Diet and snake venom evolution: can local selection alone explain intraspecific venom variation?

A symphony of destruction: Dynamic differential fibrinogenolytic toxicity by rattlesnake (Crotalus and Sistrurus) venoms

What factors influence the evolution of a heavily selected functional trait in a diverse clade? This study adopts rattlesnakes as a model group to investigate the evolutionary history of venom coagulotoxicity in the wider context of phylogenetics, natural history, and biology. Venom-induced clotting of human plasma and fibrinogen was determined and mapped onto the rattlesnake phylogenetic tree to reconstruct the evolution of coagulotoxicity across the group. Our results indicate that venom phenotype is often independent of phylogenetic relationships in rattlesnakes, suggesting the importance of diet and/or other environmental variables in driving venom evolution. Moreover, the striking inter- and intraspecific variability in venom activity on human blood highlights the considerable variability faced by physicians treating envenomation. This study is the most comprehensive effort to date to describe and characterize the evolutionary and biological aspects of coagulotoxins in rattlesnake venom. Further research at finer taxonomic levels is recommended to elucidate patterns of variation within species and lineages.

Pan-American Lancehead Pit-Vipers: Coagulotoxic Venom Effects and Antivenom Neutralisation of Bothrops asper and B. atrox Geographical Variants

The toxin composition of snake venoms and, thus, their functional activity, can vary between and within species. Intraspecific venom variation across a species’ geographic range is a major concern for antivenom treatment of envenomations, particularly for countries like French Guiana that lack a locally produced antivenom. Bothrops asper and Bothrops atrox are the most medically significant species of snakes in Latin America, both producing a variety of clinical manifestations, including systemic bleeding. These pathophysiological actions are due to the activation by the venom of the blood clotting factors Factor X and prothrombin, thereby causing severe consumptive coagulopathy. Both species are extremely wide-ranging, and previous studies have shown their venoms to exhibit regional venom variation. In this study, we investigate the differential coagulotoxic effects on human plasma of six venoms (four B. asper and two B. atrox samples) from different geographic locations, spanning from Mexico to Peru. We assessed how the venom variation of these venom samples affects neutralisation by five regionally available antivenoms: Antivipmyn, Antivipmyn-Tri, PoliVal-ICP, Bothrofav, and Soro Antibotrópico (SAB). The results revealed both inter- and intraspecific variations in the clotting activity of the venoms. These variations in turn resulted in significant variation in antivenom efficacy against the coagulotoxic effects of these venoms. Due to variations in the venoms used in the antivenom production process, antivenoms differed in their species-specific or geographical neutralisation capacity. Some antivenoms (PoliVal-ICP, Bothrofav, and SAB) showed species-specific patterns of neutralisation, while another antivenom (Antivipmyn) showed geographic-specific patterns of neutralisation. This study adds to current knowledge of Bothrops venoms and also illustrates the importance of considering evolutionary biology when developing antivenoms. Therefore, these results have tangible, real-world implications by aiding evidence-based design of antivenoms for treatment of the envenomed patient. We stress that these in vitro studies must be backed by future in vivo studies and clinical trials before therapeutic guidelines are issued regarding specific antivenom use in a clinical setting.

High levels of functional divergence in toxicity towards prey among the venoms of individual pigmy rattlesnakes

Venom is a complex molecular phenotype that shows high levels of variation in expressed proteins between individuals within and between populations. However, the functional significance of this variation in terms of toxicity towards prey is largely unknown. Here, we assessed the relative toxicity of venom from individual pygmy rattlesnakes ( Sistrurus miliarius) on brown anoles ( Anolis sagrei) using a novel assay involving tests of fixed doses of venom from individual snakes on individual lizards. We found high levels of functional variation between individual venoms within populations with individual differences (nested within population) explaining 3.6 times more variation in toxicity than population differences. Our results suggest a previously unappreciated adaptive significance to within-population variation in venom. They argue that selective mechanisms that maintain variation within populations may be of equal or greater importance to divergent selection leading to local adaption between populations as evolutionary explanations of venom variation within species.

Snake Venoms in Drug Discovery: Valuable Therapeutic Tools for Life Saving

Animal venoms are used as defense mechanisms or to immobilize and digest prey. In fact, venoms are complex mixtures of enzymatic and non-enzymatic components with specific pathophysiological functions. Peptide toxins isolated from animal venoms target mainly ion channels, membrane receptors and components of the hemostatic system with high selectivity and affinity. The present review shows an up-to-date survey on the pharmacology of snake-venom bioactive components and evaluates their therapeutic perspectives against a wide range of pathophysiological conditions. Snake venoms have also been used as medical tools for thousands of years especially in tradition Chinese medicine. Consequently, snake venoms can be considered as mini-drug libraries in which each drug is pharmacologically active. However, less than 0.01% of these toxins have been identified and characterized. For instance, Captopril^® (Enalapril), Integrilin^® (Eptifibatide) and Aggrastat^® (Tirofiban) are drugs based on snake venoms, which have been approved by the FDA. In addition to these approved drugs, many other snake venom components are now involved in preclinical or clinical trials for a variety of therapeutic applications. These examples show that snake venoms can be a valuable source of new principle components in drug discovery.

Characterization of the Venom of C. d. cumanesis of Colombia: Proteomic Analysis and Antivenomic Study

The Colombian rattlesnake Crotalus durissus cumanensis is distributed in three geographic zones of the country: the Atlantic Coast, the upper valley of the Magdalena River, and the eastern plains of the Colombian Orinoquía. Its venom induces neurological symptoms, such as eyelid ptosis, myasthenic facies, and paralysis of the respiratory muscles, which can lead to death. Identification and analysis of C. d. cumanensis showed nine groups of proteins responsible for the neurotoxic effect, of which the crotoxin complex was the most abundant (64.71%). Immunorecognition tests of C. d. cumanensis showed that the use of a commercial antivenom manufactured in Mexico resulted in immunoreactivity.

Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms

Central and South American pitvipers, belonging to the genera Bothrops and Bothriechis, have independently evolved arboreal tendencies. Little is known regarding the composition and activity of their venoms. In order to close this knowledge gap, venom proteomics and toxin activity of species of Bothriechis, and Bothrops (including Bothriopsis) were investigated through established analytical methods. A combination of proteomics and bioactivity techniques was used to demonstrate a similar diversification of venom composition between large and small species within Bothriechis and Bothriopsis. Increasing our understanding of the evolution of complex venom cocktails may facilitate future biodiscoveries.

Age-related association of venom gene expression and diet of predatory gastropods

Background

Venomous organisms serve as wonderful systems to study the evolution and expression of genes that are directly associated with prey capture. To evaluate the relationship between venom gene expression and prey utilization, we examined these features among individuals of different ages of the venomous, worm-eating marine snail Conus ebraeus. We determined expression levels of six genes that encode venom components, used a DNA-based approach to evaluate the identity of prey items, and compared patterns of venom gene expression and dietary specialization.

Results

C. ebraeus exhibits two major shifts in diet with age—an initial transition from a relatively broad dietary breadth to a narrower one and then a return to a broader diet. Venom gene expression patterns also change with growth. All six venom genes are up-regulated in small individuals, down-regulated in medium-sized individuals, and then either up-regulated or continued to be down-regulated in members of the largest size class. Venom gene expression is not significantly different among individuals consuming different types of prey, but instead is coupled and slightly delayed with shifts in prey diversity.

Conclusion

These results imply that changes in gene expression contribute to intraspecific variation of venom composition and that gene expression patterns respond to changes in the diversity of food resources during different growth stages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0592-5) contains supplementary material, which is available to authorized users.

Effects of geographical heterogeneity in species interactions on the evolution of venom genes

Geographical heterogeneity in the composition of biotic interactions can create a mosaic of selection regimes that may drive the differentiation of phenotypes that operate at the interface of these interactions. Nonetheless, little is known about effects of these geographical mosaics on the evolution of genes encoding traits associated with species interactions. Predatory marine snails of the family Conidae use venom, a cocktail of conotoxins, to capture prey. We characterized patterns of geographical variation at five conotoxin genes of a vermivorous species, Conus ebraeus, at Hawaii, Guam and American Samoa, and evaluated how these patterns of variation are associated with geographical heterogeneity in prey utilization. All populations show distinct patterns of prey utilization. Three 'highly polymorphic' conotoxin genes showed significant geographical differences in allelic frequency, and appear to be affected by different modes of selection among populations. Two genes exhibited low levels of diversity and a general lack of differentiation among populations. Levels of diversity of 'highly polymorphic' genes exhibit a positive relationship with dietary breadth. The different patterns of evolution exhibited by conotoxin genes suggest that these genes play different roles in prey capture, and that some genes are more greatly affected by differences in predator-prey interactions than others. Moreover, differences in dietary breadth appear to have a greater influence on the differentiation of venoms than differences in the species of prey.

Snake venoms are integrated systems, but abundant venom proteins evolve more rapidly

Background

While many studies have shown that extracellular proteins evolve rapidly, how selection acts on them remains poorly understood. We used snake venoms to understand the interaction between ecology, expression level, and evolutionary rate in secreted protein systems. Venomous snakes employ well-integrated systems of proteins and organic constituents to immobilize prey. Venoms are generally optimized to subdue preferred prey more effectively than non-prey, and many venom protein families manifest positive selection and rapid gene family diversification. Although previous studies have illuminated how individual venom protein families evolve, how selection acts on venoms as integrated systems, is unknown.

Results

Using next-generation transcriptome sequencing and mass spectrometry, we examined microevolution in two pitvipers, allopatrically separated for at least 1.6 million years, and their hybrids. Transcriptomes of parental species had generally similar compositions in regard to protein families, but for a given protein family, the homologs present and concentrations thereof sometimes differed dramatically. For instance, a phospholipase A₂ transcript comprising 73.4 % of the Protobothrops elegans transcriptome, was barely present in the P. flavoviridis transcriptome (<0.05 %). Hybrids produced most proteins found in both parental venoms. Protein evolutionary rates were positively correlated with transcriptomic and proteomic abundances, and the most abundant proteins showed positive selection. This pattern holds with the addition of four other published crotaline transcriptomes, from two more genera, and also for the recently published king cobra genome, suggesting that rapid evolution of abundant proteins may be generally true for snake venoms. Looking more broadly at Protobothrops, we show that rapid evolution of the most abundant components is due to positive selection, suggesting an interplay between abundance and adaptation.

Conclusions

Given log-scale differences in toxin abundance, which are likely correlated with biosynthetic costs, we hypothesize that as a result of natural selection, snakes optimize return on energetic investment by producing more of venom proteins that increase their fitness. Natural selection then acts on the additive genetic variance of these components, in proportion to their contributions to overall fitness. Adaptive evolution of venoms may occur most rapidly through changes in expression levels that alter fitness contributions, and thus the strength of selection acting on specific secretome components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1832-6) contains supplementary material, which is available to authorized users.

VTBuilder: a tool for the assembly of multi isoform transcriptomes

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

Secreted phospholipases A2 of snake venoms: effects on the peripheral neuromuscular system with comments on the role of phospholipases A2 in disorders of the CNS and their uses in industry

Neuro- and myotoxicological signs and symptoms are significant clinical features of envenoming snakebites in many parts of the world. The toxins primarily responsible for the neuro and myotoxicity fall into one of two categories--those that bind to and block the post-synaptic acetylcholine receptors (AChR) at the neuromuscular junction and neurotoxic phospholipases A2 (PLAs) that bind to and hydrolyse membrane phospholipids of the motor nerve terminal (and, in most cases, the plasma membrane of skeletal muscle) to cause degeneration of the nerve terminal and skeletal muscle. This review provides an introduction to the biochemical properties of secreted sPLA2s in the venoms of many dangerous snakes and a detailed discussion of their role in the initiation of the neurologically important consequences of snakebite. The rationale behind the experimental studies on the pharmacology and toxicology of the venoms and isolated PLAs in the venoms is discussed, with particular reference to the way these studies allow one to understand the biological basis of the clinical syndrome. The review also introduces the involvement of PLAs in inflammatory and degenerative disorders of the central nervous system (CNS) and their commercial use in the food industry. It concludes with an introduction to the problems associated with the use of antivenoms in the treatment of neuro-myotoxic snakebite and the search for alternative treatments.

Venomics of New World pit vipers: genus-wide comparisons of venom proteomes across Agkistrodon

We report a genus-wide comparison of venom proteome variation across New World pit vipers in the genus Agkistrodon. Despite the wide variety of habitats occupied by this genus and that all its taxa feed on diverse species of vertebrates and invertebrate prey, the venom proteomes of copperheads, cottonmouths, and cantils are remarkably similar, both in the type and relative abundance of their different toxin families. The venoms from all the eleven species and subspecies sampled showed relatively similar proteolytic and PLA2 activities. In contrast, quantitative differences were observed in hemorrhagic and myotoxic activities in mice. The highest myotoxic activity was observed with the venoms of A. b. bilineatus, followed by A. p. piscivorus, whereas the venoms of A. c. contortrix and A. p. leucostoma induced the lowest myotoxic activity. The venoms of Agkistrodon bilineatus subspecies showed the highest hemorrhagic activity and A. c. contortrix the lowest. Compositional and toxicological analyses agree with clinical observations of envenomations by Agkistrodon in the USA and Central America. A comparative analysis of Agkistrodon shows that venom divergence tracks phylogeny of this genus to a greater extent than in Sistrurus rattlesnakes, suggesting that the distinct natural histories of Agkistrodon and Sistrurus clades may have played a key role in molding the patterns of evolution of their venom protein genes.

BIOLOGICAL SIGNIFICANCE

A deep understanding of the structural and functional profiles of venoms and of the principles governing the evolution of venomous systems is a goal of venomics. Isolated proteomics analyses have been conducted on venoms from many species of vipers and pit vipers. However, making sense of these large inventories of data requires the integration of this information across multiple species to identify evolutionary and ecological trends. Our genus-wide venomics study provides a comprehensive overview of the toxic arsenal across Agkistrodon and a ground for understanding the natural histories of, and clinical observations of envenomations by, species of this genus.

What is the most dangerous snake?

Letter to Editor of Journal of Venomous Animals and Toxins.

Phylogeny-based comparative analysis of venom proteome variation in a clade of rattlesnakes (Sistrurus sp.)

A long-standing question in evolutionary studies of snake venoms is the extent to which phylogenetic divergence and diet can account for between-species differences in venom composition. Here we apply phylogeny-based comparative methods to address this question. We use data on venom variation generated using proteomic techniques for all members of a small clade of rattlesnakes (Sistrurus sp.) and two outgroups for which phylogenetic and diet information is available. We first complete the characterization of venom variation for all members of this clade with a “venomic” analysis of pooled venoms from two members of this genus, S. milarius streckeri and S. m. milarius. These venoms exhibit the same general classes of proteins as those found in other Sistrurus species but differ in their relative abundances of specific protein families. We then test whether there is significant phylogenetic signal in the relative abundances of major venom proteins across species and if diet (measured as percent mammals and lizards among all prey consumed) covaries with venom composition after phylogenetic divergence is accounted for. We found no evidence for significant phylogenetic signal in venom variation: K values for seven snake venom proteins and two composite venom variables [PC 1 and 2]) were all nonsignificant and lower (mean = 0.11+0.06 sd) than mean K values (>0.35) previously reported for a wide range of morphological, life history, physiological and behavioral traits from other species. Finally, analyses based on Phylogenetic Generalized Least Squares (PGLS) methods reveal that variation in abundance of some venom proteins, most strongly CRISP is significantly related to snake diet. Our results demonstrate that venom variation in these snakes is evolutionarily a highly labile trait even among very closely-related taxa and that natural selection acting through diet variation may play a role in molding the relative abundance of specific venom proteins.

Evolution of CRISPs associated with toxicoferan-reptilian venom and mammalian reproduction

Cysteine-rich secretory proteins (CRISPs) are glycoproteins found exclusively in vertebrates and have broad diversified functions. They are hypothesized to play important roles in mammalian reproduction and in reptilian venom, where they disrupt homeostasis of the prey through several mechanisms, including among others, blockage of cyclic nucleotide-gated and voltage-gated ion channels and inhibition of smooth muscle contraction. We evaluated the molecular evolution of CRISPs in toxicoferan reptiles at both nucleotide and protein levels relative to their nonvenomous mammalian homologs. We show that the evolution of CRISP gene in these reptiles is significantly influenced by positive selection and in snakes (ω = 3.84) more than in lizards (ω = 2.33), whereas mammalian CRISPs were under strong negative selection (CRISP1 = 0.55, CRISP2 = 0.40, and CRISP3 = 0.68). The use of ancestral sequence reconstruction, mapping of mutations on the three-dimensional structure, and detailed evaluation of selection pressures suggests that the toxicoferan CRISPs underwent accelerated evolution aided by strong positive selection and directional mutagenesis, whereas their mammalian homologs are constrained by negative selection. Gene and protein-level selection analyses identified 41 positively selected sites in snakes and 14 sites in lizards. Most of these sites are located on the molecular surface (nearly 76% in snakes and 79% in lizards), whereas the backbone of the protein retains a highly conserved structural scaffold. Nearly 46% of the positively selected sites occur in the cysteine-rich domain of the protein. This directional mutagenesis, where the hotspots of mutations are found on the molecular surface and functional domains of the protein, acts as a diversifying mechanism for the exquisite biological targeting of CRISPs in toxicoferan reptiles. Finally, our analyses suggest that the evolution of toxicoferan-CRISP venoms might have been influenced by the specific predatory mechanism employed by the organism. CRISPs in Elapidae, which mostly employ neurotoxins, have experienced less positive selection pressure (ω = 2.86) compared with the "nonvenomous" colubrids (ω = 4.10) that rely on grip and constriction to capture the prey, and the Viperidae, a lineage that mostly employs haemotoxins (ω = 4.19). Relatively lower omega estimates in Anguimorph lizards (ω = 2.33) than snakes (ω = 3.84) suggests that lizards probably depend more on pace and powerful jaws for predation than venom.

Snake venomics of the Central American rattlesnake Crotalus simus and the South American Crotalus durissus complex points to neurotoxicity as an adaptive paedomorphic trend along Crotalus dispersal in South America

We report a comparative venomic and antivenomic characterization of the venoms of newborn and adult specimens of the Central American rattlesnake, Crotalus simus, and of the subspecies cumanensis, durissus, ruruima, and terrificus of South American Crotalus durissus. Neonate and adult C. simus share about 50% of their venom proteome. The venom proteome of 6-week-old C. simus is predominantly made of the neurotoxic heterodimeric phospholipase A(2) (PLA(2) crotoxin) (55.9%) and serine proteinases (36%), whereas snake venom Zn(2+)-metalloproteinases (SVMPs), exclusively of class PIII, represent only 2% of the total venom proteins. In marked contrast, venom from adult C. simus comprises toxins from 7 protein families. A large proportion (71.7%) of these toxins are SVMPs, two-thirds of which belong to the PIII class. These toxin profiles correlate well with the overall biochemical and pharmacological features of venoms from adult (hemorrhagic) and newborn (neurotoxic) C. simus specimens. The venoms of the South American Crotalus subspecies belong to one of two distinct phenotypes. C. d. cumanensis exhibits high levels of SVMPs and low lethal potency (LD(50)), whereas C. d. subspecies terrificus, ruruima, and durissus have low SVMP activity and high neurotoxicity to mice. Their overall toxin compositions explain the outcome of envenomation by these species. Further, in all C. simus and C. durissus venoms, the concentration of neurotoxins (crotoxin and crotamine) is directly related with lethal activity, whereas lethality and metalloproteinase activity show an inverse relationship. The similar venom toxin profiles of newborn C. simus and adult C. durissus terrificus, ruruima, and durissus subspecies strongly suggests that the South American taxa have retained juvenile venom characteristics in the adult form (paedomorphism) along their North-South stepping-stone dispersal. The driving force behind paedomorphism is often competition or predation pressure. The increased concentration of the neurotoxins crotoxin and crotamine in South American rattlesnake venoms strongly argues that the gain of neurotoxicity and lethal venom activities to mammals may have represented the key axis along which overall venom toxicity has evolved during Crotalus durissus invasion of South America. The paedomorphic trend is supported by a decreasing LNC (lethal neurotoxicity coefficient, defined as the ratio between the average LD(50) of the venom and the crotoxin + crotamine concentration) along the North-South axis, coincident with the evolutionary dispersal pattern of the Neotropical rattlesnakes. The indistinguisable immunoreactivity patterns of Costa Rican and Venezuelan polyvalent antivenoms toward C. simus and C. durissus venoms strongly suggest the possibility of using these antivenoms indistinctly for the management of snakebites by adult C. simus and by certain C. d. cumanensis populations exhibiting a hemorrhagic venom phenotype. The antivenomic results also explain why the antivenoms effectively neutralize the hemorrhagic activity of adult C. simus venoms but does not protect against adult C. durissus sp. and newborn C. simus envenomations. The identification of evolutionary trends among tropical Crotalus, as reported here, may have an impact in defining the mixture of venoms for immunization to produce an effective pan-American anti-Crotalus antivenom.

Comparative venom gland transcriptome surveys of the saw-scaled vipers (Viperidae: Echis) reveal substantial intra-family gene diversity and novel venom transcripts

Background

Venom variation occurs at all taxonomical levels and can impact significantly upon the clinical manifestations and efficacy of antivenom therapy following snakebite. Variation in snake venom composition is thought to be subject to strong natural selection as a result of adaptation towards specific diets. Members of the medically important genus Echis exhibit considerable variation in venom composition, which has been demonstrated to co-evolve with evolutionary shifts in diet. We adopt a venom gland transcriptome approach in order to investigate the diversity of toxins in the genus and elucidate the mechanisms which result in prey-specific adaptations of venom composition.

Results

Venom gland transcriptomes were created for E. pyramidum leakeyi, E. coloratus and E. carinatus sochureki by sequencing ~1000 expressed sequence tags from venom gland cDNA libraries. A standardised methodology allowed a comprehensive intra-genus comparison of the venom gland profiles to be undertaken, including the previously described E. ocellatus transcriptome. Blast annotation revealed the presence of snake venom metalloproteinases, C-type lectins, group II phopholipases A₂, serine proteases, L-amino oxidases and growth factors in all transcriptomes throughout the genus. Transcripts encoding disintegrins, cysteine-rich secretory proteins and hyaluronidases were obtained from at least one, but not all, species. A representative group of novel venom transcripts exhibiting similarity to lysosomal acid lipase were identified from the E. coloratus transcriptome, whilst novel metallopeptidases exhibiting similarity to neprilysin and dipeptidyl peptidase III were identified from E. p. leakeyi and E. coloratus respectively.

Conclusion

The comparison of Echis venom gland transcriptomes revealed substantial intrageneric venom variation in representations and cluster numbers of the most abundant venom toxin families. The expression profiles of established toxin groups exhibit little obvious association with venom-related adaptations to diet described from this genus. We suggest therefore that alterations in isoform diversity or transcript expression levels within the major venom protein families are likely to be responsible for prey specificity, rather than differences in the representation of entire toxin families or the recruitment of novel toxin families, although the recruitment of lysosomal acid lipase as a response to vertebrate feeding cannot be excluded. Evidence of marked intrageneric venom variation within the medically important genus Echis strongly advocates further investigations into the medical significance of venom variation in this genus and its impact upon antivenom therapy.

Geographic variation in venom allelic composition and diets of the widespread predatory marine gastropod Conus ebraeus

Background

Members of the predatory gastropod genus Conus use a venom comprised of a cocktail of peptide neurotoxins, termed conotoxins or conopeptides, to paralyze prey and conotoxin gene family members diversify via strong positive selection. Because Conus venoms are used primarily to subdue prey, the evolution of venoms is likely affected by predator-prey interactions.

Methodology/Principal Findings

To identify the selective forces that drive the differentiation of venoms within species of Conus, we examined the distribution of alleles of a polymorphic O-superfamily conotoxin locus of Conus ebraeus at Okinawa, Guam and Hawaii. Previous analyses of mitochondrial cytochrome oxidase I gene sequences suggest that populations of C. ebraeus, a worm-eating Conus, are not structured genetically in the western and central Pacific. Nonetheless, because the sample size from Guam was relatively low, we obtained additional data from this location and reexamined patterns of genetic variation at the mitochondrial gene at Okinawa, Guam and Hawaii. We also utilized a DNA-based approach to identify prey items of individuals of C. ebraeus from Guam and compared this information to published data on diets at Okinawa and Hawaii. Our results show that conotoxin allelic frequencies differ significantly among all three locations, with strongest differentiation at Hawaii. We also confirm previous inferences that C. ebraeus exhibits no genetic differentiation between Okinawa, Guam and Hawaii at the mitochondrial locus. Finally, DNA-based analyses show that eunicid polychaetes comprise the majority of the prey items of C. ebraeus at Guam; while this results compares well with observed diet of this species at Okinawa, C. ebraeus preys predominantly on nereid polychaetes at Hawaii.

Conclusions/Significance

These results imply that strong selection pressures affect conotoxin allelic frequencies. Based on the dietary information, the selection may derive from geographic variation in dietary specialization and local coevolutionary arms races between Conus and their prey.

Coevolution of diet and prey-specific venom activity supports the role of selection in snake venom evolution

The processes that drive the evolution of snake venom variability, particularly the role of diet, have been a topic of intense recent research interest. Here, we test whether extensive variation in venom composition in the medically important viper genus Echis is associated with shifts in diet. Examination of stomach and hindgut contents revealed extreme variation between the major clades of Echis in the proportion of arthropod prey consumed. The toxicity (median lethal dose, LD(50)) of representative Echis venoms to a natural scorpion prey species was found to be strongly associated with the degree of arthropod feeding. Mapping the results onto a novel Echis phylogeny generated from nuclear and mitochondrial sequence data revealed two independent instances of coevolution of venom toxicity and diet. Unlike venom LD(50), the speed with which venoms incapacitated and killed scorpions was not associated with the degree of arthropod feeding. The prey-specific venom toxicity of arthropod-feeding Echis may thus be adaptive primarily by reducing venom expenditure. Overall, our results provide strong evidence that variation in snake venom composition results from adaptive evolution driven by natural selection for different diets, and underscores the need for a multi-faceted, integrative approach to the study of the causes of venom evolution.

Prey envenomation does not improve digestive performance in western diamondback rattlesnakes (Crotalus atrox)

Although the toxic properties of snake venoms have been recognized throughout history, very little is known about the adaptive significance of these powerful mixtures. This study examined the popular hypothesis that prey envenomation enhances digestion by influencing the energetic costs of digestion and assimilation, gut passage time, and apparent assimilation efficiency (ASSIM) in western diamondback rattlesnakes (Crotalus atrox), a species whose venom is recognized for its comparatively high proteolytic activities. A complete randomized block design allowed repeated measures of specific dynamic action and gut passage time to be measured in eight snakes ingesting four feeding treatments (i.e., artificially envenomated live mice, artificially envenomated prekilled mice, saline injected live mice, and saline injected prekilled mice). A second experiment measured ASSIM in eight snakes ingesting a series of six artificially envenomated or six saline injected mice meals over an 8-week period. Contrary to expectations, the results of both these experiments revealed that envenomation had no significant influence on any of the measured digestive performance variables. Gut passage time averaged 6 days and ASSIM averaged 79.1%. Twenty-one hours following ingestion, postprandial metabolic rates exhibited factorial increases that averaged 3.9-fold greater than resting metabolic rate. Specific dynamic action lasted on average 88 hr and accounted for 26% of the total ingested energy. The results of this study reinforce the need to systematically examine the potential adaptive advantages that venoms confer on the snakes that produce them.

The venom gland transcriptome of the Desert Massasauga rattlesnake (Sistrurus catenatus edwardsii): towards an understanding of venom composition among advanced snakes (Superfamily Colubroidea)

Background

Snake venoms are complex mixtures of pharmacologically active proteins and peptides which belong to a small number of superfamilies. Global cataloguing of the venom transcriptome facilitates the identification of new families of toxins as well as helps in understanding the evolution of venom proteomes.

Results

We have constructed a cDNA library of the venom gland of a threatened rattlesnake (a pitviper), Sistrurus catenatus edwardsii (Desert Massasauga), and sequenced 576 ESTs. Our results demonstrate a high abundance of serine proteinase and metalloproteinase transcripts, indicating that the disruption of hemostasis is a principle mechanism of action of the venom. In addition to the transcripts encoding common venom proteins, we detected two varieties of low abundance unique transcripts in the library; these encode for three-finger toxins and a novel toxin possibly generated from the fusion of two genes. We also observed polyadenylated ribosomal RNAs in the venom gland library, an interesting preliminary obsevation of this unusual phenomenon in a reptilian system.

Conclusion

The three-finger toxins are characteristic of most elapid venoms but are rare in viperid venoms. We detected several ESTs encoding this group of toxins in this study. We also observed the presence of a transcript encoding a fused protein of two well-characterized toxins (Kunitz/BPTI and Waprins), and this is the first report of this kind of fusion in a snake toxin transcriptome. We propose that these new venom proteins may have ancillary functions for envenomation. The presence of a fused toxin indicates that in addition to gene duplication and accelerated evolution, exon shuffling or transcriptional splicing may also contribute to generating the diversity of toxins and toxin isoforms observed among snake venoms. The detection of low abundance toxins, as observed in this and other studies, indicates a greater compositional similarity of venoms (though potency will differ) among advanced snakes than has been previously recognized.

Reappraisal of Vipera aspis venom neurotoxicity

Background

The variation of venom composition with geography is an important aspect of intraspecific variability in the Vipera genus, although causes of this variability remain unclear. The diversity of snake venom is important both for our understanding of venomous snake evolution and for the preparation of relevant antivenoms to treat envenomations. A geographic intraspecific variation in snake venom composition was recently reported for Vipera aspis aspis venom in France. Since 1992, cases of human envenomation after Vipera aspis aspis bites in south-east France involving unexpected neurological signs were regularly reported. The presence of genes encoding PLA₂ neurotoxins in the Vaa snake genome led us to investigate any neurological symptom associated with snake bites in other regions of France and in neighboring countries. In parallel, we used several approaches to characterize the venom PLA₂ composition of the snakes captured in the same areas.

Methodology/Principal Findings

We conducted an epidemiological survey of snake bites in various regions of France. In parallel, we carried out the analysis of the genes and the transcripts encoding venom PLA₂s. We used SELDI technology to study the diversity of PLA₂ in various venom samples. Neurological signs (mainly cranial nerve disturbances) were reported after snake bites in three regions of France: Languedoc-Roussillon, Midi-Pyrénées and Provence-Alpes-Côte d'Azur. Genomes of Vipera aspis snakes from south-east France were shown to contain ammodytoxin isoforms never described in the genome of Vipera aspis from other French regions. Surprisingly, transcripts encoding venom neurotoxic PLA₂s were found in snakes of Massif Central region. Accordingly, SELDI analysis of PLA₂ venom composition confirmed the existence of population of neurotoxic Vipera aspis snakes in the west part of the Massif Central mountains.

Conclusions/Significance

The association of epidemiological studies to genetic, biochemical and immunochemical analyses of snake venoms allowed a good evaluation of the potential neurotoxicity of snake bites. A correlation was found between the expression of neurological symptoms in humans and the intensity of the cross-reaction of venoms with anti-ammodytoxin antibodies, which is correlated with the level of neurotoxin (vaspin and/or ammodytoxin) expression in the venom. The origin of the two recently identified neurotoxic snake populations is discussed according to venom PLA₂ genome and transcriptome data.

Mass spectrometric analysis of the individual variability of Bothrops jararaca venom peptide fraction. Evidence for sex-based variation among the bradykinin-potentiating peptides

Variation in the snake venom proteome is well documented and it is a ubiquitous phenomenon at all taxonomical levels. However, variation in the snake venom peptidome is so far not described. In this work we used mass spectrometry [liquid chromatography/mass spectrometry (LC/MS) and matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOFMS)] to explore sex-based differences among the venom peptides of eighteen sibling specimens of Bothrops jararaca of a single litter born and raised in the laboratory. MALDI-TOFMS analyses showed individual variability among the bradykinin-potentiating peptides (BPPs), and, interestingly, four new peptides were detected only in female venoms and identified by de novo sequencing as cleaved BPPs lacking the C-terminal Q-I-P-P sequence. Similar results were obtained with venom from wild-caught adult non-sibling specimens of B. jararaca and in this case we were able to identify the gender of the specimen by analyzing the MALDI-TOF profile of the peptide fraction and finding the cleaved peptides only in female venoms. Synthetic replicates of the cleaved BPPs were less potent than the full-length BPP-10c in potentiating the bradykinin hypotensive effect, suggesting that the C-terminus is critical for the interaction of the BPPs with their mammalian molecular targets. This work represents a comprehensive mass spectrometric analysis of the peptide fraction of B. jararaca venom and shows for the first time sex-based differences in the snake venom peptidome of sibling and non-sibling snakes and suggests that the BPPs may follow distinct processing pathways in female and male individuals.

California Ground Squirrel (Spermophilus beecheyi) Defenses Against Rattlesnake Venom Digestive and Hemostatic Toxins

Previous studies have shown that some mammals are able to neutralize venom from snake predators. California ground squirrels (Spermophilus beecheyi) show variation among populations in their ability to bind venom and minimize damage from northern Pacific rattlesnakes (Crotalus oreganus), but the venom toxins targeted by resistance have not been investigated. Four California ground squirrel populations, selected for differences in local density or type of rattlesnake predators, were assayed for their ability to neutralize digestive and hemostatic effects of venom from three rattlesnake species. In Douglas ground squirrels (S. b. douglasii), we found that animals from a location where snakes are common showed greater inhibition of venom metalloprotease and hemolytic activity than animals from a location where snakes are rare. Effects on general proteolysis were not different. Douglas ground squirrels also reduced the metalloprotease activity of venom from sympatric northern Pacific rattlesnakes (C. o. oreganus) more than the activity of venom from allopatric western diamondback rattlesnakes (C. atrox), but enhanced the fibrinolysis of sympatric venom almost 1.8 times above baseline levels. Two Beechey ground squirrel (S. b. beecheyi) populations had similar inhibition of venoms from northern and southern Pacific rattlesnakes (C. o. helleri), despite differences between the populations in the locally prevalent predator. However, the venom toxins inhibited by Beechey squirrels varied among venom from Pacific rattlesnake subspecies, and between these venoms and venom from allopatric western diamondback rattlesnakes. Blood plasma from Beechey squirrels showed highest inhibition of metalloprotease activity of northern Pacific rattlesnake venom, general proteolytic activity and hemolysis of southern Pacific rattlesnake venom, and hemolysis by allopatric western diamondback venom. These results reveal previously cryptic variation in venom activity against resistant prey that suggests reciprocal adaptation at the molecular level.

California ground squirrel (Spermophilus beecheyi) defenses against rattlesnake venom digestive and hemostatic toxins

Previous studies have shown that some mammals are able to neutralize venom from snake predators. California ground squirrels (Spermophilus beecheyi) show variation among populations in their ability to bind venom and minimize damage from northern Pacific rattlesnakes (Crotalus oreganus), but the venom toxins targeted by resistance have not been investigated. Four California ground squirrel populations, selected for differences in local density or type of rattlesnake predators, were assayed for their ability to neutralize digestive and hemostatic effects of venom from three rattlesnake species. In Douglas ground squirrels (S. b. douglasii), we found that animals from a location where snakes are common showed greater inhibition of venom metalloprotease and hemolytic activity than animals from a location where snakes are rare. Effects on general proteolysis were not different. Douglas ground squirrels also reduced the metalloprotease activity of venom from sympatric northern Pacific rattlesnakes (Crotalus oreganus oreganus) more than the activity of venom from allopatric western diamondback rattlesnakes (C. atrox), but enhanced fibrinolysis of sympatric venom almost 1.8 times above baseline levels. Two Beechey ground squirrel (S. b. beecheyi) populations had similar inhibition of venoms from northern and southern Pacific rattlesnakes (C. o. helleri), despite differences between the populations in the locally prevalent predator. However, the venom toxins inhibited by Beechey squirrels did vary among venom from Pacific rattlesnake subspecies, and between these venoms and venom from allopatric western diamondback rattlesnakes. Blood plasma from Beechey squirrels showed highest inhibition of metalloprotease activity of northern Pacific rattlesnake venom, general proteolytic activity and hemolysis of southern Pacific rattlesnake venom, and hemolysis by allopatric western diamondback venom. These results reveal previously cryptic variation in venom activity against resistant prey that suggests reciprocal adaptation at the molecular level.

A multifaceted analysis of viperid snake venoms by two-dimensional gel electrophoresis: an approach to understanding venom proteomics

The complexity of Viperid venoms has long been appreciated by investigators in the fields of toxinology and medicine. However, it is only recently that the depth of that complexity has become somewhat quantitatively and qualitatively appreciated. With the resurgence of two-dimensional gel electrophoresis (2-DE) and the advances in mass spectrometry virtually all venom components can be visualized and identified given sufficient effort and resources. Here we present the use of 2-DE for examining venom complexity as well as demonstrating interesting approaches to selectively delineate subpopulations of venom proteins based on particular characteristics of the proteins such as antibody cross-reactivity or enzymatic activities. 2-DE comparisons between venoms from different species of the same genus (Bothrops) of snake clearly demonstrated both the similarity as well as the apparent diversity among these venoms. Using liquid chromatography/tandem mass spectrometry we were able to identify regions of the two-dimensional gels from each venom in which certain classes of proteins were found. 2-DE was also used to compare venoms from Crotalus atrox and Bothrops jararaca. For these venoms a variety of staining/detection protocols was utilized to compare and contrast the venoms. Specifically, we used various stains to visualize subpopulations of the venom proteomes of these snakes, including Coomassie, Silver, Sypro Ruby and Pro-Q-Emerald. Using specific antibodies in Western blot analyses of 2-DE of the venoms we have examined subpopulations of proteins in these venoms including the serine proteinase proteome, the metalloproteinase proteome, and the phospholipases A2 proteome. A functional assessment of the gelatinolytic activity of these venoms was also performed by zymography. These approaches have given rise to a more thorough understanding of venom complexity and the toxins comprising these venoms and provide insights to investigators who wish to focus on these venom subpopulations of proteins in future studies.

Snake venom hyaluronidase: an evidence for isoforms and extracellular matrix degradation

The present study attempts to establish the isoforms of hyaluronidase enzyme and their possible role in the spreading of toxins during envenomation. Screening of venoms of 15 snakes belonging to three different families revealed varied hyaluronidase activity in ELISA-like assay, but with relatively similar pH and temperature optima. The zymograms of individual venoms showed varied activity banding patterns and indicated the presence of at least two molecular forms of the enzyme. During envenomation, activity of hyaluronidase is considered crucial for the spreading of toxins and is presumed to distort the integrity of extracellular matrix through the degradation of hyaluronic acid in it. This property has been addressed through localization of hyaluronic acid in human skin and muscle tissue sections using the probe, biotinylated hyaluronic acid binding protein. Faint and discontinuous staining pattern of hyaluronidase treated tissue sections over intense staining of untreated tissue sections confirm the selective degradation of hyaluronic acid in extracellular matrix and thus provide an evidence for the spreading property of the enzyme.

Unusual neurotoxic envenomations by Vipera aspis aspis snakes in France

Vipera aspis aspis (V.a.a.) is the most dangerous poisonous snake in South-Eastern France. The clinical symptoms observed after V.a.a. envenomations involve mostly local signs (pain, edema) associated in the more severe cases with systemic symptoms (gastro-intestinal and cardiovascular manifestations). Since 1992, several unusual cases of moderate and severe 'neurotoxic' envenomations by V.a.a. snakes have been reported in a very localized area in South-Eastern France. Most of the human patients mainly suffered neurological signs owing to cephalic muscle paralysis. Drowsiness and dyspnea were observed for the most severe cases. Envenomed animals suffered respiratory distress and paralysis. The local signs were never as severe as observed after envenomations by vipers in other French regions. Human patients with moderate or severe clinical features received two intravenous injections of Viperfav antivenom, the first dose inducing the decrease of the neurological signs and the second reducing significantly the edema. Neurotoxic components immunologically cross-reacting with toxins from V. ammodytes ammodytes venom from Eastern Europe were detected in the blood of all patients suffering neurological symptoms after a V.a.a. bite. The protective efficacy of various antivenoms was evaluated in mice. The existence of geographical variations in the composition of V.a.a. venom emphasizes on the use of polyvalent antivenom in the treatment of viper envenomations in France.

Variations in biochemical and pharmacological properties of Indian cobra (Naja naja naja) venom due to geographical distribution

Indian cobra (Naja naja naja) venom obtained from three different geographical regions was studied in terms of electrophoretic pattern, biochemical and pharmacological activities. SDS-PAGE banding pattern revealed significant variation in the protein constituents of the three regional venoms. The eastern venom showed highest indirect hemolysis and hyaluronidase activity. In contrast, western and southern venoms were rich in proteolytic activity. All the three regional venoms were devoid of p-tosyl-L-arginine methyl ester hydrolysing activity. The eastern venom was found to be most lethal among the three regional venoms. The lethal potency varied as eastern > western > southern regional venoms. In addition, all the three regional venoms showed marked variations in their ability to induce symptoms/signs of neurotoxicity, myotoxicity, edema and effect on plasma coagulation process. Polyclonal antiserum prepared against the venom of eastern region cross-reacted with both southern and western regional venoms, but varied in the extent of cross-reactivity by ouchterlony immunodiffusion and ELISA.

Phospholipase A(2) with platelet aggregation inhibitor activity from Austrelaps superbus venom: protein purification and cDNA cloning

Four phospholipase A(2) (PLA(2)) enzymes (Superbins a, b, c, and d) with varying platelet aggregation inhibitor activities have been purified from Austrelaps superbus by a combination of gel filtration, ion-exchange, and reversed-phase high-pressure liquid chromatography. Purity and homogeneity of the superbins have been confirmed by high-performance capillary zone electrophoresis and mass spectrometry. The electron spray ionization mass spectrometry data showed that their molecular masses range from 13,140 to 13,236 Da. Each of the proteins has been found to be basic and exhibit varying degrees of PLA(2) activity. They also displayed different platelet aggregation inhibitory activities. Superbin a was found to possess the most potent inhibitory activity with an IC(50) of 9.0 nM, whereas Superbin d was found to be least effective with an IC(50) of 3.0 microM. Superbins b and c were moderately effective with IC(50) values of 0.05 and 0.5 microM, respectively. The amino-terminal sequencing confirmed the identity of these superbins. cDNA cloning resulted in the identification of 17 more PLA(2) isoforms in A. superbus venom. It has also provided complete information on the precursor PLA(2). The precursor PLA(2) contained a 27-amino-acid signal peptide and 117- to 125-amino-acid PLA(2) (molecular mass ranging from 13,000 to 14,000 Da). Two of these PLA(2) enzymes resembled more closely (87%) Superbin a in structure. Two unique PLA(2) enzymes containing an extra pancreatic loop also have been identified among the isoforms.