Next Article in Journal
Rapid, Sensitive and Reliable Ricin Identification in Serum Samples Using LC–MS/MS
Next Article in Special Issue
Potency Testing of Venoms and Antivenoms in Embryonated Eggs: An Ethical Alternative to Animal Testing
Previous Article in Journal
Transconjunctival versus Transcutaneous Injection of Botulinum Toxin into the Lacrimal Gland to Reduce Lacrimal Production: A Randomized Controlled Trial
Previous Article in Special Issue
A Neurotoxic Snake Venom without Phospholipase A2: Proteomics and Cross-Neutralization of the Venom from Senegalese Cobra, Naja senegalensis (Subgenus: Uraeus)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 13
Issue 2
10.3390/toxins13020078
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pan-American Lancehead Pit-Vipers: Coagulotoxic Venom Effects and Antivenom Neutralisation of Bothrops asper and B. atrox Geographical Variants

by
Lachlan A. Bourke
1,
Christina N. Zdenek
1,
Edgar Neri-Castro
2,
Melisa Bénard-Valle
2,
Alejandro Alagón
2,
José María Gutiérrez
3,
Eladio F. Sanchez
4,
Matt Aldridge
5 and
Bryan G. Fry
1,*
1
Toxin Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, QLD 4072, Australia
2
Departamento de Medicina Molecular y Bioprocesos, Instituto de Biotecnologa, Universidad Nacional Autónoma de México, Av. Universidad 2001, Cuernavaca, Morelos 62210, Mexico
3
Instituto Clodomiro Picado, Facultad de Microbiología, Universidad de Costa Rica, San José 11501, Costa Rica
4
Laboratory of Biochemistry of Proteins from Animal Venoms, Research and Development Center, Ezequiel Dias Foundation, Belo Horizonte, MG 30510-010, Brazil
5
MicroPharm Limited, Station Road Ind. Est., Station Road, Newcastle Emlyn, Carmarthenshire SA38 9BY, UK
*
Author to whom correspondence should be addressed.
Toxins 2021, 13(2), 78; https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13020078
Submission received: 21 December 2020 / Revised: 10 January 2021 / Accepted: 17 January 2021 / Published: 22 January 2021
(This article belongs to the Special Issue Interactions of Snake Venoms and Antivenoms: Prelude to Protection)
Browse Figures
Versions Notes

Abstract

:
The toxin composition of snake venoms and, thus, their functional activity, can vary between and within species. Intraspecific venom variation across a species’ geographic range is a major concern for antivenom treatment of envenomations, particularly for countries like French Guiana that lack a locally produced antivenom. Bothrops asper and Bothrops atrox are the most medically significant species of snakes in Latin America, both producing a variety of clinical manifestations, including systemic bleeding. These pathophysiological actions are due to the activation by the venom of the blood clotting factors Factor X and prothrombin, thereby causing severe consumptive coagulopathy. Both species are extremely wide-ranging, and previous studies have shown their venoms to exhibit regional venom variation. In this study, we investigate the differential coagulotoxic effects on human plasma of six venoms (four B. asper and two B. atrox samples) from different geographic locations, spanning from Mexico to Peru. We assessed how the venom variation of these venom samples affects neutralisation by five regionally available antivenoms: Antivipmyn, Antivipmyn-Tri, PoliVal-ICP, Bothrofav, and Soro Antibotrópico (SAB). The results revealed both inter- and intraspecific variations in the clotting activity of the venoms. These variations in turn resulted in significant variation in antivenom efficacy against the coagulotoxic effects of these venoms. Due to variations in the venoms used in the antivenom production process, antivenoms differed in their species-specific or geographical neutralisation capacity. Some antivenoms (PoliVal-ICP, Bothrofav, and SAB) showed species-specific patterns of neutralisation, while another antivenom (Antivipmyn) showed geographic-specific patterns of neutralisation. This study adds to current knowledge of Bothrops venoms and also illustrates the importance of considering evolutionary biology when developing antivenoms. Therefore, these results have tangible, real-world implications by aiding evidence-based design of antivenoms for treatment of the envenomed patient. We stress that these in vitro studies must be backed by future in vivo studies and clinical trials before therapeutic guidelines are issued regarding specific antivenom use in a clinical setting.
Keywords:
Bothrops; coagulotoxicity; antivenom neutralisation; antivenom; venom variation
Key Contribution: Due to differential biochemical compositions of venoms used in the immunisation process, Latin American antivenoms vary in their intraspecific and interspecific cross-reactivity.

1. Introduction

Venomous snakebites across the Americas result in an estimated 652 to 3466 deaths per year [1], with the majority of these snakebites by pit-vipers [2,3,4,5]. Bothrops is a genus of pit-viper—commonly known as lanceheads—that is distributed across South and Central America, as well as Mexico. Epidemiological data suggest that across South America, 70–96.6% of envenomings are inflicted by the Bothrops species [3,4,5,6]. Most of these bites are from B. atrox in Amazon regions [6,7] or B. asper in northern South America [8]. Similar values are also reported across Central America, with 50–80% of all snakebites attributed to B. asper [8,9]. Epidemiological data on Mexican B. asper bites are scarce, although the incidence is approximately 6 bite cases per 100,000 people per year [8]. In a more recent study, however, Chippaux [10] reported 3.3 snakebite cases per 100,000 people per year in Mexico, although the species responsible for these bites is not made clear. Thus, collectively across the Americas, B. asper and B. atrox may be considered the two species responsible for most medically significant snakebite envenomations.
Coagulopathy resulting from over-activation of the blood clotting cascade is one of the most common clinical manifestations induced by Bothrops venoms [7,8,11]. Bites from such coagulotoxic snakes are particularly difficult to treat, as they cause disruption of an integral bodily function: the body’s blood clotting mechanism (the haemostatic system). The formation of a blood clot ultimately results in the production of a fibrin mesh network that traps red blood cells to arrest bleeding [12]. Fibrin is the end-product of a process known as the coagulation cascade. In this cascading series of biochemical events in the blood, a clot is formed through sequential activation of blood coagulation factors (zymogens). In the final stages of this domino-type effect, prothrombin is converted into its active form (thrombin), which then cleaves fibrinogen into fibrin strands that cross-link to form a blood clot [12].
Numerous snake venoms immobilise prey by disrupting the coagulation cascade in various ways, including via activation or inhibition of coagulation factors [13,14,15,16]. Procoagulant venoms induce clotting by either possessing activated forms of clotting factors in the venom itself or by activating endogenous zymogens, such as prothrombin, Factor V or Factor X [13,15]. Procoagulant venoms ultimately produce large amounts of endogenous thrombin, leading to a well-ordered, fibrin clot. Snake venom metalloproteases (SVMPs) have been shown to be responsible for the prothrombin and Factor X activation produced by Bothrops venoms [7,11,17,18], with differential activity between SVMP isoforms in relative prothrombin or Factor X activation [19]. The possession of such toxins makes Bothrops venoms deadly, underscoring the need for antivenom treatment for envenomed patients. Depending on species and geographic region, a multitude of different antivenoms is used to treat Bothrops envenomations. These antivenoms include, among others, Antivipmyn (Instituto Bioclon, Mexico) [20], Central American polyspecific antivenom (PoliVal-ICP; Instituto Clodomiro Picado, Costa Rica) [21,22], Antivipmyn-Tri (Instituto Bioclon, Mexico) [21,23,24], Soro Antibotrópico (SAB; Instituto Butantan, Brazil) [25], and Bothrofav (MicroPharm, United Kingdom) [21]. Despite antivenom being available to treat envenomations, a fundamental, life-threatening issue with antivenom efficacy arises due to the natural extreme complexity of venoms that have been shaped by a myriad of selection pressures. In fact, diversity among snake venoms (both inter- and intraspecifically) is ubiquitous, with the main driver of venom variation often attributed to diet differences [26,27]. Diet-related selection pressures on venom can include prey type, prey escape potential, and prey retaliation potential [11,28]. However, neutral evolution processes, such as genetic drift, genetic bottleneck, or founder effect, are additional drivers of venom variation [29,30].
Many studies have shown intraspecific venom variation in B. asper and B. atrox [11,31,32,33]. In fact, venoms from the Bothrops genus (specifically B. asper and B. atrox) are highly malleable to intra- and interspecific venom variation, due to their evolutionary history and geographic range. During the Miocene epoch, between 23–10 million years ago (Mya), the common ancestor of the Bothrops + Bothrocophias (sister genera to Bothrops) clade colonised the then viperid-devoid South America, with the Bothrops genus splitting off subsequent to the invasion of South America, where it rapidly diversified [34]. Approximately 3.34 Mya, B. atrox and B. asper diverged from a common ancestor into separate species, and soon after, B. asper reinvaded Middle America in at least two independent dispersal events [35,36]. Today, B. atrox occurs throughout South America, from Brazil to Northern South America, while B. asper occurs across Mexico, Central, and South America [37,38]. Such a vast geographic range of both species could result in geographically isolated populations under different selection pressures, thereby differentially shaping venom evolution. This is likely the driver of B. asper and B. atrox venom variation observed in previous studies [11,31,32,36]. Furthermore, evidence of distinct lineages within B. atrox suggests the species is probably a complex of multiple species [7,39,40]. This could also be the case for B. asper, as previous work suggests [35,41]. Speciation occurs due to different selection pressures acting on each species and, thus, would undoubtedly affect venom composition.
Despite a multitude of antivenoms being available for treating envenomations by species of the same genus, there has not been a standardised comparison of them against the two widest-ranging medically important species in Latin America (B. asper and B. atrox), which thus impedes the evidence-based design of clinical management strategies. It is especially important to test the use of antivenoms made for snakes in one region for their cross-reactivity for snakebites in regions that do not have a locally produced antivenom. Use of non-local antivenom is seen in French Guiana, which, due to their lack of a local antivenom, adopted Mexico’s Antivipmyn-Tri to treat snakebites, including B. atrox snakebites, in the region [23]. However, the literature on the neutralising ability of Antivipmyn-Tri for French Guiana snakes is scarce [23,24,42,43], and, therefore, there is a knowledge gap regarding its efficacy in treating critically envenomated patients. Cross-species neutralisation has, however, been noted for other antivenoms, such as SAB from Instituto Butantan, which does not contain B. atrox in the immunising mixture but has been proven to be effective in neutralising the coagulotoxic venoms of some B. atrox populations [11]. It is therefore relevant to assess whether antivenoms available in Latin America are effective in the neutralisation of venoms of B. asper and B. atrox from various geographic locations [44].
The aim of this study, therefore, is to investigate venom in vitro clotting activity variation and antivenom efficacy across geographically distinct populations of B. asper and B. atrox using five antivenoms that are currently used to treat Bothrops bites. It is hypothesised that there will be venom variation both within and between species and that this will result in differential antivenom efficacy. This knowledge will be fundamental for improving antivenom design in Latin America.

2. Results

The spontaneous clotting control of plasma and the positive kaolin control was 459.90 ± 67.13 s (n = 4) and 63.08 ± 0.67 s (n = 4), respectively. Antivenoms did not have any effect on clotting as antivenom control values were similar to the spontaneous clotting time: Antivipmyn (n = 4) = 485.63 ± 54.41 s, PoliVal-ICP (n = 4) = 483.15 ± 27.93 s, Antivipmyn-Tri (n = 4) = 486.45 ± 20.97 s, Bothrofav (n = 4) = 521.45 ± 9.32 s, and SAB (n = 4) = 467.08 ± 30.79 s. A Brown–Forsythe ANOVA test revealed no significant differences between antivenom controls and the spontaneous clotting time (F (5.00, 9.21) = 1.13, p = 0.41).
All venoms exhibited extremely strong procoagulant effects on plasma, with clotting times at 20 µg/mL venom concentration (at which the concentration of all venoms was well into the asymptotic phase and, thus, at maximum clotting velocity), ranging from 13.50 ± 0.85 s for B. asper (Ecuador) to 47.20 ± 2.04 s for B. asper (San Andres, Tuxtla, (S.A.T.), Mexico) (Table 1, Figure 1A and Figure 2A). Both intra- and interspecific variations were observed. Based upon the area under the curve (AUC) values (Table 1), the relative ranking of the venoms from most potent (smallest AUC) to least potent (largest AUC) was B. asper (Yucatán, Mexico) 377.33 ± 9.15; B. asper (Ecuador) 384.00 ± 5.86; B. asper (Costa Rica) 481.70 ± 18.01; B. atrox (French Guiana) 669.70 ± 15.85; B. atrox (Alto Marañon (A.M.), Peru) 741.73 ± 5.47; B. asper (S.A.T., Mexico) 1,327.33 ± 26.03. Considerable intraspecies variation was observed in B. asper AUC values across localities, with B. asper (S.A.T., Mexico) AUC being 72% (100 − (377.33/1327.33) * 100) greater (less potent) than B. asper (Yucatán, Mexico) (Table 1, Figure 1A). These difference between Mexican B. asper AUC values were significant (Tukey’s multiple comparisons test, p < 0.0001; Table 2).
Amongst all the antivenoms tested, PoliVal-ICP was the best performing antivenom for all B. asper samples except B. asper (S.A.T., Mexico; Figure 1). This antivenom significantly induced the highest X-fold shift (slowing of clotting times due to venom neutralisation) among all antivenoms for these samples, with the other antivenoms having very low neutralisation compared to PoliVal-ICP (Table 3). In contrast, B. asper (S.A.T., Mexico) was best neutralised by Bothrofav and SAB (X-fold shift = 1.83 ± 0.06 and 1.74 ± 0.11, respectively; Figure 1), and the X-fold shift values were significantly higher than all other antivenoms including PoliVal-ICP (Table 3).
Across all the localities of B. atrox, SAB was the best performing antivenom, with the X-fold shift significantly higher than all other antivenoms tested (Figure 2, Table 4). Bothrofav also neutralised French Guiana and A.M. (Peru) localities in a similar way to SAB (Figure 2). All other antivenoms had lower neutralisation compared to the best performing antivenom (Figure 2). Variation in antivenom neutralisation was observed between Mexican B. asper samples and between B. atrox samples. The S.A.T (Mexico) locality was the best neutralised across all the antivenoms compared to the Yucatán (Mexico) locality (Figure 3.). A similar result was observed for the B. atrox samples, with the A.M. (Peru) locality being the best neutralised across all the antivenoms compared to the French Guiana locality (Figure 3).
A clear species-specific pattern of neutralisation was observed for PoliVal-ICP, with all B. asper samples better neutralised than B. atrox samples (Figure 3). These differences were significant: B. asper samples had significantly higher X-fold shift values than B. atrox samples, with the exception being B. asper (Yucatán, Mexico) as it was not significantly different from B. atrox (A.M., Peru; Table 5). A species-specific pattern was also observed for SAB and Bothrofav, which significantly neutralised all B. atrox samples to a greater extent than B. asper samples, except for B. asper (S.A.T., Mexico; Figure 3, Table 5). Antivipmyn displayed a possible geographically specific pattern of neutralisation: all Central American and Mexican samples were better neutralised than the South American samples (Figure 3). These differences were significant: Central American and Mexican samples had significantly higher X-fold shift values than South American samples, with the exception being B. asper (Yucatán, Mexico) as it was not significantly different from B. atrox (A.M., Peru; Table 6). No clear species-specific or geographically specific pattern of neutralisation was observed for Antivipmyn-Tri (Figure 3).

3. Discussion

The aim of this study was to investigate the venom effects of multiple B. asper and B. atrox localities on human plasma and to ascertain the relative efficacy of five regionally important antivenoms used to treat Bothrops bites. It was hypothesised that both intra- and interspecific venom variation would be observed and would result in differential antivenom efficacy. Indeed, the venoms tested differed in coagulotoxic activity both intra- and interspecifically and in antivenom neutralising potential (Figure 1, Figure 2 and Figure 3, Table 1). SAB and Bothrofav performed well against B. atrox venoms (Figure 2, Table 4), while PoliVal-ICP neutralised all B. asper samples considerably better than other antivenoms, except for B. asper (S.A.T., Mexico; Figure 1, Table 3). PoliVal-ICP, Bothrofav, and SAB showed species-specific patterns of neutralisation, while Antivipmyn showed a geographically specific pattern of neutralisation (Figure 3).
Venom from all Bothrops tested displayed potent procoagulant activity on our assays (Figure 1A and Figure 2A, Table 1). In prey animals, this procoagulant action results in rapid incapacitation by the induction of stroke. However, due to the dilution of the venom into a much large blood volume, thrombotic strokes in humans from Bothrops envenomations are rare but have been reported [45]. The more common outcome of the out-of-control activation of the clotting cascade is unclottable blood and systemic bleeding [7,8,25,46], due to venom induced consumptive coagulopathy [47]. An exception to this general trend is the venom of B. lanceolatus, endemic to the Caribbean island of Martinique, which is more likely than other Bothrops species to cause severe thrombosis and infarctions [48].
Despite the consistent procoagulant activity of all species, our study found significant differences in AUC values (representing clotting activity; smaller AUC values indicate greater procoagulant activity, and vice versa) across different localities (ANOVA: F (5, 12) = 1644.00, p < 0.0001; Figure 1A and Figure 2A, Table 1 and Table 2). This result is consistent with previous literature, which found B. asper and B. atrox venom activity to vary between different localities [11,31,32,33].
Clotting times induced by B. asper venoms were similar across Costa Rica, Ecuador, and Yucatán (Mexico) localities. The S.A.T. (Mexico) locality, however, had considerably slower clotting times (Figure 1A, Table 1). Unlike the other samples, our B. asper Mexican samples were individual venoms rather than pooled venoms (Table 7). Individual venom variation has been observed in multiple snake species [49,50,51]; thus pooled venoms are favoured as they help account for any possible variation between individuals. Therefore, in the present study, individual variation cannot be ruled out as contributing to the observed variation in venom potency and antivenom efficacy. Future work should investigate whether the differences between Mexican B. asper samples in this study are due to differences in venoms between individual snakes or differences in venom between populations of snakes from different localities. If the latter, drivers of population variation should be investigated through natural history observations.
As with procoagulant activity, variation amongst antivenom efficacy was evident, both between and within species, and across different antivenoms (Figure 1, Figure 2 and Figure 3). A possible geographical pattern of antivenom neutralisation across both Bothrops species tested was observed for Antivipmyn, whereby most Central American and Mexican venom samples, with the exception being B. asper (Yucatán, Mexico), were significantly more neutralised than South American venom samples (Figure 3, Table 6). This geographical pattern is evident, as the only South American B. asper sample (B. asper (Ecuador)) was the least neutralised B. asper sample by this antivenom. This result indicates different toxin compositions between South American and Central/North American B. asper, possibly driven by founder effects following the invasion of Central America from South America and the resulting genetic bottleneck imposed by the transit through the very narrow Isthmus of Panama. Furthermore, the venom variation observed may correspond with the Central American and Mexican B. asper samples in this study being of a different lineage than the Ecuador B. asper sample. Two divergent lineages of B. asper have been recognised: one comprising snakes from southwestern Ecuador and the Caribbean coast of Middle America (Clade A), and the other comprising snakes from northwestern South America and Pacific regions of Costa Rica and Panama (Clade B) [35]. In fact, venom variation has been observed in Costa Rican B. asper from Pacific and Caribbean regions [32,33]. Thus, differences observed between B. asper samples may reflect their geographic location and evoluntionary history.
A clear species-specific pattern was observed for PoliVal-ICP, since most B. asper venom samples were significantly more neutralised than B. atrox venom samples (Figure 3, Table 5). This suggests a high cross-reactivity of the antivenom to B. asper localities besides just the B. asper locality of Costa Rica, which is used in the antivenom immunizing mixture of PoliVal-ICP (see Table 8). PoliVal-ICP also neutralised Mexican B. asper samples better than Antivipmyn (Tukey’s multiple comparisons test, p < 0.0001; Figure 1, Table 3), which is another antivenom with B. asper venom in the immunising mixture (Table 8). However, it is important to note that only two individual venom samples from Mexico were used in the study; thus, more comprehensive sampling must be undertaken before evidence-based medical recommendations can be made. Furthermore, it is important to note that Antivipmyn significantly neutralised B. asper from S.A.T. (Mexico) better than B. asper from Yucatán (Mexico) (Tukey’s multiple comparisons test: p = 0.0158; Figure 3). Thus, not only do these venoms differ in their coagulotoxic profiles, but these differences also reflect antivenom neutralisation. To further elucidate the efficacy of Antivipmyn and the drivers of this variation, as seen amongst Mexican samples, more work needs to be done in order to better inform on treatment strategies for the envenomed patient.
SAB and Bothrofav also showed species-specific effects by neutralising all B. atrox samples to a greater extent than B. asper samples, with the exception being the high neutralisation of Bothrofav to the S.A.T. (Mexico) locality (Figure 3, Table 5). Despite SAB antivenom not having B. atrox in the immunising mixture (but using related species; see Table 8), its strong cross-reactivity with B. atrox is consistent with other studies [11]. The cross-reactivity of antivenoms similar to SAB (produced with similar species) with B. atrox has also been observed [52]. In this regard, proteomic analysis of several Bothrops venoms has been reported, such as the protein composition of the Peruvian B. atrox included in this work, and correlated with their main pharmacological properties [53]. It has become clear that a low number of toxic compounds belonging to the predominant protein families (snake venom metalloproteases, serine proteases, phospholipases A2 and l-amino acid oxidases), among others, compose the venom of Bothrops species [36,53]. Thus, data of the literature and results obtained in this study provide evidence of the degree of antigenic conservation among Bothrops snakes.
On the other hand, this is the first time the strong cross-reactivity of Bothrofav antivenom against B. atrox venoms has been shown. The high neutralisation of Bothrofav to B. asper (S.A.T., Mexico), but not B. asper (Yucatán, Mexico) (Figure 1 and Figure 3), further indicates important differences in venom composition between these samples. A previous study has described the lack of efficacy of a first-generation Bothrofav antivenom for neutralising in vitro coagulant activity of venoms of B. asper and B. atrox [54]. It is likely that changes have been introduced in the manufacturing of more recent batches of this antivenom, which would explain its efficacy in the neutralisation of procoagulant activity of B. atrox venom, as shown in our work.
Overall, these results help inform antivenom treatment plans for regions without antivenoms containing local snake venom in the immunising mixture. Antivipmyn-Tri is currently used in French Guiana to treat snakebites [23], although our results show it has low neutralising efficacy for the coagulotoxic effects of B. atrox inhabiting French Guiana and the other B. atrox venoms tested in this study. These results are congruent with a recent retrospective study that looked at the efficacy of Antivipmyn-Tri and concluded that compared to patients receiving no antivenom, it did not benefit patient recovery from coagulopathy as a result of envenomation [24]. Bothrofav and SAB perform better at neutralising the procoagulant activity of B. atrox from French Guiana. Furthermore, Bothrofav antivenom is used to treat B. lanceolatus, a species inhabiting the Caribbean island of Martinique. However, in comparison to other antivenoms used, it showed high efficacy in neutralising the coagulotoxic effect of all B. atrox studied, as well as B. asper from S.A.T. (Mexico) (Figure 1, Figure 2 and Figure 3). This finding is interesting in the light of the low procoagulant activity described for the venom of B. lanceolatus [54,55]. However, it must be noted that these studies examined the coagulation of citrated plasma without recalcifying it. It has been shown that recalcification is essential to ascertain procoagulant toxins [11]. Thus, the previous lack of demonstrated procoagulant activity of B. lanceolatus in laboratory studies is likely due to methodological conditions, especially in light of the devastating procoagulant activity that this species has been documented as having in clinical reports [45,48]. Furthermore, the nature of venom components in B. lanceolatus that elicit the antibody response capable of neutralising the procoagulant effect of B. atrox venom deserves further investigation.
It is important to note the limitations of the study so that appropriate conclusions can be drawn and future work mapped out. Firstly, all antivenom efficacy results are based on the ability of the antivenom to neutralise the coagulotoxic effects of the venom on platelet-poor plasma and not other venom effects, such as the activation and inhibitory effects on platelet aggregation [17,56,57,58], and myotoxicity [59,60]. Thus, even if an antivenom neutralises procoagulant coagulotoxic effects, other medically important venom effects may not be neutralised. Secondly, another caveat in the study is the use of some antivenoms beyond their expiry date. However, it has been documented that antivenoms, especially freeze-dried ones, maintain their neutralising effect well beyond their expiry date, with studies showing that their efficacy was retained even decades later (up to 60 years in one study) [61,62]. Furthermore, one study showed no significant change in specific antibody concentration when liquid-formulated antivenoms were stored at 4 °C vs. ambient temperature [63], hence showing how robust some antivenoms can be. Thus, while a quantitative reduction in their neutralising ability might occur, as observed in O’Leary, Kornhauser, Hodgson and Isbister [62], this does not affect the main qualitative conclusions reached in this study regarding the species-specific and geographically specific patterns of neutralisation by the antivenoms. Thirdly, proteomic work was not performed in this study; thus, differences in venom function and resultant antivenom efficacy cannot be discussed in the context of differences in proteins. Proteomic work should be investigated in futures studies to elucidate the differences in venom components and how these impact antivenom efficacy.
Another important caveat of the present study is that the experiments are in vitro. In vivo systems are far more complicated (e.g., involving multiple physiological systems and feedback loops) and further studies should confirm the efficacy of the antivenoms used in this study on in vivo models, assessing, for example, the neutralisation of hemorrhagic and defibrinogenating activities. However, while success in this study needs to be corroborated with in vivo follow-up work, it is unlikely for an antivenom to perform better in vivo than in vitro where variables are removed and the antivenom is given an ideal situation (i.e., in close proximity with no flow) to bind to the toxins. Thus, antivenoms that failed to neutralise a venom under our study conditions would be exceedingly unlikely to neutralise the venom in a more dynamic in vivo scenario. Therefore, indications of poor performance obtained under the ideal conditions are strongly suggestive of poor performance in real-world scenarios. Conversely, for antivenoms that worked well against a particular venom in the in vitro assays of this study, it is critical that such results are confirmed with in vivo assays and clinical trials before recommendations are made regarding therapeutic use. Therefore, while the overall patterns of antivenom efficacy observed in this current study are important and informative, additional work in vivo and in clinical trials are required before reaching more conclusive therapeutic recommendations.
In conclusion, the results of the present study are important not only to evolutionary biologists interested in venom variation, but also to antivenom manufacturers and clinicians treating snakebites in Latin America. The intraspecific variation of Bothrops and other snake species needs to be considered when developing and using antivenoms. Studies like the present add to the body of knowledge required for the evidence-based design of clinical management strategies.

4. Materials and Methods

4.1. Snake Venom Sample Preparation

All work was performed under University of Queensland IBSC approval #IBC134BSBS2015. All snake venom samples were sourced from either venom suppliers or the Toxin Evolution Lab collection at the University of Queensland. Samples were stored in lyophilised form until use. In this study, six venom samples from the Bothrops genus (four B. asper and two B. atrox samples), including multiple geographic localities, were used (Table 7). It was not a concern that some venom samples were from captive snakes and some from wild snakes, as studies on Bothrops have shown captivity does not influence venom variation [64]. For experiments, lyophilised venom samples were solubilised with deionised H2O and their protein concentration measured in triplicate using Thermo Fisher Scientific™ NanoDrop 2000 (Waltham, MA, USA) at 280 nm wavelength. Venom samples were then made into working stocks of 1 mg/mL 1:1 deionised H2O:glycerol and stored (frozen) at −80 °C until required, whereby they were then stored (unfrozen) at −20 °C. To avoid protein degradation, glycerol was used to stabilise the enzymes, thus preventing degradation due to autocatalysis.

4.2. Human Plasma

All human plasma work was performed under University of Queensland Biosafety Approval #IBC134BSBS2015 and Human Ethics Approval #2016000256. Surplus human platelet-poor plasma (3.2% citrated) was supplied by the Australian Red Cross (44 Musk Street, Kelvin Grove, Queensland 4059) under research approval #16-04QLD-10 and stored at −80 °C until required. Platelet-poor plasma, as opposed to platelet-rich plasma, was used due to its long-term stability. Furthermore, platelet-poor plasma is commonly used in coagulation studies, as platelets can interact with reagents and shorten clotting times [65]. For the present study, one batch of platelet-poor plasma (Batch # A540019384586, O positive, collected 3 May 2019) was defrosted at 37 °C and aliquoted into 1.5 mL tubes in a laminar flow hood to avoid contamination. Aliquots were flash-frozen with liquid nitrogen and stored at −80 °C. Platelet-poor plasma was defrosted for 5 min in a 37 °C water bath before being used in assays (see Section 4.3 and Section 4.4).

4.3. Venom-Induced Coagulation

Coagulation assays were performed on a Stago STA-R Max haemostasis analyser (Stago, Asnières sur Seine, France) using a previously established assay [66,67,68]. The clotting time of platelet-poor plasma was measured in triplicate at eight different venom concentrations (μg/mL: 0.05, 0.125, 0.25, 0.66, 1.66, 4, 10, and 20), using all venoms (Table 7). For testing, venom stock was diluted 1:10 using Owren Koller (OK) buffer (Stago Cat #00360) to make a 0.1 mg/mL solution. This working solution was added into the analyser. For the 1:1 dilution (20 μg/mL), the machine added 50 μL of the venom sample into a cuvette with cofactors (50 μL of 25 mM CaCl2 (Stago Cat# 00367), 50 μL of phospholipid (Stago Cat# 00597), and 25 µL of OK buffer). The cuvette was then incubated for 120 s at 37 °C. After incubation, 75 µL of platelet-poor plasma was added into the incubated venom mix (final cuvette volume = 250 µL) and clotting time was automatically measured via a mechanical, viscosity-based system. For subsequent venom concentrations, the volume of the venom sample added into the cuvette was adjusted and additional OK buffer was added to maintain a final cuvette volume of 250 µL. Positive and negative controls were performed to ensure the health of the plasma and reagents. A kaolin test (similar to aPTT (activated partial thromboplastin time)) was used as the positive control: 50 µL of platelet-poor plasma incubated with 50 µL kaolin (Stago kit; Cat# 00597), 50 µL phospholipid, and 50 µL OK buffer for 120 s. Then, 25 mM CaCl2 was added postincubation to initiate clotting. The negative control used was the spontaneous clotting time of plasma, in which the same volume of reagents was used as the venom coagulation assay, except a 1:10 diluted (diluent = OK buffer) 50% deionised H2O/50% glycerol sample was used instead of the venom sample. Dilution series were not performed for controls. All tests were performed in triplicate, and clotting time was measured up until 999 s (machine maximum reading time).

4.4. Antivenom Neutralising Efficacy

Five antivenoms that are used to treat Bothrops envenomings were included in the study (Table 8). Expired antivenoms were not a concern, as antivenoms have been shown to still be effective after their expiry date [61,62]. In particular, freeze-dried antivenoms are highly stable for many years after the official expiration date. Antivipmyn and Antivipmyn-Tri were both freeze-dried antivenoms, while all other antivenoms used were in liquid form. All antivenoms were centrifuged at 12,000 relative centrifugal force (RCF) on an Allegra™ X-22R Centrifuge (Lot#982501, Beckman Coulter, Brea, CA, USA) at 4 °C for 10 min; then, the supernatant was removed and filtered with a 0.45 µm filter, aliquoted (1 mL into 1.5 mL tubes), and stored at 4 °C. Antivenom coagulation assays were performed for each antivenom to determine their potential to neutralise the coagulotoxic effects of each venom. These assays were performed on a Stago STA-R Max haemostasis analyser, as previously described [66,67,68], whereby the same procedure as the venom coagulation assays (see Section 4.3) was used, except the 25 μL OK buffer was replaced with 25 μL of 2.5% (v/v) antivenom (final cuvette concentration of 0.25%). Prior to testing, the antivenom was prepared to 2.5% by diluting with OK buffer. A negative control for antivenom assays was performed, in which a 1:10 diluted (diluent = OK buffer) 50% deionised H2O/50% glycerol sample was used instead of the venom sample. This control was performed to see if the antivenom had any effect on clotting.

4.5. Statistical Analyses

All raw data is available in Supplementary Dataset S1. All statistical analyses were performed in GraphPad PRISM 8.1.1 (GraphPad Prism Inc., La Jolla, CA, USA). A p-value of less than 0.05 was considered significant. When reporting p-values and the associated information, the following abbreviations are used: p means p-value and F means F-value. When reporting F-values, the degrees of freedom are shown in brackets. All results in the study are shown as the mean ± standard deviation.
Eight-point dilution curves, showing the clotting time of each venom and venom + antivenom across a venom concentration gradient, were graphed using GraphPad PRISM 8.1.1. To more clearly view the data, the x-axis for eight-point dilution curves was presented in logarithmic view. To calculate the neutralising capacity of each venom, an area under the curve (AUC) value was obtained in GraphPad PRISM 8.1.1 for each of the venom curves and venom + antivenom curves. An X-fold shift value was calculated with the following formula: venom   +   antivenom   AUC   Venom   AUC − 1 ; a value of 0 indicates no neutralisation and a value above 0 indicates neutralisation. Bar graphs of X-fold shift values were graphed using GraphPad PRISM 8.1.1.
Venom sample AUC values were compared using an ordinary one-way analysis of variance (ANOVA; Table 2). Within each venom, X-fold shifts for different antivenoms were also statistically compared using an ordinary one-way ANOVA (Table 3 and Table 4). Assumptions of the normality of residuals and the homodescasity of residuals were checked using quantile–quantile plots, Shapiro–Wilk and Brown–Forsythe tests. The posthoc test used was Tukey’s multiple comparisons test, in which multiplicity adjusted p-values were used that accounted for multiple comparisons. The same statistical procedure was followed to statistically compare X-fold shifts between different venoms within each antivenom (Table 5 and Table 6). To test if antivenoms had any effect on clotting times, the clotting times of antivenom controls were statistically compared to the spontaneous control using the Brown–Forsythe ANOVA. A Brown–Forsythe ANOVA was used as the data violated the assumption of the homodescasity of residuals. A Welch ANOVA was also performed alongside the Brown–Forsythe test due to uncertainty as to which test is best. Despite the p-value being barely significant (p = 0.0485), all posthoc tests (Dunnett’s T3 multiple comparisons test) were extremely not-significant (ranging from p = 0.4038 to p = 0.9997); thus, it was deemed that there was no statistical difference between antivenom control and spontaneous control clotting times.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/13/2/78/s1, Dataset S1: Raw data.

Author Contributions

Conceptualization, B.G.F.; methodology, B.G.F.; formal analysis L.A.B. and B.G.F.; investigation, L.A.B., C.N.Z., E.F.S. and B.G.F.; resources, E.N.-C., M.B.-V., A.A., and J.M.G.; writing—original draft preparation, L.A.B., C.N.Z., E.F.S. and B.G.F.; writing—review and editing, L.A.B., C.N.Z., E.N.-C., M.B.-V., A.A., J.M.G., M.A., and B.G.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by a grant from Australian Discovery Project DP190100304 to B.G.F.

Institutional Review Board Statement

University of Queensland IBSC approval #IBC134BSBS2015 (5 January 2018).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Supplementary Dataset S1.

Conflicts of Interest

J.M.G. works at Instituto Clodomiro Picado (Universidad de Costa Rica), where the antivenom PoliVal-ICP is manufactured, while M.A. works at Micropharm, where BothroFav antivenom is manufactured. However, their involvement in this study was completely independent of the companies, with the companies having no oversight of this study.

References

  1. Kasturiratne, A.; Wickremasinghe, A.R.; de Silva, N.; Gunawardena, N.K.; Pathmeswaran, A.; Premaratna, R.; Savioli, L.; Lalloo, D.G.; de Silva, H.J. The global burden of snakebite: A literature analysis and modelling based on regional estimates of envenoming and deaths. PLoS Med. 2008, 5, e218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ruha, A.-M.; Kleinschmidt, K.C.; Greene, S.; Spyres, M.B.; Brent, J.; Wax, P.; Padilla-Jones, A.; Campleman, S. The epidemiology, clinical course, and management of snakebites in the North American snakebite registry. J. Med. Toxicol. 2017, 13, 309–320. [Google Scholar] [CrossRef] [PubMed]
  3. Dolab, J.A.; de Roodt, A.R.; de Titto, E.H.; García, S.I.; Funes, R.; Salomón, O.D.; Chippaux, J.-P. Epidemiology of snakebite and use of antivenom in Argentina. Trans. R. Soc. Trop. Med. Hyg. 2014, 108, 269–276. [Google Scholar] [CrossRef] [PubMed]
  4. Da Costa, M.K.B.; da Fonseca, C.S.; Navoni, J.A.; Freire, E.M.X. Snakebite accidents in Rio Grande do Norte state, Brazil: Epidemiology, health management and influence of the environmental scenario. Trop. Med. Int. Health 2019, 24, 432–441. [Google Scholar] [CrossRef]
  5. Chippaux, J.-P. Epidemiology of envenomations by terrestrial venomous animals in Brazil based on case reporting: From obvious facts to contingencies. J. Venom. Anim. Toxins Incl. Trop. Dis. 2015, 21, 1–17. [Google Scholar] [CrossRef] [Green Version]
  6. Mota-da-Silva, A.; Colombini, M.; Moura-da-Silva, A.M.; Souza, R.M.; Monteiro, W.M.; Bernarde, P.S. Epidemiological and clinical aspects of snakebites in the upper Juruá River region, western Brazilian Amazonia. Acta Amaz. 2019. [Google Scholar] [CrossRef]
  7. Monteiro, W.M.; Contreras-Bernal, J.C.; Bisneto, P.F.; Sachett, J.; Mendonça da Silva, I.; Lacerda, M.; Guimarães da Costa, A.; Val, F.; Brasileiro, L.; Sartim, M.A.; et al. Bothrops atrox, the most important snake involved in human envenomings in the amazon: How venomics contributes to the knowledge of snake biology and clinical toxinology. Toxicon X 2020, 6, 100037. [Google Scholar] [CrossRef]
  8. Otero-Patiño, R. Epidemiological, clinical and therapeutic aspects of Bothrops asper bites. Toxicon 2009, 54, 998–1011. [Google Scholar] [CrossRef]
  9. Gutiérrez, J.M. Current challenges for confronting the public health problem of snakebite envenoming in Central America. J. Venom. Anim. Toxins Incl. Trop. Dis. 2014, 20, 7. [Google Scholar] [CrossRef] [Green Version]
  10. Chippaux, J.-P. Incidence and mortality due to snakebite in the Americas. PLoS Negl. Trop. Dis. 2017, 11, e0005662. [Google Scholar] [CrossRef]
  11. Sousa, L.F.; Zdenek, C.N.; Dobson, J.S.; Op den Brouw, B.; Coimbra, F.C.P.; Gillett, A.; Del-Rei, T.H.M.; Chalkidis, H.D.M.; Sant’Anna, S.; Teixeira-da-Rocha, M.M.; et al. Coagulotoxicity of Bothrops (lancehead pit-vipers) venoms from Brazil: Differential biochemistry and antivenom efficacy resulting from prey-driven venom variation. Toxins 2018, 10, 411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Colman, R.W.; Clowes, A.W.; George, J.N.; Goldhaber, S.Z.; Marder, V.J. Overview of hemostasis. In Hemostasis and Thrombosis: Basic Principles and Clinical Practice; Colman, R.W., Clowes, A.W., George, J.N., Goldhaber, S.Z., Marder, V.J., Eds.; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2006; pp. 3–20. [Google Scholar]
  13. Kini, R.M. The intriguing world of prothrombin activators from snake venom. Toxicon 2005, 45, 1133–1145. [Google Scholar] [CrossRef] [PubMed]
  14. Kini, R.M. Anticoagulant proteins from snake venoms: Structure, function and mechanism. Biochem. J. 2006, 397, 377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kini, R.M.; Rao, V.S.; Joseph, J.S. Procoagulant proteins from snake venoms. Pathophysiol. Haemost. Thromb. 2001, 31, 218–224. [Google Scholar] [CrossRef] [PubMed]
  16. Markland, F.S. Snake venoms and the hemostatic system. Toxicon 1998, 36, 1749–1800. [Google Scholar] [CrossRef]
  17. Loría, G.D.; Rucavado, A.; Kamiguti, A.S.; Theakston, R.D.G.; Fox, J.W.; Alape, A.; Gutiérrez, J.M.A. Characterization of ‘basparin A,’ a prothrombin-activating metalloproteinase, from the venom of the snake Bothrops asper that inhibits platelet aggregation and induces defibrination and thrombosis. Arch. Biochem. Biophys. 2003, 418, 13–24. [Google Scholar] [CrossRef]
  18. Rucavado, A.; Escalante, T.; Gutiérrez, J.M.a. Effect of the metalloproteinase inhibitor batimastat in the systemic toxicity induced by Bothrops asper snake venom: Understanding the role of metalloproteinases in envenomation. Toxicon 2004, 43, 417–424. [Google Scholar] [CrossRef]
  19. Sousa, L.F.; Bernardoni, J.L.; Zdenek, C.N.; Dobson, J.; Coimbra, F.; Gillett, A.; Lopes-Ferreira, M.; Moura-da-Silva, A.M.; Fry, B.G. Differential coagulotoxicity of metalloprotease isoforms from Bothrops neuwiedi snake venom and consequent variations in antivenom efficacy. Toxicol. Lett. 2020, 333, 211–221. [Google Scholar] [CrossRef]
  20. Fuchs, J.; Faber, K.; Tuchscherer, D.T.; Tsakiris, D.A.; Weiler, S.; Hofer, K.E. Bite by a juvenile Bothrops venezuelensis (Venezuelan lancehead) resulting in severe envenomation: A case report. Toxicon 2020, 180, 39–42. [Google Scholar] [CrossRef]
  21. Resiere, D.; Mehdaoui, H.; Gutiérrez, J.M. Snakebite envenomation in the Caribbean: The role of medical and scientific cooperation. PLoS Negl. Trop. Dis. 2018, 12, e0006441. [Google Scholar] [CrossRef]
  22. Gutiérrez, J.M.; Lomonte, B.; Sanz, L.; Calvete, J.J.; Pla, D. Immunological profile of antivenoms: Preclinical analysis of the efficacy of a polyspecific antivenom through antivenomics and neutralization assays. J. Proteom. 2014, 105, 340–350. [Google Scholar] [CrossRef] [PubMed]
  23. Nadaud, A.; Perotti, F.; de Haro, L.; Boels, D. Snake envenomations in French Guiana: First clinical assessment of an antivenom imported from Mexico. Anaesth. Crit. Care Pain Med. 2019, 38, 193–194. [Google Scholar] [CrossRef] [PubMed]
  24. Heckmann, X.; Lambert, V.; Mion, G.; Ehrhardt, A.; Marty, C.; Perotti, F.; Carod, J.-F.; Jolivet, A.; Boels, D.; Lehida Andi, I.; et al. Failure of a Mexican antivenom on recovery from snakebite-related coagulopathy in French Guiana. Clin. Toxicol. 2020, 1–7. [Google Scholar] [CrossRef]
  25. De Silva Oliveira, S.; Campos Alves, E.; dos Santos Santos, A.; Freitas Nascimento, E.; Tavares Pereira, J.P.; da Mendonça Silva, I.; Sachett, J.; dos Santos Ibiapina, H.N.; Santos Sarraf, L.K.; Contreras Bernal, J.C.; et al. Bothrops snakebites in the Amazon: Recovery from hemostatic disorders after Brazilian antivenom therapy. Clin. Toxicol. 2020, 58, 266–274. [Google Scholar] [CrossRef] [PubMed]
  26. Barlow, A.; Pook, C.E.; Harrison, R.A.; Wüster, W. Coevolution of diet and prey-specific venom activity supports the role of selection in snake venom evolution. Proc. R. Soc. Biol. Sci. Ser. B 2009, 276, 2443–2449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Daltry, J.C.; Wüster, W.; Thorpe, R.S. Diet and snake venom evolution. Nature 1996, 379, 537–540. [Google Scholar] [CrossRef]
  28. Jackson, T.N.W.; Koludarov, I.; Ali, S.A.; Dobson, J.; Zdenek, C.N.; Dashevsky, D.; Op den Brouw, B.; Masci, P.P.; Nouwens, A.; Josh, P.; et al. Rapid radiations and the race to redundancy: An investigation of the evolution of Australian elapid snake venoms. Toxins 2016, 8, 309. [Google Scholar] [CrossRef]
  29. Sasa, M. Diet and snake venom evolution: Can local selection alone explain intraspecific venom variation? Toxicon 1999, 37, 253–260. [Google Scholar] [CrossRef]
  30. Williams, V.; White, J.; Schwaner, T.D.; Sparrow, A. Variation in venom proteins from isolated populations of tiger snakes (Notechis ater niger, N. scutatus) in South Australia. Toxicon 1988, 26, 1067–1075. [Google Scholar] [CrossRef]
  31. Calvete, J.J.; Sanz, L.; Pérez, A.; Borges, A.; Vargas, A.M.; Lomonte, B.; Angulo, Y.; Gutiérrez, J.M.; Chalkidis, H.M.; Mourão, R.H.V.; et al. Snake population venomics and antivenomics of Bothrops atrox: Paedomorphism along its transamazonian dispersal and implications of geographic venom variability on snakebite management. J. Proteom. 2011, 74, 510–527. [Google Scholar] [CrossRef]
  32. Alape-Girón, A.; Sanz, L.; Escolano, J.; Flores-Díaz, M.; Madrigal, M.; Sasa, M.; Calvete, J.J. Snake venomics of the lancehead pitviper Bothrops asper: Geographic, individual, and ontogenetic variations. J. Proteome Res. 2008, 7, 3556–3571. [Google Scholar] [CrossRef] [PubMed]
  33. Aragón, F.; Gubenšek, F. Bothrops asper venom from the Atlantic and Pacific zones of Costa Rica. Toxicon 1981, 19, 797–805. [Google Scholar] [CrossRef]
  34. Wüster, W.; Salomão, M.; Quijada-Mascareñas, J.A.; Thorpe, R.; BBBSP. Origins and evolution of the South American pitviper fauna: Evidence from mitochondrial DNA sequence analysis. In Biology of the Vipers; NHBS Academic & Professional Books: Eagle Mountain, UT, USA, 2002; pp. 111–128. [Google Scholar]
  35. Saldarriaga-Córdoba, M.; Parkinson, C.L.; Daza, J.M.; Wüster, W.; Sasa, M. Phylogeography of the Central American lancehead Bothrops asper (Serpentes: Viperidae). PLoS ONE 2017, 12, e0187969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Mora-Obando, D.; Salazar-Valenzuela, D.; Pla, D.; Lomonte, B.; Guerrero-Vargas, J.A.; Ayerbe, S.; Gibbs, H.L.; Calvete, J.J. Venom variation in Bothrops asper lineages from North-Western South America. J. Proteom. 2020, 229, 103945. [Google Scholar] [CrossRef] [PubMed]
  37. Nogueira, C.C.; Argôlo, A.J.S.; Arzamendia, V.; Azevedo, J.A.; Barbo, F.E.; Bérnils, R.S.; Bolochio, B.E.; Borges-Martins, M.; Brasil-Godinho, M.; Braz, H.; et al. Atlas of Brazilian Snakes: Verified point-locality maps to mitigate the Wallacean shortfall in a megadiverse snake fauna. S. Am. J. Herpetol. 2019, 14, 1–274. [Google Scholar] [CrossRef]
  38. Campbell, J.A.; Lamar, W.W. The Venomous Reptiles of the Western Hemisphere; Comstock Publishing Associates: Ithaca, NY, USA, 2004; Volume 1. [Google Scholar]
  39. Gibbs, H.L.; Sovic, M.; Amazonas, D.; Chalkidis, H.; Salazar-Valenzuela, D.; Moura-Da-Silva, A.M. Recent lineage diversification in a venomous snake through dispersal across the Amazon River. Biol. J. Linn. Soc. 2018, 123, 651–665. [Google Scholar] [CrossRef] [Green Version]
  40. Wüster, W.; Salomão, M.G.; Duckett, G.J.; Thorpe, R.S. BBBSP Mitochondrial DNA phylogeny of the Bothrops atrox species complex (Squamata: Serpentes: Viperidae). Kaupia Darmst Beitr Natgesch 1999, 8, 135–144. [Google Scholar]
  41. Salazar Valenzuela, C.D. Diversification in the Neotropics: Insights from Demographic and Phylogenetic Patterns of Lancehead Pitvipers (Bothrops spp.). Ph.D. Thesis, The Ohio State University, Columbus, OH, USA, 2016. [Google Scholar]
  42. Resiere, D.; Villalta, M.; Arias, A.S.; Kallel, H.; Nèviére, R.; Vidal, N.; Mehdaoui, H.; Gutiérrez, J.M. Snakebite envenoming in French Guiana: Assessment of the preclinical efficacy against the venom of Bothrops atrox of two polyspecific antivenoms. Toxicon 2020, 173, 1–4. [Google Scholar] [CrossRef]
  43. Resiere, D.; Houcke, S.; Pujo, J.M.; Mayence, C.; Mathien, C.; NkontCho, F.; Blaise, N.; Demar, M.P.; Hommel, D.; Kallel, H. Clinical features and management of snakebite envenoming in French Guiana. Toxins 2020, 12, 662. [Google Scholar] [CrossRef]
  44. Segura, A.; Castillo, M.C.; Núñez, V.; Yarlequé, A.; Gonçalves, L.R.C.; Villalta, M.; Bonilla, C.; Herrera, M.; Vargas, M.; Fernández, M.; et al. Preclinical assessment of the neutralizing capacity of antivenoms produced in six Latin American countries against medically-relevant Bothrops snake venoms. Toxicon 2010, 56, 980–989. [Google Scholar] [CrossRef]
  45. Thomas, L.; Chausson, N.; Uzan, J.; Kaidomar, S.; Vignes, R.; Plumelle, Y.; Bucher, B.; Smadja, D. Thrombotic stroke following snake bites by the “Fer-de-Lance” Bothrops lanceolatus in Martinique despite antivenom treatment: A report of three recent cases. Toxicon 2006, 48, 23–28. [Google Scholar] [CrossRef] [PubMed]
  46. Otero, R.; Gutiérrez, J.; Beatriz Mesa, M.a.; Duque, E.; Rodríguez, O.; Luis Arango, J.; Gómez, F.; Toro, A.; Cano, F.; María Rodríguez, L.; et al. Complications of Bothrops, Porthidium, and Bothriechis snakebites in Colombia. A clinical and epidemiological study of 39 cases attended in a university hospital. Toxicon 2002, 40, 1107–1114. [Google Scholar] [CrossRef]
  47. Isbister, G.K. Snakebite doesn’t cause disseminated intravascular coagulation: Coagulopathy and thrombotic microangiopathy in snake envenoming. Semin. Thromb. Hemost. 2010, 36, 444–451. [Google Scholar] [CrossRef] [PubMed]
  48. Resiere, D.; Mégarbane, B.; Valentino, R.; Mehdaoui, H.; Thomas, L. Bothrops lanceolatus bites: Guidelines for severity assessment and emergent management. Toxins 2010, 2, 163–173. [Google Scholar] [CrossRef] [Green Version]
  49. Currier, R.B.; Harrison, R.A.; Rowley, P.D.; Laing, G.D.; Wagstaff, S.C. Intra-specific variation in venom of the African Puff Adder (Bitis arietans): Differential expression and activity of snake venom metalloproteinases (SVMPs). Toxicon 2010, 55, 864–873. [Google Scholar] [CrossRef]
  50. Oliveira, I.S.d.; Cardoso, I.A.; Bordon, K.d.C.F.; Carone, S.E.I.; Boldrini-França, J.; Pucca, M.B.; Zoccal, K.F.; Faccioli, L.H.; Sampaio, S.V.; Rosa, J.C.; et al. Global proteomic and functional analysis of Crotalus durissus collilineatus individual venom variation and its impact on envenoming. J. Proteom. 2019, 191, 153–165. [Google Scholar] [CrossRef]
  51. Tasoulis, T.; Silva, A.; Veerati, P.C.; Baker, M.; Hodgson, W.C.; Dunstan, N.; Isbister, G.K. Intra-specific venom variation in the Australian coastal taipan Oxyuranus scutellatus. Toxins 2020, 12, 485. [Google Scholar] [CrossRef]
  52. Estevao-Costa, M.I.; Gontijo, S.S.; Correia, B.L.; Yarleque, A.; Vivas-Ruiz, D.; Rodrigues, E.; Chávez-Olortegui, C.; Oliveira, L.S.; Sanchez, E.F. Neutralization of toxicological activities of medically-relevant Bothrops snake venoms and relevant toxins by two polyvalent bothropic antivenoms produced in Peru and Brazil. Toxicon 2016, 122, 67–77. [Google Scholar] [CrossRef]
  53. Kohlhoff, M.; Borges, M.H.; Yarleque, A.; Cabezas, C.; Richardson, M.; Sanchez, E.F. Exploring the proteomes of the venoms of the Peruvian pit vipers Bothrops atrox, B. barnetti and B. pictus. J. Proteom. 2012, 75, 2181–2195. [Google Scholar] [CrossRef]
  54. Bogarín, G.; Romero, M.; Rojas, G.; Lutsch, C.; Casadamont, M.; Lang, J.; Otero, R.; Gutiérrez, J.M.a. Neutralization, by a monospecific Bothrops lanceolatus antivenom, of toxic activities induced by homologous and heterologous Bothírops snake venoms. Toxicon 1999, 37, 551–557. [Google Scholar] [CrossRef]
  55. Gutiérrez, J.M.; Sanz, L.; Escolano, J.; Fernández, J.; Lomonte, B.; Angulo, Y.; Rucavado, A.; Warrell, D.A.; Calvete, J.J. Snake venomics of the Lesser Antillean pit vipers Bothrops caribbaeus and Bothrops lanceolatus: Correlation with toxicological activities and immunoreactivity of a heterologous antivenom. J. Proteome Res. 2008, 7, 4396–4408. [Google Scholar] [CrossRef] [PubMed]
  56. Niewiarowski, S.; Kirby, E.P.; Brudzynski, T.M.; Stocker, K. Thrombocytin, a serine protease from Bothrops atrox venom. 2. Interaction with platelets and plasma-clotting factors. Biochemistry 1979, 18, 3570–3577. [Google Scholar] [CrossRef] [PubMed]
  57. Rucavado, A.; Soto, M.; Escalante, T.; Loría, G.D.; Arni, R.; Gutiérrez, J.M. Thrombocytopenia and platelet hypoaggregation induced by Bothrops asper snake venom. Thromb. Haemost. 2005, 94, 123–131. [Google Scholar] [CrossRef] [PubMed]
  58. Rucavado, A.; Soto, M.; Kamiguti, A.S.; Theakston, D.G.; Fox, J.W.; Escalante, T.; Gutiérrez, J.M. Characterization of aspercetin, a platelet aggregating component from the venom of the snake Bothrops asper which induces thrombocytopenia and potentiates metalloproteinase-induced hemorrhage. Thromb. Haemost. 2001, 85, 710–715. [Google Scholar] [CrossRef] [PubMed]
  59. Gutiérrez, J.; Ownby, C.L.; Odell, G.V. Isolation of a myotoxin from Bothrops asper venom: Partial characterization and action on skeletal muscle. Toxicon 1984, 22, 115–128. [Google Scholar] [CrossRef]
  60. Gutiérrez, J.M.; Lomonte, B. Phospholipase A2 myotoxins from Bothrops snake venoms. Toxicon 1995, 33, 1405–1424. [Google Scholar] [CrossRef]
  61. Lister, C.; Arbuckle, K.; Jackson, T.N.W.; Debono, J.; Zdenek, C.N.; Dashevsky, D.; Dunstan, N.; Allen, L.; Hay, C.; Bush, B.; et al. Catch a tiger snake by its tail: Differential toxicity, co-factor dependence and antivenom efficacy in a procoagulant clade of Australian venomous snakes. Comp. Biochem. Physiol. C Toxicol. Pharm. 2017, 202, 39–54. [Google Scholar] [CrossRef] [Green Version]
  62. O’Leary, M.A.; Kornhauser, R.S.; Hodgson, W.C.; Isbister, G.K. An examination of the activity of expired and mistreated commercial Australian antivenoms. Trans. R. Soc. Trop. Med. Hyg. 2009, 103, 937–942. [Google Scholar] [CrossRef]
  63. Al-Abdulla, I.; Casewell, N.R.; Landon, J. Long-term physicochemical and immunological stability of a liquid formulated intact ovine immunoglobulin-based antivenom. Toxicon 2013, 64, 38–42. [Google Scholar] [CrossRef]
  64. Da Galizio, N.C.; Serino-Silva, C.; Stuginski, D.R.; Abreu, P.A.E.; Sant’Anna, S.S.; Grego, K.F.; Tashima, A.K.; Tanaka-Azevedo, A.M.; de Morais-Zani, K. Compositional and functional investigation of individual and pooled venoms from long-term captive and recently wild-caught Bothrops jararaca snakes. J. Proteom. 2018, 186, 56–70. [Google Scholar] [CrossRef]
  65. Castellone, D. How to deliver quality results in the coagulation laboratory: Commonly asked questions. Lab. Med. 2004, 35, 208–213. [Google Scholar] [CrossRef]
  66. Bourke, L.A.; Youngman, N.J.; Zdenek, C.N.; op den Brouw, B.; Violette, A.; Fourmy, R.; Fry, B.G. Trimeresurus albolabris snakebite treatment implications arising from ontogenetic venom comparisons of anticoagulant function, and antivenom efficacy. Toxicol. Lett. 2020, 327, 2–8. [Google Scholar] [CrossRef] [PubMed]
  67. Debono, J.; Bos, M.H.A.; Frank, N.; Fry, B. Clinical implications of differential antivenom efficacy in neutralising coagulotoxicity produced by venoms from species within the arboreal viperid snake genus Trimeresurus. Toxicol. Lett. 2019, 316, 35–48. [Google Scholar] [CrossRef] [PubMed]
  68. Zdenek, C.N.; Hay, C.; Arbuckle, K.; Jackson, T.N.W.; Bos, M.H.A.; op den Brouw, B.; Debono, J.; Allen, L.; Dunstan, N.; Morley, T.; et al. Coagulotoxic effects by brown snake (Pseudonaja) and taipan (Oxyuranus) venoms, and the efficacy of a new antivenom. Toxicol. In Vitro 2019, 58, 97–109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Rojas, G.; Jiménez, J.; Gutiérrez, J. Caprylic acid fractionation of hyperimmune horse plasma: Description of a simple procedure for antivenom production. Toxicon 1994, 32, 351–363. [Google Scholar] [CrossRef]
Toxins 13 00078 g001 550
Figure 1. (A) Human plasma clotting times (y-axis) induced by Bothrops asper venom (red line) and venom + antivenom treatments (line colours indicated by the legend) over eight different venom concentrations (μg/mL: 20, 10, 4, 1.66, 0.66, 0.25, 0.125, and 0.05; x-axis). The percentage of antivenom in the final cuvette is 0.25%. The graphs’ x-axis is displayed in logarithmic view. Each of the eight data points per venom curve is represented by dots (mean value, n = 3), with standard deviation error bars. Some error bars are smaller than the data-point symbol. Control values (seconds ± standard deviation) were negative control (spontaneous clotting time of plasma, n = 4) = 459.90 ± 67.13, positive control (kaolin, n = 4) = 63.08 ± 0.67, and antivenom controls (Antivipmyn (n = 4) = 485.63 ± 54.41, PoliVal-ICP (n = 4) = 483.15 ± 27.93, Antivipmyn-Tri (n = 4) = 486.45 ± 20.97, Bothrofav (n = 4) = 521.45 ± 9.32, and SAB (n = 4) = 467.08 ± 30.79). See Section 4.3 and Section 4.4 for control methodologies. (B) X-fold shift in plasma clotting curves (y-axis) of Bothrops asper venom incubated with five different antivenoms (bar colours indicated by the legend). X-fold shift values were calculated for each venom by dividing the area under the curve (AUC) for the venom incubated with antivenom clotting time curve by the AUC for the venom clotting time curve, and subtracting the total by 1. Calculated values represent antivenom neutralisation: a value of 0 indicates no neutralisation, while a value of above 0 indicates neutralisation. Values are mean (n = 3) ± standard deviation. Letters above bars show statistical significance (p < 0.05). Different letters indicate a significant difference between values, while the same letter across multiple bars indicates no significance. The statistical test used was Tukey’s multiple comparisons test, following significant ordinary one-way analysis of variance (ANOVA) results. SAB = Soro Antibotrópico.
Figure 1. (A) Human plasma clotting times (y-axis) induced by Bothrops asper venom (red line) and venom + antivenom treatments (line colours indicated by the legend) over eight different venom concentrations (μg/mL: 20, 10, 4, 1.66, 0.66, 0.25, 0.125, and 0.05; x-axis). The percentage of antivenom in the final cuvette is 0.25%. The graphs’ x-axis is displayed in logarithmic view. Each of the eight data points per venom curve is represented by dots (mean value, n = 3), with standard deviation error bars. Some error bars are smaller than the data-point symbol. Control values (seconds ± standard deviation) were negative control (spontaneous clotting time of plasma, n = 4) = 459.90 ± 67.13, positive control (kaolin, n = 4) = 63.08 ± 0.67, and antivenom controls (Antivipmyn (n = 4) = 485.63 ± 54.41, PoliVal-ICP (n = 4) = 483.15 ± 27.93, Antivipmyn-Tri (n = 4) = 486.45 ± 20.97, Bothrofav (n = 4) = 521.45 ± 9.32, and SAB (n = 4) = 467.08 ± 30.79). See Section 4.3 and Section 4.4 for control methodologies. (B) X-fold shift in plasma clotting curves (y-axis) of Bothrops asper venom incubated with five different antivenoms (bar colours indicated by the legend). X-fold shift values were calculated for each venom by dividing the area under the curve (AUC) for the venom incubated with antivenom clotting time curve by the AUC for the venom clotting time curve, and subtracting the total by 1. Calculated values represent antivenom neutralisation: a value of 0 indicates no neutralisation, while a value of above 0 indicates neutralisation. Values are mean (n = 3) ± standard deviation. Letters above bars show statistical significance (p < 0.05). Different letters indicate a significant difference between values, while the same letter across multiple bars indicates no significance. The statistical test used was Tukey’s multiple comparisons test, following significant ordinary one-way analysis of variance (ANOVA) results. SAB = Soro Antibotrópico.
Toxins 13 00078 g001
Toxins 13 00078 g002 550
Figure 2. (A) Human plasma clotting times (y-axis) induced by Bothrops atrox venom (red line) and venom + antivenom treatments (line colours indicated by the legend) over eight different venom concentrations (μg/mL: 20, 10, 4, 1.66, 0.66, 0.25, 0.125, and 0.05; x-axis). The percentage of antivenom in the final cuvette is 0.25%. The graphs’ x-axis is displayed in logarithmic view. Each of the eight data points per venom curve is represented by dots (mean value, n = 3), with standard deviation error bars. Some error bars are smaller than the data-point symbol. Control values (seconds ± standard deviation) were negative control (spontaneous clotting time of plasma, n = 4) = 459.90 ± 67.13, positive control (kaolin, n = 4) = 63.08 ± 0.67, and antivenom controls (Antivipmyn (n = 4) = 485.63 ± 54.41, PoliVal-ICP (n = 4) = 483.15 ± 27.93, Antivipmyn-Tri (n = 4) = 486.45 ± 20.97, Bothrofav (n = 4) = 521.45 ± 9.32, and SAB (n = 4) = 467.08 ± 30.79). See Section 4.3 and Section 4.4 for control methodologies. (B) X-fold shift in plasma clotting curves (y-axis) of Bothrops atrox venom incubated with five different antivenoms (bar colours indicated by the legend). X-fold shift values were calculated for each venom by dividing the area under the curve (AUC) for the venom incubated with antivenom clotting time curve by the AUC for the venom clotting time curve, and subtracting the total by 1. Calculated values represent antivenom neutralisation: a value of 0 indicates no neutralisation, while a value of above 0 indicates neutralisation. Values are mean (n = 3) ± standard deviation. Letters above bars show statistical significance (p < 0.05). Different letters indicate a significant difference between values, while the same letter across multiple bars indicates no significance. The statistical test used was Tukey’s multiple comparisons test, following significant ordinary one-way analysis of variance (ANOVA) results. SAB = Soro Antibotrópico.
Figure 2. (A) Human plasma clotting times (y-axis) induced by Bothrops atrox venom (red line) and venom + antivenom treatments (line colours indicated by the legend) over eight different venom concentrations (μg/mL: 20, 10, 4, 1.66, 0.66, 0.25, 0.125, and 0.05; x-axis). The percentage of antivenom in the final cuvette is 0.25%. The graphs’ x-axis is displayed in logarithmic view. Each of the eight data points per venom curve is represented by dots (mean value, n = 3), with standard deviation error bars. Some error bars are smaller than the data-point symbol. Control values (seconds ± standard deviation) were negative control (spontaneous clotting time of plasma, n = 4) = 459.90 ± 67.13, positive control (kaolin, n = 4) = 63.08 ± 0.67, and antivenom controls (Antivipmyn (n = 4) = 485.63 ± 54.41, PoliVal-ICP (n = 4) = 483.15 ± 27.93, Antivipmyn-Tri (n = 4) = 486.45 ± 20.97, Bothrofav (n = 4) = 521.45 ± 9.32, and SAB (n = 4) = 467.08 ± 30.79). See Section 4.3 and Section 4.4 for control methodologies. (B) X-fold shift in plasma clotting curves (y-axis) of Bothrops atrox venom incubated with five different antivenoms (bar colours indicated by the legend). X-fold shift values were calculated for each venom by dividing the area under the curve (AUC) for the venom incubated with antivenom clotting time curve by the AUC for the venom clotting time curve, and subtracting the total by 1. Calculated values represent antivenom neutralisation: a value of 0 indicates no neutralisation, while a value of above 0 indicates neutralisation. Values are mean (n = 3) ± standard deviation. Letters above bars show statistical significance (p < 0.05). Different letters indicate a significant difference between values, while the same letter across multiple bars indicates no significance. The statistical test used was Tukey’s multiple comparisons test, following significant ordinary one-way analysis of variance (ANOVA) results. SAB = Soro Antibotrópico.
Toxins 13 00078 g002
Toxins 13 00078 g003 550
Figure 3. Each bar graph represents the X-fold shift value of each venom when incubated with antivenom (antivenom indicated by the graphs’ title). Note SAB = Soro Antibotrópico. X-fold shift values were calculated for each venom by dividing the area under the curve (AUC) for the venom incubated with antivenom clotting time curve by the AUC for the venom clotting time curve, and subtracting the total by 1. Calculated values represent antivenom neutralisation: a value of 0 indicates no neutralisation, while a value of above 0 indicates neutralisation. Each bar is a mean (n = 3) ± standard deviation. Red bars indicate venom samples are from Bothrops asper, while blue bars indicate venom samples are from Bothrops atrox. Venom sample information is also displayed below each bar. The percentage of antivenom in the final cuvette is 0.25%.
Figure 3. Each bar graph represents the X-fold shift value of each venom when incubated with antivenom (antivenom indicated by the graphs’ title). Note SAB = Soro Antibotrópico. X-fold shift values were calculated for each venom by dividing the area under the curve (AUC) for the venom incubated with antivenom clotting time curve by the AUC for the venom clotting time curve, and subtracting the total by 1. Calculated values represent antivenom neutralisation: a value of 0 indicates no neutralisation, while a value of above 0 indicates neutralisation. Each bar is a mean (n = 3) ± standard deviation. Red bars indicate venom samples are from Bothrops asper, while blue bars indicate venom samples are from Bothrops atrox. Venom sample information is also displayed below each bar. The percentage of antivenom in the final cuvette is 0.25%.
Toxins 13 00078 g003
Table
Table 1. Venom-induced (20 μg/mL) clotting time (seconds) and area under the curve (AUC) values for Bothrops venom samples. All values are mean (n = 3) ± standard deviation.
Table 1. Venom-induced (20 μg/mL) clotting time (seconds) and area under the curve (AUC) values for Bothrops venom samples. All values are mean (n = 3) ± standard deviation.
Venom20 μg/mL Clotting TimeAUC Value
B. asper (Costa Rica)17.73 ± 2.17481.70 ± 18.01
B. asper (Ecuador)13.50 ± 0.85384.00 ± 5.86
B. asper (Mérida, Yucatán, Mexico)14.27 ± 1.12377.33 ± 9.15
B. asper (San Andres, Tuxtla, Veracruz, Mexico)47.20 ± 2.041327.33 ± 26.03
B. atrox (French Guiana)26.80 ± 0.96669.70 ± 15.85
B. atrox (Alto Marañon, Peru)28.30 ± 0.66741.73 ± 5.47
Table
Table 2. p-values of Tukey’s multiple comparisons tests for an ordinary one-way analysis of variance (ANOVA; F (5, 12) = 1644.00, p < 0.0001) comparing area under the curve (AUC) values of venom samples. A significant result is p < 0.05. X indicates no valid comparison.
Table 2. p-values of Tukey’s multiple comparisons tests for an ordinary one-way analysis of variance (ANOVA; F (5, 12) = 1644.00, p < 0.0001) comparing area under the curve (AUC) values of venom samples. A significant result is p < 0.05. X indicates no valid comparison.
B. asper (Costa Rica)B. asper (Ecuador)B. asper (Mérida, Yucatán, Mexico)B. asper (San Andres, Tuxtla, Veracruz, Mexico)B. atrox (French
Guiana)		B. atrox (Alto Marañon, Peru)
B. asper (Costa Rica)X<0.0001<0.0001<0.0001<0.0001<0.0001
B. asper (Ecuador)XX0.9936<0.0001<0.0001<0.0001
B. asper (Mérida, Yucatán, Mexico)XXX<0.0001<0.0001<0.0001
B. asper (San Andres, Tuxtla, Veracruz, Mexico)XXXX<0.0001<0.0001
B. atrox (French Guiana)XXXXX0.0010
B. atrox (Alto Marañon, Peru)XXXXXX
Table
Table 3. p-values of Tukey’s multiple comparisons tests comparing the X-fold shift values of different antivenoms for all Bothrops asper samples. Tukey’s multiple comparisons tests were performed on each B. asper locality X-fold shift dataset following significant ordinary one-way analysis of variance (ANOVA) results: B. asper (Costa Rica; F (4, 10) = 61.52, p < 0.0001), B. asper (Ecuador; F (4, 10) = 460.40, p < 0.0001), B. asper (Mérida, Yucatán, Mexico; F (4, 10) = 25.05, p < 0.0001) and B. asper (San Andres, Tuxtla, Veracruz, Mexico; F (4, 10) = 209.90, p < 0.0001). A significant result is p < 0.05. SAB = Soro Antibotrópico.
Table 3. p-values of Tukey’s multiple comparisons tests comparing the X-fold shift values of different antivenoms for all Bothrops asper samples. Tukey’s multiple comparisons tests were performed on each B. asper locality X-fold shift dataset following significant ordinary one-way analysis of variance (ANOVA) results: B. asper (Costa Rica; F (4, 10) = 61.52, p < 0.0001), B. asper (Ecuador; F (4, 10) = 460.40, p < 0.0001), B. asper (Mérida, Yucatán, Mexico; F (4, 10) = 25.05, p < 0.0001) and B. asper (San Andres, Tuxtla, Veracruz, Mexico; F (4, 10) = 209.90, p < 0.0001). A significant result is p < 0.05. SAB = Soro Antibotrópico.
B. asper
(Costa Rica)		B. asper
(Ecuador)		B. asper (Mérida, Yucatán,
Mexico)		B. asper
(San Andres, Tuxtla, Veracruz,
Mexico)
Antivipmyn vs. PoliVal-ICP<0.0001<0.0001<0.0001<0.0001
Antivipmyn vs. Antivipmyn-Tri>0.99990.99960.99840.3096
Antivipmyn vs. Bothrofav0.97100.96670.9993<0.0001
Antivipmyn vs. SAB0.94370.01690.9525<0.0001
PoliVal-ICP vs. Antivipmyn-Tri<0.0001<0.0001<0.0001<0.0001
PoliVal-ICP vs. Bothrofav<0.0001<0.00010.0001<0.0001
PoliVal-ICP vs. SAB<0.0001<0.00010.0002<0.0001
Antivipmyn-Tri vs. Bothrofav0.98480.99160.9851<0.0001
Antivipmyn-Tri vs. SAB0.91460.01250.8598<0.0001
Bothrofav vs. SAB0.67590.00660.98820.5859
Table
Table 4. p-values of Tukey’s multiple comparisons tests comparing the X-fold shift values of different antivenoms for both Bothrops atrox samples. Tukey’s multiple comparisons tests were performed on each B. atrox locality X-fold shift dataset, following significant ordinary one-way analysis of variance (ANOVA) results: B. atrox (French Guiana; F (4, 10) = 993.80, p < 0.0001), and B. atrox (Alto Marañon, Peru; F (4, 10) = 1226.00, p < 0.0001).
Table 4. p-values of Tukey’s multiple comparisons tests comparing the X-fold shift values of different antivenoms for both Bothrops atrox samples. Tukey’s multiple comparisons tests were performed on each B. atrox locality X-fold shift dataset, following significant ordinary one-way analysis of variance (ANOVA) results: B. atrox (French Guiana; F (4, 10) = 993.80, p < 0.0001), and B. atrox (Alto Marañon, Peru; F (4, 10) = 1226.00, p < 0.0001).
B. atrox (French Guiana)B. atrox (Alto Marañon, Peru)
Antivipmyn vs. PoliVal-ICP0.0408<0.0001
Antivipmyn vs. Antivipmyn-Tri0.93180.0266
Antivipmyn vs. Bothrofav<0.0001<0.0001
Antivipmyn vs. SAB<0.0001<0.0001
PoliVal-ICP vs. Antivipmyn-Tri0.13440.0013
PoliVal-ICP vs. Bothrofav<0.0001<0.0001
PoliVal-ICP vs. SAB<0.0001<0.0001
Antivipmyn-Tri vs. Bothrofav<0.0001<0.0001
Antivipmyn-Tri vs. SAB<0.0001<0.0001
Bothrofav vs. SAB<0.0001<0.0001
Table
Table 5. p-values of Tukey’s multiple comparisons tests comparing X-fold shift values between venoms within an antivenom. Tukey’s multiple comparisons tests were performed on the PoliVal-CP, SAB (Soro Antibotrópico) and BothroFav X-fold shift datasets, following significant ordinary one-way analysis of variance (ANOVA) results: PoliVal-ICP (F (5, 12) = 71.58, p < 0.0001), SAB (F (5, 12) = 643.40, p < 0.0001), and Bothrofav (F (5, 12) = 389.90, p < 0.0001). A significant result is p < 0.05. Note that not all Tukey’s multiple comparisons are shown in this table.
Table 5. p-values of Tukey’s multiple comparisons tests comparing X-fold shift values between venoms within an antivenom. Tukey’s multiple comparisons tests were performed on the PoliVal-CP, SAB (Soro Antibotrópico) and BothroFav X-fold shift datasets, following significant ordinary one-way analysis of variance (ANOVA) results: PoliVal-ICP (F (5, 12) = 71.58, p < 0.0001), SAB (F (5, 12) = 643.40, p < 0.0001), and Bothrofav (F (5, 12) = 389.90, p < 0.0001). A significant result is p < 0.05. Note that not all Tukey’s multiple comparisons are shown in this table.
SpeciesAntivenomB. asper
(Costa Rica)		B. asper (Mérida,
Yucatán, Mexico)		B. asper (San Andres, Tuxtla,
Veracruz, Mexico)		B. asper
(Ecuador)
B. atrox (French
Guiana)		PoliVal-ICP<0.0001<0.0001<0.0001<0.0001
SAB<0.0001<0.00010.0002<0.0001
Bothrofav<0.0001<0.00010.8651<0.0001
B. atrox (Alto
Marañon,
Peru)		PoliVal-ICP<0.00010.16370.00030.0023
SAB<0.0001<0.0001<0.0001<0.0001
Bothrofav<0.0001<0.00010.0065<0.0001
Table
Table 6. p-values of Tukey’s multiple comparisons tests comparing X-fold shift values between Central American/Mexican and South American venoms within the antivenom Antivipmyn. Tukey’s multiple comparisons tests were performed on the Antivipmyn X-fold shift dataset, following significant ordinary one-way analysis of variance (ANOVA) results: Antivipmyn (F (5, 12) = 50.92, p < 0.0001). A significant result is p < 0.05. Note that not all Tukey’s multiple comparisons are shown in this table.
Table 6. p-values of Tukey’s multiple comparisons tests comparing X-fold shift values between Central American/Mexican and South American venoms within the antivenom Antivipmyn. Tukey’s multiple comparisons tests were performed on the Antivipmyn X-fold shift dataset, following significant ordinary one-way analysis of variance (ANOVA) results: Antivipmyn (F (5, 12) = 50.92, p < 0.0001). A significant result is p < 0.05. Note that not all Tukey’s multiple comparisons are shown in this table.
Central America and Mexico
B. asper (Costa Rica)B. asper (Mérida, Yucatán,
Mexico)		B. asper (San Andres, Tuxtla, Veracruz, Mexico)
South
America		B. asper (Ecuador)<0.00010.0044<0.0001
B. atrox (French Guiana)<0.00010.0213<0.0001
B. atrox (Alto Marañon, Peru)<0.00010.79040.0019
Table
Table 7. List of species, localities, abbreviations, and additional information if known (age, sex, habitat, and whether venoms were from wild-caught or captive-born snakes, single individuals or multiple individuals (pooled)). Captive-born samples originate from snakes from the specified locality. N.A. = nonapplicable.
Table 7. List of species, localities, abbreviations, and additional information if known (age, sex, habitat, and whether venoms were from wild-caught or captive-born snakes, single individuals or multiple individuals (pooled)). Captive-born samples originate from snakes from the specified locality. N.A. = nonapplicable.
Species LocalityAbbreviationAdditional Information
Bothrops asperCosta RicaN.A.Pooled from adults collected in the Pacific region (n = 40)
Bothrops asperEcuadorN.A.Pooled from captive-born adult males (n = 2)
Bothrops asperMérida, Yucatán, MexicoYucatán, MexicoAdult (individual)
Bothrops asperSan Andres, Tuxtla, Veracruz, MexicoS.A.T, MexicoAdult (individual)
Bothrops atroxFrench GuianaN.A.Pooled from captive-born and wild adults (male + females, n = 66)
Bothrops atroxAlto Marañon, PeruA.M., PeruPooled from wild adults caught in the Amazon rainforest and kept in captivity (n values not supplied)
Table
Table 8. Antivenoms used and their relevant information. Information was obtained from package inserts, scientific literature (cited in the table), and/or correspondence with researchers. Species names may have changed since production of the antivenom.
Table 8. Antivenoms used and their relevant information. Information was obtained from package inserts, scientific literature (cited in the table), and/or correspondence with researchers. Species names may have changed since production of the antivenom.
AntivenomLot # and ExpiryImmunising Mixture *PreparationNeutralising
Potency (Specified by Manufacturer)
Antivipmyn®, Instituto Bioclon, MexicoLot: B-6F-16; Exp: October 2010Bothrops asper (Mexico),
Crotalus simus simus (Mexico)		Polyspecific equine F(ab’)2 preparation1 vial neutralises >780 LD50 of Bothrops sp. and >790 of Crotalus sp.
Antivipmyn-Tri®, Instituto
Bioclon, Mexico		Lot: B-4F-13; Exp: 28 June 2009Bothrops asper (Colombia), Crotalus simus simus (Mexico), and Lachesis muta (Not stated)Polyspecific equine F(ab’)2 preparationUnknown
Soro Antibotrópico (SAB)®, Instituto Butantan, BrazilLot: 1305077; Exp: May 2016Bothrops jararaca (Brazil) (50%), Bothrops jararacussu (Brazil) (12.5%), Bothrops neuwiedi (Brazil) (12.5%), Bothrops alternatus (Brazil) (12.5%), Bothrops moojeni (Brazil) (12.5%)Polyspecific equine IgG F(ab′)2 preparation1 mL neutralises 5 mg venom of Botrops jararaca (Brazil)
Central American polyspecific antivenom (PoliVal-ICP)®, Instituto Clodomiro Picado, Universidad de Costa RicaLot: 5720416; Exp: April 2021Bothrops asper (Costa Rica), Crotalus simus (Costa Rica), and Lachesis stenophrys (Costa Rica) [69]Polyspecific equine whole IgG—purified by caprylic acid fractionation [69]1ml neutralises 3 mg venom of Bothrops asper (Costa Rica)
Bothrofav®, MicroPharm, United KingdomLot P4AP61V; Exp: October 2020Bothrops lanceolatusMonospecific equine F(ab′)2 preparationUnknown
* The species of snakes whose venom was used for antivenom production.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bourke, L.A.; Zdenek, C.N.; Neri-Castro, E.; Bénard-Valle, M.; Alagón, A.; Gutiérrez, J.M.; Sanchez, E.F.; Aldridge, M.; Fry, B.G. Pan-American Lancehead Pit-Vipers: Coagulotoxic Venom Effects and Antivenom Neutralisation of Bothrops asper and B. atrox Geographical Variants. Toxins 2021, 13, 78. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13020078

AMA Style

Bourke LA, Zdenek CN, Neri-Castro E, Bénard-Valle M, Alagón A, Gutiérrez JM, Sanchez EF, Aldridge M, Fry BG. Pan-American Lancehead Pit-Vipers: Coagulotoxic Venom Effects and Antivenom Neutralisation of Bothrops asper and B. atrox Geographical Variants. Toxins. 2021; 13(2):78. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13020078

Chicago/Turabian Style

Bourke, Lachlan A., Christina N. Zdenek, Edgar Neri-Castro, Melisa Bénard-Valle, Alejandro Alagón, José María Gutiérrez, Eladio F. Sanchez, Matt Aldridge, and Bryan G. Fry. 2021. "Pan-American Lancehead Pit-Vipers: Coagulotoxic Venom Effects and Antivenom Neutralisation of Bothrops asper and B. atrox Geographical Variants" Toxins 13, no. 2: 78. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13020078

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop