There is a paucity of robust epidemiological data on snakebite, and data available from hospitals and localized or time-limited surveys have major limitations. No study has investigated the incidence of snakebite across a whole country. We undertook a community-based national survey and model based geostatistics to determine incidence, envenoming, mortality and geographical pattern of snakebite in Sri Lanka. The survey was designed to sample a population distributed equally among the nine provinces of the country. The number of data collection clusters was divided among districts in proportion to their population. Within districts clusters were randomly selected. Population based incidence of snakebite and significant envenoming were estimated. Model-based geostatistics was used to develop snakebite risk maps for Sri Lanka. 1118 of the total of 14022 GN divisions with a population of 165665 (0. 8%of the country’s population) were surveyed. The crude overall community incidence of snakebite, envenoming and mortality were 398 (95% CI: 356–441), 151 (130–173) and 2. 3 (0. 2–4. 4) per 100000 population, respectively. Risk maps showed wide variation in incidence within the country, and snakebite hotspots and cold spots were determined by considering the probability of exceeding the national incidence. This study provides community based incidence rates of snakebite and envenoming for Sri Lanka. The within-country spatial variation of bites can inform healthcare decision making and highlights the limitations associated with estimates of incidence from hospital data or localized surveys. Our methods are replicable, and these models can be adapted to other geographic regions after re-estimating spatial covariance parameters for the particular region.
Bites of venomous snakes cause significant morbidity and mortality in the rural tropics. Despite this, snakebite did not, until very recently, receive the attention it deserves as an important public health problem. The main reason for this was the paucity of robust epidemiological data on the disease burden associated with snakebite. The most recent estimates of the global burden of snakebite highlighted the need for good quality data on snakebite, particularly from nation-wide population-based studies [1]. The paucity of reliable data is partly related to inherent methodological difficulties, which include: poorly developed reporting and recording systems in countries with the highest burden, limitations in hospital-based data that often under-estimate the problem [2][3][4][5], and seasonal and geographical variation in bite incidence [6][7], all of which make extrapolations unreliable. Sound epidemiological data are, however, important to both give credence to the magnitude of the problem and raise awareness of snakebite as an important but neglected public health issue, and to assist prioritization of resources for prevention and treatment. Sri Lanka has a mean elevation of 228 meters and consists of mainly low flat to rolling plains with mountains in the south central interior which reach 2400 meters. The country is divided into three climatic zones based on rainfall: wet (south western and central parts of the country including the Western, Central and Sabaragamuwa provinces); dry (northern and south eastern parts of the country including the Northern, Northcentral, Eastern and parts of the Southern and Uva provinces); and intermediate (the areas in between the wet and dry zones including the Northwestern, and parts of Uva and Southern provinces) [8][9]. The differences in rainfall have led to much diversity in the flora and fauna, and in land use in these zones, leading to differences in snakebite patterns [1]. Despite being a country where snakebite results in over 30000 hospital admissions annually [10], the community incidence of snakebite in Sri Lanka is unknown. Previous research has shown that hospital data underestimate deaths due to snakebites by over 50% [4], and multiple hospital admissions at different levels of care can further distort hospital statistics. Geostatistical modeling provides a method to estimate continuous spatial variation in incidence from community level surveys through modeling the occurrence of spatially correlated events in order to predict the probability or likelihood of occurrence of an event anywhere in the study-region and to identify risk-factors [11][12]. The resulting risk maps can provide useful information for healthcare decision making in resource limited settings; they include the majority of the most snakebite endemic regions, without the need for comprehensive registry data. The aims of this study were to determine the community incidence, rate of envenoming, mortality and geographical pattern of snakebite in Sri Lanka and to develop a snakebite risk map for Sri Lanka based on geo-spatial and socio-demographic predictive factors.
Ethical approval for the study was given by the Ethics Review Committee of the Faculty of Medicine, University of Kelaniya. Data collection was done by using an interviewer administered questionnaire and all interviews were conducted after obtaining informed written consent. No animals were used in the study. Permission for conducting the study was obtained from District and Divisional level public administrators before data collection. Grama Niladharis of the sampled GN divisions were informed about the study through the public administration system.
Data relating to 165665 individuals (0. 8% of the population of Sri Lanka) living in 44136 households in 1118 clusters were collected (Fig 1). 695 snakebites, 323 envenomings and 5 deaths (four of them male) were reported in the sample population in the 12 months preceding the interview. The incidence of snakebites, envenoming and deaths in the complete sample population was 398 (95% CI 356–441), 151 (95% CI 130–173) and 2. 3 (95% CI 0. 2–4. 4) per 100000 population respectively. Extrapolating this to the population of the whole country, the estimated national numbers of snakebites, envenoming and deaths were 80514 (95% CI 71774–89254), 30543 (95% CI 26203–34883) and 464 (95% CI 45–884) respectively. There was wide geographical variation in incidence within the country: the Northcentral province (rural, agricultural, dry region) had the highest incidence rate of 623 snakebites per 100000 and the Central province (which includes the high altitude areas of the country) the lowest incidence rate of 277 per 100000 (Table 1). The differences were more marked for envenoming. The Northcentral province had the highest envenoming rate of 440 per 100000 and the Central province the lowest rate of 92 per 100000 (Table 1). The proportion of snakebite victims who were envenomed varied from 20% in Western province to 70% in Northern and Northcentral province. Incidence rates of snakebite and envenoming among males were 478 (95% CI: 415–541) and 176 (95% CI: 146–207) per 100000 population. Males were approximately 1. 5 times more likely to experience both snakebites and envenoming than females (Table 2). Median age (interquartile range) of snakebite victims was 42 (31–54) years and for victims with envenoming bites was 42 (30–53). Almost 40% of the snakebites occurred between 4 pm and 8 pm. Snakebite was rarest in the early morning (4 am to 8 am). More than 55% of the male victims were involved in labour-intensive outdoor occupations such as farming and manual labour. The commonest clinical feature recalled by the snakebite victims was swelling. Symptoms representing neurotoxicity were reported in 28% of all bites. 37% of the significantly envenomed reported neurotoxic features without bleeding. Bleeding manifestations alone were reported by 8%. Neurotoxicity and bleeding together were reported by 24% (Table 3). Exploratory analysis for snakebite incidence using generalized linear models, ignoring spatial dependence, showed population density, elevation, occupation distribution and climatic zones to be significantly associated with incidence, with non-linear effects of elevation and occupation distribution. Elevation, population density and climatic zone were significantly associated with envenoming incidence with a non-linear effect of elevation (S1 Appendix). According to the geostatistical model, the median predicted incidence for Sri Lanka was 397 per 100000 (IQR: 295 to 515 per 100 000) which is similar to the estimated national snakebite incidence derived from the survey. The fitted geostatistical model for snakebite incidence is summarised in table 4. There was a positive association between elevation and incidence up to 160 meters above sea level, with incidence dropping thereafter. Intermediate and wet climatic zones had higher snakebite incidence compared to the dry zone. Snakebite incidence decreased with increasing population density. Snakebite incidence rapidly increased as the proportion of the population engaged in agriculture increased up to 9%, with a more gradual increase thereafter. Table 5 shows the results of the geostatistical model for the incidence of envenoming. There was a positive association with elevation up to 195 meters above sea level, with incidence dropping at higher elevations. The incidence of envenoming was higher in the dry zone compared to intermediate and wet climatic zones and decreased with increasing population density. Estimated snakebite incidence and envenoming incidence maps for the whole of Sri Lanka are shown in Fig 2 (A) and 2 (B) respectively. The three estimated incidence maps shown in Fig 3 include contours demarcating incidence higher than the national snakebite incidence rate (Incidence>398 per 100000) and approximately +/- 100 per 100000 from the national incidence rate for snakebites (i. e. incidence 300 per 100000 and 500 per 100000). Corresponding PCMs based on exceedance probabilities of 300 per 100000,398 per 100000 (i. e. national rate) and 500 per 100000 are shown in Fig 4. The green areas of Fig 4 (A) indicate low probability areas for snakebites (cold spots); in this region the probability of observing a snakebite incidence more than 300 per 100000 is less than 0. 3. The red islands of Fig 4 (C) show the high probability areas for snakebites (hot spots), where the probability of having a snakebite incidence more than 500 per 100000 is greater than 0. 7. The estimated envenoming incidence maps are shown in Fig 5; contours are drawn to demarcate the national rate (151 per 100000), 50 per 100000 lower than the national rate and 100 per 100000 higher than the national rate. PCMs for exceedance thresholds of 100 per 100000,151 per 100000 (i. e. national rate) and 250 per 100000 are shown in Fig 6 (A), 6 (B) and 6 (C) respectively. The green areas of Fig 6 (A) indicate low probability areas for envenoming (cold spots), in this region; the probability of observing a snakebite incidence more than 100 per 100000 is less than 0. 3. The red islands of Fig 6 (C) shows the high probability areas for snakebites (hot spots); the probability of having a snakebite incidence more than 250 per 100 000 is more than 0. 7.
In a nation-wide community-based survey on snakebite we estimated there to be over 80000 bites, 30000 envenomings and 400 deaths per year in Sri Lanka. Our results are supported by a comprehensive methodology that included a representative, large, population-based sample and associated estimation techniques, and represents the most comprehensive survey of snakebite that has ever been undertaken in any country. Previous studies have attempted to estimate the national incidence of snakebite using hospital-based or local or regional (sub-national) data as highlighted in our estimate of the global burden of snakebite [1]. Such approaches are fraught with inaccuracies as hospital-based data are dependent on health seeking behaviour of the victims and generally underestimate incidence, while localized studies are almost always conducted in regions where snakebite is common, and tend to overestimate incidence [19]. Our previous work has shown that hospital data underestimate deaths due to snakebite by as much as 60% [4]. Incidence data may also depend on seasonal variations of bites. Our study demonstrates comprehensively the considerable geographical variation in snakebite incidence, the variation in the proportion of victims that are envenomed in different regions and the flaws in official data. Consequently, an accurate assessment of the burden of snakebite and envenoming can only be achieved by performing representative community surveys. Most of our results on the demography of bites are in agreement with the current literature. The highest rates of bites and envenoming were seen in the rural and agricultural north central and north eastern regions of the country. In keeping with other parts of the world, the highest disease burden due to venomous snakebite affects areas with the population groups that are most under-served in terms of healthcare and infrastructure. This re-emphasizes the need for equitable distribution of resources to address the problem of snakebite. Other ecological factors must also be taken into account to better understand these variations. Previous studies on snakebite show that rural males of working age are at high risk of snakebite due to high exposure levels associated with their lifestyle and occupation–mainly farming. In our study, males were also predominantly affected, and for both males and females the most susceptible were the economically active age groups. We have shown that the incidence of snakebite and envenoming gradually decreases with increasing population density. This may be due to more densely populated areas being less suitable as snake habitats. Snakebite incidence increased with increasing elevation up to 160 meters and thereafter dropped with increasing elevation, likely reflecting the facts that coastal areas are not generally suitable for snake habitats and that temperature declines with increasing altitude. Agricultural work has commonly been identified as a risk factor for snakebite. We assessed this by using the percentage of all workers undertaking agricultural work in each of the clusters. Snakebite incidence is dramatically lower in non-agricultural communities (less than 9% involved in agriculture) and snakebite incidence increased as the proportion of agricultural workers in a community increased. The incidence of snakebite is higher in both the intermediate and wet zones than in the dry climatic zone. However, the overall incidence of envenoming is higher in the dry zone. This is likely to reflect variation in species and again shows the importance of capturing data on the incidence of both snakebite and envenoming. Russell’s viper is the most widely distributed snake in Sri Lanka, and can be found up to an elevation of 1800 m. Saw-scaled vipers are largely confined to the arid dry zones of the country including Northern and Eastern Provinces extending up to eastern parts of the Southern Province. Cobras are also widely distributed in the country and are found up to 1500 m. Common Kraits are mainly found in the dry zone. Hump-nosed vipers are found all over the country, but mostly in the wet and intermediate zones [20][21][22][23]. Defining geographical areas by the level of risk of snakebite and envenoming is important for decision makers; allocation and distribution of antivenom, establishing snakebite management centres, and location of emergency treatment units and intensive care units all depend on assessing snakebite risk in a given location. Decision makers require thresholds to identify high risk areas for intervention. In this study, we drew cut-off points based on the national snakebite incidence rate for the country, which is about 400 per 100000. High risk areas or hotspots and low risk populations (cold spots) were defined by adding or subtracting 100 cases per 100000 to the national rate. Similarly, cut-off points were determined for envenoming bites based on the national rate of 150 cases per 100000. Hot and cold spots were determined by adding 100 cases per 100000 and subtracting 50 cases per 100000, respectively, from the national rate. Such estimates can be further refined to improve decision making by quantifying the likelihoods of exceeding or not exceeding any pre-defined threshold of incidence at a given location. Mapping of estimated incidence may lead to some uncertainty about whether the true incidence exceeds a given cut off threshold, as precision of the estimates varies between locations. Although the point estimate maps and PCMs are qualitatively similar, quantitative differences are important in interpreting the maps and decision making. Thus PCMs provide additional information to identify confidently high risk and low risk areas with respect to any pre-defined incidence threshold, as indicated on a PCM by probabilities close to one or zero respectively. PCMs also identify areas of high uncertainty, i. e. map areas with probabilities close to 0. 5, indicating a need for additional sampling to draw firm conclusions. The methodology used in the present study can be applied to estimate snakebite incidence in other geographical locations. In the methodology, the model captures the complexity of the snakebite process of a human, snake and environment interaction through a combination of covariate adjustments and spatially correlated residual variation, representing explained and unexplained spatial variation, respectively. This flexibility allows the approach to be applied in, with re-estimation of the model parameters to take account of difference in the available covariate information and consequent re-balancing between explained and unexplained spatial variation. There are inevitably limitations to this study. The quality of the household data depended upon memory of an event and is therefore subject to recall bias. However, the main outcome of interest, a snakebite, is highly memorable and we think it unlikely that this affects our estimates although the possibility of the incidence of envenoming being increased due to inappropriate classification of envenomed bites is acknowledged. Our sampling strategy was carefully designed to make it representative of Sri Lanka. However, the use of the electoral register to identify the first household in any cluster does mean that there is a chance that the small population in slums may have been under-sampled. As slums are relatively uncommon in Sri Lanka except in big cities, this is unlikely to have affected our estimates. Urbanization, which is highly correlated to several of the district-level variables, could also have been a confounder. Finally, cluster (i. e. GN) level income statistics were not available at the Census and Statistics Department of Sri Lanka. We therefore used district mean income to estimate cluster level income. This approach means that variability in income within a district may have led to some inaccuracies in the cluster level income estimates resulting in some residual confounding. In conclusion we have, for the first time, been able to provide robust community based incidence rates of snakebite and envenoming at a national level. We also developed estimated incidence maps and probability contour maps of within-country spatial variation, which can inform local healthcare decision making. The methodology we have used is replicable, as the required software is open-source. It also addresses the methodological limitations of the many epidemiological studies on snakebite that base their results on data from hospitals or from localized or time-limited surveys.
Snakebite is a neglected tropical disease which mainly affects the rural poor in tropical countries. There is little reliable data on snakebite, which makes it difficult to estimate the true disease burden. Hospital statistics underestimate numbers of snakebites because a significant proportion of victims in tropical countries seek traditional treatments. On the other hand, time limited or localized surveys may be inaccurate as they may underestimate or overestimate numbers depending on when and where they are performed. To get a truer picture of the situation in Sri Lanka, where snakebites are an important cause of hospital admission, we undertook an island-wide community survey to determine the number of bites, envenomings and deaths due to snakebite in the previous 12 months. We found that there were more than 80,000 bites, 30,000 envenomings and 400 deaths due to snakebite, much more than claimed by official statistics. There was variation in numbers of bites and envenomings in different parts of the country and, using the data from our survey, we were able develop snakebite risk maps to identify snakebite hotspots and cold spots in the country. These maps would be useful for healthcare decision makers to allocate resources to manage snakebite in the country. We used free and open source software and replicable methods, which we believe can be adopted to other regions where snakebite is a public health problem.
cognitive science
medicine and health sciences
pathology and laboratory medicine
decision making
tropical diseases
geographical locations
vertebrates
neuroscience
animals
epidemiological methods and statistics
signs and symptoms
reptiles
cognition
neglected tropical diseases
population biology
snakebite
public and occupational health
sri lanka
epidemiology
epidemiological statistics
people and places
snakes
population metrics
diagnostic medicine
asia
squamates
hemorrhage
biology and life sciences
population density
vascular medicine
amniotes
organisms
Mapping the Risk of Snakebite in Sri Lanka - A National Survey with Geospatial Analysis