Unveiling Novel Kunitz- and Waprin-Type Toxins in the Micrurus mipartitus Coral Snake Venom Gland: An In Silico Transcriptome Analysis.

Kunitz-type peptide expression has been described in the venom of snakes of the Viperidae, Elapidae and Colubridae families. This work aimed to identify these peptides in the venom gland transcriptome of the coral snake Micrurus mipartitus. Transcriptomic analysis revealed a high diversity of venom-associated Kunitz serine protease inhibitor proteins (KSPIs). A total of eight copies of KSPIs were predicted and grouped into four distinctive types, including short KSPI, long KSPI, Kunitz-Waprin (Ku-WAP) proteins, and a multi-domain Kunitz-type protein. From these, one short KSPI showed high identity with Micrurus tener and Austrelaps superbus. The long KSPI group exhibited similarity within the Micrurus genus and showed homology with various elapid snakes and even with the colubrid Pantherophis guttatus. A third group suggested the presence of Kunitz domains in addition to a whey-acidic-protein-type four-disulfide core domain. Finally, the fourth group corresponded to a transcript copy with a putative 511 amino acid protein, formerly annotated as KSPI, which UniProt classified as SPINT1. In conclusion, this study showed the diversity of Kunitz-type proteins expressed in the venom gland transcriptome of M. mipartitus.

Protein composition and biochemical characterization of venom from Sonoran Coral Snakes (Micruroides euryxanthus).

The elapid genus, Micruroides, is considered the sister clade of all New World coral snakes (Genus Micrurus), is monotypic, and is represented by Sonoran Coral Snakes, Micruroides euryxanthus. Coral snakes of the genus Micrurus have been reported to have venoms that are predominantly composed of phospholipases A2 (PLA2) or three finger toxins (3FTx), but the venoms of the genus Micruroides are almost completely unstudied. Here, we present the first description of the venom of M. euryxanthus including identification of some proteins as well as transcriptomic, and biological activity assays. The most abundant components within M. euryxanthus venom are 3FTxs (62.3%) and there was relatively low proportion of PLA2s (14.2%). The venom phenotype supports the hypothesis that the common ancestor of Micrurus and Micruroides had a 3FTx-dominated venom. Within the venom, there were two nearly identical α-neurotoxins (α-Ntx), one of which was designated Eurytoxin, that account for approximately 60% of the venom's lethality to mice. Eurytoxin was cloned, expressed in a soluble and active form, and used to produce rabbit hyperimmune serum. This allowed the analysis of its immunochemical properties, showing them to be different from the recombinant αNTx D.H., present in the venoms of some species of Micrurus. Finally, we observed that the commercial antivenom produced in Mexico for coral snake envenomation is unable to neutralize the lethality from M. euryxanthus venom. This work allowed the classification of Micruroides venom into the 3FTx-predominant group and identified the main components responsible for toxicity to mice.

Functional, proteomic and transcriptomic characterization of the venom from Micrurus browni browni: Identification of the first lethal multimeric neurotoxin in coral snake venom.

Proteomic characterization of Micrurus browni browni venom showed approximately 41 components belonging to 9 protein families, mainly phospholipases A2 (PLA2s) and three-finger toxins (3FTxs). Venom gland transcriptome yielded 39 venom transcripts belonging to 10 protein families. Functional characterization identified a multimeric toxin, here designated Brownitoxin-1, which comprises at least one PLA2 and one 3FTx. Its components have no or very low lethality individually but become extremely lethal when combined; both were partially characterized. Other two lethal components were identified: A neurotoxic PLA2, and a postsynaptic α-neurotoxin. LD50s as well as PLA2 and nAChR-blocking activities were determined for whole venom and isolated components. Application of venom to murine neuromuscular preparations caused a progressive decrease of twitch force that was irreversible after washing. Inhibition of PLA2 activity with p-bromophenacyl bromide (pBPB) showed that approximately 90% of toxicity is dependent on this activity. Non-lethal components include diverse 3FTxs, at least three enzymatically active PLA2s and the nociceptive toxin MitTx. No evidence of specificity towards prey was observed. This work is one of the most complete characterizations of a coral snake venom so far and its findings highlight the relevance of protein complexes in venom function. SIGNIFICANCE: This study represents a profound analysis of the venom of the coral snake Micrurus browni browni, including a venom proteome, venom gland transcriptomic data and functional studies of whole venom and isolated toxins. It significantly contributes to the understanding of North American coral snake venoms, which are currently largely unknown. It includes characterization of relevant venom components, one of which represents the first description of a lethal multimeric neurotoxin in coral snake venom. This work highlights the importance of protein complexes in coral snake venom and could serve as a basis for the finding of several other multimeric toxins. Finally, we report the absence of taxon specificity, which has been previously reported in the venoms of other snakes of the same genus.

Genome-wide SNPs clarify lineage diversity confused by coloration in coralsnakes of the Micrurus diastema species complex (Serpentes: Elapidae).

New world coralsnakes of the genus Micrurus are a diverse radiation of highly venomous and brightly colored snakes that range from North Carolina to Argentina. Species in this group have played central roles in developing and testing hypotheses about the evolution of mimicry and aposematism. Despite their diversity and prominence as model systems, surprisingly little is known about species boundaries and phylogenetic relationships within Micrurus, which has substantially hindered meaningful analyses of their evolutionary history. Here we use mitochondrial genes together with thousands of nuclear genomic loci obtained via ddRADseq to study the phylogenetic relationships and population genomics of a subclade of the genus Micrurus: The M. diastema species complex. Our results indicate that prior species and species-group inferences based on morphology and color pattern have grossly misguided taxonomy, and that the M. diastema complex is not monophyletic. Based on our analyses of molecular data, we infer the phylogenetic relationships among species and populations, and provide a revised taxonomy for the group. Two non-sister species-complexes with similar color patterns are recognized, the M. distans and the M. diastema complexes, the first being basal to the monadal Micrurus and the second encompassing most North American monadal taxa. We examined all 13 species, and their respective subspecies, for a total of 24 recognized taxa in the M. diastema species complex. Our analyses suggest a reduction to 10 species, with no subspecific designations warranted, to be a more likely estimate of species diversity, namely, M. apiatus, M. browni, M. diastema, M. distans, M. ephippifer, M. fulvius, M. michoacanensis, M. oliveri, M. tener, and one undescribed species.

Unlocking the secrets of banded coral snake (Calliophis intestinalis, Malaysia): A venom with proteome novelty, low toxicity and distinct antigenicity.

The Asiatic coral snakes are basal in the phylogeny of coral snakes. Although envenoming by the Asiatic coral snakes is rarely fatal, little is known about their venom properties and variability from the American coral snakes. Integrating reverse-phase high performance liquid chromatography and nano-liquid chromatography-tandem mass spectrometry, we showed that the venom proteome of the Malaysian banded or striped coral snake (Calliophis intestinalis) was composed of mainly phospholipases A2 (PLA2, 43.4%) and three-finger toxins (3FTx, 20.1%). Within 3FTx, the cytotoxins or cardiotoxins (CTX) dominated while the neurotoxins' content was much lower. Its subproteomic details contrasted with the 3FTx profile of most Micrurus sp., illustrating a unique dichotomy of venom phenotype between the Old and the New World coral snakes. Calliophis intestinalis venom proteome was correlated with measured enzymatic activities, and in vivo it was myotoxic but non-lethal to mice, frogs and geckos at high doses (5-10 µg/g). The venom contains species-specific toxins with distinct sequences and antigenicity, and the antibodies raised against PLA2 and CTX of other elapids showed poor binding toward its venom antigens. The unique venom proteome of C. intestinalis unveiled a repertoire of novel toxins, and the toxicity test supported the need for post-bite monitoring of myotoxic complication. SIGNIFICANCE: Malaysian banded or striped coral snake (Calliophis intestinalis) has a cytotoxin (CTX)-predominating venom proteome, a characteristic shared by its congener, the Malayan blue coral snake (Calliophis bivirgata). With little neurotoxins (NTX), it illustrates a CTX/NTX dichotomy of venom phenotype between the Old World and the New World coral snakes. The low toxicity of the venom imply that C. intestinalis bite envenoming can be managed via symptomatic relief of the mild to moderate pain with appropriate analgesia. Systemically, the serum creatine kinase level of patients should be monitored serially for potential complication of myotoxicity. The distinct antigenicity of the venom proteins implies that the empirical use of heterologous antivenom is mostly inappropriate and not recommended.

Cloning and sequencing of three-finger toxins from the venom glands of four Micrurus species from Mexico and heterologous expression of an alpha-neurotoxin from Micrurus diastema.

The three-finger toxins (3FTxs) represent an extremely diverse protein family in elapid venoms, where the short chain α-neurotoxins are the most relevant toxin group from the clinical point of view. Essentially, the 3FTxs variability and the low proportions of α-neurotoxins in the venoms of North American coral snakes make it difficult to obtain effective elapid antivenoms against the envenomation symptoms caused mainly by these α-neurotoxins. In this work, thirty 3FTx transcript sequences were obtained from the venom glands of four coral snake species from Mexico (M. diastema, M. laticollaris, M. browni and M. tener). The transcripts were mined using a forward oligonucleotide based on the highly conserved signal peptide from the 3FTxs, and four of these transcripts, named MlatA1, B.D, B.E and D.H, encoded for short-chain α-neurotoxins. Additionally, one isoform of the D.H α-neurotoxin transcript was identified in the venom of M. diastema. The mature α-neurotoxin coded in the D.H transcript was heterologously expressed, and it was found soluble (4.2 mg/l) in the cytoplasm of a bacterial system. The recombinant D.H (rD.H) had an IC50 value of 31.5 ±â€¯4.4 nM on nicotinic acetylcholine receptors of the muscular type expressed in rhabdomyosarcoma cells (TE671). The rDH also had an LD50 of 0.15 mg/kg mice, and it was evaluated as a potential immunogen in New Zealand rabbits. The protective capacity of rabbit sera was tested against two native coral snake α-neurotoxins, and against rD.H. One of the anti-rD.H rabbit sera was able to neutralize the lethality of all three neurotoxins when tested in groups of CD1 mice. This work contributes to the increasing understanding of the high diversity of 3FTxs, and shows that recombinant coral snake α-neurotoxins are promising supplements for hyperimmunization protocols for coral snake antivenom production.

Ancient Diversification of Three-Finger Toxins in Micrurus Coral Snakes.

Coral snakes, most notably the genus Micrurus, are the only terrestrial elapid snakes in the Americas. Elapid venoms are generally known for their potent neurotoxicity which is usually caused by Three-Finger Toxin (3FTx) proteins. These toxins can have a wide array of functions that have been characterized from the venom of other elapids. We examined publicly available sequences from Micrurus 3FTx to show that they belong to 8 monophyletic clades that diverged as deep in the 3FTx phylogenetic tree as the other clades with characterized functions. Functional residues from previously characterized clades of 3FTx are not well conserved in most of the Micrurus toxin clades. We also analyzed the patterns of selection on these toxins and find that they have been diversifying at different rates, with some having undergone extreme diversifying selection. This suggests that Micrurus 3FTx may contain a previously underappreciated functional diversity that has implications for the clinical outcomes of bite victims, the evolution and ecology of the genus, as well as the potential for biodiscovery efforts focusing on these toxins.

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa.

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A2 (PLA2s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA2s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the ß-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A2s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.

Molecular cloning and structural modelling of gamma-phospholipase A2 inhibitors from Bothrops atrox and Micrurus lemniscatus snakes.

Phospholipases A2 inhibitors (PLIs) produced by venomous and non-venomous snakes play essential role in this resistance. These endogenous inhibitors may be classified by their fold in PLIα, PLIß and PLIγ. Phospholipases A2 (PLA2s) develop myonecrosis in snake envenomation, a consequence that is not efficiently neutralized by antivenom treatment. This work aimed to identify and characterize two PLIs from Amazonian snake species, Bothrops atrox and Micrurus lemniscatus. Liver tissues RNA of specimens from each species were isolated and amplified by RT-PCR using PCR primers based on known PLIγ gene sequences, followed by cloning and sequencing of amplified fragments. Sequence similarity studies showed elevated identity with inhibitor PLIγ gene sequences from other snake species. Molecular models of translated inhibitors' gene sequences resemble canonical three finger fold from PLIγ and support the hypothesis that the decapeptide (residues 107-116) may be responsible for PLA2 inhibition. Structural studies and action mechanism of these PLIs may provide necessary information to evaluate their potential as antivenom or as complement of the current ophidian accident treatment.

Coral snakes predict the evolution of mimicry across New World snakes.

Batesian mimicry, in which harmless species (mimics) deter predators by deceitfully imitating the warning signals of noxious species (models), generates striking cases of phenotypic convergence that are classic examples of evolution by natural selection. However, mimicry of venomous coral snakes has remained controversial because of unresolved conflict between the predictions of mimicry theory and empirical patterns in the distribution and abundance of snakes. Here we integrate distributional, phenotypic and phylogenetic data across all New World snake species to demonstrate that shifts to mimetic coloration in nonvenomous snakes are highly correlated with coral snakes in both space and time, providing overwhelming support for Batesian mimicry. We also find that bidirectional transitions between mimetic and cryptic coloration are unexpectedly frequent over both long- and short-time scales, challenging traditional views of mimicry as a stable evolutionary 'end point' and suggesting that insect and snake mimicry may have different evolutionary dynamics.