Trends in snakebite deaths in India from 2000 to 2019 in a nationally representative mortality study
Abstract
The World Health Organization call to halve global snakebite deaths by 2030 will require substantial progress in India. We analyzed 2833 snakebite deaths from 611,483 verbal autopsies in the nationally representative Indian Million Death Study from 2001 to 2014, and conducted a systematic literature review from 2000 to 2019 covering 87,590 snakebites. We estimate that India had 1.2 million snakebite deaths (average 58,000/year) from 2000 to 2019. Nearly half occurred at ages 30–69 years and over a quarter in children < 15 years. Most occurred at home in the rural areas. About 70% occurred in eight higher burden states and half during the rainy season and at low altitude. The risk of an Indian dying from snakebite before age 70 is about 1 in 250, but notably higher in some areas. More crudely, we estimate 1.11–1.77 million bites in 2015, of which 70% showed symptoms of envenomation. Prevention and treatment strategies might substantially reduce snakebite mortality in India.
Introduction
The World Health Organization (WHO) estimates that 81,000–138,000 people die each year from snakebites worldwide, and about three times that number survive and but are left with amputations and permanent disabilities (World Health Organization (WHO), 2019a). Bites by venomous snakes can cause acute medical emergencies involving shock, paralysis, hemorrhage, acute kidney injury and severe local tissue destruction that can prove fatal or lead to permanent disability if left untreated. Most deaths and serious consequences from snakebite envenomation (exposure to venom toxins from the bite) are avoidable by timely access to safe and effective antivenoms (Gutiérrez et al., 2017). Snakebite deaths and envenomation are largely neglected topics in global health. However, in 2017, the WHO included snakebite envenoming in the priority list of neglected tropical diseases (World Health Organization (WHO), 2019b) and launched in 2019 a strategy for prevention and control of snakebite, aiming to halve the numbers of deaths and cases of serious disability by 2030 as compared to 2015 baseline (World Health Organization (WHO), 2019c). Achieving this goal will require substantial progress in India, which is home to approximately half of global snakebite deaths. Snakebite deaths and disability remain a major public health challenge also for poor rural communities in many parts of Asia, Africa, Latin America and Oceania.
Direct estimation of 46,000 annual snakebite deaths in India in 2005 (Mohapatra et al., 2011) prompted a revision of the WHO’s global total, which had estimated about that number for the entire world. The 2005 Indian estimate relied upon analyses of about 123,000 verbal autopsy records from 2001 to 2003 in the Registrar General of India’s (RGI) Million Death Study (MDS), one of the largest nationally representative mortality surveys. Now the MDS has reported cause-specific mortality patterns on over 600,000 deaths from 2001 to 2014 for the whole of India. Here, we report seasonal and temporal trends in snakebite mortality over the last two decades in India and its spatial distribution. We provide estimates of total snakebite deaths for the 20-year period 2000–2019 by age and sex. Our earlier report estimated a crude ratio of about one death to 20 envenomations. We now further quantify the levels of envenomations based on a systematic review of 88,000 snakebites in the published literature. The literature also provides details on the specific causes, bite locations, and treatment of envenomations. Finally, enhanced surveillance including facility-based tracking will be central to the Government of India’s strategies to reduce snakebite deaths. Thus, we provide estimates on the degree to which snakebites and deaths are reported adequately in public facilities. Appendix 1—figure 1 shows the overall study design, data sources, input resources and outcomes.
Results
Trends in snakebite mortality and its geographic and temporal patterns
From 2001 to 2014, the MDS reported deaths with causes classified by physicians who examined verbal autopsy records collected from over 3.6 million households in three distinct nationally representative sampling frames (1993–2003; 2004–13; and 2014–23). Two of 404 independent physicians coded each death to the International Classification of Diseases-10th revision (ICD-10), reconciling (anonymously) any coding differences with a senior physician adjudicating any persistent disagreements (Gomes et al., 2017; Aleksandrowicz et al., 2014; Menon et al., 2019). Among 611,483 available records, 2833 deaths were assigned to snakebites (ICD-10 code X20). The two physicians agreed on the diagnosis 92% of the time. About 94% of snakebite deaths occurred in rural areas, and 77% occurred out of hospital (Appendix 1—table 1).
We applied the age- and sex-specific proportion of snakebite deaths to total deaths as estimated by the United Nations Population Division (UN) for India (United Nations, 2019) to estimate national death rates by age and sex, as well as absolute totals for each year (Table 1). The UN totals are based on careful demographic review of census and other data sources. The fieldwork procedures of the Sample Registration System (SRS the underlying demographic survey on which the MDS is based) leads to some undercounts (of about 5–10%) of expected deaths (Gerland, 2014). The SRS is representative at the state and rural/urban strata, and has a large, distributed sampling covering over 7000 small areas in the whole of the country (Registrar General of India, 2017). Hence, any missing deaths are generally randomly distributed across states, and not clustered in one state or one key sub-group, such as in rural areas (Dhingra et al., 2010; Aleksandrowicz et al., 2014; Menon et al., 2019). Thus, the proportion of snakebite deaths is not likely an underestimate. However, total snakebite deaths might be underestimated. The use of the UN death totals adjusts for these possible undercounts and provides a plausible national total for each year.
Total snakebite deaths in India from 2001 to 2014 totaled about 808,000, with reasonably narrow uncertainty range of 738,000 to 833,000, based on both physicians immediately assigning snakebites or only one physician doing so. Some age-specific death rates fell, but as population growth averaged 1.1% annually, the application of annual age-specific rates to the UN death totals for that year showed that the overall number of snakebite deaths grew from about 55,000 in 2001 to about 61,000 in 2014. During the 2001–2014 MDS study period, the average age-standardized snakebite death rate (using the Indian census population of 2001 to take into account the minor change in age structure) was 4.8 per 100,000 population, falling annually by 0.8%.
Declines in the age-specific snakebite death rate were fastest for children aged 0–14 years (declining by about 1.6% annually), with slower declines in young adults aged 15–29 years (1.2% annually) and no declines among middle-aged adults (30–69 years). Before 2010, snakebite death rates were higher in boys than girls but from 2010 to 2014, death rates in girls exceeded those for boys (Appendix 1—table 1). The age-specific risks translate to a probability of 0.37% (uncertainty range 0.34–0.38%) of dying from snakebite before age 70 years in the absence of competing mortality (Table 1). This suggests that the average risk of an Indian dying from snakebite prematurely before age 70 is approximately 1 in 250.
Because the risk of dying from snakebites has been stable from 2001 to 2014, we can make reasonably reliable forward projections from 2015 to 2019 and backward projections from 2001 to 2000 (Table 2). This reveals that 1.2 million snakebite deaths occurred over this 20-year period. Of these deaths, 602,000 occurred among males and 565,000 occurred among females. With both sexes combined, about 543,000 (47%) occurred in middle-age (30–69 years), 325,000 (28%) among children below 15 years, 197,000 (17%) among adults aged 15–29 years, and 102,000 (9%) among those over age 70 years. Using the agreement of one or two physicians on the cause yielded generally narrow uncertainty estimates for each sex and age groups.
From 2001 to 2014, just under 70% of these snakebite deaths occurred in eight states with about 55% of the population: Bihar, Jharkhand, Madhya Pradesh, Odisha, Uttar Pradesh, Andhra Pradesh (which includes Telangana, a recently defined state), Rajasthan and Gujarat (Table 3). In these high-burden states, the age-standardized death rate was about six per 100,000. Snakebite death rates generally rose over time in most high-burden states, particularly in Bihar, but fell in Andhra Pradesh. The remaining lower-burden states began the study period with age-standardized death rates of about 3.7, which fell over time. Figure 1 shows the absolute risk of dying from snakebite using data from 7400 small areas (the small sampling units used in the RGI’s Sample Registration System for the MDS) from 2004 to 2013. The absolute risks were calculated applying spatially smoothed predictive relative risks from a spatial Poisson model to the overall national risk before age 70 years (of about 0.4%, Table 1) after adjusting for any differences in rural/urban status, female illiteracy levels, temperature, and altitude of local areas. We observed greater than 0.6% (1 in 167) mortality risk before age 70 years in the highest risk sub-areas of Andhra Pradesh, Odisha, Bihar, Uttar Pradesh, Madhya Pradesh, Chhattisgarh, and Rajasthan. About 260 million people lived in these areas in 2015, including about 4 million people living in hot spots that had a 1% or greater risk of death from snakebite. Appendix 3 provides statistical details and credible intervals of these risk estimates.
Half of all snakebite deaths occurred during the southwest monsoon seasons from June to September. Seasonality was similar in each of the study years (Appendix 1—figure 2) and was similar in higher-burden and lower-burden states (data not shown). We used a Poisson time series model for snakebite deaths from 2001 to 2014 to predict the average daily snakebite mortality in India. The peak (294 deaths per day) was in mid-July and the trough (78 deaths per day) was in mid-February (Figure 2A). We also aggregated the deaths from 2001 to 2013 by every 100 m of altitude above sea level (Figure 2B). The crude death rates in areas below 400 m were about three times those in areas at about 1000 m. Over 80% of snakebite deaths occurred below 400 m and 50% occurred below 200 m.
Characteristics of snakebites from a systematic review of the literature
A systematic literature review yielded 87,590 snakebite cases in India (both fatal and non-fatal) from 2000 to 2019 based on screening 1417 papers and including 78 studies from 24 states or union territories in India (Figure 3, Appendix 2—figure 1, Appendix 2—table 2). In the published studies, snakebites were more common in males (59%), at ages 30–69 years (57%), from June to September (48%), and occurring outdoors (64%). These results match the relevant results for the MDS. However, MDS snakebite deaths were equal between males and females. The leg was the dominant site of bite (77%), and the time of reported bite was throughout the day. Of the treated cases, nearly two-thirds (66%) were seen within 1–6 hr, with the remainder seen after six hours. The proportion treated within 1–6 hr improved over time (data not shown). In the fewer studies that attempted to identify the snake species, Russell’s viper (Daboia russelii) constituted 43%, followed by unknown species (21%), krait (Bungarus species) (18%), and cobra (Naja species) (12%).
Snakebite and mortality surveillance
The Government of India relies on reporting via public hospitals to track snakebites and deaths (Government of India, 2015). We examined the total bites and deaths available from 2003 to 2015 in Government hospitals and compared these deaths to the MDS in-hospital deaths (Table 4). Over this 13-year period, the MDS estimated about 154,000 snakebite deaths in public and private hospitals, and the Government reported 15,500 deaths in hospitals, meaning that the routine reporting system captured only 10% of the expected hospital-based deaths. The most complete reporting was in Karnataka which captured 26% of expected hospital snakebite deaths.
Snakebite prevalence and envenomation
Among the 87,590 snakebites reported in the literature, there were 3329 reported deaths (Appendix 2—table 1). We fitted death and bite data from each study to an ordinary least square regression to calculate a case-fatality rate, after removing the extreme outliers. We estimated a crude case-fatality rate of 3.2% for in-hospital cases. Based on mostly cautious assumptions about the ratio of in-hospital to out-of-hospital prevalence of snakebites (Appendix 1—table 2), we estimate the total number of snakebites to range from 1.11 to 1.77 million in 2015. Based on 44 hospital studies where 70% of patients sought treatment, were diagnosed with systemic envenomation, and received antivenom, we estimate that the annual number of envenomations is about 0.77 to 1.24 million with the remainder being ‘dry bites’ or bites by non-venomous species (0.33 to 0.53 million).
Discussion
Our nationally representative mortality study documents about 1.2 million snakebite deaths from 2000 to 2019. Most occurred at home in the rural areas. About 70% occurred in eight higher-burden states and half occurred during the rainy season and in low altitude rural areas. While rates of childhood and young adult snakebite mortality have fallen, those in middle age have not. Thus, the average risk of an Indian dying from snakebite before age 70 is approximately 1 in 250, but in some areas, this risk approaches 1 in 100. Over 260 million Indians live in areas of moderate risk of about 1 in 167. More crudely, approximately 1.11–1.77 million bites occur annually with about 70% representing envenomation, and 58,000 dying. While snakebite deaths represent only about 0.5% of the approximately 10 million deaths that occurred in India in 2015, they are nonetheless important, as they are nearly all avoidable.
Many of the features of snakebites and deaths were known or suspected, but few were quantified reliably (Mohapatra et al., 2011). Our study’s novel contributions are to quantify some of these features, and identify other findings that are relevant to improved epidemiological understanding and to prevention and treatment in snakebite control programs. The map of the snakebite mortality risk (Figure 1) highlights ‘hot spots’ in each state, which are at lower altitude. This reflects not only the more highly populated and the more extensive and intensively farmed arable land at lower altitudes, but also the species and population densities of snake species of medical importance. These snake densities are sometimes very high, particularly in grain agriculture which attracts the largest rodent and amphibian populations that are eaten by snakes (Whitaker and Captain, 2004; Mise et al., 2016). Focusing on agrarian communities in specific areas which carry the highest risk of mortality, especially during the monsoon seasons, could reduce mortality and morbidity attributable to snakebites. Targeting these areas with education about simple methods, such as ‘snake-safe’ harvest practices, wearing rubber boots and gloves and using rechargeable torches (or mobile phone flashlights) could reduce the risk of snakebites. Mass distribution of mosquito nets (which also protect against scorpion sting and mosquito-borne diseases) is a relevant strategy that could build upon the National Vector Borne Disease Control Program’s efforts to control malaria, kala-azar, and arboviral infections.
Our study has implications for better treatment, particularly in the distribution of effective antivenom to the areas and populations in greatest need. Increased use of antivenom would require tactful cooperation with local traditional healers and ayurvedic practitioners to persuade them to refer severely ill patients for treatment with antivenom, and raising awareness of the effectiveness of antivenom. Government hospitals can make antivenom freely available to snakebite victims (Whitaker and Whitaker, 2012). Health services could monitor adverse reactions to antivenom and improve distribution and cold-chain storage, matching supply to places and times of greatest need. Training of local medical staff and emergency responders should be improved so that they can administer antivenom by intravenous injection and also identify and treat early anaphylactic reactions. India has sufficient manufacturing capacity to make large amounts of snake antivenom. Better understanding of the distribution of India’s many venomous snake species could help in the development of more appropriate antivenoms. The current Indian polyvalent antivenoms neutralize venom from only common cobra (Naja naja) (there are three other Indian cobra species), common krait (Bungarus caeruleus) (seven other krait species), Russell’s viper (Daboia russelii) and saw-scaled viper (Echis carinatus) (World Health Organization (WHO), 2010; Warrell, 2011; Warrell et al., 2013, Whitaker and Martin, 2015; Senji Laxme et al., 2019). The few studies from healthcare facilities found that antivenom treatment reduced deaths by over 90% (Appendix 2—table 1).
We estimate, crudely, that the in-hospital case-fatality rate based on the literature was about 3%. This in part reflects delay in reaching medical care, with only about half of the cases doing so within 6 hr. Public-private partnerships for ambulance services are possible. In some states, an emergency ambulance service equips vehicles with lifesaving equipment and drugs, including antivenom. Ambulances can be summoned in 15 states of India by calling a toll free number (Gimkala et al., 2016). In 2014, of 27,509 snakebite patients transported to hospitals within 6 hr, 359 patients died within 48 hr of follow-up. This represents a crude case-fatality rate of 1.3%, below the rate we estimate for in-hospital bites. This ambulance model is relevant to other parts of the country, especially the more remote areas in Bihar, Jharkhand and Odisha, which are not currently covered.
Finally, improved surveillance is required of venomous snake species as well as the human consequences of bites. An enhanced snake species database, hosted in collaboration with agricultural and forest departments, of their habitat details, clear photographs, and geographical distributions is now available as a downloadable Google app (Indian Snakes, 2019; World Health Organization (WHO), 2019d). There are at least 15 species of snake in India responsible for human deaths, and better information about them would aid control (Whitaker and Whitaker, 2006; Whitaker and Martin, 2015). We show that public facility-based reporting of deaths captures only 10% of expected deaths in public and private hospitals. While much care in India occurs in private hospitals, disease reporting and surveillance by private facilities is likely to be similarly or even more deficient (Jha and Laxminarayan, 2009). The Government of India could designate and enforce snakebite as a ‘Notifiable Disease’ within the Integrated Disease Surveillance Program. However, since most deaths occur at home, community death tracking through ongoing mortality surveillance will be needed. Both community and facility-based surveillance data are essential to better align interventions to prevent and treat snakebites to their heterogeneously distributed burden. For both hospital and community bites, it would be invaluable to document the circumstances and consequences of snakebites, including morbidity sequelae, with simple questions to investigate the circumstances of each bite or death/disability. This could include data such as use of boots, walking in the dark, sleeping patterns, and other questions. Similar investigative tools for HIV/AIDS have recently been proposed to add to the WHO’s standard verbal autopsy (Bogoch et al., 2018).
Limitations of our study
The major source of uncertainty in our estimates of snakebite deaths at national level arises not from random errors, as the MDS has a large sample size and the vital rates used as underlying denominators are reasonably complete (Menon et al., 2019), but from the misclassification of causes of death in the verbal autopsy. Earlier evaluations of the MDS showed strong reproducibility of the dual physician-coded verbal autopsies, generally low rates of misclassification in children and young and middle-age adults, and high consistency with relevant hospital or clinical data (Gomes et al., 2017; Aleksandrowicz et al., 2014; Menon et al., 2019). Moreover, two independent physicians agreed about 92% of the time on a diagnosis of snakebite deaths. Snakebite mortality may be under-estimated because the phenomenon of painless, unsuspected nocturnal krait (Bungarus) bites resulting in ‘early morning paralysis’ may not be attributed to snakebite (Saini et al., 1986; Ariaratnam et al., 2008), but this is likely to be a small bias. Our estimates for some states are uncertain due to small number of deaths recorded annually, which also prevented us from examining yearly spatial clustering patterns.
Our estimates of case-fatality rate and envenomations based on the systematic literature review has obvious limitations. First, the exact species of snake cannot be easily identified, and indeed, there are four species of cobra and eight species of kraits in the country (Whitaker and Captain, 2004). In addition, each species varies in the circumstances, seasonal and diurnal variation and types of terrain where bites most often occur. For example, anecdotal experience indicate that most bites from common kraits (Bungarus caeruleus) occur at night while people are sleeping on mats on the floor or ground, in or near home, unprotected by tucked-in mosquito nets (Ariaratnam et al., 2008; Kularatne, 2002; Bawaskar and Bawaskar, 2002). Most bites from saw-scaled vipers (Echis carinatus) happen either when the snake is stepped upon in bare or sandaled feet at night or when cutting grass by hand with a short sickle. Bites by cobras are divided into circumstances such as defensive bites while stepping on them during planting/harvesting crops, reaching into piles of straw or firewood and predatory bites when the cobra mistakes a human hand or foot for a prey item (Alirol et al., 2010). Russell's viper (D. russelii) bites occur during the day, inflicted on farmers in the paddy fields or while hand harvesting peanut plants, or at night when someone walks without using a light and steps on the snake (Whitaker and Captain, 2004).
Another major uncertainty from the literature review are in the hospital-based data, given expected problems with publication biases and in the differences between patients seeking or not seeking hospital care. For example, serious cases are more likely to be hospitalized, raising the observed case-fatality rate. However, under various scenarios, our report of 20 to 40 bites per death is greater than the crude estimate of 20 envenomations per death we made earlier (Mohapatra et al., 2011). Additional hospital-based surveillance, including tracking the severity of treated cases, could further refine the actual ratio of envenomations to deaths. These uncertainties demand appropriate caution in interpreting our basic estimate of the number of envenomations.
Conclusion
We conclude that snakebite deaths in India are concentrated largely within limited geographical areas, and involve particular communities during specific seasons. Our identification of the focused geographic and temporal spread of snakebites allows targeted prevention and treatment strategies that could help India to achieve the WHO’s goal of halving snakebite death and morbidity rates by 2030. Further use of nation-wide, representative epidemiological studies will be essential to review the success of such control programs.
Materials and methods
Data sources
Request a detailed protocolTo derive comprehensive and up-to-date estimates of snakebite mortality and prevalence, we collected all possible statistics related to snakebites in India from 2000 to 2015. The main data sources for this study were snakebite mortality data from the Indian Million Death Study (MDS), a systematic review of studies published in the scientific literature, and chronological statistics published by the Ministry of Health and Family Welfare of the Government of India.
Nationally representative mortality data
Request a detailed protocolThe methods, strengths, and limitations of the MDS and key results for various diseases have been extensively reviewed and published (Aleksandrowicz et al., 2014; Gomes et al., 2017; Menon et al., 2019). Briefly, in collaboration with the Registrar General of India, the MDS monitored approximately 23 million people in 3.6 million nationally representative households in India from 1998 to 2014. The Registrar General of India’s Sample Registration System (SRS) established three sampling frames for the MDS, which covered years 1993-2003, 2004-2013, and 2014-2023. The SRS randomly selects sampling units based on the 1991, 2001, and 2011 censuses for the respective sampling frames (Registrar General of India, 2017). Mortality data used in this study were from 2001-2003, 2004-2013, and 2014, generated from these sampling frames. Every six months, about 900 non-medical surveyors recorded the details of each death that occurred in these households during the preceding six months using a well-validated verbal autopsy instrument (based on the 2012 WHO instrument and including a half-page local language narrative). Each record is converted to an electronic form and randomly assigned to two of 404 trained physicians, who each assign a cause of death using ICD-10 codes. Disagreements in assignment undergo anonymous reconciliation, and persisting differences undergo adjudication by a third physician. We included 2833 snakebite deaths in our study by carefully examining 3020 probable snakebite deaths that either of the two physicians had coded as X20 (venomous snakes), X27 (venomous animals) or X29 (not specific). We followed the same inclusion/exclusion method described in our earlier analyses (Mohapatra et al., 2011). Out of the 3020 possible snakebite deaths, there were 2779 (92%) deaths in which both coders initially coded to X20. Review of these yielded no misclassified deaths. Re-examination of the symptoms and physician keywords for 105 deaths that one coder had coded as X20 and other coder as X27 or X29 revealed that 54 (2%) were snakebite deaths. No misclassified deaths were found in the 136 deaths that one coder had coded X20, X27 or X29 and the other coder had assigned a different ICD-10 code.
Statistical methods
Geospatial mapping of deaths
Request a detailed protocolThe SRS provided population data for the sampling units for 2004–2013 (Registrar General of India, 2017). These population values were partitioned into single-year ages by applying 2011 Census Registrar General of India, 2011 district-level single-year age structure proportions (by sex and rural/urban setting). District codes from 2011 Census (Registrar General of India, 2011) were converted to 2001 codes prior to linking with sampling units given that the 2004–2013 sampling frame used 2001 codes. Sampling units that belonged to districts that split in 2011 used 2001 district codes from the parent districts. All sampling units that belonged to the same district, rural/urban setting and sex shared the same age structure proportions. We further linked the MDS data to these population data at the sampling unit level. Statistical analyses were based on data from 7377 geocoded sampling units (out of 7597 sampling units), after exclusion of sampling units from the islands.
We derived the spatially-smoothed absolute risks of snakebite mortality in India for 2004–2013. First, using mortality data of ages 0–69 years (by 5-year age group) as the outcome and sampling unit population of the same age range as the offset, we fitted a Bayesian Poisson model to obtain the age- and sex-specific snakebite death rates at the national level. We did not include an intercept in the model, but included age-sex interaction term and time trend (using year 2010 as the reference value) as covariates. This formulation allowed us to obtain the estimated age- and sex-specific national death rates for year 2010. We then used a geostatistical Bayesian Poisson model to estimate the spatially smoothed relative risks of snakebite mortality, by comparing the observed snakebite death rates at each sampling unit versus the national death rate (see Appendix 3 for details). The geostatistical models adjusted for time trends, urban/rural status, female illiteracy in rural areas, altitude, and average of long-term monthly mean temperature. We adjusted for urban/rural status of the sampling unit due to the higher risks of snakebite in rural areas compared to urban areas (Chaves et al., 2015). We included female illiteracy in rural areas as a proxy of poverty effects on snakebite mortality, since the poor have higher risk of snakebite (Harrison et al., 2009). We used sub-district-level female illiteracy data from the 2011 Indian census (Registrar General of India, 2011). We included altitude and long-term monthly temperature as covariates since they affect the occurrence of snakebite (Chaves et al., 2015). Altitude data came from the NASA Shuttle Radar Topographic Mission’s digital elevation data version 4 (Jarvis et al., 2008). Long-term monthly mean temperature came from the University of Delaware’s air temperature gridded dataset V5.01, with a 0.5 degree latitude/longitude grid resolution (Willmott and Matsuura, 2001). We also included spatial random effects and sampling unit-level random effects in the model. Thus, the spatially smoothed relative risks were the predicted relative risks of snakebite mortality across India. These relative risks were assigned to grid cells that covered the country (see Appendix 3: statistical supplement). Finally, we calculated the absolute risks of snakebite mortality across India by multiplying the spatially-smoothed relative risks (in each grid cell) by the national risks of dying before age 70 years from Table 1. National risks of dying used the average of annual risks of dying for 2004–13. The absolute risks represent the risk of dying from snakebite before age 70 years at the grid cell location. Population estimates in high-risk areas were obtained by overlaying the absolute risk surface on the Gridded Population of the World version 4 for year 2015 (Center for International Earth Science Information Network - CIESIN - Columbia University, 2015). Further technical explanation on the geostatistical Bayesian model is published (Brown, 2015). Appendix 3 provides the model form, equations and implementation.
Systematic literature review
Request a detailed protocolWe performed a systematic review of snakebite studies in India. We searched the literature using a combination of keywords related to the study setting, metrics, treatment, snake species, and geography in Ovid MEDLINE(R), PubMed, Web of Science, and Scopus electronic databases. We selected relevant studies published in the English language from January 1 2000 to September 1 2019, to collect data for understanding case-fatality patterns and important snakebite characteristics in India. In addition, we hand searched articles that had cited our 2011 publication (Mohapatra et al., 2011). Appendix 2 provides the keywords and inclusion and exclusion criteria. The search initially found 1417 snakebite mortality and morbidity studies. After a careful review of titles, abstracts and quality of the studies by three independent reviewers (MB, KP and WS), 78 of the 95 possible studies were included in our analysis. We categorized the 78 studies to four study types: autopsy (seven studies: only deaths), emergency medical services (EMS) (one study of prevalence), hospital (66 studies: both prevalence and deaths) and community (four studies: both prevalence and deaths).
Mortality rates
Request a detailed protocolWe applied the SRS probability of selection sampling weights to the snakebite death frequencies to address urban and rural differences. We calculated snakebite mortality fractions using three-year backward moving averages of weighted snakebite death frequencies for each age, sex, and urban/rural stratum for each state of India. We interpolated the mortality fractions, using standard statistical methods for strata with zero death count (SAS Institute, 2014). We applied these mortality fractions to SRS and the India census demographic framework to obtain the snakebite death rates. We then adjusted the death rates (usually upward by slight amounts) to the United Nations Population Division (United Nations, 2019) estimated India death totals (United Nations, 2019) to obtain the numbers of national and sub-national snakebite deaths. To address the remaining noise from crude death rates, we fitted cubic spline regressions to 2003 to 2014 cause-specific death rates while adjusting snakebite mortality to other causes of death to obtain the final estimates. We obtained estimates for the beginning and end of the period, including years 2001 and 2002 where data used for moving averages were less than three years, by extrapolating the spline curves to cover the overall period of 2000–2015. For comparison of rates across the years, we standardized the death rates to the 2001 census population. We calculated the number of in-hospital and out-of-hospital deaths by multiplying estimated deaths by percentages of study deaths that occurred in-hospital and out-of-hospital, as reported in the MDS. Appendix 1—figure 1 shows the data sources, data inputs and outcomes.
Snakebite prevalence estimates
Request a detailed protocolWe used an indirect method to estimate snakebite prevalence (which given very short duration of each bite effectively represents incidence) measured in terms of in-hospital and out-of-hospital prevalences. This involved using the estimated MDS hospital deaths divided by case-fatality rate from systematic review of the literature to estimate the in-hospital prevalence and apply a hypothetical relationship between in-hospital and out-of-hospital prevalence to estimate the out-of-hospital prevalence. This is, by necessity, crude but provides some reasonable ranges to estimate the numbers of bites and envenomations in India in recent years. We applied all such prevalence estimates for 2015, as that was the closest year to the last MDS round of 2014. Appendix 1—table 2 provides details of the calculation and the assumptions.
Appendix 1
The role of data sources and snakebite mortality/prevalence estimates for India.
Appendix 2
Snakebite systematic review search strategy and selection criteria of studies
Search strategy
We searched Ovid MEDLINE(R), PubMed, Web of Science and Scopus databases for relevant studies from their date of inception from January 1, 2000 to September 1, 2019.
List of Keyword Categories
Keyword category | ||||
---|---|---|---|---|
Setting | Metrics | Treatment | Snake / Species | Geography / Location |
Hospital | Death | Antivenom | Venomous | India |
Survey | Mortality | Treatment | Snakebites | not *Indians, North American |
Urban or Rural | Incidence | Treated | Bites or Stings | List (of Indian States and Union Territories) |
Community | Prevalence | Untreated | Envenomation | |
Cases | List (of India venomous snakes from www.indiansnakes.org) | |||
Case fatality | ||||
Mortality rate | ||||
Prevalence rate | ||||
Incidence rate |
Inclusion Criteria
We included studies conducted in India that described:
Snakebite cases, both fatal and non-fatal, either exclusively or as a subset of acute poisoning cases,
Demographic distribution, complications, treatment and outcomes of the snake bite cases.
Exclusion Criteria
We excluded studies that were:
Secondary reviews,
Animal or pharmacological studies on therapy drugs or herbal therapies,
Injury studies on bites and stings by animals other than snake,
Studies beyond the Indian context,
Biomedical studies that described only the biological mechanism and manifestations of a snakebite,
Poor quality, failed at the quality assessment.
List of articles included in systematic review:
Autopsy studies
Community studies
Emergency Medical Services (EMS) studies
Hospital studies
Appendix 3
Statistical supplement to snakebite analysis
Introduction
Here, we provide further details on distribution of MDS sampling units, the Bayesian statistical methods, and model validations of the estimated absolute risk of the study.
The sampling units
The India’s Sample Registration System (SRS) is based on a system of dual recording of births and deaths in fairly representative sample units spread all over the country. In the 2004–13 sampling frame, there were 7597 sampling units, and 7416 of them were geocoded. Excluding the islands, there were 7377 geocoded sampling units (Appendix 3—figure 1).
Statistical models
The statistical models were fitted using Bayesian inference, and posterior distributions were computed using the INLA methodology from inla Rue et al., 2009 as explained in the following sections.
Spatially smoothed relative risks of snakebite mortality
A Poisson regression model with a spatial random effect was used for the spatial analyses, a model which is a form of Generalized Linear Geostatistical Model described by Diggle and Ribeiro, 2006. The approach is broadly similar to the methods used by Bhatt et al., 2015 for estimating malaria prevalence in Africa, and by Jiang et al., 2014 for assessing the contribution to cancer incidence of ambient radiation near a nuclear generating facility in Canada.
Writing as the outcome in sampling unit (located at the spatial coordinates ) at time , the model is as follows.
is the expected number of snakebite death (described below). The (log-transformed) relative risk depends on the spatially referenced covariates , a non-linear time trend f, the spatial random effect , and sampling unit-level effect .
The spatial explanatory variable includes an intercept, a rural-urban indicator variable, the sub-district-level proportion of female illiteracy in rural areas from the 2011 census, altitude, and long-term monthly mean temperature. We applied change points for altitude at 400 meters and long-term monthly mean temperature at 20°C based upon exploratory analysis of the relationship between these covariates and the relative risks of snakebite mortality.
The time trend is a non-parametric effect (a wiggly line) modelled as a second order random walk, with year 2010 as the reference value.
The spatial random effect is a Gaussian random field with a Matern spatial correlation function specifying how correlation decreases with distance , depending on the value of a range parameter . The range parameter determines how rough or smooth the relative risk surface is.
The unit-level effect is spatially independent, and unlike the spatial effect , two sampling units in close proximity will have unrelated values of . The can be thought of as accounting for short-scale spatial variation or sampling-unit-level risk factors not included in the model as covariates.
The absolute risks mapped in main text Figure 1 is the posterior median of the spatial relative risk in 2010 (since f2010 = 0) multiplied by the national value for risk of dying before age 70 years.
The expected counts were obtained by computing age-sex-specific rates (by 5 year age group) using a non-spatial Poisson regression model with the age-sex specific death counts as the response, the population of the age-sex group in the study as an offset, and a linear time trend (with reference at year 2010). The estimated rates were multiplied by each sampling unit’s age-sex population in each year and summed over groups to produce each unit’s yearly expected count.
Implementation
As there are 7377 geocoded sampling units, and 6587 different spatial locations , model fitting is computationally intensive and an approximation to the spatial covariance matrix is necessary. The Markov random field approximation from Lindgren et al., 2011 is used here, and implemented in the geostatsp package (Brown, 2015) for the R statistical programming language (R Development Core Team, 2018), which in turns calls the inla software (Rue et al., 2009). A full description of the methodology is published Brown, 2016. Although there are other methods available for fitting models of this type, the task is complex and computationally demanding and there is currently no rival to the Bayesian methodology in the inla software for fitting spatial models with non-Normal responses.
Bayesian inference requires specifying prior distributions, and spatial models are particularly susceptible to producing spurious results from ill-chosen priors for the spatial parameters and . Here the penalized complexity prior distributions from Simpson et al., 2017 are used, priors which discourage a spatial effect (wanting flat and close to zero) unless the data indicate a clear preference for a spatial model. Following Simpson et al., 2017, the standard deviations and have exponential priors, as does the scale parameter . The prior median for , the distance beyond which the correlation between two locations is under 10%, is or roughly one sixth of the distance across India. The prior medians of and are both , a value at which a one standard deviation increase in or doubles mortality risk.
Model validation
Snakebite crude death rates
A map of the crude death rates of snakebite at the sampling units is provided in Appendix 3—figure 2. Overall, areas with higher crude death rates correspond to areas with higher absolute snakebite risks (in Main text: Figure 1). It should be noted that the crude death rates have not accounted for the other factors that were included as covariates in the spatial Poisson regression model.
Uncertainty of absolute snakebite risks
Uncertainty of the absolute snakebite risks can be depicted by the 95% credible interval of the posterior distribution for each grid cell (Appendix 3—figure 3a-b). For comparison, the estimated median absolute snakebite risks is also provided (Appendix 3—figure 3c).
Covariate parameter estimates
Here, we provide the parameter estimates for covariates included in the Bayesian spatial Poisson model (Appendix 3—table 1). Urban/rural status, female illiteracy in rural areas, altitude above the change point of 400 m, and temperature below the change point of 20°C showed statistically significant association with snakebite mortality risk. Female illiteracy in rural areas and increased temperature below 20°C were associated with higher risks of snakebite mortality, and vice versa for urban status and increased altitude at attitude level above 400 m. For time trend, snakebite mortality risks showed a decreasing trend over time (in terms of year, Appendix 3—figure 3).
Data availability
Data from the Million Death Study (MDS) India cannot be redistributed outside of the Centre for Global Health Research due to legal agreement with the Registrar General of India. However, access to MDS data can be granted via data transfer agreements, upon request to the Office of the Registrar General, RK Puram, New Delhi, India (rgi[dot]rgi[at]nic[dot]in). The public census reports and Sample Registration data can be accessed thorough http://censusindia.gov.in/. Source data for Figure 1 is already explained in Appendix 3. Source data files have been provided for Figure 2 and for Appendix 1-figure 2. The systematic review data used for Figure 3 is already presented in Appendix 2-table 1.
References
-
Snakebite envenomation in India: a rural medical emergencyIndian Pediatric 43:553–554.
-
Snake bites profile in Haldia, Purba Midnapur, West BengalIndian Journal of Forensic Medicine and Toxicology 7:243–245.
-
Retrospective analysis of snake victims in Northern India admitted in a tertiary level instituteJournal of Anaesthesiology Clinical Pharmacology 28:45.https://doi.org/10.4103/0970-9185.92434
-
Study of seasonal variation in snake bite patients admitted to a Mangalore based tertiary care hospitalIndian Journal of Critical Care Medicine 18:S37.
-
Snake bite in South Asia: a reviewPLOS Neglected Tropical Diseases 4:e603.https://doi.org/10.1371/journal.pntd.0000603
-
Role of neostigmine and polyvalent antivenom in Indian common krait (Bungarus caeruleus) biteJournal of Infection and Public Health 3:83–87.https://doi.org/10.1016/j.jiph.2010.01.002
-
Distinctive epidemiologic and clinical features of common krait (Bungarus caeruleus) bites in Sri LankaThe American Journal of Tropical Medicine and Hygiene 79:458–462.https://doi.org/10.4269/ajtmh.2008.79.458
-
Retrospective study of case profile amongst snake bite patients managed at a tertiary care centre situated in a tribal district of Maharashtra during a span of three consecutive yearsInternational Journal of Scientific Research 7:10.
-
Acute renal failure in snake envenomation: a large prospective studySaudi Journal of Kidney Diseases and Transplantation : An Official Publication of the Saudi Center for Organ Transplantation, Saudi Arabia 19:404–410.
-
Snake bites in rural area of Maharashtra state, IndiaTropical Doctor 29:104–105.https://doi.org/10.1177/004947559902900215
-
Epidemiological aspects of acute childhood poisoning among patients attending a hospital at KolkataIndian Journal of Public Health 49:25–26.
-
Profile of snakebite envenoming in rural Maharashtra, IndiaThe Journal of the Association of Physicians of India 56:88–95.
-
Profile of snakebite envenoming in Western Maharashtra, IndiaTransactions of the Royal Society of Tropical Medicine and Hygiene 96:79–84.https://doi.org/10.1016/S0035-9203(02)90250-6
-
Snake bite: our experience at a tertiary care centre in North Western IndiaClinical Toxicology 47:727.
-
A study of clinical profile of snake bite at a tertiary care centreToxicology International 21:203–208.https://doi.org/10.4103/0971-6580.139811
-
Model-Based geostatistics the easy wayJournal of Statistical Software 63:1–24.https://doi.org/10.18637/jss.v063.i12
-
BookGeostatistics in small-area health applicationsIn: Lawson A. B, Banerjee S, Haining R. P, Ugarte M. D, editors. Handbook of Spatial Epidemiology. Chapman and Hall/CRC. pp. 211–224.
-
Epidemiological profile of snake-bite cases from Andhra Pradesh using immunoanalytical approachThe Indian Journal of Medical Research 125:661.
-
SoftwareGridded population of the world, version 4 (GPWv4)Population count adjusted to match 2015 revision of UN WPP country totals, beta release.
-
Surgical implications of snakebitesThe Indian Journal of Pediatrics 71:397–399.https://doi.org/10.1007/BF02725623
-
A profile of fatal snake bite cases in the Bankura district of West BengalJournal of Forensic and Legal Medicine 18:18–20.https://doi.org/10.1016/j.jflm.2010.11.007
-
Predictors of mortality in patients of poisonous snake bite: experience from a tertiary care hospital in Central IndiaInternational Journal of Critical Illness and Injury Science 4:101.https://doi.org/10.4103/2229-5151.134145
-
Pre-hospital treatment of snake envenomation in patients presented AT a tertiary care hospital in Northwestern IndiaJournal of Venomous Animals and Toxins Including Tropical Diseases 11:.https://doi.org/10.1590/S1678-91992005000300006
-
High or low- a trial of low dose anti snake venom in the treatment of poisonous snakebitesThe Journal of the Association of Physicians of India 61:387–396.
-
Profile of acute poisoning cases at Pravara rural hospital, LoniJournal of Indian Academy of Forensic Medicine 37:400.https://doi.org/10.5958/0974-0848.2015.00104.9
-
Adverse drug reaction profile of anti-snake venom in a rural tertiary care teaching hospitalJournal of Young Pharmacists 5:41–45.https://doi.org/10.1016/j.jyp.2013.02.003
-
Incidence of fatal snake bite in loni, Maharashtra: an autopsy based retrospective study (2004-2014)Journal of Forensic and Legal Medicine 39:61–64.https://doi.org/10.1016/j.jflm.2016.01.013
-
ConferenceUN population division’s methodology in preparing base population for projections: case study for IndiaAsian Population Studies. pp. 274–303.
-
Transporting snake bite victims to appropriate health facility within golden hour through toll free emergency ambulance service in India, save livesInternational Journal of TROPICAL DISEASE & Health 17:1–12.https://doi.org/10.9734/IJTDH/2016/25793
-
A study on the clinico-epidemiological profile and outcome of snake bite cases at a tertiary care hospitalIndian Journal of Pharmacology 45:S178–S179.
-
ReportHandbook on Treatment Guidelines for Snake Bite and Scorpion StingTamil Nadu Health Systems Project, Health and Family Welfare Department, Chennai, India.
-
Epidemiological profile and management of snakebite cases – A cross sectional study from Himachal Pradesh, IndiaClinical Epidemiology and Global Health 3:S96–S100.https://doi.org/10.1016/j.cegh.2015.11.007
-
Profile of poisoning cases at Belgaum, Karnataka: a cross sectional studyMedico-Legal Update 11:24–27.
-
A study on the clinico-epidemiological profile and the outcome of snake bite victims in a tertiary care centre in Southern IndiaJournal of Clinical and Diagnostic Research 7:e122.
-
Snake envenoming: a disease of povertyPLOS Neglected Tropical Diseases 3:e569.https://doi.org/10.1371/journal.pntd.0000569
-
A study on the acute kidney injury in snake bite victims in a tertiary care centreJournal of Clinical and Diagnostic Research : JCDR 7:853–856.https://doi.org/10.7860/JCDR/2013/5495.2957
-
Coagulation profile in snake bite victimsIndian Journal of Hematology and Blood Transfusion 31:S14.
-
WebsiteHole-filled SRTM for the globe version 4, available from the CGIAR-CSI SRTM 90m DatabaseAccessed February 12, 2020.
-
Snake bite mortality in children: beyond bite to needle timeArchives of Disease in Childhood 102:445–449.https://doi.org/10.1136/archdischild-2016-311142
-
Geostatistical survival models for environmental risk assessment with large retrospective cohortsJournal of the Royal Statistical Society: Series A 177:679–695.https://doi.org/10.1111/rssa.12041
-
Record based study of snake bite cases admitted at Shri Vasantrao Naik government medical college & hospital, Yavatmal (Maharashtra)Indian Journal of Public Health 50:35–37.
-
A multidisciplinary approach with hyperbaric oxygen therapy improve outcome in snake bite injuriesToxicology International 22:104–109.https://doi.org/10.4103/0971-6580.172287
-
Common krait (Bungarus caeruleus) bite in Anuradhapura, Sri Lanka: a prospective clinical study, 1996-98Postgraduate Medical Journal 78:276–280.https://doi.org/10.1136/pmj.78.919.276
-
An explicit link between gaussian fields and gaussian markov random fields: the stochastic partial differential equation approachJournal of the Royal Statistical Society: Series B 73:423–498.https://doi.org/10.1111/j.1467-9868.2011.00777.x
-
Human, snake, and environmental factors in human - snake conflict in North Bihar - A one-year descriptive studyChristian Journal for Global Health 3:36–45.https://doi.org/10.15566/cjgh.v3i1.77
-
Diet & nutrition profile of Chenchu population - A vulnerable tribe in Telangana & Andhra pradesh, IndiaIndian Journal of Medical Research, 141:688–696.
-
A comparative study between retrospective and prospective assessment of rationality and effectiveness of snake bite management in a tertiary care teaching hospitalInternational Journal of Basic & Clinical Pharmacology 8:1547.https://doi.org/10.18203/2319-2003.ijbcp20192571
-
Agriculture and snakebite in Bahia, Brazil - An ecological studyAnnals of Agricultural and Environmental Medicine 23:416–419.https://doi.org/10.5604/12321966.1219179
-
A study on snake bite poisoning in a tertiary care hospital in rural Odisha, IndiaJournal of Punjab Academy of Forensic Medicine & Toxicology 19:108.https://doi.org/10.5958/0974-083X.2019.00021.9
-
Management pattern of poisoning cases in a tertiary care center: what can be done at health care and national levels to improve outcomes?Indian Journal of Physiology and Pharmacology 56:210–211.
-
Snakebite mortality in India: a nationally representative mortality surveyPLOS Neglected Tropical Diseases 5:e1018.https://doi.org/10.1371/journal.pntd.0001018
-
Clinico-epidemiological features of viper bite envenomation: a study from Manipal, South IndiaSingapore Medical Journal 53:203–207.
-
Survey on pattern of snake bite cases admitted in South Indian tertiary care hospitalsInternational Journal of Pharmaceutical Sciences and Research 6:4362–4367.
-
Snake Bite Envenomation in a Tertiary Care CentreThe Journal of the Association of Physicians of India 66:55.
-
Pattern of use of antibiotics following snake bite in a tertiary care hospitalJournal of Clinical and Diagnostic Research 9:05–09.https://doi.org/10.7860/JCDR/2015/14753.6322
-
Study of poisoning cases in an Indian tertiary care teaching hospitalJournal of Indian Academy of Forensic Medicine 37:165.https://doi.org/10.5958/0974-0848.2015.00040.8
-
A Clinico-Epidemiological profile of neuroparalytic snake bite: using low dose ASV in a tertiary care centre from North IndiaThe Journal of the Association of Physicians of India 64:16–20.
-
Clinical profile and outcome of envenomous snake-bite at tertiary care centre in Western MaharashtraInternational Journal of Medicine and Public Health 1:28–38.https://doi.org/10.5530/ijmedph.4.2011.7
-
Profile of snake bite envenomation in rural Bangalore presenting to a tertiary care hospital, rural Bangalore, South IndiaIndian Journal of Critical Care Medicine 17:26.
-
Management of snake-bite in rural Maharashtra: a 10-year experienceThe National Medical Journal of India 18:71–75.
-
SoftwareR: a language and environment for statistical computingR Foundation for Statistical Computing, Vienna, Austria.
-
Snakebite profile from a medical college in rural setting in the hills of Himachal Pradesh, IndiaIndian Journal of Critical Care Medicine 18:134–138.https://doi.org/10.4103/0972-5229.128702
-
Study the assessment of poisoning cases in a rural tertiary care teaching hospital by a clinical pharmacistAsian Journal of Pharmaceutical and Clinical Research 5:138–141.
-
Pattern and outcome of acute poisoning cases in a tertiary care hospital in Karnataka, IndiaIndian Journal of Critical Care Medicine 13:152–155.https://doi.org/10.4103/0972-5229.58541
-
WebsiteDistrict census handbook (DCHB) - Part B - Primary census abstract 2011Accessed March 1, 2019.
-
SoftwareSample Registration System statistical report 2017New Delhi Office of the Registrar General of India.
-
Approximate bayesian inference for latent gaussian models by using integrated nested Laplace approximationsJournal of the Royal Statistical Society: Series B 71:319–392.https://doi.org/10.1111/j.1467-9868.2008.00700.x
-
Snake bite poisoning presenting as early morning neuroparalytic syndrome in Jhuggi dwellersThe Journal of the Association of Physicians of India 34:415–417.
-
Management of snake bite victims in a tertiary care intensive care unit in North IndiaIndian Journal of Critical Care Medicine 18:544–545.https://doi.org/10.4103/0972-5229.138164
-
Snakebite envenomation: a comprehensive evaluation of severity, treatment and outcome in a tertiary care South Indian hospitalThe Internet Journal of Emergency Medicine 5:.
-
Clinical profile, species-specific severity grading, and outcome determinants of snake envenomation: an Indian tertiary care hospital-based prospective studyIndian Journal of Critical Care Medicine 16:187–192.https://doi.org/10.4103/0972-5229.106499
-
Epidemiological profile of snake bite at a tertiary care hospital bengaluru, IndiaArchives of Pharmacy Practice 6:58–63.https://doi.org/10.4103/2045-080X.160991
-
Trends of poisoning in a tertiary care centre of North West Uttar PradeshJournal of Indian Academy of Forensic Medicine 37:396.https://doi.org/10.5958/0974-0848.2015.00103.7
-
Assessment of poisoning cases in tertiary care hospital - VIMS, BellaryIndian Journal of Pharmacology 45:S273.
-
Profile of fatal poisoning in and around Bidar, KarnatakaIndian Journal of Forensic Medicine and Toxicology 6:127–129.
-
Human fatalities caused by animal attacks: a six-year autopsy studyMedico-Legal Journal 85:194–199.https://doi.org/10.1177/0025817217707166
-
An epidemiological study of animal bites and envenomings in a rural district of Tamil Nadu, IndiaOnline Journal of Health and Allied Sciences.
-
Demography, clinical profile, morbidity and mortality pattern of snake bite cases in children: a study at tertiary teaching hospitalInternational Journal of Contemporary Pediatrics 6:1472.
-
Snake bite: a neglected problem in twenty-first century IndiaThe National Medical Journal of India 24:321–324.
-
New approaches & technologies of venomics to meet the challenge of human envenoming by snakebites in IndiaThe Indian Journal of Medical Research 138:38–59.
-
BookGopalakrishnakoneIn: Faiz P, Fernando A, Gnanathasan R, editors. Clinical Toxinology in Asia Pacific and Africa. Springer. pp. 115–136.
-
Venom, antivenom production and the medically important snakes of IndiaCurrent Science 103:635–643.
-
WebsiteGuidelines for the production control and regulation of snake antivenom immunoglobulinsWHO Press. Accessed August 1, 2019.
-
WebsiteVenomous snakes distribution and species risk categoriesWHO Press. Accessed August 1, 2019.
Article and author information
Author details
Funding
University of Toronto
- Prabhat Jha
Canadian Institutes of Health Research (FDN154277)
- Prabhat Jha
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank David Lightfoot for assistance with the literature search and Peter Rodriguez and Leslie Newcombe for data support.
Copyright
© 2020, Suraweera et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 48,114
- views
-
- 1,583
- downloads
-
- 168
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.