Biochemical composition, lethality and pathophysiology of venom from two cobras-- Naja naja and N. kaouthia

The cobras Naja naja and N. kaouthia are abundant in eastern and north-eastern India, accounting for maximum snakebite deaths. Here we report on variation in the composition of Naja kaouthia and N. naja venom from eastern India on corresponding differences in the severity of pathogenesis. These two venoms differ in chromatographic elution profile through Sephadex G-50 and enzyme activity, protein and carbohydrate contents associated with each fraction. The presence of greater amounts of basic phospholipase A2, L-amino acid oxidase and low molecular weight membrane active polypeptides in the N. naja venom makes it more toxic than N. kaouthia venom. A commercial polyvalent antivenom raised against N. naja venom inactivates lethality and variety of toxic effects of homologous venom more effectively than N. kaouthia venom.     

Comparative characterisation of Russell's viper (Daboia/Vipera russelli) venoms from different regions of the Indian peninsula.

Russell's viper (Daboia/Vipera russelli) venom from different regions of India was subjected to chromatographic, electrophoretic, biochemical and immunological analysis. The elution profiles from ion-exchange chromatography and protein banding pattern from SDS-PAGE showed a significant variation in the constituents of venoms. The acidic proteins are found to be predominant in the venoms of eastern and western regions while basic proteins are the major contributors of the northern and southern regional venoms. The major variation of phospholipases A(2) in the venom samples of India may be described as: southern regional venom is rich in basic, toxic PLA(2) while this activity showed a dramatic decrease as one moves towards west, north and eastern regions of India. In addition, the caseinolytic, TAME-hydrolytic, anticoagulant, oedema-inducing and haemorrhagic activities of the venoms have also varied from one region to another. The muscle specimens of mice injected with venoms of different regions showed variable change in the muscle fibre damage and cell morphology. The eastern regional venom is most lethal among all the venoms. The lethal potencies for four regional venoms vary as: eastern>western>southern>northern. The polyclonal antibodies prepared against the venom of southern region showed cross-reaction with the venoms of other regions, but the extent of cross-reaction and diffusion patterns are different. However, the polyclonal antibodies prepared against southern regional venom showed no protection against lethal toxicity of other regional venoms.     

Purification,Characterization, and Chemical Modification of Neurotoxic Peptide from Daboia russelii Snake Venom of India

Comprehensive knowledge of venom composition is very important for effective management of snake envenomation and antivenom preparation. Daboia russelii venom from the eastern region of India is the most neurotoxic among the four venom samples investigated. From the eastern D. russelii venom sample, neurotoxic peptide has been purified by combined method of ion exchange gel permeation chromatography and reversed phase high performance liquid chromatography. Molecular weight of Daboia neurotoxin III (DNTx‐III) found to be 6,849 Da (as measured on matrix‐assisted laser desorption/ionisation‐time of flight mass spectrometer), and N‐terminal amino acid sequences is I K C F I T P D U T S Q A. Approximate LD₅₀ dosage was 0.24 mg/kg body weight. It produced concentration‐ and time‐dependent inhibition of indirectly stimulated twitches of Rana hexadactyla sciatic nerve gastrocnemius muscle preparations. Chemical modification of DNTx‐III tryptophan residue(s) reduced the twitch height inhibition property of toxin, signifying the importance of tryptophan residues for the neurotoxic function. This type of neurotoxic peptide is unique to east Indian regional D. russelii venom. © 2013 Wiley Periodicals, Inc. J BiochemMol Toxicol 27:295‐304, 2013; View this article online at wileyonlinelibrary.com . DOI 10.1002/jbt.21486     

Deconstructing a Complex Molecular Phenotype: Population-Level Variation in Individual Venom Proteins in Eastern Massasauga Rattlesnakes (<Emphasis Type="Italic">Sistrurus c. catenatus</Emphasis>)

Identifying the molecular basis for complex adaptations such as the toxic proteins used by venomous snakes to subdue and digest prey is an important step in understanding the evolutionary and functional basis for such traits. Recent proteomics-based analyses have made possible the identification of all constituent proteins in whole venom samples. Here we exploit this advance to study patterns of population-level variation in venom proteins from 254 adult eastern massasauga rattlesnakes (Sistrurus c. catenatus) collected from 10 populations. Analysis of presence–absence variation in specific proteins from 1D PAGE gels shows that: (1) The frequency spectra for individual protein bands is U-shaped with a large number of specific proteins either being consistently “common” or “rare” across populations possibly reflecting functional differences. (2) Multivariate axes which summarize whole venom variation consist of bands from all major types of proteins implying the integration of functionally distinct components within the overall venom phenotype. (3) There is significant differentiation in venom proteins across populations and the specific classes of proteins contributing to this differentiation have been identified. (4) Levels of population differentiation in venom proteins are not correlated with levels of neutral genetic differentiation, or genetically effective population sizes which argues that patterns of venom variation are not simply a consequence of population structure but leaves open the role of selection in generating population differences in venom. Our results identify particular classes of venom proteins and their associated genes as being fruitful targets for future studies of the molecular and functional basis for this complex adaptive phenotype.     

Variation 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.     

Population systematics of Russell's viper: a multivariate study

A multivariate analysis of the population systematics of Russell's viper, based on scalation and colour pattern characters, reveals that the populations of this viper constitute two well-defined taxa: a western form, comprising all populations from the Indian subcontinent, and an eastern form, comprising all populations from east of the Bay of Bengal. The two forms could be considered either as subspecies of one species, or as two separate species, depending on the species concept used. Within the western form, there is no clear pattern of geographic variation. Within the eastern form, the populations from the Lesser Sunda Islands are clearly divergent from the populations of mainland Asia and Java. The conventionally recognized subspecies of Vipera russelli fail to portray this pattern of geographic variation. There is no clear relationship between the pattern of geographic variation in morphology and the pattern of geographic variation in the clinical effects of the venom in human bite victims: some populations with considerable differences in venom effects are equally distinct morphologically, whereas other populations with equally strong venom differences are morphologically very similar. The distribution of Russell's viper can be attributed to Pleistocene changes in climate and sea level, coupled with the viper's ecological requirements, which appear to include a seasonally dry climate.     

Maternal footprints of Southeast Asians in North India

We have analyzed 7,137 samples from 125 different caste, tribal and religious groups of India and 99 samples from three populations of Nepal for the length variation in the COII/tRNA(Lys) region of mtDNA. Samples showing length variation were subjected to detailed phylogenetic analysis based on HVS-I and informative coding region sequence variation. The overall frequencies of the 9-bp deletion and insertion variants in South Asia were 1.9 and 0.6%, respectively. We have also defined a novel deep-rooting haplogroup M43 and identified the rare haplogroup H14 in Indian populations carrying the 9-bp deletion by complete mtDNA sequencing. Moreover, we redefined haplogroup M6 and dissected it into two well-defined subclades. The presence of haplogroups F1 and B5a in Uttar Pradesh suggests minor maternal contribution from Southeast Asia to Northern India. The occurrence of haplogroup F1 in the Nepalese sample implies that Nepal might have served as a bridge for the flow of eastern lineages to India. The presence of R6 in the Nepalese, on the other hand, suggests that the gene flow between India and Nepal has been reciprocal.     

P-III hemorrhagic metalloproteinases from Russell's viper venom: cloning, characterization, phylogenetic and functional site analyses

Two homologous P-III hemorrhagic metalloproteinases were purified from Russell's viper venoms from Myanmar and Kolkata (eastern India), and designated as daborhagin-M and daborhagin-K, respectively. They induced severe dermal hemorrhage in mice at a minimum hemorrhagic dose of 0.8-0.9 microg. Daborhagin-M specifically hydrolyzed an Aalpha-chain of fibrinogen, fibronectin, and type IV collagen in vitro. Analyses of its cleavage sites on insulin chain B and kinetic specificities toward oligopeptides suggested that daborhagin-M prefers hydrophobic residues at the P(1), P(1)', and P(2)' positions on the substrates. Of the eight Daboia geographic venom samples analyzed by Western blotting, only those from Myanmar and eastern India showed a strong positive band at 65kDa, which correlated with the high risk of systemic hemorrhagic symptoms elicited by Daboia envenoming in both regions. The full sequence of daborhagin-K was determined by cDNA cloning and sequencing, and then confirmed by peptide mass fingerprinting. Furthermore, molecular phylogenetic analyses based on 27 P-IIIs revealed the co-evolution of two major P-III classes with distinct hemorrhagic potencies, and daborhagin-K belongs to the most hemorrhagic subclass. By comparing the absolute complexity profiles between these two classes, we identified four structural motifs probably responsible for the phylogenetic subtyping and hemorrhagic potencies of P-III SVMPs.     

Proteomic analysis of ontogenetic and diet-related changes in venom composition of juvenile and adult Dusky Pigmy rattlesnakes (Sistrurus miliarius barbouri)

Snake venom proteins show high levels of variation at the level of the individual yet the environmental and molecular mechanisms that generate this diversity remain unclear. Here we report the results of a controlled feeding experiment combined with proteomic analyses of periodically collected venom samples to assess the roles of ontogenetic and diet-related effects on venom composition of captive juvenile and adult Dusky Pigmy rattlesnakes (Sistrurus miliarius barbouri). Juvenile snakes fed from birth with mice, lizards, or frogs showed little evidence for an ontogenetic shift in venom composition from 5 to 26 months in terms of substantial changes in the relative abundance of major classes of venom toxins. However, there were fine-scale changes in the relative abundance of D49-PLA? 15, PI-SVMPs, and PIII-SVMP 28, and a decline in the abundance of other PIII-SVMPs. Although juveniles raised on different prey exhibited distinct relative toxin compositional change rates, at 26 months old, their venoms showed similar patterns of protein composition suggesting little effect of diet on the overall make-up of venom in snakes this age or younger. In contrast, adult females raised on different prey over a 26 month period showed prey-related changes in the relative abundance of major protein families from initial to final samples. Specifically, mouse-fed females showed substantial increases in the relative abundance of total PLA?s and serine proteinases of 95% and >100%, respectively, whereas comparable values for lizard- (42% and -22%) and frog-fed females (2% and 11%) were distinctly smaller in magnitude. Venom from adult snakes fed on different prey also showed distinct changes in the abundance of PLA? molecules 15, 19a, and 19, which were, respectively, (>100%, 33%, 63%), (>100%, 0%, 35%), and (71%, 20%, -4%) for the mice-, lizard-, and frog-diet. Venom from snakes raised on frogs contained a small (1.1%) but consistent amount of a PLA? molecule (15a) not present in snakes fed on mice or lizards. This work provides evidence that venom composition is somewhat plastic in both juvenile and adult S. m. barbouri and that, at least in adults, prey consumed may influence the relative abundance of possibly functionally-distinct classes of venom proteins.     

Genetic and Ecological Correlates of Intraspecific Variation in Pitviper Venom Composition Detected Using Matrix-Assisted Laser Desorption Time-of-Flight Mass Spectrometry (MALDI-TOF-MS) and Isoelectric Focusing

The ability to detect biochemical diversity in animal venoms has wide-ranging implications for a diverse array of scientific disciplines. Matrix-assisted laser desorption time-of-flight mass spectrometry (and, for comparative purposes, isoelectric focusing) were used to characterize venoms from a geographically diverse sample of Trimeresurus stejnegeri ( n < 229) from Taiwan. Previously unrealized levels of heterogeneity were detected in venom phospholipase A(2) isoforms (PLA(2)) and in whole venom profiles. Geographic variation in venom was primarily between Taiwan and two Pacific islets. Despite the common assumption that venom variation is a product of neutral molecular evolution, statistical testing failed to link venom variation with phylogenetic descent convincingly. Instead, pronounced differences in venom composition may be the product of natural selection for regional diets or of independent founder effects. More data are required on the functional differences between the isoforms to distinguish between these alternatives.     

Regional variation of biochemical characteristics and antigeneity in Great Basin rattlesnake (Crotalus viridis lutosus) venom

1. Three pooled and 20 individual venom samples of Crotalus viridis lutosus from different localities in Utah and Arizona were screened and fractionated with HPLC-anion exchange. 2. Pooled venom samples and fractions were tested for hemorrhagic, collagenase, and phospholipase activities, and reactivity to a monoclonal antibody against a hemorrhagin from C. atrox venom (CA-P-8) using ELISA. 3. The 20 individual samples were organized into four groups based on their HPLC profiles. 4. ELISA results and specific hemorrhagic activity of the venom samples displayed a variation in latitidinal distribution although from the same species.     

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₅₀) 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₅₀, 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.     

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.     

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age‐specific selection is constrained to act on the same genome by altering age‐specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype–phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.     

Biogeographical venom variation in the Indian spectacled cobra (Naja naja) underscores the pressing need for pan-India efficacious snakebite therapy

BackgroundSnake venom composition is dictated by various ecological and environmental factors, and can exhibit dramatic variation across geographically disparate populations of the same species. This molecular diversity can undermine the efficacy of snakebite treatments, as antivenoms produced against venom from one population may fail to neutralise others. India is the world’s snakebite hotspot, with 58,000 fatalities and 140,000 morbidities occurring annually. Spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are known to cause the majority of these envenomations, in part due to their near country-wide distributions. However, the impact of differing ecologies and environment on their venom compositions has not been comprehensively studied.MethodsHere, we used a multi-disciplinary approach consisting of venom proteomics, biochemical and pharmacological analyses, and in vivo research to comparatively analyse N. naja venoms across a broad region (>6000 km; seven populations) covering India’s six distinct biogeographical zones.FindingsBy generating the most comprehensive pan-Indian proteomic and toxicity profiles to date, we unveil considerable differences in the composition, pharmacological effects and potencies of geographically-distinct venoms from this species and, through the use of immunological assays and preclinical experiments, demonstrate alarming repercussions on antivenom therapy. We find that commercially-available antivenom fails to effectively neutralise envenomations by the pan-Indian populations of N. naja, including a complete lack of neutralisation against the desert Naja population.ConclusionOur findings highlight the significant influence of ecology and environment on snake venom composition and potency, and stress the pressing need to innovate pan-India effective antivenoms to safeguard the lives, limbs and livelihoods of the country’s 200,000 annual snakebite victims.     

Identifying Intraspecific Variation in Venom Yield of Chinese Cobra(Naja atra) from Ten Populations in Mainland China

Detailed information on venom yield is helpful in preparing antivenoms and treating snakebites, but such information is lacking for many species of venomous snakes. The Chinese cobra(Naja atra) is a large sized, venomous snake commonly found in southeastern China, where it causes a heavy burden of snakebites. To examine the effects of various factors(morphology, sex, age, season, and geographical origin) on the venom yield in this snake, we collected venom samples of 446 individuals(426 adults and 20 neonates) from 10 populations of N. atra over an eightyear period. We used two variables, lyophilized venom mass(venom yield) and solid content of venom(% solids), to quantify the venom yield. We used linear regression analysis to check if venom yield was related to morphological factors, one-way ANOVA and one-way ANCOVA to detect the sexual, ontogenetic, and geographic variation in venom yield, and repeated-measures ANOVA to examine seasonal shifts in venom yield. Our results indicate that venom yield of N. atra is positively related to the morphological traits examined, with male snakes expelling more venom than females. Venom yield in N. atra was age-related, with elder snakes always expelling more venom than younger ones. Geographic variation in venom yield was also observed, while seasonal variation was not. The solid content of venom was lower in males than in females, but this was not related to morphology, season, age, or geography. Our findings suggest that venom yield in N. atra is influenced by multiple factors, as well as by the interactions among these factors.     

Proteomics analysis to compare the venom composition between Naja naja and Naja kaouthia from the same geographical location of eastern India: Correlation with pathophysiology of envenomation and immunological cross-reactivity towards commercial polyantivenom

ABSTRACT

Background : Cobra bite is frequently reported across the Indian subcontinent and is associated with a high rate of death and morbidity. In eastern India (EI) Naja naja and Naja kaouthia are reported to be the two most abundant species of cobra.

Research design and methods : The venom proteome composition of N. naja (NnV) and N. kaouthia (NkV) from Burdwan districts of EI were compared by separation of venom proteins by 1D-SDS-PAGE followed by LC-MS/MS analysis of protein bands. The potency of commercial polyantivenom (PAV) was assessed by neutralization, ELISA, immuno-blot and venom-PAV immunoaffinity chromatography studies.

Results : Proteomic analysis identified 52 and 55 proteins for NnV and NkV, respectively, when searched against the Elapidae database. A small quantitative difference in venom composition between these two species of cobra was observed. PAVs exhibited poor cross-reactivity against low molecular mass toxins (<20 kDa) of both cobra venoms, which was substantiated by a meager neutralization of their phospholipase A₂ activity. Phospholipase A₂ and 3FTx, the two major classes of nonenzymatic and enzymatic proteins, respectively, were partially recognized by PAVs.

Conclusions : Efforts must be made to improve immunization protocols and supplement existing antivenoms with antibodies raised against the major toxins of these venoms.     

Alteration of the protein composition of bothrops jararaca venom and venom gland by isoproterenol treatment

• 1.1. The protein composition of Bothrops jararaca venom and venom gland was analyzed through SDS-PAGE, after isoproterenol (IPR) treatment.
• 2.2. Some proteins (47, 48, 57 and 72 kDa) were detected in the gland homogenate from the control but not from the IPR-treated samples.
• 3.3. Three proteins (26.5, 44.5 and 53 kDa) were detected in the venom gland from IPR-treated snakes but not from the venom gland from the control.
• 4.4. In the venom samples proteins of 41 and 74 kDa were detected only in the IPR treated samples, while proteins of 17 and 28 kDa were detected only in the control.
• 5.5. The biological activity of the venom did not change with IPR treatment.

     

Phylogenetic affinities of the Iraqi populations of Saw-scaled vipers of the genus Echis (Serpentes: Viperidae), revealed by sequences of mtDNA genes

The Saw-scaled vipers of the species Echis carinatus range from Sri Lanka and India westwards to Iraq, including the eastern Arabian Peninsula. We collected the species in southern Iraq and compared two mtDNA genes (16S and Cyt b) with other populations of this species and with other species of the genus. Analyses of both Maximum Likelihood and Bayesian Inference confirmed E. carinatus as a monophyletic species. The samples from Iraq cluster with populations of the species from Pakistan and UAE. Populations from India, however, are situated in a separate phylogenetic lineage. This can be explained by the geographic barriers between western (Iraq, Pakistan and UAE) and eastern (India) populations of the species. Soleyman Mountain in southern Pakistan is the main barrier between them and its role is reflected in the genetic distance between populations.