ASSESSING TARGET SPECIFICITY OF THE SMALL MOLECULE INHIBITOR MARIMASTAT TO SNAKE VENOM TOXINS: A NOVEL APPLICATION OF THERMAL PROTEOME PROFILING

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation, and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including Thermal Proteome Profiling (TPP) and Proteome Integral Solubility Alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply TPP and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Identification of an Antimicrobial Peptide from the Venom of the Trinidad Thick-Tailed Scorpion Tityus trinitatis with Potent Activity against ESKAPE Pathogens and Clostridioides difficile

Envenomation by the Trinidad thick-tailed scorpion Tityus trinitatis may result in fatal myocarditis and there is a high incidence of acute pancreatitis among survivors. Peptidomic analysis (reversed-phase HPLC followed by MALDI-TOF mass spectrometry and automated Edman degradation) of T. trinitatis venom led to the isolation and characterization of three peptides with antimicrobial activity. Their primary structures were established asTtAP-1 (FLGSLFSIGSKLLPGVFKLFSRKKQ.NH2), TtAP-2 (IFGMIPGLIGGLISAFK.NH2) and TtAP-3 (FFSLIPSLIGGLVSAIK.NH2). In addition, potassium channel and sodium channel toxins, present in the venom in high abundance, were identified by CID-MS/MS sequence analysis. TtAP-1 was the most potent against a range of clinically relevant Gram-positive and Gram-negative aerobes and against the anaerobe Clostridioides difficile (MIC = 3.1-12.5 µg/mL). At a concentration of 1× MIC, TtAP-1 produced rapid cell death (<15 min against Acinetobacter baumannii and Staphylococcus aureus). The therapeutic potential of TtAP-1 as an anti-infective agent is limited by its high hemolytic activity (LC50 = 18 µg/mL against mouse erythrocytes) but the peptide constitutes a template for the design of analogs that maintain the high bactericidal activity against ESKAPE pathogens but are less toxic to human cells. It is suggested that the antimicrobial peptides in the scorpion venom facilitate the action of the neurotoxins by increasing the membrane permeability of cells from either prey or predator.

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.

Proteomic analysis of the mandibular glands from the Chinese crocodile lizard, Shinisaurus crocodilurus - Another venomous lizard?

Based on its phylogenetic relationship to monitor lizards (Varanidae), Gila monsters (Heloderma spp.), and the earless monitor Lanthanotus borneesis, the Chinese crocodile lizard, Shinisaurus crocodilurus, has been assigned to the Toxicofera clade, which comprises venomous reptiles. However, no data about composition and biological activities of its oral secretion have been reported. In the present study, a proteomic analysis of the mandibular gland of S. crocodilurus and, for comparison, of the herbivorous Solomon Island skink Corucia zebrata, was performed. Scanning electron microscopy (SEM) of the teeth from S. crocodilurus revealed a sharp ridge on the anterior surface, but no grooves, whereas those of C. zebrata possess a flattened crown with a pointed cusp. Proteomic analysis of their gland extracts provided no evidence of venom-derived peptides or proteins, strongly supporting the non-venomous character of these lizards. Data are available via ProteomeXchange with identifier PXD039424.

Monoclonal-Based Antivenomics Reveals Conserved Neutralizing Epitopes in Type I PLA2 Molecules from Coral Snakes

For over a century, polyclonal antibodies have been used to treat snakebite envenoming and are still considered by the WHO as the only scientifically validated treatment for snakebites. Nevertheless, moderate innovations have been introduced to this immunotherapy. New strategies and approaches to understanding how antibodies recognize and neutralize snake toxins represent a challenge for next-generation antivenoms. The neurotoxic activity of Micrurus venom is mainly due to two distinct protein families, three-finger toxins (3FTx) and phospholipases A₂ (PLA₂). Structural conservation among protein family members may represent an opportunity to generate neutralizing monoclonal antibodies (mAbs) against family-conserved epitopes. In this work, we sought to produce a set of monoclonal antibodies against the most toxic components of M. altirostris venom. To this end, the crude venom was fractionated, and its major toxic proteins were identified and used to generate a panel of five mAbs. The specificity of these mAbs was characterized by ELISA and antivenomics approaches. Two of the generated mAbs recognized PLA₂ epitopes. They inhibited PLA₂ catalytic activity and showed paraspecific neutralization against the myotoxicity from the lethal effect of Micrurus and Naja venoms' PLA_2s. Epitope conservation among venom PLA₂ molecules suggests the possibility of generating pan-PLA₂ neutralizing antibodies.

Combined Molecular and Elemental Mass Spectrometry Approaches for Absolute Quantification of Proteomes: Application to the Venomics Characterization of the Two Species of Desert Black Cobras, Walterinnesia aegyptia and Walterinnesia morgani

We report a novel hybrid, molecular and elemental mass spectrometry (MS) setup for the absolute quantification of snake venom proteomes shown here for two desert black cobra species within the genus Walterinnesia, Walterinnesia aegyptia and Walterinnesia morgani. The experimental design includes the decomplexation of the venom samples by reverse-phase chromatography independently coupled to four mass spectrometry systems: the combined bottom-up and top-down molecular MS for protein identification and a parallel reverse-phase microbore high-performance liquid chromatograph (RP-μHPLC) on-line to inductively coupled plasma (ICP-MS/MS) elemental mass spectrometry and electrospray ionization quadrupole time-of-flight mass spectrometry (ESI-QToF MS). This allows to continuously record the absolute sulfur concentration throughout the chromatogram and assign it to the parent venom proteins separated in the RP-μHPLC-ESI-QToF parallel run via mass profiling. The results provide a locus-resolved and quantitative insight into the three desert black cobra venom proteome samples. They also validate the units of measure of our snake venomics strategy for the relative quantification of snake venom proteomes as % of total venom peptide bonds as a proxy for the % by weight of the venom toxins/toxin families. In a more general context, our work may pave the way for broader applications of hybrid elemental/molecular MS setups in diverse areas of proteomics.

Venomics of the poorly studied hognosed pitvipers Porthidium arcosae and Porthidium volcanicum

We report the first proteomics analyses of the venoms of two poorly studied snakes, the Manabi hognosed pitviper Porthidium arcosae endemic to the western coastal province of Manabí (Ecuador), and the Costa Rican hognosed pitviper P. volcanicum with distribution restricted to South Pacific Costa Rica and western Panamá. These venom proteomes share a conserved compositional pattern reported in four other congeneric species within the clade of South American Porthidium species, P. nasutum, P. lansbergii, P. ophryomegas, and P. porrasi. The paraspecific immunorecognition profile of antivenoms produced in Costa Rica (ICP polyvalent), Perú (Instituto Nacional de Salud) and Brazil (soro antibotrópico pentavalente, SAB, from Instituto Butantan) against the venom of P. arcosae was investigated through a third-generation antivenomics approach. The maximal venom-binding capacities of the investigated antivenoms were 97.1 mg, 21.8 mg, and 25.7 mg of P. arcosae venom proteins per gram of SAB, ICP, and INS-PERU antibody molecules, respectively, which translate into 28.4 mg, 13.1 mg, and 15.2 mg of total venom proteins bound per vial of SAB, ICP, and INS-PERU AV. The antivenomics results suggest that 21.8%, 7.8% and 6.1% of the SAB, ICP, and INS-PERU antibody molecules recognized P. arcosae venom toxins. The SAB antivenom neutralized P. arcosae venom's lethality in mice with an ED₅₀ of 31.3 mgV/g SAB AV. This preclinical neutralization paraspecificity points to Brazilian SAB as a promising candidate for the treatment of envenomings by Ecuadorian P. arcosae. BIOLOGICAL SIGNIFICANCE: Assessing the preclinical efficacy profile of antivenoms against homologous and heterologous medically relevant snake venoms represents an important goal towards defining the biogeographic range of their clinical utility. This is particularly relevant in regions, such as Mesoamerica, where a small number of pharmaceutical companies produce antivenoms against the venoms of a small number of species of maximum medical relevance among the local rich herpetofauna, leaving a wide range of snakes of secondary medical relevance, but also causing life-threatening human envenomings without nominal clinical coverage. This work is part of a larger project aiming at mapping the immunological characteristics of antivenoms generated in Latin American countries towards venoms of such poorly studied snakes of the local and neighboring countries' herpetofauna. Here we report the proteomics characterization of the Manabi hognosed pitviper Porthidium arcosae endemic to the western coastal province of Manabí (Ecuador), and the Costa Rican hognosed pitviper P. volcanicum with distribution restricted to southwestern Costa Rica, the antivenomics assessment of three bothropoid commercial antivenoms produced in Costa Rica, Perú, and Brazil against the venom components of P. arcosae, and the in vivo capacity of the Brazilian soro antibotrópico pentavalente (SAB) from Instituto Butantan to neutralize the murine lethality of P. arcosae venom. The preclinical paraspecific ED₅₀ of 31.3 mg of P. arcosae venom per gram of antivenom points to Brazilian SAB as a promising candidate for the treatment of envenomings by the Manabi hognosed pitviper P. arcosae.

Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a "cocktail", as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an "ecological stance". We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies - "venomics" and "antivenomics" - can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.

Preclinical validation of a repurposed metal chelator as an early-intervention therapeutic for hemotoxic snakebite

Snakebite envenoming causes 138,000 deaths annually, and ~400,000 victims are left with permanent disabilities. Envenoming by saw-scaled vipers (Viperidae: Echis) leads to systemic hemorrhage and coagulopathy and represents a major cause of snakebite mortality and morbidity in Africa and Asia. The only specific treatment for snakebite, antivenom, has poor specificity and low affordability and must be administered in clinical settings because of its intravenous delivery and high rates of adverse reactions. This requirement results in major treatment delays in resource-poor regions and substantially affects patient outcomes after envenoming. Here, we investigated the value of metal ion chelators as prehospital therapeutics for snakebite. Among the tested chelators, dimercaprol (British anti-Lewisite) and its derivative 2,3-dimercapto-1-propanesulfonic acid (DMPS) were found to potently antagonize the activity of Zn²⁺-dependent snake venom metalloproteinases in vitro. Moreover, DMPS prolonged or conferred complete survival in murine preclinical models of envenoming against a variety of saw-scaled viper venoms. DMPS also considerably extended survival in a "challenge and treat" model, where drug administration was delayed after venom injection and the oral administration of this chelator provided partial protection against envenoming. Last, the potential clinical scenario of early oral DMPS therapy combined with a delayed, intravenous dose of conventional antivenom provided prolonged protection against the lethal effects of envenoming in vivo. Our findings demonstrate that the safe and affordable repurposed metal chelator DMPS can effectively neutralize saw-scaled viper venoms in vitro and in vivo and highlight the promise of this drug as an early, prehospital, therapeutic intervention for hemotoxic snakebite envenoming.

Seminal Plasma: Relevant for Fertility?

Seminal plasma (SP), the non-cellular component of semen, is a heterogeneous composite fluid built by secretions of the testis, the epididymis and the accessory sexual glands. Its composition, despite species-specific anatomical peculiarities, consistently contains inorganic ions, specific hormones, proteins and peptides, including cytokines and enzymes, cholesterol, DNA and RNA-the latter often protected within epididymis- or prostate-derived extracellular vesicles. It is beyond question that the SP participates in diverse aspects of sperm function pre-fertilization events. The SP also interacts with the various compartments of the tubular genital tract, triggering changes in gene function that prepares for an eventual successful pregnancy; thus, it ultimately modulates fertility. Despite these concepts, it is imperative to remember that SP-free spermatozoa (epididymal or washed ejaculated) are still fertile, so this review shall focus on the differences between the in vivo roles of the SP following semen deposition in the female and those regarding additions of SP on spermatozoa handled for artificial reproduction, including cryopreservation, from artificial insemination to in vitro fertilization. This review attempts, including our own results on model animal species, to critically summarize the current knowledge of the reproductive roles played by SP components, particularly in our own species, which is increasingly affected by infertility. The ultimate goal is to reconcile the delicate balance between the SP molecular concentration and their concerted effects after temporal exposure in vivo. We aim to appraise the functions of the SP components, their relevance as diagnostic biomarkers and their value as eventual additives to refine reproductive strategies, including biotechnologies, in livestock models and humans.

Comparative venomics and preclinical efficacy evaluation of a monospecific Hemachatus antivenom towards sub-Saharan Africa cobra venoms

Cobras are the most medically important elapid snakes in Africa. The African genera Naja and Hemachatus include snakes with neurotoxic and cytotoxic venoms, with shared biochemical, toxinological and antigenic characteristics. We have studied the antigenic cross-reactivity of four sub-Saharan Africa cobra venoms against an experimental monospecific Hemachatus haemachatus antivenom through comparative proteomics, preclinical assessment of neutralization, and third generation antivenomics. The venoms of H. haemachatus, N. annulifera, N. mossambica and N. nigricollis share an overall qualitative family toxin composition but depart in their proportions of three-finger toxin (3FTxs) classes, phospholipases A₂ (PLA₂s), snake venom metalloproteinases (SVMPs), and cysteine-rich secretory proteins (CRISPs). A monospecific anti-Hemachatus antivenom produced by Costa Rican Instituto Clodomiro Picado neutralized the lethal activity of the homologous and heterologous neuro/cytotoxic (H. haemachatus) and cyto/cardiotoxic (N. mossambica and N. nigricollis) venoms of the three spitting cobras sampled, while it was ineffective against the lethal and toxic activities of the neurotoxic venom of the non-spitting snouted cobra N. annulifera. The ability of the anti-Hemachatus-ICP antivenom to neutralize toxic (dermonecrotic and anticoagulant) and enzymatic (PLA₂) activities of spitting cobra venoms suggested a closer kinship of H. haemachatus and Naja subgenus Afrocobra spitting cobras than to Naja subgenus Uraeus neurotoxic taxa. These results were confirmed by third generation antivenomics. BIOLOGICAL SIGNIFICANCE: African Naja species represent the most widespread medically important elapid snakes across Africa. To gain deeper insight into the spectrum of medically relevant toxins, we compared the proteome of three spitting cobras (Hemachatus haemachatus, Naja mossambica and N. nigricollis) and one non-spitting cobra (N. annulifera). Three finger toxins and phospholipases A₂ are the two major protein families among the venoms analyzed. The development of antivenoms of broad species coverage is an urgent need in sub-Saharan Africa. An equine antivenom raised against H. haemachatus venom showed cross-reactivity with the venoms of H. haemachatus, N. mossambica and N. nigricollis, while having poor recognition of the venom of N. annulifera. This immunological information provides clues for the design of optimum venom mixtures for the preparation of broad spectrum antivenoms.

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

Convergent evolution provides insights into the selective drivers underlying evolutionary change. Snake venoms, with a direct genetic basis and clearly defined functional phenotype, provide a model system for exploring the repeated evolution of adaptations. While snakes use venom primarily for predation, and venom composition often reflects diet specificity, three lineages of cobras have independently evolved the ability to spit venom at adversaries. Using gene, protein, and functional analyses, we show that the three spitting lineages possess venoms characterized by an up-regulation of phospholipase A2 (PLA2) toxins, which potentiate the action of preexisting venom cytotoxins to activate mammalian sensory neurons and cause enhanced pain. These repeated independent changes provide a fascinating example of convergent evolution across multiple phenotypic levels driven by selection for defense.

The earless monitor lizard Lanthanotus borneensis - A venomous animal?

Based on its mandibular gland secretion, the earless monitor lizard, Lanthanotus borneensis, has been considered a venomous animal like other members of the Toxicofera group, including Heloderma. In the present study, the gland structure and teeth of L. borneensis were examined by micro-tomography (μCT) and scanning electron microscopy (SEM), respectively, and proteomic analysis of the gland extract was performed. The mandibular gland consists of six compartments with separate ducts. The pleurodont teeth of the lower and upper jaw are not grooved but possess a sharp ridge on the anterior surface. Proteomic analysis of the gland extract confirmed previous studies that kallikrein enzymes are the major biologically active components. In view of the lizard's biology, its mandibular gland secretion is obviously not needed for prey capture or defence. It seems not justified the labelling of L. borneensis as a venomous animal. However, definitively answering this question requires toxinological studies on natural prey.

The molecular basis of venom resistance in a rattlesnake-squirrel predator-prey system

Understanding how interspecific interactions mould the molecular basis of adaptations in coevolving species is a long-sought goal of evolutionary biology. Venom in predators and venom resistance proteins in prey are coevolving molecular phenotypes, and while venoms are highly complex mixtures it is unclear if prey respond with equally complex resistance traits. Here, we use a novel molecular methodology based on protein affinity columns to capture and identify candidate blood serum resistance proteins ("venom interactive proteins" [VIPs]) in California Ground Squirrels (Otospermophilus beecheyi) that interact with venom proteins from their main predator, Northern Pacific Rattlesnakes (Crotalus o. oreganus). This assay showed that serum-based resistance is both population- and species-specific, with serum proteins from ground squirrels showing higher binding affinities for venom proteins of local snakes compared to allopatric individuals. Venom protein specificity assays identified numerous and diverse candidate prey resistance VIPs but also potential targets of venom in prey tissues. Many specific VIPs bind to multiple snake venom proteins and, conversely, single venom proteins bind multiple VIPs, demonstrating that a portion of the squirrel blood serum "resistome" involves broad-based inhibition of nonself proteins and suggests that resistance involves a toxin scavenging mechanism. Analyses of rates of evolution of VIP protein homologues in related mammals show that most of these proteins evolve under purifying selection possibly due to molecular constraints that limit the evolutionary responses of prey to rapidly evolving snake venom proteins. Our method represents a general approach to identify specific proteins involved in co-evolutionary interactions between species at the molecular level.

Comparative proteomic profiling and functional characterization of venom pooled from captive Crotalus durissus terrificus specimens and the Brazilian crotalic reference venom

The South American rattlesnake Crotalus durissus spp has a wide geographic distribution in Brazil. Although responsible for only a low proportion of ophidian accidents, it is considered one of the most medically important species of venomous snakes due to the high mortality rate (1.87%). Snake venom is a complex phenotype commonly subjected to individual intraspecific, ontogenetic and geographic variability. Compositional differences in pooled venom used in the immunization process may impact the efficacy of the antivenom. In order to assure standardized high-quality antivenom, the potency of each Brazilian crotalic antivenom batch is determined against the 'Brazilian Crotalic Reference Venom' (BCRV). BCRV is produced by Instituto Butantan using venom obtained from the first milking of recently wild-caught C. d. terrificus specimens brought to the Institute. The decrease in the number of snake donations experienced in recent years can become a threat to the production of future batches of BCRV. To evaluate the feasibility of using venom from long-term captive animals in the formulation of BCRV, we have compared the proteomic, biochemical and biological profiles of C. d. terrificus venom pooled from captive specimens (CVP- captive venom pool) and BCRV. Electrophoretic and venomics analyses revealed a very similar venom composition profile, but also certain differences in toxins abundance, with some low abundant protein families found only in BCRV. Enzymatic (L-amino acid oxidase, phospholipase A₂ and proteolytic) and biological (myotoxic and coagulant) activities showed higher values in CVP than in BCRV. CVP also possessed slightly higher lethal effect, although the Instituto Butantan crotalic antivenom showed equivalent potency neutralizing BCRV and CVP. Our results strongly suggest that venom from long-term captive C. d. terrificus might be a valid alternative to generate an immunization mixture of equivalent quality to the currently in use reference venom.

Comparative characterization of Viperidae snake venoms from Perú reveals two compositional patterns of phospholipase A2 expression

Snake species within the Bothrops complex (sensu lato) are of medical relevance in Latin America, but knowledge on their venom characteristics is limited, or even unavailable, for some taxa. Perú harbors 17 species of pit vipers, within the genera Bothrops, Bothriechis, Bothrocophias, Porthidium, Crotalus, and Lachesis. This study compared the venoms of twelve species, through chromatographic and electrophoretic profiles, as well as proteolytic and phospholipase A₂ (PLA₂) activities. Also, proteomic profiles were analyzed for nine of the venoms using a shotgun approach. Results unveiled conspicuous differences in the expression of venom PLA₂s among species, six of them presenting scarce levels as judged by RP-HPLC profiles. Since most species within the bothropoid lineage possess venoms with high to intermediate abundances of this protein family, our findings suggest the existence of a phenotypic duality in the expression of venom PLA₂s within the Bothrops (sensu lato) complex. Bothrops barnetti and Bothrocophias andianus venoms, very scarce in PLA₂s, were shown to lack significant myotoxic activity, highlighting that the observed variability in PLA₂ expression bears toxicological correlations with effects attributed to these proteins. Finally, an attempt to identify phylogenetic relationships of bothropoid species from Perú presenting low- or high-PLA₂ venom phenotypes showed an interspersed pattern, thus precluding a simple phylogenetic interpretation of this venom compositional dichotomy.

•

Venoms from 12 viperids of Perú were compared.

•

Conspicuous differences in the expression of PLA₂ were found.

•

Venoms presenting scarce levels of PLA₂ lack myotoxicity.

•

A new phenotypic dichotomy in venom PLA₂ expression is described within Bothrops (sensu lato).

Dagestan blunt-nosed viper, Macrovipera lebetina obtusa (Dwigubsky, 1832), venom. Venomics, antivenomics, and neutralization assays of the lethal and toxic venom activities by anti-Macrovipera lebetina turanica and anti-Vipera berus berus antivenoms

We have applied a combination of venomics, in vivo neutralization assays, and in vitro third-generation antivenomics analysis to assess the preclinical efficacy of the monospecific anti-Macrovipera lebetina turanica (anti-Mlt) antivenom manufactured by Uzbiopharm® (Uzbekistan) and the monospecific anti-Vipera berus berus antivenom from Microgen® (Russia) against the venom of Dagestan blunt-nosed viper, Macrovipera lebetina obtusa (Mlo). Despite their low content of homologous (anti-Mlt, 5–10%) or para-specific (anti-Vbb, 4–9%) F(ab')₂ antibody fragments against M. l. obtusa venom toxins, both antivenoms efficiently recognized most components of the complex venom proteome's arsenal, which is made up of toxins derived from 11 different gene families and neutralized, albeit at different doses, key toxic effects of M. l. obtusa venom, i.e., in vivo lethal and hemorrhagic effects in a murine model, and in vitro phospholipase A₂, proteolytic and coagulant activities. The calculated lethality neutralization potencies for Uzbiopharm® anti-Mlt and anti-Vbb Microgen® antivenoms were 1.46 and 1.77 mg/mL, indicating that 1 mL of Uzbiopharm® and Microgen® antivenoms may protect mice from 41 to 50 LD₅₀s of Mlo venom, respectively. The remarkable degree of conservation of immunogenic determinants between species of the clades of European and Oriental viper, which evolved geographically segregated since the early Miocene, suggests an eventual window of opportunity for the treatment of envenomings by Eurasian snakes. Clearly, the rational use of heterologous antivenoms requires establishing their para-specificity landscapes. This paper illustrates the analytical power of combining in vitro and in vivo preclinical quantitative assays toward this goal.

•

Efficacy against M. l. obtusa venom by two antivenoms was investigated.

•

Quantification of lethality neutralizing antibodies was assessed.

•

Anti-Vipera berus antivenom showed paraspecificity against M. l. obtusa venom.

•

This study provides hints as how to improve the potency of the antivenoms sampled.

Venomics and antivenomics of the poorly studied Brazil's lancehead, Bothrops brazili (Hoge, 1954), from the Brazilian State of Pará

Background:

The Brazil’s lancehead, Bothrops brazili, is a poorly studied pit viper distributed in lowlands of the equatorial rainforests of southern Colombia, northeastern Peru, eastern Ecuador, southern and southeastern Venezuela, Guyana, Suriname, French Guiana, Brazil, and northern Bolivia. Few studies have been reported on toxins isolated from venom of Ecuadorian and Brazilian B. brazili. The aim of the present study was to elucidate the qualitative and quantitative protein composition of B. brazili venom from Pará (Brazil), and to carry out a comparative antivenomics assessment of the immunoreactivity of the Brazilian antibothropic pentavalent antivenom [soro antibotrópico (SAB) in Portuguese] against the venoms of B. brazili and reference species, B. jararaca.

Methods:

We have applied a quantitative snake venomics approach, including reverse-phase and two-dimensional electrophoretic decomplexation of the venom toxin arsenal, LC-ESI-MS mass profiling and peptide-centric MS/MS proteomic analysis, to unveil the overall protein composition of B. brazili venom from Pará (Brazil). Using third-generation antivenomics, the specific and paraspecific immunoreactivity of the Brazilian SAB against homologous (B. jararaca) and heterologous (B. brazili) venoms was investigated.

Results:

The venom proteome of the Brazil’s lancehead (Pará) is predominantly composed of two major and three minor acidic (19%) and two major and five minor basic (14%) phospholipase A₂ molecules; 7-11 snake venom metalloproteinases of classes PI (21%) and PIII (6%); 10-12 serine proteinases (14%), and 1-2 L-amino acid oxidases (6%). Other toxins, including two cysteine-rich secretory proteins, one C-type lectin-like molecule, one nerve growth factor, one 5'-nucleotidase, one phosphodiesterase, one phospholipase B, and one glutaminyl cyclase molecule, represent together less than 2.7% of the venom proteome. Third generation antivenomics profile of the Brazilian pentabothropic antivenom showed paraspecific immunoreactivity against all the toxin classes of B. brazili venom, with maximal binding capacity of 132.2 mg venom/g antivenom. This figure indicates that 19% of antivenom's F(ab')₂ antibodies bind B. brazili venom toxins.

Conclusion:

The proteomics outcome contribute to a deeper insight into the spectrum of toxins present in the venom of the Brazil’s lancehead, and rationalize the pathophysiology underlying this snake bite envenomings. The comparative qualitative and quantitative immunorecognition profile of the Brazilian pentabothropic antivenom toward the venom toxins of B. brazili and B. jararaca (the reference venom for assessing the bothropic antivenom's potency in Brazil), provides clues about the proper use of the Brazilian antibothropic polyvalent antivenom in the treatment of bites by the Brazil’s lancehead.

Isolation and characterization of cytotoxic and insulin-releasing components from the venom of the black-necked spitting cobra Naja nigricollis (Elapidae)

Four peptides with cytotoxic activity against BRIN-BD11 rat clonal β-cells were purified from the venom of the black-necked spitting cobra Naja nigricollis using reversed-phase HPLC. The peptides were identified as members of the three-finger superfamily of snake toxins by ESI-MS/MS sequencing of tryptic peptides. The most potent peptide (cytotoxin-1N) showed strong cytotoxic activity against three human tumor-derived cell lines (LC₅₀ = 0.8 ± 0.2 μM for A549 non-small cell lung adenocarcinoma cells; LC₅₀ = 7 ± 1 μM for MDA-MB-231 breast adenocarcinoma cells; and LC₅₀ = 9 ± 1 μM for HT-29 colorectal adenocarcinoma cells). However, all the peptides were to varying degrees cytotoxic against HUVEC human umbilical vein endothelial cells (LC₅₀ in the range 2–22 μM) and cytotoxin-2N was moderately hemolytic (LC₅₀ = 45 ± 3 μM against mouse erythrocytes). The lack of differential activity against cells derived from non-neoplastic tissue limits their potential for development into anti-cancer agents. In addition, two proteins in the venom, identified as isoforms of phospholipase A₂, effectively stimulated insulin release from BRIN-BD11 cells (an approximately 6-fold increase in rate compared with 5.6 mM glucose alone) at a concentration (1 μM) that was not cytotoxic to the cells suggesting possible application in therapy for Type 2 diabetes.

•

Four members of the three-finger superfamily of toxins were isolated from N. nigricollis venom.

•

All peptides were cytotoxic to human tumor-derived cells.

•

The peptides were also cytotoxic to non-neoplastic HUVEC cells.

•

Two isoforms of phospholipase A₂ effectively stimulated insulin release from rat clonal β-cells.

When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

Phylovenomics of Daboia russelii across the Indian subcontinent. Bioactivities and comparative in vivo neutralization and in vitro third-generation antivenomics of antivenoms against venoms from India, Bangladesh and Sri Lanka

Russell's viper (Daboia russelii) is, together with Naja naja, Bungarus caeruleus and Echis carinatus, a member of the medically important 'Big Four' species responsible for causing a large number of morbidity and mortality cases across the Indian subcontinent. Despite the wide distribution of Russell's viper and the well-documented ubiquity of the phenomenon of geographic variability of intraspecific snake venom composition, Indian polyvalent antivenoms against the "Big Four" venoms are raised against venoms sourced mainly from Chennai in the southeastern Indian state of Tamil Nadu. Biochemical and venomics investigations have consistently revealed notable compositional, functional, and immunological differences among geographic variants of Russell's viper venoms across the Indian subcontinent. However, these studies, carried out by different laboratories using different protocols and involving venoms from a single geographical region, make the comparison of the different venoms difficult. To bridge this gap, we have conducted bioactivities and proteomic analyses of D. russelii venoms from the three corners of the Indian subcontinent, Pakistan, Bangladesh, and Tamil Nandu (India) and Sri Lanka, along with comparative in vivo neutralization and in vitro third-generation antivenomics of antivenoms used in India, Bangladesh and Sri Lanka. These analyses let us to propose two alternative routes of radiation for Russell's viper in the Indian subcontinent. Both radiations, towards the northeast of India and Bangladesh and towards south India and Sri Lanka, have a common origin in Pakistan, and provide a phylovenomics ground for rationalizing the geographic variability in venom composition and their distinct immunoreactivity against available antivenoms. BIOLOGICAL SIGNIFICANCE: Russell's viper (Daboia russelii), the Indian cobra (Naja naja), the common krait (Bungarus caeruleus), and the saw-scaled viper (Echis carinatus) constitute the 'Big Four' snake species responsible for most snakebite envenomings and deaths in the Indian subcontinent. Despite the medical relevance of Daboia russelii, and the well documented variations in the clinical manifestations of envenomings by this wide distributed species, which are doubtless functionally related to differences in venom composition of its geographic variants, antivenoms for the clinical treatment of envenomings by D. russelii across the Indian subcontinent are invariably raised using venom sourced mainly from the southeastern Indian state of Tamil Nadu. We have applied a phylovenomics approach to compare the venom proteomes of Russell's vipers from the three corners of the Indian subcontinent, Pakistan, Bangladesh, and South India/Sri Lanka, and have assessed the in vitro (third-generation antivenomics) and in vivo preclinical efficacy of a panel of homologous antivenoms. The identification of two dispersal routes of ancestral D. russelii into the Indian subcontinent provides the ground for rationalizing the variability in composition and immunoreactivity of the venoms of extant geographic variants of Russell's viper. Such knowledge is relevant for envisioning strategies to improve the clinical coverage of anti- D. russelii antivenoms.

Third-generation antivenomics analysis of the preclinical efficacy of Bothrofav® antivenom towards Bothrops lanceolatus venom

Bothrops lanceolatus inflicts severe envenomings in the Lesser Caribbean island of Martinique. Bothrofav®, a monospecific antivenom against B. lanceolatus venom, has proven highly effective at the preclinical and clinical levels. Here, we report a detailed third-generation antivenomics quantitative analysis of Bothrofav®. With the exception of poorly-immunogenic peptides, Bothrofav® immunocaptured all the major protein components. These results, along with previous preclinical and clinical observations, underscore the high neutralizing efficacy of the antivenom against B. lanceolatus venom.

•

The immunoreactivity of the monospecific antivenom Bothrofav® was evaluated against Bothrops lanceolatus venom.

•

Third Generation Antivenomics approach was used in this study.

•

Antivenom immunocaptured all major venom components.

•

Results agree with previous preclinical studies and demonstrate the efficacy of this antivenom.

The Harderian gland transcriptomes of Caraiba andreae, Cubophis cantherigerus and Tretanorhinus variabilis, three colubroid snakes from Cuba

The Harderian gland is a cephalic structure, widely distributed among vertebrates. In snakes, the Harderian gland is anatomically connected to the vomeronasal organ via the nasolacrimal duct, and in some species can be larger than the eyes. The function of the Harderian gland remains elusive, but it has been proposed to play a role in the production of saliva, pheromones, thermoregulatory lipids and growth factors, among others. Here, we have profiled the transcriptomes of the Harderian glands of three non-front-fanged colubroid snakes from Cuba: Caraiba andreae (Cuban Lesser Racer); Cubophis cantherigerus (Cuban Racer); and Tretanorhinus variabilis (Caribbean Water Snake), using Illumina HiSeq2000 100 bp paired-end. In addition to ribosomal and non-characterized proteins, the most abundant transcripts encode putative transport/binding, lipocalin/lipocalin-like, and bactericidal/permeability-increasing-like proteins. Transcripts coding for putative canonical toxins described in venomous snakes were also identified. This transcriptional profile suggests a more complex function than previously recognized for this enigmatic organ.

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.

Clues for Understanding the Structure and Function of a Prototypic Human Integrin: The Platelet Glycoprotein IIb/IIIa Complex

The glycoprotein (GP) IIb/IIIa, a Ca2+-dependent heterodimer, is the major integrin on the platelet plasma membrane. On resting platelets GPIIb/IIIa is maintained in an inactive conformation and serves as a low affinity adhesion receptor for surface-coated fibrinogen, whereas upon platelet activation signals within the cytoplasma alter the receptor function of GPIIb/IIIa (inside-out signalling), which undergoes a measurable conformational change within its exoplasmic domains, and becomes a competent receptor for soluble fibrinogen and some other RGD sequence-containing plasma adhesive proteins. Upon ligand binding, further structural alterations trigger the association of receptor-occupied GPIIb/IIIa complexes with themselves within the plane of the membrane. The simultaneous binding of dimeric fibrinogen molecules to GPIIb/IIIa clusters on adjacent platelets leads to platelet aggregation, which promotes attachment of fibrinogen-GPIIb/IIIa clusters to the cytoskeleton (outside-in signalling). This, in turn, provides the necessary physical link for clot retraction to occur, and generates a cascade of intracellular biochemical reactions which result in the formation of a multiprotein signalling complex at the cytoplasmic domains of GPIIb/IIIa. Glycoprotein IMIIa, also called αIIbβ3 in the integrin nomenclature, plays thus a primary role in both platelet adhesion and thrombus formation at the site of vascular injury. In addition, the human glycoprotein Ilb/IIIa complex is the most thoroughly studied integrin receptor, its molecular biology and major features of its primary structure having been elucidated mainly during the last six years. Furthermore, localization of functionally relevant monoclonal antibody epitopes, determination of the cross-linking sites of inhibitory peptide ligands, proteolytic dissection of the isolated integrin, and analysis of natural and artificial GPIIb/IIIa mutants have recently provided a wealth of information regarding structure-function relationships of human GPIIb/IIIa. The aim of this review is to summarize these many structural and functional data in the perspective of an emerging model. Although most of the interpretations based on structural elements of this initial biochemical model require independent confirmation, they may help us to understand the structure-function relationship of this major platelet receptor, and of other members of the integrin superfamily, as well as to perform further investigations in order to test current hypotheses.

Identification of a Glycoprotein III a Dimer in Polyacrylamide Gel Separations of Human Platelet Membranes

Membrane glycoproteins IIb and IIIa play a major role in human blood platelet aggregation. The absence or the severe reduction of these two membrane glycoproteins, as observed in platelets of Glanzmann’s thrombasthenic patients, is related to a lack of platelet aggregation. Separation of Glanzmann’s thrombasthenic platelet samples by two-dimensional polyacrylamide O’Farrell gels show the absence of a high and several low molecular mass glycoproteins, in addition to the loss of glycoproteins IIb and IIIa (McGregor J. L. et al. Eur. J. Biochem. 1981; 116: 379-388). The aim of this study was to identify the nature of the high molecular mass component, absent in thrombasthenic platelets. A high molecular mass glycoprotein (200 kDa), present in two-dimensional SDS-polyacrylamide O-Farrell gel separations, was recognized by a monoclonal antibody (MP37) directed against glycoprotein IIIa. Moreover, the tryptic peptide map of this high molecular mass glycoprotein was nearly identical to that of glycoprotein IIIa. These results indicate that this high molecular mass glycoprotein present in SDS-polyacrylamide gels is a dimer of glycoprotein IIIa. This work raises the possibility that the high molecular mass glycoprotein, absent in two-dimensional O’Farrell gel separations of thrombasthenic platelets, is a dimer of glycoprotein IIIa.

Toxin-resolved antivenomics-guided assessment of the immunorecognition landscape of antivenoms

Snakebite envenoming represents a major issue in rural areas of tropical and subtropical regions across sub-Saharan Africa, South to Southeast Asia, Latin America and Oceania. Antivenoms constitute the only scientifically validated therapy for snakebite envenomings, provided they are safe, effective, affordable, accessible and administered appropriately. However, the lack of financial incentives in a technology that has remained relatively unchanged for more than a century, has contributed to some manufacturers leaving the market and others downscaling production or increasing the prices, leading to a decline in the availability and accessibility for these life-saving antidotes to millions of rural poor most at risk from snakebites in low income countries. The shortage of antivenoms can be significantly alleviated by optimizing the use of current antivenoms (through the assessment of their specific and paraspecific efficacy against the different medically relevant homologous and heterologous snake venoms) and by generating novel polyspecific antivenoms exhibiting broad clinical spectrum and wide geographic distribution range. Research on venoms has been continuously enhanced by advances in technology. Particularly, the last decade has witnessed the development of omics strategies for unravelling the toxin composition of venoms ("venomics") and to assess the immunorecognition profile of antivenoms ("antivenomics"). Here, we review recent developments and reflect on near future innovations that promise to revolutionize the mutually enlightening relationship between evolutionary and translational venomics.

The paraspecific neutralisation of snake venom induced coagulopathy by antivenoms

Snake envenoming causes several potentially lethal pathologies. The specific pathology is dictated by the toxin composition of venom, which varies by species, geography and ontogeny. This variation severely restricts the paraspecific efficacy of antivenoms used to treat snakebite victims. With a view to devising pathology-specific snakebite treatments, we assessed the procoagulant activity of 57 snake venoms and investigated the efficacy of various antivenoms. We find that procoagulant venoms act differentially on key steps of the coagulation cascade, and that certain monospecific antivenoms work in a previously unrecognised paraspecific manner to neutralise this activity, despite conventional assumptions of congener-restricted efficacy. Moreover, we demonstrate that the metal chelator EDTA is also capable of neutralising venom-induced lethality in vivo. This study illustrates the exciting potential of developing new, broad-spectrum, toxin-targeting antivenoms capable of treating key snakebite pathologies, and advocates a thorough re-examination of enzyme inhibiting compounds as alternative therapies for treating snakebite victims.

A synthetic biology approach for consistent production of plant-made recombinant polyclonal antibodies against snake venom toxins

Antivenoms developed from the plasma of hyperimmunized animals are the only effective treatment available against snakebite envenomation but shortage of supply contributes to the high morbidity and mortality toll of this tropical disease. We describe a synthetic biology approach to affordable and cost-effective antivenom production based on plant-made recombinant polyclonal antibodies (termed pluribodies). The strategy takes advantage of virus superinfection exclusion to induce the formation of somatic expression mosaics in agroinfiltrated plants, which enables the expression of complex antibody repertoires in a highly reproducible manner. Pluribodies developed using toxin-binding genetic information captured from peripheral blood lymphocytes of hyperimmunized camels recapitulated the overall binding activity of the immune response. Furthermore, an improved plant-made antivenom (plantivenom) was formulated using an in vitro selected pluribody against Bothrops asper snake venom toxins and has been shown to neutralize a wide range of toxin activities and provide protection against lethal venom doses in mice.

Transcriptomics-guided bottom-up and top-down venomics of neonate and adult specimens of the arboreal rear-fanged Brown Treesnake, Boiga irregularis, from Guam

The Brown Treesnake (Boiga irregularis) is an arboreal, nocturnal, rear-fanged venomous snake native to northern and eastern regions of Australia, Papua New Guinea and the Solomon Islands. It was inadvertently introduced onto the island of Guam during the late 1940's to early 1950's, and it has caused massive declines and extirpations of the native bird, lizard, and mammal populations. In the current study, we report the characterization of the venom proteome of an adult and a neonate B. irregularis specimens from Guam by a combination of venom gland transcriptomic and venomic analyses. Venom gland transcriptomic analysis of an adult individual identified toxins belonging to 18 protein families, with three-finger toxin isoforms being the most abundantly expressed transcripts, comprising 94% of all venom protein transcript reads. Transcripts for PIII-metalloproteinases, C-type lectins, cysteine-rich secretory proteins, acetylcholinesterases, natriuretic peptides, ficolins, phospholipase A₂ (PLA₂) inhibitors, PLA₂s, vascular endothelial growth factors, Kunitz-type protease inhibitors, cystatins, phospholipase Bs, cobra venom factors, waprins, SVMP inhibitors, matrix metalloproteinases, and hyaluronidases were also identified, albeit, at very low abundances ranging from 0.05% to 1.7% of the transcriptome. The venom proteomes of neonate and adult B. irregularis were also both overwhelmingly (78 and 84%, respectively) dominated by monomeric and dimeric 3FTxs, followed by moderately abundant (21% (N) and 13% (A)) CRISPs, low abundance (1% (N), 3% (A)) PIII-SVMPs, and very low abundance (<0.01%) PLA₂ and SVMP inhibitors. The differences in relative toxin abundances identified between neonate and adult snakes likely correlates to shifts in prey preference between the two age classes, from nearly-exclusively lizards to lizards, birds and small mammals. Immunoaffinity antivenomics with experimentally designed rabbit anti-Brown Treesnake (anti-BTS) venom IgGs against homologous venom from adult snakes demonstrated that CRISPs, PIII-SVMPs, and 60-70% of 3FTxs were effectively immunocaptured. Western blot analysis showed that all venom proteins were recognized by anti-BTS IgGs, and cross-reactivity with other rear-fanged snake venoms was also observed. Incubation of anti-BTS venom IgGs with crude B. irregularis venom resulted in a significant decrease in proteolytic (SVMP) activity against azocasein. These results provide the first comparative venomic and anti-venomic analysis of neonate and adult B. irregularis from Guam, further highlighting evolutionary trends in venom composition among rear-fanged venomous snakes.

SIGNIFICANCE PARAGRAPH

The Brown Treesnake (Boiga irregularis) has caused extensive ecological and economic damage to the island of Guam where it has become a classic example of the negative impacts of invasive species. In the current study, we report the first combined transcriptomic and proteomic analysis of B. irregularis venom of Guam origin. The transcriptome of an adult snake contained toxin sequences belonging to 18 protein families, with three-finger toxin (3FTx) isoforms being the most abundant and representing 94% of all venom protein transcript reads. Our bottom-up and top-down venomic analyses confirmed that 3FTxs are the major components of B. irregularis venom, and a comparative analysis of neonate and adult venoms demonstrate a clear ontogenetic shift in toxin abundance, likely driven by dietary variation between the two age classes. Second-generation antivenomics and Western blot analysis using purified anti-Brown Treesnake rabbit serum IgGs (anti-BTS IgGs) showed strong immunoreactivity toward B. irregularis venom. Interestingly, our anti-BTS IgGs did not cross-react with 3FTxs found in several other rear-fanged snake venoms, or against 3FTxs in the venom of the elapid Ophiophagus hannah, indicating that epitopes in these 3FTx molecules are quite distinct.

Snakebite envenoming

Snakebite envenoming is a neglected tropical disease that kills >100,000 people and maims >400,000 people every year. Impoverished populations living in the rural tropics are particularly vulnerable; snakebite envenoming perpetuates the cycle of poverty. Snake venoms are complex mixtures of proteins that exert a wide range of toxic actions. The high variability in snake venom composition is responsible for the various clinical manifestations in envenomings, ranging from local tissue damage to potentially life-threatening systemic effects. Intravenous administration of antivenom is the only specific treatment to counteract envenoming. Analgesics, ventilator support, fluid therapy, haemodialysis and antibiotic therapy are also used. Novel therapeutic alternatives based on recombinant antibody technologies and new toxin inhibitors are being explored. Confronting snakebite envenoming at a global level demands the implementation of an integrated intervention strategy involving the WHO, the research community, antivenom manufacturers, regulatory agencies, national and regional health authorities, professional health organizations, international funding agencies, advocacy groups and civil society institutions.

The medical threat of mamba envenoming in sub-Saharan Africa revealed by genus-wide analysis of venom composition, toxicity and antivenomics profiling of available antivenoms

Mambas (genus Dendroaspis) are among the most feared of the medically important elapid snakes found in sub-Saharan Africa, but many facets of their biology, including the diversity of venom composition, remain relatively understudied. Here, we present a reconstruction of mamba phylogeny, alongside genus-wide venom gland transcriptomic and high-resolution top-down venomic analyses. Whereas the green mambas, D. viridis, D. angusticeps, D. j. jamesoni and D. j. kaimosae, express 3FTx-predominant venoms, black mamba (D. polylepis) venom is dominated by dendrotoxins I and K. The divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas makes it plausible that this major difference in venom composition is due to dietary variation. The pattern of intrageneric venom variability across Dendroaspis represented a valuable opportunity to investigate, in a genus-wide context, the variant toxicity of the venom, and the degree of paraspecific cross-reactivity between antivenoms and mamba venoms. To this end, the immunological profiles of the five mamba venoms were assessed against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview of which available antivenoms may be more efficacious in neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the notably different potency and pharmacological profiles of Dendroaspis venoms at locus resolution. This understanding will allow selection and design of toxin immunogens with a view to generating a safer and more efficacious pan-specific antivenom against any mamba envenomation.

BIOLOGICAL SIGNIFICANCE

The mambas (genus Dendroaspis) comprise five especially notorious medically important venomous snakes endemic to sub-Saharan Africa. Their highly potent venoms comprise a high diversity of pharmacologically active peptides, including extremely rapid-acting neurotoxins. Previous studies on mamba venoms have focused on the biochemical and pharmacological characterisation of their most relevant toxins to rationalize the common neurological and neuromuscular symptoms of envenomings caused by these species, but there has been little work on overall venom composition or comparisons between them. Only very recently an overview of the composition of the venom of two Dendroaspis species, D. angusticeps and D. polylepis, has been unveiled through venomics approaches. Here we present the first genus-wide transcriptomic-proteomic analysis of mamba venom composition. The transcriptomic analyses described in this paper have contributed 29 (D. polylepis), 23 (D. angusticeps), 40 (D. viridis), 25 (D. j. jamesoni) and 21 (D. j. kaimosae), novel full-length toxin sequences to the non-redundant Dendroaspis sequence database. The mamba genus-wide venomic analysis demonstrated that major D. polylepis venom components are Kunitz-fold family toxins. This feature is unique in relation to the relatively conserved three-finger toxin (3FTx)-dominated venom compositions of the green mambas. Venom variation was interpreted in the context of dietary variation due to the divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas. Additionally, the degree of cross-reactivity conservation of mamba venoms was assessed by antivenomics against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview to infer which available antivenoms may be capable of neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the pharmacological profiles of mamba venoms at locus resolution. This understanding will contribute to the generation of a safer and more efficacious pan-Dendroaspis therapeutic antivenom against any mamba envenomation.

Integrated Venomics and Venom Gland Transcriptome Analysis of Juvenile and Adult Mexican Rattlesnakes Crotalus simus, C. tzabcan, and C. culminatus Revealed miRNA-modulated Ontogenetic Shifts

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA₂ crotoxin and acidic PLA₂ mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

Absolute venomics: Absolute quantification of intact venom proteins through elemental mass spectrometry

We report the application of a hybrid element and molecular MS configuration for the parallel absolute quantification of μHPLC-separated intact sulfur-containing venom proteins, via ICP triple quadrupole MS and ³²S/³⁴S isotope dilution analysis, and identification by ESI-QToF-MS of the toxins of the medically important African black-necked spitting cobra, Naja nigricollis (Tanzania); New Guinea small-eyed snake, Micropechis ikaheka; and Papuan black snake, Pseudechis papuanus. The main advantage of this approach is that only one generic sulfur-containing standard is required to quantify each and all intact Cys- and/or Met-containing toxins of the venom proteome. The results of absolute quantification are in reasonably good agreement with previously reported relative quantification of the most abundant protein families. However, both datasets depart in the quantification of the minor ones, showing a tendency for this set of proteins to be underestimated in standard peptide-centric venomics approaches. The molecular identity, specific toxic activity, and concentration in the venom, are the pillars on which the toxicovenomics-aimed discovery of the most medically-relevant venom toxins, e.g. those that need to be neutralized by an effective therapeutic antivenom, should be based. The pioneering venom proteome-wide absolute quantification shown in this paper represents thus a significant advance towards this goal. The potential of ICP triple quadrupole MS in proteomics in general, and venomics in particular, is critically discussed.

BIOLOGICAL SIGNIFICANCE

Animal venoms provide excellent model systems for investigating interactions between predators and prey, and the molecular mechanisms that contribute to adaptive protein evolution. On the other hand, numerous cases of snake bites occur yearly by encounters of humans and snakes in their shared natural environment. Snakebite envenoming is a serious global public health issue that affects the most impoverished and geopolitically disadvantaged rural communities in many tropical and subtropical countries. Unveiling the temporal and spatial patterns of venom variability is of fundamental importance to understand the molecular basis of envenoming, a prerequisite for developing therapeutic strategies against snakebite envenoming. Research on venoms has been continuously enhanced by advances in technology. The combined application of next-generation transcriptomic and venomic workflows has demonstrated unparalleled capabilities for venom characterization in unprecedented detail. However, mass spectrometry is not inherently quantitative, and this analytical limitation has sparked the development of methods to determine absolute abundance of proteins in biological samples. Here we show the potential of a hybrid element and molecular MS configuration for the parallel ESI-QToF-MS and ICP-QQQ detection and absolute quantification of intact sulfur-containing venom proteins via ³²S/³⁴S isotope dilution analysis. This configuration has been applied to quantify the toxins of the medically important African snake Naja nigricollis (Tanzania), and the Papuan species Micropechis ikaheka and Pseudechis papuanus.

Venom On-a-Chip: A Fast and Efficient Method for Comparative Venomics

Venom research has attracted an increasing interest in disparate fields, from drug development and pharmacology, to evolutionary biology and ecology, and rational antivenom production. Advances in “-omics” technologies have allowed the characterization of an increasing number of animal venoms, but the methodology currently available is suboptimal for large-scale comparisons of venom profiles. Here, we describe a fast, reproducible and semi-automated protocol for investigating snake venom variability, especially at the intraspecific level, using the Agilent Bioanalyzer on-chip technology. Our protocol generated a phenotype matrix which can be used for robust statistical analysis and correlations of venom variation with ecological correlates, or other extrinsic factors. We also demonstrate the ease and utility of combining on-chip technology with previously fractionated venoms for detection of specific individual toxin proteins. Our study describes a novel strategy for rapid venom discrimination and analysis of compositional variation at multiple taxonomic levels, allowing researchers to tackle evolutionary questions and unveiling the drivers of the incredible biodiversity of venoms.