Sequences, geographic variations and molecular phylogeny of venom phospholipases and threefinger toxins of eastern India Bungarus fasciatus and kinetic analyses of its Pro31 phospholipases A2

Identification of poorly immunodepleted phospholipase A2 (PLA2) proteins of Bungarus fasciatus venom from Assam, India and evaluation of Indian polyvalent antivenom using third-generation antivenomics

Bungarus fasciatus also referred to as the Banded krait is a snake which possesses venom and belongs to the Elapidae family. It is widely distributed across the Indian subcontinent and South East Asian countries and is responsible for numerous snakebites in the population. B. fasciatus possesses a neurotoxic venom and envenomation by the snake results in significant morbidity and occasional morbidity in the victim if not treated appropriately. In this study, the efficacy of Indian polyvalent antivenom (Premium Serums polyvalent antivenom) was evaluated against the venom of B. fasciatus from Guwahati, Assam (India) employing the Third-generation antivenomics technique followed by identification of venom proteins from three poorly immunodepleted peaks (P5, P6 and P7) using LC-MS/MS analysis. Seven proteins were identified from the three peaks and all these venom proteins belonged to the phospholipase A2 (PLA2) superfamily. The identified PLA2 proteins were corroborated by the in vitro enzymatic activities (PLA2 and Anticoagulant activity) exhibited by the three peaks and previous reports of pathological manifestation in the envenomated victims. Neutralization of enzymatic activities by Premium Serums polyvalent antivenom was also assessed in vitro for crude venom, P5, P6 and P7 which revealed moderate to poor inhibition. Inclusion of venom proteins/peptides, which are non-immunodepleted or poorly immunodepleted, into the immunization mixture of venom used for antivenom production may help in enhancing the efficacy of the polyvalent antivenom.

Molecular phylogeny reveals distinct evolutionary lineages of the banded krait, Bungarus fasciatus (Squamata, Elapidae) in Asia

The banded krait, Bungarus fasciatus is a widespread elapid snake, likely to comprise several distinct species in different geographic regions of Asia. Therefore, based on molecular phylogenetics and comparative morphology data, we present an overview of the systematic composition of the species to delimit potential biogeographic boundaries. Our phylogenetic analyses, based on four mitochondrial genes, reveal the existence of at least three evolutionary lineages within B. fasciatus, corresponding to Indo-Myanmar, Sundaic and eastern Asian lineages. We are convinced that there are at least three taxonomic entities within the nomen B. fasciatus and restrict the distribution of B. fasciatus sensu stricto to the Indo-Myanmar region. We also provide additional natural history data of the taxon from eastern India. Finally, we advocate further studies to establish the degree of reproductive isolation among these diverging evolutionary lineages and to reassess the systematic status of this species complex especially the Sundaic and eastern Asian lineages.

De novo venom gland transcriptomics of Calliophis bivirgata flaviceps: uncovering the complexity of toxins from the Malayan blue coral snake

Background:

The Malayan blue coral snake, Calliophis bivirgata flaviceps, is a medically important venomous snake in Southeast Asia. However, the complexity and diversity of its venom genes remain little explored.

Methods:

To address this, we applied high-throughput next-generation sequencing to profile the venom gland cDNA libraries of C. bivirgata flaviceps. The transcriptome was de novo assembled, followed by gene annotation, multiple sequence alignment and analyses of the transcripts.

Results:

A total of 74 non-redundant toxin-encoding genes from 16 protein families were identified, with 31 full-length toxin transcripts. Three-finger toxins (3FTx), primarily delta-neurotoxins and cardiotoxin-like/cytotoxin-like proteins, were the most diverse and abundantly expressed. The major 3FTx (Cb_FTX01 and Cb_FTX02) are highly similar to calliotoxin, a delta-neurotoxin previously reported in the venom of C. bivirgata. This study also revealed a conserved tyrosine residue at position 4 of the cardiotoxin-like/cytotoxin-like protein genes in the species. These variants, proposed as Y-type CTX-like proteins, are similar to the H-type CTX from cobras. The substitution is conservative though, preserving a less toxic form of elapid CTX-like protein, as indicated by the lack of venom cytotoxicity in previous laboratory and clinical findings. The ecological role of these toxins, however, remains unclear. The study also uncovered unique transcripts that belong to phospholipase A₂ of Groups IA and IB, and snake venom metalloproteinases of PIII subclass, which show sequence variations from those of Asiatic elapids.

Conclusion:

The venom gland transcriptome of C. bivirgata flaviceps from Malaysia was de novo assembled and annotated. The diversity and expression profile of toxin genes provide insights into the biological and medical importance of the species.

Convergent evolution of pain-inducing defensive venom components in spitting cobras

Convergent evolution provides insights into the selective drivers underlying evolutionary change. Snake venoms, with a direct genetic basis and clearly defined functional phenotype, provide a model system for exploring the repeated evolution of adaptations. While snakes use venom primarily for predation, and venom composition often reflects diet specificity, three lineages of cobras have independently evolved the ability to spit venom at adversaries. Using gene, protein, and functional analyses, we show that the three spitting lineages possess venoms characterized by an up-regulation of phospholipase A2 (PLA2) toxins, which potentiate the action of preexisting venom cytotoxins to activate mammalian sensory neurons and cause enhanced pain. These repeated independent changes provide a fascinating example of convergent evolution across multiple phenotypic levels driven by selection for defense.

Variations in neurotoxicity and proteome profile of Malayan krait (Bungarus candidus) venoms

Malayan krait (Bungarus candidus) is a medically important snake species found in Southeast Asia. The neurotoxic effects of envenoming present as flaccid paralysis of skeletal muscles. It is unclear whether geographical variation in venom composition plays a significant role in the degree of clinical neurotoxicity. In this study, the effects of geographical variation on neurotoxicity and venom composition of B. candidus venoms from Indonesia, Malaysia and Thailand were examined. In the chick biventer cervicis nerve-muscle preparation, all venoms abolished indirect twitches and attenuated contractile responses to nicotinic receptor agonists, with venom from Indonesia displaying the most rapid neurotoxicity. A proteomic analysis indicated that three finger toxins (3FTx), phospholipase A2 (PLA2) and Kunitz-type serine protease inhibitors were common toxin groups in the venoms. In addition, venom from Thailand contained L-amino acid oxidase (LAAO), cysteine rich secretory protein (CRISP), thrombin-like enzyme (TLE) and snake venom metalloproteinase (SVMP). Short-chain post-synaptic neurotoxins were not detected in any of the venoms. The largest quantity of long-chain post-synaptic neurotoxins and non-conventional toxins was found in the venom from Thailand. Analysis of PLA2 activity did not show any correlation between the amount of PLA2 and the degree of neurotoxicity of the venoms. Our study shows that variation in venom composition is not limited to the degree of neurotoxicity. This investigation provides additional insights into the geographical differences in venom composition and provides information that could be used to improve the management of Malayan krait envenoming in Southeast Asia.

Phospholipase A2 from krait Bungarus fasciatus venom induces human cancer cell death in vitro

Background

Snake venoms are the complex mixtures of different compounds manifesting a wide array of biological activities. The venoms of kraits (genus Bungarus, family Elapidae) induce mainly neurological symptoms; however, these venoms show a cytotoxicity against cancer cells as well. This study was conducted to identify in Bungarus fasciatus venom an active compound(s) exerting cytotoxic effects toward MCF7 human breast cancer cells and A549 human lung cancer cells.

Methods

The crude venom of B. fasciatus was separated by gel-filtration on Superdex HR 75 column and reversed phase HPLC on C18 column. The fractions obtained were screened for cytotoxic effect against MCF7, A549, and HK2 cell lines using colorimetric assay with the tetrazolium dye MTT- 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide. The primary structure of active protein was established by ultra high resolution LC-MS/MS. The molecular mechanism of the isolated protein action on MCF7 cells was elucidated by flow cytometry.

Results

MTT cell viability assays of cancer cells incubated with fractions isolated from B. fasciatus venom revealed a protein with molecular mass of about 13 kDa possessing significant cytotoxicity. This protein manifested the dose and time dependent cytotoxicity for MCF7 and A549 cell lines while showed no toxic effect on human normal kidney HK2 cells. In MCF7, flow cytometry analysis revealed a decrease in the proportion of Ki-67 positive cells. As Ki-67 protein is a cellular marker for proliferation, its decline indicates the reduction in the proliferation of MCF7 cells treated with the protein. Flow cytometry analysis of MCF7 cells stained with propidium iodide and Annexin V conjugated with allophycocyanin showed that a probable mechanism of cell death is apoptosis. Mass spectrometric studies showed that the cytotoxic protein was phospholipase A₂. The amino acid sequence of this enzyme earlier was deduced from cloned cDNA, and in this work it was isolated from the venom as a protein for the first time. It is also the first krait phospholipase A₂ manifesting the cytotoxicity for cancer cells.

Therapeutic potential of krait venom

Kraits belong to Elapideae and are widely distributed in East and South-East Asian countries. Krait venom possesses neurotoxins, membrane toxins, cardiotoxins, three finger toxins, metalloproteinases, cholinesterases, L-amino acid oxidases and serine proteases. The therapeutic potential of krait venom in pathophysiological conditions such as microbial and parasitic infections, cancer, arthritis, inflammation and blood coagulation disorder is discussed in this review. More intensive new research ventures are required to establish the therapeutic potential of krait venom in complex and emerging diseases.

Inhibition of nicotinic acetylcholine receptors, a novel facet in the pleiotropic activities of snake venom phospholipases A2

Phospholipases A2 represent the most abundant family of snake venom proteins. They manifest an array of biological activities, which is constantly expanding. We have recently shown that a protein bitanarin, isolated from the venom of the puff adder Bitis arietans and possessing high phospholipolytic activity, interacts with different types of nicotinic acetylcholine receptors and with the acetylcholine-binding protein. To check if this property is characteristic to all venom phospholipases A2, we have studied the capability of these enzymes from other snakes to block the responses of Lymnaea stagnalis neurons to acetylcholine or cytisine and to inhibit α-bungarotoxin binding to nicotinic acetylcholine receptors and acetylcholine-binding proteins. Here we present the evidence that phospholipases A2 from venoms of vipers Vipera ursinii and V. nikolskii, cobra Naja kaouthia, and krait Bungarus fasciatus from different snake families suppress the acetylcholine- or cytisine-elicited currents in L. stagnalis neurons and compete with α-bungarotoxin for binding to muscle- and neuronal α7-types of nicotinic acetylcholine receptor, as well as to acetylcholine-binding proteins. As the phospholipase A2 content in venoms is quite high, under some conditions the activity found may contribute to the deleterious venom effects. The results obtained suggest that the ability to interact with nicotinic acetylcholine receptors may be a general property of snake venom phospholipases A2, which add a new target to the numerous activities of these enzymes.

Proteomic characterization and comparison of Malaysian Bungarus candidus and Bungarus fasciatus venoms

Kraits (Bungarus spp.) are highly venomous elapids that are only found in Asia. In the current study, 103 and 86 different proteins were identified from Bungarus candidus and Bungarus fasciatus venoms, respectively. These proteins were classified into 18 different venom protein families. Both venoms were found to contain a high percentage of three finger toxins, phospholipase A2 enzymes and Kunitz-type inhibitors. Smaller number of high molecular weight enzymes such as L-amino acid oxidase, hyaluronidases, and acetylcholinesterase were also detected in the venoms. We also detected some unique proteins that were not known to be present in these venoms. The presence of a natriuretic peptide, vespryn, and serine protease families was detected in B. candidus venom. We also detected the presence of subunit A and B of β-bungarotoxin and α-bungarotoxin which had not been previously found in B. fasciatus venom. Understanding the proteome composition of Malaysian krait species will provide useful information on unique toxins and proteins which are present in the venoms. This knowledge will assist in the management of krait envenoming. In addition, these proteins may have potential use as research tools or as drug-design templates.

BIOLOGICAL SIGNIFICANCE

This study has revealed the proteome composition of Malaysian B. candidus and B. fasciatus venoms, two medically important snake species in Asia. Information on the venom proteome of these species will provide useful information for krait bite management and aid in antivenom selection. Venom proteome profiles of these venoms showed that there are significant differences in the venom protein family compositions. Detection of proteins and peptides that have not been documented in these species such as natriuretic peptides, vespryn and serine proteases provides new knowledge on the composition of these venoms. The roles of these new proteins and peptides in krait envenoming are still unknown. Discovery of these proteins and peptides may also be useful for future research tool and therapeutic development.

Isolation and characterization of a presynaptic neurotoxin, P-elapitoxin-Bf1a from Malaysian Bungarus fasciatus venom

Presynaptic neurotoxins are one of the major components in Bungarus venom. Unlike other Bungarus species that have been studied, β-bungarotoxin has never been isolated from Bungarus fasciatus venom. It was hypothesized that the absence of β-bungarotoxin in this species was due to divergence during evolution prior to evolution of β-bungarotoxin. In this study, we have isolated a β-bungarotoxin isoform we named P-elapitoxin-Bf1a by using gel filtration, cation-exchange and reverse-phase chromatography from Malaysian B. fasciatus venom. The toxin consists of two heterogeneous subunits, subunit A and subunit B. LCMS/MS data showed that subunit A was homologous to acidic phospholipase A2 subunit A3 from Bungarus candidus and B. multicinctus venoms, whereas subunit B was homologous with subunit B1 from B. fasciatus venom that was previously detected by cDNA cloning. The toxin showed concentration- and time-dependent reduction of indirect-twitches without affecting contractile responses to ACh, CCh or KCl at the end of experiment in the chick biventer preparation. Toxin modification with 4-BPB inhibited the neurotoxic effect suggesting the importance of His-48. Tissue pre-incubation with monovalent B. fasciatus (BFAV) or neuro-polyvalent antivenom (NPV), at the recommended titer, was unable to inhibit the twitch reduction induced by the toxin. This study indicates that Malaysian B. fasciatus venom has a unique β-bungarotoxin isoform which was not neutralized by antivenoms. This suggests that there might be other presynaptic neurotoxins present in the venom and there is a variation in the enzymatic neurotoxin composition in venoms from different localities.

In-vitro neurotoxicity of two Malaysian krait species (Bungarus candidus and Bungarus fasciatus) venoms: neutralization by monovalent and polyvalent antivenoms from Thailand

Bungarus candidus and Bungarus fasciatus are two species of krait found in Southeast Asia. Envenoming by these snakes is often characterized by neurotoxicity and, without treatment, causes considerable morbidity and mortality. In this study, the in vitro neurotoxicity of each species, and the effectiveness of two monovalent antivenoms and a polyvalent antivenom, against the neurotoxic effects of the venoms, were examined in a skeletal muscle preparation. Both venoms caused concentration-dependent inhibition of indirect twitches, and attenuated responses to exogenous nicotinic receptor agonists, in the chick biventer preparation, with B. candidus venom being more potent than B. fasciatus venom. SDS-PAGE and western blot analysis indicated different profiles between the venoms. Despite these differences, most proteins bands were recognized by all three antivenoms. Antivenom, added prior to the venoms, attenuated the neurotoxic effect of the venoms. Interestingly, the respective monovalent antivenoms did not neutralize the effects of the venom from the other Bungarus species indicating a relative absence of cross-neutralization. Addition of a high concentration of polyvalent antivenom, at the t₉₀ time point after addition of venom, partially reversed the neurotoxicity of B. fasciatus venom but not B. candidus venom. The monovalent antivenoms had no significant effect when added at the t₉₀ time point. This study showed that B. candidus and B. fasciatus venoms display marked in vitro neurotoxicity in the chick biventer preparation and administration of antivenoms at high dose is necessary to prevent or reverse neurotoxicity.

Exact solutions for species tree inference from discordant gene trees

Phylogenetic analysis has to overcome the grant challenge of inferring accurate species trees from evolutionary histories of gene families (gene trees) that are discordant with the species tree along whose branches they have evolved. Two well studied approaches to cope with this challenge are to solve either biologically informed gene tree parsimony (GTP) problems under gene duplication, gene loss, and deep coalescence, or the classic RF supertree problem that does not rely on any biological model. Despite the potential of these problems to infer credible species trees, they are NP-hard. Therefore, these problems are addressed by heuristics that typically lack any provable accuracy and precision. We describe fast dynamic programming algorithms that solve the GTP problems and the RF supertree problem exactly, and demonstrate that our algorithms can solve instances with data sets consisting of as many as 22 taxa. Extensions of our algorithms can also report the number of all optimal species trees, as well as the trees themselves. To better asses the quality of the resulting species trees that best fit the given gene trees, we also compute the worst case species trees, their numbers, and optimization score for each of the computational problems. Finally, we demonstrate the performance of our exact algorithms using empirical and simulated data sets, and analyze the quality of heuristic solutions for the studied problems by contrasting them with our exact solutions.

[New weak toxins from the cobra venom]

A protein with M 7485 Da containing five disulfide bonds was isolated from the venom of cobra Naja oxiana using various types of liquid chromatography. The complete amino acid sequence of the protein was determined by protein chemistry methods, which permitted us to assign it to the group of weak toxins. This is the first weak toxin isolated from the venom of N. oxiana. In a similar way, two new toxins with M 7628 and 7559 Da, which fall into the range of weak toxin masses, were isolated from the venom of the cobra N. kaouthia. The characterization of these proteins using Edman degradation and MALDI mass spectrometry has shown that one of these proteins is a novel weak toxin and the other is the known weak toxin WTX with an oxidized methionine residue in position 9. Such a modification was detected in weak toxins for the first time. A study of the biological activity of the toxin from N. oxiana showed that, like other weak toxins, it can be bound by muscle nicotinic cholinoreceptors and alpha7 nicotinic cholinoreceptors.