Sequences and structural organization of phospholipase A2 genes from Vipera aspis aspis, V. aspis zinnikeri and Vipera berus berus venom. Identification of the origin of a new viper population based on ammodytin I1 heterogeneity

Envenomation by Vipera aspis in Piedmont (Italy): A report of three cases, including one case with neurological symptoms

This report describes three cases of human envenomation by the asp viper (Vipera aspis) in Piedmont, north-west Italy. A woman was bitten on the ankle while she was hiking and two herpetologists received bites on the hand while they were manipulating the animals. In the first case, the victim presented severe systemic symptoms (abdominal pain, vomiting, diarrhea) that required treatment with two vials of antivenom and hospitalization for one week. In the second case, the patient manifested neurological symptoms (blepharoptosis, ophtalmoplegia); he was treated with antivenom and discharged after five days. In the third case, the patient was bitten by a juvenile viper and showed only local symptoms (edema and bruising). All patients reported prolonged functional impairment after discharge from hospital. Although uncommon, envenomation by Vipera aspis can cause severe consequences that require immediate management and antivenom administration. These cases highlight the importance of obtaining better knowledge of the intraspecific variability of venoms and its clinical significance, as well as of the factors that determine the severity of snakebite injuries.

Biochemistry and toxicology of proteins and peptides purified from the venom of Vipera berus berus

The isolation and characterization of individual snake venom components is important for a deeper understanding of the pathophysiology of envenomation and for improving the therapeutic procedures of patients. It also opens possibilities for the discovery of novel toxins that might be useful as tools for understanding cellular and molecular processes. The variable venom composition, toxicological and immunological properties of the common vipers (Vipera berus berus) have been reviewed. The combination of venom gland transcriptomics, bottom-up and top-down proteomics enabled comparison of common viper venom proteomes from multiple individuals. V. b. berus venom contains proteins and peptides belonging to 10–15 toxin families: snake venom metalloproteinase, phospholipases A₂ (PLA₂), snake venom serine proteinase, aspartic protease, L-amino acid oxidase (LAAO), hyaluronidase, 5′-nucleotidase, glutaminyl-peptide cyclotransferase, disintegrin, C-type lectin (snaclec), nerve growth factor, Kunitz type serine protease inhibitor, snake venom vascular endothelial growth factor, cysteine-rich secretory protein, bradykinin potentiating peptide, natriuretic peptides. PLA₂ and LAAO from V. b. berus venom produce more pronounced cytotoxic effects in cancer cells than normal cells, via induction of apoptosis, cell cycle arrest and suppression of proliferation. Proteomic data of V. b. berus venoms from different parts of Russia and Slovakian Republic have been compared with analogous data for Vipera nikolskii venom. Proteomic studies demonstrated quantitative differences in the composition of V. b. berus venom from different geographical regions. Differences in the venom composition of V. berus were mainly driven by the age, sex, habitat and diet of the snakes. The venom variability of V. berus results in a loss of antivenom efficacy against snakebites. The effectiveness of antibodies is discussed. This review presents an overview with a special focus on different toxins that have been isolated and characterized from the venoms of V. b. berus. Their main biochemical properties and toxic actions are described.

•

Vipera berus berus venom composition is variable among different populations.

•

Venom contains about 15 protein/peptide families.

•

It disturbs blood coagulation inducing pro- or anticoagulant effects.

•

Venom contains different types of blood factor X activators.

•

PLA₂ and L-amino acid oxidase produce cytotoxic effects in cancer cells.

Individual variability of venom from the European adder (Vipera berus berus) from one locality in Eastern Hungary

We have revealed intra-population variability among venom samples from several individual European adders (Vipera berus berus) within a defined population in Eastern Hungary. Individual differences in venom pattern were noticed, both gender-specific and age-related, by one-dimensional electrophoresis. Gelatin zymography demonstrated that these individual venoms have different degradation profiles indicating varying protease activity in the specimens from adders of different ages and genders. Some specimens shared a conserved region of substrate degradation, while others had lower or extremely low protease activity. Phospholipase A₂ activity of venoms was similar but not identical. Interspecimen diversity of the venom phospholipase A₂-spectra (based on the components' molecular masses) was detected by MALDI-TOF MS. The lethal toxicity of venoms (LD₅₀) also showed differences among individual snakes. Extracted venom samples had varying neuromuscular paralysing effect on chick biventer cervicis nerve-muscle preparations. The paralysing effect of venom was lost when calcium in the physiological salt solution was replaced by strontium; indicating that the block of twitch responses to nerve stimulation is associated with the activity of a phospholipase-dependent neurotoxin. In contrast to the studied V. b. berus venoms from different geographical regions so far, this is the first V. b. berus population discovered to have predominantly neurotoxic neuromuscular activity. The relevance of varying venom yields is also discussed. This study demonstrates that individual venom variation among V. b. berus living in particular area of Eastern Hungary might contribute to a wider range of clinical manifestations of V. b. berus envenoming than elsewhere in Europe.

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

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

BIOLOGICAL SIGNIFICANCE

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

Is Hybridization a Source of Adaptive Venom Variation in Rattlesnakes? A Test, Using a Crotalus scutulatus × viridis Hybrid Zone in Southwestern New Mexico

Venomous snakes often display extensive variation in venom composition both between and within species. However, the mechanisms underlying the distribution of different toxins and venom types among populations and taxa remain insufficiently known. Rattlesnakes (Crotalus, Sistrurus) display extreme inter- and intraspecific variation in venom composition, centered particularly on the presence or absence of presynaptically neurotoxic phospholipases A₂ such as Mojave toxin (MTX). Interspecific hybridization has been invoked as a mechanism to explain the distribution of these toxins across rattlesnakes, with the implicit assumption that they are adaptively advantageous. Here, we test the potential of adaptive hybridization as a mechanism for venom evolution by assessing the distribution of genes encoding the acidic and basic subunits of Mojave toxin across a hybrid zone between MTX-positive Crotalus scutulatus and MTX-negative C. viridis in southwestern New Mexico, USA. Analyses of morphology, mitochondrial and single copy-nuclear genes document extensive admixture within a narrow hybrid zone. The genes encoding the two MTX subunits are strictly linked, and found in most hybrids and backcrossed individuals, but not in C. viridis away from the hybrid zone. Presence of the genes is invariably associated with presence of the corresponding toxin in the venom. We conclude that introgression of highly lethal neurotoxins through hybridization is not necessarily favored by natural selection in rattlesnakes, and that even extensive hybridization may not lead to introgression of these genes into another species.

Hybrid origin of European Vipers (Vipera magnifica and Vipera orlovi) from the Caucasus determined using genomic scale DNA markers

Background

Studying patterns of introgression can illuminate the role of hybridization in speciation, and help guide decisions relevant to the conservation of rare taxa. Vipera magnifica and Vipera orlovi are small vipers that have high conservation status due to their rarity and restricted distributions in an area of the Caucasus region where two other related species are present – V. kaznakovi and V. renardi. Despite numerous observations of hybridization between different species of small vipers, and the potential of a hybrid origin for V. magnifica and V. orlovi based on their distribution with respect to V. kaznakovi and V. renardi, hypotheses of a hybrid origin have not been formally tested. Here we generate genomic-scale data by performing next generation sequencing of double digest restriction-site associated DNA libraries, and use these multilocus data to test whether these two species are of hybrid origin.

Results

We generated over nine hundred loci for 38 specimens of six taxa, and analysed the dataset using Bayesian clustering and multivariate methods, as well as Patterson D-statistics, which can distinguish between incomplete lineage sorting and introgression as explanations for shared polymorphism. The results demonstrate a pattern of historical admixture in the two purported hybrids that is consistent with past gene flow from V. renardi into V. kaznakovi. The average admixture proportion in individuals was low (6.39 %) in the case of V. magnifica, but was higher in V. orlovi (19.02 %). We also show that the specific individual samples used in D-statistic tests can have a significant impact on inferences regarding the magnitude of introgression, suggesting the importance of including multiple individuals in these analyses.

Conclusions

Our results support the conclusion that both V. orlovi and V. magnifica had formed through a hybridization event between V. kaznakovi and V. renardi. Given a low proportion of admixture and absence of clear ecological and morphological differences V. magnifica should be treated as a marginal population of V. kaznakovi. Further studies that include analyses of ecological segregation of V. orlovi from parental taxa and search for evolutionary consequences of hybridisation would clarify if V. orlovi is a distinct hybrid species. Until this we recommend preserving the current taxonomy and protection status of V. orlovi.

Electronic supplementary material

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

Management of snakebites in France

Although not a major health problem in Europe, snakebite in the old continent was the focus of recent studies to evaluate their overall incidence and define management techniques. The purpose of this three-part report is to present the experience of the Marseille Poison Centre with snakebite in France. The first section deals with viper envenomation that now benefits from a validated therapeutic protocol using of purified antivenom of proven efficacy and tolerance in patients showing grade 2 and 3 symptoms. The second section describes the highly variable snakebite situation in French overseas territories that include areas where local species require specialized management, e.g. Martinique and French Guiana. The third section involves the emerging problems associated with the keeping of exotic snakes as pets with problems related to the use of antivenoms from foreign countries. The exotic-snake pets fashion was at the origin of the creation of a national antivenom bank by two French poison centers (Angers and Marseille) to ensure prompt delivery of antivenoms for exotic snake envenomation anywhere in mainland France.

cDNA cloning, structural, and functional analyses of venom phospholipases A₂ and a Kunitz-type protease inhibitor from steppe viper Vipera ursinii renardi

Snake venom phospholipases A₂ (PLA₂s) display a wide array of biological activities and are each characteristic to the venom. Here, we report on the cDNA cloning and characterization of PLA₂s from the steppe viper Vipera ursinii renardi venom glands. Among the five distinct PLA₂ cDNAs cloned and sequenced, the most common were the clones encoding a basic Ser-49 containing PLA₂ (Vur-S49). Other clones encoded either ammodytin analogs I1, I2d and I2a (designated as Vur-PL1, Vur-PL2 and Vur-PL3, respectively) or an ammodytoxin-like PLA₂ (Vurtoxin). Additionally, a novel Kunitz-type trypsin inhibitor for this venom species was cloned and sequenced. Comparison of these PLA₂ and Kunitz inhibitor sequences with those in the sequence data banks suggests that the viper V. u. renardi is closely related to Vipera ammodytes and Vipera aspis. Separation of V. u. renardi venom components by gel-filtration and ion-exchange chromatography showed the presence of many PLA₂ isoforms. Remarkably, the most abundant PLA₂ isolated was Vur-PL2 while Vur-S49 analog was in very low yield. There are great differences between the proportion of cDNA clones and that of the proteins isolated. Two Vur-PL2 isoforms (designated as Vur-PL2A and Vur-PL2B) indistinguishable by masses, peptide mass fingerprinting, N-terminal sequences and CD spectroscopy were purified from the pooled venom. However, when rechromatographed on cation-exchanger, Vur-PL2A showed only one peak corresponding to Vur-PL2B, suggesting the existence of conformers for Vur-PL2. Vur-PL2B was weakly cytotoxic to rat pheochromocytoma PC12 cells and showed both strong anticoagulant and anti-platelet activities. This is the first case of a strong anticoagulating ammodytin I analog in Vipera venom.

Asp Viper (Vipera aspis) envenomation: experience of the Marseille Poison Centre from 1996 to 2008

A retrospective case review study of viper envenomations collected by the Marseille's Poison Centre between 1996 and 2008 was performed.

RESULTS

174 cases were studied (52 grade 1 = G1, 90 G2 and 32 G3). G1 patients received symptomatic treatments (average hospital stay 0.96 day). One hundred and six (106) of the G2/G3 patients were treated with the antivenom Viperfav* (2.1+/-0.9 days in hospital), while 15 of them received symptomatic treatments only (plus one immediate death) (8.1+/-4 days in hospital, 2 of them died). The hospital stay was significantly reduced in the antivenom treated group (p < 0.001), and none of the 106 antivenom treated patients had immediate (anaphylaxis) or delayed (serum sickness) allergic reactions.

CONCLUSION

Viperfav* antivenom was safe and effective for treating asp viper venom-induced toxicity.

Phospholipase A2 diversity and polymorphism in European viper venoms: Paradoxical molecular evolution in Viperinae

We report the diversity and polymorphism of phospholipase A(2) (PLA(2)) transcripts from snakes belonging to nine European viper subspecies. This diversity results in the expression of a combination of six PLA(2) species--ammodytin I1, ammodytin I2, ammodytin L, ammodytoxin, vaspin A and vaspin B--with 19 known isoforms of the first five of these species. Most of the European viper venoms studied contained either a myotoxin or a neurotoxin, and all contained ammodytin I1 and ammodytin I2. There is no evidence that a given pattern of PLA(2) species constitutes a taxonomic criterion, and isoform analysis would be required for such discrimination. Analysis of the phylogenetic relationships between PLA(2) species from European vipers and those of other members of the Viperinae revealed a strong correlation between the geographical source of the viper and the clustering seen for the different isoforms, for each PLA(2) species. The K(a)/K(s) values calculated for the mature protein-coding region of paralogous genes showed that ratios for pairs including vaspin B or one ammodytoxin isoform were greater than 1.09, whereas those for most of the remaining pairs were less than 1. Different patterns of mutation were observed in comparisons of the different PLA(2) isoforms. The mechanisms directing a mutation toward a precise exon remain unresolved.

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.

Venom phospholipases of Russell's vipers from Myanmar and eastern India--cloning, characterization and phylogeographic analysis

Venoms of Russell's vipers (genus Daboia) are known for their deadly coagulopathic and other effects. We herein studied various isoforms of venom phospholipases A(2) (PLAs) from two Daboia species at their geographic boundary. From Myanmar Daboia siamensis venom (designated as DsM), four PLAs (designated DsM-aI, aI', aII' and bI') were purified, and the cDNAs encoding two acidic (DsM-aI and aII) and two basic PLAs (DsM-bI and S1) were also cloned from its venom-glands. DsM-S1 is identical to the major venom PLA of southern India Daboia russelii, but the protein is absent from the venom. Additionally, four PLAs (designated DrK-aI, aII, bI and bII) were cloned from cDNA obtained from venom glands of a Kolkata D. russelii, and the PLAs were purified from the pooled venom (designated as DrK). The acidic DrK-aI is the most neurotoxic and lethal among these PLAs; DsM-aI which differs from DrK-aI by only the Phe2 substitution shows greatly reduced enzymatic activity and lethality. Both acidic PLAs do not form dimeric complex with basic PLAs in the same venoms. DsM-bI' is neurotoxic and lethal but its orthologous DrK-bI (97% identical to DsM-bI') is a much weaker toxin. Given the fact that most of the orthologous PLAs of DrK and DsM share 97-100% sequence identity, Daboia vipers of Myanmar and Kolkata must be closely related. Molecular phylogenetic analyses on 30 venom PLAs of Eurasian vipers' revealed co-evolution of five subtypes of venom PLAs in both Daboia and Vipera genera. Our results shed light on the intra- and inter-species variations and structure-function relationships of viperid venom PLAs.

Sequence analysis and phylogenetic relationship of genes encoding heterodimeric phospholipases A2 from the venom of the scorpion Anuroctonus phaiodactylus

Some scorpion venom contain heterodimeric phospholipases A2. They were shown to be toxic to insects and to cause edema and/or hemolysis of mammalian erythrocytes. This manuscript describes the results of cDNA cloning of five different heterodimeric phospholipases from the venomous glands of the Mexican scorpion Anuroctonus phaiodactylus. The amino acid sequence deduced from the heterodimeric phospholipases open reading frames corresponds in each case to a different isoform. The nucleotide sequences corresponding to two of these genes were also obtained by directly sequencing genomic DNA. The cDNA isoforms show high similarity with the heterodimeric phospholipase Phaiodactylipin purified from the same scorpion. However, similar phospholipases were also found in scorpions from other species and the sequences available were used to construct a phylogenetic tree. In order to understand better the gene structure and phylogeny of these enzymes we analyzed their sequences and compared them with secretory phospholipases of other sources from groups I, II and III. The genomic DNA sequence of a similar phospholipase from bee venomous glands was also cloned. The information available on a Drosophila phospholipase was included in this analysis. The phospholipases of groups I and II contain a conserved exon-intron structure (four or five exons of the mature segment of the enzyme are separated by three or four introns). Also, the gene structure of the phospholipases from A. phaiodactylus and that of the bee venom, belonging to group III phospholipases, are interrupted by three introns. The mature peptide of the bee enzyme is a single polypeptide chain, coded by four exons, whereas those from the scorpion studied here although having four exons, showed the presence of two different polypeptides in its native state. The mature protein is processed after synthesis, producing the heterodimeric structure: a long and a short-peptide chain, linked by a disulfide bridge. The small subunit is the one coded by the fourth exon. The human phospholipase A2 and that of Drosophila, also classified into the group III phospholipases, have a considerably different exon-intron organization.

Influence of phospholipases A2 from snake venoms on survival and neurite outgrowth in pheochromocytoma cell line PC12

To determine whether the ability to induce neurite outgrowth in rat pheochromocytoma cell line PC12 is characteristic of phospholipases of different types, we have studied the influence of phospholipase A(2) (PLA2) from cobra Naja kaouthia venom and two PLA2s from viper Vipera nikolskii venom on PC12 cells. Phospholipases from the viper venom are heterodimers in which only one of the subunits is enzymatically active, while PLA2 from the cobra venom is a monomer. It was found that all three PLA2s induce neurite outgrowth in PC12. The PLA2 from cobra venom exhibits this effect at higher concentrations as compared to the viper enzymes. We have not observed such an activity for isolated subunits of viper PLA2s, since the enzymatically active subunits have very high cytotoxicity, while the other subunits are not active at all. However, co-incubation of active and inactive subunits before addition to the cells leads to a marked decrease in cytotoxicity and to restoration of the neurite-inducing activity. It has also been shown that all enzymatically active PLA2s are cytotoxic, the PLA2 from cobra venom being the least active. Thus, for the first time we have shown that PLA2s from snake venoms can induce neurite outgrowth in PC12 cells.

Characterisation of some pharmacological effects of the venom from Vipera lebetina

Vipera lebetina is one of the most venomous snakes on the Iran plateau. Serious clinical problems such as edema, hemorrhage and tissue necrosis are observed in humans following V. lebetina envenomating. However, little information on the pharmacological properties of the venom is available. To determine haemodynamic actions of the venom of V. lebetina, the changes in the mean arterial blood pressure of anaesthetised rats following the administration of the venom were recorded. Venom (1 mg/kg, i.v.) produced rapid cardiovascular collapse, while 0.3 mg/kg (i.v.) caused only a small transient decrease in mean arterial blood pressure. Effects of the venom on perfusion pressure in the isolated rat mesenteric bed, and on contractions of the isolated rat right atrium and the isolated guinea-pig ileum, were also studied. Exposure of the isolated rat right atrium to venom (0.1-1 mg/ml) caused a transient increase followed by a sustained reduction in the amplitude and frequency of spontaneous contractions. The transient positive inotropic and chronotropic effects were abolished when the preparation was preincubated with propranolol, but not with tolazoline. N(G)-nitro-l-arginine methyl ester pretreatment attenuated the vascular hyporeactivity to phenylephrine induced by the venom in the isolated rat mesenteric vascular bed. This suggests that nitric oxide (NO) or NO-like compounds may be present in the venom and involved in its hypotensive effect. The venom (0.3-1 mg/ml) caused concentration-dependant blockade of isolated guinea-pig ileum contractions induced by electrical field stimulation, acetylcholine or KCl. This inhibitory effect of the venom was significantly reduced by prior incubation of the venom with manoalide (1 microM) indicating involvement of a phospholipase A(2) component. Further, investigation is required to identify specific toxins responsible for the above pharmacological effects.