Theo Tasoulis
Tara L. Pukala
Geoffrey K. Isbister- 1Clinical Toxicology Research Group, University of Newcastle, Callaghan, NSW, Australia
- 2Department of Chemistry, University of Adelaide, Adelaide, SA, Australia
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%, PLA2–27%, SVMP–2.8%, and SVSP–0.1%, and for vipers: 3FTx–0.5%, PLA2–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.
1 Introduction
Venom proteomes of snakes have been published with increasing frequency since 2004 (Juárez et al., 2004), with a review published in late 2017 of 132 snake species characterised up to that point (Tasoulis and Isbister, 2017). Individual snake venoms consist of three to approximately 20 different recognised toxin protein families (Sanz et al., 2019; Calvette et al., 2021; Sunagar et al., 2021), made up from a total pool of 57 protein families so far reported in published snake venom proteomes. This total includes a number of regulatory proteins and low abundance protein families of unknown functional or biological significance. Each of these toxin protein families are believed to be monophyletic, meaning that each family is the result of a single recruitment event into the venom proteome, followed by up-regulation of expression and orthologous diversification (Casewell, 2012; Casewell et al., 2020; Giorgianni et al., 2020). When averaged across all snake species, the majority of snake venoms are composed of four dominant protein families; phospholipase A2 (PLA2), three-finger toxins (3FTx), snake venom serine protease (SVSP) and snake venom metalloprotease (SVMP) (Tasoulis and Isbister, 2017). A further six secondary protein families make up most of the remaining composition of snake venoms and include; cysteine-rich secretory protein (CRiSP), kuntiz peptides (KUN), l-amino acid oxidase (LAO), natriuretic peptides (NP), C-type lectins (CTL), and disintegrins (DIS) (Tasoulis and Isbister, 2017). Some of these protein families may contain up to 80 different toxins/proteoforms in a single snake species (Tan et al., 2017; Giorgianni et al., 2020; Wang et al., 2020; Sunagar et al., 2021). This amazing diversity of toxins present in some protein families, combined with their synteny and tandemly arrayed location on the genes, is most likely a result of multiple gene duplication events (Wong and Belov, 2012; Hargreaves et al., 2014; Giorgianni et al., 2020).
Characterisation of snake venom proteomes now allows a much more refined way to investigate the evolutionary history of snake venoms. The evolutionary origins of toxins can be approximated by comparisons with phylogeny, providing information on when each of the toxin protein families were ancestrally recruited into the venom proteome. The evolutionary history (orthologous divergence and functional diversification), of these toxin families can then be traced through selected lineages. Venom composition data can be overlaid onto phylogenetic cladograms to reveal if venom composition among related species/genera co-varies with phylogenetic distance. If examples are found in which venom composition diverges from phylogeny, this alerts us to possible unusual evolutionary events that may have occurred and caused changes to venom phenotypes. Comparative venom composition data is also useful for investigating aspects of molecular biology such as gene regulation, gene loss and gene duplication. Used in conjunction with genomic and transcriptomic data, proteomics can provide a wealth of information on the processes governing gene expression, alternative splicing, and other molecular mechanisms responsible for phenotypic divergence. Finally, venom composition data is essential for determining the extent of intra-specific variation, and the influence of ecology and diet versus gene flow on venom evolution (Rautsaw et al., 2019; Schield et al., 2019).
There are two major parts to characterising snake venom proteomes; peptide and protein identification, followed by quantification of the different peptides/proteins. This gives us an estimate of protein diversity and protein family (or individual toxin) abundance. Several previous studies have aimed to establish methods for accurately separating and quantifying protein families in snake venoms (Calvete, 2011; Calvete, 2014; Eichberg et al., 2015; Calderón-Celis et al., 2017; Ghezellou et al., 2019). Some of these studies have emphasised the importance of de-complexing steps prior to mass spectrometry (MS). This workflow usually incorporates reverse-phase high-performance liquid chromatography (RP-HPLC) and 1D sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). However, size-exclusion chromatography (SEC), and 2D electrophoresis (2DE), have also been used with favourable results (Eichberg et al., 2015; Choudhury et al., 2017). Recently, researchers have demonstrated that using a single method to estimate snake venom composition can result in incomplete proteomic coverage and that combining methods may be preferred (Choudhury et al., 2017; Chanda and Mukherjee, 2020). Some recent studies have compared different methodological steps and demonstrated a concerning degree of variation in the results (Choudhury et al., 2017; Modahl et al., 2018).
Protein identification is most often achieved by matching MS derived peptide sequences with public databases. These databases contain protein sequences of previously investigated snake species that have been uploaded, usually in a non-systematic manner. This introduces further errors into the process, most importantly that toxins from new snake species may not be present in the database and so will not be detected. This limits the ability of studies to identify new toxins. The greater the difference in sequence homology of a toxin the greater the risk of it not being detected. Although the corollary to this is that toxins with high sequence similarity may not be distinguished and hence diversity will be underrepresented. A major step forward has been the advent of transcriptomics technology. This now means that species-specific protein databases can be constructed, and toxins in a venom proteome can be directly matched with the same species transcriptome. This allows the identification of new toxins and a better understanding of the true diversity of snake venom proteomes.
To explore the current state of snake venom proteomics and transcriptomics we aimed to establish which venom proteomic methods are most commonly used, and systematically compare the advantages and disadvantages of different techniques, including methods of isolating, identifying and quantifying toxins, matching to databases and workflow patterns. We aimed to investigate the association between methodological (experimental) variability and biological (i.e. intra-specific) variability to determine if some apparent intra-specific variation may be the result of different experimental techniques. Finally, we compared the diversity and abundance of the identified protein families.
2 Materials and Methods
We conducted a systematic review of all studies on snake venom proteomes published between September 2017 and April 2021. We searched the Ovid database > Medline and PubMed, using the search keywords; “snake”, “snake venoms”, “toxins”, “snake venomics”, “elapid”, “viper”, “colubrid”, “proteomics”, “proteomes” and “snake venom proteome”. We identified 143 studies, but only 41 of these met our inclusion criteria for complete venom proteome coverage. We supplemented this search by screening the contents of the following journals for the same time period: Toxins, Toxicon, Toxicon X, Journal of Proteomics, Journal of Proteome Research, International Journal of Biological Macromolecules, Comparative Biochemistry and Physiology Part D. Genomics and Proteomics, Comparative Biochemistry and Physiology Part C. Toxicology and Pharmacology, Proceedings of the Royal Society B, and Acta Tropica. This yielded a further 40 studies to give a total of 81 studies (Supplementary Material S1).
From the included studies we extracted data on the types and amounts of each toxin, to determine toxin diversity and toxin abundance for each snake species. We defined toxin diversity as the total number of different toxins identified in a snake venom. We defined relative toxin abundance as the proportion (percentage) of the whole venom made up by each toxin protein family.
Next, we extracted information on the types of experimental methods and the workflow pattern for each study. Proteomic methodologies were then classified into five different stages: Prefractionation: size-exclusion chromatography (SEC), reverse-phase high performance liquid chromatography (RP-HPLC). Analytical separation: One and two dimensional sodium dodecyl sulfate-polyacrylamide gel electrophoresis (1D SDS-PAGE and 2D SDS-PAGE), liquid chromatography-tandem mass spectrometry (LC-MS/MS). Sample preparation: in-solution trypsin digestion, in-gel trypsin digestion. Protein identification: electrospray ionisation (ESI), liquid chromatography tandem mass spectrometry (LC-MS/MS), matrix-assisted laser desorption ionisation-time of flight mass spectrometry (MALDI-TOF MS). Data analysis; MS derived peptide sequences matched to transcriptome, and MS derived peptide sequences matched to a public data base.
3 Results
3.1 Proteomic Methodology
We identified 81 studies that investigated the composition of snake venoms (Supplementary Table S1), 71 investigated relative toxin abundance (and sometimes diversity), while ten investigated only toxin diversity.
3.1.1 A Comparison of Proteomic Methodologies for Determining Toxin Abundance
Of the 71 studies investigating toxin abundance, 50 started the workflow with RP-HPLC (Figure 1), 20 commenced the workflow with crude venom without any prior chromatographic prefractionation (Figure 2), 13 studies used both RP-HPLC prefractionation workflow in conjunction with a crude venom workflow, and three studies commenced with size-exclusion chromatography (SEC) alone, as a prefractionation step (Figure 3). Only one study used centrifugation as a pre-fractionation step, and one study used hydrophilic interaction liquid chromatography (HILIC) (Figure 3).
FIGURE 1. Summary of the workflows for the 50 studies investigating toxin abundance that commenced with the RP-HPLC prefractionation step (13 of these studies additionally used a crude venom pathway, purple). 23 of the 50 studies (46%), followed the 1D SDS-PAGE →in-gel trypsin digestion → ESI -LC-MS/MS pathway.
FIGURE 2. Summary of the workflows for the 33 studies investigating toxin abundance that commenced with crude venom without a prefractionation step (13 of these studies additionally used the RP-HPLC pathway–Figure 1).
FIGURE 3. Other methods used for characterising snake venom proteomes; hydrophilic interaction liquid chromatography (HILIC), centrifugation, and a summary of the workflows used for the three studies investigating toxin abundance that commenced with a SEC (also called gel filtration), pathway.
For the second step in the RP-HPLC pathway (Figure 1), 43 of these 50 studies used 1D SDS-PAGE, although 11 of these only used the gels for visualisation of the venom. The remaining 32 followed the classic bottom-up proteomics pathway utilising a workflow of RP-HPLC > 1D SDS-PAGE > in-gel tryptic digestion of the peak fractions > LC-MS/MS. 17 studies used RP-HPLC as a prefractionation step for in-solution trypsin digestion prior to MS, and 16 of the studies used RP-HPLC as the prefractionation step for intact MS (top-down proteomics, black lines on right Figure 1).
Thirty three studies used workflows that commenced with crude venom without any prior chromatographic prefractionation (Figure 2). However, 13 of these studies additionally used a chromatographic prefractionation workflow in conjunction with a crude venom workflow. 16 of the 33 crude venom studies commenced with either 1D or 2D electrophoresis, followed by trypsin digestion and then MS. Fourteen used a crude venom → in-solution trypsin digestion → MS pathway.
Of the three studies that used SEC pre-fractionation (Figure 3), two continued with in-solution trypsin digestion of the fractions followed by MS, and two studies also used 1D SDS-PAGE as a visualisation tool.
Protein identification was achieved by matching peptide sequences with a public database in 64 studies, with only 17 studies using transcriptomics (24%) (supplementary Material S2).
For the 71 studies which quantified toxin abundance, the most common method was a three-step process of RP-HPLC peak integration/SDS-PAGE densitometry/mass spectrometry, with MS1 spectral intensity being used more than twice as commonly as MS2 spectral counting (35–16 studies respectively–supplementary Material S2).
3.1.2 Increased Protein Identification by Preliminary Decomplexation of Crude Venom
Several studies employed multiple methods on the same snake species to compare the outcomes. Additionally, some snake species were the focus of several independent studies allowing for a comparison of different techniques (Table 1). One study (Choudhury et al., 2017), characterised the venoms of Indian cobra Naja naja and common krait Bungarus caeruleus, comparing three different methods; tryptic digestion of crude venom, gel-filtration (SEC) fractionation prior to proteomic analysis, and SDS-PAGE fractionation (without prior fractionation using RP-HPLC), followed by LC-MS/MS. For both species, prior fractionation by SEC gave the most comprehensive proteome coverage (Table 1). For N. naja, SEC recovered 75 toxins, crude venom digestion 40 toxins, and SDS-PAGE was least efficient, only recovering 25 toxins. For B. caeruleus, the results showed a similar trend, SEC–34 toxins, crude venom digestion–30 toxins, and SDS-PAGE 23 toxins (Table 1; Supplementary Table S3). Another study (Ghezellou et al., 2021) compared tryptic digestion of crude venom versus preliminary fractionation by SEC in three populations of saw-scaled vipers Echis carinatus sochureki from Iran. For population 1, crude venom digestion recovered 54 proteins compared to 101 proteins using SEC, for population 2, crude venom digestion identified 49 proteins versus 99 proteins using SEC, and for population 3, crude venom digestion identified 53 toxins versus 91 proteins using SEC (Table 1; Supplementary Table S3). A third study (Kunalan et al., 2018), compared trypsin digestion of crude venom with preliminary SEC for the venoms of king cobra Ophiophagus hannah, and Malayan pit viper Calloselasma rhodostoma. They found C. rhodostoma venom contained 47 toxins when using trptic digestion of the crude venom and 99 toxins when venom was prefractionated using SEC. For O. hannah, 76 toxins were identified using trypsin digestion of crude venom and 150 toxins by prior fractionation with SEC (Table 1; Supplementary Table S3). These findings of the importance of decomplexation of venom prior to MS support the conclusions of earlier researchers (Calvete, 2011; Eichberg et al., 2015).
TABLE 1. A comparison of three different methods of proteome coverage of snake venoms with the number of toxins identified for each method: 1. tryptic digestion of crude venom 2. Prior fractionation of crude venom using SEC, and 3. SDS-PAGE separation of venom mass fractions prior to LC-MS/MS. The studies compared three populations of the saw-scaled viper Echis carinatus (Ghezellou et al., 2021), Indian cobra Naja naja, common krait Bungarus caeruleus (Choudhury et al., 2017), king cobra Ophiophagus hannah, and Malayan pit viper Calloselasma rhodostoma (Kunalan et al., 2018).
3.1.3 Increased Protein Identification by Using a Species-Specific Transcriptome Reference and a Public Database Versus a Public Database Only
Two studies (Choudhury et al., 2017; Sunagar et al., 2021), compared the venoms of the common krait Bungarus caeruleus and produced different results for venom diversity. The major difference in methodology between the two studies was in protein identification from MS data. One study identified 46 toxins using three different workflow methodologies, but matching MS peptides to the NCBI (public), database, and the other study used venom from a single snake from Maharashtra North-western India and identified 225 toxins by matching MS pepides to a Bungarus transcriptome and the NCBI database.
Another study characterised the venom proteomes of two species of colubrids, Ahaetulla prasina and Borikenophis portoricensis, using species-specific transcriptome references versus a public database only (Modahl et al., 2018). This study found that the number of peptide fragments/spectra mapped was far greater when species-specific transcriptome references were used. Moreover, the pie charts generated for protein family abundance using species-specific transcriptomes were more in agreement with their SDS-PAGE results than the pie charts obtained from using a public database. However, the study also emphasised the importance of combining species-specific transcriptome references with public databases, due to the risk of mRNA degradation or issues with venom gland transcriptome assemblies, resulting in the loss of some toxin transcripts.
3.1.4 Comparison of Top-Down and Bottom-Up Proteomics
“Bottom-up” (BU) proteomics, centred on peptide-based identification from protein digests, has been invaluable for high-throughput identification and quantification of venom proteomes and is demonstrably the most widely adopted approach in venomics (Figures 1, 2). However, identification on the peptide level is suboptimal for characterizing post-translational modifications and sequence variants due to the “peptide-to-protein” inference problem (Nesvizhskii and Aebersold, 2005). This is a particular challenge for venoms in which the same peptide is often present in multiple different proteoforms, leading to ambiguity when determining the identity of toxins, and in turn complicating the determination of the total number and relative abundance of proteoforms present in the venom (Melani et al., 2017; Ghezellou et al., 2019).
First coined in 2013 (Verano-Braga et al., 2013), “TD venomics” has to date principally focussed on mass measurement and sequence analysis of denatured intact proteins, offering the ability to directly identify proteoforms and localize modifications. Despite the advantages of TD strategies, a number of methodological and technical challenges remain, including issues associated with efficient ionisation, fragmentation and isotope resolution of large analytes, which currently limit routine analysis to proteins of approximately 30 kDa (Brown et al., 2020). Therefore, investigation of viperid venom proteomes containing high molecular weight proteins, such as those from SVMPs LAAO and SVSP toxin families, may suffer from reduced TD proteomics detection. Similarly, a recent study of Echis carinatus sochureki found that TD proteomics failed to detect some low abundance protein families (NGF and AChE) (Ghezellou et al., 2021). Nevertheless, Calvete and co-workers first reported an integrated analysis of the Indonesian King Cobra venom, with a workflow coupling LC-MS/MS, intact mass measurement of reduced and non-reduced proteins, and BU protein identification with attention to locus-specificity (Ainsworth et al., 2018). Given the complementarity of BU and TD proteomics approaches, multiple recent examples have since demonstrated that, together, they facilitate the most thorough identification of the diversity of snake venom components (Göçmen et al., 2015; Melani et al., 2016; Damm et al., 2018).
Continued developments in intact protein analysis further contribute to venom characterisation. For example, the classical TD venomics approach was extended to elucidate the profile of Vipera anatolica senliki venom by an alternative in-source decay (ISD) MALDI based proteomic workflow (Hempel et al., 2020). Although limited somewhat by sensitivity and to sample fractions with minimal complexity, the method was able to overcome challenges with TD identification of high molecular mass venom components. Denaturing fractionation and traditional TD methodology also destroys non-covalent protein-protein and protein-ligand interactions that contribute to biological activity. Therefore employing native MS and native TD proteomics can provide more detailed characterization of protein complexes, resolving combinatorial variations arising from proteolytic protein processing, sequence polymorphism, complexation and post-translational modification (Melani et al., 2016). Finally, beyond identification, a promising workflow combining RP-LC with inductively coupled plasma MS and denaturing MS was also recently reported for absolute quantitation and mass assignment of intact venom proteins (Calderón-Celis et al., 2016).
3.2 Venom Composition
3.2.1 Toxin Abundance
The relative mean (standard deviation) abundance of the four dominant protein families was for elapids; 3FTx -52% (25%), PLA2–27% (24%), SVMP–2.8% (3.1%), and SVSP–0.11% (0.30%), and for vipers, 3FTx–0.5% (2.7%), PLA2–24% (18%), SVMP -27% (16%), and SVSP–12% (7%) (Figure 4). Elapid venoms are dominated by 3FTx and PLA2, and viper venoms are dominated by SVMP, PLA2 and SVSP. Major differences in the venom composition between elapids and vipers include SVMP being a large component in the venom of vipers but not elapids, and the reverse being the case for 3FTx. The secondary protein families LAO/CRiSP/NP/CTL/DIS also made up a larger component of viper venoms than elapids, with a combined mean of 24% in vipers versus 10% in elapids. Colubrid, non-front fanged snake (NFFS), venoms were overwhelmingly dominated by SVMP - 38% (37%), and 3FTx - 42% (34%), with the remainder of their venom being made up of CRiSP -11% (5.6%), and small amounts of CTL -2.8% (5.7%), PLA2 - 2.% (5%) and LAO -0.9% (2.1%) (Figure 4B) (Supplementary Material S3).
FIGURE 4. Relative abundance (as % of the whole venom), of the four dominant and six secondary protein families found in the venoms of elapids (A), colubrids–non front-fanged snakes (Non Front-fanged snake) (B), true vipers (C), and pit vipers (D), showing all 70 species for which abundance data was available.
3.2.2 Toxin Diversity
Total toxin diversity within the venoms of each species was readily recoverable from studies of 18 species of elapids, seven species of true vipers and ten species of pit vipers. The median number of different toxins in elapid venoms was 46–[inter quartile range (IQR): 39 to 79, range: 27–225]. For true vipers, the median was 56 (IQR: 25 to 80, range: 25–99). For pit vipers, the median was 39 (IQR: 31 to 68, range: 22–97).
Toxin diversity within each protein family was readily recoverable for studies of 23 species of elapids, and nine species of vipers (Figure 5). Toxin diversity within protein families for elapids was; 3FTx: median 15, range: 0 to 80. PLA2: median 4, range: 0 to 49. SVSP: median 0, range: 0 to 7. SVMP: median 4, range: 0 to 42. For vipers; 3FTx: range: 0 to 1. PLA2: median 2, range: 0 to 14. SVSP, median 7, range: 0 to 39. SVMP: median 5, range: 2 to 29.
FIGURE 5. Toxin diversity (number of different toxins), for the four dominant and six secondary protein families found in the venoms of elapids [panel (A)], and vipers [panel (B)]. Total number of toxins for each protein family in the venom for 23 species of elapids and nine species of vipers are indicated, in some cases averaged from multiple populations/studies of a species.
3.2.3 Relationship Between Toxin Diversity and Protein Family Abundance
In elapids, there was no correlation between increasing toxin diversity and higher protein family abundance (Figure 6). This was particularly the case for the dominant protein families (Figures 6A–C). The secondary protein families in elapid venoms showed indications of a linear relationship but it was inconsistent. The strongest correlation was for LAO and CRiSP in cobras and kraits (Figures 6D,E), and KUN in mambas (Figure 6F). However, there were, many examples of phylogenetic groupings in the datasets, with species from the same genera forming clusters e.g. mambas for 3FTx, PLA2, KUN, cobras for 3FTx, PLA2, SVMP, sea snakes for 3FTx and PLA2.
FIGURE 6. Scatter plots showing the lack of a clear relationship between toxin diversity (X-axis), and protein family abundance (Y-axis) for the dominant and secondary protein families in elapid venoms. SVSP occurred at such low abundance in elapids it was omitted (see Figure 4A). Colours denote different genera; blue squares = cobras (Naja - seven species), red triangles = mambas (Dendroaspis - four species), black crosses = kraits (Bungarus - five species), purple circles = sea snakes (Hydrophis - two species), black triangle = Collett’s snake Pseudechis colletti, and green open circle = Malayan striped coral snake Calliophis intestinalis.
A comparison of toxin diversity versus protein family abundance in vipers similarly showed no correlation in the dominant protein families (Figure 7).
FIGURE 7. Scatter plots showing the lack of clear relationship between toxin diversity (X-axis), and protein family abundance (Y-axis) for the three dominant protein families in viper venoms (PLA2, SVMP, and SVSP); These are the seven species from Figure 5B green open circle = Trimeresurus nebularis, blue squares = Russell’s vipers Daboia russelii and D. siamensis, black circle = saw-scaled viper Echis carinatus, black cross = puff adder Bitis arietans, and red triangle = Malayan pit viper Calloselasma rhodostoma, black diamond = Protobothrops flavoviridis.
4 Discussion
We compared different methodologies used in 67 studies investigating snake venom proteomics and found that the most common workflow used in just over half of studies to determine protein abundance was prefractionation by RP-HPLC, followed by 1D SDS-PAGE, in-gel trypsin digestion and MS (Figure 1). Our analysis showed that preliminary decomplexation of crude venom by chromatography leads to increased protein identification (Choudhury et al., 2017; Kunalan et al., 2018; Ghezellou et al., 2021). Utilizing different methodological strategies such as combining bottom-up and top-down proteomics, and integrating venom gland transcriptomics and public databases, will maximize proteome coverage (Petras et al., 2019; Ghezellou et al., 2021). Limited work has been carried out on individual versus pooled venom. Pooled venom provides a typical representation of the quantitative toxin abundance, but may give an inflated picture of individual venom complexity (Sanz et al., 2020).
With continued technological advancement, TD proteomics is expected to play an increasing role in characterizing the full diversity of venom protein/peptide families, particularly with regard to more detailed description of post-translation modifications and protein complexes (Melani et al., 2017). Furthermore, it is clear from this survey of recent publications, quantitative profiling of venom proteins is increasingly being described to add further informative value beyond qualitative cataloguing of components. Notably, several variants of label free proteomic quantification approaches are known, and largely operate under the assumption that a linear relationship exists between protein abundance and measured MS-based parameters (either the number of spectral counts per protein or MS1 peak intensities) (Ahrné et al., 2013). However, these approaches still suffer from major drawbacks resulting from differential peptide ionizability, or the problem with missing values (Adams and Collyer, 2015). Such strategies for label-free quantification have been developed using organisms for which comprehensive genomic or transcriptomic databases are available. In these cases, it is assumed that the likelihood of ion selection for MS/MS sequencing is higher for abundant peptides, and that the number of peptide identifications (normalised to protein size, since larger proteins generally give rise to more tryptic peptides), can provide a surrogate measure for the abundance of the parent protein. However, of particular importance in the venomics field, where a comprehensive sequence database is missing, quantification is biased toward successful peptide identifications (Calvette et al., 2021).
In this context, although laborious, Hus et al. (2020), recently demonstrated that the traditional multi-stage venomics protocol for quantitative analysis remains the most accurate method available (Calvete, 2018). These results also highlighted the need for caution in the interpretation of data from such quantitative experiments, particularly for comparisons between different studies where label-free proteomics strategies have been employed.
The proteoform complexity and wide dynamic range of venom proteins means that low abundant protein components can often be excluded from identification on the basis of typical data-dependent acquisition (DDA) modes, whereby the highest abundance precursor ions are preferentially targeted for fragmentation and sequence analysis. In the proteomics field broadly, data-independent acquisition (DIA) strategies have emerged to increase reproducibility and depth of coverage. In brief, for DIA modes, all ions present in a wide mass isolation window are fragmented without selection. In this case, parallel MS/MS sequence ions are in principle generated for all peptide precursors within the mass range of interest. These complex ion patterns therefore have to be deconvoluted to reliably relate the observed fragments to a known peptide sequence which is challenging using a conventional genome-wide species-specific database.
Typically, a project-specific spectral library is first required to be generated from multiple fractionated DDA analyses of the same sample searched against a protein sequence database. Matching of the peptide elution time and fragment ion patterns from DIA data to the spectral library then aids in detection. However, given the limitations already discussed for venomics analysis in producing DDA datasets, particularly in cases where protein sequence databases are incomplete, DIA methods have found limited application in venomics analyses.
Elapid venoms consisted of predominately two protein families, with 3FTx making up about half of the venoms on average and PLA2 making up about one quarter. However, the relative abundance of PLA2 in elapids could be either higher or lower, as many elapid species exhibit a 3FTx/PLA2 dichotomy, with one family making up most of the venom (Lomonte et al., 2016). In contrast to elapids, viper venoms were predominately made up of three protein families, with SVMP and PLA2 making up a quarter each and SVSP an eighth. Venoms of true and pit vipers showed a high degree of similarity, with only minor differences in the abundance of some of the secondary protein families. Among Viperidae, kunitz peptides were only recorded in true vipers, and natriuretic peptides were more abundant in pit vipers.
Interestingly, colubrid (NFFS) venoms showed features of both elapid and viper venoms, sharing a high abundance of 3FTxs with elapids and a high abundance of SVMPs with vipers. Although only a small number of NFFS were investigated, the studies suggest a 3FTx/SVMP dichotomy, with NFFS venoms being either 3FTx dominant or SVMP dominant. NFFS as a group are not generally considered dangerous, with bites to humans usually only presenting painful local inflammation, and rarely systemic effects (Weinstein and Keyler, 2009; Medeiros et al., 2019; Ineich et al., 2020; Castro et al., 2021). The few species that have caused human fatalities appear to all belong to the group containing predominately SVMP [e.g. boomslang (Pla et al., 2017) and African vine snakes (Debono et al., 2020)], but data is still lacking for Asian Nactricines. 3FTx dominant species are yet to be responsible for human fatalities, despite numerous documented bites (e.g. Boiga) (Fritts et al., 1994).
Most venom proteomics studies have focused on the abundance of various toxins and toxin families, and diversity is often not reported or investigated. Toxin diversity is likely to provide significant insights into the biology and evolution of snake venoms. Diversity in snake venoms can be considered as either the number of different protein families in a snake venom, or toxin diversity, being the number of toxins within these protein families. A much greater proportion of viper venoms are made up of secondary protein families compared to elapids (mean values 24–10%; Figure 4). In contrast, elapids appear to have greater toxin diversity within the protein families, compared to vipers. Examples of this are snouted cobra Naja annulifera in which most of the venom is made up of the single protein family 3FTx (78%), but within this family there are 18 different 3FTxs (Tan et al., 2020), or Eastern green mamba, which has been recorded to express 80 different 3FTx proteoforms (Ainsworth et al., 2018).
We compared toxin abundance and toxin diversity within different families and groups of snakes within each family. It appears that there is no obvious association between toxin diversity and toxin abundance. In fact, venom composition more closely followed phylogenetic groupings (e.g. cobras and mambas–Figure 6), than any simple relationship between abundance and diversity. This highlights the strong influence of phylogeny on venom composition.
For accurate comparisons of toxin diversity to be made between snake species, definitions of toxin nomenclature and toxin classification need to be agreed upon and adopted. Determining whether two toxins are simply proteoforms or different toxins is not simple. These issues have already begun to be addressed and rational solutions offered (King et al., 2008; Smith and Kelleher, 2013). Recommendations include referring to all protein variations arising from a single gene as proteoforms (in preference to the terms isoforms and protein species), with products of different genes (proteins), being identified as toxins. The term proteoforms also includes variations caused by post-translational modifications. Additionally, a rational nomenclature for peptide toxins has been proposed by King et al. (2008). Further refinements would include establishing guidelines for minimum reporting requirements in publications, such as standardised presentation of data. For example, venom proteome pie charts should state the percentage of whole venom for each protein family, followed by the number of toxins for each protein family in brackets. Papers should also include a summary table of the supplementary toxin tables in the results, which many studies have already adopted. One concerning finding of this review was examples of studies in which there was mislabelling of toxin diversity with relative abundance, or presenting diversity results using vague or ambiguous language. Adoption of the terms toxin diversity and toxin abundance will eliminate this confusion and should be provided in all studies.
A further complication when trying to evaluate toxin diversity is that different results can be achieved depending on whether toxins are matched to a public database or a transcriptome. Public databases may not include a full complement of the toxins present in the venom, resulting in an underestimation of true toxin diversity. Conversely, when matching to a transcriptome, the same peptide sequence can match to several RNA splice variants. This can result in an inflated picture of diversity in the venom proteome. Without a reference genome or the latest technology, single-molecule real time (SMRT) RNA sequencing, it may be difficult to determine if these are different toxins arising from different genes or are merely different mRNA isoforms (toxin proteoforms - some of which may not possess toxic properties). However, these issues will not affect abundance determinations, because this reflects the total amount of the toxin family, not the number of different toxins for each family. We found in this review that transcriptome databases are now used in approximately 25% of all studies (17/71 studies). A final issue which may require consideration is investigating the current status of the curation of public toxin databases, to determine if greater oversight required.
Venomics is essentially a tool that allows us to better understand snake venoms from an evolutionary, biological and clinical perspective. Not only does it provide a simple measure of the abundance of different protein families within snake venoms, it provides information on the diversity of these different protein families and the diversity of individual toxins within the protein families. In the development of antivenoms to treat snake envenoming, an understanding of the difference between protein family abundance/diversity and actual toxin diversity is essential. Antibodies are likely to cross-neutralise toxins within the same toxin family, potentially those in the same toxin family but in different snake venoms (Isbister et al., 2010; Silva et al., 2016). In this way, identifying the toxin families that are the most medically important and focusing antivenoms on these will improve the efficacy of antivenom treatment.
There has been an exponential increase in the characterisation of snake venom proteomes in the last 2 decades (Calvette et al., 2021), and a corresponding increase in the use of this data for inter-species and inter-generic comparisons to unravel evolutionary histories/processes, and medical implications (Petras et al., 2011; Lomonte et al., 2016; Strickland et al., 2018; Barua and Mikheyev, 2019; Jackson et al., 2019; Barua and Mikheyev, 2020; Damm et al., 2021). Given these developments, it is essential that equivalent comparisons are being made, so there is an important need for clarification of the best methodologies, as well as an acceptance of accurate definitions, reporting and nomenclature, in this important and rapidly developing field of research.
Author Contributions
Conceptualisation—GI, TT, and TP, Data extraction and collation—TT, Data curation—TT, GI, and TP, Data analysis—GI, TP, and TT, Manuscript original draft—TT, Manuscript editing—GI and TP, Funding—GI.
Funding
GI is supported by an NHMRC Senior Research Fellowship 1061041. The study was funded by an NHMRC Centre for Research Excellence 1110343.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.768015/full#supplementary-material
References
Adams, D. C., and Collyer, M. L. (2015). Permutation Tests for Phylogenetic Comparative Analyses of High-Dimensional Shape Data: what You Shuffle Matters. Evolution 69 (3), 823–829. doi:10.1111/evo.12596
Ahrné, E., Molzahn, L., Glatter, T., and Schmidt, A. (2013). Critical Assessment of Proteome-wide Label-free Absolute Abundance Estimation Strategies. PROTEOMICS 13 (17), 2567–2578. doi:10.1002/pmic.201300135
Ainsworth, S., Petras, D., Engmark, M., Süssmuth, R. D., Whiteley, G., Albulescu, L. O., et al. (2018). 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. J. Proteomics 172, 173–189. doi:10.1016/j.jprot.2017.08.016
Barua, A., and Mikheyev, A. S. (2019). Many Options, Few Solutions: Over 60 My Snakes Converged on a Few Optimal Venom Formulations. Mol. Biol. Evol. 36 (9), 1964–1974. doi:10.1093/molbev/msz125
Barua, A., and Mikheyev, A. S. (2020). Toxin Expression in Snake Venom Evolves Rapidly with Constant Shifts in Evolutionary Rates. Proc. Biol. Sci. 287 (1926), 20200613. doi:10.1098/rspb.2020.0613
Brown, K. A., Melby, J. A., Roberts, D. S., and Ge, Y. (2020). Top-down Proteomics: Challenges, Innovations, and Applications in Basic and Clinical Research. Expert Rev. Proteomics 17 (10), 719–733. doi:10.1080/14789450.2020.1855982
Calderón-Celis, F., Diez-Fernández, S., Costa-Fernández, J. M., Encinar, J. R., Calvete, J. J., and Sanz-Medel, A. (2016). Elemental Mass Spectrometry for Absolute Intact Protein Quantification without Protein-Specific Standards: Application to Snake Venomics. Anal. Chem. 88 (19), 9699–9706. doi:10.1021/acs.analchem.6b02585
Calderón-Celis, F., Cid-Barrio, L., Encinar, J. R., Sanz-Medel, A., and Calvete, J. J. (2017). Erratum to "Absolute Venomics: Absolute Quantification of Intact Venom Proteins through Elemental Mass Spectrometry". J. Proteomics 165, 138–140. doi:10.1016/j.jprot.2017.07.009
Calvette, J. J., Lomonte, B., Saviola, A. J., Bonilla, F., Sasa, M., Williams, D. J., et al. (2021). Mutual Enlightenment: A Toolbox of Concepts and Methods for Integrating Evolutionary and Clinical Toxinology via Snake Venomics and the Contextual Stance. Toxicon X 9-10, 100070–100110. doi:10.1016/j.toxcx.2021.100070
Calvete, J. J. (2011). Proteomic Tools against the Neglected Pathology of Snake Bite Envenoming. Expert Rev. Proteomics 8 (6), 739–758. doi:10.1586/epr.11.61
Calvete, J. J. (2014). Next-generation Snake Venomics: Protein-Locus Resolution through Venom Proteome Decomplexation. Expert Rev. Proteomics 11 (3), 315–329. doi:10.1586/14789450.2014.900447
Calvete, J. J. (2018). Snake Venomics - from Low-Resolution Toxin-Pattern Recognition to Toxin-Resolved Venom Proteomes with Absolute Quantification. Expert Rev. Proteomics 15 (7), 555–568. doi:10.1080/14789450.2018.1500904
Casewell, N. R., Jackson, T. N. W., Laustsen, A. H., and Sunagar, K. (2020). Causes and Consequences of Snake Venom Variation. Trends Pharmacol. Sci. 41 (8), 570–581. doi:10.1016/j.tips.2020.05.006
Casewell, N. R. (2012). On the Ancestral Recruitment of Metalloproteinases into the Venom of Snakes. Toxicon 60 (4), 449–454. doi:10.1016/j.toxicon.2012.02.006
Castro, F. C., Nogueira de Souza, S., Maria de Almeida-Santos, S., Miyaji, K. T., and Roberto de Medeiros, C. (2021). Bites by Philodryas Olfersii (Lichtenstein, 1823) and Philodryas Aestiva (Duméril, Bibron and Duméril, 1854) (Serpentes, Dipsadidae) in São Paulo, Brazil: A Retrospective Observational Study of 155 Cases. Toxicon 197, 55–64. doi:10.1016/j.toxicon.2021.04.014
Chanda, A., and Mukherjee, A. K. (2020). Quantitative Proteomics to Reveal the Composition of Southern India Spectacled Cobra (Naja naja) Venom and its Immunological Cross-Reactivity towards Commercial Antivenom. Int. J. Biol. Macromol 160, 224–232. doi:10.1016/j.ijbiomac.2020.05.106
Choudhury, M., McCleary, R. J. R., Kesherwani, M., Kini, R. M., and Velmurugan, D. (2017). Comparison of Proteomic Profiles of the Venoms of Two of the 'Big Four' Snakes of India, the Indian Cobra (Naja naja) and the Common Krait (Bungarus Caeruleus), and Analyses of Their Toxins. Toxicon 135, 33–42. doi:10.1016/j.toxicon.2017.06.005
Damm, M., Hempel, B. F., Nalbantsoy, A., and Süssmuth, R. D. (2018). Comprehensive Snake Venomics of the Okinawa Habu Pit Viper, Protobothrops Flavoviridis, by Complementary Mass Spectrometry-Guided Approaches. Molecules 23 (8), 1893. doi:10.3390/molecules23081893
Damm, M., Hempel, B. F., and Süssmuth, R. D. (2021). Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations. Toxins (Basel) 13 (6), 427. doi:10.3390/toxins13060427
Debono, J., Dashevsky, D., Nouwens, A., and Fry, B. G. (2020). The Sweet Side of Venom: Glycosylated Prothrombin Activating Metalloproteases from Dispholidus Typus (Boomslang) and Thelotornis Mossambicanus (Twig Snake). Comp. Biochem. Physiol. C Toxicol. Pharmacol. 227, 108625. doi:10.1016/j.cbpc.2019.108625
Eichberg, S., Sanz, L., Calvete, J. J., and Pla, D. (2015). Constructing Comprehensive Venom Proteome Reference Maps for Integrative Venomics. Expert Rev. Proteomics 12 (5), 557–573. doi:10.1586/14789450.2015.1073590
Fritts, T. H., McCoid, M. J., and Haddock, R. L. (1994). Symptoms and Circumstances Associated with Bites by the Brown Tree Snake (Colubridae: Boiga Irregularis) on Guam. J. Herpetology 28 (1), 27–33. doi:10.2307/1564676
Ghezellou, P., Garikapati, V., Kazemi, S. M., Strupat, K., Ghassempour, A., and Spengler, B. (2019). A Perspective View of Top-Down Proteomics in Snake Venom Research. Rapid Commun. Mass. Spectrom. 33 Suppl 1 (S1), 20–27. doi:10.1002/rcm.8255
Ghezellou, P., Albuquerque, W., Garikapati, V., Casewell, N. R., Kazemi, S. M., Ghassempour, A., et al. (2021). Integrating Top-Down and Bottom-Up Mass Spectrometric Strategies for Proteomic Profiling of Iranian Saw-Scaled Viper, Echis carinatus Sochureki, Venom. J. Proteome Res. 20 (1), 895–908. doi:10.1021/acs.jproteome.0c00687
Giorgianni, M. W., Dowell, N. L., Griffin, S., Kassner, V. A., Selegue, J. E., and Carroll, S. B. (2020). The Origin and Diversification of a Novel Protein Family in Venomous Snakes. Proc. Natl. Acad. Sci. U S A. 117 (20), 10911–10920. doi:10.1073/pnas.1920011117
Göçmen, B., Heiss, P., Petras, D., Nalbantsoy, A., and Süssmuth, R. D. (2015). Mass Spectrometry Guided Venom Profiling and Bioactivity Screening of the Anatolian Meadow Viper, Vipera Anatolica. Toxicon 107, 163–174. doi:10.1016/j.toxicon.2015.09.013
Hargreaves, A. D., Swain, M. T., Hegarty, M. J., Logan, D. W., and Mulley, J. F. (2014). Restriction and Recruitment-Gene Duplication and the Origin and Evolution of Snake Venom Toxins. Genome Biol. Evol. 6 (8), 2088–2095. doi:10.1093/gbe/evu166
Hempel, B. F., Damm, M., Mrinalini, , Göçmen, B., Karış, M., Nalbantsoy, A., et al. (2020). Extended Snake Venomics by Top-Down In-Source Decay: Investigating the Newly Discovered Anatolian Meadow Viper Subspecies, Vipera Anatolica Senliki. J. Proteome Res. 19 (4), 1731–1749. doi:10.1021/acs.jproteome.9b00869
Hus, K. K., Marczak, Ł., Petrilla, V., Petrillová, M., Legáth, J., and Bocian, A. (2020). Different Research Approaches in Unraveling the Venom Proteome of Naja Ashei. Biomolecules 10 (9), 1282. doi:10.3390/biom10091282
Ineich, I., Girard, F., and Weinstein, S. A. (2020). Local Envenoming by the Schokari Sand Racer, Psammophis Schokari Forskål, 1775 (Serpentes, Psammophiidae) and a Brief Review of Reported Bites by Sand Racers (Psammophis spp.). Toxicon 185, 72–75. doi:10.1016/j.toxicon.2020.06.020
Isbister, G. K., O'Leary, M. A., Hagan, J., Nichols, K., Jacoby, T., Davern, K., et al. (2010). Cross-neutralisation of Australian Brown Snake, Taipan and Death Adder Venoms by Monovalent Antibodies. Vaccine 28 (3), 798–802. doi:10.1016/j.vaccine.2009.10.055
Jackson, T. N. W., Jouanne, H., and Vidal, N. (2019). Snake Venom in Context: Neglected Clades and Concepts. Front. Ecol. Evol. 7, 332. doi:10.3389/fevo.2019.00332
Juárez, P., Sanz, L., and Calvete, J. J. (2004). Snake Venomics: Characterization of Protein Families in Sistrurus Barbouri Venom by Cysteine Mapping, N-Terminal Sequencing, and Tandem Mass Spectrometry Analysis. Proteomics 4 (2), 327–338. doi:10.1002/pmic.200300628
King, G. F., Gentz, M. C., Escoubas, P., and Nicholson, G. M. (2008). A Rational Nomenclature for Naming Peptide Toxins from Spiders and Other Venomous Animals. Toxicon 52 (2), 264–276. doi:10.1016/j.toxicon.2008.05.020
Kunalan, S., Othman, I., Syed Hassan, S., and Hodgson, W. C. (2018). Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma Rhodostoma (Malayan Pit Viper) and Ophiophagus Hannah (King Cobra). Toxins (Basel) 10 (11), 434. doi:10.3390/toxins10110434
Lomonte, B., Rey-Suárez, P., Fernández, J., Sasa, M., Pla, D., Vargas, N., et al. (2016). Venoms of Micrurus Coral Snakes: Evolutionary Trends in Compositional Patterns Emerging from Proteomic Analyses. Toxicon 122, 7–25. doi:10.1016/j.toxicon.2016.09.008
Medeiros, C. R., Souza, S. N., Silva, M. Cd., Ventura, J. S., Piorelli, R. O., and Puorto, G. (2019). Bites by Tomodon Dorsatus (Serpentes, Dipsadidae): Clinical and Epidemiological Study of 86 Cases. Toxicon 162, 40–45. doi:10.1016/j.toxicon.2019.03.005
Melani, R. D., Skinner, O. S., Fornelli, L., Domont, G. B., Compton, P. D., and Kelleher, N. L. (2016). Mapping Proteoforms and Protein Complexes From King Cobra Venom Using Both Denaturing and Native Top-Down Proteomics. Mol. Cel Proteomics 15 (7), 2423–2434. doi:10.1074/mcp.M115.056523
Melani, R. D., Nogueira, F. C. S., and Domont, G. B. (2017). It Is Time for Top-Down Venomics. J. Venom Anim. Toxins Incl Trop. Dis. 23 (1), 44. doi:10.1186/s40409-017-0135-6
Modahl, C. M., Frietze, S., and Mackessy, S. P. (2018). Transcriptome-facilitated Proteomic Characterization of Rear-Fanged Snake Venoms Reveal Abundant Metalloproteinases with Enhanced Activity. J. Proteomics 187, 223–234. doi:10.1016/j.jprot.2018.08.004
Nesvizhskii, A. I., and Aebersold, R. (2005). Interpretation of Shotgun Proteomic Data: the Protein Inference Problem. Mol. Cel Proteomics 4 (10), 1419–1440. doi:10.1074/mcp.R500012-MCP200
Petras, D., Sanz, L., Segura, A., Herrera, M., Villalta, M., Solano, D., et al. (2011). Snake Venomics of African Spitting Cobras: Toxin Composition and Assessment of Congeneric Cross-Reactivity of the Pan-African EchiTAb-Plus-ICP Antivenom by Antivenomics and Neutralization Approaches. J. Proteome Res. 10 (3), 1266–1280. doi:10.1021/pr101040f
Petras, D., Hempel, B. F., Göçmen, B., Karis, M., Whiteley, G., Wagstaff, S. C., et al. (2019). Intact Protein Mass Spectrometry Reveals Intraspecies Variations in Venom Composition of a Local Population of Vipera Kaznakovi in Northeastern Turkey. J. Proteomics 199, 31–50. doi:10.1016/j.jprot.2019.02.004
Pla, D., Sanz, L., Whiteley, G., Wagstaff, S. C., Harrison, R. A., Casewell, N. R., et al. (2017). What Killed Karl Patterson Schmidt? Combined Venom Gland Transcriptomic, Venomic and Antivenomic Analysis of the South African green Tree Snake (The Boomslang), Dispholidus Typus. Biochim. Biophys. Acta Gen. Subj 1861 (4), 814–823. doi:10.1016/j.bbagen.2017.01.020
Rautsaw, R. M., Hofmann, E. P., Margres, M. J., Holding, M. L., Strickland, J. L., Mason, A. J., et al. (2019). Intraspecific Sequence and Gene Expression Variation Contribute Little to Venom Diversity in Sidewinder Rattlesnakes ( Crotalus cerastes). Proc. Biol. Sci. 286 (1906), 20190810. doi:10.1098/rspb.2019.0810
Sanz, L., de Freitas-Lima, L. N., Quesada-Bernat, S., Graça-de-Souza, V. K., Soares, A. M., Calderón, L. A., et al. (2019). Comparative Venomics of Brazilian Coral Snakes: Micrurus Frontalis, Micrurus Spixii Spixii, and Micrurus Surinamensis. Toxicon 166, 39–45. doi:10.1016/j.toxicon.2019.05.001
Sanz, L., Quesada-Bernat, S., Pérez, A., De Morais-Zani, K., SantˈAnna, S. S., Hatakeyama, D. M., et al. (2020). Danger in the Canopy. Comparative Proteomics and Bioactivities of the Venoms of the South American Palm Pit Viper Bothrops Bilineatus Subspecies Bilineatus and Smaragdinus and Antivenomics of B. B. Bilineatus (Rondônia) Venom against the Brazilian Pentabothropic Antivenom. J. Proteome Res. 19 (8), 3518–3532. doi:10.1021/acs.jproteome.0c00337
Schield, D. R., Perry, B. W., Adams, R. H., Card, D. C., Jezkova, T., Pasquesi, G. I. M., et al. (2019). Allopatric Divergence and Secondary Contact with Gene Flow: a Recurring Theme in Rattlesnake Speciation. Biol. J. Linn. Soc. 128 (1), 149–169. doi:10.1093/biolinnean/blz077
Silva, A., Hodgson, W. C., and Isbister, G. K. (2016). Cross-Neutralisation of In Vitro Neurotoxicity of Asian and Australian Snake Neurotoxins and Venoms by Different Antivenoms. Toxins (Basel) 8 (10), 302. doi:10.3390/toxins8100302
Smith, L. M., and Kelleher, N. L. (2013). Proteoform: a Single Term Describing Protein Complexity. Nat. Methods 10 (3), 186–187. doi:10.1038/nmeth.2369
Strickland, J. L., Mason, A. J., Rokyta, D. R., and Parkinson, C. L. (2018). Phenotypic Variation in Mojave Rattlesnake (Crotalus Scutulatus) Venom Is Driven by Four Toxin Families. Toxins (Basel) 10 (4), 135. doi:10.3390/toxins10040135
Sunagar, K., Khochare, S., Senji Laxme, R. R., Attarde, S., Dam, P., Suranse, V., et al. (2021). A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India. Toxins 13 (1), 69. doi:10.3390/toxins13010069
Tan, N. H., Wong, K. Y., and Tan, C. H. (2017). Venomics of Naja Sputatrix, the Javan Spitting Cobra: A Short Neurotoxin-Driven Venom Needing Improved Antivenom Neutralization. J. Proteomics 157, 18–32. doi:10.1016/j.jprot.2017.01.018
Tan, K. Y., Wong, K. Y., Tan, N. H., and Tan, C. H. (2020). Quantitative Proteomics of Naja Annulifera (Sub-Saharan Snouted Cobra) Venom and Neutralization Activities of Two Antivenoms in Africa. Int. J. Biol. Macromol 158, 605–616. doi:10.1016/j.ijbiomac.2020.04.173
Tasoulis, T., and Isbister, G. K. (2017). A Review and Database of Snake Venom Proteomes. Toxins (Basel) 9 (9), 290. doi:10.3390/toxins9090290
Verano-Braga, T., Dutra, A. A., León, I. R., Melo-Braga, M. N., Roepstorff, P., Pimenta, A. M., et al. (2013). Moving Pieces in a Venomic Puzzle: Unveiling Post-translationally Modified Toxins from Tityus Serrulatus. J. Proteome Res. 12 (7), 3460–3470. doi:10.1021/pr4003068
Wang, C. R., Bubner, E. R., Jovcevski, B., Mittal, P., and Pukala, T. L. (2020). Interrogating the Higher Order Structures of Snake Venom Proteins Using an Integrated Mass Spectrometric Approach. J. Proteomics 216, 103680. doi:10.1016/j.jprot.2020.103680
Weinstein, S. A., and Keyler, D. E. (2009). Local Envenoming by the Western Hognose Snake (Heterodon Nasicus): A Case Report and Review of Medically Significant Heterodon Bites. Toxicon 54 (3), 354–360. doi:10.1016/j.toxicon.2009.04.015
Keywords: snake, venom, proteomics, toxin, protein family classification, transcriptomics, mass spectrometry
Citation: Tasoulis T, Pukala TL and Isbister GK (2022) Investigating Toxin Diversity and Abundance in Snake Venom Proteomes. Front. Pharmacol. 12:768015. doi: 10.3389/fphar.2021.768015
Received: 31 August 2021; Accepted: 10 December 2021;
Published: 14 January 2022.
Edited by:
Choo Hock Tan, University of Malaya, MalaysiaReviewed by:
Juan J. Calvete, Spanish National Research Council (CSIC), SpainKaren Zani, Butantan Institute, Brazil
Bhargab Kalita, Amrita Vishwa Vidyapeetham (kochi campus), India
Copyright © 2022 Tasoulis, Pukala and Isbister. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Geoffrey K. Isbister, Z2VvZmYuaXNiaXN0ZXJAZ21haWwuY29t