Venom gland EST analysis of the saw-scaled viper, Echis ocellatus, reveals novel α9β1 integrin-binding motifs in venom metalloproteinases and a new group of putative toxins, renin-like aspartic proteases

Dabsylated Bradykinin Is Cleaved by Snake Venom Proteases from Echis ocellatus

The vasoactive peptide bradykinin (BK) is an important member of the renin-angiotensin system. Its discovery is tightly interwoven with snake venom research, because it was first detected in plasma following the addition of viper venom. While the fact that venoms liberate BK from a serum globulin fraction is well described, its destruction by the venom has largely gone unnoticed. Here, BK was found to be cleaved by snake venom metalloproteinases in the venom of Echis ocellatus, one of the deadliest snakes, which degraded its dabsylated form (DBK) in a few minutes after Pro7 (RPPGFSP↓FR). This is a common cleavage site for several mammalian proteases such as ACE, but is not typical for matrix metalloproteinases. Residual protease activity < 5% after addition of EDTA indicated that DBK is also cleaved by serine proteases to a minor extent. Mass spectrometry-based protein analysis provided spectral proof for several peptides of zinc metalloproteinase-disintegrin-like Eoc1, disintegrin EO4A, and three serine proteases in the venom.

Toxicological analyses of the venoms of Nigerian vipers Echis ocellatus and Bitis arietans

BACKGROUND

Among the medically important snakes in Nigeria, Echis ocellatus and Bitis arietans have the most lethal venom. These venoms were classified according to the presence of snake venom metalloproteinases (SVMPs), snake venom phospholipase A2 (PLA2s), and snake venom serine proteases (SVSPs). Toxicological analyzes were performed to understand the significance of different protein families in venoms.

METHODS

Proteins were separated from venom using column chromatography. The skin and footpad of mice were used to determine hemorrhagic and edematogenic activities. Caprine blood plasma was used to test fibrinolytic activity in vitro.

RESULTS

The results showed that, compared to the crude venom, the SVMP fraction induced hemorrhagic effects with a diameter of 26.00 ± 1.00 mm in E. ocellatus and 21.33 ± 1.52 mm in B. arietans. Both SVSP and SVMP had anticoagulant effects; however, the SVSP fraction had a stronger effect, with a longer anticoagulation time of 30.00 ± 3.00 min in E. ocellatus and 26.00 ± 2.00 min in B. arietans. These main venom toxins, SVMPs, SVSPs, and PLA2, were found to have edema-forming effects that were optimal at 2 h after envenomation. PLA2s had the highest edema-inducing activity, with onset 30 min after envenomation.

CONCLUSIONS

Given the importance of SVMPs in altering the integrity of the membrane structure and impairing the blood coagulation system, an antivenom that can specifically neutralize its activity could inhibit the hemorrhage effects of the venoms.

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933–2022)

Simple Summary

Around the world, snake envenomation poses a serious health risk. Proteins with pharmacological effects are present in snake venom. Recent studies elaborate snake venom and its potential application, including as a cancer drug and antibacterial substances. Our study aimed to analyze the global profile of the literature in snake venom research from documents indexed in the Scopus database between 1933 and 2022. In total, 2999 documents were published with Brazil showing the highest productivity. Antivenom, proteomics, and transcriptomics are emerging as hot topics on a global scale. The present study offers a distinctive overview of snake venom research conducted worldwide.

Abstract

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

De Novo Genome Assembly Highlights the Role of Lineage-Specific Duplications in the Evolution of Venom in Fea’s Viper

Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea’s Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A₂ (PLA₂) proteins expressed by genes that likely arose from lineage-specific PLA₂ gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.

Antivenom: An immunotherapy for the treatment of snakebite envenoming in sub-Saharan Africa

Snakebite envenoming (SBE) leads to significant morbidity and mortality, resulting in over 90,000 deaths and approximately 400,000 amputations annually. In sub-Saharan Africa (SSA) alone, SBE accounts for over 30,000 deaths per annum. Since 2017, SBE has been classified as a priority Neglected Tropical Disease (NTD) by the World Health Organisation (WHO). The major species responsible for mortality from SBE within SSA are from the Bitis, Dendroaspis, Echis and Naja genera. Pharmacologically active toxins such as metalloproteinases, serine proteinases, 3-finger toxins, kunitz-type toxins, and phospholipase A₂s are the primary snake venom components. These toxins induce cytotoxicity, coagulopathy, hemorrhage, and neurotoxicity in envenomed victims. Antivenom is currently the only available venom-specific treatment for SBE and contains purified equine or ovine polyclonal antibodies, collected from donor animals repeatedly immunized with low doses of adjuvanted venom. The resulting plasma or serum contains a high titre of specific antibodies, which can then be collected and stored until required. The purified antibodies are either whole IgG, monovalent fragment antibody (Fab) or divalent fragment antibody F(ab')₂. Despite pharmacokinetic and pharmacodynamic differences, all three are effective in the treatment of SBE. No antivenom is without adverse reactions but, the level of their impact and severity varies from benign early adverse reactions to the rarely occurring fatal anaphylactic shock. However, the major side effects are largely reversible with immediate administration of adrenaline and corticosteroids. There are 16 different antivenoms marketed within SSA, but the efficacy and safety profiles are only published for less than 50% of these products.

Venom proteomic analysis of medically important Nigerian viper Echis ocellatus and Bitis arietans snake species

Snakebite envenoming remains a neglected tropical disease which poses severe health hazard, especially for the rural inhabitants in Africa. In Nigeria, vipers are responsible for the highest number of deaths. Hydrophilic interaction liquid chromatography coupled with LC-MS/MS was used to analyze the crude venoms of Echis ocellatus (Carpet viper) and Bitis arietans (Puff adder) in order to understand their venom proteomic identities. Results obtained revealed that gel-free proteomic analysis of the crude venoms led to the identification of 85 and 79 proteins, respectively. Seventy-eight (78) proteins were common between the two snake species with a 91.8% similarity score. The identified proteins belong to 18 protein families in E. ocellatus and 14 protein families in B. arietans. Serine proteases (22.31%) and metalloproteinases (21.06%) were the dominant proteins in the venom of B. arietans; while metalloproteinases (34.84%), phospholipase A₂s (21.19%) and serine proteases (15.50%) represent the major toxins in the E. ocellatus venom. Other protein families such as three-finger toxins and cysteine-rich venom proteins were detected in low proportions. This study provides an insight into the venom proteomic analysis of the two Nigerian viper species, which could be useful in identifying the toxin families to be neutralized in case of envenomation.

•

Venom proteomic of Nigeria's most medically important snakes is presented.

•

SVMP, SVSP and PLA₂ were the major toxin families in E. ocellatus and B. arietans.

•

The venom proteomes of these vipers displayed 91.8% similarity in composition.

Moringa oleifera Extract Extenuates Echis ocellatus Venom-Induced Toxicities, Histopathological Impairments and Inflammation via Enhancement of Nrf2 Expression in Rats

Echis ocellatus snakebite causes more fatalities than all other African snake species combined. Moringa oleifera reportedly possesses an antivenom property. Therefore, we evaluated the effectiveness of M. oleifera ethanol extract (MOE) against E. ocellatus venom (EOV) toxicities. Thirty male rats were grouped as follows (n = 5): Group 1 (normal control received saline), groups 2 to 6 were administered intraperitoneally, 0.22 mg/kg (LD50) of EOV. Group 2 was left untreated while group 3 to 6 were treated post-envenoming with 0.2 mL of polyvalent antivenom, 200, 400, and 600 mg/kg of MOE respectively. MOE significantly (p < 0.05) normalized the altered haematological indices and blood electrolytes profiles. MOE attenuated venom-induced cellular dysfunctions, characterized by a significant increase in NRF2, and concomitant downregulation of increased antioxidant enzymes (SOD and CAT) activities in the serum and heart of the treated rats. MOE normalized the elevated TNF-α and IL-1β in serum and heart tissues. Furthermore, the IgG titre value was significantly (p < 0.5) higher in the envenomed untreated group compared to the MOE-treated groups. Hemorrhagic, hemolytic and coagulant activities of the venom were strongly inhibited by the MOE dose, dependently. Lesions noticed on tissues of vital organs of untreated rats were abolished by MOE. Our findings substantiate the effectiveness of MOE as a potential remedy against EOV toxicities.

High throughput screening and identification of coagulopathic snake venom proteins and peptides using nanofractionation and proteomics approaches

Snakebite is a neglected tropical disease that results in a variety of systemic and local pathologies in envenomed victims and is responsible for around 138,000 deaths every year. Many snake venoms cause severe coagulopathy that makes victims vulnerable to suffering life-threating haemorrhage. The mechanisms of action of coagulopathic snake venom toxins are diverse and can result in both anticoagulant and procoagulant effects. However, because snake venoms consist of a mixture of numerous protein and peptide components, high throughput characterizations of specific target bioactives is challenging. In this study, we applied a combination of analytical and pharmacological methods to identify snake venom toxins from a wide diversity of snake species that perturb coagulation. To do so, we used a high-throughput screening approach consisting of a miniaturised plasma coagulation assay in combination with a venom nanofractionation approach. Twenty snake venoms were first separated using reversed-phase liquid chromatography, and a post-column split allowed a small fraction to be analyzed with mass spectrometry, while the larger fraction was collected and dispensed onto 384-well plates. After fraction collection, any solvent present in the wells was removed by means of freeze-drying, after which it was possible to perform a plasma coagulation assay in order to detect coagulopathic activity. Our results demonstrate that many snake venoms simultaneously contain both procoagulant and anticoagulant bioactives that contribute to coagulopathy. In-depth identification analysis from seven medically-important venoms, via mass spectrometry and nanoLC-MS/MS, revealed that phospholipase A2 toxins are frequently identified in anticoagulant venom fractions, while serine protease and metalloproteinase toxins are often associated with procoagulant bioactivities. The nanofractionation and proteomics approach applied herein seems likely to be a valuable tool for the rational development of next-generation snakebite treatments by facilitating the rapid identification and fractionation of coagulopathic toxins, thereby enabling specific targeting of these toxins by new therapeutics such as monoclonal antibodies and small molecule inhibitors.

Rapid Identification of Phospholipase A₂ Transcripts from Snake Venoms

Phospholipase A₂ (PLA₂) is a major component in snake venoms and it is found in many different isoforms. To identify transcripts encoding different PLA₂ isoforms, we developed a simple, rapid procedure. Total RNA was extracted from the venoms of three cottonmouth snakes and two diamondback rattlesnakes, and further reverse-transcribed into complementary DNA (cDNA). Using one pair of cottonmouth PLA₂-specific primers and a Reverse Transcription Polymerase Chain Reaction (RT-PCR) technique, we identified 27 unique full-length PLA₂ transcripts, including nine sequences identical to the previously documented ones in the database and one novel GIII-like PLA₂. Two common transcripts respectively encoding Asp49 and Lys49 PLA₂ isoforms were identified in all three cottonmouth venoms, that contain more PLA₂ transcripts than the diamondback rattlesnake venoms. The placement of cloned PLA₂ transcripts in snake venom PLA₂s was further discussed by phylogenetic analysis. The procedure developed in this study paves the way for accelerated acquisition of transcriptome data on any other venom toxin families. The results obtained are crucial for insight into the structure and function of PLA₂ isoforms for scientific and potential therapeutic purposes.

Does size matter? Venom proteomic and functional comparison between night adder species (Viperidae: Causus) with short and long venom glands

Night adders (Causus species within the Viperidae family) are amphibian specialists and a common source of snakebite in Africa. Some species are unique in that they have the longest venom glands of any viper, extending approximately 10% of the body length. Despite their potential medical importance and evolutionary novelty, their venom has received almost no research attention. In this study, venoms from a short-glanded species (C. lichtensteinii) and from a long-glanded species (C. rhombeatus) were compared using a series of proteomic and bioactivity testing techniques to investigate and compare the toxin composition and functioning of the venoms of these two species. Both C. rhombeatus and C. lichtensteinii were similar in overall venom composition and inhibition of blood coagulation through non-clotting proteolytic cleavage of fibrinogen. While the 1D gel profiles were very similar to each other in the toxin types present, 2D gel analyses revealed isoformic differences within each toxin classes. This variation was congruent with differential efficacy of South African Institute for Medical Research snake polyvalent antivenom, with C. lichtensteinii unaffected at the dose tested while C. rhombeatus was moderately but significantly neutralized. Despite the variation within toxin classes, the similarity in overall venom biochemistry suggests that the selection pressure for the evolution of long glands served to increase venom yield in order to subjugate proportionally large anurans as a unique form of niche partitioning, and is not linked to significant changes in venom function. These results not only contribute to the body of venom evolution knowledge but also highlight the limited clinical management outcomes for Causus envenomations.

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.

Isolation and characterization of renin-like aspartic-proteases from Echis ocellatus venom

Three aspartic proteases (SVAPs) have been isolated from venom of the saw-scaled viper, Echis ocellatus. In confirmation of prior transcriptomic predictions, all three forms match to sequences of either of the two SVAP transcripts (EOC00051 and EOC00123), have a molecular weight of 42 kDa and possess a single N-glycan. The SVAPs act in a renin-like manner, specifically cleaving human and porcine angiotensinogen into angiotensin-1 and possess no general protease activity. Their activity is completely inhibited by the aspartyl protease inhibitor Pepstatin A.

Comparative transcriptome analysis of venom glands from Cotesia vestalis and Diadromus collaris, two endoparasitoids of the host Plutella xylostella

Venoms secreted by the venom gland (VG) of parasitoid wasp help ensure successful parasitism by host immune suppression and developmental regulation. Cotesia vestalis, a larval endoparasitoid, and Diadromus collaris, a pupal endoparasitoid, parasitize the diamondback moth (DBM), Plutella xylostella. To explore and compare the venom components of two endoparasitoids, we sequenced transcriptomes of the VGs and wasp bodies without VGs (BWVGs) of the two endoparasitoids. Statistically enriched GO terms and KEGG pathways of the two VGs compared to respective whole-body background were similar and reflected active protein biosynthesis activities in the two VGs. 1,595 VG specific genes of the D. collaris VG and 1,461 VG specific genes of the C. vestalis VG were identified by comparative transcript profiling. A total of 444 and 513 genes encoding potential secretory proteins were identified and defined as putative venom genes in D. collaris VG and C. vestalis VG, respectively. The putative venom genes of the two wasps showed no significant similarity or convergence. More venom genes were predicted in D. collaris VG than C. vestalis VG, especially hydrolase-coding genes. Differences in the types and quantities of putative venom genes shed light on different venom functions.

Venomics of Vipera berus berus to explain differences in pathology elicited by Vipera ammodytes ammodytes envenomation: Therapeutic implications

Vipera berus berus (Vbb) is the most widely distributed and Vipera ammodytes ammodytes (Vaa) the most venomous viper in Europe. In particular areas of the Old continent their toxic bites constitute a considerable public health problem. To make the current envenomation therapy more effective we have analysed the proteome of Vbb venom and compared it with that of Vaa. We found the proteome of Vbb to be much less complex and to contain smaller levels of particularly snaclecs and sPLA2s. Snaclecs are probably responsible for thrombocytopenia. The neurotoxic sPLA2s, ammodytoxins, are responsible for the most specific feature of the Vaa venom poisoning - induction of signs of neurotoxicity in patients. These molecules were not found in Vbb venom. Both venoms induce haemorrhage and coagulopathy in man. As Vaa and Vbb venoms possess homologous P-III snake venom metalloproteinases, the main haemorrhagic factors, the severity of the haemorrhage is dictated by concentration and specific activity of these molecules. The much greater anticoagulant effect of Vaa venom than that of Vbb venom lies in its higher extrinsic pathway coagulation factor-proteolysing activity and content of ammodytoxins which block the prothrombinase complex formation.

BIOLOGICAL SIGNIFICANCE

Envenomations by venomous snakes constitute a considerable public health problem worldwide, and also in Europe. In the submitted work we analysed the venom proteome of Vipera berus berus (Vbb), the most widely distributed venomous snake in Europe and compared it with the venom proteome of the most venomous viper in Europe, Vipera ammodytes ammodytes (Vaa). We have offered a possible explanation, at the molecular level, for the differences in clinical pictures inflicted by the Vbb and Vaa venoms. We have provided an explanation for the effectiveness of treatment of Vbb envenomation by Vaa antiserum and explained why full protection of Vaa venom poisoning by Vbb antiserum should not be always expected, especially not in cases of severe poisoning. The latter makes a strong case for Vaa antiserum production as we are faced with its shortage due to ceasing of production of two most frequently used products.

Stabilising the Integrity of Snake Venom mRNA Stored under Tropical Field Conditions Expands Research Horizons

BACKGROUND

Snake venoms contain many proteinaceous toxins that can cause severe pathology and mortality in snakebite victims. Interestingly, mRNA encoding such toxins can be recovered directly from venom, although yields are low and quality is unknown. It also remains unclear whether such RNA contains information about toxin isoforms and whether it is representative of mRNA recovered from conventional sources, such as the venom gland. Answering these questions will address the feasibility of using venom-derived RNA for future research relevant to biomedical and antivenom applications.

METHODOLOGY/PRINCIPAL FINDINGS

Venom was extracted from several species of snake, including both members of the Viperidae and Elapidae, and either lyophilized or immediately added to TRIzol reagent. TRIzol-treated venom was incubated at a range of temperatures (4-37°C) for a range of durations (0-48 hours), followed by subsequent RNA isolation and assessments of RNA quantity and quality. Subsequently, full-length toxin transcripts were targeted for PCR amplification and Sanger sequencing. TRIzol-treated venom yielded total RNA of greater quantity and quality than lyophilized venom, and with quality comparable to venom gland-derived RNA. Full-length sequences from multiple Viperidae and Elapidae toxin families were successfully PCR amplified from TRIzol-treated venom RNA. We demonstrated that venom can be stored in TRIzol for 48 hours at 4-19°C, and 8 hours at 37°C, at minimal cost to RNA quality, and found that venom RNA encoded multiple toxin isoforms that seemed homologous (98-99% identity) to those found in the venom gland.

CONCLUSIONS/SIGNIFICANCE

The non-invasive experimental modifications we propose will facilitate the future investigation of venom composition by using venom as an alternative source to venom gland tissue for RNA-based studies, thus obviating the undesirable need to sacrifice snakes for such research purposes. In addition, they expand research horizons to rare, endangered or protected snake species and provide more flexibility to performing fieldwork on venomous snakes in tropical conditions.

Distal M domain of cobra ADAM-like metalloproteinase mediates the binding of positively charged cysteine-rich domain to αvβ3 integrin in the suppression of cell migration

We have previously identified two new P-III type ADAM-like snake venom metalloproteinases (SVMPs), i.e., atragin and kaouthiagin-like, from Taiwan cobra venom and determined their 3D structures with a distinct C- and I-shaped metalloproteinase/disintegrin-like/cysteine-rich (MDC) modular architecture. Herein, we investigated their functional targets to elucidate the role of cobra SVMPs in perturbing wound healing in snakebite victims. We showed that the non-RGD (Arg-Gly-Asp) C-shaped SVMP atragin binds about ten-fold stronger than the RGD-containing I-shaped SVMP kaouthiagin-like to αvβ3 integrin in the surface-immobilized form. Atragin binds to αvβ3 integrin through a novel interaction mode involving distal M and C domains via the RRN sequence motif in the hyper variable loop. In a cell adhesion assay, the adhesion of fibroblasts to atragin was mediated by αvβ3 integrin. Furthermore, atragin inhibited wound healing and suppressed cell migration in a αvβ3 integrin-dependent manner. These results, together with our previous demonstration of non-cytotoxic cobra CTX A5 in targeting αvβ3 integrin, suggest that cobra venom consists of several non-RGD toxins with integrin-binding specificity that could perturb wound healing in snakebite victims.

Gene expression profiling of the venom gland from the Venezuelan mapanare (Bothrops colombiensis) using expressed sequence tags (ESTs)

BACKGROUND

Bothrops colombiensis is a highly dangerous pit viper and responsible for over 70% of snakebites in Venezuela. Although the composition in B. colombiensis venom has been identified using a proteome analysis, the venom gland transcriptome is currently lacking.

RESULTS

We constructed a cDNA library from the venom gland of B. colombiensis, and a set of 729 high quality expressed sequence tags (ESTs) was identified. A total number of 344 ESTs (47.2% of total ESTs) was related to toxins. The most abundant toxin transcripts were metalloproteinases (37.5%), phospholipases A2s (PLA2, 29.7%), and serine proteinases (11.9%). Minor toxin transcripts were linked to waprins (5.5%), C-type lectins (4.1%), ATPases (2.9%), cysteine-rich secretory proteins (CRISP, 2.3%), snake venom vascular endothelium growth factors (svVEGF, 2.3%), L-amino acid oxidases (2%), and other putative toxins (1.7%). While 160 ESTs (22% of total ESTs) coded for translation proteins, regulatory proteins, ribosomal proteins, elongation factors, release factors, metabolic proteins, and immune response proteins. Other proteins detected in the transcriptome (87 ESTs, 11.9% of total ESTs) were undescribed proteins with unknown functions. The remaining 138 (18.9%) cDNAs had no match with known GenBank accessions.

CONCLUSION

This study represents the analysis of transcript expressions and provides a physical resource of unique genes for further study of gene function and the development of novel molecules for medical applications.

Transcriptomics and venomics: implications for medicinal chemistry

Over the last three decades, transcriptomic studies of venom gland cells have continuously evolved, opening up new possibilities for exploring the molecular diversity of animal venoms, a prerequisite for the discovery of new drug candidates and molecular phylogenetics. The molecular complexity of animal venoms is much greater than initially thought. In this review, we describe the different technologies available for transcriptomic studies of venom, from the original individual cloning approaches to the more recent global Next Generation Sequencing strategies. Our understanding of animal venoms is evolving, with the discovery of complex and diverse bio-optimized cocktails of compounds, including mostly peptides and proteins, which are now beginning to be studied by academic and industrial researchers.

Interrogating the Venom of the Viperid Snake Sistrurus catenatus edwardsii by a Combined Approach of Electrospray and MALDI Mass Spectrometry

The complete sequence characterization of snake venom proteins by mass spectrometry is rather challenging due to the presence of multiple isoforms from different protein families. In the present study, we investigated the tryptic digest of the venom of the viperid snake Sistrurus catenatus edwardsii by a combined approach of liquid chromatography coupled to either electrospray (online) or MALDI (offline) mass spectrometry. These different ionization techniques proved to be complementary allowing the identification a great variety of isoforms of diverse snake venom protein families, as evidenced by the detection of the corresponding unique peptides. For example, ten out of eleven predicted isoforms of serine proteinases of the venom of S. c. edwardsii were distinguished using this approach. Moreover, snake venom protein families not encountered in a previous transcriptome study of the venom gland of this snake were identified. In essence, our results support the notion that complementary ionization techniques of mass spectrometry allow for the detection of even subtle sequence differences of snake venom proteins, which is fundamental for future structure-function relationship and possible drug design studies.

Testing the Toxicofera: comparative transcriptomics casts doubt on the single, early evolution of the reptile venom system

The identification of apparently conserved gene complements in the venom and salivary glands of a diverse set of reptiles led to the development of the Toxicofera hypothesis - the single, early evolution of the venom system in reptiles. However, this hypothesis is based largely on relatively small scale EST-based studies of only venom or salivary glands and toxic effects have been assigned to only some putative Toxicoferan toxins in some species. We set out to examine the distribution of these proposed venom toxin transcripts in order to investigate to what extent conservation of gene complements may reflect a bias in previous sampling efforts. Our quantitative transcriptomic analyses of venom and salivary glands and other body tissues in five species of reptile, together with the use of available RNA-Seq datasets for additional species, shows that the majority of genes used to support the establishment and expansion of the Toxicofera are in fact expressed in multiple body tissues and most likely represent general maintenance or "housekeeping" genes. The apparent conservation of gene complements across the Toxicofera therefore reflects an artefact of incomplete tissue sampling. We therefore conclude that venom has evolved multiple times in reptiles.

VTBuilder: a tool for the assembly of multi isoform transcriptomes

Background

Within many research areas, such as transcriptomics, the millions of short DNA fragments (reads) produced by current sequencing platforms need to be assembled into transcript sequences before they can be utilized. Despite recent advances in assembly software, creating such transcripts from read data harboring isoform variation remains challenging. This is because current approaches fail to identify all variants present or they create chimeric transcripts within which relationships between co-evolving sites and other evolutionary factors are disrupted. We present VTBuilder, a tool for constructing non-chimeric transcripts from read data that has been sequenced from sources containing isoform complexity.

Results

We validated VTBuilder using reads simulated from 54 Sanger sequenced transcripts (SSTs) expressed in the venom gland of the saw scaled viper, Echis ocellatus. The SSTs were selected to represent genes from major co-expressed toxin groups known to harbor isoform variants. From the simulated reads, VTBuilder constructed 55 transcripts, 50 of which had a greater than 99% sequence similarity to 48 of the SSTs. In contrast, using the popular assembler tool Trinity (r2013-02-25), only 14 transcripts were constructed with a similar level of sequence identity to just 11 SSTs. Furthermore VTBuilder produced transcripts with a similar length distribution to the SSTs while those produced by Trinity were considerably shorter. To demonstrate that our approach can be scaled to real world data we assembled the venom gland transcriptome of the African puff adder Bitis arietans using paired-end reads sequenced on Illumina’s MiSeq platform. VTBuilder constructed 1481 transcripts from 5 million reads and, following annotation, all major toxin genes were recovered demonstrating reconstruction of complex underlying sequence and isoform diversity.

Conclusion

Unlike other approaches, VTBuilder strives to maintain the relationships between co-evolving sites within the constructed transcripts, and thus increases transcript utility for a wide range of research areas ranging from transcriptomics to phylogenetics and including the monitoring of drug resistant parasite populations. Additionally, improving the quality of transcripts assembled from read data will have an impact on future studies that query these data. VTBuilder has been implemented in java and is available, under the GPL GPU V0.3 license, from http:// http://www.lstmed.ac.uk/vtbuilder.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0389-8) contains supplementary material, which is available to authorized users.

Venom gland transcriptomics for identifying, cataloging, and characterizing venom proteins in snakes

Snake venoms are cocktails of protein toxins that play important roles in capture and digestion of prey. Significant qualitative and quantitative variation in snake venom composition has been observed among and within species. Understanding these variations in protein components is instrumental in interpreting clinical symptoms during human envenomation and in searching for novel venom proteins with potential therapeutic applications. In the last decade, transcriptomic analyses of venom glands have helped in understanding the composition of various snake venoms in great detail. Here we review transcriptomic analysis as a powerful tool for understanding venom profile, variation and evolution.

Novel apigenin based small molecule that targets snake venom metalloproteases

The classical antivenom therapy has appreciably reduced snakebite mortality rate and thus is the only savior drug available. Unfortunately, it considerably fails to shield the viper bite complications like hemorrhage, local tissue degradation and necrosis responsible for severe morbidity. Moreover, the therapy is also tagged with limitations including anaphylaxis, serum sickness and poor availability. Over the last decade, snake venom metalloproteases (SVMPs) are reported to be the primary component responsible for hemorrhage and tissue degradation at bitten site. Thus, antivenom inability to offset viper venom-induced local toxicity has been a basis for an insistent search for SVMP inhibitors. Here we report the inhibitory effect of compound 5d, an apigenin based molecule against SVMPs both in silico and in vivo. Several apigenin analogues are synthesized using multicomponent Ugi reactions. Among them, compound 5d effectively abrogated Echis carinatus (EC) venom-induced local hemorrhage, tissue necrosis and myotoxicity in a dose dependant fashion. The histopathological study further conferred effective inhibition of basement membrane degradation, and accumulation of inflammatory leucocytes at the site of EC venom inoculation. The compound also protected EC venom-induced fibrin and fibrinogen degradation. The molecular docking of compound 5d and bothropasin demonstrated the direct interaction of hydroxyl group of compound with Glu146 present in hydrophobic pocket of active site and does not chelate Zn²⁺. Hence, it is concluded that compound 5d could be a potent agent in viper bite management.

Medically important differences in snake venom composition are dictated by distinct postgenomic mechanisms

Variation in venom composition is a ubiquitous phenomenon in snakes and occurs both interspecifically and intraspecifically. Venom variation can have severe outcomes for snakebite victims by rendering the specific antibodies found in antivenoms ineffective against heterologous toxins found in different venoms. The rapid evolutionary expansion of different toxin-encoding gene families in different snake lineages is widely perceived as the main cause of venom variation. However, this view is simplistic and disregards the understudied influence that processes acting on gene transcription and translation may have on the production of the venom proteome. Here, we assess the venom composition of six related viperid snakes and compare interspecific changes in the number of toxin genes, their transcription in the venom gland, and their translation into proteins secreted in venom. Our results reveal that multiple levels of regulation are responsible for generating variation in venom composition between related snake species. We demonstrate that differential levels of toxin transcription, translation, and their posttranslational modification have a substantial impact upon the resulting venom protein mixture. Notably, these processes act to varying extents on different toxin paralogs found in different snakes and are therefore likely to be as important as ancestral gene duplication events for generating compositionally distinct venom proteomes. Our results suggest that these processes may also contribute to altering the toxicity of snake venoms, and we demonstrate how this variability can undermine the treatment of a neglected tropical disease, snakebite.

Intraspecific venom variation in the medically significant Southern Pacific Rattlesnake (Crotalus oreganus helleri): biodiscovery, clinical and evolutionary implications

Due to the extreme variation of venom, which consequently results in drastically variable degrees of neutralization by CroFab antivenom, the management and treatment of envenoming by Crotalus oreganus helleri (the Southern Pacific Rattlesnake), one of the most medically significant snake species in all of North America, has been a clinician's nightmare. This snake has also been the subject of sensational news stories regarding supposed rapid (within the last few decades) evolution of its venom. This research demonstrates for the first time that variable evolutionary selection pressures sculpt the intraspecific molecular diversity of venom components in C. o. helleri. We show that myotoxic β-defensin peptides (aka: crotamines/small basic myotoxic peptides) are secreted in large amounts by all populations. However, the mature toxin-encoding nucleotide regions evolve under the constraints of negative selection, likely as a result of their non-specific mode of action which doesn't enforce them to follow the regime of the classic predator-prey chemical arms race. The hemorrhagic and tissue destroying snake venom metalloproteinases (SVMPs) were secreted in larger amounts by the Catalina Island and Phelan rattlesnake populations, in moderate amounts in the Loma Linda population and in only trace levels by the Idyllwild population. Only the Idyllwild population in the San Jacinto Mountains contained potent presynaptic neurotoxic phospholipase A2 complex characteristic of Mohave Rattlesnake (Crotalus scutulatus) and Neotropical Rattlesnake (Crotalus durissus terrificus). The derived heterodimeric lectin toxins characteristic of viper venoms, which exhibit a diversity of biological activities, including anticoagulation, agonism/antagonism of platelet activation, or procoagulation, appear to have evolved under extremely variable selection pressures. While most lectin α- and β-chains evolved rapidly under the influence of positive Darwinian selection, the β-chain lectin of the Catalina Island population appears to have evolved under the constraint of negative selection. Both lectin chains were conspicuously absent in both the proteomics and transcriptomics of the Idyllwild population. Thus, we not only highlight the tremendous biochemical diversity in C. o. helleri's venom-arsenal, but we also show that they experience remarkably variable strengths of evolutionary selection pressures, within each toxin class among populations and among toxin classes within each population. The mapping of geographical venom variation not only provides additional information regarding venom evolution, but also has direct medical implications by allowing prediction of the clinical effects of rattlesnake bites from different regions. Such information, however, also points to these highly variable venoms as being a rich source of novel toxins which may ultimately prove to be useful in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results have direct implications for the treatment of envenomed patients. The variable venom profile of Crotalus oreganus helleri underscores the biodiscovery potential of novel snake venoms.

Pharmacological and biochemical studies on the venom of a clinically important viper snake (Echis carinatus) of Pakistan

Echis carinatus (saw-scaled viper) has been the major culprit responsible for serious envenomation casualties throughout the subcontinent. The present study describes the electrophoretic and zymographic characterization of E. carinatus venom and its effect on mammalian smooth muscle. Crude venom showed the presence of disintegrin, PLA2, C-type lectin/lectin-like components, CRISP, Serine protease, l-amino acid oxidase and very high concentrations of SVMPs. E. carinatus venom (1, 10, 30, 50, 100 μg/ml) inhibited the active tension/force of muscle contraction in a time and concentration dependent manner. The observed effects abolished when the venom was heated at 100 °C for 5 min. However, a decrease in bath temperature from 37 °C to 26 °C or an increase in CaCl2 concentration to 5 mM did not prevent the inhibition of contractile activity. The contractile response elicited by exogenous application of 50 mM KCl and 1 μM acetylcholine (ACh) was also significantly inhibited by all venom concentrations. Prior administration of commercially available polyvalent anti-venom partially neutralized and prevented the effect of E. carinatus venom whereas addition of anti-venom at t50 failed to reverse the inhibitory effect. Studies on isolated intestinal muscle indicate involvement of myotoxic and apoptotic components in E. carinatus venom for irreversible damage to muscle tissue.

Venom down under: dynamic evolution of Australian elapid snake toxins

Despite the unparalleled diversity of venomous snakes in Australia, research has concentrated on a handful of medically significant species and even of these very few toxins have been fully sequenced. In this study, venom gland transcriptomes were sequenced from eleven species of small Australian elapid snakes, from eleven genera, spanning a broad phylogenetic range. The particularly large number of sequences obtained for three-finger toxin (3FTx) peptides allowed for robust reconstructions of their dynamic molecular evolutionary histories. We demonstrated that each species preferentially favoured different types of α-neurotoxic 3FTx, probably as a result of differing feeding ecologies. The three forms of α-neurotoxin [Type I (also known as (aka): short-chain), Type II (aka: long-chain) and Type III] not only adopted differential rates of evolution, but have also conserved a diversity of residues, presumably to potentiate prey-specific toxicity. Despite these differences, the different α-neurotoxin types were shown to accumulate mutations in similar regions of the protein, largely in the loops and structurally unimportant regions, highlighting the significant role of focal mutagenesis. We theorize that this phenomenon not only affects toxin potency or specificity, but also generates necessary variation for preventing/delaying prey animals from acquiring venom-resistance. This study also recovered the first full-length sequences for multimeric phospholipase A₂ (PLA₂) ‘taipoxin/paradoxin’ subunits from non-Oxyuranus species, confirming the early recruitment of this extremely potent neurotoxin complex to the venom arsenal of Australian elapid snakes. We also recovered the first natriuretic peptides from an elapid that lack the derived C-terminal tail and resemble the plesiotypic form (ancestral character state) found in viper venoms. This provides supporting evidence for a single early recruitment of natriuretic peptides into snake venoms. Novel forms of kunitz and waprin peptides were recovered, including dual domain kunitz-kunitz precursors and the first kunitz-waprin hybrid precursors from elapid snakes. The novel sequences recovered in this study reveal that the huge diversity of unstudied venomous Australian snakes are of considerable interest not only for the investigation of venom and whole organism evolution but also represent an untapped bioresource in the search for novel compounds for use in drug design and development.

Predicting function from sequence in a large multifunctional toxin family

Venoms contain active substances with highly specific physiological effects and are increasingly being used as sources of novel diagnostic, research and treatment tools for human disease. Experimental characterisation of individual toxin activities is a severe rate-limiting step in the discovery process, and in-silico tools which allow function to be predicted from sequence information are essential. Toxins are typically members of large multifunctional families of structurally similar proteins that can have different biological activities, and minor sequence divergence can have significant consequences. Thus, existing predictive tools tend to have low accuracy. We investigated a classification model based on physico-chemical attributes that can easily be calculated from amino-acid sequences, using over 250 (mostly novel) viperid phospholipase A₂ toxins. We also clustered proteins by sequence profiles, and carried out in-vitro tests for four major activities on a selection of isolated novel toxins, or crude venoms known to contain them. The majority of detected activities were consistent with predictions, in contrast to poor performance of a number of tested existing predictive methods. Our results provide a framework for comparison of active sites among different functional sub-groups of toxins that will allow a more targeted approach for identification of potential drug leads in the future.

The genesis of an exceptionally lethal venom in the timber rattlesnake (Crotalus horridus) revealed through comparative venom-gland transcriptomics

Background

Snake venoms generally show sequence and quantitative variation within and between species, but some rattlesnakes have undergone exceptionally rapid, dramatic shifts in the composition, lethality, and pharmacological effects of their venoms. Such shifts have occurred within species, most notably in Mojave (Crotalus scutulatus), South American (C. durissus), and timber (C. horridus) rattlesnakes, resulting in some populations with extremely potent, neurotoxic venoms without the hemorrhagic effects typical of rattlesnake bites.

Results

To better understand the evolutionary changes that resulted in the potent venom of a population of C. horridus from northern Florida, we sequenced the venom-gland transcriptome of an animal from this population for comparison with the previously described transcriptome of the eastern diamondback rattlesnake (C. adamanteus), a congener with a more typical rattlesnake venom. Relative to the toxin transcription of C. adamanteus, which consisted primarily of snake-venom metalloproteinases, C-type lectins, snake-venom serine proteinases, and myotoxin-A, the toxin transcription of C. horridus was far simpler in composition and consisted almost entirely of snake-venom serine proteinases, phospholipases A₂, and bradykinin-potentiating and C-type natriuretic peptides. Crotalus horridus lacked significant expression of the hemorrhagic snake-venom metalloproteinases and C-type lectins. Evolution of shared toxin families involved differential expansion and loss of toxin clades within each species and pronounced differences in the highly expressed toxin paralogs. Toxin genes showed significantly higher rates of nonsynonymous substitution than nontoxin genes. The expression patterns of nontoxin genes were conserved between species, despite the vast differences in toxin expression.

Conclusions

Our results represent the first complete, sequence-based comparison between the venoms of closely related snake species and reveal in unprecedented detail the rapid evolution of snake venoms. We found that the difference in venom properties resulted from major changes in expression levels of toxin gene families, differential gene-family expansion and loss, changes in which paralogs within gene families were expressed at high levels, and higher nonsynonymous substitution rates in the toxin genes relative to nontoxins. These massive alterations in the genetics of the venom phenotype emphasize the evolutionary lability and flexibility of this ecologically critical trait.

Complex cocktails: the evolutionary novelty of venoms

Venoms have evolved on numerous occasions throughout the animal kingdom. These 'biochemical weapon systems' typically function to facilitate, or protect the producing animal from, predation. Most venomous animals remain unstudied despite venoms providing model systems for investigating predator-prey interactions, molecular evolution, functional convergence, and novel targets for pharmaceutical discovery. Through advances in 'omic' technologies, venom composition data have recently become available for several venomous lineages, revealing considerable complexity in the processes responsible for generating the genetic and functional diversity observed in many venoms. Here, we review these recent advances and highlight the ecological and evolutionary novelty of venom systems.

Unusual stability of messenger RNA in snake venom reveals gene expression dynamics of venom replenishment

Venom is a critical evolutionary innovation enabling venomous snakes to become successful limbless predators; it is therefore vital that venomous snakes possess a highly efficient venom production and delivery system to maintain their predatory arsenal. Here, we exploit the unusual stability of messenger RNA in venom to conduct, for the first time, quantitative PCR to characterise the dynamics of gene expression of newly synthesised venom proteins following venom depletion. Quantitative PCR directly from venom enables real-time dynamic studies of gene expression in the same animals because it circumvents the conventional requirement to sacrifice snakes to extract mRNA from dissected venom glands. Using qPCR and proteomic analysis, we show that gene expression and protein re-synthesis triggered by venom expulsion peaks between days 3-7 of the cycle of venom replenishment, with different protein families expressed in parallel. We demonstrate that venom re-synthesis occurs very rapidly following depletion of venom stores, presumably to ensure venomous snakes retain their ability to efficiently predate and remain defended from predators. The stability of mRNA in venom is biologically fascinating, and could significantly empower venom research by expanding opportunities to produce transcriptomes from historical venom stocks and rare or endangered venomous species, for new therapeutic, diagnostic and evolutionary studies.

First draft of the genomic organization of a PIII-SVMP gene

The evolutionary pathway of highly toxic proteins expressed in snake venom glands from proteins without toxic function and expressed in non-parotid tissues remains poorly understood. Here we examine gene structure of a representative of a venom protein with an ADAMs metalloproteinase evolutionary origin. The structure of the 15,652 bp Echis ocellatus pre-pro EOC00089-like PIII-SVMP gene was assembled from PCR-amplified sequences of overlapping genomic fragments. The gene comprises 12 exons interrupted by 11 introns. In a homology model of the EOC00089-like protein, the insertion of introns interrupting coding regions lie just after or between secondary structure elements. Long interspersed nuclear retroelements (LINE) L2/CR1 and RTE/Bov-B, short interspersed nuclear retroelements SINE/Sauria, and a hobo-activator DNA (Charlie, hAT) transposon were identified within introns 1, 3, 7 and 8. Pairwise amino acid sequence comparisons between EOC00089-like PIII-SVMP and its closest orthologs, ADAM28, from a mammal, Homo sapiens, and the lizard, Anolis carolinensis, showed that the ORFs of these three proteins share 42%/59%, 49%/69%, and 48%/65% (identity/similarity), respectively. The protein-coding positions interrupted by each of the 11 introns of the Echis PIII-SVMP gene are entirely conserved in the A. carolinensis and human ADAM28 genes. However, the lizard and the human ADAM28 genes contain 5 introns not present in the E. ocellatus gene. Furthermore, Echis and Anolis introns exhibit quantitatively and qualitatively distinctions in their inserted retroelements. These findings identify introns as possible key elements in the recruitment and amplification process of SVMPs into the venom gland of extant snakes. Ongoing reptile genome sequencing projects may shed light on this intriguing aspect of the emergence and evolution of venom toxin genes. Furthermore, the organization of the PIII-SVMP reported here provides a genomic explanation for the emergence of dimeric disintegrin subunits encoded by short messengers.

The structural and functional diversification of the Toxicofera reptile venom system

The evolutionary origin and diversification of the reptilian venom system is described. The resolution of higher-order molecular phylogenetics has clearly established that a venom system is ancestral to snakes. The diversification of the venom system within lizards is discussed, as is the role of venom delivery in the behavioural ecology of these taxa (particularly Varanus komodoensis). The more extensive diversification of the venom system in snakes is summarised, including its loss in some clades. Finally, we discuss the contentious issue of a definition for "venom", supporting an evolutionary definition that recognises the homology of both the venom delivery systems and the toxins themselves.

Research strategies to improve snakebite treatment: challenges and progress

Antivenom is an effective treatment of snakebite but, because of the complex interplay of fiscal, epidemiological, therapeutic efficacy and safety issues, the mortality of snakebite remains unacceptably high. Efficiently combating this high level of preventable death amongst the world's most disadvantaged communities requires the globally-coordinated action of multiple intervention programmes. This is the overall objective of the Global Snakebite Initiative. This paper describes the challenges facing the research community to develop snakebite treatments that are more efficacious, safe and affordable than current therapy.

Evolutionary analysis of novel serine proteases in the venom gland transcriptome of Bitis gabonica rhinoceros

BACKGROUND

Serine proteases are major components of viper venom and target various stages of the blood coagulation system in victims and prey. A better understanding of the diversity of serine proteases and other enzymes present in snake venom will help to understand how the complexity of snake venom has evolved and will aid the development of novel therapeutics for treating snake bites.

METHODOLOGY AND PRINCIPAL FINDINGS

Four serine protease-encoding genes from the venom gland transcriptome of Bitis gabonica rhinoceros were amplified and sequenced. Mass spectrometry suggests the four enzymes corresponding to these genes are present in the venom of B. g. rhinoceros. Two of the enzymes, rhinocerases 2 and 3 have substitutions to two of the serine protease catalytic triad residues and are thus unlikely to be catalytically active, though they may have evolved other toxic functions. The other two enzymes, rhinocerases 4 and 5, have classical serine protease catalytic triad residues and thus are likely to be catalytically active, however they have glycine rather than the more typical aspartic acid at the base of the primary specificity pocket (position 189). Based on a detailed analysis of these sequences we suggest that alternative splicing together with individual amino acid mutations may have been involved in their evolution. Changes within amino acid segments which were previously proposed to undergo accelerated change in venom serine proteases have also been observed.

CONCLUSIONS AND SIGNIFICANCE

Our study provides further insight into the diversity of serine protease isoforms present within snake venom and discusses their possible functions and how they may have evolved. These multiple serine protease isoforms with different substrate specificities may enhance the envenomation effects and help the snake to adapt to new habitats and diets. Our findings have potential for helping the future development of improved therapeutics for snake bites.

Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing

BACKGROUND

A long term research goal of venomics, of applied importance for improving current antivenom therapy, but also for drug discovery, is to understand the pharmacological potential of venoms. Individually or combined, proteomic and transcriptomic studies have demonstrated their feasibility to explore in depth the molecular diversity of venoms. In the absence of genome sequence, transcriptomes represent also valuable searchable databases for proteomic projects.

RESULTS

The venom gland transcriptomes of 8 Costa Rican taxa from 5 genera (Crotalus, Bothrops, Atropoides, Cerrophidion, and Bothriechis) of pitvipers were investigated using high-throughput 454 pyrosequencing. 100,394 out of 330,010 masked reads produced significant hits in the available databases. 5.165,220 nucleotides (8.27%) were masked by RepeatMasker, the vast majority of which corresponding to class I (retroelements) and class II (DNA transposons) mobile elements. BLAST hits included 79,991 matches to entries of the taxonomic suborder Serpentes, of which 62,433 displayed similarity to documented venom proteins. Strong discrepancies between the transcriptome-computed and the proteome-gathered toxin compositions were obvious at first sight. Although the reasons underlaying this discrepancy are elusive, since no clear trend within or between species is apparent, the data indicate that individual mRNA species may be translationally controlled in a species-dependent manner. The minimum number of genes from each toxin family transcribed into the venom gland transcriptome of each species was calculated from multiple alignments of reads matched to a full-length reference sequence of each toxin family. Reads encoding ORF regions of Kazal-type inhibitor-like proteins were uniquely found in Bothriechis schlegelii and B. lateralis transcriptomes, suggesting a genus-specific recruitment event during the early-Middle Miocene. A transcriptome-based cladogram supports the large divergence between A. mexicanus and A. picadoi, and a closer kinship between A. mexicanus and C. godmani.

CONCLUSIONS

Our comparative next-generation sequencing (NGS) analysis reveals taxon-specific trends governing the formulation of the venom arsenal. Knowledge of the venom proteome provides hints on the translation efficiency of toxin-coding transcripts, contributing thereby to a more accurate interpretation of the transcriptome. The application of NGS to the analysis of snake venom transcriptomes, may represent the tool for opening the door to systems venomics.

Ending the drought: new strategies for improving the flow of affordable, effective antivenoms in Asia and Africa

The development of snake antivenoms more than a century ago should have heralded effective treatment of the scourge of snakebite envenoming in impoverished, mostly rural populations around the world. That snakebite still exists today, as a widely untreated illness that maims, kills and terrifies men, women and children in vulnerable communities, is a cruel anachronism. Antivenom can be an effective, safe and affordable treatment for snakebites, but apathy, inaction and the politicisation of public health have marginalised both the problem (making snakebite arguably the most neglected of all neglected tropical diseases) and its solution. For lack of any coordinated approach, provision of antivenoms has been pushed off the public health agenda, leading to an incongruous decline in demand for these crucial antidotes, excused and fed by new priorities, an absence of epidemiological data, and a poor regulatory framework. These factors facilitated the infiltration of poor quality products that degrade user confidence and undermine legitimate producers. The result is that tens of thousands are denied an essential life-saving medicine, allowing a toll of human suffering that is a summation of many individual catastrophes. No strategy has been developed to address this problem and to overcome the intransigence and inaction responsible for the global tragedy of snakebite. Attempts to engage with the broader public health community through the World Health Organisation (WHO), GAVI, and other agencies have failed. Consequently, the toxinology community has taken on a leadership role in a new approach, the Global Snakebite Initiative, which seeks to mobilise the resources, skills and experience of scientists and clinicians for whom venoms, toxins, antivenoms, snakes and snakebites are already fields of interest. Proteomics is one such discipline, which has embraced the potential of using venoms in bio-discovery and systems biology. The fields of venomics and antivenomics have recently evolved from this discipline, offering fresh hope for the victims of snakebites by providing an exciting insight into the complexities, nature, fundamental properties and significance of venom constituents. Such a rational approach brings with it the potential to design new immunising mixtures from which to raise potent antivenoms with wider therapeutic ranges. This addresses a major practical limitation in antivenom use recognised since the beginning of the 20th century: the restriction of therapeutic effectiveness to the specific venom immunogen used in production. Antivenomic techniques enable the interactions between venoms and antivenoms to be examined in detail, and if combined with functional assays of specific activity and followed up by clinical trials of effectiveness and safety, can be powerful tools with which to evaluate the suitability of current and new antivenoms for meeting urgent regional needs. We propose two mechanisms through which the Global Snakebite Initiative might seek to end the antivenom drought in Africa and Asia: first by establishing a multidisciplinary, multicentre, international collaboration to evaluate currently available antivenoms against the venoms of medically important snakes from specific nations in Africa and Asia using a combination of proteomic, antivenomic and WHO-endorsed preclinical assessment protocols, to provide a validated evidence base for either recommending or rejecting individual products; and secondly by bringing the power of proteomics to bear on the design of new immunising mixtures to raise Pan-African and Pan-Asian polyvalent antivenoms of improved potency and quality. These products will be subject to rigorous clinical assessment. We propose radically to change the basis upon which antivenoms are produced and supplied for the developing world. Donor funding and strategic public health alliances will be sought to make it possible not only to sustain the financial viability of antivenom production partnerships, but also to ensure that patients are relieved of the costs of antivenom so that poverty is no longer a barrier to the treatment of this important, but grossly neglected public health emergency.

A high-throughput venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes

Despite causing considerable human mortality and morbidity, animal toxins represent a valuable source of pharmacologically active macromolecules, a unique system for studying molecular adaptation, and a powerful framework for examining structure-function relationships in proteins. Snake venoms are particularly useful in the latter regard as they consist primarily of a moderate number of proteins and peptides that have been found to belong to just a handful of protein families. As these proteins and peptides are produced in dedicated glands, transcriptome sequencing has proven to be an effective approach to identifying the expressed toxin genes. We generated a venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) using Roche 454 sequencing technology. In the current work, we focus on transcripts encoding toxins. We identified 40 unique toxin transcripts, 30 of which have full-length coding sequences, and 10 have only partial coding sequences. These toxins account for 24% of the total sequencing reads. We found toxins from 11 previously described families of snake-venom toxins and have discovered two putative, previously undescribed toxin classes. The most diverse and highly expressed toxin classes in the C. adamanteus venom-gland transcriptome are the serine proteinases, metalloproteinases, and C-type lectins. The serine proteinases are the most abundant class, accounting for 35% of the toxin sequencing reads. Metalloproteinases are the most diverse; 11 different forms have been identified. Using our sequences and those available in public databases, we detected positive selection in seven of the eight toxin families for which sufficient sequences were available for the analysis. We find that the vast majority of the genes that contribute directly to this vertebrate trait show evidence for a role for positive selection in their evolutionary history.

A rapid and sensitive fluorometric method for the quantitative analysis of snake venom metalloproteases and their inhibitors

Metalloproteases are responsible for the hemorrhagic effects of many snake venoms and contribute to other pathways that lead to local tissue damage. Methods that quantify snake venom metalloproteases (SVMP) are therefore valuable tools in research on the clinical, physiological, and biochemical effects of envenomation. Comparative analysis of individual, population, and species differences requires screening of large numbers of samples and treatments, and therefore require a method of quantifying SVMP activity that is simple, rapid, and sensitive. This paper demonstrates the properties of a new fluorometric assay of SVMP activity that can provide a measure of metalloprotease activity in 1 h. The assay is reliable, with variation among replicates sufficiently small to reliably detect differences in between species (F(19,60) = 2924, p < 0.001), even for those venoms with low overall activity. It is also sensitive enough to detect differences among venoms using <2 ng of whole venom protein. We provide an example use of this assay to detect the presence of natural SVMP inhibitors in minute samples of blood plasma from rock squirrels (S. variegatus), a natural prey species for North American rattlesnakes. We propose this assay is a useful addition to the set of tools used to characterize venoms, as well as high-throughput screening of natural or synthetic inhibitors, or other novel therapeutic agents against SVMP effects.

The venom composition of the parasitic wasp Chelonus inanitus resolved by combined expressed sequence tags analysis and proteomic approach

Background

Parasitic wasps constitute one of the largest group of venomous animals. Although some physiological effects of their venoms are well documented, relatively little is known at the molecular level on the protein composition of these secretions. To identify the majority of the venom proteins of the endoparasitoid wasp Chelonus inanitus (Hymenoptera: Braconidae), we have randomly sequenced 2111 expressed sequence tags (ESTs) from a cDNA library of venom gland. In parallel, proteins from pure venom were separated by gel electrophoresis and individually submitted to a nano-LC-MS/MS analysis allowing comparison of peptides and ESTs sequences.

Results

About 60% of sequenced ESTs encoded proteins whose presence in venom was attested by mass spectrometry. Most of the remaining ESTs corresponded to gene products likely involved in the transcriptional and translational machinery of venom gland cells. In addition, a small number of transcripts were found to encode proteins that share sequence similarity with well-known venom constituents of social hymenopteran species, such as hyaluronidase-like proteins and an Allergen-5 protein.

An overall number of 29 venom proteins could be identified through the combination of ESTs sequencing and proteomic analyses. The most highly redundant set of ESTs encoded a protein that shared sequence similarity with a venom protein of unknown function potentially specific of the Chelonus lineage. Venom components specific to C. inanitus included a C-type lectin domain containing protein, a chemosensory protein-like protein, a protein related to yellow-e3 and ten new proteins which shared no significant sequence similarity with known sequences. In addition, several venom proteins potentially able to interact with chitin were also identified including a chitinase, an imaginal disc growth factor-like protein and two putative mucin-like peritrophins.

Conclusions

The use of the combined approaches has allowed to discriminate between cellular and truly venom proteins. The venom of C. inanitus appears as a mixture of conserved venom components and of potentially lineage-specific proteins. These new molecular data enrich our knowledge on parasitoid venoms and more generally, might contribute to a better understanding of the evolution and functional diversity of venom proteins within Hymenoptera.

A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu)

Background

The genus Bothrops is widespread throughout Central and South America and is the principal cause of snakebite in these regions. Transcriptomic and proteomic studies have examined the venom composition of several species in this genus, but many others remain to be studied. In this work, we used a transcriptomic approach to examine the venom gland genes of Bothrops alternatus, a clinically important species found in southeastern and southern Brazil, Uruguay, northern Argentina and eastern Paraguay.

Results

A cDNA library of 5,350 expressed sequence tags (ESTs) was produced and assembled into 838 contigs and 4512 singletons. BLAST searches of relevant databases showed 30% hits and 70% no-hits, with toxin-related transcripts accounting for 23% and 78% of the total transcripts and hits, respectively. Gene ontology analysis identified non-toxin genes related to general metabolism, transcription and translation, processing and sorting, (polypeptide) degradation, structural functions and cell regulation. The major groups of toxin transcripts identified were metalloproteinases (81%), bradykinin-potentiating peptides/C-type natriuretic peptides (8.8%), phospholipases A₂ (5.6%), serine proteinases (1.9%) and C-type lectins (1.5%). Metalloproteinases were almost exclusively type PIII proteins, with few type PII and no type PI proteins. Phospholipases A₂ were essentially acidic; no basic PLA₂ were detected. Minor toxin transcripts were related to L-amino acid oxidase, cysteine-rich secretory proteins, dipeptidylpeptidase IV, hyaluronidase, three-finger toxins and ohanin. Two non-toxic proteins, thioredoxin and double-specificity phosphatase Dusp6, showed high sequence identity to similar proteins from other snakes. In addition to the above features, single-nucleotide polymorphisms, microsatellites, transposable elements and inverted repeats that could contribute to toxin diversity were observed.

Conclusions

Bothrops alternatus venom gland contains the major toxin classes described for other Bothrops venoms based on trancriptomic and proteomic studies. The predominance of type PIII metalloproteinases agrees with the well-known hemorrhagic activity of this venom, whereas the lower content of serine proteases and C-type lectins could contribute to less marked coagulopathy following envenoming by this species. The lack of basic PLA₂ agrees with the lower myotoxicity of this venom compared to other Bothrops species with these toxins. Together, these results contribute to our understanding of the physiopathology of envenoming by this species.