Nucleotide sequence of phospholipase A(2) gene expressed in snake pancreas reveals the molecular evolution of toxic phospholipase A(2) genes

Comparative venomics and preclinical efficacy evaluation of a monospecific Hemachatus antivenom towards sub-Saharan Africa cobra venoms

Cobras are the most medically important elapid snakes in Africa. The African genera Naja and Hemachatus include snakes with neurotoxic and cytotoxic venoms, with shared biochemical, toxinological and antigenic characteristics. We have studied the antigenic cross-reactivity of four sub-Saharan Africa cobra venoms against an experimental monospecific Hemachatus haemachatus antivenom through comparative proteomics, preclinical assessment of neutralization, and third generation antivenomics. The venoms of H. haemachatus, N. annulifera, N. mossambica and N. nigricollis share an overall qualitative family toxin composition but depart in their proportions of three-finger toxin (3FTxs) classes, phospholipases A₂ (PLA₂s), snake venom metalloproteinases (SVMPs), and cysteine-rich secretory proteins (CRISPs). A monospecific anti-Hemachatus antivenom produced by Costa Rican Instituto Clodomiro Picado neutralized the lethal activity of the homologous and heterologous neuro/cytotoxic (H. haemachatus) and cyto/cardiotoxic (N. mossambica and N. nigricollis) venoms of the three spitting cobras sampled, while it was ineffective against the lethal and toxic activities of the neurotoxic venom of the non-spitting snouted cobra N. annulifera. The ability of the anti-Hemachatus-ICP antivenom to neutralize toxic (dermonecrotic and anticoagulant) and enzymatic (PLA₂) activities of spitting cobra venoms suggested a closer kinship of H. haemachatus and Naja subgenus Afrocobra spitting cobras than to Naja subgenus Uraeus neurotoxic taxa. These results were confirmed by third generation antivenomics. BIOLOGICAL SIGNIFICANCE: African Naja species represent the most widespread medically important elapid snakes across Africa. To gain deeper insight into the spectrum of medically relevant toxins, we compared the proteome of three spitting cobras (Hemachatus haemachatus, Naja mossambica and N. nigricollis) and one non-spitting cobra (N. annulifera). Three finger toxins and phospholipases A₂ are the two major protein families among the venoms analyzed. The development of antivenoms of broad species coverage is an urgent need in sub-Saharan Africa. An equine antivenom raised against H. haemachatus venom showed cross-reactivity with the venoms of H. haemachatus, N. mossambica and N. nigricollis, while having poor recognition of the venom of N. annulifera. This immunological information provides clues for the design of optimum venom mixtures for the preparation of broad spectrum antivenoms.

How the Toxin got its Toxicity

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules—“toxins”—share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as “exochemicals.” There has been debate within the field of toxinology concerning the evolutionary pathways leading to the “recruitment” of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel “environment,” i.e., a range of novel partners for molecular interaction. After stressing the distinction between “activity” and “function,” we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of “neofunctionalization”—simply the origins of a novel molecular function by any process—and “recruitment”—the “weaponization” of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of “birth-and-death.” Thus, recruitment leads to a form of “downwards causation,” in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

Toxinology provides multidirectional and multidimensional opportunities: A personal perspective

In nature, toxins have evolved as weapons to capture and subdue the prey or to counter predators or competitors. When they are inadvertently injected into humans, they cause symptoms ranging from mild discomfort to debilitation and death. Toxinology is the science of studying venoms and toxins that are produced by a wide variety of organisms. In the past, the structure, function and mechanisms of most abundant and/or most toxic components were characterized to understand and to develop strategies to neutralize their toxicity. With recent technical advances, we are able to evaluate and determine the toxin profiles using transcriptomes of venom glands and proteomes of tiny amounts of venom. Enormous amounts of data from these studies have opened tremendous opportunities in many directions of basic and applied research. The lower costs for profiling venoms will further fuel the expansion of toxin database, which in turn will provide greater exciting and bright opportunities in toxin research.

Venom Systems as Models for Studying the Origin and Regulation of Evolutionary Novelties

A central goal in biology is to determine the ways in which evolution repeats itself. One of the most remarkable examples in nature of convergent evolutionary novelty is animal venom. Across diverse animal phyla, various specialized organs and anatomical structures have evolved from disparate developmental tissues to perform the same function, that is, produce and deliver a cocktail of potent molecules to subdue prey or predators. Venomous organisms therefore offer unique opportunities to investigate the evolutionary processes of convergence of key adaptive traits, and the molecular mechanisms underlying the emergence of novel genes, cells, and tissues. Indeed, some venomous species have already proven to be highly amenable as models for developmental studies, and recent work with venom gland organoids provides manipulatable systems for directly testing important evolutionary questions. Here, we provide a synthesis of the current knowledge that could serve as a starting point for the establishment of venom systems as new models for evolutionary and molecular biology. In particular, we highlight the potential of various venomous species for the study of cell differentiation and cell identity, and the regulatory dynamics of rapidly evolving, highly expressed, tissue-specific, gene paralogs. We hope that this review will encourage researchers to look beyond traditional study organisms and consider venom systems as useful tools to explore evolutionary novelties.

Clinical Condition and Management of 114 Mamushi (Gloydius blomhoffii) Bites in a General Hospital in Japan

Objective Mamushi (Gloydius blomhoffii) snakebite is the most common type of snake injury in Japan and is also seen in China and Korea. Although the components of Mamushi venom have been investigated, epidemiological and clinical descriptions still remain limited in the English literature. The aim of this study was to review the clinical features and management of patients with injuries related to Mamushi snakebites. Methods We conducted a retrospective examination of 114 Mamushi snakebite cases encountered at a general hospital in Japan from January 2004 to November 2016. Data were collected from the medical records. Results We found that Mamushi snakebites commonly occurred during summer and the daytime, with elderly men typically being affected. The symptom grade at initial consultation was significantly worse in the walk-in group than in the ambulance admission group, probably due to treatment delay. The number of fangs that pierced the skin was not related to the severity of the symptoms. The group treated with a tourniquet more frequently exhibited exacerbation of symptoms than those that received other treatments (p<0.001). Conclusion The delay between patients being bitten and arriving at hospital as well as the number of fangs that pierced the skin did not affect the duration of hospitalization; however, proximal tourniquation should be avoided in such cases, as significant exacerbation of local symptoms was observed when this procedure was applied.

Structural characteristics and evolution of the Protobothrops elegans pancreatic phospholipase A2 gene in contrast with those of Protobothrops genus venom phospholipase A2 genes

The nucleotide sequence of the gene encoding Protobothrops elegans (Crotalinae) pancreatic phospholipase A(2) (PLA(2)), abbreviated PePancPLA(2), was determined by means of inverted PCR techniques. Since its deduced amino acid sequence contains a pancreatic loop and shows high similarity to that of Laticauda semifasciata (Elapinae) group IB pancreatic PLA(2), PePancPLA(2) is classified into group IB PLA(2). The nucleotide sequences of the PePancPLA(2) gene, the L. semifasciata group IB pancreatic PLA(2) gene, and the L. semifasciata group IA venom PLA(2) gene are similar to one another but greatly dissimilar to those of Protobothrops genus (Crotalinae) group II venom PLA(2) genes, suggesting that the Elapinae group IB PLA(2) gene and the group IA PLA(2) gene appeared after Elapinae was established, and that the Crotalinae group II venom PLA(2) genes came into existence before Elapinae and Crotalinae diverged. A phylogenetic analysis of their amino acid sequences confirms this.

Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes

BACKGROUND

Kraits (genus Bungarus) and cobras (genus Naja) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for Bungarus multicinctus and Naja atra, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; i.e., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP).

RESULTS

In total, 1092 valid expressed sequences tags (ESTs) for B. multicinctus and 1166 ESTs for N. atra were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and β bungarotoxin (25.1%) comprise the main toxin classes in B. multicinctus, while 3FTx (95.8%) is the dominant toxin in N. atra. We also observed several less abundant venom families in B. multicinctus and N. atra, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in B. multicinctus. A total of 71 positive toxin BAC clones in B. multicinctus and N. atra were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-positive BACs were sequenced to reveal gene structures of 3FTx toxin genes.

CONCLUSIONS

Based on the toxin ESTs and 3FTx gene sequences, the major components of B. multicinctus venom transcriptome are neurotoxins, including long chain alpha neurotoxins (α-ntx) and the recently originated β bungarotoxin, whereas the N. atra venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain α-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (dN/dS) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.

Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom

BACKGROUND

Micrurus corallinus (coral snake) is a tropical forest snake belonging to the family Elapidae. Its venom shows a high neurotoxicity associated with pre- and post-synaptic toxins, causing diaphragm paralysis, which may result in death. In spite of a relatively small incidence of accidents, serum therapy is crucial for those bitten. However, the adequate production of antiserum is hampered by the difficulty in obtaining sufficient amounts of venom from a small snake with demanding breeding conditions. In order to elucidate the molecular basis of this venom and to uncover possible immunogens for an antiserum, we generated expressed sequences tags (ESTs) from its venom glands and analyzed the transcriptomic profile. In addition, their immunogenicity was tested using DNA immunization.

RESULTS

A total of 1438 ESTs were generated and grouped into 611 clusters. Toxin transcripts represented 46% of the total ESTs. The two main toxin classes consisted of three-finger toxins (3FTx) (24%) and phospholipases A(2) (PLA(2)s) (15%). However, 8 other classes of toxins were present, including C-type lectins, natriuretic peptide precursors and even high-molecular mass components such as metalloproteases and L-amino acid oxidases. Each class included an assortment of isoforms, some showing evidence of alternative splicing and domain deletions. Five antigenic candidates were selected (four 3FTx and one PLA(2)) and used for a preliminary study of DNA immunization. The immunological response showed that the sera from the immunized animals were able to recognize the recombinant antigens.

CONCLUSION

Besides an improvement in our knowledge of the composition of coral snake venoms, which are very poorly known when compared to Old World elapids, the expression profile suggests abundant and diversified components that may be used in future antiserum formulation. As recombinant production of venom antigens frequently fails due to complex disulfide arrangements, DNA immunization may be a viable alternative. In fact, the selected candidates provided an initial evidence of the feasibility of this approach, which is less costly and not dependent on the availability of the venom.

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.

Inventing an arsenal: adaptive evolution and neofunctionalization of snake venom phospholipase A2 genes

BACKGROUND

Gene duplication followed by functional divergence has long been hypothesized to be the main source of molecular novelty. Convincing examples of neofunctionalization, however, remain rare. Snake venom phospholipase A2 genes are members of large multigene families with many diverse functions, thus they are excellent models to study the emergence of novel functions after gene duplications.

RESULTS

Here, I show that positive Darwinian selection and neofunctionalization is common in snake venom phospholipase A2 genes. The pattern of gene duplication and positive selection indicates that adaptive molecular evolution occurs immediately after duplication events as novel functions emerge and continues as gene families diversify and are refined. Surprisingly, adaptive evolution of group-I phospholipases in elapids is also associated with speciation events, suggesting adaptation of the phospholipase arsenal to novel prey species after niche shifts. Mapping the location of sites under positive selection onto the crystal structure of phospholipase A2 identified regions evolving under diversifying selection are located on the molecular surface and are likely protein-protein interactions sites essential for toxin functions.

CONCLUSION

These data show that increases in genomic complexity (through gene duplications) can lead to phenotypic complexity (venom composition) and that positive Darwinian selection is a common evolutionary force in snake venoms. Finally, regions identified under selection on the surface of phospholipase A2 enzymes are potential candidate sites for structure based antivenin design.

Molecular evolution of toxin genes in Elapidae snakes

The venom of the sea krait, Laticauda semifasciata, consists primarily of two toxic proteins, phospholipase A(2) (PLA(2)) and a three-finger-structure toxin. We have cloned both toxic protein genes, including the upstream region. PLA(2) genes contain three types of inserted sequences: an AG-rich region, a chicken repeat 1-like long interspersed nucleotide element sequence and an intron II 3' side repeat sequence. The molecular divergence of L. semifasciata PLA(2) genes was defined on the basis of the inserted sequences and their sequence homology. The length of intron I in the three-finger-structure toxin genes differs from species to species. The alignment analysis of intron I of the three-finger-structure toxin genes revealed that the intron I sequence of the ancestral gene comprised ten genetic regions. A hypothetical evolutionary process for the three-finger-structure toxin genes has also been developed.

Snake fangs from the Lower Miocene of Germany: evolutionary stability of perfect weapons

There is a general consensus that most of today's nonvenomous snakes are descendants of venomous snakes that lost their venomous capabilities secondarily. This implies that the evolutionary history of venomous snakes and their venom apparatus should be older than the current evidence from the fossil record. We compared some of the oldest-known fossil snake fangs from the Lower Miocene of Germany with those of modern viperids and elapids and found their morphology to be indistinguishable from the modern forms. The primary function of recent elapid and viperid snake fangs is to facilitate the extremely rapid, stab-like application of highly toxic venoms. Our findings therefore indicate that the other components of the venom-delivery system of Early Miocene vipers and elapids were also highly developed, and that these snakes used their venom in the same way as their modern relatives. Thus, the fossil record supports the view that snakes used their venoms to rapidly subdue prey long before the mid-Tertiary onset of the global environmental changes that seem to have supported the successful radiation of venomous snakes.

Identification of crotasin, a crotamine-related gene of Crotalus durissus terrificus

Crotamine is a cationic peptide (4.9 kDa, pI 9.5) of South American rattlesnake, Crotalus durissus terrificus' venom. Its presence varies according to the subspecies or the geographical locality of a given species. At the genomic level, we observed the presence of 1.8 kb gene, Crt-p1, in crotamine-positive specimens and its absence in crotamine-negative ones. In this work, we described a crotamine-related 2.5 kb gene, crotasin (Cts-p2), isolated from crotamine-negative specimens. Reverse transcription coupled to polymerase chain reaction indicates that Cts-p2 is abundantly expressed in several snake tissues, but scarcely expressed in the venom gland. The genome of crotamine-positive specimen contains both Crt-p1 and Cts-p2 genes. The present data suggest that both crotamine and crotasin have evolved by duplication of a common ancestor gene, and the conservation of their three disulfide bonds indicates that they might adopt the same fold as beta-defensin. The physiological function of the crotasin is not yet known.

Biochemical characterization of phospholipase A2 (trimorphin) from the venom of the Sonoran Lyre Snake Trimorphodon biscutatus lambda (family Colubridae)

Phospholipases A(2) (PLA(2)), common venom components and bioregulatory enzymes, have been isolated and sequenced from many snake venoms, but never from the venom (Duvernoy's gland secretion) of colubrid snakes. We report for the first time the purification, biochemical characterization and partial sequence of a PLA(2) (trimorphin) from the venom of a colubrid snake, Trimorphodon biscutatus lambda (Sonoran Lyre Snake). Specific phospholipase activity of the purified PLA(2) was confirmed by enzyme assays. The molecular weight of the enzyme has been determined by SDS-PAGE and mass spectrometry to be 13,996 kDa. The sequence of 50 amino acid residues from the N-terminal has been identified and shows a high degree of sequence homology to the type IA PLA(2)s, especially the Asp-49 enzymes. The Cys-11 residue, characteristic of the group IA PLA(2)s, and the Ca(2+) binding loop residues (Tyr-28, Gly-30, Gly-32, and Asp-49) are conserved. In addition, the His-48 residue, a key component of the active site, is also conserved in trimorphin. The results of phylogenetic analysis on the basis of amino acid sequence homology demonstrate that trimorphin belongs to the type IA family, and it appears to share a close evolutionary relationship with the PLA(2)s from hydrophiine elapid snakes (sea snakes and Australian venomous snakes).

Comparative analysis of gene expression mechanisms between group IA and IB phospholipase A2 genes from sea snake Laticauda semifasciata

Phospholipase A(2) (PLA(2)) genes expressed in the venom glands of the sea snake, Laticauda semifasciata, were investigated. Both mRNAs, encoding group IA (without a pancreatic loop) and group IB (with pancreatic loop), were detected from venom glands by Northern blot hybridization analysis and RT-PCR. The results of quantitative PCR analysis indicated that the expression amount of group IA genes was around 100-300 times greater than that of group IB genes. Sequence analysis of 5'-upstream regions and a reporter gene assay of the genes (groups IA and IB) previously cloned showed that the functional sequence (411 bp) was inserted in the 5'-flanking region of the group IA PLA(2) genes. It seemed that the contribution of the inserted sequence to the amount of transcribed mRNAs was greater than that of number of genes present in the genome. Comparative analysis of the 5'-flanking sequences from several snake genes encoding toxic PLA(2)s revealed that this sequence was probably inserted into an ancestral gene of PLA(2) with a pancreatic loop. After the duplication of the gene, which contained the inserted sequence, the PLA(2) gene without a pancreatic loop evolved from one of the duplicate genes. This inserted sequence might determine the future of the genes expressed in the venom glands.

Identification of alpha-bungarotoxin (A31) as the major postsynaptic neurotoxin, and complete nucleotide identity of a genomic DNA of Bungarus candidus from Java with exons of the Bungarus multicinctus alpha-bungarotoxin (A31) gene

The Malayan krait (Bungarus candidus) is one of the most medically significant snake species in Southeast Asia. No specific antivenom exists to treat envenoming by this species. Death within 30 min after its bite has been reported from Java, suggesting the presence of highly lethal postsynaptic neurotoxins in the venom of these snakes. We purified and identified the major postsynaptic toxin in the venom of B. candidus from Java. The toxin was indistinguishable from alpha-bungarotoxin (A31), a toxin originally isolated from Bungarus multicinctus, in its mass (7983.75 Da), LD50 (0.23 microg/g in mice i.p.), affinity to nicotinic acetylcholine receptors, and by its 40 N-terminal amino acid residues as determined by Edman degradation. Identity with alpha-bungarotoxin was confirmed by cloning and sequencing a genomic DNA from B. candidus which encodes the 74 amino acid sequence of alpha-bungarotoxin (A31) and part of its signal peptide, revealing complete identity to the alpha-bungarotoxin (A31) gene in exon and 98.9% identity in intron sequences. The entire mitochondrial cytochrome b gene of the krait species B. candidus from Java and B. multicinctus from Taiwan was sequenced for comparison, suggesting that these snakes are phylogenetically closely related. alpha-Bungarotoxin appears to be widely present and conserved in Southeast and East Asian black-and-white kraits across populations and taxa.