Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest.

Bothrops jararaca is a slender and semi-arboreal medically relevant pit viper species endemic to tropical and subtropical forests in southern Brazil, Paraguay, and northern Argentina (Misiones). Within its geographic range, it is often abundant and is an important cause of snakebite. Although no subspecies are currently recognized, geographic analyses have revealed the existence of two well-supported B. jararaca clades that diverged during the Pliocene ~3.8Mya and currently display a southeastern (SE) and a southern (S) Atlantic rainforest (Mata Atlântica) distribution. The spectrum, geographic variability, and ontogenetic changes of the venom proteomes of snakes from these two B. jararaca phylogroups were investigated applying a combined venom gland transcriptomic and venomic analysis. Comparisons of the venom proteomes and transcriptomes of B. jararaca from the SE and S geographic regions revealed notable interpopulational variability that may be due to the different levels of population-specific transcriptional regulation, including, in the case of the southern population, a marked ontogenetic venom compositional change involving the upregulation of the myotoxic PLA2 homolog, bothropstoxin-I. This population-specific marker can be used to estimate the proportion of venom from the southern population present in the B. jararaca venom pool used for the Brazilian soro antibotrópico (SAB) antivenom production. On the other hand, the southeastern population-specific D49-PLA2 molecules, BinTX-I and BinTX-II, lend support to the notion that the mainland ancestor of Bothrops insularis was originated within the same population that gave rise to the current SE B. jararaca phylogroup, and that this insular species endemic to Queimada Grande Island (Brazil) expresses a pedomorphic venom phenotype. Mirroring their compositional divergence, the two geographic B. jararaca venom pools showed distinct bioactivity profiles. However, the SAB antivenom manufactured in Vital Brazil Institute neutralized the lethal effect of both venoms to a similar extent. In addition, immobilized SAB antivenom immunocaptured most of the venom components of the venoms of both B. jararaca populations, but did not show immunoreactivity against vasoactive peptides. The Costa Rican bothropic-crotalic-lachesic (BCL) antivenom showed the same lack of reactivity against vasoactive peptides but, in addition, was less efficient immunocapturing PI- and PIII-SVMPs from the SE venom, and bothropstoxin-I, a CRISP molecule, and a D49-PLA2 from the venom of the southern B. jararaca phylogroup. The remarkable paraspecificity exhibited by the Brazilian and the Costa Rican antivenoms indicates large immunoreactive epitope conservation across the natural history of Bothrops, a genus that has its roots in the middle Miocene. This article is part of a Special Issue entitled: Omics Evolutionary Ecolog.

Snake venomics and venom gland transcriptomic analysis of Brazilian coral snakes, Micrurus altirostris and M. corallinus.

The venom proteomes of Micrurus altirostris and M. corallinus were analyzed by combining snake venomics and venom gland transcriptomic surveys. In both coral snake species, 3FTx and PLA(2) were the most abundant and diversified toxin families. 33 different 3FTxs and 13 PLA(2) proteins, accounting respectively for 79.5% and 13.7% of the total proteins, were identified in the venom of M. altirostris. The venom of M. corallinus comprised 10 3FTx (81.7% of the venom proteome) and 4 (11.9%) PLA(2) molecules. Transcriptomic data provided the full-length amino acid sequences of 18 (M. altirostris) and 10 (M. corallinus) 3FTxs, and 3 (M. altirostris) and 1 (M. corallinus) novel PLA(2) sequences. In addition, venom from each species contained single members of minor toxin families: 3 common (PIII-SVMP, C-type lectin-like, L-amino acid oxidase) and 4 species-specific (CRISP, Kunitz-type inhibitor, lysosomal acid lipase in M. altirostris; serine proteinase in M. corallinus) toxin classes. The finding of a lipase (LIPA) in the venom proteome and in the venom gland transcriptome of M. altirostris supports the view of a recruitment event predating the divergence of Elapidae and Viperidae more than 60 Mya. The toxin profile of both M. altirostris and M. corallinus venoms points to 3FTxs and PLA(2) molecules as the major players of the envenoming process. In M. altirostris venom, all major, and most minor, 3FTxs display highest similarity to type I α-neurotoxins, suggesting that these postsynaptically acting toxins may play the predominant role in the neurotoxic effect leading to peripheral paralysis, respiratory arrest, and death. M. corallinus venom posesses both, type I α-neurotoxins and a high-abundance (26% of the venom proteome) protein of subfamily XIX of 3FTxs, exhibiting similarity to bucandin from Malayan krait, Bungarus candidus, venom, which enhances acetylcholine release presynaptically. This finding may explain the presynaptic neurotoxicity of M. corallinus venom and the lack of this effect in M. altirostris venom. The anti-Micrurus (corallinus and frontalis) antivenom produced by Instituto Butantan quantitatively immunodepleted the minor toxins from M. altirostris and M. corallinus venoms but showed impaired crossreactivity towards their major 3FTx and PLA(2) molecules. The structural diversity of 3FTxs among Micrurus sp. may underlay the impaired cross-immunoreactivity of the Butantan antivenom towards M. altirostris and M. corallinus toxins, hampering the possibility to raise an antivenom against a simple venom mixture exhibiting paraspecific neutralization of other Micrurus venoms.

Snake venomics and venom gland transcriptomic analysis of Brazilian coral snakes, Micrurus altirostris and M. corallinus

The venom proteomes of Micrurus altirostris and M. corallinus were analyzed by combining snake venomics and venom gland transcriptomic surveys. In both coral snake species, 3FTx and PLA2 were the most abundant and diversified toxin families. 33 different 3FTxs and 13 PLA2 proteins, accounting respectively for 79.5% and 13.7% of the total proteins, were identified in the venom of M. altirostris. The venom of M. corallinus comprised 10 3FTx (81.7% of the venom proteome) and 4 (11.9%) PLA2 molecules. Transcriptomic data provided the full-length amino acid sequences of 18 (M. altirostris) and 10 (M. corallinus) 3FTxs, and 3 (M. altirostris) and 1 (M. corallinus) novel PLA2 sequences. In addition, venom from each species contained single members of minor toxin families: 3 common (PIII-SVMP, C-type lectin-like, L-amino acid oxidase) and 4 species-specific (CRISP, Kunitz-type inhibitor, lysosomal acid lipase in M. altirostris; serine proteinase in M. corallinus) toxin classes.