Open Access
This article is
- freely available
- re-usable
Toxins 2016, 8(10), 303; doi:10.3390/toxins8100303
Article
The Snake with the Scorpion’s Sting: Novel Three-Finger Toxin Sodium Channel Activators from the Venom of the Long-Glanded Blue Coral Snake (Calliophis bivirgatus)
1
Department of Pharmacology, Biomedicine Discovery Institute, Monash University, Clayton 3168, Australia
2
Venom Evolution Lab, School of Biological Sciences, University of Queensland, St. Lucia 4072, Australia
3
Institute for Molecular Bioscience, University of Queensland, St. Lucia 4072, Australia
4
Bejing Genomics Institute-Shenzhen, Shenzhen 518083, China
5
School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia 4072, Australia
6
Lee Kong Chian Natural History Museum, National University of Singapore, 2 Conservatory Drive, Singapore 117377, Singapore
7
Mtoxins, 1111 Washington ave, Oshkosh, WI 54901, USA
8
School of Pharmacy, University of Queensland, Woolloongabba 4102, Australia
*
Correspondence: Tel: +61-7-3346-2660 (I.V.); +61-4-0019-3182 (B.G.F.)
†
These authors contributed equally to this work.
Academic Editor:
Jan Tygat
Received: 15 September 2016 / Accepted: 10 October 2016 / Published: 18 October 2016
Abstract
:Millions of years of evolution have fine-tuned the ability of venom peptides to rapidly incapacitate both prey and potential predators. Toxicofera reptiles are characterized by serous-secreting mandibular or maxillary glands with heightened levels of protein expression. These glands are the core anatomical components of the toxicoferan venom system, which exists in myriad points along an evolutionary continuum. Neofunctionalisation of toxins is facilitated by positive selection at functional hotspots on the ancestral protein and venom proteins have undergone dynamic diversification in helodermatid and varanid lizards as well as advanced snakes. A spectacular point on the venom system continuum is the long-glanded blue coral snake (Calliophis bivirgatus), a specialist feeder that preys on fast moving, venomous snakes which have both a high likelihood of prey escape but also represent significant danger to the predator itself. The maxillary venom glands of C. bivirgatus extend one quarter of the snake’s body length and nestle within the rib cavity. Despite the snake’s notoriety its venom has remained largely unstudied. Here we show that the venom uniquely produces spastic paralysis, in contrast to the flaccid paralysis typically produced by neurotoxic snake venoms. The toxin responsible, which we have called calliotoxin (δ-elapitoxin-Cb1a), is a three-finger toxin (3FTx). Calliotoxin shifts the voltage-dependence of NaV1.4 activation to more hyperpolarised potentials, inhibits inactivation, and produces large ramp currents, consistent with its profound effects on contractile force in an isolated skeletal muscle preparation. Voltage-gated sodium channels (NaV) are a particularly attractive pharmacological target as they are involved in almost all physiological processes including action potential generation and conduction. Accordingly, venom peptides that interfere with NaV function provide a key defensive and predatory advantage to a range of invertebrate venomous species including cone snails, scorpions, spiders, and anemones. Enhanced activation or delayed inactivation of sodium channels by toxins is associated with the extremely rapid onset of tetanic/excitatory paralysis in envenomed prey animals. A strong selection pressure exists for the evolution of such toxins where there is a high chance of prey escape. However, despite their prevalence in other venomous species, toxins causing delay of sodium channel inhibition have never previously been described in vertebrate venoms. Here we show that NaV modulators, convergent with those of invertebrates, have evolved in the venom of the long-glanded coral snake. Calliotoxin represents a functionally novel class of 3FTx and a structurally novel class of NaV toxins that will provide significant insights into the pharmacology and physiology of NaV. The toxin represents a remarkable case of functional convergence between invertebrate and vertebrate venom systems in response to similar selection pressures. These results underscore the dynamic evolution of the Toxicofera reptile system and reinforces the value of using evolution as a roadmap for biodiscovery.
Keywords:
toxicofera; venom; evolution; neurotoxin; sodium channel; pharmacology1. Introduction
Squamate reptiles are a research area of great controversy and debate at all levels, including higher order organismal relationships [1]. Traditional frameworks based upon morphology [2,3,4,5,6] were dramatically over-turned by genetic studies which revealed that the morphological plasticity of these animals had obscured their authentic evolutionary relationships [7,8,9,10,11,12,13,14,15,16]. While the genetics are now well-resolved, a form of scientific culture-war persists, with some adherents of morphologically based taxonomy not modifying their position in the face of genetic evidence supporting the paradigm shifts [17]. Attempts at reconstructing the evolutionary history by using a combination of morphology, fossils, and genetics resulted in poorly resolved trees [18] that had significantly different frameworks from those reconstructed using only genetic evidence [7,8,9,10,11,12,13,14,15,16]. The utilisation of ambiguous morphological characters (which may be arbitrarily scored as divergent or convergent) creates errors relative to the use of non-ambiguous genetic data. Use of such skewed trees results in flawed conclusions when used to reconstruct the evolutionary history not only of the organisms but also their associated venom systems, such as in [17,19]. Another study which examined tissue pattern expressions of toxin and non-toxin homologues in an attempt to cast doubt on the shared ancestry of reptile venom systems was methodologically flawed due to expression values being calculated using replicate averages of values up to nearly 6,000-fold apart, only averaging N = 2 for many of the experiments, and including failed experiments with zero values in N = 3 calculations (Supplementary Tables S5–S9 of [20]).
Squamate reptiles are an extremely morphologically plastic clade, where dramatic variations in size and shape [21] may not be an accurate reflection of the organismal genetic relationships [22,23,24]. One particular clade (Toxicofera) evolved a morphological feature not otherwise described in extant squamates. The basic substrate for the evolution of reptile venom was the exaptation of existing salivary glands in the common ancestor of toxicoferan reptiles (Anguimorpha and Iguania lizards plus snakes), whereby the protein secreting region was increased in size and further segregated from the mucus secreting region [25]. These glands remained in an apparently incipient state in most Iguania except for vertebrate feeding lineages such as Anolis equestris, Pogona sp., and Chlamydosaurus Kingii. In snakes and anguimorph lizards, however, extensive diversification has occurred [19,25,26,27,28,29,30]. In particular, this key evolutionary innovation under-pinned the explosive diversification of the advanced snakes [10,11,12,13,14,15,27,30,31,32,33,34,35,36]. While initial venom proteins were pre-existing salivary proteins that were heightened in expression, such as CRiSP and kallikrein [27,37,38,39], the promiscuity venom glands share with other secretory tissues allowed for lineage-specific expression of other protein types [20,25,26,27,28,29,40,41,42,43,44,45,46,47,48] in addition to low-level reversal of expression of venom-gland-derived proteins in other tissues [40]. For many toxin types, particularly those with high disulphide-bond ratios relative to total amino acid content [37], evolutionary selection pressures resulted in neofunctionalisations that produced a myriad of derived activities [41].
As with many other modes of evolution, periodically the venom proteins undergo punctuated molecular evolution, dynamically diversifying into new functional forms early in their evolution, with periods of stasis after that interrupted by further periods of Gouldian punctuated evolution [49,50,51,52,53]. Examples include elapids first arriving in Australia when it was a continent largely devoid of snakes other than slow moving pythons and burrowing blind snakes [54]. Subsequent explosive organismal diversification was paralleled by dramatic shifts in venom chemistry even between populations with minimal organismal genetic variation [55,56].
Venom is a key evolutionary innovation, in myriad lineages, that diversifies in concert with organismal morphological diversification [11,13,14,15,25,26,27,28,29,30,33,34,37,41,43,44,45,46,47,55,56,57,58,59,60,61,62,63]. This mechanical and chemical interrelationship continues throughout the evolutionary history of a venomous lineage. In some cases, explosive morphological changes apparently occur in an evolutionarily short period of time, with the venom changing along with it in response to emergent selection pressures. Considerations of the venom delivery system must also include the morphological features associated with prey contact. Ambush feeding snakes such as most viperids rely on camouflage, which itself may be both morphological (pattern) and chemical (control of smell emission) [64]. They typically have muscular builds that facilitate fast and powerful strikes capable of driving their large, flexible fangs deep into a prey item and delivering copious amounts of venom from their bulky venom glands. Elapid snakes, however, are generally more agile as they are usually active foragers that deliver smaller amounts of very toxic venom through short, rigid fangs connected to the relatively small venom glands.
As proteinaceous toxins are large in size, they necessitate delivery via a wound, however small [26,27,58]. Thus, in addition to variation in the types and activities of proteins secreted by the venom glands, there has been extensive evolutionary tinkering of all elements of this integrated weapons system, including the glands themselves and the associated dentition. The morphological plasticity of the glands has resulted in diverse forms. Plesiomorphic anguimorph lizard glands are relatively unstructured despite containing large lumens in the well-defined protein-secreting region, but two independent lineages of anguimorph lizards have evolved completely segregated, membrane-encapsulated protein- and mucus-secreting regions [25,27,29]. In snakes there is considerable variation in oral gland form and function, while in the advanced snakes the maxillary venom glands have been fully segregated into distinct protein- and mucus-secreting tissues, accompanied by extensive variation of features such as the relative presence, size, and shape of liquid venom storage lumens [27,44,65].
While a venom delivering wound can be created by any form of tooth capable of penetrating skin, in the advanced snakes enlarged rear teeth have independently evolved on multiple occasions and also display extensive variation [44], in some cases accompanied by grooving that is convergent with that seen in other extant and extinct venomous lineages including archosauriforms [66], conodonts [67], sphenodons [68], insectivorous mammals such as shrews and solenodons [69,70,71,72], and bird-like dinosaurs hypothesised to specialise in feeding upon on early birds [73].
Refinements of the ancestral snake venom system have included, on three independent occasions, the evolution of high-pressure delivery systems terminating in hollow teeth, with the teeth in each case being derived from enlarged rear teeth from within three different non-front-fanged lineages [27,74]. Within each of these front-fanged lineages, on at least one occasion (twice within the elapid snakes) venom glands have become elongated to extend down into the neck region: Atractaspis within the Lamprophiidae; Calliophis and Toxicocalamus within the Elapidae; and Causus within the Viperidae [27]. While Atractaspis and Causus venom glands are elongations of the venom glands extending directly down the neck, Calliophis elongated glands differ sharply in that a long duct extends down the neck and then inside the ribs, with the glands located now inside the body cavity.
The influence of these elongations upon venom composition within each lineage has remained uninvestigated. For example, the sarafotoxin from Atractaspis venoms are an example of a lineage specific recruitment of protein type for use as a venom component, but it is unknown whether this precedes or is subsequent to the diversification of the venom glands into an apomorphic long-glanded state found only in one clade of Atractaspis, but not all. Similarly, Asian elapid snakes in the Calliophis genus have plesiomorphic short-glanded (C. gracilis) and apomorphic long-glanded forms (C. bivirgatus and C. intestinalis).
Coral snakes in the genus Calliophis feed upon other snakes, including other snake-eating venomous species of Elapidae such as kraits (Bungarus) and king cobras (Ophiophagus) [75,76]. A unique evolutionary scenario ensues, a chemical arms race between predator and prey in which the risk of role reversal becomes a key selection pressure driving the evolution of toxins that rapidly render prey incapable of retaliation or escape. Snakes that hunt animals capable of inflicting serious retaliatory wounds often release their intended prey after envenomation. In this situation, selection may favour the evolution of toxins that rapidly disable prey, either to prevent it moving too far to be recovered or to prevent the possibility of it attacking and injuring the snake.
With its combination of electric blue dorsolateral stripes and neon red head, tail, and ventral scales, the blue coral snake, Calliophis bivirgatus, is arguably one of the world’s most striking species of snake (Figure 1a). An encounter with one is high on the list for many reptile enthusiasts and nature photographers visiting southern Thailand, Malaysia, Singapore, and western Indonesia. The species is of additional interest to anatomists and toxinologists studying the evolution and diversification of the snake venom system as it (along with its congener C. intestinalis) possesses novel elongated venom glands that extend for up to one quarter of the length of its body [77] (Figure 1b). It is also of medical significance as, in spite of only a small handful of confirmed bites, it has been responsible for at least one human fatality [76], is suspected of causing at least one more [78], and has no known antivenom. In spite of these high levels of interest, the venom has been subject to relatively few studies [79,80,81]. Those studies that examined the toxin content of the venom concluded that all the three-finger toxins present were exclusively cytotoxic in their effects [80,81]. However, this limited scope of activity attributed to the venom was reflective of the very narrow scope of assays performed and cytotoxicity was largely assumed based on similarity of partial sequences to other toxin types from other snakes rather than full activity characterisation. One study, which examined the usefulness of Taiwan antivenom, preincubated the venom with antivenom (a clinically unrealistic situation) and even then required very high doses to exert any meaningful level of inhibition [81].
2. Results and Discussion
In the present study, the pharmacology of blue coral snake venom was investigated. A form of neurotoxicity, previously known from cone snail and scorpion venoms, was identified for the first time from the venom of a snake. In the indirectly stimulated chick biventer cervicis nerve-muscle assay, which is a skeletal muscle preparation, C. bivirgatus venom (10 μg/mL) produced large muscle contractions and fasciculations (Figure 1c) which were significantly inhibited by the addition of the sodium channel antagonist tetrodotoxin (Figure 1d; TTX; 0.1 μM) (Figure 1c). Activity-guided fractionation using SH-SY5Y human neuroblastoma cells identified a peak dominated by a single isotopic mass of 6725.89 Da as the active component (Figure 1d–g). Consistent with the effects of crude venom in the chick biventer cervicis nerve-muscle assay, toxin-induced responses in neuroblastoma cells were abolished by TTX (Control response, 3.12 ± 0.04 AFU; TTX (1 μM), 0.12 ± 0.03 AFU), suggesting direct effects on voltage-gated sodium channels (Figure 1e). The amino acid sequence of the active component, which we called calliotoxin or δ-elapitoxin-Cb1a, was determined as: LECYDTIFKWHTMTCPEGQNLCFYYFTWRIFLVRGCTATCPVGYSHTHCCDTDKCNN using a combination of Edman degradation sequencing and venom gland transcriptome analysis (Figure 2).
The sequence calculated monotopic molecular weight is 6725.91 Da, which is in accord with the mass spectrometry monoisotopic molecular weight of 6725.9 Da. Calliotoxin belongs to the class of three-finger toxins (3FTx) but has low sequence homology to other known toxins in the family (Figure 2).
3FTxs have diversified into forms with novel pharmacology that include post-synaptic (muscarinic and nicotinic receptors), synaptic (acetylcholinesterase inhibitors), and presynaptic (L-type calcium channels) activities [61,82]. However, none to date have been reported with activity at voltage-gated sodium channels. Therefore, we further examined the pharmacological effect of calliotoxin on HEK293 cells heterologously expressing NaV1.4, the NaV subtype essential for skeletal muscle function in mammals, using whole-cell patch-clamp recordings. Calliotoxin enhanced peak inward current and delayed inactivation (Figure 3a,b), causing a small but significant hyperpolarizing shift in the V1/2 of activation (control, −35.0 ± 0.1 mV; calliotoxin, −38.1 ± 0.3 mV; p < 0.0001) and a significant depolarizing shift in the V1/2 of fast inactivation (control, −66.1 ± 0.5 mV; calliotoxin, −61.6 ± 0.7 mV; p < 0.0001) (Figure 3c). Calliotoxin voltage-dependently delayed the inactivation time constant (−20 mV pulse: control, 0.24 ± 0.02 ms; calliotoxin, 5.22 ± 0.9 ms; p < 0.05; Figure 3d), which resulted in a persistent current (−20 mV pulse: control, −6 ± 4 pA; calliotoxin, −536 ± 33 pA; p = 0.015; Figure 3e). In line with the above NaV activator activity, calliotoxin also significantly increased inward ramp currents (control, −143 ± 12 pA; calliotoxin, −728 ± 52 pA; p = 0.0004; Figure 3f,g).
These results demonstrate that calliotoxin acts directly at NaV in a manner reminiscent of other NaV activators from the venom of cone snails, scorpions, spiders, wasps, and anemones [26,83,84,85,86,87,88,89,90]. The potency is similar to that characterised for Australian hexathelid spiders such as funnel webs (Atrax and Hadronyche species) and the eastern mouse spider, which produce similar fasciculations in the neuromuscular organ bath assay (Missulena bradleyi) [91,92] (Figure 1c). Thus, calliotoxin is the first identified NaV activator from snake venom and represents a structurally novel class of NaV gating modifiers. While crotamine peptides from South American Crotalus species (rattlesnakes) had been previously thought to interact with sodium channels [93,94,95], in addition to being potent myotoxins, it has been recently observed that the neurotoxic activity is guided by potassium channel interactions [96,97,98,99]. A PLA2 toxin previously isolated from the Asian pit-viper Gloydius ussurensis showed voltage-dependence of activation of sodium channels in sensory neurons [100]. However, in contrast to calliotoxin, this toxin caused predominant effects on the voltage-dependence of activation and decreased peak current. In addition, this toxin also affected voltage-gated potassium channels. It remains to be determined whether calliotoxin also has effects on other members of the voltage-gated ion channel family.
There is a general evolutionary trend for an inverse relationship between mechanical forms of prey subjugation and chemical forms. In some predatory lineages, once a mechanical form of prey capture has evolved, it is co-opted for use as a delivery system for venom. Spiders are one such example, in which the ancestral mechanical form of predation using chelicerae to inflict fatal wounds in prey items resulted in a selection pressure for the derived use as a venom delivery system, with a subsequent reduction in chelicerae size [101,102]. Similarly, the elongated dentition used by cleaner wrasse mimics from the Plagiotremus fangblenny genus for their unique parasitic feeding strategy was co-opted for venom delivery in the Meiacanthus genus of fangblenny [103,104]. In other cases, a structural feature has been derived to become a venom-delivering apparatus, such as the fin-supporting spines in fish [105]. In early toxicoferan venomous reptiles, the pre-existing teeth were sufficient for the generation of a wound, allowing for low-pressure venom delivery via chewing, with the array of extant dentition types subsequently evolving as refinements [27,58]. An alternate scenario is one in which a hypertrophic mechanical structure is atrophied subsequent to the evolution of a different part of the anatomy to serve as a venom-delivering structure. Examples include the large claws of plesiomorphic scorpions accompanied by small stingers (telsons), relative to the gracile claws of apomorphic scorpions which have large telsons [106], and octopus beaks and glands, with species having large beaks having smaller glands than those which have small beaks [107,108].
The primary function of predatory venoms is the subjugation of prey items rather than lethality [27]. From an evolutionary and practical perspective, there is little difference between an unconscious or completely immobilized prey item and a dead one as the helpless former would simply suffocate in the stomach after being swallowed. Thus, the primary shaping pressure in predatory venom evolution is rapid prey subjugation rather than rapid lethality [1].
Chemical forms of prey subjugation allow for a decoupling of the physical interaction between predator and prey. Snakes that hunt animals capable of inflicting serious retaliatory wounds often release their intended prey after envenomation. In this situation, selection may favour the evolution of toxins that rapidly disable prey, either to prevent it from moving too far to be recovered or to prevent the possibility of it attacking and injuring the snake. For example, Northern Pacific rattlesnakes (Crotalus oreganus) and inland taipans (Oxyuranus microlepidotus) both feed on rodents, which are capable of inflicting life-threatening bite wounds upon snakes. Both species release their prey after biting. Crotamine, a neurotoxin in the venom of the rattlesnake, paralyses the hind limbs of the rodent in seconds, ensuring it is unable to travel far from the heavy bodied ambush-hunting snake in the approximately two minutes it takes for it to die [109,110]. Inland taipans, which are fast moving active hunters, face their rodent prey in confined underground spaces and thus are at considerable risk of sustaining retaliatory injuries. Their venom is the most toxic (to rodents) of any snake [111] and contains an exceptionally fast acting fXa:fVa coagulation enzyme complex, which causes small clots to form throughout the rat’s circulatory system resulting in rapid knockdown via stroke injury [27,112,113,114]. In human bite victims the same amount of venom is diluted into a much larger blood volume, producing countless microthrombi which by themselves are too small to cause a stroke but instead consume all of the clotting factors, with death the result of internal bleeding such as cerebral hemorrhage [115].
Rapid knockdown effects may also evolve in situations where the cost of prey recovery is particularly high due to the distance that might be traversed by dying prey. Mambas (genus Dendroaspis) are arboreal elapid snakes that feed on birds, which they subdue with venom containing unique kunitz peptide neurotoxins that rapidly cause excitatory paralysis through the inhibition of voltage-dependent potassium channels, resulting in the sustained release of acetylcholine, while 3FTx inhibit acetylcholinesterase, further elongating the acetylcholine action [27,82,116,117,118,119,120]. Another species of arboreal elapid snake, the Stephens’ banded snake (Hoplocephalus stephensi), uses natriuretic peptide toxins to cause a precipitous drop in the blood pressure of prey animals [121]. In the case of cone snails, the venom was first evolved for worm-hunting, with the subsequent evolution of fish-specific venom for defence against predators, which was then co-opted on at least two occasions for fish-specialist prey preference [122,123,124]. Piscivorus cone snails also risk being unable to recover their prey if the fish are able to escape after envenomation—even if the fish only swims a small distance before dying, the slow-moving snail faces a high risk of losing its meal to an opportune scavenger. The risk of going hungry thus forms the selection pressure driving the evolution of the devastating chemical weaponry of venoms.
NaV are widely expressed across animal phyla including in the invertebrate and vertebrate nervous systems and are critically important for neuromuscular action potential conduction. Accordingly, toxins that target NaV are widely found in venoms, with known families including the µ-conotoxins, µ-theraphotoxins, δ-conotoxins, α- and β- scorpion toxins, and δ-theraphotoxins [26,83,84,85,86,87,88,89,90]. Although the mechanism of action of these toxins is distinct, all impair normal NaV function, which critically depends on voltage-dependent activation and inactivation to allow regulated Na+ influx. As a consequence, these toxins lead to catastrophic disruption of neuronal or neuromuscular physiology. Although it initially seems surprising that a snake toxin would converge on the same activity as that of toxins from cone snails, spiders, wasps, and scorpions, convergence in molecular targeting is actually common amongst venomous organisms [26].
On a deeper theoretical level, these results have relevance to the debate regarding whether evolution is contingent or predictable (e.g., [125,126]). In fact, as Daniel Dennett has pointed out [127], it is both. Targeting NaV is, in Dennett’s terminology, a “Good Trick”. Since there is likely to be a limited number of good tricks available to venomous organisms with the need to rapidly render prey incapable of escape or retaliation, it is predictable that convergence will occur [1]. The possibilities for evolutionary innovation, however, are constrained by historical contingency [128]. It is historical contingency that makes it more likely (thus more predictable) that the novel activity observed in the present study should emerge within the 3FTx, a toxin class typically under strong positive selection within the Elapidae [129]. Thus we see a predictable convergence of molecular targeting arising (predictably) within a toxin class “selected” by historical contingency—evolution is not predictable rather than contingent, it is predictable because it is contingent.
In this case, the selection pressure driving this contingent and predictable evolutionary trajectory is involvement with a prey animal that relies upon extreme coordination in its fight or flight response. Cone snails have such an interaction with fish, mamba snakes with birds, and long-glanded blue coral snakes with fast moving snakes, some of which are capable of venomous retaliatory actions potentially lethal to the coral snakes. Thus, a novel, extreme selection pressure has driven the venom evolution in a new direction. The hallmark of other elapid snake venoms are α-neurotoxins that antagonistically bind to the post-synaptic nicotinic acetylcholine receptors to produce flaccid paralysis, as is also the case in other venomous snakes lineages including gracile colubrid snakes with enlarged rear-fangs [19,27,44,47,59,62,130,131,132,133,134,135,136]. Potent α-neurotoxicity preceded the evolution of the high-pressure front-fanged system of elapids, as was revealed by our discovery of α-colubritoxin, the first 3FTx isolated and characterized from a non-front-fanged lineage [59]. The toxin was of the same form as had been studied in elapid snakes previously, where they had been called ‘weak neurotoxins’ since they were only weakly potent on mice [137,138,139,140,141,142,143,144,145,146,147,148,149]. A problem that had also confounded studies of Boiga irregularis (brown tree snake), which had thus concluded they were non-venomous based on their effects in a murine model [150,151,152,153,154,155]. However, all such plesiotypic 3FTx rich venoms were later shown to be much more potent in a taxon-specific manner on diapsid (bird/reptile) than on synapsid (mammal) post-synaptic nicotinic acetylcholine receptors, with some non-front-fanged snakes being as potent as elapid snakes in such a taxon specific manner and thus were not ‘weak’ against natural prey items [59,82,130,156,157,158,159,160].
In contrast, C. bivirgatus venom produces a new action with a net effect that is diametrically opposed to the typical elapid snake mode of action: spastic paralysis that is potent on both the avian and murine assays models used in this study. The early evolving elapid snake neurotoxic effect is one of respiratory failure from a paralysed diaphragm because the snake toxins prevent acetylcholine facilitating the diaphragm’s contraction, thus leaving it stuck in the non-contracted resting state. In contrast, C. bivirgatus venom keeps sodium channel toxins open, blocking the closing of the channel that would terminate the nerve transmission and allow the muscle to go back to the resting state. In this case, the net effect is that muscles are stuck in the contracted activated state instead of the typical elapid effect of being stuck in the non-contracted resting state. This is convergent with the rigid paralysis cone snails produce when feeding upon fish or the muscle fasciculations characteristic of mamba envenomations, both of which have been shaped by prey with high escape potential. The unique selection pressure operating with the niche occupied by C. bivirgatus has resulted in the flipping of the neurotoxic polarity, with this novel effect evolving due to the selection pressure of fast moving prey with high escape potential but also the potential for lethal retaliatory actions.
While the speed of action, and thus prey subjugation, are rapid, the question of relative toxicity must still be considered. Why else would sodium channel delay of inactivation toxins not be more common and only selected for in such extreme circumstances? Why has this only emerged once in snakes? It very well may be that while the α-neurotoxins are slower to take effect, they may be ultimately more effective as a consequence of their ability for sustained effect due to the pathology being produced due to simple steric inhibition. This is reflective in C. bivirgatus having a murine intravenous LD50 of 0.7–0.8 mg/kg [81] in comparison to that of a venom almost exclusively containing classic α-neurotoxins such as Acanthophis antarcticus which has a murine intravenous LD50 of 0.25 mg/kg [159,160]. Consistent with this, the black mamba (Dendroaspis polylepis), which consume a higher percentage of rodents than birds in their diet, has a murine intravenous LD50 of 0.5 mg/kg, while the green mambas, which have a higher percentage of birds in their diets, and thus would benefit from less toxic but more rapidly acting knock-down venoms, have LD50s of over 2.5 mg/kg for the Eastern green mamba (Dendroaspis angusticeps), 0.9 mg/kg for the Jameson’s mamba (Dendroaspis jamesoni), and 1 mg/kg for the Western green mamba (Dendroaspis viridis) [161]. Such tests reinforce that α-neurotoxicity is likely to be selected for due to its sustained paralytic effect, which is ultimately lethal due to its persistence; whereby complete immobilization also results in death but with the former being the outcome selected for by the shaping evolutionary pressures as there is no functional difference between a prey in suspended animation and one that is deceased.
It is interesting to note that unlike plesiotypic 3FTx which retain all 10 ancestral cysteines and are more potent to avian/reptilian post-synaptic nicotinic acetylcholine receptors than those of mammals [59,61,82,130,132,134,135], the C. bivirgatus plesiotypic presynaptic sodium channel toxins do not display taxon specific effects, as the crude venom and pure toxins were potent on both avian neuromuscular organ bath preparations and mammalian ion-channels FLIPR assays (Figure 1 and Figure 3) . Calliotoxin lacks the second and third plesiotypic cysteines (Figure 2), similar to that of the classic Type I (aka: short-chain) α-neurotoxins, which also do not display a taxon specific effect [61,160]. This is consistent with the fact that while there are very few human C. bivirgatus envenomations on record, the majority were lethal. Thus, while the venom may be three times less potent than that of α-neurotoxic snakes, this is off-set by the massive venom yields produce by the elongated venom glands, with dry weight single venom extraction quantities reaching 150 mg in this study for large (>1.2 m) specimens. This massive venom yield proportional to the length of its slender body is produced by the innovatively slender venom glands connected to the fang by a venom duct that extends until the venom gland is inside the rib cavity (Figure 1b) with the venom gland compressing musculature now wrapped completely around the gland, squeezing it from the back third in a manner analogous to that of the cone snail venom system. Thus the venom diversification in this lineage is paralleled by a morphological uniqueness, the evolution of which may have been driven by the same selection pressures as that of the toxins themselves. Studies between this species and the short-glanded C. gracilis and the long-glanded C. intestinalis will be revealing in regards to timing the venom gland elongation versus the timing of the molecular evolution of this extremely novel new diversification of the 3FTx framework. Regardless of the relative timing, the rapid interplay underscores that reptiles as a general character are extremely morphologically plastic and that their venom systems also dynamically evolve at all trophic levels.
A crucial difference between the two paralytic strategies (inhibitory/flaccid/limp or excitatory/contracting/spastic) rests upon how long the effect can be sustained. In the case of nerve transmission, certain ions move along an ion-gradient whenever a particular channel is opened at any number of steps in the transmission of a nerve impulse. Over-stimulation has inbuilt rate-limiting steps as the ion-pumps can only work at a finite rate to reset the ion-gradient. Thus, keeping the sodium channel open by delaying its inactivation, the long-glanded blue coral snake toxins will sustain this action only so long as there is a sufficient sodium gradient. Once isotonic levels for this ion are reached, there will be no more spontaneous action until the ion pumps re-establish a sufficient gradient. Thus, the delay of presynaptic sodium channel inactivation is a resource which may be exhausted. In contrast, α-neurotoxins bind to and block the post-synaptic nicotinic acetylcholine receptors; an action that may be sustained for a much longer period as it does not require the involvement of other molecules or chemical gradients. It is a case of antagonistic binding to the receptors, with fundamental interactions such as electrostatic charge or hydrophobic attractions determining the interaction strength. Once bound to the receptor, some toxins are virtually irreversible. Such immobilizing paralysis would subjugate the prey items quite efficiently for an indefinite period, thus having a marked degree of lethality. Therefore, there would exist a strong selection pressure to maintain such a predatory weapon. This underscores just how radical a change in functionality has occurred in C. bivirgatus venom.
The shift also parallels a shift in the metabolic state of the prey. Many Australian elapids snakes feed on reptilian prey that is inactive/dormant at the time of the snake’s attack, and thus suffocation through lethal lethargy inducing flaccid paralysis is an effective predatory strategy in elapid snakes [56]. Such a predatory strategy is also widely used in a myriad of venomous snakes that lack the apomorphic high-pressure delivery systems [44,47,59,62,130,131,132,134,135,136]. In contrast, the long-glanded blue coral snake feeds on other snakes that are at their highest activity level. Thus, the rapidity of effect becomes the primary selection pressure [1].
Due to its extraordinary venom glands and extremely novel venom chemistry, this enigmatic species may be considered as the epitome of toxicoferan reptile venom system derivation. Study of this species increases our knowledge of venom evolution and demonstrates the relevance of studying toxins, proteins unconstrained by endophysiological functions and under extreme selection pressures [56], in order to gain insight into more general patterns of evolution. In addition, the results are of considerable interest to those seeking to understand the mechanisms underlying NaV gating. These results underscore the dynamic evolution of the Toxicofera reptile system at all trophic levels and reinforces the value of using evolution as a roadmap for biodiscovery.
3. Materials and Methods
3.1. Materials
Specimens of the long-glanded blue coral snake (Calliophis bivirgatus) were captive animals of Malaysia stock. All other reagents were from Sigma Aldrich (Castle Hill, Sydney, NSW, Australia) unless otherwise specified.
3.2. Animal Ethics
All animal experiments used in this study were approved by the SOBS-B Monash University Animal Ethics Committee MARP/2014/97 (1 December 2014).
3.3. Neurotoxicity Studies
Male chicks (4–10 days) were euthanised by CO2 and exsanguination. Both chick biventer cervicis nerve muscle preparations were isolated and mounted on wire tissue holders under 1 g resting tension in 5 mL organ baths containing physiological salt solution (NaCl, 118.4 mM; KCl, 4.7 mM; MgSO4, 1.2 mM; KH2PO4, 1.2 mM; CaCl2, 2.5 mM; NaHCO3, and 25 mM glucose, 11.1 mM), maintained at 34 °C and bubbled with 95% O2/5% CO2. Indirect twitches were evoked by electrical stimulation of the motor nerve (supramaximal voltage, 0.2 ms, 0.1 Hz) using a Grass S88 stimulator (Grass Instruments, Quincy, MA, USA). d-Tubocurarine (10 μM) was added, and subsequent abolition of twitches confirmed selective stimulation of the motor nerve, after which thorough washing with physiological salt solution was applied to re-establish twitches. The preparation was equilibrated for 30 min before the addition of venom, which was left in contact with the preparation for a maximum of 3 h to test for slow developing effects. Efficacy of tetrodotoxin (TTX; 0.1 μM) was assessed via a 10 min pre-incubation in the organ bath.
3.3.1. Assay-Guided Fractionation
Crude Calliophis bivirgatus venom (150 µg) was fractionated into thirty second fractions on a ThermoScientific Hypersil (Scoresby, Vic, Australia) BDS C18 column (200 mm, 5 µm) using a linear gradient from 5% to 70% solvent B over 60 min (solvent A, H2O/0.05% TFA; solvent B, 90% acetonitrile/0.05% TFA) with a flow rate of 1 mL/min. Orthogonal purification of active fractions was carried out using a Phenomenex Synergi Hydro-RP 80A HPLC column (Phenomenex, Torrance, CA, USA) (250 × 4.6 mm, 4 µm) with flow rate of 0.7 mL/min and linear gradient from 30% to 50% B over 40 min (solvent A, H2O/0.05% TFA; solvent B, 90% acetonitrile/0.05% TFA). Activity testing was conducted as previously described [162] in human neuroblastoma SH-SY5Y cells (ECACC, Salisbury, Wiltshire, UK) maintained at 37 °C/5% CO2 in Roswell Park Memorial Institute (RPMI) medium containing 15% foetal bovine serum (FBS) and 2 mM l-glutamine. Cells were loaded with Calcium 4 No-Wash dye (Molecular Devices, Sunnvale, CA, USA) in 384-well black-walled imaging plates (Sigma Aldrich, Castle Hill, NSW, Australia) 48 h after plating at a density of 50,000 cells/well and incubated for 30 min at 37 °C. Freeze-dried fractions were resuspended in physiological salt solution (PSS; composition: 140 mM NaCl, 11.5 mM glucose, 5.9 mM KCl, 1.4 mM MgCl2, 1.2 mM NaH2PO4, 5 mM NaHCO3, 1.8 mM CaCl2, 10 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), pH 7.4) and added to dye-loaded cells using a FLIPRTETRA (Molecular Devices, Sunnyvale, CA, USA) plate reader while measuring real-time fluorescence responses (excitation, 470–495 nm; emission, 515–575 nm). Selective activation of endogenously expressed NaV channels was confirmed through inhibition of calliotoxin-induced responses by TTX (1 μM) and lack of effect of the nAChR antagonist d-tubocurarine (10 μM), the M1 mAChR antagonist pirenzepine (100 μM) and the α1-adrenoreceptor antagonist prazosin (10 μM) (Figure 3h).
3.3.2. Mass Spectrometry
Samples were separated using reversed-phase chromatography on a Dionex Ultimate 3000 RSLC nano-system. Using a flow rate of 30 µL/min, samples were desalted on a Thermo PepMap 100 C18 trap (0.3 × 5 mm, 5 µm) for 5 min, followed by separation on a Acclaim PepMap RSLC C18 (150 mm × 75 µm) column at a flow rate of 300 nL/min with a gradient of 10%–95% buffer B over 60 min where buffer A = 1% ACN/0.1% FA and buffer B = 80% ACN/0.1% FA. Eluted molecules were directly analysed on an Orbitap Elite mass spectrometer (Thermo, Brisbane, Australia) using an NSI electrospray interface. Source parameters included a capillary temperature of 275 °C; S-Lens RF level at 60%; source voltage of 2 kV and maximum injection times of 200 ms for MS. Data were deconvoluted using Protein Deconvolution software (Thermo).
3.3.3. Edman Degradation and Venom Gland Transcriptomics
Edman degradation was carried out by the Australian Proteome Analysis Facility (APAF, Sydney, Australia). Purified native calliotoxin was solubilised in ammonium bicarbonate (25 mM)/10% ACN and reduced using DTT (25 mM) at 56 °C for 0.5 h, followed by alkylation using iodoacetamide (55 mM) at room temperature for 0.5 h. The reaction mix was then desalted/purified by RP-HPLC using a Zorbax 300SB-C18 column (3 × 150 mm, Agilent, Santa Clara, CA, USA). The volume was reduced under vacuum and loaded onto a precycled, Biobrene-treated disc and subjected to 60 cycles of Edman N-terminal sequencing using an Applied Biosystems 494 Procise Protein Sequencing System (Applied Biosystems, Foster City, CA, USA), resulting in unambiguous identification of 47 amino acid residues. Venom gland transcriptomics were conducted by the IMB Sequencing Facility (Institute for Molecular Bioscience, The University of Queensland, St Lucia, Qld, Australia). Libraries were prepared with the TruSeq Stranded mRNA kit (Illumina, San Diego, CA, USA), and were sequenced on the Illumina NextSeq (Illumina) 500 using 2 × 150 bp reads and V2 chemistry. To identify the full sequence of δ-elapitoxin-Cb1a, forward and reverse sequences were merged using MacQIIME (Werner Lab, SUNY Cortland, NY, USA) join_paired_end.py and matched to the sequence determined by Edman degradation using standalone BLAST.
3.3.4. Electrophysiology
HEK-293 cells stably expressing hNaV1.4 (SB Drug Discovery, Glasgow, UK) were cultured in MEM containing 10% v/v foetal bovine serum, supplemented with l-glutamine 2 mM and selection antibiotics, as recommended by the manufacturer. Cells were grown in a humidified 5% CO2 incubator at 37 °C, grown to 70%–80% confluence, and passaged every three to four days using TrypLE Express (Invitrogen, Scoresby, VIC, Australia). For electrophysiology experiments, cells were dissociated by incubating with Detachin (Bio-Scientific, Kirrawee, NSW, Australia) at 37 °C for 5 min, then resuspended in Ex-Cell ACF CHO Medium with 25 mM HEPES (Sigma-Aldrich, Castle Hill, NSW, Australia) and allowed to recover with stirring for 30 min.
Whole-cell patch-clamp experiments were performed on a QPatch-16 automated electrophysiology platform (Sophion Bioscience, Ballerup, Denmark) using 16-channel planar patch chip plates (QPlates; Sophion Bioscience, Ballerup, Denmark) with a patch hole diameter of 1 µm and resistance of 2 ± 0.03 MΩ. Cell positioning and sealing parameters were set as follows: positioning pressure −60 mbar, minimum seal resistance 0.1 GΩ, holding potential −100 mV, holding pressure −20 mbar. Whole-cell currents were filtered at 5 kHz and acquired at 25 kHz.
The extracellular solution contained in mM: NaCl 145, KCl 4, CaCl2 2, MgCl2 1, HEPES 10, and glucose 10; pH 7.4; osmolarity 305 mOsm. The intracellular solution contained in mM: CsF 140, EGTA/CsOH 1/5, HEPES 10 and NaCl 10; pH 7.3; osmolarity 320 mOsm. Purified native calliotoxin was diluted in extracellular solution to a concentration of 200 nM and incubated for 5 min. The effects of calliotoxin were compared to pre-toxin control parameters in the same cell.
Current (I)-voltage (V) curves were obtained with a holding potential of −80 mV followed by a pre-pulse of −100 mV for 50 ms and a series of 50 ms step pulses that ranged from −80 to +60 mV in 5 mV increments before returning to a holding potential of −80 mV (repetition interval 5 s). Conductance-voltage curves were obtained by calculating the conductance (G) at each voltage (V) using the equation G = I/(V − Vrev), where Vrev is the reversal potential, and fitted with a Boltzmann equation: GNa = GNa,max/1 + exp[(Vm − V1/2)/k], where GNa is the voltage-dependent sodium conductance, GNa,max is the maximal sodium conductance, V1/2 is the potential at which activation is half-maximal, Vm is the membrane potential, and k is the slope factor. Fast inactivation time constants were calculated by fitting current decay traces obtained from the above I-V protocol with a single exponential function, and persistent current was determined as the average current 40–50 ms after pulse onset.
Voltage dependence of steady-state fast inactivation was measured using a series of 500 ms pre-pulses, ranging from −120 to −10 mV in 10 mV increments, followed by a 20 ms pulse of −20 mV to assess the available non-inactivated channels (repetition interval 30 s). Peak inward currents (I) were normalized to the maximal inward current (Imax) and fitted using a Boltzmann equation: I/Imax = 1/(1 + exp[(Vm − V1/2)/k)], where Imax is the maximal inward current, V1/2 is the half-maximal sodium current, Vm is the pre-conditioning pulse potential, and k is the slope factor.
Ramp currents were evoked by a depolarization from a holding potential of −100 to +20 mV at a rate of 2.4 mV/ms.
3.4. Data Analysis and Statistics
Unless otherwise stated, all data are expressed as the mean ± standard error of the mean (SEM) determined from at least n = 3 replicates. FLIPRTetra results were converted to response over baseline using ScreenWorks 3.2.0.14 (Molecular Devices, Sunnyvale, CA, USA) and plotted using GraphPad Prism 6 (GraphPad Software, San Diego, CA, USA). Twitch heights were measured from the baseline in two minute intervals. Responses were expressed as a percentage of twitch height prior to the addition of the peptide. Contractile responses to agonists obtained at the conclusion of the experiment were measured and expressed as a percentage of the response obtained prior to the addition of peptide. Where indicated, a one-way analysis of variance (ANOVA) followed by a Bonferroni-corrected post-hoc test was used to determine statistical significance of responses. Statistical analysis was performed using the GraphPad Prism 5.
Acknowledgments
Daryl C. Yang was supported by a Monash University postgraduate scholarship. Daniel Dashevsky, Timothy N. W. Jackson and Ivan Koludarov were supported by University of Queensland postgraduate scholarships. Irina Vetter and Bryan G. Fry were supported by the Australian Research Council.
Author Contributions
Bryan G. Fry, Iwan Hendrikx and Irina Vetter conceived and designed the experiments; Daryl C. Yang, Jennifer R. Deuis, Daniel Dashevsky, James Dobson, Timothy N. W. Jackson, Andreas Brust, Bing Xie, Ivan Koludarov, Jordan Debono, Iwan Hendrikx, Wayne C. Hodgson, Peter Josh, Amanda Nouwens, Gregory J. Baillie, Timothy J. C. Bruxner, Paul F. Alewood, Kelvin Kok Peng Lim, Nathaniel Frank, Irina Vetter, Bryan G. Fry performed the experiments, analysed the data and wrote the paper; Nathaniel Frank contributed essential reagents/materials/analysis tools.
Conflicts of Interest
The authors declare no conflict of interest and the founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.
References
- Jackson, T.N.; Fry, B.G. A Tricky Trait: Applying the Fruits of the “Function Debate” in the Philosophy of Biology to the “Venom Debate” in the Science of Toxinology. Toxins (Basel) 2016, 8. [Google Scholar] [CrossRef] [PubMed]
- Losos, J.B.; Hillis, D.M.; Greene, H.W. Evolution. Who speaks with a forked tongue? Science 2012, 338, 1428–1429. [Google Scholar] [CrossRef] [PubMed]
- Schwenk, K.; Wagner, G.P. Function and the evolution of phenotypic stability: Connecting pattern to process. Am. Zool. 2001, 41, 552–563. [Google Scholar] [CrossRef]
- Vitt, L.J. Walking the Natural-History Trail. Herpetologica 2013, 69, 105–117. [Google Scholar] [CrossRef]
- Gauthier, J.A.; Kearney, M.; Maisano, J.A.; Rieppel, O.; Behlke, A.D.B. Assembling the Squamate tree of life: Perspectives from the phenotype and the fossil record. Bull. Peabody Mus. Nat. Hist. 2012, 53, 3–308. [Google Scholar] [CrossRef]
- Estes, R.; de Queiroz, K.; Gauthier, J. Phylogenetic relationships of squamate reptiles. In Phylogenetic Relationships of the Lizard Families; Stanford University Press: Stanford, CA, USA, 1988; pp. 119–280. [Google Scholar]
- Pyron, R.A.; Burbrink, F.T. Extinction, ecological opportunity, and the origins of global snake diversity. Evolution 2012, 66, 163–178. [Google Scholar] [CrossRef] [PubMed]
- Reeder, T.W.; Townsend, T.M.; Mulcahy, D.G.; Noonan, B.P.; Wood, P.L., Jr.; Sites, J.W., Jr.; Wiens, J.J. Integrated analyses resolve conflicts over squamate reptile phylogeny and reveal unexpected placements for fossil taxa. PLoS ONE 2015, 10. [Google Scholar] [CrossRef] [PubMed]
- Townsend, T.; Larson, A.; Louis, E.; Macey, J.R. Molecular phylogenetics of squamata: The position of snakes, amphisbaenians, and dibamids, and the root of the squamate tree. Syst. Biol. 2004, 53, 735–757. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; David, P. New insights into the early history of snakes inferred from two nuclear genes. Mol. Phylogenet. Evol. 2004, 31, 783–787. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Hedges, S.B. Higher-level relationships of caenophidian snakes inferred from four nuclear and mitochondrial genes. C. R. Biol. 2002, 325, 987–995. [Google Scholar] [CrossRef]
- Vidal, N.; Hedges, S.B. Molecular evidence for a terrestrial origin of snakes. Proc. Biol. Sci. R. Soc. 2004, 271 (Suppl. S4), S226–S229. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Hedges, S.B. The phylogeny of squamate reptiles (lizards, snakes, and amphisbaenians) inferred from nine nuclear protein-coding genes. C. R. Biol. 2005, 328, 1000–1008. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Hedges, S.B. The molecular evolutionary tree of lizards, snakes, and amphisbaenians. C. R. Biol. 2009, 332, 129–139. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Rage, J.C.; Couloux, A.; Hedges, S.B. Snakes (Serpentes). In The Timetree of Life; Oxford University Press: New York, NY, USA, 2009. [Google Scholar]
- Wiens, J.J.; Hutter, C.R.; Mulcahy, D.G.; Noonan, B.P.; Townsend, T.M.; Sites, J.W., Jr.; Reeder, T.W. Resolving the phylogeny of lizards and snakes (Squamata) with extensive sampling of genes and species. Biol. Lett. 2012, 8, 1043–1046. [Google Scholar] [CrossRef] [PubMed]
- Sweet, S.S. Chasing Flamingos: Toxicofera and the Misinterpretation of Venom in Varanid Lizards; Institute for Research and Development, Suan Sunandha Rajabhat University: Bangkok, Thailand, 2016. [Google Scholar]
- Hsiang, A.Y.; Field, D.J.; Webster, T.H.; Behlke, A.D.; Davis, M.B.; Racicot, R.A.; Gauthier, J.A. The origin of snakes: Revealing the ecology, behavior, and evolutionary history of early snakes using genomics, phenomics, and the fossil record. BMC Evol. Biol. 2015, 15, 87. [Google Scholar] [CrossRef] [PubMed]
- Mackessy, S.P.; Saviola, A.J. Understanding biological roles of venoms among the Caenophidia: The importance of rear-fanged snakes. Integr. Comp. Biol. 2016. [Google Scholar] [CrossRef] [PubMed]
- Hargreaves, A.D.; Swain, M.T.; Logan, D.W.; Mulley, J.F. Testing the Toxicofera: Comparative transcriptomics casts doubt on the single, early evolution of the reptile venom system. Toxicon 2014, 92, 140–156. [Google Scholar] [CrossRef] [PubMed]
- Sweet, S.; Pianka, E.R. Monitors, Mammals, and Wallace’s Line. Mertensiella 2007, 16, 79–99. [Google Scholar]
- Ast, J.C. Mitochondrial DNA evidence and evolution in Varanoidea (Squamata). Cladistics 2001, 17, 211–226. [Google Scholar] [CrossRef]
- Thompson, G.G.; Clemente, C.J.; Withers, P.C.; Fry, B.G.; Norman, J.A. Is body shape of varanid lizards linked with retreat choice? Aust. J. Zool. 2008, 56, 351–362. [Google Scholar] [CrossRef]
- Vidal, N.; Marin, J.; Sassi, J.; Battistuzzi, F.U.; Donnellan, S.; Fitch, A.J.; Fry, B.G.; Vonk, F.J.; Rodriguez de la Vega, R.C.; Couloux, A.; et al. Molecular evidence for an Asian origin of monitor lizards followed by Tertiary dispersals to Africa and Australasia. Biol. Lett. 2012, 8, 853–855. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Vidal, N.; Norman, J.A.; Vonk, F.J.; Scheib, H.; Ramjan, S.F.; Kuruppu, S.; Fung, K.; Hedges, S.B.; Richardson, M.K.; et al. Early evolution of the venom system in lizards and snakes. Nature 2006, 439, 584–588. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Roelants, K.; Champagne, D.E.; Scheib, H.; Tyndall, J.D.; King, G.F.; Nevalainen, T.J.; Norman, J.A.; Lewis, R.J.; Norton, R.S.; et al. The toxicogenomic multiverse: Convergent recruitment of proteins into animal venoms. Annu. Rev. Genom. Hum. Genet. 2009, 10, 483–511. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Sunagar, K.; Casewell, N.R.; Kochva, E.; Roelants, K.; Scheib, H.; Wüster, W.; Vidal, N.; Young, B.; Burbrink, F.; et al. The origin and evolution of the Toxicofera reptile venom system. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 1–31. [Google Scholar]
- Fry, B.G.; Undheim, E.A.; Ali, S.A.; Jackson, T.N.; Debono, J.; Scheib, H.; Ruder, T.; Morgenstern, D.; Cadwallader, L.; Whitehead, D.; et al. Squeezers and leaf-cutters: Differential diversification and degeneration of the venom system in toxicoferan reptiles. Mol. Cell. Proteom. 2013, 12, 1881–1899. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Wroe, S.; Teeuwisse, W.; van Osch, M.J.; Moreno, K.; Ingle, J.; McHenry, C.; Ferrara, T.; Clausen, P.; Scheib, H.; et al. A central role for venom in predation by Varanus komodoensis (Komodo dragon) and the extinct giant Varanus (Megalania) priscus. Proc. Natl. Acad. Sci. USA 2009, 106, 8969–8974. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N. Colubroid systematics: Evidence for an early appearance of the venom apparatus followed by extensive evolutionary tinkering. J. Toxicol. Toxin Rev. 2002, 21, 21–41. [Google Scholar] [CrossRef]
- Pyron, R.A.; Burbrink, F.T.; Colli, G.R.; de Oca, A.N.; Vitt, L.J.; Kuczynski, C.A.; Wiens, J.J. The phylogeny of advanced snakes (Colubroidea), with discovery of a new subfamily and comparison of support methods for likelihood trees. Mol. Phylogenet. Evol. 2011, 58, 329–342. [Google Scholar] [CrossRef] [PubMed]
- Pyron, R.A.; Burbrink, F.T.; Wiens, J.J. A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes. BMC Evol. Biol. 2013, 13, 93. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Branch, W.R.; Pauwels, O.S.G.; Hedges, S.B.; Broadley, D.G.; Wink, M.; Cruaud, C.; Joger, U.; Nagy, Z.T. Dissecting the major African snake radiation: A molecular phylogeny of the Lamprophiidae Fitzinger (Serpentes, Caenophidia). Zootaxa 2008, 1945, 51–66. [Google Scholar]
- Vidal, N.; Delmas, A.S.; David, P.; Cruaud, C.; Couloux, A.; Hedges, S.B. The phylogeny and classification of caenophidian snakes inferred from seven nuclear protein-coding genes. C. R. Biol. 2007, 330, 182–187. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Dewynter, M.; Gower, D.J. Dissecting the major American snake radiation: A molecular phylogeny of the Dipsadidae Bonaparte (Serpentes, Caenophidia). C. R. Biol. 2010, 333, 48–55. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Kindl, S.G.; Wong, A.; Hedges, S.B. Phylogenetic relationships of xenodontine snakes inferred from 12S and 16S ribosomal RNA sequences. Mol. Phylogenet. Evol. 2000, 14, 389–402. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G. From genome to “venome”: Molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins. Genome Res. 2005, 15, 403–420. [Google Scholar] [CrossRef] [PubMed]
- Vaiyapuri, S.; Sunagar, K.; Gibbins, J.M.; Jackson, T.N.W.; Reeks, T.; Fry, B.G. Kallikrein Enzymes. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 267–280. [Google Scholar]
- Sunagar, K.; Jackson, T.N.W.; Reeks, T.; Fry, B.G. Cysteine-rich secretory proteins. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 239–246. [Google Scholar]
- Casewell, N.R.; Huttley, G.A.; Wuster, W. Dynamic evolution of venom proteins in squamate reptiles. Nat. Commun. 2012, 3, 1066. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G. Venomous Reptiles & Their Toxins. In Venomous Reptiles & Their Toxins; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015. [Google Scholar]
- Fry, B.G.; Richards, R.; Earl, S.; Cousin, X.; Jackson, T.N.W.; Weise, C.; Sunagar, K. Lesser-Known or Putative Reptile Toxins. In Venomous Reptiles & Their Toxins; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 364–407. [Google Scholar]
- Fry, B.G.; Roelants, K.; Winter, K.; Hodgson, W.C.; Griesman, L.; Kwok, H.F.; Scanlon, D.; Karas, J.; Shaw, C.; Wong, L.; et al. Novel venom proteins produced by differential domain-expression strategies in beaded lizards and gila monsters (genus Heloderma). Mol. Biol. Evol. 2010, 27, 395–407. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Scheib, H.; van der Weerd, L.; Young, B.; McNaughtan, J.; Ramjan, S.F.; Vidal, N.; Poelmann, R.E.; Norman, J.A. Evolution of an arsenal: Structural and functional diversification of the venom system in the advanced snakes (Caenophidia). Mol. Cell. Proteom. 2008, 7, 215–246. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Vidal, N.; van der Weerd, L.; Kochva, E.; Renjifo, C. Evolution and diversification of the Toxicofera reptile venom system. J. Proteom. 2009, 72, 127–136. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Winter, K.; Norman, J.A.; Roelants, K.; Nabuurs, R.J.; van Osch, M.J.; Teeuwisse, W.M.; van der Weerd, L.; McNaughtan, J.E.; Kwok, H.F.; et al. Functional and structural diversification of the Anguimorpha lizard venom system. Mol. Cell. Proteom. 2010, 9, 2369–2390. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Wuster, W. Assembling an arsenal: Origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences. Mol. Biol. Evol. 2004, 21, 870–883. [Google Scholar] [CrossRef] [PubMed]
- Reyes-Velasco, J.; Card, D.C.; Andrew, A.L.; Shaney, K.J.; Adams, R.H.; Schield, D.R.; Casewell, N.R.; Mackessy, S.P.; Castoe, T.A. Expression of venom gene homologs in diverse python tissues suggests a new model for the evolution of snake venom. Mol. Biol. Evol. 2015, 32, 173–183. [Google Scholar] [CrossRef] [PubMed]
- Eldredge, N.; Gould, S.J. On punctuated equilibria. Science 1997, 276, 338–341. [Google Scholar] [CrossRef] [PubMed]
- Gould, S.J. Tempo and mode in the macroevolutionary reconstruction of Darwinism. Proc. Natl. Acad. Sci. USA 1994, 91, 6764–6771. [Google Scholar] [CrossRef] [PubMed]
- Gould, S.J.; Eldredge, N. Punctuated equilibrium comes of age. Nature 1993, 366, 223–227. [Google Scholar] [CrossRef] [PubMed]
- Gould, S.J. Punctuated equilibrium: Empirical response. Science 1986, 232, 439. [Google Scholar] [CrossRef] [PubMed]
- Gould, S.J. Punctuated equilibrium and the fossil record. Science 1983, 219, 439–440. [Google Scholar] [CrossRef] [PubMed]
- Sunagar, K.; Moran, Y. The Rise and Fall of an Evolutionary Innovation: Contrasting Strategies of Venom Evolution in Ancient and Young Animals. PLoS Genet. 2015, 11. [Google Scholar] [CrossRef] [PubMed]
- Jackson, T.N.; Sunagar, K.; Undheim, E.A.; Koludarov, I.; Chan, A.H.; Sanders, K.; Ali, S.A.; Hendrikx, I.; Dunstan, N.; Fry, B.G. Venom down under: Dynamic evolution of Australian elapid snake toxins. Toxins (Basel) 2013, 5, 2621–2655. [Google Scholar] [CrossRef] [PubMed]
- Jackson, T.N.W.; Koludarov, I.; Ali, S.A.; Dobson, J.; Zdenek, C.N.; Dashevsky, D.; op den Brouw, B.; Masci, P.; Nouwens, A.; Josh, P.; et al. Rapid radiations and the race to redundancy: An investigation of the evolution of Australian elapid snake venoms. Toxins 2016, in press. [Google Scholar]
- Fry, B.G. Structure-function properties of venom components from Australian elapids. Toxicon 1999, 37, 11–32. [Google Scholar] [CrossRef]
- Fry, B.G.; Casewell, N.R.; Wuster, W.; Vidal, N.; Young, B.; Jackson, T.N. The structural and functional diversification of the Toxicofera reptile venom system. Toxicon 2012, 60, 434–448. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Lumsden, N.G.; Wuster, W.; Wickramaratna, J.C.; Hodgson, W.C.; Kini, R.M. Isolation of a neurotoxin (alpha-colubritoxin) from a nonvenomous colubrid: Evidence for early origin of venom in snakes. J. Mol. Evol. 2003, 57, 446–452. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Roelants, K.; Norman, J.A. Tentacles of venom: Toxic protein convergence in the Kingdom Animalia. J. Mol. Evol. 2009, 68, 311–321. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Wuster, W.; Kini, R.M.; Brusic, V.; Khan, A.; Venkataraman, D.; Rooney, A.P. Molecular evolution and phylogeny of elapid snake venom three-finger toxins. J. Mol. Evol. 2003, 57, 110–129. [Google Scholar] [CrossRef] [PubMed]
- Fry, B.G.; Wuster, W.; Ramjan, S.F.R.; Jackson, T.; Martelli, P.; Kini, R.M. Analysis of Colubroidea snake venoms by liquid chromatography with mass spectrometry: Evolutionary and toxinological implications. Rapid Commun. Mass Spectrom. 2003, 17, 2047–2062. [Google Scholar] [CrossRef] [PubMed]
- Vidal, N.; Lecointre, G. Weighting and congruence: A case study based on three mitochondrial genes in pitvipers. Mol. Phylogenet. Evol. 1998, 9, 366–374. [Google Scholar] [CrossRef] [PubMed]
- Miller, A.K.; Maritz, B.; McKay, S.; Glaudas, X.; Alexander, G.J. An ambushers arsenal: Chemical crypsis in the puff adder (Bitis arietans). Proc. Biol. Sci. 2015, 282, 20152182. [Google Scholar] [CrossRef] [PubMed]
- Jackson, T.N.W.; Young, B.; McCarthy, C.J.; Kochva, E.; Vidal, N.; Underwood, G.; Fry, B.G. Endless forms most beautiful: The evolution of ophidian oral glands, including the venom system, and the use of appropriate terminology for homologous structures. Zoomorphology 2016, in press. [Google Scholar]
- Mitchell, J.S.; Heckert, A.B.; Sues, H.D. Grooves to tubes: Evolution of the venom delivery system in a Late Triassic “reptile”. Naturwissenschaften 2010, 97, 1117–1121. [Google Scholar] [CrossRef] [PubMed]
- Szaniawski, H. The earliest known venomous animals recognized among conodonts. Acta Palaeontol. Pol. 2009, 54, 669–676. [Google Scholar] [CrossRef]
- Reynoso, V.H. Possible evidence of a venom apparatus in a Middle Jurassic sphenodontian from the Huizachal red beds of Tamaulipas, Mexico. J. Vertebr. Paleontol. 2005, 25, 646–654. [Google Scholar] [CrossRef]
- Cuenca-Bescos, G.; Rofes, J. First evidence of poisonous shrews with an envenomation apparatus. Naturwissenschaften 2007, 94, 113–116. [Google Scholar] [CrossRef] [PubMed]
- Dufton, M.J. Venomous mammals. Pharmacol. Ther. 1992, 53, 199–215. [Google Scholar] [CrossRef]
- Ligabue-Braun, R.; Verli, H.; Carlini, C.R. Venomous mammals: A review. Toxicon 2012, 59, 680–695. [Google Scholar] [CrossRef] [PubMed]
- Rofes, J.; Cuenca-Bescos, G. First record of Beremendia fissidens (Mammalia, Soricidae) in the Pleistocene of the Iberian Peninsula, with a review of the biostratigraphy, biogeography and palaeoecology of the species. C. R. Palevol 2009, 8, 21–37. [Google Scholar] [CrossRef]
- Gong, E.; Martin, L.D.; Burnham, D.A.; Falk, A.R. The birdlike raptor Sinornithosaurus was venomous. Proc. Natl. Acad. Sci. USA 2010, 107, 766–768. [Google Scholar] [CrossRef] [PubMed]
- Vonk, F.J.; Admiraal, J.F.; Jackson, K.; Reshef, R.; de Bakker, M.A.G.; Vanderschoot, K.; van den Berge, I.; van Atten, M.; Burgerhout, E.; Beck, A.; et al. Evolutionary origin and development of snake fangs. Nature 2008, 454, 630–633. [Google Scholar] [CrossRef] [PubMed]
- Stuebing, R.B.; Inger, R.F. Field Guide to the Snakes of Borneo; Natural History Publications: Borneo, Malaysia, 1993. [Google Scholar]
- Tweedie, M.W.F. The Snakes of Malaya; Singapore National Printers: Singapore, 1983. [Google Scholar]
- Gopalakrishnakone, P.; Kochva, E. Unusual aspects of the venom apparatus of the blue coral snake Maticora bivirgata. Arch. Histol. Cytol. 1990, 53, 199–210. [Google Scholar] [CrossRef] [PubMed]
- Lim, K.P.L.; Lim, F.L.K. A Guide to the Amphibians and Reptiles of Singapore; Singapore Science Centre: Singapore, 1992. [Google Scholar]
- Durkin, J.P.; Pickwell, G.V.; Trotter, J.T.; Shier, W.T. Phospholipase A2 electrophoretic variants in reptile venoms. Toxicon 1981, 19, 535–546. [Google Scholar] [CrossRef]
- Takasaki, C.; Yoshida, H.; Shimazu, T.; Teruuchi, T.; Toriba, M.; Tamiya, N. Studies on the venom components of the long-glanded coral snake, Maticora bivirgata. Toxicon 1991, 29, 191–200. [Google Scholar] [CrossRef]
- Tan, C.H.; Fung, S.Y.; Yap, M.K.; Leong, P.K.; Liew, J.L.; Tan, N.H. Unveiling the elusive and exotic: Venomics of the Malayan blue coral snake (Calliophis bivirgata flaviceps). J. Proteom. 2016, 132, 1–12. [Google Scholar] [CrossRef] [PubMed]
- Utkin, Y.; Sunagar, K.; Jackson, T.N.W.; Reeks, T.; Fry, B.G. Three-Finger Toxins (3FTxs). In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 215–227. [Google Scholar]
- Eitan, M.; Fowler, E.; Herrmann, R.; Duval, A.; Pelhate, M.; Zlotkin, E. A scorpion venom neurotoxin paralytic to insects that affects sodium current inactivation: Purification, primary structure, and mode of action. Biochemistry 1990, 29, 5941–5947. [Google Scholar] [CrossRef] [PubMed]
- Grieco, T.M.; Raman, I.M. Production of resurgent current in NaV1.6-null Purkinje neurons by slowing sodium channel inactivation with beta-pompilidotoxin. J. Neurosci. 2004, 24, 35–42. [Google Scholar] [CrossRef] [PubMed]
- Gur Barzilai, M.; Kahn, R.; Regev, N.; Gordon, D.; Moran, Y.; Gurevitz, M. The specificity of Av3 sea anemone toxin for arthropods is determined at linker DI/SS2-S6 in the pore module of target sodium channels. Biochem. J. 2014, 463, 271–277. [Google Scholar] [CrossRef] [PubMed]
- Jouiaei, M.; Sunagar, K.; Federman Gross, A.; Scheib, H.; Alewood, P.F.; Moran, Y.; Fry, B.G. Evolution of an ancient venom: Recognition of a novel family of cnidarian toxins and the common evolutionary origin of sodium and potassium neurotoxins in sea anemone. Mol. Biol. Evol. 2015, 32, 1598–1610. [Google Scholar] [CrossRef] [PubMed]
- Liu, P.; Jo, S.; Bean, B.P. Modulation of neuronal sodium channels by the sea anemone peptide BDS-I. J. Neurophysiol. 2012, 107, 3155–3167. [Google Scholar] [CrossRef] [PubMed]
- Moran, Y.; Gordon, D.; Gurevitz, M. Sea anemone toxins affecting voltage-gated sodium channels—Molecular and evolutionary features. Toxicon 2009, 54, 1089–1101. [Google Scholar] [CrossRef] [PubMed]
- Schiavon, E.; Stevens, M.; Zaharenko, A.J.; Konno, K.; Tytgat, J.; Wanke, E. Voltage-gated sodium channel isoform-specific effects of pompilidotoxins. FEBS J. 2010, 277, 918–930. [Google Scholar] [CrossRef] [PubMed]
- Shon, K.J.; Grilley, M.M.; Marsh, M.; Yoshikami, D.; Hall, A.R.; Kurz, B.; Gray, W.R.; Imperial, J.S.; Hillyard, D.R.; Olivera, B.M. Purification, characterization, synthesis, and cloning of the lockjaw peptide from Conus purpurascens venom. Biochemistry 1995, 34, 4913–4918. [Google Scholar] [CrossRef] [PubMed]
- Graudins, A.; Wilson, D.; Alewood, P.F.; Broady, K.W.; Nicholson, G.M. Cross-reactivity of Sydney funnel-web spider antivenom: Neutralization of the in vitro toxicity of other Australian funnel-web (Atrax and Hadronyche) spider venoms. Toxicon 2002, 40, 259–266. [Google Scholar] [CrossRef]
- Rash, L.D.; Birinyi-Strachan, L.C.; Nicholson, G.M.; Hodgson, W.C. Neurotoxic activity of venom from the Australian eastern mouse spider (Missulena bradleyi) involves modulation of sodium channel gating. Br. J. Pharmacol. 2000, 130, 1817–1824. [Google Scholar] [CrossRef] [PubMed]
- Brazil, O.V.; Prado-Franceschi, J.; Laure, C.J. Repetitive muscle responses induced by crotamine. Toxicon 1979, 17, 61–67. [Google Scholar] [CrossRef]
- Chang, C.C.; Tseng, K.H. Effect of crotamine, a toxin of South American rattlesnake venom, on the sodium channel of murine skeletal muscle. Br. J. Pharmacol. 1978, 63, 551–559. [Google Scholar] [CrossRef] [PubMed]
- Matavel, A.C.; Ferreira-Alves, D.L.; Beirao, P.S.; Cruz, J.S. Tension generation and increase in voltage-activated Na+ current by crotamine. Eur. J. Pharmacol. 1998, 348, 167–173. [Google Scholar] [CrossRef]
- Peigneur, S.; Orts, D.J.; Prieto da Silva, A.R.; Oguiura, N.; Boni-Mitake, M.; de Oliveira, E.B.; Zaharenko, A.J.; de Freitas, J.C.; Tytgat, J. Crotamine pharmacology revisited: Novel insights based on the inhibition of KV channels. Mol. Pharmacol. 2012, 82, 90–96. [Google Scholar] [CrossRef] [PubMed]
- Rizzi, C.T.; Carvalho-de-Souza, J.L.; Schiavon, E.; Cassola, A.C.; Wanke, E.; Troncone, L.R. Crotamine inhibits preferentially fast-twitching muscles but is inactive on sodium channels. Toxicon 2007, 50, 553–562. [Google Scholar] [CrossRef] [PubMed]
- Yount, N.Y.; Kupferwasser, D.; Spisni, A.; Dutz, S.M.; Ramjan, Z.H.; Sharma, S.; Waring, A.J.; Yeaman, M.R. Selective reciprocity in antimicrobial activity versus cytotoxicity of hBD-2 and crotamine. Proc. Natl. Acad. Sci. USA 2009, 106, 14972–14977. [Google Scholar] [CrossRef] [PubMed]
- Prieto da Silva, A.R.B.; Fry, B.G.; Sunagar, K.; Scheib, H.; Jackson, T.N.W.; Rádis-Baptista, G.; Zaharenko, A.; de Sá, P.L., Jr.; Pereira, A.; Oguiura, N.; et al. Beta-Defensins. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 228–238. [Google Scholar]
- Bao, Y.; Bu, P.; Jin, L.; Hongxia, W.; Yang, Q.; An, L. Purification, characterization and gene cloning of a novel phospholipase A2 from the venom of Agkistrodon blomhoffii ussurensis. Int. J. Biochem. Cell Biol. 2005, 37, 558–565. [Google Scholar] [CrossRef] [PubMed]
- Haupt, J. The Mesothelae—A Monograph of an Exceptional Group of Spiders (Araneae: Mesothelae). Zoologica 2003, 154, 1–102. [Google Scholar]
- Foelix, R.; Erb, B. Mesothelae have venom glands. J. Arachnol. 2010, 38, 596–598. [Google Scholar] [CrossRef]
- Martini, F. The venom apparatus of the fanged blenny, Meiacanthus atrodorsalis. Am. Zool. 1988, 28, A76. [Google Scholar]
- Fishelson, L. Histology and ultrastructure of the recently found buccal toxic gland in the fish Meiacanthus nigrolineatus (Blenniidae). Copeia 1974, 1974, 386–392. [Google Scholar] [CrossRef]
- Smith, W.L.; Stern, J.H.; Girard, M.G.; Davis, M.P. Evolution of venomous cartilaginous and ray-finned fishes. Integr. Comp. Biol. 2016. [Google Scholar] [CrossRef] [PubMed]
- Polis, G.A. The Biology of Scorpions; Stanford University Press: Palo Alto, CA, USA, 1990. [Google Scholar]
- Undheim, E.A.; Georgieva, D.N.; Thoen, H.H.; Norman, J.A.; Mork, J.; Betzel, C.; Fry, B.G. Venom on ice: First insights into Antarctic octopus venoms. Toxicon 2010, 56, 897–913. [Google Scholar] [CrossRef] [PubMed]
- Allcock, A.L.; Hochberg, F.G.; Rodhouse, P.G.K.; Thorpe, J.P. Adelieledone, a new genus of octopodid from the Southern Ocean. Antarct. Sci. 2003, 15, 415–424. [Google Scholar] [CrossRef]
- Goncalves, J.M. Estudos sobre venenos de serpents brasileiras. II. Crotalus terrificus crotaminicus, subespecie biologica. Ann. Acad. Bras. Cienc. 1956, 28, 365–367. [Google Scholar]
- Gonçalves, J.M. Purification and properties of crotamine. In Venoms; Buckley, E.E., Porges, N., Eds.; American Association for the Advancement of Science: Washington, DC, USA, 1956. [Google Scholar]
- Broad, A.J.; Sutherland, S.K.; Coulter, A.R. Lethality in mice of dangerous Australian and other snake venom. Toxicon 1979, 17, 661–664. [Google Scholar] [CrossRef]
- Trabi, M.; Sunagar, K.; Jackson, T.N.W.; Fry, B.G. Factor Xa Enzymes. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 261–266. [Google Scholar]
- Earl, S.; Sunagar, K.; Jackson, T.N.W.; Reeks, T.; Fry, B.G. Factor Va Proteins. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 255–260. [Google Scholar]
- Herrera, M.; Fernandez, J.; Vargas, M.; Villalta, M.; Segura, A.; Leon, G.; Angulo, Y.; Paiva, O.; Matainaho, T.; Jensen, S.D.; et al. Comparative proteomic analysis of the venom of the taipan snake, Oxyuranus scutellatus, from Papua New Guinea and Australia: Role of neurotoxic and procoagulant effects in venom toxicity. J. Proteom. 2012, 75, 2128–2140. [Google Scholar] [CrossRef] [PubMed]
- Isbister, G.K. Procoagulant Snake Toxins: Laboratory Studies, Diagnosis, and Understanding Snakebite Coagulopathy. Semin. Thromb. Hemost. 2009, 35, 93–103. [Google Scholar] [CrossRef] [PubMed]
- Harvey, A.L.; Anderson, A.J. Dendrotoxins: Snake toxins that block potassium channels and facilitate neurotransmitter release. Pharmacol. Ther. 1985, 31, 33–55. [Google Scholar] [CrossRef]
- Harvey, A.L.; Karlsson, E. Dendrotoxin from the venom of the green mamba, Dendroaspis angusticeps. A neurotoxin that enhances acetylcholine release at neuromuscular junctions. Naunyn-Schmiedeberg’s Arch. Pharmacol. 1980, 312, 1–6. [Google Scholar] [CrossRef]
- Karlsson, E.; Mbugua, P.M.; Rodriguez-Ithurralde, D. Fasciculins, anticholinesterase toxins from the venom of the green mamba Dendroaspis angusticeps. J. Physiol. (Paris) 1984, 79, 232–240. [Google Scholar]
- Osman, O.H.; Ismail, M.; El-Asmar, M.F. Pharmacological studies of snake (Dendroaspis angusticeps) venom. Toxicon 1973, 11, 185–192. [Google Scholar] [CrossRef]
- Eng, W.S.; Fry, B.G.; Sunagar, K.; Takacs, Z.; Jackson, T.N.W.; Guddat, L.W. Kunitz Peptides. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; pp. 281–290. [Google Scholar]
- Fry, B.G.; Jackson, T.N.W.; Takacs, Z.; Reeks, T.; Sunagar, K. C-Type Natriuretic Peptides. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Oxford University Press: New York, NY, USA, 2015; pp. 318–326. [Google Scholar]
- Aman, J.W.; Imperial, J.S.; Ueberheide, B.; Zhang, M.M.; Aguilar, M.; Taylor, D.; Watkins, M.; Yoshikami, D.; Showers-Corneli, P.; Safavi-Hemami, H.; et al. Insights into the origins of fish hunting in venomous cone snails from studies of Conus tessulatus. Proc. Natl. Acad. Sci. USA 2015, 112, 5087–5092. [Google Scholar] [CrossRef] [PubMed]
- Dutertre, S.; Jin, A.H.; Vetter, I.; Hamilton, B.; Sunagar, K.; Lavergne, V.; Dutertre, V.; Fry, B.G.; Antunes, A.; Venter, D.J.; et al. Evolution of separate predation- and defence-evoked venoms in carnivorous cone snails. Nat. Commun. 2014, 5. [Google Scholar] [CrossRef] [PubMed]
- Jin, A.H.; Israel, M.R.; Inserra, M.C.; Smith, J.J.; Lewis, R.J.; Alewood, P.F.; Vetter, I.; Dutertre, S. delta-Conotoxin SuVIA suggests an evolutionary link between ancestral predator defence and the origin of fish-hunting behaviour in carnivorous cone snails. Proc. Biol. Sci. R. Soc. 2015, 282. [Google Scholar] [CrossRef] [PubMed]
- Gould, S.J. Wonderful LIFE: The Burgess Shale and the Nature of History; W.W. Norton & Company: New York, NY, USA, 1990. [Google Scholar]
- Morris, S.C.; Gould, S.J. Showdown on the Burgess Shale. Nat. Hist. Mag. 1998, 107, 48–55. [Google Scholar]
- Dennett, D.C. Darwin’s Dangerous Idea: Evolution and the Meaning of Life; Simon & Schuster: New York, NY, USA, 1996. [Google Scholar]
- Ferrada, E.; Wagner, A. Evolutionary innovations and the organization of protein functions in genotype space. PLoS ONE 2010, 5. [Google Scholar] [CrossRef] [PubMed]
- Sunagar, K.; Jackson, T.; Undheim, E.; Ali, S.; Antunes, A.; Fry, B. Three-Fingered RAVERs: Rapid Accumulation of Variations in Exposed Residues of Snake Venom Toxins. Toxins 2013, 5, 2172–2208. [Google Scholar] [CrossRef] [PubMed]
- Heyborne, W.H.; Mackessy, S.P. Identification and characterization of a taxon-specific three-finger toxin from the venom of the Green Vinesnake (Oxybelis fulgidus; family Colubridae). Biochimie 2013, 95, 1923–1932. [Google Scholar] [CrossRef] [PubMed]
- Hill, R.E.; Mackessy, S.P. Characterization of venom (Duvernoy’s secretion) from twelve species of colubrid snakes and partial sequence of four venom proteins. Toxicon 2000, 38, 1663–1687. [Google Scholar] [CrossRef]
- Mackessy, S.P.; Sixberry, N.A.; Heyborne, W.H.; Fritts, T. Venom of the Brown Treesnake, Boiga irregularis: Ontogenetic shifts and taxa-specific toxicity. Toxicon 2006, 47, 537–548. [Google Scholar] [CrossRef] [PubMed]
- Margres, M.J.; McGivern, J.J.; Seavy, M.; Wray, K.P.; Facente, J.; Rokyta, D.R. Contrasting modes and tempos of venom expression evolution in two snake species. Genetics 2015, 199, 165–176. [Google Scholar] [CrossRef] [PubMed]
- Pawlak, J.; Mackessy, S.P.; Fry, B.G.; Bhatia, M.; Mourier, G.; Fruchart-Gaillard, C.; Servent, D.; Menez, R.; Stura, E.; Menez, A.; et al. Denmotoxin, a three-finger toxin from the colubrid snake Boiga dendrophila (Mangrove Catsnake) with bird-specific activity. J. Biol. Chem. 2006, 281, 29030–29041. [Google Scholar] [CrossRef] [PubMed]
- Pawlak, J.; Mackessy, S.P.; Sixberry, N.M.; Stura, E.A.; Le Du, M.H.; Menez, R.; Foo, C.S.; Menez, A.; Nirthanan, S.; Kini, R.M. Irditoxin, a novel covalently linked heterodimeric three-finger toxin with high taxon-specific neurotoxicity. FASEB J. 2009, 23, 534–545. [Google Scholar] [CrossRef] [PubMed]
- Weldon, C.L.; Mackessy, S.P. Biological and proteomic analysis of venom from the Puerto Rican Racer (Alsophis portoricensis: Dipsadidae). Toxicon 2010, 55, 558–569. [Google Scholar] [CrossRef] [PubMed]
- Almehdar, H.A.; Adel-Sadek, M.A.; Redwan, E.M. Immunoreactivity and two-dimensional gel-electrophoresis characterization of Egyptian cobra venom proteome. Pak. J. Pharm. Sci. 2015, 28, 59–64. [Google Scholar] [PubMed]
- Chang, L.S.; Chung, C.; Liou, J.C.; Chang, C.W.; Yang, C.C. Novel neurotoxins from Taiwan banded krait (Bungarus multicinctus) venom: Purification, characterization and gene organization. Toxicon 2003, 42, 323–330. [Google Scholar] [CrossRef]
- Eletskii, A.V.; Maslennikov, I.V.; Kukhtina, V.V.; Utkin Iu, N.; Tsetlin, V.I.; Arsen’ev, A.S. Structure and conformational heterogeneity of the weak toxin from the cobra Naja kaouthia venom. Bioorg. Khim. 2001, 27, 89–101. [Google Scholar] [PubMed]
- Kulkeaw, K.; Chaicumpa, W.; Sakolvaree, Y.; Tongtawe, P.; Tapchaisri, P. Proteome and immunome of the venom of the Thai cobra, Naja kaouthia. Toxicon 2007, 49, 1026–1041. [Google Scholar] [CrossRef] [PubMed]
- Leong, P.K.; Fung, S.Y.; Tan, C.H.; Sim, S.M.; Tan, N.H. Immunological cross-reactivity and neutralization of the principal toxins of Naja sumatrana and related cobra venoms by a Thai polyvalent antivenom (Neuro Polyvalent Snake Antivenom). Acta Trop. 2015, 149, 86–93. [Google Scholar] [CrossRef] [PubMed]
- Li, J.; Zhang, H.; Liu, J.; Xu, K. Novel genes encoding six kinds of three-finger toxins in Ophiophagus hannah (king cobra) and function characterization of two recombinant long-chain neurotoxins. Biochem. J. 2006, 398, 233–242. [Google Scholar] [CrossRef] [PubMed]
- Lyukmanova, E.N.; Shenkarev, Z.O.; Shulepko, M.A.; Paramonov, A.S.; Chugunov, A.O.; Janickova, H.; Dolejsi, E.; Dolezal, V.; Utkin, Y.N.; Tsetlin, V.I.; et al. Structural Insight into Specificity of Interactions between Nonconventional Three-finger Weak Toxin from Naja kaouthia (WTX) and Muscarinic Acetylcholine Receptors. J. Biol. Chem. 2015, 290, 23616–23630. [Google Scholar] [CrossRef] [PubMed]
- Malih, I.; Ahmad rusmili, M.R.; Tee, T.Y.; Saile, R.; Ghalim, N.; Othman, I. Proteomic analysis of Moroccan cobra Naja haje legionis venom using tandem mass spectrometry. J. Proteom. 2014, 96, 240–252. [Google Scholar]
- Nirthanan, S.; Gopalakrishnakone, P.; Gwee, M.C.; Khoo, H.E.; Kini, R.M. Non-conventional toxins from Elapid venoms. Toxicon 2003, 41, 397–407. [Google Scholar] [CrossRef]
- Ogay, A.Y.; Rzhevsky, D.I.; Murashev, A.N.; Tsetlin, V.I.; Utkin, Y.N. Weak neurotoxin from Naja kaouthia cobra venom affects haemodynamic regulation by acting on acetylcholine receptors. Toxicon 2005, 45, 93–99. [Google Scholar] [CrossRef] [PubMed]
- Rzhevskii, D.I.; Murashev, A.N.; Kukhtina, V.V.; Tsetlin, V.I.; Utkin Iu, N. The weak neurotoxin from Naja kaouthia cobra venom decreases the arterial blood pressure in rats. Bioorg. Khim. 2001, 27, 221–223. [Google Scholar] [PubMed]
- Starkov, V.G.; Poliak Iu, L.; Vul’fius, E.A.; Kriukova, E.V.; Tsetlin, V.I.; Utkin Iu, N. New weak toxins from the cobra venom. Bioorg. Khim. 2009, 35, 15–24. [Google Scholar] [CrossRef] [PubMed]
- Tan, K.Y.; Tan, C.H.; Fung, S.Y.; Tan, N.H. Venomics, lethality and neutralization of Naja kaouthia (monocled cobra) venoms from three different geographical regions of Southeast Asia. J. Proteom. 2015, 120, 105–125. [Google Scholar] [CrossRef] [PubMed]
- Hayes, W.K.; Lavin-Murcio, P.; Kardong, K.V. Delivery of Duvernoy’s secretion into prey by the brown tree snake, Boiga irregularis (Serpentes:Colubridae). Toxicon 1993, 31, 881–887. [Google Scholar] [CrossRef]
- Kardong, K.V. Evolutionary patterns in advanced snakes. Am. Zool. 1980, 20, 269–282. [Google Scholar] [CrossRef]
- Kardong, K.V. The evolution of the venom apparatus in snakes from colubrids to viperids and elapids. Mem. Inst. Butanan 1982, 46, 105–118. [Google Scholar]
- Kardong, K.V. Replies to Fry et al. (Toxicon 2012, 60/4, 434–448). Part B. Properties and biological roles of squamate oral products: The “venomous lifestyle” and preadaptation. Toxicon 2012, 60, 964–966. [Google Scholar] [CrossRef] [PubMed]
- Kardong, K.W. Colubrid snakes and Duvernoy’s “venom” glands. J. Toxicol. Toxin Rev. 2002, 21, 1–19. [Google Scholar] [CrossRef]
- Rochelle, M.J.; Kardong, K.V. Constriction vs envenomation in prey capture by brown tree snakes (Boiga irregularis) (Squamata, Colubridae). Herpetologica 1993, 49, 301–304. [Google Scholar]
- Lumsden, N.G.; Fry, B.G.; Kini, R.M.; Hodgson, W.C. In vitro neuromuscular activity of ‘colubrid’ venoms: Clinical and evolutionary implications. Toxicon 2004, 43, 819–827. [Google Scholar] [CrossRef] [PubMed]
- Lumsden, N.G.; Fry, B.G.; Ventura, S.; Kini, R.M.; Hodgson, W.C. The in vitro and in vivo pharmacological activity of Boiga dendrophila (mangrove catsnake) venom. Auton. Autacoid Pharmacol. 2004, 24, 107–113. [Google Scholar] [CrossRef] [PubMed]
- Lumsden, N.G.; Fry, B.G.; Ventura, S.; Kini, R.M.; Hodgson, W.C. Pharmacological characterisation of a neurotoxin from the venom of Boiga dendrophila (mangrove catsnake). Toxicon 2005, 45, 329–334. [Google Scholar] [CrossRef] [PubMed]
- Sutherland, S.K.; Tibbals, J. Australian Animal Toxins: The Creatures, Their Toxins and Care of the Poisoned Patient, 2nd ed.; Oxford University Press: New York, NY, USA, 2001. [Google Scholar]
- Modahl, C.M.; Mukherjee, A.K.; Mackessy, S.P. An analysis of venom ontogeny and prey-specific toxicity in the Monocled Cobra (Naja kaouthia). Toxicon 2016, 119, 8–20. [Google Scholar] [CrossRef] [PubMed]
- an, N.H.; Arunmozhiarasi, A.; Ponnudurai, G. A comparative study of the biological properties of Dendroaspis (mamba) snake venoms. Comp. Biochem. Physiol. Part C: Comp. Pharmacol. 1991, 99, 463–466. [Google Scholar]
- Vetter, I.; Dekan, Z.; Knapp, O.; Adams, D.J.; Alewood, P.F.; Lewis, R.J. Isolation, characterization and total regioselective synthesis of the novel muO-conotoxin MfVIA from Conus magnificus that targets voltage-gated sodium channels. Biochem. Pharmacol. 2012, 84, 540–548. [Google Scholar] [CrossRef] [PubMed]
Figure 1.
Isolation of calliotoxin (δ-elapitoxin-Cb1a), the first NaV activator from snake venom. (a) Specimen of Calliophis bivirgatus, the blue coral snake (photo by Tom Charlton); (b) Dissected preserved 112 cm Calliophis bivirgatus specimen with 29 cm elongated venom glands (arrows); (c) C. bivirgatus crude venom (10 μg/mL) elicits rising contractions of the skeletal muscle in the chick biventer cervicis nerve-muscle preparation while abolishing nerve mediated contractions; (d) Pre-treatment with tetrodotoxin (0.1 μM) prevents the increase in baseline tension elicited by crude C. bivirgatus venom. * p Value = 0.0272, significantly different from venom alone, unpaired t-test; (e) RP-HPLC fractionation of crude C. bivirgatus venom (150 μg) on a BDS Hypersil C18 column. Dotted line; gradient (25.6%–56.2% solvent B). Solid grey line shows elution of purified calliotoxin (green peak) under identical conditions. Arrow indicates active fraction; (f) Left axis (green line): orthogonal purification of RP-HPLC fraction containing calliotoxin on a Synergi-Hydro RP column. Dotted line; gradient (36.5%–42.1% solvent B). Right axis (black line/filled circles): response of SH-SY5Y cells to corresponding 20 s-fractions. Fractions indicated by blue circles were collected for sequencing and pharmacological analysis; (g) The fraction containing purified native calliotoxin is dominated by a single isotopic mass of 6725.9 Da.
Figure 1.
Isolation of calliotoxin (δ-elapitoxin-Cb1a), the first NaV activator from snake venom. (a) Specimen of Calliophis bivirgatus, the blue coral snake (photo by Tom Charlton); (b) Dissected preserved 112 cm Calliophis bivirgatus specimen with 29 cm elongated venom glands (arrows); (c) C. bivirgatus crude venom (10 μg/mL) elicits rising contractions of the skeletal muscle in the chick biventer cervicis nerve-muscle preparation while abolishing nerve mediated contractions; (d) Pre-treatment with tetrodotoxin (0.1 μM) prevents the increase in baseline tension elicited by crude C. bivirgatus venom. * p Value = 0.0272, significantly different from venom alone, unpaired t-test; (e) RP-HPLC fractionation of crude C. bivirgatus venom (150 μg) on a BDS Hypersil C18 column. Dotted line; gradient (25.6%–56.2% solvent B). Solid grey line shows elution of purified calliotoxin (green peak) under identical conditions. Arrow indicates active fraction; (f) Left axis (green line): orthogonal purification of RP-HPLC fraction containing calliotoxin on a Synergi-Hydro RP column. Dotted line; gradient (36.5%–42.1% solvent B). Right axis (black line/filled circles): response of SH-SY5Y cells to corresponding 20 s-fractions. Fractions indicated by blue circles were collected for sequencing and pharmacological analysis; (g) The fraction containing purified native calliotoxin is dominated by a single isotopic mass of 6725.9 Da.
![Toxins 08 00303 g001](/web/20161031184148im_/http://www.mdpi.com/toxins/toxins-08-00303/article_deploy/html/images/toxins-08-00303-g001.png)
Figure 2.
Sequence of calliotoxin (δ-elapitoxin-Cb1a) and the most closely related three-finger toxins. Homology to known three-finger toxins (3FTxs) was determined using BLAST and multiple sequence alignment conducted with CLUSTAL Omega (1.2.2). Calliotoxin (δ-elapitoxin-Cb1a; blue sequence) is most closely related to ρ-EPTX-Dp1b (P25518.1; Dendroaspis polylepis polylepis; 50% sequence identity), ρ-EPTX-Dp1a (P80495.1; Dendroaspis polylepis polylepis; 50% sequence identity), ρ-EPTX-Da1a (P85092.1; Dendroaspis angusticeps; 49% sequence identity), and ρ-EPTX-Da1b (P86419.1; Dendroaspis angusticeps; 53% sequence identity), three-finger toxins with activity at α1-adrenoreceptors, as well as muscarinic toxin 3 (P81031.2; Dendroaspis angusticeps; 49% sequence identity) targeting the mAChR and Toxin CM-12 (P62394.1; Naja haje haje; 42% sequence identity), Cardiotoxin 7 (P49122.1; Naja atra; 42% sequence identity), and Three-finger toxin 6 (JAS05190.1; Micrurus tener; 47% sequence identity) with unknown pharmacological activity. Conserved ancestral cysteines shown in bold, disulphide-bond pattern shown in black lines.
Figure 2.
Sequence of calliotoxin (δ-elapitoxin-Cb1a) and the most closely related three-finger toxins. Homology to known three-finger toxins (3FTxs) was determined using BLAST and multiple sequence alignment conducted with CLUSTAL Omega (1.2.2). Calliotoxin (δ-elapitoxin-Cb1a; blue sequence) is most closely related to ρ-EPTX-Dp1b (P25518.1; Dendroaspis polylepis polylepis; 50% sequence identity), ρ-EPTX-Dp1a (P80495.1; Dendroaspis polylepis polylepis; 50% sequence identity), ρ-EPTX-Da1a (P85092.1; Dendroaspis angusticeps; 49% sequence identity), and ρ-EPTX-Da1b (P86419.1; Dendroaspis angusticeps; 53% sequence identity), three-finger toxins with activity at α1-adrenoreceptors, as well as muscarinic toxin 3 (P81031.2; Dendroaspis angusticeps; 49% sequence identity) targeting the mAChR and Toxin CM-12 (P62394.1; Naja haje haje; 42% sequence identity), Cardiotoxin 7 (P49122.1; Naja atra; 42% sequence identity), and Three-finger toxin 6 (JAS05190.1; Micrurus tener; 47% sequence identity) with unknown pharmacological activity. Conserved ancestral cysteines shown in bold, disulphide-bond pattern shown in black lines.
![Toxins 08 00303 g002](/web/20161031184148im_/http://www.mdpi.com/toxins/toxins-08-00303/article_deploy/html/images/toxins-08-00303-g002.png)
Figure 3.
Activity of calliotoxin (200 nM) in HEK-293 cells heterologously expressing hNaV1.4 assessed by automated patch clamping. (a) Representative trace of sodium currents before (upper) and after addition of calliotoxin (lower) elicited by depolarizing steps between −80 and +60 mV in 10 mV increments. The red trace highlights the depolarizing step to −20 mV; (b) Current-voltage relationship before and after the addition of calliotoxin; (c) Voltage-dependence of activation (circles) and voltage-dependence of steady-state fast inactivation (squares) before and after addition of calliotoxin; (d) Voltage-dependence of fast inactivation time constants before and after addition of calliotoxin; (e) Average persistent current elicited 40–50 ms after a −20 mV depolarizing step before and after addition of calliotoxin; (f) Representative trace of ramp sodium current before and after addition of calliotoxin; (g) Peak inward current elicited by a depolarizing ramp (2.4 mV/ms) before and after addition of calliotoxin. Data are shown as mean ± SEM with n = 3 replicates; (h) Calliotoxin-induced effects are mediated through NaV channels. Calliotoxin-induced Ca2+ responses in SH-SY5Y cells were inhibited by tetrodotoxin (TTX) (1 μM) but not the nAChR antagonist d-tubocurarine (10 μM), the M1 mAChR antagonist pirenzepine (100 μM), or the α1-adrenoreceptor antagonist prazosin (10 μM).
Figure 3.
Activity of calliotoxin (200 nM) in HEK-293 cells heterologously expressing hNaV1.4 assessed by automated patch clamping. (a) Representative trace of sodium currents before (upper) and after addition of calliotoxin (lower) elicited by depolarizing steps between −80 and +60 mV in 10 mV increments. The red trace highlights the depolarizing step to −20 mV; (b) Current-voltage relationship before and after the addition of calliotoxin; (c) Voltage-dependence of activation (circles) and voltage-dependence of steady-state fast inactivation (squares) before and after addition of calliotoxin; (d) Voltage-dependence of fast inactivation time constants before and after addition of calliotoxin; (e) Average persistent current elicited 40–50 ms after a −20 mV depolarizing step before and after addition of calliotoxin; (f) Representative trace of ramp sodium current before and after addition of calliotoxin; (g) Peak inward current elicited by a depolarizing ramp (2.4 mV/ms) before and after addition of calliotoxin. Data are shown as mean ± SEM with n = 3 replicates; (h) Calliotoxin-induced effects are mediated through NaV channels. Calliotoxin-induced Ca2+ responses in SH-SY5Y cells were inhibited by tetrodotoxin (TTX) (1 μM) but not the nAChR antagonist d-tubocurarine (10 μM), the M1 mAChR antagonist pirenzepine (100 μM), or the α1-adrenoreceptor antagonist prazosin (10 μM).
![Toxins 08 00303 g003](/web/20161031184148im_/http://www.mdpi.com/toxins/toxins-08-00303/article_deploy/html/images/toxins-08-00303-g003.png)
© 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).