Skip to main content
Advertisement
  • Loading metrics

Use of Molecular Diagnostic Tools for the Identification of Species Responsible for Snakebite in Nepal: A Pilot Study

  • Sanjib Kumar Sharma,

    Affiliation B.P. Koirala Institute of Health Sciences, Dharan, Nepal

  • Ulrich Kuch,

    Affiliation Institute of Occupational Medicine, Social Medicine and Environmental Medicine, Goethe University, Frankfurt am Main, Germany

  • Patrick Höde,

    Affiliation Institute of Occupational Medicine, Social Medicine and Environmental Medicine, Goethe University, Frankfurt am Main, Germany

  • Laura Bruhse,

    Affiliation Institute of Occupational Medicine, Social Medicine and Environmental Medicine, Goethe University, Frankfurt am Main, Germany

  • Deb P. Pandey,

    Affiliations Senckenberg Biodiversity and Climate Research Centre (BiK-F), Frankfurt am Main, Germany, Senckenberg Forschungsinstitut, Frankfurt am Main, Germany

  • Anup Ghimire,

    Affiliation B.P. Koirala Institute of Health Sciences, Dharan, Nepal

  • François Chappuis,

    Affiliation Division of Tropical and Humanitarian Medicine, University Hospitals of Geneva, Geneva, Switzerland

  • Emilie Alirol

    emilie.alirol@hcuge.ch

    Affiliations Division of Tropical and Humanitarian Medicine, University Hospitals of Geneva, Geneva, Switzerland, Médecins Sans Frontières UK, London, United Kingdom

Abstract

Snakebite is an important medical emergency in rural Nepal. Correct identification of the biting species is crucial for clinicians to choose appropriate treatment and anticipate complications. This is particularly important for neurotoxic envenoming which, depending on the snake species involved, may not respond to available antivenoms. Adequate species identification tools are lacking. This study used a combination of morphological and molecular approaches (PCR-aided DNA sequencing from swabs of bite sites) to determine the contribution of venomous and non-venomous species to the snakebite burden in southern Nepal. Out of 749 patients admitted with a history of snakebite to one of three study centres, the biting species could be identified in 194 (25.9%). Out of these, 87 had been bitten by a venomous snake, most commonly the Indian spectacled cobra (Naja naja; n = 42) and the common krait (Bungarus caeruleus; n = 22). When both morphological identification and PCR/sequencing results were available, a 100% agreement was noted. The probability of a positive PCR result was significantly lower among patients who had used inadequate “first aid” measures (e.g. tourniquets or local application of remedies). This study is the first to report the use of forensic genetics methods for snake species identification in a prospective clinical study. If high diagnostic accuracy is confirmed in larger cohorts, this method will be a very useful reference diagnostic tool for epidemiological investigations and clinical studies.

Author Summary

Snakebite is an important medical problem in sub-tropical and tropical regions, including Nepal where tens of thousands of people are bitten every year. Snakebite can result in life-threatening envenoming, and correct identification of the biting species is crucial for care providers to choose appropriate treatment and anticipate complications. This paper explores a number of methods, including molecular techniques, to assist care providers in identifying the species responsible for bites in rural Nepal. Out of 749 patients with a history of snakebite, the biting species could be identified in 194 (25.9%). Out of these, 87 had been bitten by a venomous snake, most commonly cobras (n = 42) and kraits (n = 22). This study is the first to report the use of molecular techniques for snake species identification. The diagnostic accuracy of this method appears high but needs to be confirmed in larger studies.

Introduction

In rural Nepal snakebite is an important public health problem. A survey conducted in the 1980s showed that about 20’000 people were bitten each year, resulting in over 1’000 deaths [1]. These official figures, however, significantly underestimate the true burden [27]. Annual incidence and mortality figures of 1’162/100’000 and 162/100’000, respectively, have been reported in some regions of Nepal [4]. Children are among the primary victims [2,8,9]. Most snakebites occur in the southern plains of Terai, a region with a hot tropical climate and high population density. Most bites occur during the rainy season, from June to September, which corresponds to the peak period for agricultural work.

Eighty nine snake species have been recorded in Nepal, of which 17 are known to be venomous [10]. The Elapidae family includes two species of cobra (Naja naja and Naja kaouthia), the king cobra (Ophiophagus hannah), one species of coral snake (Sinomicrurus macclellandii) and six species of krait (Bungarus bungaroides, Bungarus caeruleus, Bungarus fasciatus, Bungarus lividus, Bungarus niger, and Bungarus walli). The Viperidae are represented by seven species. The most dangerous of these, Russell’s viper (Daboia russelii), appears to be rare in Nepal. Pitvipers (Gloydius himalayanus, Himalayophis tibetanus, Ovophis monticola, Protobothrops spp. Trimeresurus albolabris and Trimeresurus septentrionalis), on the other hand, are widespread from the lowlands to the high mountains. Non-venomous species are also very common and may be involved in bites. Some of these non-venomous species are easily mistaken for venomous ones [10]. For example, rat snakes (Ptyas and Coelognathus species) may be confused with cobras, while wolf snakes, which are common inside and around houses, have a colour pattern similar to that of kraits [1113]. Bites can therefore be inflicted by a variety of species, in all kinds of environments. Neither the geographical distribution of these species nor their contribution to snakebite mortality and morbidity have been systematically studied in Nepal.

Snakebite envenoming can be life-threatening, and recognizing early signs of systemic envenoming is crucial for the optimal management of patients. Depending on the species of snakes, different organs and tissues can be affected [14]. In envenoming following the bites of elapid snakes, neurotoxicity with progressive descending paralysis is characteristic, and patients usually die of respiratory failure if not adequately ventilated [1517]. Clinical prognosis and response to antivenom depend on the species, hence knowing which snake is responsible for the bite is of primary importance [14]. In South Asia appropriate tools for species identification are not available. The snake is rarely seen, and if it is, its description by the victim is often misleading [18]. Hospital personnel are generally not trained to properly identify the biting species, even when it is killed and brought by the victim [1921]. The morphological resemblance of venomous by non-venomous species complicates this task.

Several approaches have been investigated to improve snake species identification in peripheral health centres. Immunodiagnosis of circulating venom antigen can be used to identify the biting species or ‘immunogroups’ of antigenically similar, cross-reacting venoms [2225], but commercial point-of-care venom detection kits have been marketed only in Australia, for species occurring on that continent [26]. In Sri Lanka, syndromic approaches have been proposed [20] and clinical scores developed based on a systematic analysis of features of envenom envenoming in cases where snake specimens brought to hospital were expertly identified [27]. However, the extent to which these approaches can also be applied to the Nepal setting is not known. The present study aimed at investigating the accuracy, feasibility and usefulness of different, complementary approaches to species identification in patients bitten by snakes in rural Nepal. The main objectives were (1) to clarify the contribution of different snake species causing envenoming and non-envenoming bites through morphological analysis of preserved snake specimens and the use of forensic methods for DNA-based identification, and (2) to explore the utility, diagnostic performance and feasibility of DNA-based identification by forensic methods.

Methods

Study setting and population

The study was conducted in three centres, namely the Snake Bite Treatment Centre of Damak Red Cross Society and the Snake Bite Management Centre of Charali, both in Jhapa district, and Bharatpur District Hospital, Bharatpur, Chitwan district. All three centres are located in the Terai lowland region of Nepal.

All snakebite victims presenting to the Damak Red Cross Centre and to the Charali Snake Bite Management Centre during the study period were offered to participate in the study, irrespective of the presence and nature of envenoming signs. Hence, victims of bites by elapids, viperids and non-venomous snakes were included. Patients were excluded from the study if they were below 5 years of age or if they had already received antivenom prior to admission. Patients unable or unwilling to give consent were also excluded. In Bharatpur District Hospital, additional exclusion criteria were applied due to the contemporary conduct of a randomized controlled trial on neurotoxic envenoming (clinicaltrials.gov number NCT01284855). In this centre, only snakebite patients presenting with signs of systemic neurotoxic envenoming were included in the study. Those presenting more than 24 hours after the bite, pregnant or breastfeeding women and patients with proven viper bite envenoming were excluded. The recruitment period in Bharatpur was shorter and extended from 1st April 2011 to 31st October 2012.

Species identification

Dead snakes brought by bite victims were systematically labelled with patient number, initials, date of birth and date of admission, and preserved in 70% ethanol. Morphological identification was conducted by taxonomic experts (UK, DP) who remained blind to the circumstances of bites and to envenoming status of the victims. For each preserved specimen, morphological characters of scalation and dentition were analysed by comparison to relevant reference specimens in museum collections and an existing database [28,29]. Additionally, genetic information from tissue samples of killed snakes were obtained in the form of DNA ‘barcodes’ of the mitochondrial cytochrome b gene [28,29].

Molecular analysis

Whenever the bite site could be located, trace DNA of the biting snake was collected by rubbing the cotton swab of a Prionix evidence collection tube on the bite site (see Standard Operating Procedure in supplementary material). The sample was left to dry at room temperature. For DNA extraction, one half of the cotton bud was cut off and subjected to lysis in CTAB (cetyl trimethylammonium bromide) buffer with proteinase K at 56°C, followed by one extraction with phenol-chloroform-isoamyl alcohol and then with chloroform-isoamyl alcohol. DNA was precipitated with 98% ethanol and 3M sodium acetate (pH 5,2) and stored overnight at -20°C before washing with 70% ethanol, drying and dissolving in TE-buffer. For PCR we used primers flanking a 400 bp sequence of the mitochondrial cytochrome b gene (cytb) which preferentially amplify snake rather than human cytb sequences. For nested PCR we designed primers with cytb binding sites upstream and downstream of the PCR primers. PCR was performed in 30 μl volumes containing 6 μl of template DNA solution, 0.3 pM of each primer, 83 pM of each dNTP, 3 μl 10 × PCR buffer without MgCl2 (Fermentas), 3.25 mM MgCl2, 1 unit TrueStart Hot Start Taq polymerase (Fermentas), and 13.7 μl H2O. PCR products were visualized by SYBRGreen (Invitrogen) staining in 1% agarose gels and ultraviolet (302 nm) light illumination. Their lengths were estimated using a 100 bp ladder (Fermentas). Cycle sequencing of both strands was performed with 0.16 μl of BigDyeTerminator 3.1 reaction mix (Applied Biosystems), 1 μl primer, 1.92 μl 5 × sequencing buffer (Applied Biosystems) and 5.92 μl milliQ H2O (Millipore). Products of the sequencing reaction were separated on an ABI 3730 sequencer (Applied Biosystems) operated with a 50 cm capillary, 8 sec injection time and 1500 V injection voltage. After visual checks of electropherograms, correction of base-calls and comparison of complementary strands using BioEdit and SequenceScanner software, sequences were submitted to BLAST searches for comparison with the sequences in the GenBank, ENA and DDBJ nucleotide sequence database to determine the snake species. New DNA sequences obtained in the course of this study were deposited in this database. A full description of the methods is included in Melaun and Kuch [30].

Features of envenoming

Local envenoming was defined as the presence of one or more of the following: (1) necrosis, (2) bullae or blisters, (3) enlarged regional lymph nodes plus either local bleeding or ecchymosis or swelling and (4) swelling extending at least halfway between two articulations. Snakebite victims presenting with moderate local swelling (i.e., not extending farther than one articulation) were not considered as locally envenomed. Systemic envenoming was defined as the presence of one or more of the following: (1) Incoagulable blood as indicated by the 20 Minutes Whole Blood Clotting Test (20’WBCT), (2) spontaneous and continuous bleeding from the bite site, IV line, gums or old wound, (3) gastrointestinal bleeding or blood in urine, (4) neurotoxic sign(s) including inability to frown, bilateral ptosis, inability to open the mouth, inability to protrude the tongue beyond incisors, inability to clear secretions, broken neck sign, skeletal muscle weakness, gag reflex loss and paradoxical breathing. Neurotoxic envenoming was defined as the presence of one or more of the above-mentioned neurotoxic signs. Envenoming status was determined retrospectively, by one of the authors (EA) who was blinded to the results of snake species identification.

Clinical management of patients

The clinical management of snakebite victims followed the Nepal national protocol and WHO SEARO guidelines [14]. All patients presenting with systemic envenoming received antivenom, as did patients presenting with severe local envenoming (e.g., edema extending over full limb). Antivenoms available in Nepal are all manufactured in India and target four species, the Indian spectacled cobra (Naja naja), the common Indian krait (Bungarus caeruleus), the saw-scaled viper (Echis carinatus) and Russell’s viper (Daboia russelii). Additional treatments included anticholinesterases and assisted ventilation for those patients experiencing respiratory paralysis. Patients who were asymptomatic on admission were kept under observation for 24 hours.

Data collection and analyses

A standard Case Report Form (CRF) was designed to prospectively collect data on the circumstances of the bite, the participants’ demographic characteristics, and the clinical features on admission. First aid measures, whether appropriate or not, were recorded, as well as features of clinical management after admission. To ensure harmonization of data collection, study staff were trained on CRF completion and were instructed to follow Standard Operating Procedures (SOP). Separate forms were used for the results of morphological identification and molecular analysis. As these analyses were performed well after recruitment was over, their results did not influence the assessment of snakebite victims by the care-providers. The researchers performing the PCR and DNA sequence analyses were also blind to the results of morphological identification and patient data.

Characteristics on admission were described using percentages for categorical variables and by calculating means and standard deviations or median and Inter Quartile Ranges (IQR) for continuous variables. The Case Fatality Rate (CFR) was calculated as the percentage of patients who died among those showing signs of envenoming. Descriptive analyses were restricted to participants with snake species identified by PCR or morphology. Determinants of PCR positivity were analysed for all patients for whom a PCR was done, as follows: continuous variables were compared using Student’s t test, or the Wilcoxon- rank sum test for non-normally distributed variables. Categorical variables were compared using the Χ2 test with continuity correction or the Fisher exact test, as appropriate. For categorical variables with more than one category, a test for trend was used. Associations were examined at the p<0.05 level of significance. To assess the probability of a positive PCR, Risk Ratios (RR) with 95% Confidence Interval (95% CI) were calculated.

Ethical aspects

The study was conducted in accordance with the Declaration of Helsinki 1964, as revised in Seoul, 2008, and in compliance with the protocol, Good Clinical Practices (GCP) and Nepal regulatory requirements. The B.P. Koirala Institute of Health Sciences (BPKIHS) Ethics Committee, the Nepal National Health Research Council (NHRC) and Geneva University Hospital Ethics Committee approved the study prior to its start. All participants provided written informed consent before being included in the study.

Results

Between the 1st of April 2010 and the 31st of October 2012, a total of 749 patients were found eligible to be included in the study (Fig 1).

thumbnail
Fig 1. Flow diagram showing numbers of individuals screened and included in each study centre.

Between 01/04/2010 and the 31/10/2012, 749 victims of snakebite were included in the study and the snake species responsible for the bite could be ascertained in 194 cases.

https://doi.org/10.1371/journal.pntd.0004620.g001

The responsible snake species could be ascertained for 194 (25.9%) patients. In total 722 snakebite victims presented with fang marks, and swabs were collected and analysed from 565 of them. PCR and DNA sequencing yielded a positive result in 153 out of these (27.1%). Sixty two snakes were brought by the victim and all were morphologically identified. Collectively the snake species could be determined in 87 (32.9%) of 264 envenomed patients and in 107 (22.1%) of 485 non-envenomed patients. Table 1 lists the species identified in all three centers. Among the non-venomous species, the checkered keelback (Xenochrophis piscator) was responsible for the majority of bites. Among the venomous species identified the most common were the spectacled cobra (N. naja) and common krait (B. caeruleus). There was a very clear difference between study sites. Most B. caeruleus specimens (20 out of 22) were brought in Bharatpur, while most N. naja (35 out of 42) came from either Charali or Damak. In 21 cases, both morphological identification and a snake DNA sequence were available. In these cases there was a 100% agreement between the two identification methods.

thumbnail
Table 1. List of species responsible for snakebite in southern Nepal between April 2010 and October 2012 (n = 194).

Species were identified either by morphological examination of preserved specimen or through PCR and DNA sequencing performed using bite-site swabs.

https://doi.org/10.1371/journal.pntd.0004620.t001

Epidemiological characteristics

The baseline characteristics and circumstances of the bites of the 194 snakebite victims for whom a species could be ascertained, are summarized in Table 2. The majority of snakebite victims were male (52.6%) and the median age was 30.8 years (IQR = 18–40). There were 28 (14.4%) children (5–15 years). A clear seasonal pattern was observed in all three study sites, with the incidence of bites being higher during the rainy season (June to September).

thumbnail
Table 2. Circumstances of the bite and baseline characteristics of victims of snakebite in southern Nepal between April 2010 and October 2012 (n = 194).

https://doi.org/10.1371/journal.pntd.0004620.t002

The median time needed to reach the treatment centre was 85 minutes (IQR = 50 to 120 minutes), with a longer time observed in Bharatpur (194 minutes; IQR: 11 to 269 minutes). Only four (2.1%) snakebite victims visited a traditional healer before reaching the study centre. Use of first aid measures was more common in Damak (n = 91/93; 97.8%) and Charali (n = 72; 98.6%) compared to Bharatpur (n = 11; 39.5%). Tourniquets were by far the most common first aid method applied.

S1 Table compares baseline characteristics, circumstances of the bite and first aid measures between snakebite victims with identified and unidentified snake species.

Clinical features on admission

Among snakebite victims presenting to the treatment centres in Damak or Charali during the study period (n = 676), 191 (27.4%) showed signs of envenoming. Twenty three (24.7%) patients in Damak and four (5.5%) in Charali presented with neurotoxic signs. As expected by local inclusion criteria, all 73 patients recruited in Bharatpur had signs of systemic neurotoxic envenoming on admission, and among those, seven also had local signs (see S1 Text).

The signs and symptoms presented on admission by the 194 snakebite victims for whom a species identity could be ascertained are described in Table 3. Among 87 victims for whom a venomous species was identified, 16 (18.4%) had not developed any sign of envenoming. Among the 73 patients for whom an elapid snakebite could be ascertained, 64 (87.7%) developed signs of envenoming. As expected, none of the 107 patients bitten by non-venomous species exhibited local or systemic signs of envenoming.

thumbnail
Table 3. Clinical features on admission of victims of snakebite in southern Nepal between April 2010 and October 2012 (n=194).

Initial assessment of snakebite victims included vital signs, patients’ complaints and standardized evaluation of envenoming signs.

https://doi.org/10.1371/journal.pntd.0004620.t003

Swelling of the bitten limb was present in 43 out of 194 patients (22.2%), and bleeding from the bite site was present in 60 (30.9%). Of note, both signs also occurred in patients bitten by non-venomous species, in particular 25 out of 88 X. piscator bite cases (28.4%) presented with bleeding from the bite site.

The 20 Minutes Whole Blood Clotting test was performed on all victims presenting to one of the study centres during the study period. Out of those with incoagulable blood on admission, the snake species could be determined in 7 cases: 5 were N. naja and 2 were white-lipped pit vipers (Trimeresurus cf. albolabris). No snakebite victim presented with gum bleeding, gut bleeding, or blood in urine on admission.

S2 Table compares clinical features on admission of snakebite victims with identified and unidentified snake species.

Factors affecting PCR positive rate

As the number of cases in which the snake species could be identified by PCR sequencing was low, we investigated baseline characteristics and circumstances of the bite that may influence the sensitivity of this method. Results are summarized in Table 4. The median time to reach the center was significantly longer in patients with a negative PCR. The probability of a positive PCR result was significantly lower among patients who had used first aid measures. In particular, applying local remedies (e.g., herbs, honey, etc.) was associated with a 2-fold decrease in the probability of a positive PCR. Patients bitten on the upper limb were 60% less likely to have a positive PCR compared to those bitten on the lower limb. Local bleeding also increased the probability of a positive PCR by 1.45 fold. Whereas bite by a venomous species did not affect the chances of a positive PCR (RR = 1.004, 95%CI: 0.893–1.128), showing signs of envenoming significantly increased the chances of a positive result (RR = 1.444, 95% CI: 1.091–1.912, p = 0.017).

thumbnail
Table 4. Factors associated with a positive PCR among 565 snake bite victims.

Unadjusted Risk Ratio (RR) and their 95% Confidence Interval (95% CI) were calculated with respect to the baseline category, i.e., absence of the risk factor (RR = 1).

https://doi.org/10.1371/journal.pntd.0004620.t004

Management and outcome

A total of 67 (34.5%) snakebite victims received antivenom. Eight were put under mechanical ventilation and nine were transferred to a tertiary care centre. Among patients treated with antivenom the median dose was 2 vials (inter-quartile-range = 1–3).

In total there were 5 deaths (Case Fatality Rate = 5/194 = 2.6%): one in Damak, one in Charali and three in Bharatpur). Of these 5 deaths, 4 had been bitten by a B. caeruleus and 1 by a N. naja.

Discussion

Between the 1st of April 2010 and the 31st of October 2012, 749 individuals with a history of snakebite presented to one of three study centres in southern Nepal. In 194 (25.9%) patients, the responsible snake species could be ascertained, by either morphological identification or PCR plus DNA sequencing. Most species identified were non-venomous ones. The non-venomous checkered keelback (X. piscator) was the most frequently identified species, followed by the spectacled cobra (N. naja) and the common krait (B. caeruleus). Other venomous species contributing to the snakebite burden in this study comprised several pitvipers (O. monticola, Trimeresurus sp., T. albolabris, and T. popeiorum) as well as various additional elapid snakes, including the first cases of envenoming by the greater black krait (Bungarus niger) and the king cobra (O. hannah) ever reported in Nepal.

As few victims (11.6%) brought the dead snakes to the study centres, the use of PCR to amplify snake trace DNA from bite-site swabs increased to 25.9% the proportion of victims for whom the snake species could be ascertained. Morphological identification of preserved specimens by a qualified herpetologist is the gold standard for species identification. However, this method is seldom used, as snakes are rarely captured and preserved, and as care-providers working in snakebite treatment centres generally lack the appropriate expertise [1921]. Alternative approaches must therefore be developed to complement morphological identification. Molecular techniques have shown promising results in animal models [3033], and the present study shows that PCR amplification of a mitochondrial gene region from snake trace DNA is feasible in field setting. Sampling is straight forward and requires minimal training (see SOP in supplementary material), and the storage and transport conditions for the Prionix evidence collection tubes (room temperature and protected from light) can easily be met.

The PCR yielded a positive result in 27.1% of the cases only. This could probably be improved if several factors shown in this study to be associated with a lower sensitivity, such as inappropriate first-aid measures (use of tourniquet or application of local remedies on the bite site) or prolonged time to reach the treatment centre, were corrected or improved by public health interventions. Although the sensitivity of the PCR did not seem to be affected by whether the species was venomous or not, envenoming status did have an effect. This may be due to confounding effects of venom injection. Bites associated with venom injection may go deeper in tissues, and snakes may deposit more trace DNA on bite sites during envenoming bites compared to ‘dry’ bites or bites by non-venomous snakes e.g., snake DNA can also be recovered from snake venom [33].

In all cases where both snake morphological identification and DNA sequence from the bite-site swab were available (n = 21), a 100% agreement was observed between the two methods, suggesting a high specificity of molecular identification. An on-going study conducted by the authors is expected to complement this preliminary data and validate PCR-aided sequencing of snake trace DNA on a larger cohort of patients. Although molecular tools are not yet appropriate for point-of-care (POC) testing and hence cannot be used to guide clinical management, the encouraging results presented here, if confirmed in larger studies, suggest that they could be used as reference tests in future epidemiological and clinical studies. Progress in the development and validation of POC tests for snake species has indeed been hindered by the difficult implementation of the diagnostic gold standard (morphological identification) in rural regions where most bites occur. Besides, molecular tools could be very useful in clarifying the contribution of different snake species to the snakebite burden, and help identify new medically important species.

The species found to have caused neurotoxicity in this study were the two species of cobra (N. naja and N. kaouthia) and three different species of krait (B. caeruleus, B. lividus, and B. niger). Of these, only two (N. naja and B. caeruleus) are included in the production of the Indian polyvalent antivenoms that are available in Nepal. No pre-clinical data exist on the efficacy of Indian polyvalent antivenoms to neutralize the venoms of these two species in Nepal, and no clinical data have been published so far. In this study, patients bitten by a krait had significantly higher chances of being put under mechanical ventilation and being transferred to an intensive care unit compared to those being bitten by a cobra (see supplementary information). This is consistent with published literature which suggests that krait bites often result in poorer outcomes for patients, and high mortality rates [34,35]. The efficacy of the available Indian antivenom in reversing envenoming by kraits is increasingly being questioned, and several case series have reported little or no benefit of immunotherapy [17,3638].

The frequency of neurotoxicity observed in our study is consistent with elapid snake species being most commonly involved in envenoming bites in Nepal. The fact that we did not observe late-appearing signs such as broken neck sign, muscle weakness and loss of gag reflex may be due to the early presentation of victims to the health centre and the prompt initiation of antivenom therapy. Interestingly, among those patients who presented with incoagulable blood on admission were five victims of N. naja bites. This is not the first report of apparent coagulopathy following bites by this species [39], and a few in vitro studies have reported anticoagulant activities in N. naja venoms [4042], however, these results need to be interpreted with caution. In fact, the 20 Minutes Whole Blood Clotting Test can give erroneous results if performed incorrectly, in particular if the tubes used bear traces of detergent [14].

The present study has several limitations, the principal one being that the study population differed between study sites. Snakebite victims admitted to Bharatpur District Hospital were only included if they presented with signs of neurotoxic envenoming. It is therefore not surprising that all snakes identified in this centre were venomous, resulting in an overestimation of the contribution of venomous species (and in particular elapids) to the snakebite burden. When Bharatpur was excluded from the analysis, venomous species accounted for only 35.6% of identified bites. The checkered keelback and the spectacled cobra remained the most common species identified (see supplementary information). Another limitation relates to the geographical coverage of the study. The list of species identified here is not representative of all snakes causing bites in Nepal. The three study centres are located in the lowlands of the Terai region where most snakebites occur, but their catchment areas do not cover all of Nepal’s great biogeographical diversity. In particular, species found in mountain regions (although present, e.g., O. monticola) were probably under-represented in our study.

The fact that the biting species could be identified only in a relatively small proportion of patients (25.9%) could in theory lead to bias. Although we cannot exclude that some selection bias occurred in the present study, its impact is likely to be minimal. We compared bite circumstances and baseline characteristics of snakebite victims with or without identification of snake species (S1 and S2 Tables). Differences were seen with regard to season of bite, location and activity at the time of bite and consultation of a traditional healer. However, the magnitude of these differences was minimal. Moreover, the epidemiological characteristics of our study population are consistent with other published reports [25,8,43,44], further ruling out the possibility of selection bias and reaffirming the external validity of our findings.

Finally, morphological identification and molecular analysis results were both available in only 21 cases, limiting our ability to evaluate the diagnostic performance of the molecular diagnosis method. Findings presented here thus need to be interpreted with caution. A follow-up prospective validation study is ongoing in Nepal and Myanmar to address this issue.

Snakebite envenoming is an important health problem in Nepal, accounting for up to 39.7% of poisoning cases admitted to emergency units in certain regions [6]. Neurotoxicity following the bites of elapid snakes is of particular concern. This study for the first time addresses the distribution and medical importance of snake species contributing to the burden of snakebite in Nepal. It provides crucial information for clinicians and health workers involved in the management of snakebite victims in Nepal. It notably highlights that the majority of bites are caused by non-venomous snakes, and that the diversity of venomous snake species involved in bites is greater than previously believed. Finally, this study provides initial evidence on the utility of forensic DNA-based methods in the identification of biting snake species.

Supporting Information

S2 Text. Standard Operating Procedure (SOP) for the collection, storage and shipment of swab samples.

https://doi.org/10.1371/journal.pntd.0004620.s003

(PDF)

S1 Table. Comparison of baseline characteristics, circumstances of the bite and first-aid measures between snake bite victims with identified (n = 194) and unidentified (n = 555) snake species.

https://doi.org/10.1371/journal.pntd.0004620.s004

(DOCX)

S2 Table. Comparison of clinical features on admission of victims of snake bite with identified (n = 194) and unidentified (n = 555) snake species.

https://doi.org/10.1371/journal.pntd.0004620.s005

(DOCX)

Acknowledgments

We are grateful to Mamit Rai for his kind help with project coordination and data entry. We also thank all paramedical staff involved in managing the snakebite patients in the three study centres.

Author Contributions

Conceived and designed the experiments: EA FC UK SKS. Performed the experiments: SKS PH LB DPP AG. Analyzed the data: EA PH LB UK. Contributed reagents/materials/analysis tools: SKS PH LB UK. Wrote the paper: EA FC SKS. Critical review of the manuscript: UK PH LB DPP AG.

References

  1. 1. Zoonotic disease control: Baseline epidemiological study on snake-nite treatment and management. Wkly Epidemiol Rec. 1987;42: 319–320.
  2. 2. Sharma SK, Khanal B, Pokhrel P, Khan A, Koirala S. Snakebite-reappraisal of the situation in Eastern Nepal. Toxicon. 2003;41: 285–9. Available: http://www.ncbi.nlm.nih.gov/pubmed/12565750 pmid:12565750
  3. 3. Sharma SK, Koirala S, Dahal G, Sah C. Clinico-epidemiological features of snakebite: a study from Eastern Nepal. Trop Doct. 2004;34: 20–2. Available: http://www.ncbi.nlm.nih.gov/pubmed/14959965 pmid:14959965
  4. 4. Sharma SK, Chappuis F, Jha N, Bovier PA, Loutan L, Koirala S. Impact of snake bites and determinants of fatal outcomes in Southeastern Nepal. Am J Trop Med Hyg. 2004;71: 234–238. pmid:15306717
  5. 5. Pandey DP. Epidemiology of snakebites based on field survey in Chitwan and Nawalparasi districts, Nepal. J Med Toxicol. 2007;3: 164–8. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3550019&tool=pmcentrez&rendertype=abstract pmid:18072170
  6. 6. Gurung CK, Dahal R, Khanal P, Nepal S, Jaiswal AK. Pattern of Poisoning Cases in a Hospital in a Terai District of Central Nepal. Nepal Med Coll J. 2011;13: 160–163. pmid:22808805
  7. 7. Pyakurel R, Sharma N, Paudel D, Coghill A, Sinden L, Bost L, et al. Cause of Death in Women of Reproductive Age in Rural Nepal Obtained Through Community-Based Surveillance: Is Reducing Maternal Mortality the Right Priority for Women’s Health Programs? Health Care Women Int. 2014; 1–8.
  8. 8. Magar CT, Devkota K, Gupta R, Shrestha RK, Sharma SK, Pandey DP. A hospital based epidemiological study of snakebite in Western Development Region, Nepal. Toxicon. 2013;69: 98–102. pmid:23624195
  9. 9. Paudel KM, Sharma S. Study of Clinico-Epidemiological Profile and Outcome of Poisonous Snake Bites in Children. 2012;32: 47–52.
  10. 10. Sharma SK, Pandey DP, Tillack F, Chappuis F, Magar CT, Alirol E, et al. Venomous snakes of Nepal: A photographic guide [Internet]. Lalitpur, Nepal: BP Koirala Institute of Health Sciences; 2013. Available: www.tinyurl.com/nepalsnakes
  11. 11. de Silva HA. Coulour guide to the snakes of Sri Lanka. Portishead, UK: R & A Publishing Ltd; 1990.
  12. 12. Fernando B. Kraits and wolfsnakes—the karawalas. Loris. 1977;14: 157.
  13. 13. Harris JB, Faiz MA, Rahman MR, Jalil MMA, Ahsan MF, Theakston RDG, et al. Snake bite in Chittagong Division, Bangladesh: a study of bitten patients who developed no signs of systemic envenoming. Trans R Soc Trop Med Hyg. 2010;104: 320–7. pmid:20096910
  14. 14. Warrell DA. Guidelines for the management of snake-bites. New Delhi, India: WHO SEARO; 2010.
  15. 15. Kuch U, Sharma SK, Alirol E, Chappuis F. Fatal neurotoxic envenomation from the bite of a Lesser Black Krait (Bungarus lividus) in Nepal. Southeast Asian J Trop Med Public Health. 2011;42: 960–964. Available: http://www.tm.mahidol.ac.th/en/seameo/publication.htm\nhttp://ovidsp.ovid.com/ovidweb.cgi?T=JS&CSC=Y&NEWS=N&PAGE=fulltext&D=cagh&AN=20113273840\nhttp://lshtmsfx.hosted.exlibrisgroup.com/lshtm?sid=OVID:caghdb&id=pmid:&id=&issn=0125-1562&isbn=&volume=42&issue=4&spage=960&pages=960-964&date=2011&title=Southeast+Asian+Journal+of+Tropical+Medicine+and+Public+Health&atitle=Fatal+neurotoxic+envenomation+from+the+bite+of+a+Lesser+Black+Krait+(Bungarus+lividus)+in+Nepal.&aulast=Kuch&pid=<author pmid:22299479
  16. 16. Pillai L V, Ambike D, Husainy S, Khaire A, Captain A, Kuch U. Severe Neurotoxic Envenoming and Cardiac Complications after the Bite of a “Sind Krait” (Bungarus cf. sindanus) in Maharashtra, India. Trop Med Health. 2012;40: 103–8. pmid:23264729
  17. 17. Ariaratnam C a, Sheriff MHR, Theakston RDG, Warrell D a. Distinctive epidemiologic and clinical features of common krait (Bungarus caeruleus) bites in Sri Lanka. Am J Trop Med Hyg. 2008;79: 458–62. Available: http://www.ncbi.nlm.nih.gov/pubmed/18784244 pmid:18784244
  18. 18. Kularatne AM, Silva A, Maduwage K, Ratnayake I, Walathara C, Ratnayake C, et al. Victims’ response to snakebite and socio-epidemiological factors of 1018 snakebites in a tertiary care hospital in Sri Lanka. Wilderness Environ Med. Elsevier; 2014;25: 35–40. pmid:24412659
  19. 19. Viravan C, Looareesuwan S, Kosakarn W, Wuthiekanun V, Mccarthy CJ, Stimson AF, et al. A national hospital- based survey of snakes responsible for bites in Thailand. Trans R Soc Trop Med Hyg. 1992;86: 100–106. pmid:1566285
  20. 20. Ariaratnam CA, Sheriff MHR, Arambepola C, Theakston RDG, Warrell DA. Syndromic approach to treatment of snake bite in Sri Lanka based on results of a prospective national hospital-based survey of patients envenomed by identified snakes. Am J Trop Med Hyg. 2009;81: 725–31. pmid:19815895
  21. 21. Joseph JK, Simpson ID, Menon NCS, Jose MP, Kulkarni KJ, Raghavendra GB, et al. First authenticated cases of life-threatening envenoming by the hump-nosed pit viper (Hypnale hypnale) in India. Trans R Soc Trop Med Hyg. 2007;101: 85–90. pmid:16839578
  22. 22. Theakston RD, Lloyd-Jones MJ, Reid HA. Micro-ELISA for detecting and assaying snake venom and venom-antibody. Lancet. 1977;2: 639–41. Available: http://www.ncbi.nlm.nih.gov/pubmed/71452 pmid:71452
  23. 23. Brunda G, Sashidhar RB, Sarin RK. Use of egg yolk antibody (IgY) as an immunoanalytical tool in the detection of Indian cobra (Naja naja naja) venom in biological samples of forensic origin. Toxicon. 2006;48: 183–94. pmid:16846624
  24. 24. Selvanayagam ZE, Gnanavendhan SG, Ganesh K a, Rajagopal D, Rao P V. ELISA for the detection of venoms from four medically important snakes of India. Toxicon. 1999;37: 757–70. Available: http://www.ncbi.nlm.nih.gov/pubmed/10219987 pmid:10219987
  25. 25. Dong L Van, Selvanayagam ZE, Gopalakrishnakone P, Eng KH. A new avidin-biotin optical immunoassay for the detection of beta-bungarotoxin and application in diagnosis of experimental snake envenomation. J Immunol Methods. 2002;260: 125–36. Available: http://www.ncbi.nlm.nih.gov/pubmed/11792383 pmid:11792383
  26. 26. Chandler HM, Hurrell JG. A new enzyme immunoassay system suitable for field use and its application in a snake venom detection kit. Clin Chim Acta. 1982;121: 225–30. Available: http://www.ncbi.nlm.nih.gov/pubmed/7046996 pmid:7046996
  27. 27. Pathmeswaran a, Kasturiratne a, Fonseka M, Nandasena S, Lalloo DG, de Silva HJ. Identifying the biting species in snakebite by clinical features: an epidemiological tool for community surveys. Trans R Soc Trop Med Hyg. 2006;100: 874–8. pmid:16412486
  28. 28. Kuch U, Keogh JS, Weigel J, Smith L a, Mebs D. Phylogeography of Australia’s king brown snake (Pseudechis australis) reveals Pliocene divergence and Pleistocene dispersal of a top predator. Naturwissenschaften. 2005;92: 121–7. pmid:15688185
  29. 29. Kuch U, Kizirian D, Truong NQ, Lawson R, Donnelly MA, Mebs D. A New Species of Krait (Squamata: Elapidae) from the Red River System of Northern Vietnam. Copeia. 2005;4: 818–833. Available: http://www.asihcopeiaonline.org/doi/abs/10.1643/0045-8511(2005)005[0818:ANSOKS]2.0.CO;2
  30. 30. Melaun C, Kuch U. Snake species identification by PCR amplification and sequence analysis of mitochondrial DNA from bite-site swabs. PLoS Negl Trop Dis.
  31. 31. Fernandes C a., Ginja C, Pereira I, Tenreiro R, Bruford MW, Santos-Reis M. Species-specific mitochondrial DNA markers for identification of non-invasive samples from sympatric carnivores in the Iberian Peninsula. Conserv Genet. 2007;9: 681–690.
  32. 32. Suntrarachun S, Pakmanee N, Tirawatnapong T, Chanhome L, Sitprija V. Development of a polymerase chain reaction to distinguish monocellate cobra (Naja khouthia) bites from other common Thai snake species, using both venom extracts and bite-site swabs. Toxicon. 2001;39: 1087–90. Available: http://www.ncbi.nlm.nih.gov/pubmed/11223099 pmid:11223099
  33. 33. Pook CE, McEwing R. Mitochondrial DNA sequences from dried snake venom: a DNA barcoding approach to the identification of venom samples. Toxicon. 2005;46: 711–5. pmid:16157361
  34. 34. Warrell DA. Clinical toxicology of Snake bites in Asia. In: White MA, editor. Handbook of clinical toxicology of animal venoms and poisons. CRC Press; 1995. pp. 493–588.
  35. 35. Bawaskar HS, Bawaskar PH. Envenoming by the Common Krait (Bungarus caeruleus) and Asian Cobra (Naja naja): Clinical Manifestations and Their Management in a Rural Setting. Wilderness Environ Med. Elsevier; 2004;15: 257–266. pmid:15636376
  36. 36. Theakston RD, Phillips RE, Warrell D a, Galagedera Y, Abeysekera DT, Dissanayaka P, et al. Envenoming by the common krait (Bungarus caeruleus) and Sri Lankan cobra (Naja naja naja): efficacy and complications of therapy with Haffkine antivenom. Trans R Soc Trop Med Hyg. 1990;84: 301–8. Available: http://www.ncbi.nlm.nih.gov/pubmed/2389328 pmid:2389328
  37. 37. Kularatne SAM. Common krait (Bungarus caeruleus) bite in Anuradhapura, Sri Lanka: a prospective clinical study, 1996–98. Postgrad Med J. 2002;78: 276–80. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1742360&tool=pmcentrez&rendertype=abstract pmid:12151569
  38. 38. Anil a, Singh S, Bhalla A, Sharma N, Agarwal R, Simpson ID. Role of neostigmine and polyvalent antivenom in Indian common krait (Bungarus caeruleus) bite. J Infect Public Health. King Saud Bin Abdulaziz University for Health Sciences; 2010;3: 83–7. pmid:20701896
  39. 39. Kularatne SAM, Budagoda BDSS, Gawarammana IB, Kularatne WKS. Epidemiology, clinical profile and management issues of cobra (Naja naja) bites in Sri Lanka: first authenticated case series. Trans R Soc Trop Med Hyg. 2009;103: 924–30. pmid:19439335
  40. 40. Kumar MS, Devaraj VR, Vishwanath BS, Kemparaju K. Anti-coagulant activity of a metalloprotease: further characterization from the Indian cobra (Naja naja) venom. J Thromb Thrombolysis. 2010;29: 340–8. pmid:19629641
  41. 41. Sundell IB, Rånby M, Zuzel M, Robinson KA, Theakston RDG. In vitro procoagulant and anticoagulant properties of Naja naja naja venom. Toxicon. 2003;42: 239–47. Available: http://www.ncbi.nlm.nih.gov/pubmed/14559074 pmid:14559074
  42. 42. Asad MHH Bin, Razi MT, Khan T, Najam-Us-saqib Q, Murtaza G, Hussain MS, et al. Coagulopathies in Naja naja karachiensis (black Pakistan cobra) bites and its effect on coagulation tests upon storage of platelet-poor plasma. Acta Pol Pharm. 69: 1031–4. Available: http://www.ncbi.nlm.nih.gov/pubmed/23285662 pmid:23285662
  43. 43. Sharma SK, Bovier P, Jha N, Alirol E, Loutan L, Chappuis F. Effectiveness of rapid transport of victims and community health education on snake bite fatalities in rural Nepal. Am J Trop Med Hyg. 2013;89: 145–50. pmid:23568287
  44. 44. Chappuis F, Sharma SK, Jha N, Loutan L, Bovier P a. Protection against snake bites by sleeping under a bed net in southeastern Nepal. Am J Trop Med Hyg. 2007;77: 197–9. Available: http://www.ncbi.nlm.nih.gov/pubmed/17620654 pmid:17620654