Molecular cloning of the major lethal toxins from two kraits (Bungarus flaviceps and Bungarus candidus)

On characterizing the Red-headed Krait (Bungarus flaviceps) venom: Decomplexation proteomics, immunoreactivity and toxicity cross-neutralization by hetero-specific antivenoms

The Red-headed Krait (Bungarus flaviceps) is a medically important venomous snake species in Southeast Asia, while there is no specific antivenom available for its envenoming. This study investigated the venom composition through a decomplexation proteomic approach, and examined the immunoreactivity as well as cross-neutralization efficacy of two hetero-specific krait antivenoms, Bungarus candidus Monovalent Antivenom (BcMAV) and Bungarus fasciatus Monovalent Antivenom (BfMAV), against the venom of B. flaviceps from Peninsular Malaysia. A total of 43 non-redundant proteoforms belonging to 10 toxin families were identified in the venom proteome, which is dominated by phospholipases A₂ including beta-bungarotoxin lethal subunit (56.20 % of total venom proteins), Kunitz-type serine protease inhibitors (19.40 %), metalloproteinases (12.85 %) and three-finger toxins (7.73 %). The proteome varied in quantitative aspect from the earlier reported Indonesian (Sumatran) sample, suggesting geographical venom variation. BcMAV and BfMAV were immunoreactive toward the B. flaviceps venom, with BcMAV being more efficacious in immunological binding. Both antivenoms cross-neutralized the venom lethality with varying efficacy, where BcMAV was more potent than BfMAV by ~13 times (normalized potency: 38.04 mg/g vs. 2.73 mg/g, defined as the venom amount completely neutralized by one-gram antivenom protein), supporting the potential utility of BcMAV for para-specific neutralization against B. flaviceps venom.

Complex target SELEX-based identification of DNA aptamers against Bungarus caeruleus venom for the detection of envenomation using a paper-based device

Complex target SELEX always have been an intriguing approach to the scientific community, as it offers the potential discovery of novel biomarkers. We herein successfully performed SELEX on Bungarus caeruleus venom to develop a panel of highly affine aptamers that specifically recognizes the B. caeruleus (common krait) venom and was able to discriminate the B. caeruleus venom from Cobra, Russell's, and Saw-scaled viper's venom. The aptamers generated against the crude venom also lead to the identification of the specific component of the venom, which is β-Bungarotoxin, a toxin uniquely present in the B. caeruleus venom. The best performing aptamer candidates were used as a molecular recognition element in a paper-based device and were able to detect as low as 2 ng krait venom in human serum background. The developed aptamer-based paper device can be used for potential point-of-care venom detection applications due to its simplicity and affordability.

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.

Kinetic analysis of effects of temperature and time on the regulation of venom expression in Bungarus multicinctus

Venom gland is a highly efficient venom production system to maintain their predatory arsenal. Venom toxins mRNA has been shown to increase abruptly in snake after venom expenditure, while the dynamics of venom accumulation during synthesis are poorly understood. Here, PacBio long-read sequencing, Illumina RNA sequencing (RNA-seq), and label-free proteome quantification were used to investigate the composition landscape and time- and temperature-dependent dynamics changes of the Bungarus multicinctus venom gland system. Transcriptome data (19.5223 Gb) from six adult B. multicinctus tissues were sequenced using PacBio RS II to generate a reference assembly, and average 7.28 Gb of clean RNA-seq data was obtained from venom glands by Illumina sequencing. Differentially expressed genes (DEGs) mainly were protein processing rather than venom toxins. Post-translational modifications provided the evidence of the significantly different proportions of toxins in the venom proteome with the changing of replenishment time and temperature, but constant of venom toxins mRNA in the venom gland transcriptome. Dynamic of toxins and genes involved in venom gland contraction suggesting the formation of the mature venom gland system would take at least 9 days. In addition, 59 toxin processing genes were identified, peptidylprolyl isomerase B of which underwent positive selection in Toxicofera. These results provide a reference for determining the extraction time of venom, production of polyclonal and monoclonal antibody for precise treatment plans of venomous snakebites, and construction of an in vitro synthesis system for snake venom protein.

Comparison of proteomic profiles of the venoms of two of the 'Big Four' snakes of India, the Indian cobra (Naja naja) and the common krait (Bungarus caeruleus), and analyses of their toxins

Snake venoms are mixtures of biologically-active proteins and peptides, and several studies have described the characteristics of some of these toxins. However, complete proteomic profiling of the venoms of many snake species has not yet been done. The Indian cobra (Naja naja) and common krait (Bungarus caeruleus) are elapid snake species that are among the 'Big Four' responsible for the majority of human snake envenomation cases in India. As understanding the composition and complexity of venoms is necessary for successful treatment of envenomation in humans, we utilized three different proteomic profiling approaches to characterize these venoms: i) one-dimensional SDS-PAGE coupled with in-gel tryptic digestion and electrospray tandem mass spectrometry (ESI-LC-MS/MS) of individual protein bands; ii) in-solution tryptic digestion of crude venoms coupled with ESI-LC-MS/MS; and iii) separation by gel-filtration chromatography coupled with tryptic digestion and ESI-LC-MS/MS of separated fractions. From the generated data, 81 and 46 different proteins were identified from N. naja and B. caeruleus venoms, respectively, belonging to fifteen different protein families. Venoms from both species were found to contain a variety of phospholipases A₂ and three-finger toxins, whereas relatively higher numbers of snake venom metalloproteinases were found in N. naja compared to B. caeruleus venom. The analyses also identified less represented venom proteins including L-amino acid oxidases, cysteine-rich secretory proteins, 5'-nucleotidases and venom nerve growth factors. Further, Kunitz-type serine protease inhibitors, cobra venom factors, phosphodiesterases, vespryns and aminopeptidases were identified in the N. naja venom, while acetylcholinesterases and hyaluronidases were found in the B. caeruleus venom. We further analyzed protein coverage (Lys/Arg rich and poor regions as well as potential glycosylation sites) using in-house software. These studies expand our understanding of the proteomes of the venoms of these two medically-important species.

Codon usage pattern and prediction of gene expression level in Bungarus species

Codon bias study in an organism gains significance in understanding the molecular mechanism as well as the functional conservation of gene expression during the course of evolution. The prime focus in this study is to compare the codon usage patterns among the four species belonging to the genus Bungarus (B. multicinctus, B. fasciatus, B. candidus and B. flaviceps) using several codon bias parameters. Our results suggested that relatively low codon bias exists in the coding sequences of the selected species. The compositional constraints together with gene expression level might influence the patterns of codon usage among the genes of Bungarus species. Both natural selection and mutation pressure affect the codon usage pattern in Bungarus species as evident from correspondence analysis. Neutrality plot indicates that natural selection played a major role while mutation pressure played a minor role in codon usage pattern of the genes in Bungarus species.

Proteomic characterization and comparison of venoms from two elapid snakes (Bungarus multicinctus and Naja atra) from China

Bungarus multicinctus (many-banded krait) and Naja atra (Chinese cobra) are widely distributed and medically important venomous snakes in China; however, their venom proteomic profiles have not been fully compared. Here, we fractionated crude venoms and analyzed them using a combination of proteomic techniques. Three-finger toxins (3-FTx) and phospholipase A2 (PLA2) were most abundant in both species, respectively accounting for 32.6% and 66.4% of total B. multicinctus venom, and 84.3% and 12.2% of total N. atra venom. Venoms from these two species contained one common protein family and six less abundant species-specific protein families. The proteomic profiles of B. multicinctus and N. atra venoms and analysis of toxicological activity in mice suggested that 3-FTx and PLA2 are the major contributors to clinical symptoms caused by envenomation. The venoms differed in enzymatic activity, likely the result of inter-specific variation in the amount of related venom components. Antivenomics assessment revealed that a small number of venom components (3-FTxs and PLA2s in B. multicinctus, and 3-FTxs in N. atra) could not be immunocaptured completely, suggesting that we should pay attention to enhancing the immune response of these components in designing commercial antivenoms for B. multicinctus and N. atra.

BIOLOGICAL SIGNIFICANCE

The proteomic profiles of venoms from two medically important snake species - B. multicinctus and N. atra - have been explored. Quantitative and qualitative differences are evident in both venoms when proteomic profiles and transcriptomic results are compared; this is a reminder that combined approaches are needed to explore the precise composition of snake venom. Two protein families (3-FTx and PLA2) of high abundance in these snake venoms are major players in the biochemical and pharmacological effects of envenomation. Elucidation of the proteomic profiles of these snake venoms is helpful in understanding composition-function relationships and will facilitate the clinical application of antivenoms.

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.

Transcriptomic analysis of the venom gland of the red-headed krait (Bungarus flaviceps) using expressed sequence tags

Background

The Red-headed krait (Bungarus flaviceps, Squamata: Serpentes: Elapidae) is a medically important venomous snake that inhabits South-East Asia. Although the venoms of most species of the snake genus Bungarus have been well characterized, a detailed compositional analysis of B. flaviceps is currently lacking.

Results

Here, we have sequenced 845 expressed sequence tags (ESTs) from the venom gland of a B. flaviceps. Of the transcripts, 74.8% were putative toxins; 20.6% were cellular; and 4.6% were unknown. The main venom protein families identified were three-finger toxins (3FTxs), Kunitz-type serine protease inhibitors (including chain B of β-bungarotoxin), phospholipase A₂ (including chain A of β-bungarotoxin), natriuretic peptide (NP), CRISPs, and C-type lectin.

Conclusion

The 3FTxs were found to be the major component of the venom (39%). We found eight groups of unique 3FTxs and most of them were different from the well-characterized 3FTxs. We found three groups of Kunitz-type serine protease inhibitors (SPIs); one group was comparable to the classical SPIs and the other two groups to chain B of β-bungarotoxins (with or without the extra cysteine) based on sequence identity. The latter group may be functional equivalents of dendrotoxins in Bungarus venoms. The natriuretic peptide (NP) found is the first NP for any Asian elapid, and distantly related to Australian elapid NPs. Our study identifies several unique toxins in B. flaviceps venom, which may help in understanding the evolution of venom toxins and the pathophysiological symptoms induced after envenomation.

The evolution of two partner LINE/SINE families and a full-length chromodomain-containing Ty3/Gypsy LTR element in the first reptilian genome of Anolis carolinensis

Transposable elements have been characterized in a number of vertebrates, including whole genomes of mammals, birds, and fishes. The Anolis carolinensis draft assembly provides the first opportunity to study retroposons in a reptilian genome. Here, we identified and reconstructed a number of retroposons based on database searches: Five Sauria short interspersed element (SINE) subfamilies, 5S-Sauria SINE chimeras, Anolis Bov-B long interspersed element (LINE), Anolis SINE 2, Anolis LINE 2, Anolis LINE 1, Anolis CR 1, and a chromodomain-containing Ty3/Gypsy LTR element. We focused on two SINE families (Anolis Sauria SINE and Anolis SINE 2) and their partner LINE families (Anolis Bov-B LINE and Anolis LINE 2). We demonstrate that each SINE/LINE pair is distributed similarly and predict that the retrotransposition of evolutionarily younger Sauria SINE members is via younger Bov-B LINE members while a correlation also exists between their respective evolutionarily older SINE/LINE members. The evolutionarily youngest Sauria SINE sequences evolved as part of novel rolling-circle transposons. The evolutionary time frame when Bov-B LINEs and Sauria SINEs were less active in their retrotransposition is characterized by a high retrotransposition burst of Anolis SINE 2 and Anolis LINE 2 elements. We also characterized the first full-length chromoviral LTR element in amniotes (Amn-ichi). This newly identified chromovirus is widespread in the Anolis genome and has been very well preserved, indicating that it is still active. Transposable elements in the Anolis genome account for approximately 20% of the total DNA sequence, whereas the proportion is more than double that in many mammalian genomes in which such elements have important biological functions. Nevertheless, 20% transposable element coverage is sufficient to predict that Anolis retroposons and other mobile elements also may have biologically and evolutionarily relevant functions. The new SINEs and LINEs and other ubiquitous genomic elements characterized in the Anolis genome will prove very useful for studies in comparative genomics, phylogenetics, and functional genetics.

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

Eight phospholipases A2 (PLAs) and four three-finger-toxins (3FTx) from the pooled venom of Bungarus fasciatus (Bf) were previously studied and sequenced, but their expression pattern in individual Bf venom and possible geographic variations remained to be investigated. We herein analyze the individual venom of two Bf specimens from Kolkata (designated as KBf) to address this question. Seven PLAs and five 3FTx were purified from the KBf venoms, and respective cDNAs were cloned from venom glands of one of the snakes. Comparison of their mass and N-terminal sequence revealed that all the PLAs were conserved in both KBf venoms, but that two of their 3FTx isoforms were variable. When comparing the sequences of these KBf-PLAs with those published, only one was found to be identical to that of Bf Vb-2, and the other five were 94-98% identical to those of Bf II, III, Va, VI and XI-2, respectively. Notably, the most abundant PLA isoforms of Bf and KBf venoms contain Pro31 substitution. They were found to have abnormally low k(cat) values but high affinity for Ca2+. Phylogenetic analysis based on the sequences of venom group IA PLAs showed a close relationship between Bungarus and Australian and marine Elapidae. As the five deduced sequences of KBf-3FTx are only 62-82% identical to the corresponding Bf-3FTx from the pooled venom, the 3FTx apparently have higher degree of individual and geographic variations than the PLAs. None of the KBf-3FTx was found to be neurotoxic or very lethal; phylogenetic analyses of the 3FTx also revealed the unique evolution of Bf as compared with other kraits.