Venom proteomics and antivenom neutralization for the Chinese eastern Russell’s viper, Daboia siamensis from Guangxi and Taiwan

The eastern Russell’s viper (Daboia siamensis) causes primarily hemotoxic envenomation. Applying shotgun proteomic approach, the present study unveiled the protein complexity and geographical variation of eastern D. siamensis venoms originated from Guangxi and Taiwan. The snake venoms from the two geographical locales shared comparable expression of major proteins notwithstanding variability in their toxin proteoforms. More than 90% of total venom proteins belong to the toxin families of Kunitz-type serine protease inhibitor, phospholipase A2, C-type lectin/lectin-like protein, serine protease and metalloproteinase. Daboiasiamensis Monovalent Antivenom produced in Taiwan (DsMAV-Taiwan) was immunoreactive toward the Guangxi D. siamensis venom, and effectively neutralized the venom lethality at a potency of 1.41 mg venom per ml antivenom. This was corroborated by the antivenom effective neutralization against the venom procoagulant (ED = 0.044 ± 0.002 µl, 2.03 ± 0.12 mg/ml) and hemorrhagic (ED50 = 0.871 ± 0.159 µl, 7.85 ± 3.70 mg/ml) effects. The hetero-specific Chinese pit viper antivenoms i.e. Deinagkistrodon acutus Monovalent Antivenom and Gloydius brevicaudus Monovalent Antivenom showed negligible immunoreactivity and poor neutralization against the Guangxi D. siamensis venom. The findings suggest the need for improving treatment of D. siamensis envenomation in the region through the production and the use of appropriate antivenom.

Immunoreactivity of antivenoms to D. siamensis venoms. The four antivenoms (DsMAV-Taiwan, DsMAV-Thai, DaMAV and GbMAV) and an additional combination of DaMAV and GbMAV in a ratio of 1:1 were tested for their immunological binding activity toward Ds-Guangxi and Ds-Taiwan venoms (Fig. 2). The highest immunoreactivity (reflecting antigen binding activity) was observed in the reaction of DsMAV-Taiwan with Ds-Taiwan venom (Fig. 2). The Ds-Guangxi venom showed a lower immunoreactivity (approximately 50% lower) when reacting with the homologous antivenom from Taiwan. Both venoms also showed immunoreactivity for DsMAV-Thai but the relative magnitude of reactivity was apparently lower, found in the range of 15% (Ds-Guangxi) to 30% (Ds-Taiwan). The heterologous antivenoms (DaMAV, GbMAV and a combination of DaMAV:GbMAV in a ratio of 1:1) were generally low in immunoreactivity (<10%) to both Ds-Guangxi and Ds-Taiwan venoms.

Procoagulant activity of D. siamensis venoms and antivenom neutralization. Both D. siamensis
venoms exhibited potent procoagulant effect with minimal coagulation dose (MCD) of 0.23 ± 0.06 µg/ml and 0.15 ± 0.04 µg/ml for Ds-Guangxi and Ds-Taiwan venoms, respectively (Table 3). DsMAV-Taiwan neutralized the procoagulant effect of Ds-Guangxi and Ds-Taiwan venoms in a dose-dependent manner (Fig. 3A), and the efficacy values of neutralization (defined as ED, effective dose) were comparable for both venoms (ED = ~1.4 to 2.0 mg venom neutralized per millilitre of antivenom, Table 3). The neutralization by GbMAV was extremely poor, at least 30-fold less effective compared with DsMAV-Taiwan. On the other hand, DaMAV was totally ineffective in neutralizing the procoagulant effect of Ds-Guangxi and Ds-Taiwan venoms. For comparison, the effective doses of the antivenoms were normalized by the respective protein concentrations ( Table 3). The normalized effective dose (n-ED) for procoagulant effect of DsMAV-Taiwan (in mg/g, milligram of venom neutralized per gram of antivenom protein) was at least 250-fold higher than those of GbMAV and DaMAV.

Hemorrhagic activity of Daboia siamensis venoms and antivenom neutralization. The Ds-Guangxi
venom showed a lower minimal hemorrhagic dose (MHD) of 3.42 ± 0.12 µg/mouse compared with the Ds-Taiwan venom (MHD = 8.21 ± 0.31 µg/mouse). DsMAV-Taiwan neutralized the hemorrhagic effect of the venoms dose-dependently (Fig. 3B), but the neutralization was much more effective against the hemorrhagic effect of Ds-Taiwan venom than that of Ds-Guangxi venom based on the normalized median effective doses (n-ED 50 ) ( Table 4). In comparison with DsMAV-Taiwan, GbMAV and DaMAV were extremely weak in neutralizing the hemorrhagic effect induced by the D. siamensis venoms. DsMAV-Taiwan was at least 70-300 folds more effective in neutralizing the hemorrhagic effect.
Lethality of D. siamensis venoms and in vivo neutralization in mice. When administrated intravenously, Ds-Guangxi and Ds-Taiwan venoms had median lethal doses (LD 50 ) of 0.18 µg/g and 0.09 µg/g,  Taiwan (island); KSPI, Kunitz-type serine protease inhibitor; PLA 2 , phospholipase A 2 ; Snaclec, snake venom C-type lectin/lectin-like protein; SVSP, snake venom serine protease; SVMP, snake venom metalloproteinase; LAAO, L-amino acid oxidase; svVEGF, snake venom vascular endothelial growth factor; svNGF, snake venom nerve growth factor; 5′NUC, 5′-nucleotidase; CRiSP, cysteinerich secretory protein; PDE, phosphodiesterase. Note: SDS-PAGE of Ds-Guangxi and Ds-Taiwan venoms were cropped from the same gel for display purpose. The full-length gel is provided in the Supplementary Information File S1.  (Table 5). DsMAV-Taiwan showed dose-dependent neutralization effects against the lethality of the venoms (constituting 5 LD 50 ) (Fig. 3C). As the two venoms were different in their LD 50 and the total amount of venom injected, the lethality neutralization was expressed in terms of neutralization potency (P) ( Table 5). The neutralization of DsMAV-Taiwan was slightly more potent on Ds-Taiwan venom (P = 1.62 mg/ml, or n-P = 83.9 mg/g) than on Ds-Guangxi venom (P = 1.41 mg/ml, or normalized potency, n-P = 73.1 mg/g). On the other hand, GbMAV had a much lower potency of cross-neutralization against the venom lethality (P = 0.17 mg/ml or n-P = 1.0 mg/g), at least 70-fold lesser when comparing with the n-P values of DsMAV-Taiwan. DaMAV was totally ineffective to cross-neutralize the venom lethality at the maximal dose of antivenom (200 µl) administered intravenously into the mice.

Discussion
Thorpe et al. 2 suggested that the eastern Russell's viper underwent an almost simultaneous rapid divergence 2-3 mybp. The Myanmese or Cambodian specimen is basal in the phylogenetic tree, while the Javan branch is sister to the geographically distant Chinese and Taiwanese branches 2 . This was likely to be at a time of mainland range expansion that enabled rapid overland colonization of the snake into the Taiwan once physically joined to the mainland of China over the Pleistocene 43 . Nonetheless, the clinical manifestations of snakebite envenoming often show marked geographical variations that are broadly unrelated to the phylogeny, and this phenomenon has been elucidated by venom proteomic and transcriptomic studies of a number of snake species 6,10,11,14,44 . Hence, the observed incongruence between the envenoming effects and the primary phylogenetic division of the Russell's viper complex does not allow the prediction of Russell's viper envenoming effects and antivenom responses in areas where these have not been studied directly 2,10,12,14 . In this context, D. siamensis venoms from the southern China and Taiwan were reported to cause hemotoxic envenomation [19][20][21] ; however, knowledge on the venom composition was largely limited to isolated toxins of specimens from either the insular Taiwan 31,32,34,45 or mainland China [46][47][48] . Using a comparative approach, this present work successfully unveiled the proteomes and toxic properties of the venom sourced from two locales (Guangxi and Taiwan). The findings revealed a greater number of protein families in the proteomes of the two D. siamensis venoms than that of Myanmese D. siamensis reported previously where only six protein families were identified (i.e. phospholipase A 2 , PLA 2 ; C-type lectin/lectin-like protein, snaclecs; L-amino acid oxidase. LAAO; snake venom serine protease, SVSP; snake venom metalloproteinase, SVMP and snake venom vascular endothelial growth factor; svVEGF) 37 . The discrepancy observed could be probably due to different techniques used, as the current approach (integrating whole venom in-solution shotgun proteomics with recent database mining) might be a more sensitive venomic tool for studying the complexity  and diversity of venom proteins [49][50][51][52][53] . The major protein families (Kunitz-type serine protease inhibitor or KSPI, PLA 2 , snaclec, SVMP and SVSP) expressed in both Ds-Guangxi and Ds-Taiwan venoms were essentially the same, suggesting that these paralogous proteins from specimens of the two locales were highly conserved with evolutionary significance. The variations in the proteoforms probably implied further molecular adaptation to different ecological niches. The predominance of KSPI in both Ds-Guangxi and Ds-Taiwan venoms is in agreement with the high abundance of this non-enzymatic protein reported in the venom proteomes of the Russell's viper from Pakistan (28.4%) 12 and western India (32.5%) 14 . KSPI however was not reported in the venom proteome of the Myanmese D. siamensis 37 , although more recently this protein has been isolated from the Myanmese and Chinese D. siamensis venoms 29,30 . The KSPI proteoforms detected in the present study showed sequences matched to those reported previously for D. siamensis of Myanmar 30 , China 29 and an unreported locale 54 . In general, KSPI are protease inhibitors with approximately 60-66 amino acids (~7 kDa) and are homologous to the conserved Kunitz motif present in the bovine pancreatic trypsin inhibitor 55 . Besides its serine protease inhibitory activity 29,30 which could be important for venom protein storage in the venom glands, KSPI of Russell's viper venom has also been shown to exhibit anticoagulant effect 56,57 . The pathophysiological role of KSPI in Ds-Guangxi and Ds-Taiwan venoms hence may be coagulopathy-related but the properties of the purified toxin await further investigation.
The high abundance of PLA 2 in both D. siamensis venoms from Guangxi and Taiwan is in line with the finding of high content of PLA 2 (~35%) in the Myanmese D. siamensis venom 31 . PLA 2 is a dominant venom protein family in virtually all reported D. russelii venom proteomes, including those from Pakistan (32.8-63.8%) 12,13 , Western India (32.5%) 14 and Sri Lanka (35%) 15 . However, the activities of the different PLA 2 subtypes are diverse, and the effects can be ranging from no toxicity to high lethality 31,58,59 . The PLA 2 detected in Ds-Guangxi and Ds-Taiwan venoms were apparently of different subtypes, implying that the PLA 2 toxic activities could be diverse. In Ds-Taiwan venom, the presence of PLA 2 RV-4 and RV-7 (in a ratio of 1:1) is consistent with the PLA 2 isolated from the Taiwanese Russell's viper venom reported earlier, where the PLA 2 formed a heterodimeric complex that induced presynaptic neurotoxicity 60 . Putative neurotoxic PLA 2 were also detected in Ds-Guangxi venom (daboiatoxin A, B chain and DsM-b1) 61,62 . However, neurotoxicity is mainly reported in envenomation by D. russelii in Sri Lanka and some parts of southern India 8,26,63,64 ; it is not a commonly observed clinical feature in envenomation by the eastern Russell's vipers in Southeast Asia 26 , Taiwan 21 or China 17,65 . The neurotoxicity induced by D. siamensis PLA 2 in laboratory animals probably reflects the complex interactions between toxins and the neurons of different specificity in animals 60 , where the natural prey such as rodents appear to be more susceptible to the PLA 2 -induced neurotoxicity than human beings are.    [17][18][19][20][21] . The hemotoxicity of D. siamensis venom is collectively caused by a number of toxins. Viperid PLA 2 , including a neutral PLA 2 purified from D. russelii venom, were known to exhibit    66,67 . Other venom proteins detected in substantial amounts in Ds-Guangxi and Ds-Taiwan venoms were mainly hemotoxins that can induce consumptive coagulopathy and hemorrhage, such as snaclec, SVSP and SVMP. Snaclecs (comprising C-type lectins and C-type lectin-like proteins of snake venom) are non-enzymatic toxins that can modulate thrombosis and hemostasis [68][69][70] . Various proteoforms of snaclec were detected in both venom proteomes, including dabocetin, a heterodimer consisting alpha and beta subunits that inhibit ristocetin-induced platelet aggregation 32 . Two RVV-X light chains which are homologous to snaclecs were also detected in the venoms. As part of the RVV-X metalloproteinase, these C-type lectin-like proteins were suggested to play a regulatory role in the calcium-dependant activation of factor X, probably through the recognition of specific sites of the zymogen factor X 68,71 .
Snake venom serine protease (SVSP) is another important protein family of viperid venoms that can cause venom-induced consumptive coagulopathy 38,72,73 . Several SVSP proteoforms were detected in the D. siamensis venom proteomes, one of which is Factor V activating enzyme (RVV-V), a serine protease that specifically activates Factor V (through cleavage at Arg1545-Ser1546 bond) to induce prothrombinase complex in a calcium-dependent manner 74 . The presence of alpha and beta fibrinogenases in the venoms also suggests that fibrinogenolytic activity may contribute to systemic coagulopathy in envenomation. Besides, the SVMP VLAIP-A detected in the Ds-Guangxi and Ds-Taiwan venoms is another hemotoxin known to induce fibrinogenolysis and coagulopathy 75 . More importantly, the Factor X activating enzyme (RVV-X), present in both D. siamensis venoms, is a potent procoagulant enzyme unique to Russell's viper. RVV-X is a heterotrimeric metalloproteinase (93 kDa) composed of a heavy chain (from PIII-SVMP containing metalloproteinase, disintegrin-like and cysteine rich domains, 58 kDa) and two light chains of snaclecs (beta-chain, ~19 kDa and gamma-chain, ~16 kDa) 68,71,76 . In this study, the higher abundance of RVV-X in the Ds-Taiwan venom corroborated the stronger procoagulant effect of the venom on human citrated plasma. On the other hand, the hemorrhagic PIII-SVMP daborhagin-K was detected only in Ds-Guangxi venom; this PIII-SVMP was similar to the potent hemorrhagin SVMP purified from the Myanmese D. siamensis venom 34 . The substantial amount of PIII-SVMP detected in the venom proteomes hence supported the clinical presentation of hemorrhages in D. siamensis envenomation. Nevertheless, the distinctive presence of daborhagin-K in Ds-Guangxi venom correlated with the more potent hemorrhagic effect of the venom (shown in this study), and this venom property might be associated with the prominent local and systemic bleeding reported in Chinese D. siamensis envenoming 20 . Clinically, the potential renal complication (acute kidney injury) of Russell's viper envenoming could be due to the nephrotoxic effect of SVMP and PLA 2 mediated through cytotoxic activity or secondary to renal hypoperfusion in severe bleeding 26,77,78 .
In the current study, LAAO was detected only in Ds-Guangxi venom, consistent with the observation of a more intense protein band around 60 kDa on the SDS−PAGE under reducing conditions. The absence of LAAO in the proteome of Ds-Taiwan venom is puzzling as this enzyme is present in the venoms of most Viperidae and Elapidae 11,79-81 , including the Russell's vipers of South Asia (Sri Lanka, India, Pakistan) [12][13][14][15][16] and Myanmar 37 . In fact, the colors of the two venoms studied in the present study exhibited marked differences: the yellow coloration of Ds-Guangxi venom was likely due to the presence of flavin-containing LAAO, while the lyophilized Ds-Taiwan venom was white in color, implying the lack of this enzyme in the venom. The finding indicated that the amount of LAAO in Ds-Taiwan venom was low or negligible, a characteristic perhaps influenced by ecology and the condition of diet. LAAO may exhibit cytotoxicity and anti-microbial activity to facilitate prey digestion 82 , but its pathophysiological role in D. siamensis envenoming remains unclear.
The minor venom proteins (<10% of total venom proteins) i.e. snake venom vascular endothelial growth factor (svVEGF), nerve growth factor (svNGF), 5′nucleotidase (5′NUC), cysteine-rich secretory protein (CRiSP) and phosphodiesterase (PDE) may play a role in the predatory or digestive function of the venoms. Some of these venom components were known to induce hypotensive or proinflammatory effects, thereby facilitating the subduing of prey. For instance, svVEGF has been shown to increase capillary permeability 83 and induces hypotensive effect 84 . Meanwhile, svNGF through the release of nitric oxide and/or histamine, 5′NUC through the release of purines (adenosine) and PDE through the alteration of extracellular levels of adenosine and other purines 85   contribute to the venom-induced hypotensive effect or venom spread in the prey. Furthermore, the release of adenosine may cause platelet aggregation inhibition, thereby worsening the hemostatic derangement in Russell's viper envenomation 85 . Also, it has been shown that PDE could strongly inhibit ADP-induced platelet aggregation in human plasma, hence potentiating the hemotoxic effect of the venom 86 . On the other hand, CRiSP (detected in Ds-Guangxi venom) was similar to ablomin of Gloydius blomhoffii, a minor venom protein that targets L-type voltage-gated calcium channels (Ca v ) and blocks smooth muscle contraction 87 . The toxicological property of this protein in D. siamensis venom remains to be further studied. In mainland China, there are two widely distributed medically important pit vipers namely D. acutus (five-pace snake, also known as sharp-nosed pit viper) and G. brevicaudus (short-tail pit viper/Chinese mamushi). Specific antivenoms effective against these two pit vipers have been developed for clinical treatment [88][89][90] . These antivenoms for pit vipers (DaMAV for D. acutus and GbMAV for G. brevicaudus) were anecdotally reported to have been used as alternative antidote to treat D. siamensis envenomation in mainland China in the absence of the specific D. siamensis antivenom. In this study, the weak immunoreactivity of DaMAV and GbMAV against the D. siamensis venoms indicated that the two heterologous antivenoms had very limited binding activity toward the venom antigens of D. siamensis. This is likely due to the substantial differences in the compositions and antigenicity of venom proteins between D. siamensis and the two Chinese pit vipers. For instance, the PLA 2 subtypes reported from the venoms of G. brevicaudus and D. acutus have amino acid sequences varied from those of D. siamensis in this study 91,92 . KSPI which formed the major component of D. siamensis venom have not been reported from the two pit viper venoms. Furthermore, the PIII-SVMP subtype was the main SVMP form expressed in both Ds-Guangxi and Ds-Taiwan venoms. However, in the Chinese G. brevicaudus venom, close to 65% of the total venom proteins were made up of a mixture of PII and PIII-SVMP 91 , while approximately only 10% of PI and PIII-SVMP was reported in the Taiwanese D. acutus venom 93 . The variability of protein families and proteoforms in these venoms imply that the toxins are likely diverse in their antigenicity. Compared with Ds-Taiwan venom, Ds-Guangxi venom was slightly less immunoreactive toward DsMAV-Taiwan, indicating that the antigenicity of some venom proteins varied between the two D. siamensis. On the other hand, the weak binding activity of DsMAV-Thai toward both Ds-Guangxi and Ds-Taiwan venom proteins implied that the Thai D. siamensis venom could be more diverse antigenically from the venoms of the geographically distant Ds-Guangxi and Ds-Taiwan.
The potent procoagulant and hemorrhagic effects of D. siamensis venoms correlated with the hemotoxic syndrome of D. siamensis envenomation in the region. In China, D. siamensis envenoming has been reported to cause severe bleeding including cerebral hemorrhage 20 ; hence, it is essential to ensure that the antivenom used clinically is able to neutralize this hemorrhagic effect of the venom, besides neutralizing the procoagulant and lethal effects. The present study demonstrated that the heterologous DaMAV and GbMAV were rather ineffective in cross-neutralizing the hemotoxic (procoagulant and hemorrhagic) effects of Ds-Guangxi venom, even though GbMAV showed weak cross-neutralizing capability against the venom lethality when given at a very high dosage, judging from its low potency. The finding in this study therefore suggests that both DaMAV and GbMAV may not be the appropriate treatment of D. siamensis envenomation. On the other hand, Ds-Guangxi and Ds-Taiwan venoms shared conserved antigenicity of key toxins, thus enabling DsMAV-Taiwan to neutralize the procoagulant, hemorrhagic and lethal effects of Ds-Guangxi venom effectively.

Conclusion
Shotgun proteomics showed that the principal toxins in the venoms of D. siamensis from Guangxi and Taiwan were comparable. The venom proteins within each protein family, however, varied between the two D. siamensis venoms. The subproteomic variation between the two could be reflective of ecological adaptation to diet since the mainland and the insular populations are physically long separated by the Taiwan Strait. However, the venom divergence could also be the result of random fixation of neutral alleles, with sequence differences that relate to neither dietary nor ecological factors on fitness. Both D. siamensis venoms exhibited potent hemotoxicity and lethality. Ds-Guangxi venom was comparatively more potent in hemorrhagic effect while Ds-Taiwan venom was a stronger procoagulant. Immunoreactivity and neutralization studies further revealed that the antigenicity of the major toxins of D. siamensis were relatively well conserved across the Strait. The study further showed that the heterologous DaMAV and GbMAV did not confer effective cross-neutralization against the venom hemotoxicity and lethality of Ds-Guangxi.

Materials and Methods
Chemicals and materials. All chemicals and reagents used in the study were of analytical grade.
Ammonium bicarbonate, dithiothreitol (DTT) and iodoacetamide were purchased from Sigma-Aldrich (USA). Mass spectrometry sequencing grade of trypsin proteases, and HPLC grade solvents used in the studies were purchased from Thermo Scientific ™ Pierce ™ (USA). Millipore ZipTip ® C 18 Pipette Tips were purchased from Merck (USA).

Whole venom in-solution tryptic digestion and protein identification by tandem mass spectrometry (nano-ESI-LCMS/MS).
Whole venom in-solution tryptic digestion was carried out in three technical replicates for each of the venom. Twenty micrograms for each sample of D. siamensis venoms (Ds-Guangxi and Ds-Taiwan) were subjected to reduction with DTT, alkylation with iodoacetamide, and digested in-solution with mass-spectrometry grade trypsin proteases as described previously 49 . Briefly, the digested peptides eluates were reconstituted in 7 µl of 0.1% formic acid in water. Peptides separation were performed by 1260 Infinity Nanoflow LC system (Agilent, Santa Clara, CA, USA) connected to Accurate-Mass Q-TOF 6550 series with a nano electrospray ionization source. The eluate was subjected to HPLC Large-Capacity Chip Column Zorbax 300-SB-C18 (160 nl enrichment column, 75 µm × 150 mm analytical column and 5 µm particles) (Agilent, Santa Clara, CA, USA). Injection volume was adjusted to 1 µl per sample, using a flow rate of 0.4 µl/min, with linear gradient of 5-70% of solvent B (0.1% formic acid in 100% acetonitrile). Drying gas flow was 11 L/min and drying gas temperature was 290 °C. Fragmentor voltage was 175 V and the capillary voltage was set to 1800 V. Mass spectra was acquired using Mass Hunter acquisition software (Agilent, Santa Clara, CA, USA) in a MS/MS mode with an MS scan range of 200-3000 m/z and MS/MS scan range of 50-3200 m/z. Data were extracted with MH+ mass range between 50 and 3200 Da and processed with Agilent Spectrum Mill MS Proteomics Workbench software packages version B.04.00 against merged database incorporating both non-redundant NCBI databases of Serpentes (taxid: 8570) and in-house transcripts database. Carbamidomethylation was specified as a fixed modification and oxidized methionine as a variable modification. The identified proteins or peptides were validated with the following filters: protein score > 20, peptide score > 10 and scored peak intensity (SPI) >70%. Identified proteins were filtered to achieve False discovery rate (FDR) <1% for the peptide-spectrum matches. The proteins identified were classified as toxins or non-toxins according to their putative functions. The abundance of individual venom toxin was estimated based on its mean spectral intensity (MSI) relative to the total MSI of all proteins identified through the in-solution mass spectrometry 49 .

Sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE). SDS-polyacrylamide
gel electrophoresis (SDS-PAGE) was conducted according to method of Laemmli 95 calibrated with the Thermo Scientific Spectra Multicolor Broad Range Protein Ladder (10-260 kDa). The venoms of both D. siamensis (Ds-Guangxi and Ds-Taiwan) was loaded onto a 15% gel and the electrophoresis was performed under reducing condition at 80 V for 2.5 h. Proteins were stained with Coomassie Brilliant Blue R-250 for gel visualization.
Immunological binding assay. Immunological binding activities between venom antigens and antivenoms were examined with an indirect enzyme-linked immunosorbent assay (ELISA). The immunoplate wells were precoated with 10 ng of venom antigens at 4 °C overnight. The plate was then flicked dry and rinsed four times with phosphate-buffered saline with 0.5% Tween ® 20 (PBST). Antivenoms were prepared at a protein concentration of 10 mg/ml each, and 100 µl appropriately diluted antivenom (1:3000) was added to each antigen-coated well, followed by incubation for 1 h at room temperature. After washing the plate four times with PBST, 100 µl of appropriately diluted horseradish peroxidase-conjugated antihorse-IgG (Jackson ImmunoResearch Inc., USA) in PBST (1:8000) was added to the well and incubated for another hour at room temperature. The excess components were removed by washing four times with PBST. A hundred microliters of freshly prepared substrate solution (0.5 mg/mL o-phenylenediamine and 0.006% hydrogen peroxide in 0.1 M citrate-phosphate buffer, pH 5.0) was added per well. The enzymatic reaction was allowed to take place in the dark for 30 min at room temperature. The reaction was subsequently terminated by adding 50 µL of 12.5% sulphuric acid, and the absorbance at 492 nm was read using Tecan Infinite M1000 Pro plate reader (Tecan Laboratories, Switzerland). Immunological binding activity was expressed as percentage of relative absorbance (highest binding activity set as 100%) between two comparing antivenoms in immunological binding toward the both D. siamensis venoms (Ds-Guangxi and Ds-Taiwan) 96 . Values were means of triplicates ± S.E.M., and the significance of difference was analyzed using unpaired t-test with p value <0.05.  96,97 . The absorbance at 405 nm was monitored every 30 s over 30 min. This produced a plot of absorbance versus time (min) with an initial lag time at which absorbance began to increase drastically due to cloudiness of clot formation, followed by a hyperbolic curve to a plateau. The clotting time was determined as the time when the absorbance became 0.02 U greater than the mean of the first two absorbance measurements. The minimal clotting dose (MCD) was defined as the dose of venom that induces coagulation in 5 min.
In neutralization study, venom concentration equivalent to 2 times minimal coagulation dose (2 MCD) was pre-incubated with various doses of the antivenoms (DsMAV-Taiwan, GbMAV and DaMAV) at 37 °C for 30 min before the addition of 100 µl citrated human plasma. The determination of clotting time in the neutralization assay was performed as described above. The effective dose (ED) was defined as the dose of antivenom that prolonged the clotting time of human plasma 3 times that of the control (2 MCD, without antivenom).
The effective dose (ED) was calculated using the following formulation: (1) Effective dose (ED, mg/ml) 2 times minimal coagulation dose (2 MCD, mg) Dose of antivenom that prolonged the clotting time 3 times that of the control (ml) = For comparative purpose, ED values of antivenoms were normalized (n-ED) by their respective protein amount and expressed as milligram of venom neutralized per gram of antivenom protein (mg/g).
Hemorrhagic activity of D. siamensis venom and antivenom neutralization. Hemorrhagic effect was assessed by intradermal venom injection into the dorsal skin of ICR mice (20-25 g, n = 3) as described by Gutiérrez et al. 98 . The animals were euthanized with urethane 90 min after venom exposure and the skins were removed. Minimal hemorrhagic dose (MHD) was defined as the amount of venom that induces a skin hemorrhagic lesion of 10 mm diameter. For neutralization assays, various doses of antivenom (DsMAV-Taiwan, GbMAV and DaMAV) were pre-incubated with a constant amount of venom challenge dose (2MHD) at 37 °C for 30 min prior to intradermal injection into the animals. The neutralization of hemorrhagic effects was expressed as median effective dose (ED 50 ), defined as the amount of reconstituted antivenom in µl or the ratio of mg venom/ ml reconstituted antivenom in which the venom activity was reduced by 50%.
The median effective dose (ED 50 ) was calculated using the following formulation: = Median effective dose (ED , mg/ml) 2 times minimal hemorhagic dose (2 MHD, mg) Dose of antivenom in which the venom activity was reduced by 50% (ml) (2) 50 For comparative purpose, ED 50 values of antivenoms were normalized (n-ED 50 ) by their respective protein amount and expressed as milligram of venom neutralized per gram of antivenom protein (mg/g).

Determination of venom lethality and neutralization by antivenom.
Median lethal doses (LD 50 ) of venoms were determined by intravenous injection (via caudal veins) into ICR mice (n = 4 per dose, 20-25 g). The survival ratio was recorded after 48 h. In antivenom neutralization assay, pre-incubation of venom and antivenom was conducted as described by Tan et al. 11 . A challenge dose at 5 times LD 50 of the venom dissolved in normal saline was pre-incubated with various dilutions of antivenoms (DsMAV-Taiwan, GbMAV and DaMAV) at 37 °C for 30 min and the mixture was injected intravenously into the mice (n = 4 per dose, 20-25 g). The mice were allowed free access to food and water ad libitum, and the ratio of survival was recorded at 48 h post injection. Neutralizing capacity was expressed as ED 50 , defined as the amount of reconstituted antivenom that gives 50% survival in the venom-challenged animals. These parameters were calculated according to the Probit analysis method 99 using BioStat 2009 analysis software (AnalystSoft Inc., Canada). Neutralization capacity was also expressed in term of 'neutralization potency' (P, defined as the amount of venom in milligram neutralized completely by a unit volume of antivenom in millilitre, mg/ml) 96,100 . The neutralization potency (P) is a more direct indicator of antivenom neutralizing capacity, and is theoretically unaffected by the number of LD 50 in the challenge dose. For comparative purpose, P values of antivenoms were normalized (n-P) by their respective protein amount and expressed as milligram of venom neutralized per gram of antivenom protein (mg/g).