Integrated Venomics and Venom Gland Transcriptome Analysis of Juvenile and Adult Mexican Rattlesnakes Crotalus simus, C. tzabcan, and C. culminatus Revealed miRNA-modulated Ontogenetic Shifts.

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA2 crotoxin and acidic PLA2 mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

Constructing comprehensive venom proteome reference maps for integrative venomics.

OBJECTIVE: Understanding the molecular basis of complex adaptive traits, such as snake venom, demands qualitative and quantitative comparisons of the temporal and spatial patterns of venom variation. Here, we assessed the proof-of-concept that locus-resolved reference venom proteome maps can be achieved through efficient pre-MS venom proteome decomplexation, peptide-centric MS/MS analysis and species-specific database searching. METHODS: Venom proteome components were fractionated and quantified by RP-HPLC, SDS-PAGE and 2DE prior to LC-MS/MS matching against a species-specific transcriptomic dataset. RESULTS: Combination of RP-HPLC/SDS-PAGE and 2DE followed by LC-MS/MS showed the existence of ∼178-180 venom protein species generated from ∼48 unique transcripts. CONCLUSIONS: Our results underscore that if sufficient pre-MS and MS efforts are applied, comprehensive venom maps can be achieved. And - equally important - dissociating the venom decomplexing steps from the protein identification process represents the key to achieving a quantitative and locus-resolved insight of the venom proteome.

Venom Proteomics of Trimeresurus gracilis, a Taiwan-Endemic Pitviper, and Comparison of Its Venom Proteome and VEGF and CRISP Sequences with Those of the Most Related Species.

Trimeresurus gracilis is an endemic alpine pitviper in Taiwan with controversial phylogeny, and its venom proteome remains unknown. In this study, we conducted a proteomic analysis of T. gracilis venom using high-performance liquid chromatography-tandem mass spectrometry and identified 155 toxin proteoforms that belong to 13 viperid venom toxin families. By searching the sequences of trypsin-digested peptides of the separated HPLC fractions against the NCBI database, T. gracilis venom was found to contain 40.3% metalloproteases (SVMPs), 15.3% serine proteases, 6.6% phospholipases A2, 5.0% L-amino acid oxidase, 4.6% Cys-rich secretory proteins (CRISPs), 3.2% disintegrins, 2.9% vascular endothelial growth factors (VEGFs), 1.9% C-type lectin-like proteins, and 20.2% of minor toxins, nontoxins, and unidentified peptides or compounds. Sixteen of these proteoforms matched the toxins whose full amino-acid sequences have been deduced from T. gracilis venom gland cDNA sequences. The hemorrhagic venom of T. gracilis appears to be especially rich in PI-class SVMPs and lacks basic phospholipase A2. We also cloned and sequenced the cDNAs encoding two CRISP and three VEGF variants from T. gracilis venom glands. Sequence alignments and comparison revealed that the PI-SVMP, kallikrein-like proteases, CRISPs, and VEGF-F of T. gracilis and Ovophis okinavensis are structurally most similar, consistent with their close phylogenetic relationship. However, the expression levels of some of their toxins were rather different, possibly due to their distinct ecological and prey conditions.

Biological and proteomic characterization of the venom from Peruvian Andes rattlesnake Crotalus durissus.

The Peruvian rattlesnake Crotalus durissus is a venomous species that is restricted to the Peruvian Departments of Puno and Madre de Dios. Although clinically meaningful in this region, Crotalus durissus venom composition remains largely elusive. In this sense, this work aimed to provide a primary description of Peruvian C. durissus venom (PCdV). The enzymatic activities (SVMP, SVSP, LAAO, Hyaluronidase and PLA2) of PCdV were analyzed and compared to Brazilian Crotalus durissus terrificus venom (BCdtV). PCdV showed higher PLA2 activity when compared to the Brazilian venom. PCdV also showed cytotoxicity in VERO cells. For proteomic analysis, PCdV proteins were separated by HPLC, followed by SDS-PAGE. Gel bands were excised and tryptic digested for MALDI-TOF/TOF identification. Approximately 21 proteins were identified, belonging to 7 families. Phospholipases A2 (PLA2, 66.63%) were the most abundant proteins of the venom, followed by snake venom serine proteinases (SVSPs, 13.37%), C-type lectins (Snaclec, 8.98%) and snake venom metalloproteinases (SVMPs, 7.13%), crotamine (2.98%) and phosphodiesterase (PDE, 0.87%). Moreover, antivenom recognition assays indicated that both Brazilian and Peruvian antivenoms recognize PCdV, indicating the presence of antigenically related proteins in crotalic venoms. The results reported here, may impact in the venom selection for the production of effective Pan-American crotalic antivenom.

Biochemical and proteomic analyses of venom from a new pit viper, Protobothrops kelomohy.

Background: A new pit viper, Protobothrops kelomohy, has been recently discovered in northern and northwestern Thailand. Envenoming by the other Protobothrops species across several Asian countries has been a serious health problem since their venom is highly hematotoxic. However, the management of P. kelomohy bites is required as no specific antivenom is available. This study aimed to investigate the biochemical properties and proteomes of P. kelomohy venom (PKV), including the cross-neutralization to its lethality with antivenoms available in Thailand. Methods: PKV was evaluated for its neutralizing capacity (ER50), lethality (LD50), procoagulant and hemorrhagic effects with three monovalent antivenoms (TAAV, DSAV, and CRAV) and one polyvalent (HPAV) hematotoxic antivenom. The enzymatic activities were examined in comparison with venoms of Trimeresurus albolabris (TAV), Daboia siamensis (DSV), Calloselasma rhodostoma (CRV). Molecular mass was separated on SDS-PAGE, then the specific proteins were determined by western blotting. The venom protein classification was analyzed using mass spectrometry-based proteomics. Results: Intravenous LD50 of PKV was 0.67 µg/g. ER50 of HPAV, DSAV and TAAV neutralize PKV at 1.02, 0.36 and 0.12 mg/mL, respectively. PKV exhibited procoagulant effect with a minimal coagulation dose of 12.5 ± 0.016 µg/mL and hemorrhagic effect with a minimal hemorrhagic dose of 1.20 ± 0.71 µg/mouse. HPAV was significantly effective in neutralizing procoagulant and hemorrhagic effects of PKV than those of TAAV, DSAV and CRAV. All enzymatic activities among four venoms exhibited significant differences. PKV proteome revealed eleven classes of putative snake venom proteins, predominantly metalloproteinase (40.85%), serine protease (29.93%), and phospholipase A2 (15.49%). Conclusions: Enzymatic activities of PKV are similarly related to other viperid venoms in this study by quantitatively hematotoxic properties. Three major venom toxins were responsible for coagulopathy in PKV envenomation. The antivenom HPAV was considered effective in neutralizing the lethality, procoagulant and hemorrhagic effects of PKV.

Venomics and antivenomics of Indian spectacled cobra (Naja naja) from the Western Ghats.

Venom proteome profiling of Naja naja from the Western Ghats region in Kerala was achieved through SDS-PAGE and RP-HPLC followed by Q-TOF LC-MS/MS analysis, incorporating PEAKS and Novor assisted de novo sequencing methodologies. A total of 115 proteins distributed across 17 different enzymatic and non-enzymatic venom protein families were identified through conventional and 39 peptides through homology-driven proteomics approaches. Fourteen peptides derived through de novo complements the Mascot data indicating the importance of homology-driven approaches in improving protein sequence information. Among the protein families identified, glutathione peroxidase and endonuclease were reported for the first time in the Indian cobra venom. Immunological cross-reactivity assessed using Indian polyvalent antivenoms suggested that VINS showed better EC50 (2.48 µg/mL) value than that of PSAV (6.04 µg/mL) and Virchow (6.03 µg/mL) antivenoms. Western blotting experiments indicated that all the antivenoms elicited poor binding specificities, especially towards low molecular mass proteins. Second-generation antivenomics studies revealed that VINS antivenom was less efficient to detect many low molecular mass proteins such as three-finger toxins and Kunitz-type serine protease Inhibitors. Taken together, the present study enabled a large-scale characterization of the venom proteome of Naja naja from the Western Ghats and emphasized the need for developing more efficient antivenoms.

A comparison of two different analytical workflows to determine the venom proteome composition of Naja kaouthia from North-East India and immunological profiling of venom against commercial antivenoms.

The Indian monocled cobra (Naja kaouthia) is one of the most prevalent venomous snakes in northeast India (NEI) and is the cause of many fatalities. The composition of NEI N. kaouthia venom (NkV) was deciphered using two different proteomic approaches: (i) 1D SDS-PAGE coupled to label-free quantification of protein bands using stringent identification criteria and (ii) reversed-phase high-performance liquid chromatography (RP-HPLC) followed by quantification based on area under the RP-HPLC peaks. The proteomic data from both strategies were compared. Proteomic analyses from both workflows identified 32 proteins (toxins) distributed over 10-14 snake venom protein families in NEI NkV. The relative abundances of the venom proteins determined from the analytical workflows coincided with the densitometry band intensities of the NEI NkV. Phospholipase A2 (13.1-16.0%) and three-finger toxins (58.5-64.2%) represented the most abundant enzymatic and non-enzymatic proteins in NEI NkV, respectively. Immuno-cross-reactivity studies by enzyme-linked immunoassay and immunoblot analyses pointed to the poor efficacy of commercial PAVs in recognizing the low molecular mass (<15 kDa) toxins of NEI NkV. Spectrofluorometric titration determined the presence of NEI NkV-specific antibodies in commercial PAV, at a level that was higher than that previously reported for eastern India NkV-specific antibodies in commercial antivenom.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India.

The Common Krait (Bungarus caeruleus) shares a distribution range with many other 'phenotypically-similar' kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus' krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer's claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.

Venomics and Cellular Toxicity of Thai Pit Vipers (Trimeresurus macrops and T. hageni).

The two venomous pit vipers, Trimeresurus macrops and T. hageni, are distributed throughout Thailand, although their abundance varies among different areas. No species-specific antivenom is available for their bite victims, and the only recorded treatment method is a horse antivenom raised against T. albolabris crude venom. To facilitate assessment of the cross-reactivity of heterologous antivenoms, protein profiles of T. macrops and T. hageni venoms were explored using mass-spectrometry-based proteomics. The results show that 185 and 216 proteins were identified from T. macrops and T. hageni venoms, respectively. Two major protein components in T. macrops and T. hageni venoms were snake venom serine protease and metalloproteinase. The toxicity of the venoms on human monocytes and skin fibroblasts was analyzed, and both showed a greater cytotoxic effect on fibroblasts than monocytic cells, with toxicity occurring in a dose-dependent rather than a time-dependent manner. Exploring the protein composition of snake venom leads to a better understanding of the envenoming of prey. Moreover, knowledge of pit viper venomics facilitates the selection of the optimum heterologous antivenoms for treating bite victims.

Delineating the venom toxin arsenal of Malabar pit viper (Trimeresurus malabaricus) from the Western Ghats of India and evaluating its immunological cross-reactivity and in vitro cytotoxicity.

The venom protein components of Malabar pit viper (Trimeresurus malabaricus) were identified by combining SDS-PAGE and ion-exchange chromatography pre-fractionation techniques with LC-MS/MS incorporating Novor and PEAKS-assisted de novo sequencing strategies. Total 97 proteins that belong to 16 protein families such as L-amino acid oxidase, metalloprotease, serine protease, phospholipase A2, 5'-nucleotidase, C-type lectins/snaclecs and disintegrin were recognized from the venom of a single exemplar species. Of the 97 proteins, eighteen were identified through de novo approaches. Immunological cross-reactivity assessed through ELISA and western blot indicate that the Indian antivenoms binds less effectively to Malabar pit viper venom components compared to that of Russell's viper venom. The in vitro cell viability assays suggest that compared to the normal cells, MPV venom induces concentration dependent cell death in various cancer cells. Moreover, crude venom resulted in chromatin condensation and apoptotic bodies implying the induction of apoptosis. Taken together, the present study enabled in dissecting the venom proteome of Trimeresurus malabaricus and revealed the immuno-cross-reactivity profiles of commercially available Indian polyvalent antivenoms that, in turn, is expected to provide valuable insights on the need in improving antivenom preparations against its bite.

Phylovenomics of Daboia russelii across the Indian subcontinent. Bioactivities and comparative in vivo neutralization and in vitro third-generation antivenomics of antivenoms against venoms from India, Bangladesh and Sri Lanka.

Russell's viper (Daboia russelii) is, together with Naja naja, Bungarus caeruleus and Echis carinatus, a member of the medically important 'Big Four' species responsible for causing a large number of morbidity and mortality cases across the Indian subcontinent. Despite the wide distribution of Russell's viper and the well-documented ubiquity of the phenomenon of geographic variability of intraspecific snake venom composition, Indian polyvalent antivenoms against the "Big Four" venoms are raised against venoms sourced mainly from Chennai in the southeastern Indian state of Tamil Nadu. Biochemical and venomics investigations have consistently revealed notable compositional, functional, and immunological differences among geographic variants of Russell's viper venoms across the Indian subcontinent. However, these studies, carried out by different laboratories using different protocols and involving venoms from a single geographical region, make the comparison of the different venoms difficult. To bridge this gap, we have conducted bioactivities and proteomic analyses of D. russelii venoms from the three corners of the Indian subcontinent, Pakistan, Bangladesh, and Tamil Nandu (India) and Sri Lanka, along with comparative in vivo neutralization and in vitro third-generation antivenomics of antivenoms used in India, Bangladesh and Sri Lanka. These analyses let us to propose two alternative routes of radiation for Russell's viper in the Indian subcontinent. Both radiations, towards the northeast of India and Bangladesh and towards south India and Sri Lanka, have a common origin in Pakistan, and provide a phylovenomics ground for rationalizing the geographic variability in venom composition and their distinct immunoreactivity against available antivenoms. BIOLOGICAL SIGNIFICANCE: Russell's viper (Daboia russelii), the Indian cobra (Naja naja), the common krait (Bungarus caeruleus), and the saw-scaled viper (Echis carinatus) constitute the 'Big Four' snake species responsible for most snakebite envenomings and deaths in the Indian subcontinent. Despite the medical relevance of Daboia russelii, and the well documented variations in the clinical manifestations of envenomings by this wide distributed species, which are doubtless functionally related to differences in venom composition of its geographic variants, antivenoms for the clinical treatment of envenomings by D. russelii across the Indian subcontinent are invariably raised using venom sourced mainly from the southeastern Indian state of Tamil Nadu. We have applied a phylovenomics approach to compare the venom proteomes of Russell's vipers from the three corners of the Indian subcontinent, Pakistan, Bangladesh, and South India/Sri Lanka, and have assessed the in vitro (third-generation antivenomics) and in vivo preclinical efficacy of a panel of homologous antivenoms. The identification of two dispersal routes of ancestral D. russelii into the Indian subcontinent provides the ground for rationalizing the variability in composition and immunoreactivity of the venoms of extant geographic variants of Russell's viper. Such knowledge is relevant for envisioning strategies to improve the clinical coverage of anti- D. russelii antivenoms.

Proteomic profile, biological activities and antigenic analysis of the venom from Bothriopsis bilineata smaragdina ("loro machaco"), a pitviper snake from Peru.

In order to determine Bothriopsis bilineata smaragdina venom (BbsV) composition, proteomic approaches were performed. Venom components were analyzed by RP-HPLC, SDS- PAGE and nano LC on line with LTQ Orbitrap XL. Results showed a total of 189 identified proteins, grouped into 11 different subgroups, which include snake venom metalloproteinases (SVMPs, 54.67%), snake C-type lectins (Snaclecs, 15.78%), snake venom serine proteinases (SVSPs, 14.69%), cystein-rich secretory proteins (CRISP, 2.61%), phospholipases A2 (PLA2, 1.14%), phosphodiesterase (PDE, 1.17%), venom endothelial growth factor (VEGF, 1.06%) 5'nucleotidases (0.33%), L-amino acid oxidases (LAAOs, 0.28%) and other proteins. In vitro enzymatic activities (SVMP, SVSP, LAAO, Hyal and PLA2) of BbsV were also analyzed. BbsV showed high SVSP activity but low PLA2 activity, when compared to other Bothrops venoms. In vivo, BbsV induced hemorrhage and edema in mice and showed intraperitoneal median lethal dose (LD50) of 92.74 (± 0.15) µg/20 g of mice. Furthermore, BbsV reduced cell viability when incubated with VERO cells. Peruvian and Brazilian bothropic antivenoms recognize BbsV proteins, as detected by ELISA and Western Blotting. Both antivenoms were able to neutralize in vivo edema and hemorrhage. SIGNIFICANCE: In Peru, snakebite is a public health problem, especially in the rain forest, as a result of progressive colonization of this geographical area. This country is the second in Latin America, after Brazil, to exhibit the largest variety of venomous snakes. B. atrox and B. b. smaragdina snakes are sympatric species in Peruvian Amazon region and are responsible for approximately 95% of the envenomings reported in this region. B. b. smaragdina may cause a smaller share (3 to 38%) of those accidents, due to its arboreal habits, that make human encounters with these snakes less likely to happen. Despite B. b. smaragdina recognized medical importance, its venom composition and biological activities have been poorly studied. Furthermore, BbsV is not a component of the antigenic pool used to produce the corresponding Peruvian bothropic antivenom (P-BAV). Our results not only provide new insights on BbsV composition and biological activity, but also demonstrate that both P-BAV and B-BAV polyvalent antivenoms have a considerable recognition of proteins from BbsV and, more importantly, neutralized hemorrhage and edema, the main local effects of bothropic envenomation.

Mass spectrometry-assisted venom profiling of Hypnale hypnale found in the Western Ghats of India incorporating de novo sequencing approaches.

Hypnale hypnale (hump-nosed pit viper) is considered to be one among the medically important venomous snake species of India and Sri Lanka. In the present study, venom proteome profiling of a single Hypnale hypnale from Western Ghats of India was achieved using SDS-PAGE based protein separation followed by LC-MS/MS analysis. The identities of the proteins that were not established using the Mascot search were determined through de novo sequencing tools such as Novor followed by MS-BLAST based sequence similarity search algorithm and PEAKS proteomics software. The combined proteomics analysis revealed a total of 37 proteins belonging to nine different snake venom families, in which 7 proteins were exclusively identified through de novo strategies. The enzymatic and non-enzymatic venom protein families identified include serine proteases, metalloproteases, phospholipase A2, thrombin-like enzymes, phospholipase B, C-type lectins/snaclecs, disintegrins, cysteine rich secretory proteins and nerve growth factor. Among these, disintegrins, nerve growth factor, phospholipase B and cysteine rich secretory protein families were identified for the first time in HPV venom. This could possibly explain the regiospecific venom variation seen across snake species. Taken together, the venom proteome profiling on Indian Hypnale hypnale venom correlates with the clinical manifestations often seen in the envenomed victims.

Snake Venomics Display: An online toolbox for visualization of snake venomics data.

With the introduction of powerful mass spectrometry equipment into the field of snake venom proteomics, a large body of venomics data is accumulating. To allow for better comparison between venom compositions from different snake species and to provide an online database containing this data, we devised the Snake Venomics Display toolbox for visualization of snake venomics data on linear scales. This toolbox is freely available to be used online at https://tropicalpharmacology.com/tools/snake-venomics-display/ and allows researchers to visualize venomics data in a Relative Abundance (%) visualization mode and in an Absolute Abundance (mg) visualization mode, the latter taking venom yields into account. The curated venomics data for all snake species included in this database is also made available in a downloadable Excel file format. The Snake Venomics Display toolbox represents a simple way of handling snake venomics data, which is better suited for large data sets of venom compositions from multiple snake species.

Transcriptomics-guided bottom-up and top-down venomics of neonate and adult specimens of the arboreal rear-fanged Brown Treesnake, Boiga irregularis, from Guam.

The Brown Treesnake (Boiga irregularis) is an arboreal, nocturnal, rear-fanged venomous snake native to northern and eastern regions of Australia, Papua New Guinea and the Solomon Islands. It was inadvertently introduced onto the island of Guam during the late 1940's to early 1950's, and it has caused massive declines and extirpations of the native bird, lizard, and mammal populations. In the current study, we report the characterization of the venom proteome of an adult and a neonate B. irregularis specimens from Guam by a combination of venom gland transcriptomic and venomic analyses. Venom gland transcriptomic analysis of an adult individual identified toxins belonging to 18 protein families, with three-finger toxin isoforms being the most abundantly expressed transcripts, comprising 94% of all venom protein transcript reads. Transcripts for PIII-metalloproteinases, C-type lectins, cysteine-rich secretory proteins, acetylcholinesterases, natriuretic peptides, ficolins, phospholipase A2 (PLA2) inhibitors, PLA2s, vascular endothelial growth factors, Kunitz-type protease inhibitors, cystatins, phospholipase Bs, cobra venom factors, waprins, SVMP inhibitors, matrix metalloproteinases, and hyaluronidases were also identified, albeit, at very low abundances ranging from 0.05% to 1.7% of the transcriptome. The venom proteomes of neonate and adult B. irregularis were also both overwhelmingly (78 and 84%, respectively) dominated by monomeric and dimeric 3FTxs, followed by moderately abundant (21% (N) and 13% (A)) CRISPs, low abundance (1% (N), 3% (A)) PIII-SVMPs, and very low abundance (<0.01%) PLA2 and SVMP inhibitors. The differences in relative toxin abundances identified between neonate and adult snakes likely correlates to shifts in prey preference between the two age classes, from nearly-exclusively lizards to lizards, birds and small mammals. Immunoaffinity antivenomics with experimentally designed rabbit anti-Brown Treesnake (anti-BTS) venom IgGs against homologous venom from adult snakes demonstrated that CRISPs, PIII-SVMPs, and 60-70% of 3FTxs were effectively immunocaptured. Western blot analysis showed that all venom proteins were recognized by anti-BTS IgGs, and cross-reactivity with other rear-fanged snake venoms was also observed. Incubation of anti-BTS venom IgGs with crude B. irregularis venom resulted in a significant decrease in proteolytic (SVMP) activity against azocasein. These results provide the first comparative venomic and anti-venomic analysis of neonate and adult B. irregularis from Guam, further highlighting evolutionary trends in venom composition among rear-fanged venomous snakes. SIGNIFICANCE PARAGRAPH: The Brown Treesnake (Boiga irregularis) has caused extensive ecological and economic damage to the island of Guam where it has become a classic example of the negative impacts of invasive species. In the current study, we report the first combined transcriptomic and proteomic analysis of B. irregularis venom of Guam origin. The transcriptome of an adult snake contained toxin sequences belonging to 18 protein families, with three-finger toxin (3FTx) isoforms being the most abundant and representing 94% of all venom protein transcript reads. Our bottom-up and top-down venomic analyses confirmed that 3FTxs are the major components of B. irregularis venom, and a comparative analysis of neonate and adult venoms demonstrate a clear ontogenetic shift in toxin abundance, likely driven by dietary variation between the two age classes. Second-generation antivenomics and Western blot analysis using purified anti-Brown Treesnake rabbit serum IgGs (anti-BTS IgGs) showed strong immunoreactivity toward B. irregularis venom. Interestingly, our anti-BTS IgGs did not cross-react with 3FTxs found in several other rear-fanged snake venoms, or against 3FTxs in the venom of the elapid Ophiophagus hannah, indicating that epitopes in these 3FTx molecules are quite distinct.

Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis).

Bothriechis is a genus of eleven currently recognized slender and arboreal venomous snakes, commonly called palm-pitvipers that range from southern Mexico to northern South America. Despite dietary studies suggesting that palm-pitvipers are generalists with an ontogenetic shift toward endothermic prey, venom proteomic analyses have revealed remarkable divergence between the venoms of the Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. To achieve a more complete picture of the venomic landscape across Bothriechis, the venom proteomes of biodiversity of the northern Middle American highland palm-pitvipers, B. thalassinus, B. aurifer, and B. bicolor from Guatemala, B. marchi from Honduras, and neonate Costa Rican B. lateralis and B. schlegelii, were investigated. B. thalassinus and B. aurifer venoms are comprised by similar toxin arsenals dominated by SVMPs (33-39% of the venom proteome), CTLs (11-16%), BPP-like molecules (10-13%), and CRISPs (5-10%), and are characterized by the absence of PLA2 proteins. Conversely, the predominant (35%) components of B. bicolor are D49-PLA2 molecules. The venom proteome of B. marchi is similar to B. aurifer and B. thalassinus in that it is rich in SVMPs and BPPs, but also contains appreciable amounts (14.3%) of PLA2s. The major toxin family found in the venoms of both neonate B. lateralis and B. schlegelii, is serine proteinase (SVSP), comprising about 20% of their toxin arsenals. The venom of neonate B. schlegelii is the only palm-pitviper venom where relative high amounts of Kunitz-type (6.3%) and γPLA2 (5.2%) inhibitors have been identified. Despite notable differences between their proteomes, neonate venoms are more similar to each other than to adults of their respective species. However, the ontogenetic changes taking place in the venom of B. lateralis strongly differ from those that occur in the venom of B. schlegelii. Thus, the ontogenetic change in B. lateralis produces a SVMP-rich venom, whereas in B. schlegelii the age-dependent compositional shift generates a PLA2-rich venom. Overall, genus-wide venomics illustrate the high evolvability of palm-pitviper venoms. The integration of the pattern of venom variation across Bothriechis into a phylogenetic and biogeographic framework may lay the foundation for assessing, in future studies, the evolutionary path that led to the present-day variability of the venoms of palm-pitvipers. SIGNIFICANCE: Bothriechis represents a monophyletic basal genus of eleven arboreal palm-pitvipers that range from southern Mexico to northern South America. Despite palm-pitvipers' putative status as diet generalists, previous proteomic analyses have revealed remarkable divergence between the venoms of Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. Our current proteomic study of Guatemalan species, B. thalassinus, B. aurifer, and B. bicolor, Honduran B. marchi, and neonate B. lateralis and B. schlegelii from Costa Rica was undertaken to deepen our understanding of the evolutionary pattern of venom proteome diversity across Bothriechis. Ancestral characters are often, but not always, preserved in an organism's development. Venoms of neonate B. lateralis and B. schlegelii are more similar to each other than to adults of their respective species, suggesting that the high evolvability of palm-pitviper venoms may represent an inherent feature of Bothriechis common ancestor. Our genus-wide data identified four nodes of venom phenotype differentiation across the phylogeny of Bothriechis. Integrated into a phylogenetic and biogeographic framework, the pattern of venom variation across Bothriechis may lay the groundwork to establish whether divergence was driven by selection for efficient resource exploitation in arboreal 'islands', thereby contributing to the ecological speciation of the genus.

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.

Venomics and antivenomics of Bothrops erythromelas from five geographic populations within the Caatinga ecoregion of northeastern Brazil.

The Caatinga lancehead, Bothrops erythromelas, is a medically relevant species, responsible for most of the snakebite accidents in most parts of its distribution range in northeastern Brazil. The spectrum and geographic variability of its venom toxins were investigated applying a venomics approach to venom pools from five geographic areas within the Caatinga ecoregion. Despite its wide habitat, populations of B. erythromelas from Ceará, Pernambuco, Juazeiro, Paraiba, and Ilha de Itaparica exhibit highly conserved venom proteomes. Mirroring their compositional conservation, the five geographic venom pools also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated in Vital Brazil (BR) and Clodomiro Picado (CR) Institutes, using different venoms in the immunization mixtures. The paraspecificity exhibited by the Brazilian SAB and the Costa Rican BCL antivenoms against venom toxins from B. erythromelas indicates large immunoreactive epitope conservation across genus Bothrops during the last ~14 million years, thus offering promise for the possibility of generating a broad-spectrum bothropic antivenom. Biological Significance Accidental snakebite envenomings represent an important public health hazard in Brazil. Ninety per cent of the yearly estimated 20-30,000 snakebite accidents are caused by species of the Bothrops genus. Bothrops erythromelas, a small, moderately stocky terrestrial venomous snake, is responsible for most of the snakebite accidents in its broad distribution range in the Caatinga, a large ecoregion in northeastern Brazil. To gain a deeper insight into the spectrum of medically important toxins present in the venom of the Caatinga lancehead, we applied a venomics approach to define the proteome and geographic variability of adult B. erythromelas venoms from five geographic regions. Although intraspecific compositional variation between venoms among specimens from different geographic regions has long been appreciated by herpetologists and toxinologists as a general feature of highly adaptable and widely distributed snake species, the five B. erythromelas populations investigated exhibit highly conserved venom proteomes. The overall toxin profile of the Caatinga lancehead's venom explains the local and systemic effects observed in envenomations by B. erythromelas. The five geographic venom pools sampled also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated using different bothropic venoms in the immunization mixtures. The large immunoreactive epitope conservation across genus Bothrops offers promise for the generation of a broad-spectrum bothropic antivenom.

Unveiling the nature of black mamba (Dendroaspis polylepis) venom through venomics and antivenom immunoprofiling: Identification of key toxin targets for antivenom development.

The venom proteome of the black mamba, Dendroaspis polylepis, from Eastern Africa, was, for the first time, characterized. Forty- different proteins and one nucleoside were identified or assigned to protein families. The most abundant proteins were Kunitz-type proteinase inhibitors, which include the unique mamba venom components 'dendrotoxins', and α-neurotoxins and other representatives of the three-finger toxin family. In addition, the venom contains lower percentages of proteins from other families, including metalloproteinase, hyaluronidase, prokineticin, nerve growth factor, vascular endothelial growth factor, phospholipase A2, 5'-nucleotidase, and phosphodiesterase. Assessment of acute toxicity revealed that the most lethal components were α-neurotoxins and, to a lower extent, dendrotoxins. This venom also contains a relatively high concentration of adenosine, which might contribute to toxicity by influencing the toxin biodistribution. ELISA immunoprofiling and preclinical assessment of neutralization showed that polyspecific antivenoms manufactured in South Africa and India were effective in the neutralization of D. polylepis venom, albeit showing different potencies. Antivenoms had higher antibody titers against α-neurotoxins than against dendrotoxins, and displayed high titers against less toxic proteins of high molecular mass. Our results reveal the complexity of D. polylepis venom, and provide information for the identification of its most relevant toxins to be neutralized by antivenoms. BIOLOGICAL SIGNIFICANCE: The black mamba, D. polylepis, is one of the most feared snakes in the world, owing to the potency of its venom, the severity and rapid onset of clinical manifestations of envenomings, and its ability to strike fast and repeatedly. The present study reports the first proteomic analysis of this venom. Results revealed a complex venom constituted predominantly by proteins belonging to the Kunitz-type proteinase inhibitor family, which comprises the dendrotoxins, and to α-neurotoxins of the three-finger toxin family. The proteins showing highest acute toxicity were α-neurotoxins, which induce post-synaptic blockade of the neuromuscular junctions, followed by dendrotoxins, which inhibit the voltage-dependent potassium channels. The combination of these two types of toxins in the venom underscores the presence of a dual strategy that results in a highly effective mechanism for prey subduction. This complex toxic arsenal is likely to provide D. polylepis with high trophic versatility. The rapid onset and severity of neurotoxic clinical manifestations in envenomings by D. polylepis demand the rapid administration of effective and safe antivenoms. Preclinical tests showed that an antivenom from South Africa and two antivenoms from India were effective in the neutralization of this venom, albeit differing in their potency. Moreover, ELISA immunoprofiling of these antivenoms against all venom fractions revealed that antivenoms have higher titers against α-neurotoxins than against dendrotoxins, thus underscoring the need to develop improved immunization strategies. The results of this investigation identified the most relevant toxins present in D. polylepis venom, which need to be targeted by antivenoms or other type of inhibitors.

Comparative venomics of the Prairie Rattlesnake (Crotalus viridis viridis) from Colorado: Identification of a novel pattern of ontogenetic changes in venom composition and assessment of the immunoreactivity of the commercial antivenom CroFab®.

Here we describe and compare the venomic and antivenomic characteristics of both neonate and adult Prairie Rattlesnake (Crotalus viridis viridis) venoms. Although both neonate and adult venoms contain unique components, similarities among protein family content were seen. Both neonate and adult venoms consisted of myotoxin, bradykinin-potentiating peptide (BPP), phospholipase A2 (PLA2), Zn(2+)-dependent metalloproteinase (SVMP), serine proteinase, L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRISP) and disintegrin families. Quantitative differences, however, were observed, with venoms of adults containing significantly higher concentrations of the non-enzymatic toxic compounds and venoms of neonates containing higher concentrations of pre-digestive enzymatic proteins such as SVMPs. To assess the relevance of this venom variation in the context of snakebite and snakebite treatment, we tested the efficacy of the common antivenom CroFab® for recognition of both adult and neonate venoms in vitro. This comparison revealed that many of the major protein families (SVMPs, CRISP, PLA2, serine proteases, and LAAO) in both neonate and adult venoms were immunodepleted by the antivenom, whereas myotoxins, one of the major toxic components of C. v. viridis venom, in addition to many of the small peptides, were not efficiently depleted by CroFab®. These results therefore provide a comprehensive catalog of the venom compounds present in C. v. viridis venom and new molecular insight into the potential efficacy of CroFab® against human envenomations by one of the most widely distributed rattlesnake species in North America. BIOLOGICAL SIGNIFICANCE: Comparative proteomic analysis of venoms of neonate and adult Prairie Rattlesnake (Crotalus viridis viridis) from a discrete population in Colorado revealed a novel pattern of ontogenetic shifts in toxin composition for viperid snakes. The observed stage-dependent decrease of the relative content of disintegrins, catalytically active D49-PLA2s, L-amino acid oxidase, and SVMPs, and the concomitant increase of the relative abundance of paralytic small basic myotoxins and ohanin-like toxin, and hemostasis-disrupting serine proteinases, may represent an age-dependent strategy for securing prey and avoiding injury as the snake switches from small ectothermic prey and newborn rodents to larger endothermic prey. Such age-dependent shifts in venom composition may be relevant for antivenom efficacy and treatment of snakebite. However, applying a second-generation antivenomics approach, we show that CroFab®, developed against venom of three Crotalus and one Agkistrodon species, efficiently immunodepleted many, but not all, of the major compounds present in neonate and adult C. v. viridis venoms.