Constructing comprehensive venom proteome reference maps for integrative venomics

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.

The Royal Armoury: Venomics and antivenomics of king cobra (Ophiophagus hannah) from the Indian Western Ghats

Despite being famous as 'the king' of the snake world, the king cobra (Ophiophagus hannah) has remained a mysterious species, particularly with respect to its venom ecology. In contrast, venom research has largely focussed on the 'big four' snakes that are greatly responsible for the burden of snakebite in the Indian subcontinent. This study aims to bridge the current void in our understanding of the O. hannah venom by investigating its proteomic, biochemical, pharmacological, and toxinological profiles via interdisciplinary approaches. Considering their physical resemblance, the king cobra is often compared to the spectacled cobra (Naja naja). Comparative venomics of O. hannah and N. naja in this study provided interesting insights into their venom compositions, activities, and potencies. Our findings suggest that the O. hannah venom, despite being relatively less complex than the N. naja venom, is equally potent. Finally, our in vitro and in vivo assays revealed that both Indian polyvalent and Thai Red Cross monovalent antivenoms completely fail to neutralise the O. hannah venom. Our findings provide guidelines for the management of bites from this clinically important yet neglected snake species in India.

One Size Fits All-Venomics of the Iberian Adder (Vipera seoanei, Lataste 1878) Reveals Low Levels of Venom Variation across Its Distributional Range

European vipers (genus Vipera) are medically important snakes displaying considerable venom variation, occurring at different levels in this group. The presence of intraspecific venom variation, however, remains understudied in several Vipera species. Vipera seoanei is a venomous snake endemic to the northern Iberian Peninsula and south-western France, presenting notable phenotypic variation and inhabiting several diverse habitats across its range. We analysed the venoms of 49 adult specimens of V. seoanei from 20 localities across the species' Iberian distribution. We used a pool of all individual venoms to generate a V. seoanei venom reference proteome, produced SDS-PAGE profiles of all venom samples, and visualised patterns of variation using NMDS. By applying linear regression, we then assessed presence and nature of venom variation between localities, and investigated the effect of 14 predictors (biological, eco-geographic, genetic) on its occurrence. The venom comprised at least 12 different toxin families, of which five (i.e., PLA2, svSP, DI, snaclec, svMP) accounted for about 75% of the whole proteome. The comparative analyses of the SDS-PAGE venom profiles showed them to be remarkably similar across the sampled localities, suggesting low geographic variability. The regression analyses suggested significant effects of biological and habitat predictors on the little variation we detected across the analysed V. seoanei venoms. Other factors were also significantly associated with the presence/absence of individual bands in the SDS-PAGE profiles. The low levels of venom variability we detected within V. seoanei might be the result of a recent population expansion, or of processes other than directional positive selection.

Interpopulational variation and ontogenetic shift in the venom composition of Lataste's viper (Vipera latastei, Boscá 1878) from northern Portugal

Lataste's viper (Vipera latastei) is a venomous European viper endemic to the Iberian Peninsula, recognised as medically important by the World Health Organization. To date, no comprehensive characterisation of this species' venom has been reported. Here, we analysed the venoms of juvenile and adult specimens of V. latastei from two environmentally different populations from northern Portugal. Using bottom-up venomics, we produced six venom proteomes (three per population) from vipers belonging to both age classes (i.e., two juveniles and four adults), and RP-HPLC profiles of 54 venoms collected from wild specimens. Venoms from juveniles and adults differed in their chromatographic profiles and relative abundances of their toxins, suggesting the occurrence of ontogenetic changes in venom composition. Specifically, snake venom metalloproteinase (SVMP) was the most abundant toxin family in juvenile venoms, while snake venom serine proteinases (SVSPs), phospholipases A₂ (PLA₂s), and C-type lectin-like (CTLs) proteins were the main toxins comprising adult venoms. The RP-HPLC venom profiles were found to vary significantly between the two sampled localities, indicating geographic variability. Furthermore, the presence/absence of certain peaks in the venom chromatographic profiles appeared to be significantly correlated also to factors like body size and sex of the vipers. Our findings show that V. latastei venom is a variable phenotype. The intraspecific differences we detected in its composition likely mirror changes in the feeding ecology of this species, taking place during different life stages and under different environmental pressures. SIGNIFICANCE: Lataste's viper (Vipera latastei) is a medically important viper endemic to the Iberian Peninsula, inhabiting different habitats and undergoing a marked ontogenetic dietary shift. In the current study, we report the first proteomic analysis of V. latastei venom from two environmentally different localities in northern Portugal. Our bottom-up venomic analyses show that snake venom serine proteinases (SVSPs), phospholipases A2 (PLA2s), and C-type lectin-like (CTLs) proteins are the major components of adult V. latastei venom. The comparative analysis of young and adult venoms suggests the occurrence of ontogenetic shift in toxin abundances, with snake venom metalloproteinases (SVMPs) being the predominant toxins in juvenile venoms. Moreover, geographic venom variation between the two studied populations is also detected, with our statistical analyses suggesting that factors like body size and sex of the vipers are possibly at play in its determination. Our work represents the first assessment of the composition of V. latastei venom, and the first step towards a better understanding of the drivers behind its variability.

Investigating Toxin Diversity and Abundance in Snake Venom Proteomes

Understanding snake venom proteomes is becoming increasingly important to understand snake venom biology, evolution and especially clinical effects of venoms and approaches to antivenom development. To explore the current state of snake venom proteomics and transcriptomics we investigated venom proteomic methods, associations between methodological and biological variability and the diversity and abundance of protein families. We reviewed available studies on snake venom proteomes from September 2017 to April 2021. This included 81 studies characterising venom proteomes of 79 snake species, providing data on relative toxin abundance for 70 species and toxin diversity (number of different toxins) for 37 species. Methodologies utilised in these studies were summarised and compared. Several comparative studies showed that preliminary decomplexation of crude venom by chromatography leads to increased protein identification, as does the use of transcriptomics. Combining different methodological strategies in venomic approaches appears to maximize proteome coverage. 48% of studies used the RP-HPLC →1D SDS-PAGE →in-gel trypsin digestion → ESI -LC-MS/MS pathway. Protein quantification by MS1-based spectral intensity was used twice as commonly as MS2-based spectral counting (33–15 studies). Total toxin diversity was 25–225 toxins/species, with a median of 48. The relative mean abundance of the four dominant protein families was for elapids; 3FTx–52%, PLA₂–27%, SVMP–2.8%, and SVSP–0.1%, and for vipers: 3FTx–0.5%, PLA₂–24%, SVMP–27%, and SVSP–12%. Viper venoms were compositionally more complex than elapid venoms in terms of number of protein families making up most of the venom, in contrast, elapid venoms were made up of fewer, but more toxin diverse, protein families. No relationship was observed between relative toxin diversity and abundance. For equivalent comparisons to be made between studies, there is a need to clarify the differences between methodological approaches and for acceptance of a standardised protein classification, nomenclature and reporting procedure. Correctly measuring and comparing toxin diversity and abundance is essential for understanding biological, clinical and evolutionary implications of snake venom composition.

Recombinant Production and Characterization of a New Toxin from Cryptops iheringi Centipede Venom Revealed by Proteome and Transcriptome Analysis

Among the Chilopoda class of centipede, the Cryptops genus is one of the most associated with envenomation in humans in the metropolitan region of the state of São Paulo. To date, there is no study in the literature about the toxins present in its venom. Thus, in this work, a transcriptomic characterization of the Cryptops iheringi venom gland, as well as a proteomic analysis of its venom, were performed to obtain a toxin profile of this species. These methods indicated that 57.9% of the sequences showed to be putative toxins unknown in public databases; among them, we pointed out a novel putative toxin named Cryptoxin-1. The recombinant form of this new toxin was able to promote edema in mice footpads with massive neutrophils infiltration, linking this toxin to envenomation symptoms observed in accidents with humans. Our findings may elucidate the role of this toxin in the venom, as well as the possibility to explore other proteins found in this work.

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.

Panacea within a Pandora's box: the antiparasitic effects of phospholipases A2 (PLA2s) from snake venoms

Parasitic diseases affect millions of individuals worldwide, mainly in low-income regions. There is no cure for most of these diseases, and the treatment relies on drugs that have side effects and lead to drug resistance, emphasizing the urgency to find new treatments. Snake venom has been gaining prominence as a rich source of molecules with antiparasitic potentials, such as phospholipases A₂ (PLA₂s). Here, we compile the findings involving PLA₂s with antiparasitic activities against helminths, Plasmodium, Toxoplasma, and trypanosomatids. We indicate their molecular features, highlighting the possible antiparasitic mechanisms of action of these proteins. We also demonstrate interactions between PLA₂s and some parasite membrane components, shedding light on potential targets for drug design that may provide better treatment for the illnesses caused by parasites.

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.

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.

Antivenomics and in vivo preclinical efficacy of six Latin American antivenoms towards south-western Colombian Bothrops asper lineage venoms

Background

Bothrops asper represents the clinically most important snake species in Central America and Northern South America, where it is responsible for an estimated 50–80% of snakebites. Compositional variability among the venom proteomes of B. asper lineages across its wide range mirrors clinical differences in their envenomings. Bothropic antivenoms generated in a number of Latin American countries commonly exhibit a certain degree of paraspecific effectiveness in the neutralization of congeneric venoms. Defining the phylogeographic boundaries of an antivenom's effectivity has implications for optimizing its clinical use. However, the molecular bases and impact of venom compositions on the immune recognition and neutralization of the toxic activities of across geographically disparate populations of B. asper lineages has not been comprehensively studied.

Methodology/Principal findings

Third-generation antivenomics was applied to quantify the cross-immunorecognizing capacity against the individual components of venoms of three B. asper lineages (B. asper (sensu stricto), B. ayerbei and B. rhombeatus) distributed in south-western (SW) Colombia, of six Latin American antivenoms, produced against homologous (Colombia, INS-COL and PROBIOL) and Costa Rica (ICP)), and heterologous (Argentina (BIOL), Perú (INS-PERU) and Venezuela (UCV)) bothropic venoms. In vivo neutralization assays of the lethal, hemorrhagic, coagulant, defibrinogenating, myotoxic, edematogenic, indirect hemolytic, and proteolytic activities of the three SW Colombian B. asper lineage venoms were carried to compare the preclinical efficacy of three (Colombian INS-COL and PROBIOL, and Costa Rican ICP) antivenoms frequently used in Colombia. Antivenomics showed that all the six antivenom affinity matrices efficiently immunoretained most of the B. asper lineages venom proteins and exhibited impaired binding towards the venoms' peptidomes. The neutralization profile of the INS-COL, PROBIOL and ICP antivenoms towards the biological activities of the venoms of SW Colombian B. asper (sensu stricto), B. ayerbei and B. rhombeatus lineages was coherent with the antivenomics outcome. In addition, the combination of in vitro (antivenomics) and in vivo neutralization results allowed us to determine their toxin-specific and venom neutralizing antibody content. Noteworthy, heterologous INS-PERU, BIOL, and UCV bothropic antivenoms had equal or higher binding capacity towards the venoms components of SW Colombian B. asper lineages that the homologous Colombian and Costa Rican antivenoms.

Conclusions/Significance

The combined in vitro and in vivo preclinical outcome showed that antivenoms manufactured in Colombia and Costa Rica effectively neutralize the major toxic activities of SW Colombian B. asper lineage venoms. The antivenomics profiles of the heterologous antivenoms manufactured in Argentina, Venezuela, and Perú strongly suggests their (pre)clinical adequacy for the treatment of B. asper lineage envenomings in SW Colombia. However, their recommendation in the clinical setting is pending on in vivo neutralization testing and clinical testing in humans. Bothrops asper is a highly adaptable snake species complex, which is considered the most dangerous snake throughout much of its distribution range from the Atlantic lowland of eastern México to northwestern Perú. Antivenoms are the only scientifically validated treatment of snakebite envenomings. Venom variation is particularly common in wide ranging species, such as B. asper, and may result in variable clinical presentations of envenomings, as is the case for the B. asper species complex, potentially undermining the efficacy of snakebite treatments depending on the immunization mixture used in the generation of the antivenom. Conversely, phylogenetic conservation of antigenic determinants confers an unpredictable degree of paraspecificity to homologous antivenoms produced for a geographic area, but also to heterologous congeneric antivenoms, towards the venom components of allopatric conspecific populations. This work aimed at comparing the preclinical profile of a panel of Latin American homologous and heterologous antivenoms against the venoms of B. asper lineages distributed in SW Colombia. The outcome of this study strongly suggests the suitability of considering the heterologous antivenoms BIOL (Argentina), UCV (Venezuela) and INS-PERU (Perú) as alternatives to homologous Colombian INS-COL and PROBIOL and Costa Rican ICP antivenoms for the treatment of envenomings by B. asper (sensu stricto) in W Colombia and Ecuador, B. ayerbei in Cauca and Nariño (Colombia), and B. rhombeatus in Cauca river valley, SW Colombia.

Snakebite envenoming is an important occupational health problem, particularly in rural areas of developing countries. The timely administration of an effective antivenom remains the mainstay of snakebite management. However, the use of antivenoms is often limited by non-availability due to high cost or by lack of effectiveness. Antivenom shortage can be addressed through the generation of novel polyspecific antivenoms of wide clinical efficacy against the venoms of the medically-relevant snake species within the geographical range where these antivenoms are intended to be deployed, but also by optimizing the paraspecific use of current antivenoms. In Colombia, antivenoms are supplied by two manufacturers, one public, the Instituto Nacional de Salud (INS), and one private, Laboratorios Probiol (PROBIOL). However, the antivenom supply in Colombia has traditionally been insufficient, a circumstance that has led the Colombian Ministerio de Salud y Protección Social to issue several resolutions and decrees to announce this health emergency in the country, and to import antivenoms produced in México and Costa Rica. Contrary to these countries, where B. asper represents the only species of the genus, in SW Colombia three close phylogenetically related B. asper lineages, B. asper (sensu stricto), B. rhombeatus, and B. ayerbei, are responsible for most severe cases of snakebite accidents and exhibit remarkable differences in the physiopathological profile of their envenomings. This work aimed to assess the immunorecognition characteristics of a panel of antivenoms manufactured in Colombia, Costa Rica, Argentina, Perú and Venezuela towards the venoms of the three SW Colombian B. asper lineages. Additionally, combined quantitative in vitro and in vivo data show that the homologous antivenoms produced in Colombia (INS-COL, PROBIOL) and Costa Rica (ICP) effectively neutralize the lethality and the major toxic activities tested of the three SW Colombian B. asper lineage venoms. Heterologous Argentinian (BIOL), Venezuelan (UCV) and Peruvian (INS-PERU) antivenoms also showed comparable, even higher, effective immunocapturing ability towards the venom proteomes of SW Colombian B. asper (sensu stricto), B. rhombeatus, and B. ayerbei, than the Colombian and Costa Rican antivenoms. These results are in line with previous studies highlighting the notable conservation of paraspecific antigenic determinants across the phylogeny of genus Bothrops, and advocate for considering the heterologous Argentinian, Venezuelan and Peruvian antivenoms as further therapeutic alternatives for the treatment of B. asper spp. snakebites in Colombia.

Deadly Proteomes: A Practical Guide to Proteotranscriptomics of Animal Venoms

Animal venoms are renowned for their toxicity, biochemical complexity, and as a source of compounds with potential applications in medicine, agriculture, and industry. Polypeptides underlie much of the pharmacology of animal venoms, and elucidating these arsenals of polypeptide toxins-known as the venom proteome or venome-is an important step in venom research. Proteomics is used for the identification of venom toxins, determination of their primary structure including post-translational modifications, as well as investigations into the physiology underlying their production and delivery. Advances in proteomics and adjacent technologies has led to a recent upsurge in publications reporting venom proteomes. Improved mass spectrometers, better proteomic workflows, and the integration of next-generation sequencing of venom-gland transcriptomes and venomous animal genomes allow quicker and more accurate profiling of venom proteomes with greatly reduced starting material. Technologies such as imaging mass spectrometry are revealing additional insights into the mechanism, location, and kinetics of venom toxin production. However, these numerous new developments may be overwhelming for researchers designing venom proteome studies. Here, the field of venom proteomics is reviewed and some practical solutions for simplifying mass spectrometry workflows to study animal venoms are offered.

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 proteomics of geographically distinct saw-scaled viper (Echis carinatus) venoms from India

Snakebite is a socio-economic problem in tropical countries and it is exacerbated by geographical venom variation of snakes. We investigated on venom variation in geographically distinct populations of Echis carinatus from three ecologically distinct regions: Tamil Nadu (ECVTN), Goa (ECVGO), and Rajasthan (ECVRAJ). Venom was fractionated by RP-HPLC, combined with SDS-PAGE, and subjected to tandem mass spectrometry. Toxins were identified, and their relative abundance was estimated. Using NCBI database of Echis genus, we queried the MS/MS spectra, and found 69, 38 and 38 proteins in ECVTN, ECVGO and ECVRAJ respectively, belonging to 8-10 different toxin families. The differences in the venom profiles were due to change in the relative composition of the toxin families. Snake venom metalloproteinase (svMP), Snaclecs and Phospholipase A2 (PLA2) were the major venom components in all the venoms. Heteromeric Disintegrins were found in ECVTN and absent in other venoms. ECVRAJ showed higher abundance of low-molecular-weight (>30 kDa) proteins than ECVTN and ECVGO. Cysteine-rich venom protein (CRISP) was highest in ECVRAJ (7.34%), followed by ECVTN (0.01%) and in ECVGO, it was not detected. These findings highlight the need for evaluating the efficacy of the polyvalent anti-venom to neutralize the toxins from geographically distinct venoms of E. carinatus.

Bothrops atrox, the most important snake involved in human envenomings in the amazon: How venomics contributes to the knowledge of snake biology and clinical toxinology

Bothrops atrox snakes are mostly endemic of the Amazon rainforest and is certainly the South American pit viper responsible for most of the snakebites in the region. The composition of B. atrox venom is significantly known and has been used to trace the relevance of the venom phenotype for snake biology and for the impacts in the clinics of human patients involved in accidents by B. atrox. However, in spite of the wide distribution and the great medical relevance of B. atrox snakes, B. atrox taxonomy is not fully resolved and the impacts of the lack of taxonomic resolution on the studies focused on venom or envenoming are currently unknown. B. atrox venom presents different degrees of compositional variability and is generally coagulotoxic, inducing systemic hematological disturbances and local tissue damage in snakebite patients. Antivenoms are the effective therapy for attenuating the clinical signs. This review brings a comprehensive discussion of the literature concerning B. atrox snakes encompassing from snake taxonomy, diet and venom composition, towards clinical aspects of snakebite patients and efficacy of the antivenoms. This discussion is highly supported by the contributions that venomics and antivenomics added for the advancement of knowledge of B. atrox snakes, their venoms and the treatment of accidents they evoke.

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 Zn2+-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.

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.

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.

Vipera berus berus Venom from Russia: Venomics, Bioactivities and Preclinical Assessment of Microgen Antivenom

The common European adder, Vipera berus berus, is a medically relevant species, which is widely distributed in Russia and thus, is responsible for most snakebite accidents in Russia. We have investigated the toxic and enzymatic activities and have determined the proteomic composition of its venom. Phospholipases A₂ (PLA₂, 25.3% of the venom proteome), serine proteinases (SVSP, 16.2%), metalloproteinases (SVMP, 17.2%), vasoactive peptides (bradykinin-potentiating peptides (BPPs), 9.5% and C-type natriuretic peptides (C-NAP, 7.8%), cysteine-rich secretory protein (CRISP, 8%) and L-amino acid oxidase (LAO, 7.3%) represent the major toxin classes found in V. b. berus (Russia) venom. This study was also designed to assess the in vivo and in vitro preclinical efficacy of the Russian Microgen antivenom in neutralizing the main effects of V. b. berus venom. The results show that this antivenom is capable of neutralizing the lethal, hemorrhagic and PLA₂ activities. Third-generation antivenomics was applied to quantify the toxin-recognition landscape and the maximal binding capacity of the antivenom for each component of the venom. The antivenomics analysis revealed that 6.24% of the anti-V. b. berus F(ab’)₂ molecules fraction are toxin-binding antibodies, 60% of which represent clinically relevant antivenom molecules.

Defining the pathogenic threat of envenoming by South African shield-nosed and coral snakes (genus Aspidelaps), and revealing the likely efficacy of available antivenom

While envenoming by the southern African shield-nosed or coral snakes (genus Aspidelaps) has caused fatalities, bites are uncommon. Consequently, this venom is not used in the mixture of snake venoms used to immunise horses for the manufacture of regional SAIMR (South African Institute for Medical Research) polyvalent antivenom. Aspidelaps species are even excluded from the manufacturer's list of venomous snakes that can be treated by this highly effective product. This leaves clinicians, albeit rarely, in a therapeutic vacuum when treating envenoming by these snakes. This is a significantly understudied small group of nocturnal snakes and little is known about their venom compositions and toxicities. Using a murine preclinical model, this study determined that the paralysing toxicity of venoms from Aspidelaps scutatus intermedius, A. lubricus cowlesi and A. l. lubricus approached that of venoms from highly neurotoxic African cobras and mambas. This finding was consistent with the cross-genus dominance of venom three-finger toxins, including numerous isoforms which showed extensive interspecific variation. Our comprehensive analysis of venom proteomes showed that the three Aspidelaps species possess highly similar venom proteomic compositions. We also revealed that the SAIMR polyvalent antivenom cross-reacted extensively in vitro with venom proteins of the three Aspidelaps. Importantly, this cross-genus venom-IgG binding translated to preclinical (in a murine model) neutralisation of A. s. intermedius venom-induced lethality by the SAIMR polyvalent antivenom, at doses comparable with those that neutralise venom from the cape cobra (Naja nivea), which the antivenom is directed against. Our results suggest a wider than anticipated clinical utility of the SAIMR polyvalent antivenom, and here we seek to inform southern African clinicians that this readily available antivenom is likely to prove effective for victims of Aspidelaps envenoming. BIOLOGICAL SIGNIFICANCE: Coral and shield-nosed snakes (genus Aspidelaps) comprise two species and several subspecies of potentially medically important venomous snakes distributed in Namibia, Botswana, Zimbabwe, Mozambique and South Africa. Documented human fatalities, although rare, have occurred from both A. lubricus and A. scutatus. However, their venom proteomes and the pathological effects of envenomings by this understudied group of nocturnal snakes remain uncharacterised. Furthermore, no commercial antivenom is made using venom from species of the genus Aspidelaps. To fill this gap, we have conducted a transcriptomics-guided comparative proteomics analysis of the venoms of the intermediate shield-nose snake (A. s. intermedius), southern coral snake (A. l. lubricus), and Cowle's shield snake (A. l. cowlesi); investigated the mechanism of action underpinning lethality by A. s. intermedius in the murine model; and assessed the in vitro immunoreactivity of the SAIMR polyvalent antivenom towards the venom toxins of A. l. lubricus and A. l. cowlesi, and the in vivo capability of this antivenom at neutralising the lethal effect of A. s. intermedius venom. Our data revealed a high degree of conservation of the global composition of the three Aspidelaps venom proteomes, all characterised by the overwhelming predominance of neurotoxic 3FTxs, which induced classical signs of systemic neurotoxicity in mice. The SAIMR polyvalent antivenom extensively binds to Aspidelaps venom toxins and neutralised, with a potency of 0.235 mg venom/mL antivenom, the lethal effect of A. s. intermedius venom. Our data suggest that the SAIMR antivenom could be a useful therapeutic tool for treating human envenomings by Aspidelaps species.

Translational Venomics: Third-Generation Antivenomics of Anti-Siamese Russell's Viper, Daboia siamensis, Antivenom Manufactured in Taiwan CDC's Vaccine Center

The venom proteome of Siamese Russell’s viper from Taiwan, alongside complementary in vivo lethality neutralization assay and in vitro third-generation antivenomics assessment of the preclinical efficacy of the homologous antivenom manufactured in Taiwan CDC’s Vaccine Center, are here reported. Taiwanese Russell’s viper venom proteome comprised 25 distinct gene products, with the heterodimeric PLA₂ viperotoxin-F representing the most abundant toxin (47.5% of total venom proteome). Coagulation FV-activating serine proteinase (RVV-V, 14%), the PIV-SVMP activator of FX (RVV-FX, 8.5%), and less abundant toxins from nine protein families, make up its venom proteome. Venom composition-pathology correlations of D. siamensis envenomings in Taiwan are discussed. The lethal effect of Taiwanese D. siamensis venom was 0.47 mg/g mouse. Antivenomics-guided assessment of the toxin recognition landscape of the Taiwanese Russell’s viper antivenom, in conjunction with complementary in vivo neutralization analysis, informed the antivenom’s maximal toxin immunorecognition ability (14 mg total venom proteins/vial), neutralization capacity (6.5 mg venom/vial), and relative content of lethality neutralizing antibodies (46.5% of the toxin-binding F(ab’)₂ antibodies). The antivenomics analysis also revealed suboptimal aspects of the CDC-Taiwan antivenom. Strategies to improve them are suggested.

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.

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.

Proteomics and antivenomics of Echis carinatus carinatus venom: Correlation with pharmacological properties and pathophysiology of envenomation

The proteome composition of Echis carinatus carinatus venom (ECV) from India was studied for the first time by tandem mass spectrometry analysis. A total of 90, 47, and 22 distinct enzymatic and non-enzymatic proteins belonging to 15, 10, and 6 snake venom protein families were identified in ECV by searching the ESI-LC-MS/MS data against non-redundant protein databases of Viperidae (taxid 8689), Echis (taxid 8699) and Echis carinatus (taxid 40353), respectively. However, analysis of MS/MS data against the Transcriptome Shotgun Assembly sequences (87 entries) of conger E. coloratus identified only 14 proteins in ECV. Snake venom metalloproteases and snaclecs, the most abundant enzymatic and non-enzymatic proteins, respectively in ECV account for defibrinogenation and the strong in vitro pro-coagulant activity. Further, glutaminyl cyclase, aspartic protease, aminopeptidase, phospholipase B, vascular endothelial growth factor, and nerve growth factor were reported for the first time in ECV. The proteome composition of ECV was well correlated with its biochemical and pharmacological properties and clinical manifestations observed in Echis envenomed patients. Neutralization of enzymes and pharmacological properties of ECV, and immuno-cross-reactivity studies unequivocally point to the poor recognition of <20 kDa ECV proteins, such as PLA₂, subunits of snaclec, and disintegrin by commercial polyvalent antivenom.

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.

Rapid screening and identification of ACE inhibitors in snake venoms using at-line nanofractionation LC-MS

This study presents an analytical method for the screening of snake venoms for inhibitors of the angiotensin-converting enzyme (ACE) and a strategy for their rapid identification. The method is based on an at-line nanofractionation approach, which combines liquid chromatography (LC), mass spectrometry (MS), and pharmacology in one platform. After initial LC separation of a crude venom, a post-column flow split is introduced enabling parallel MS identification and high-resolution fractionation onto 384-well plates. The plates are subsequently freeze-dried and used in a fluorescence-based ACE activity assay to determine the ability of the nanofractions to inhibit ACE activity. Once the bioactive wells are identified, the parallel MS data reveals the masses corresponding to the activities found. Narrowing down of possible bioactive candidates is provided by comparison of bioactivity profiles after reversed-phase liquid chromatography (RPLC) and after hydrophilic interaction chromatography (HILIC) of a crude venom. Additional nanoLC-MS/MS analysis is performed on the content of the bioactive nanofractions to determine peptide sequences. The method described was optimized, evaluated, and successfully applied for screening of 30 snake venoms for the presence of ACE inhibitors. As a result, two new bioactive peptides were identified: pELWPRPHVPP in Crotalus viridis viridis venom with IC₅₀ = 1.1 μM and pEWPPWPPRPPIPP in Cerastes cerastes cerastes venom with IC₅₀ = 3.5 μM. The identified peptides possess a high sequence similarity to other bradykinin-potentiating peptides (BPPs), which are known ACE inhibitors found in snake venoms.

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.

Combined venom profiling and cytotoxicity screening of the Radde's mountain viper (Montivipera raddei) and Mount Bulgar Viper (Montivipera bulgardaghica) with potent cytotoxicity against human A549 lung carcinoma cells

Here we report the first characterization of the endemic Mount Bulgar Viper (Montivipera bulgardaghica) and Radde's mountain viper (Montivipera raddei) venom by a combined approach using intact mass profiling and bottom-up proteomics. The cytotoxicity screening of crude venom as well as isolated serine proteases revealed a high activity against A549 human lung carcinoma cells. By means of intact mass profiling of native and reduced venom we observed basic and acidic phospholipases type A₂. Moreover, the analysis revealed snake venom metalloproteases, cysteine-rich secretory proteins, disintegrins, snake venom serine proteases, C-type lectins, a vascular endothelial growth factor and an L-amino acid oxidase.

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.

Preclinical Evaluation of the Efficacy of Antivenoms for Snakebite Envenoming: State-of-the-Art and Challenges Ahead

Animal-derived antivenoms constitute the mainstay in the therapy of snakebite envenoming. The efficacy of antivenoms to neutralize toxicity of medically-relevant snake venoms has to be demonstrated through meticulous preclinical testing before their introduction into the clinical setting. The gold standard in the preclinical assessment and quality control of antivenoms is the neutralization of venom-induced lethality. In addition, depending on the pathophysiological profile of snake venoms, the neutralization of other toxic activities has to be evaluated, such as hemorrhagic, myotoxic, edema-forming, dermonecrotic, in vitro coagulant, and defibrinogenating effects. There is a need to develop laboratory assays to evaluate neutralization of other relevant venom activities. The concept of the 3Rs (Replacement, Reduction, and Refinement) in Toxinology is of utmost importance, and some advances have been performed in their implementation. A significant leap forward in the study of the immunological reactivity of antivenoms against venoms has been the development of “antivenomics”, which brings the analytical power of mass spectrometry to the evaluation of antivenoms. International partnerships are required to assess the preclinical efficacy of antivenoms against snake venoms in different regions of the world in order to have a detailed knowledge on the neutralizing profile of these immunotherapeutics.

Protein-species quantitative venomics: looking through a crystal ball

In this paper we discuss recent significant developments in the field of venom research, specifically the emergence of top-down proteomic applications that allow achieving compositional resolution at the level of the protein species present in the venom, and the absolute quantification of the venom proteins (the term “protein species” is used here to refer to all the different molecular forms in which a protein can be found. Please consult the special issue of Jornal of Proteomics “Towards deciphering proteomes via the proteoform, protein speciation, moonlighting and protein code concepts” published in 2016, vol. 134, pages 1-202). Challenges remain to be solved in order to achieve a compact and automated platform with which to routinely carry out comprehensive quantitative analysis of all toxins present in a venom. This short essay reflects the authors’ view of the immediate future in this direction for the proteomic analysis of venoms, particularly of snakes.

Strategies in 'snake venomics' aiming at an integrative view of compositional, functional, and immunological characteristics of venoms

This work offers a general overview on the evolving strategies for the proteomic analysis of snake venoms, and discusses how these may be combined through diverse experimental approaches with the goal of achieving a more comprehensive knowledge on the compositional, toxic, and immunological characteristics of venoms. Some recent developments in this field are summarized, highlighting how strategies have evolved from the mere cataloguing of venom components (proteomics/venomics), to a broader exploration of their immunological (antivenomics) and functional (toxicovenomics) characteristics. Altogether, the combination of these complementary strategies is helping to build a wider, more integrative view of the life-threatening protein cocktails produced by venomous snakes, responsible for thousands of deaths every year.

The genetics of venom ontogeny in the eastern diamondback rattlesnake (Crotalus adamanteus)

The same selective forces that give rise to rapid inter- and intraspecific divergence in snake venoms can also favor differences in venoms across life-history stages. Ontogenetic changes in venom composition are well known and widespread in snakes but have not been investigated to the level of unambiguously identifying the specific loci involved. The eastern diamondback rattlesnake was previously shown to undergo an ontogenetic shift in venom composition at sexual maturity, and this shift accounted for more venom variation than geography. To characterize the genetics underlying the ontogenetic venom compositional change in C. adamanteus, we sequenced adult/juvenile pairs of venom-gland transcriptomes from five populations previously shown to have different adult venom compositions. We identified a total of 59 putative toxin transcripts for C. adamanteus, and 12 of these were involved in the ontogenetic change. Three toxins were downregulated, and nine were upregulated in adults relative to juveniles. Adults and juveniles expressed similar total levels of snake-venom metalloproteinases but differed substantially in their featured paralogs, and adults expressed higher levels of Bradykinin-potentiating and C-type natriuretic peptides, nerve growth factor, and specific paralogs of phospholipases A₂ and snake venom serine proteinases. Juvenile venom was more toxic to mice, indicating that the expression differences resulted in a phenotypically, and therefore potentially ecologically, significant difference in venom function. We also showed that adult and juvenile venom-gland transcriptomes for a species with known ontogenetic venom variation were equally effective at individually providing a full characterization of the venom genes of a species but that any particular individual was likely to lack several toxins in their transcriptome. A full characterization of a species’ venom-gene complement therefore requires sequencing more than one individual, although the ages of the individuals are unimportant.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

What killed Karl Patterson Schmidt? Combined venom gland transcriptomic, venomic and antivenomic analysis of the South African green tree snake (the boomslang), Dispholidus typus

Background

Non-front-fanged colubroid snakes comprise about two-thirds of extant ophidian species. The medical significance of the majority of these snakes is unknown, but at least five species have caused life-threatening or fatal human envenomings. However, the venoms of only a small number of species have been explored.

Methods

A combined venomic and venom gland transcriptomic approach was employed to characterise of venom of Dispholidus typus (boomslang), the snake that caused the tragic death of Professor Karl Patterson Schmidt. The ability of CroFab™ antivenom to immunocapture boomslang venom proteins was investigated using antivenomics.

Results

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families: three-finger toxin (3FTx); phospholipase A₂ (PLA₂); cysteine-rich secretory proteins (CRISP); snake venom (SV) serine proteinase (SP); C-type lectin-like (CTL); SV metalloproteinases (SVMPs); and disintegrin-like/cysteine-rich (DC) proteolytic fragments. CroFab™ antivenom efficiently immunodepleted some boomslang SVMPs.

Conclusions

The present work is the first to address the overall proteomic profile of D. typus venom. This study allowed us to correlate the toxin composition with the toxic activities of the venom. The antivenomic analysis suggested that the antivenom available at the time of the unfortunate accident could have exhibited at least some immunoreactivity against the boomslang SVMPs responsible for the disseminated intravascular coagulation syndrome that caused K.P. Schmidt's fatal outcome.

General significance

This study may stimulate further research on other non-front-fanged colubroid snake venoms capable of causing life-threatening envenomings to humans, which in turn should contribute to prevent fatal human accidents, such as that unfortunately suffered by K.P. Schmidt.

•

The venom proteome of Dispholidus typus (boomslang) is reported.

•

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families.

•

Boomslang venom proteome is dominated (75%) by snake venom PIII-metalloproteinases (PIII-SVMPs).

•

CroFab™ antivenom efficiently immunodepleted some boomslang PIII-SVMPs.

Small Packages, Big Returns: Uncovering the Venom Diversity of Small Invertebrate Conoidean Snails

Venomous organisms used in research were historically chosen based on size and availability. This opportunity-driven strategy created a species bias in which snakes, scorpions, and spiders became the primary subjects of venom research. Increasing technological advancements have enabled interdisciplinary studies using genomics, transcriptomics, and proteomics to expand venom investigation to animals that produce small amounts of venom or lack traditional venom producing organs. One group of non-traditional venomous organisms that have benefitted from the rise of -omic technologies is the Conoideans. The Conoidean superfamily of venomous marine snails includes, the Terebridae, Turridae (s.l), and Conidae. Conoidea venom is used for both predation and defense, and therefore under strong selection pressures. The need for conoidean venom peptides to be potent and specific to their molecular targets has made them important tools for investigating cellular physiology and bioactive compounds that are beneficial to improving human health. A convincing case for the potential of Conoidean venom is made with the first commercially available conoidean venom peptide drug Ziconotide (Prialt®), an analgesic derived from Conus magus venom that is used to treat chronic pain in HIV and cancer patients. Investigation of conoidean venom using -omics technology provides significant insights into predator-driven diversification in biodiversity and identifies novel compounds for manipulating cellular communication, especially as it pertains to disease and disorders.