Venom Variation Research Articles

High throughput identification of human monoclonal antibodies and heavy-chain-only antibodies to treat snakebite

Snakebite envenoming is a priority Neglected Tropical Disease that causes an estimated 81,000–135,000 fatalities each year. The development of a new generation of safer, affordable, and accessible antivenom therapies is urgently needed. With this goal in mind, rigorous characterisation of the specific toxins in snake venom is key to generating novel therapies for snakebite. Monoclonal antibodies directed against venom toxins are emerging as potentially strong candidates in the development of new snakebite diagnostics and treatment. Venoms comprise many different toxins of which several are responsible for their pathological effects. Due to the large variability of venoms within and between species, formulations of combinations of human antibodies are proposed as the next generation antivenoms. Here a high-throughput screening method employing antibody-based ligand fishing of venom toxins in 384 filter-well plate format has been developed to determine the antibody target/s The approach uses Protein G beads for antibody capture followed by exposure to a full venom or purified toxins to bind their respective ligand toxin(s). This is followed by a washing/centrifugation step to remove non-binding toxins and an in-well tryptic digest. Finally, peptides from each well are analysed by nanoLC-MS/MS and subsequent Mascot database searching to identify the bound toxin/s for each antibody under investigation. The approach was successfully validated to rapidly screen antibodies sourced from hybridomas, derived from venom-immunised mice expressing either regular human antibodies or heavy-chain-only human antibodies (HCAbs).

Predatory and Defensive Strategies in Cone Snails.

Cone snails are carnivorous marine animals that prey on fish (piscivorous), worms (vermivorous), or other mollusks (molluscivorous). They produce a complex venom mostly made of disulfide-rich conotoxins and conopeptides in a compartmentalized venom gland. The pharmacology of cone snail venom has been increasingly investigated over more than half a century. The rising interest in cone snails was initiated by the surprising high human lethality rate caused by the defensive stings of some species. Although a vast amount of information has been uncovered on their venom composition, pharmacological targets, and mode of action of conotoxins, the venom-ecology relationships are still poorly understood for many lineages. This is especially important given the relatively recent discovery that some species can use different venoms to achieve rapid prey capture and efficient deterrence of aggressors. Indeed, via an unknown mechanism, only a selected subset of conotoxins is injected depending on the intended purpose. Some of these remarkable venom variations have been characterized, often using a combination of mass spectrometry and transcriptomic methods. In this review, we present the current knowledge on such specific predatory and defensive venoms gathered from sixteen different cone snail species that belong to eight subgenera: Pionoconus, Chelyconus, Gastridium, Cylinder, Conus, Stephanoconus, Rhizoconus, and Vituliconus. Further studies are needed to help close the gap in our understanding of the evolved ecological roles of many cone snail venom peptides.

Investigating snake venom variation to mitigate snakebite envenomation in Indonesia

Snakebite envenoming is one of the neglected tropical diseases and is still rarely studied in Indonesia. The high diversity of venomous snakes in Indonesia is one of the challenges. To overcome this, we consider that snake venomic research guided by phylogenetic relationships can serve valuable information that may contribute to snakebite mitigation. Here, we briefly introduce our recent study using HPLC-MS/MS to analyse the venom composition and variation across four Indonesian Trimeresurus. This adds more information on venom variation among Trimeresurus species within the close geographic origin, which congruent with their phylogenetic relationships. We also conducted a preliminary study to detect intraspecific variation between the venom of males and females of T. puniceus using SDS-PAGE. There is an addition of protein with a molecular weight of 13 kDa in the venom of male T. puniceus, while a protein with a molecular weight of 16 kDa is only detected in female venom. We summarize recent studies showing different factors that can affect venom variation between and within snake species. Finally, we discuss the importance of transdisciplinary research to understand snake venom variation and suggest future directions, particularly from a herpetological view, to mitigate human-snake conflict in Indonesia.

Analysis of the Snake Venom Peptidome.

Snake venom peptidomes are known to be a large source of molecules with different pharmacological properties. The complexity and variability of snake venoms, the presence of proteinases, and the lack of complete species-specific genome sequences make snake venom peptidome profiling a challenging task that requires especial technical strategies for sample processing and mass spectrometric analysis. Here, we describe a method for assessing the content of snake venom peptides and highlight the importance of sampling procedures, as they substantially influence the peptidomic complexity of snake venoms.

Immunoreactivity and neutralization efficacy of Pakistani Viper Antivenom (PVAV) against venoms of Saw-scaled Vipers (Echis carinatus subspp.) and Western Russell's Vipers (Daboia russelii) from the Indian subcontinent

Snakebite envenoming (SBE) is a priority Neglected Tropical Disease listed by the World Health Organization. South Asia is heavily affected, and virtually all countries in the region import polyvalent antivenom products from India for clinical use. The imported antivenoms, however, have suboptimal effectiveness due to geographical venom variation. Recently, a domestic bivalent product, named Pakistani Viper Antivenom (PVAV) has been developed specifically for Pakistani vipers, Echis carinatus sochureki and Daboia russelii. As a bivalent viperid antivenom, it is unknown yet if PVAV exhibits higher immunological binding and neutralization activities against viper venoms from distant locales compared with polyvalent antivenoms manufactured in India. This study thus examined the preclinical efficacy of PVAV against venoms of Western Russell's Vipers and Saw-scaled Viper subspecies from selected locales in the Indian subcontinent. PVAV generally outperformed the commonly used VINS polyvalent antivenom (VPAV, manufactured in India) in binding toward venoms, and showed superior or comparable neutralization efficacy against the venom procoagulant and hemorrhagic effects of Saw-scaled Vipers as well as Russell's Vipers from Pakistan and Sri Lanka. Based on normalized potency values, PVAV is far more potent than VPAV in neutralizing the lethality of all viper venoms, except that of the Indian Russell's Viper. The study shows conserved antigenicity of toxins responsible for major toxicity across these viperid venoms, and suggests the feasible production of a viper-specific antivenom with higher potency and broader geographical utility for the region.

Venomics of Leiurus abdullahbayrami, the most lethal scorpion in the Levant region of the Middle East

The scorpion Leiurus abdullahbayrami has been associated with severe/lethal envenomings throughout the Levant region of the Middle East, encompassing Turkey, Syria, and Lebanon, and only scarce information is available on its venom composition, activity, and antigenicity. We report on the composition of L. abdullahbayrami venom collected from Lebanese specimens using nESI-MS/MS, MALDI-TOF MS, SDS-PAGE and RP-HPLC. Venom lethality, through LD50 determination in mice (intraperitoneal), was also assessed (0.75 (0.16–1.09) mg/kg), confirming L.abdullahbayrami venom vertebrate toxicity. Fifty-four peaks were detected using RP-HPLC, half of which eluted in the gradient region between 20 and 40% acetonitrile. In reducing SDS-PAGE, most predominant components were <10 kDa, with minor components at higher molecular masses of 24.4, 43.1, and 48.9 kDa. Venom mass fingerprint by MALDI-TOF detected 21 components within the 1000–12,000 m/z range. Whole venom ‘shotgun’ bottom-up nLC-MS/MS approach, combined with in-gel tryptic digestion of SDS-PAGE bands, identified at least 113 different components belonging to 15 venom families and uncharacterized proteins, with ion channel-active components (K+ channel toxins (28); Na+ channel toxins (42); Cl− channel toxins (4); Ca+2 toxins (2)) being predominant. A single match for a L. adbullahbayrami NaTx was found in the UniProt database with other congeneric species, toxin h3.1 from Leiurus hebraeus, suggesting this might be an indication of venom divergence within Leiurus, eventhough this warrants further investigation involving venom proteomics and transcriptomics of relevant species. Considering such potential interspecific venom variation, future work should address whether preparation of a specific anti-L. abdullahbayrami antivenom is justified.

Varying Modes of Selection Among Toxin Families in the Venoms of the Giant Desert Hairy Scorpions (Hadrurus).

Venoms are primarily believed to evolve under strong diversifying selection resulting from persistent coevolution between predator and prey. Recent research has challenged this hypothesis, proposing that venoms from younger venomous lineages (e.g., snakes and cone snails) are governed predominantly by diversifying selection, while venoms from older venomous lineages (e.g., centipedes, scorpions, and spiders) are under stronger purifying selection. However, most research in older lineages has tested selection at more diverse phylogenetic scales. Although these tests are important for evaluating broad macroevolutionary trends underlying venom evolution, they are less equipped to detect species-level evolutionary trends, which likely have large impacts on venom variation seen at more diverse phylogenetic scales. To test for selection among closely related species from an older venomous lineage, we generated high-throughput venom-gland transcriptomes and venom proteomes for four populations of Giant Desert Hairy Scorpions (Hadrurus), including three Hadrurus arizonensis populations and one Hadrurus spadix population. We detected significant episodic and pervasive diversifying selection across a highly abundant toxin family that likely has a major role in venom function ([Formula: see text]KTxs), providing a contrast to the stronger purifying selection identified from other studies on scorpion venoms. Conversely, we detected weak episodic diversifying and/or stronger purifying selection in four toxin families (non-disulfide bridged peptides, phospholipase A2s, scorpine-like antimicrobial peptides, and serine proteases), most of which were less abundant and likely have ancillary functional roles. Finally, although we detected several major toxin families at disproportionate transcriptomic and/or proteomic abundances, we did not identify significant sex-based variation in Hadrurus venoms.

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

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

Global parameter optimisation and sensitivity analysis of antivenom pharmacokinetics and pharmacodynamics

In recent years it has become possible to design snakebite antivenoms with diverse pharmacokinetic properties. Owing to the pharmacokinetic variability of venoms, the choice of antivenom scaffold may influence a treatment's neutralisation coverage. Computation offers a useful medium through which to assess the pharmacokinetics and pharmacodynamics of envenomation-treatment systems, as antivenoms with identical neutralising capacities can be simulated. In this study, we simulate envenomation and treatment with a variety of antivenoms, to define the properties of effective antivenoms. Systemic envenomation and treatment were described using a two-compartment pharmacokinetic model. Treatment of Naja sumatrana and Cryptelytrops purpureomaculatus envenomation was simulated with a set of 200,000 theoretical antivenoms across 10 treatment time delays. These two venoms are well-characterised and have differing pharmacokinetic properties. The theoretical antivenom set varied across molecular weight, dose, kon, koff, and valency. The best and worst treatments were identified using an area under the curve metric, and a global sensitivity analysis was performed to quantify the influence of the input parameters on treatment outcome. The simulations show that scaffolds of diverse molecular formats can be effective. Molecular weight and valency have a negligible direct impact on treatment outcome, however low molecular weight scaffolds offer more flexibility across the other design parameters, particularly when treatment is delayed. The simulations show kon to primarily mediate treatment efficacy, with rates above 105 M−1s−1 required for the most effective treatments. koff has the greatest impact on the performance of less effective scaffolds. While the same scaffold preferences for improved treatment are seen for both model snakes, the parameter bounds for C. purpureomaculatus envenomation are more constrained. This paper establishes a computational framework for the optimisation of antivenom design.

Snakes on a plain: biotic and abiotic factors determine venom compositional variation in a wide-ranging generalist rattlesnake

BackgroundSnake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom.ResultsUsing shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition.ConclusionsOur findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium.

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

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.

Phylogeny-Related Variations in Venomics: A Test in a Subset of Habu Snakes (Protobothrops).

We conducted a comparative analysis to unveil the divergence among venoms from a subset of Old World habu snakes (Protobothrops) in terms of venomic profiles and toxicological and enzymatic activities. A total of 14 protein families were identified in the venoms from these habu snakes, and 11 of them were shared among these venoms. The venoms of five adult habu snakes were overwhelmingly dominated by SVMP (32.56 ± 13.94%), PLA2 (22.93 ± 9.26%), and SVSP (16.27 ± 4.79%), with a total abundance of over 65%, while the subadult P. mangshanensis had an extremely low abundance of PLA2 (1.23%) but a high abundance of CTL (51.47%), followed by SVMP (22.06%) and SVSP (10.90%). Apparent interspecific variations in lethality and enzymatic activities were also explored in habu snake venoms, but no variations in myotoxicity were found. Except for SVSP, the resemblance of the relatives within Protobothrops in other venom traits was estimated to deviate from Brownian motion evolution based on phylogenetic signals. A comparative analysis further validated that the degree of covariation between phylogeny and venom variation is evolutionarily labile and varies among clades of closely related snakes. Our findings indicate a high level of interspecific variation in the venom proteomes of habu snakes, both in the presence or absence and the relative abundance of venom protein families, and that these venoms might have evolved under a combination of adaptive and neutral mechanisms.

Venom phenotype conservation suggests integrated specialization in a lizard-eating snake

Biological specialization reduces the size of niche space while increasing efficiency in the use of available resources. Specialization often leads to phenotypic changes via natural selection aligning with niche space constraints. Commonly observed changes are in size, shape, behavior, and traits associated with feeding. One often selected trait for dietary specialization is venom, which, in snakes, often shows variation dependent on diet across and within species. The Neotropical Blunt-headed Treesnake (Imantodes cenchoa) is a highly specialized, rear-fanged, arboreal, lizard hunter that displays a long thin body, enlarged eyes, and a large Duvernoy's gland. However, toxin characterization of I. cenchoa has never been completed. Here, we use RNA-seq and mass spectrometry to assemble, annotate, and analyze the venom gland transcriptomes of four I. cenchoa from across their range. We find a lack of significant venom variation at the sequence and expression levels, suggesting venom conservation across the species. We propose this conservation provides evidence of a specialized venom repertoire, adapted to maximize efficiency of capturing and processing lizards. Importantly, this study provides the most complete venom gland transcriptomes of I. cenchoa and evidence of venom specialization in a rear-fanged snake, giving insight into selective pressures of venom across all snake species.

High-Throughput Venomics.

In this study, we present high-throughput (HT) venomics, a novel analytical strategy capable of performing a full proteomic analysis of a snake venom within 3 days. This methodology comprises a combination of RP-HPLC-nanofractionation analytics, mass spectrometry analysis, automated in-solution tryptic digestion, and high-throughput proteomics. In-house written scripts were developed to process all the obtained proteomics data by first compiling all Mascot search results for a single venom into a single Excel sheet. Then, a second script plots each of the identified toxins in so-called Protein Score Chromatograms (PSCs). For this, for each toxin, identified protein scores are plotted on the y-axis versus retention times of adjacent series of wells in which a toxin was fractionated on the x-axis. These PSCs allow correlation with parallel acquired intact toxin MS data. This same script integrates the PSC peaks from these chromatograms for semiquantitation purposes. This new HT venomics strategy was performed on venoms from diverse medically important biting species; Calloselasma rhodostoma, Echis ocellatus, Naja pallida, Bothrops asper, Bungarus multicinctus, Crotalus atrox, Daboia russelii, Naja naja, Naja nigricollis, Naja mossambica, and Ophiophagus hannah. Our data suggest that high-throughput venomics represents a valuable new analytical tool for increasing the throughput by which we can define venom variation and should greatly aid in the future development of new snakebite treatments by defining toxin composition.

Individual variation in venom from the left and right glands of snakes of the species Bothrops atrox

Individual variation in venom from the left and right glands of snakes of the species Bothrops atrox

Bungarus fasciatus venom from eastern and north-east India: venom variation and immune cross-reactivity with Indian polyvalent antivenoms

Bungarus fasciatus venom from eastern and north-east India: venom variation and immune cross-reactivity with Indian polyvalent antivenoms

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.).

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Differentiation of snake venom using Raman spectroscopic analysis.

Snake venoms are complex mixtures of different substances, proteins being their predominant components. To study the composition of venoms, methods based on chromatographic separation and mass spectrometric analysis are currently used, requiring the application of a number of sophisticated instruments. To assess the composition of snake venoms, we propose an alternative method based on Raman spectroscopy, which is an express method to study the structural features of different substances, including proteins. The method does not require preliminary preparation of the samples, which are used in small quantities; this makes Raman spectroscopy extremely attractive for venom research. In this work, we have carried out Raman spectroscopic studies on a number of dry venoms from various venomous snakes. Based on the obtained Raman spectra, with the help of mathematical methods of dimensionality reduction and clustering, differentiation of venoms reflecting their composition and the assignment of the venom producing snake to the corresponding family or even genus were performed. The proposed method can be used to analyze both the composition of and variations in venoms of different snake species, including rare and endangered ones.

A current perspective on snake venom composition and constituent protein families.

Snake venoms are heterogeneous mixtures of proteins and peptides used for prey subjugation. With modern proteomics there has been a rapid expansion in our knowledge of snake venom composition, resulting in the venom proteomes of 30% of vipers and 17% of elapids being characterised. From the reasonably complete proteomic coverage of front-fanged snake venom composition (179 species-68 species of elapids and 111 species of vipers), the venoms of vipers and elapids contained 42 different protein families, although 18 were only reported in < 5% of snake species. Based on the mean abundance and occurrence of the 42 protein families, they can be classified into 4 dominant, 6 secondary, 14 minor, and 18 rare protein families. The dominant, secondary and minor categories account for 96% on average of a snake's venom composition. The four dominant protein families are: phospholipase A2 (PLA2), snake venom metalloprotease (SVMP), three-finger toxins (3FTx), and snake venom serine protease (SVSP). The six secondary protein families are: L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRiSP), C-type lectins (CTL), disintegrins (DIS), kunitz peptides (KUN), and natriuretic peptides (NP). Venom variation occurs at all taxonomic levels, including within populations. The reasons for venom variation are complex, as variation is not always associated with geographical variation in diet. The four dominant protein families appear to be the most important toxin families in human envenomation, being responsible for coagulopathy, neurotoxicity, myotoxicity and cytotoxicity. Proteomic techniques can be used to investigate the toxicological profile of a snake venom and hence identify key protein families for antivenom immunorecognition.