RESEARCH ARTICLE Estimating the Global Burden of Endemic Canine Rabies Katie Hampson1*, Laurent Coudeville2, Tiziana Lembo1, Maganga Sambo3, Alexia Kieffer2, Michaël Attlan2, Jacques Barrat4, Jesse D. Blanton5, Deborah J. Briggs6, Sarah Cleaveland1, Peter Costa6, Conrad M. Freuling7, Elly Hiby8, Lea Knopf6, Fernando Leanes9, François-Xavier Meslin10, Artem Metlin11, Mary Elizabeth Miranda6, Thomas Müller7, Louis H. Nel12, Sergio Recuenco5, Charles E. Rupprecht6,13, Carolin Schumacher14, Louise Taylor6, Marco Antonio Natal Vigilato9, Jakob Zinsstag15, Jonathan Dushoff16, on behalf of the Global Alliance for Rabies Control Partners for Rabies Prevention OPEN ACCESS Citation: Hampson K, Coudeville L, Lembo T, Sambo M, Kieffer A, Attlan M, et al. (2015) Estimating the Global Burden of Endemic Canine Rabies. PLoS Negl Trop Dis 9(4): e0003709. doi:10.1371/journal. pntd.0003709 Editor: Marilia Sá Carvalho, Oswaldo Cruz Foundation, BRAZIL Received: March 31, 2014 1 Boyd Orr Centre for Population and Ecosystem Health, Institute for Biodiversity, Animal Health & Comparative Medicine, University of Glasgow, Glasgow, United Kingdom, 2 Sanofi Pasteur, Lyon, France, 3 Ifakara Health Institute, Ifakara, Tanzania, 4 ANSES—French Agency for Food, Environmental and Occupational Health and Safety, Rabies and Wildlife laboratory of Nancy, Atton, France, 5 Centers for Disease Control and Prevention, Atlanta, Georgia, United States of America, 6 Global Alliance for Rabies Control, Manhattan, Kansas, United States of America, 7 Friedrich-Loeffler-Institute—Federal Research Institute for Animal Health, Greifswald—Insel Riems, Germany, 8 International Companion Animal Management Coalition, Cambridge, United Kingdom, 9 Pan-American Health Organization, Rio de Janeiro, Brazil, 10 Independent rabies scientist, Geneva, Switzerland, 11 Institut Pasteur in Cambodia, Phnom Penh, Cambodia, 12 University of Pretoria, Pretoria, South Africa, 13 Ross University School of Veterinary Medicine, Basseterre, St. Kitts, West Indes, 14 Merial, Lyon, France, 15 Swiss Tropical and Public Health Institute, Basel, Switzerland, 16 McMaster University, Hamilton, Ontario, Canada *katie.hampson@glasgow.ac.uk, katiehampson@gmail.com Abstract Accepted: March 18, 2015 Published: April 16, 2015 Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication. Funding: This study was funded by the UBS Optimus Foundation (http://www.ubs.com/optimusfoundation) and the Wellcome Trust (095787/Z/11/Z). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: LC, AK and MA are employed by Sanofi Pasteur and CS is employed by Merial. This does not alter our adherence to all PLOS NTDs policies on sharing data and materials. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the Centers for Disease Control and Prevention. Background Rabies is a notoriously underreported and neglected disease of low-income countries. This study aims to estimate the public health and economic burden of rabies circulating in domestic dog populations, globally and on a country-by-country basis, allowing an objective assessment of how much this preventable disease costs endemic countries. Methodology/Principal Findings We established relationships between rabies mortality and rabies prevention and control measures, which we incorporated into a model framework. We used data derived from extensive literature searches and questionnaires on disease incidence, control interventions and preventative measures within this framework to estimate the disease burden. The burden of rabies impacts on public health sector budgets, local communities and livestock economies, with the highest risk of rabies in the poorest regions of the world. This study estimates that globally canine rabies causes approximately 59,000 (95% Confidence Intervals: 25-159,000) human deaths, over 3.7 million (95% CIs: 1.6-10.4 million) disability-adjusted life years (DALYs) and 8.6 billion USD (95% CIs: 2.9-21.5 billion) economic losses annually. The largest component of the economic burden is due to premature death (55%), followed PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 1 / 20 The Global Burden of Endemic Canine Rabies by direct costs of post-exposure prophylaxis (PEP, 20%) and lost income whilst seeking PEP (15.5%), with only limited costs to the veterinary sector due to dog vaccination (1.5%), and additional costs to communities from livestock losses (6%). Conclusions/Significance This study demonstrates that investment in dog vaccination, the single most effective way of reducing the disease burden, has been inadequate and that the availability and affordability of PEP needs improving. Collaborative investments by medical and veterinary sectors could dramatically reduce the current large, and unnecessary, burden of rabies on affected communities. Improved surveillance is needed to reduce uncertainty in burden estimates and to monitor the impacts of control efforts. Author Summary Rabies is a fatal viral disease largely transmitted to humans from bites by infected animals —predominantly from domestic dogs. The disease is entirely preventable through prompt administration of post-exposure prophylaxis (PEP) to bite victims and can be controlled through mass vaccination of domestic dogs. Yet, rabies is still very prevalent in developing countries, affecting populations with limited access to health care. The disease is also grossly underreported in these areas because most victims die at home. This leads to insufficient prioritization of rabies prevention in public health agendas. To address this lack of information on the impacts of rabies, in this study, we compiled available data to provide a robust estimate of the health and economic implications of dog rabies globally. The most important impacts included: loss of human lives (approximately 59,000 annually) and productivity due to premature death from rabies, and costs of obtaining PEP once an exposure has occurred. The greatest risk of developing rabies fell upon the poorest regions of the world, where domestic dog vaccination is not widely implemented and access to PEP is most limited. A greater focus on mass dog vaccination could eliminate the disease at source, reducing the need for costly PEP and preventing the large and unnecessary burden of mortality on at-risk communities. Introduction Rabies is a fatal viral infection that can infect all mammals, but domestic dogs cause over 99% of all human deaths from rabies [1]. Human rabies can be prevented through prompt administration of post-exposure prophylaxis (PEP) to victims of bites by rabid animals [2], and infection can be eliminated at source through sustained mass vaccination of reservoir populations [3]. Most industrialized countries have eliminated rabies from domestic dog populations. However, in the majority of developing countries, rabies remains endemic in domestic dog populations and poorly controlled [4]. Our focus is on the impacts of canine-adapted variants of the rabies virus, sustained predominantly or entirely by transmission in domestic dogs (it is unclear whether independent transmission in wildlife might be sufficient for maintenance in some areas [5–7]). Our definition therefore includes rabies cases or exposures caused by canine variants of rabies virus also transmitted from wildlife, though these are negligible compared to those transmitted by domestic dogs. PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 2 / 20 The Global Burden of Endemic Canine Rabies The human and economic costs of canine rabies are poorly known [1]. A major challenge to estimating the burden of rabies is the absence of reliable surveillance data for countries where the disease is most prevalent. Basic information on how many lives are lost to rabies and the economic costs of preventing disease amongst those exposed are needed to advocate for sustainable control programmes. Official reporting of incidence data on rabies and rabies exposures remains desperately poor in most canine rabies-endemic countries, and is increasingly recognized to grossly underestimate the true number of cases [8, 9]. Active surveillance studies highlight the disparities between officially recorded and likely occurring rabies deaths. These include recent studies from both Asia and Africa based on probability decision tree modelling [10, 11]; extensive verbal autopsy surveys [12]; community surveys [13, 14] and contact tracing [15], which all show much higher mortality than officially reported. Specific features of rabies contribute to the problem of underreporting. Death is inevitable following clinical onset and therefore a large number of rabies victims never report to health facilities and are never diagnosed. Misdiagnosis to other neurological syndromes is frequent, especially in malaria endemic regions [16]. Shortages of life-saving PEP [15] and centres that provide PEP for bite victims [17] and poorly monitored sales of PEP to private suppliers all complicate counting the number of rabies diagnoses made and the number of treatments given. These problems of PEP provision particularly increase the risks of disease among the rural poor, an already marginalized sector of society. Moreover, poor infrastructure and a lack of personnel and facilities for rabies surveillance and diagnosis in most developing countries means that only very limited data of questionable reliability are available. In the absence of either reliable mortality reporting systems or more widespread active surveillance studies, extrapolations are required to estimate the global burden of rabies. Predictive methods have been developed to overcome the underreporting of disease, including a probability decision-tree method to determine the likelihood of the onset of clinical rabies in humans following the bite of a suspect rabid dog [8]. Using this method and drawing on data from a limited number of countries, Knobel et al. estimated that canine rabies caused approximately 55,000 human deaths annually across Africa and Asia [18]. Since this 2005 study, more data have become available and the disease situation has changed, with concerted control efforts in some parts of the world [19], increased incidence in others [20], as well as emergence in previously rabies-free areas [21, 22]. An updated assessment of the global rabies burden is therefore required. The rabies burden is made up of different components. Societal costs include mortality and lost productivity from premature death, and morbidity from adverse events (AE) of vaccination using nerve tissue vaccines (NTVs) and psychological effects of exposure to this fatal disease, expressed as disability-adjusted life years (DALYs). Direct costs of PEP (depending on the use of rabies immunoglobulin (RIG), and the type of vaccine and regimen, for example intramuscular (IM) versus intradermal (ID) administration) and indirect costs of seeking PEP (travel and accommodation for multiple clinic visits and lost income) fall upon the medical sector and affected communities, whilst the veterinary sector typically incurs costs related to dog vaccination. Veterinary and medical sectors both have responsibility for surveillance costs. Livestock losses depend on the size of at-risk livestock populations and preventative measures taken, and impact both national economies and households. The goal of this study was to make the best possible estimate of the burden of rabies, both globally and on a country-by-country basis, by combining all available data sources in a modelling framework that allows us to estimate missing components. We built upon earlier model frameworks [8, 18] to assess the current status of canine rabies globally and provide countryspecific estimates of disease burden and associated economic costs. Our model relied on data collected from many different sources including published studies, international databases, PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 3 / 20 The Global Burden of Endemic Canine Rabies market data for vaccine use and expert opinion surveyed for this study. We established relationships between rabies mortality and rabies control measures, which we incorporated into our estimation methods. These relationships indicated how interventions could affect the future burden of disease. Methods Model Framework We adapted the probability decision-tree framework developed by Cleaveland et al. [8] for Tanzania and used by Knobel et al. [18] to estimate the burden of rabies in Africa and Asia. The model uses the product of bite incidence, the probabilities of (i) a biting animal being rabid, RP, (ii) a bite victim receiving PEP, PP, and (iii) in the absence of PEP, developing rabies, DP, to extrapolate human rabies deaths and DALYs. An economic component is included to calculate the costs of rabies prevention and control, such as PEP administration, surveillance and livestock losses from rabies (Fig 1 and Table 1). We parameterized the model using country-specific data or aggregated cluster estimates as described below. In endemic areas we assumed that dog rabies incidence depends on vaccination coverage in the dog population, and that the probability that a bite is by a rabid animal depends on Fig 1. Model framework used in this study for estimating the global burden of canine rabies. Probability steps correspond to the probability that a bite is by a rabid animal (RP); that the victim received post-exposure prophylaxis, PEP, (PP); and, in the absence of PEP, that the bite victim developed rabies (DP). Data inputs (Table 1) are shown in red and model outputs in blue. Red arrows show estimated relationships (Fig 2). DALY = Disability-Adjusted Life Year; HDI = Human Development Index. Further details are given in Table 1. doi:10.1371/journal.pntd.0003709.g001 PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 4 / 20 The Global Burden of Endemic Canine Rabies Table 1. Data sources and inference methods within framework used for estimating the burden of canine rabies. Input Dependencies Inference method Data source Persons seeking PEP for dog bite exposures, Bites None Cluster values applied where missing data Literature, surveys with Delphi process Probability that bite is due to a rabid animal, RP RPmax (in absence of dog vaccination) and dog vaccination coverage, VC ML fit between dog vaccination coverage and dog rabies incidence time series (Fig 2). Extrapolation based on country and cluster estimates of VC RPmax from [23]. Coverage and incidence data from REDIPRA reports and [50]. Vaccination coverage in the dog population, VC None Cluster values applied where missing data Literature, surveys and market data, with Delphi process. Dog population, Dogs Human population data Reported dog population estimates or projected from dog: human ratios using human population data UN human population data (esa.un. org/wpp), literature and surveys Dog Population Control (culling and sterilization) None Cluster values applied where missing data Literature and surveys Laboratory Surveillance for human and animal rabies cases None Cluster values applied where missing data Databases (www.who-rabies-bulletin. org/, www.oie.int/wahis_2/public/ wahid.php,siepi.panaftosa.org.br), literature and surveys Livestock losses (deaths of cattle, sheep and goats) Livestock, RP and VC ML fit between cross-sectional vaccination coverage and livestock incidence data (Fig 2) Literature [14, 23–25] Populations of cattle, sheep and goats, Livestock None FAO database country values FAO (kids.fao.org/glipha/) Probability of receiving PEP following exposure by a rabid animal, PP HDI, officially reported deaths and bites MLE fitted relationship between HDI and probability of receiving PEP (Fig 2) UN (hdr.undp.org/en/statistics/hdi), literature and surveys Probability of developing rabies in the absence of PEP following a rabies exposure, DP None Literature [15] PEP regimen including whether RIG administered, type of vaccine, route of administration and clinic visits required None Literature, surveys and market data Literature, surveys and market data Life tables for DALY calculations None Global Burden of Disease Study, 2010 [28] Disability Weightings for DALY calculations None Literature [18] PEP = post-exposure prophylaxis, HDI = Human Development Index, FAO = Food and Agriculture Organization of the United Nations, REDIPRA = Directors of National Programs to Control Rabies in the Americas, ML = Maximum Likelihood, Bites = Bite Incidence (the same annotation is used in Fig 1). doi:10.1371/journal.pntd.0003709.t001 incidence. Using concurrent time series from the Americas we identified the best fitting relationship between incidence, I, and vaccination coverage, VC, estimating the asymptote, Imax and exponent, S, of this relationship by maximum likelihood (Fig 2A): I ¼ Imax à ð1ÀVCÞ s We assume that the probability that a bite is by a rabid animal is proportional to incidence: RP ¼ RPmax à ð1ÀVCÞ s where RPmax (0.74) is the estimated proportion of bites due to rabid animals in countries with negligible vaccination coverage [23]. Using this relationship we generated country and cluster estimates of RP. We similarly inferred the relationship between livestock rabies incidence and PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 5 / 20 The Global Burden of Endemic Canine Rabies Fig 2. Functional relationships estimated between A) rabies incidence in domestic dogs and average biannual dog vaccination coverage and B) probability of receiving post-exposure prophylaxis (PEP) and the Human Development Index. Upper inset in A shows the relationship on a log scale and the lower inset shows the relationship between rabies incidence in livestock and vaccination coverage in domestic dogs. Grey shading shows the 95% confidence intervals of the fitted relationships. doi:10.1371/journal.pntd.0003709.g002 vaccination coverage, based on cross-sectional data from published studies [14, 23–25] (Fig 2A inset). We assumed that in the absence of PEP, 19% of victims bitten by rabid dogs develop rabies and die (DP = 0.19, see [26]). The probability that a bite victim receives PEP (anti-rabies vaccine, sometimes supplemented by RIG for severe exposures) determines the likelihood of progression to rabies and death, however, we know of only one study quantifying this probability, PP [14]. We therefore used an alternative inference method. Reporting of PEP use and of rabies deaths varies according to health infrastructure, however both can be very poor in developing countries. If we assume equivalent reporting rates of rabies deaths, D, and PEP use, T, we can use the ratio D/T, as well as the probability that a reported bite was by a rabid dog, RP, and the probability of developing rabies in the absence of PEP, DP, to calculate the probability that a bite victim receives PEP: PP ¼ RP à DP = ððRP à DPÞ þ D=TÞ We examined whether a relationship exists between country economic/welfare measures (the Human Development Index, HDI) and these estimates of the probability of receiving PEP, PP (Fig 2B) using a generalized linear mixed-effects model [27] with country as a random effect. We generated country-specific parameter estimates of PP and used bootstrap resampling from the fitted mixed model for sensitivity analyses and to generate prediction intervals. We used these parameter estimates together with data on bite incidence (see below) to generate estimates of rabies exposures and deaths. Disease burden was expressed in terms of standard DALYs and calculated in accordance with methods developed by the World Health Organization (WHO). DALY calculations involve two components: Years of Life Lost (YLL), which, for rabies, captures deaths due to infection, and the Years of Life lived with Disability (YLD), which, for rabies, captures disability following AEs due to use of nerve tissue vaccines, which are still used in some parts of the world. Total YLL due to rabies were estimated using the reference-standard life table from the 2010 Global Burden of Disease Study [28] and the age distribution of rabies cases and PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 6 / 20 The Global Burden of Endemic Canine Rabies exposures from previous research [18](Table 1). We provide burden estimates based on ageand country-specific mortality rates for comparison (S1 Table). The disability weighting used to calculate AEs from nerve tissue vaccines was based on previous research [18]. We explored an additional component in YLD as the anxiety associated with dog bites that may develop into rabies, assuming that the disability level associated with anxiety was 0.108 and that anxiety lasted for 60 days. However, we did not include this in the total disease burden calculation, due to a lack of data to validate these assumptions. The data, parameter values and code to replicate all the analyses are provided in the Supporting Information (S1 compressed file archive). Clustering of Countries We aggregated countries into clusters on the basis of similar rabies epidemiological situations, socioeconomic conditions and geographical proximity. Countries were classified as canine rabies-free based on historical freedom or literature reporting canine rabies elimination. Oceania, Western Europe, the US and Canada and canine rabies-free countries in Asia (Japan, Malaysia, Singapore, and Republic of Korea) comprised four canine rabies-free clusters that were not included in further analysis. The disease burden was estimated in endemic countries within clusters. Where possible, country level data was used to parameterize the model, but, in the absence of country-specific information, we applied the average estimate from countries within the cluster. Data Sources Data were obtained from: 1. Surveys involving the medical, veterinary and laboratory sectors. The information gathered included reported clinical and laboratory-confirmed rabies cases in humans and animals, protocols for and expenditure related to PEP and costs of control efforts (vaccinating, sterilizing and killing dogs). Surveys were translated into French, Spanish, Portuguese and Russian and made available online. We solicited responses from country representatives, particularly those attending regional rabies meetings or identified through regional networks (such as the Southern and Eastern Africa Rabies Group, SEARG, and Directors of National Programs to Control Rabies in the Americas, REDIPRA [19]), responsible for reporting on rabies surveillance and diagnosis in their countries. Data was collected from 136 respondents (spanning all sectors and 45 countries). These surveys provided valuable data on rabies prevention and control practices, including PEP protocols and unit costs. However, most quantitative data on incidence were incomplete and we therefore used published data instead (see point 2). 2. An extensive literature search for estimates of human and animal rabies, dog bite incidence, control efforts and associated economic costs. We searched Web of Knowledge and PubMed for publications from 2000 to 2013 using ‘rabies’ AND ‘dogs’ as key words and resulting papers were reviewed to determine their relevance. In addition to scientific publications within the standard search, we collated technical reports and presentations from regional meetings, soliciting data through members of the Global Alliance for Rabies Control Partners for Rabies Prevention (PRP, rabiesalliance.org/about-us/partners). We identified 551 articles addressing canine rabies and its control, but useable quantitative information was only available from a smaller subset (113, see the supporting bibliography, S1 text), with considerable underreporting evident in official reports of bite incidence and rabies deaths from low and middle-income countries. We used estimates of bite incidence from empirical studies involving active surveillance and only PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 7 / 20 The Global Burden of Endemic Canine Rabies incorporated official data where no other sources were available and where these were deemed valid by the PRP group. We used the most recent data available since 2000. Searches were last updated on 1 June 2013. 3. International databases for country-specific estimates of human populations, demographic rates, economic indicators and livestock (detailed in Table 1). GDP and 2010 exchange rates were from the International Monetary Fund (IMF) databases; 2010 human population estimates and HDI estimates were from the United Nations; health costs were from the WHO (WHO-CHOICE database; CHOosing Interventions that are CostEffective), and livestock populations were from the Food and Agriculture Organization (FAO) Global Livestock Production and Health Atlas. The World Organization for Animal Health (OIE) World Animal Health Information System (WAHID) and surveillance databases for specific regions were used to obtain recent country-specific reports of laboratory confirmed rabies cases and diagnostic tests performed. 4. Estimates of regional markets for dog rabies vaccines from Merial and for human postexposure vaccines and RIG from Sanofi Pasteur. These were used to validate dog vaccination coverage estimates or to provide estimates of coverage for areas where no other source of information was available. If market estimates differed from other sources by >10%, experts were consulted (amongst the PRP group) and a 2-round delphi process [29] used to obtain consensus on values used in the analysis. Economic Costs The economic cost of deaths due to rabies was estimated using the human capital approach based on productivity losses. For each rabies death, the number of discounted life years lost was based on the age distribution of rabies deaths and remaining life expectancy using the reference-standard life table from the 2010 Global Burden of Disease Study [28]. Productivity losses were calculated by weighting the life years lost by the country-specific GDP per capita without discounting. Estimates were also calculated using age-weighting and time discounting for comparison with other studies. Time lost by victims and accompanying care-givers (assuming all minors were accompanied by one adult) whilst seeking PEP was incorporated into economic losses. Country estimates of unit costs for delivering PEP, dog vaccination and surveillance were largely obtained from surveys (detailed in 1). The total economic cost of rabies was obtained by combining data on unit cost per case, livestock losses, and costs of control and prevention. We updated reported costs to 2010 US dollars (USD) using IMF statistics. We corrected for international differences in medical costs using the WHO-CHOICE database. Indirect costs were corrected for differences in income using the ratio of income per capita (expressed in International dollars, I$) calculated using IMF statistics. Direct non-medical costs were corrected only for differences in purchasing power. Livestock losses due to dog rabies were extrapolated from the inferred relationship between coverage and incidence (Fig 2A inset) and using livestock population estimates from FAO. Incidence of rabies in cattle was multiplied by the cost per head of cattle. Using the fitted relationship between dog vaccination coverage and reported rabies incidence in livestock [14, 23, 25, 30] described above, and country and cluster values for dog vaccination coverage we estimated livestock losses. We converted the costs of different livestock into the costs of cattle using FAO livestock unit measures [31]. PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 8 / 20 The Global Burden of Endemic Canine Rabies Sensitivity Analyses Uncertainty was modelled by drawing from distributions for each parameter estimate. Bite incidence and dog vaccination coverage were modelled using triangular distributions, with maxima and minima set according to cluster ranges agreed from the 2-round Delphi process. Uncertainty in PEP probability was modelled using bootstrap resampling from the fitted mixed model described above. Uncertainty in the other probabilities (RP, and DP) was modelled using permutation-based resampling based on the original binomial sample. The range of variation of the results of our analysis was assessed using 1,000 Monte Carlo simulations. Parameters that varied across countries or clusters used the same quantile draw globally for each realization. Results Data Validation and Model Fitting Rabies incidence in dogs was best described by the fitted model (Fig 2A): I ¼ 0:002 à ð1ÀVCÞ 1:9 with dog vaccination coverage corresponding to the average from the previous two years. We used this relationship to infer the probability that a bite was due to a rabid animal, RP. We inferred a similar relationship between rabies incidence in livestock, IL and dog vaccination coverage (Fig 2A inset): IL ¼ 0:0017 à ð1ÀVCÞ 9 We found a strong significant relationship between inferred probability of receiving PEP and both country GDP and HDI (p<0.001). We used the relationship with HDI, which better captures inequalities within countries (Fig 2B) to generate estimates of the probability of receiving PEP, PP. The parameter estimates RP and PP are detailed by cluster in Table 2 and by country in S1 Table. Using these relationships and the data described above, we implemented the probability decision tree model to generate burden estimates. Human Rabies Deaths and DALY Estimates Results of predicted rabies mortality, morbidity and DALYs are provided by cluster in Table 2 and by country in S1 Table. We estimated that around 59,000 [95% CIs: 25,000–159,200] human rabies deaths occur annually globally, with the vast majority of these in Africa (36.4%) and Asia (59.6%). Less than 0.05% of estimated deaths occurred in the Americas [182, 95% CIs: 84–428], of which over 70% were from Haiti. India, with 35% of human rabies deaths, accounted for more deaths than any other country, but the estimated per-person death rate was highest in the poorest countries in sub-Saharan Africa. The global distribution of estimated human deaths and death rates due to rabies is illustrated in Fig 3. The parameters that had the greatest impact on variation in estimated human rabies mortality were bite incidence, followed by the probability of receiving PEP and the probability of developing rabies after a rabid animal bite in the absence of PEP (Fig 4). Globally, around 3.7 [1.6–10.4] million DALYs were estimated to be lost due to rabies, with over 95% lost in Africa (36.2%) and Asia (59.9%) and less than 0.5% (11,950 DALYs) in the Americas. The vast majority (>99%) of DALYs lost (3.68 million) were due to the premature death of rabies victims (YLL). A very small part (0.8%) of the DALY score (30,400) was due to AEs (in terms of YLD) of outdated, mostly locally produced, nerve tissue vaccines still in use in at least 10 countries in 2010. Anxiety due to suspect rabid dog bites potentially accounts for a PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 9 / 20 7,196 [4–16] Congo Basin PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 15,694,905 [4501– 28455] 29,175,105 [9907– 51132] 233,883 [116–359] 672,177 [353–1531] 308,522 [96–596] 37,854 [31–80] 128,538 [88–706] 73,557 [65–82] 427,604 [282–485] 168,104 [99–416] 517,409 [200–1227] 350,374 [317–566] 116,433 [86–189] 403,632 [120–699] 242,725 [271–3013] 8,209,470 [2832– 16149] 14,943,066 [3811– 19223] 862,641 [502–2123] 335,740 [220–1612] 1,143,377 [418– 2076] PEP 12,311 [12–198] 2,923,041 [846– 5892] 21,960 [10–41] 47,458 [23–127] 12,421 [4–29] 3,434 [1–8] 1,564 [1–9] 3,878 [2–6] 6,492 [3–10] 7,763 [2–20] 45,738 [19–147] 43,080 [30–82] 15,548 [10–29] 3,683,974 [1304– 10433] 14,310 [6–39] 117,116 [46–368] 1,948 [0–5] 189 [0–1] 472 [0–3] 8,581 [4–17] 1,023 [0–2] 1,108 [0–3] 395,297 [155– 1702] 375,023 [206–971] 449,382 [244– 1031] 123,074 [38–467] 20,384 [22–293] 35,124 [11–80] 1,301,865 [377– 3436] 374,851 [60–674] 16,521 [10–83] 152,263 [69–799] 338,639 [76–635] YLL rabies 849,658 [293– 1974] 1,574,518 [313– 2345] 49,675 [32–191] 44,068 [26–227] 140,277 [44–275] Prevented deaths 30,359 [12– 86] [0–0] [0–0] [0–0] 81 [0–0] 23 [0–0] [0–0] [0–0] 473 [0–1] 2,867 [2– 11] 3,714,333 [1316– 10519] 14,310 [6–39] 117,116 [46–368] 1,948 [0–5] 270 [0–1] 495 [0–3] 8,581 [4–17] 1,023 [0–2] 1,582 [0–4] 398,164 [157– 1713] 375,023 [206–971] 449,382 [244– 1031] [0–0] [0–0] 123,074 [38–467] 12,311 [12–198] 1,301,865 [377– 3436] 374,851 [60–674] 16,521 [10–83] 160,801 [75–853] 357,015 [80–655] DALYs total [0–0] [0–0] [0–0] [0–0] [0–0] 8,538 [6– 54] 18,376 [4– 19] DALY NTVs 32% 22% 62% 52% 58% 41% 69% 57% 23% 10% 9% 10% 24% 15% 14% 36% 5% 9% Dog Vaccination Coverage % 0.476 0.474 0.218 0.176 0.166 0.191 0.080 0.194 0.661 0.718 0.701 0.562 0.442 0.545 0.555 0.337 0.673 0.671 RP 0.991 0.980 0.999 0.998 0.994 0.992 0.997 0.997 0.913 0.837 0.829 0.892 0.972 0.976 0.996 0.995 0.952 0.973 PP doi:10.1371/journal.pntd.0003709.t002 Kyrgyzstan, Mongolia, the Russian Federation, Turkmenistan, Tajikistan, and Uzbekistan. and Democratic People’s Republic of Korea; SADC comprises countries in the Southern African Development Community, Eurasia comprises Afghanistan, Kazakhstan, PEP = post-exposure prophylaxis, YLL = Years of life lost, DALY = Disability Adjusted Life Year, NTVs = adverse events from Nerve Tissue Vaccines. Asia 4 comprises the Philippines, Sri Lanka, Thailand (High PEP use); Asia 3 comprises Bhutan, Nepal, Bangladesh, Pakistan (Himalayan region); Asia 2 comprises Cambodia, Myanmar, Laos, Vietnam *95% Confidence Intervals (Thousands) Estimates by country are provided in S1 Table including which cluster countries were assigned to. 58,991 [21– 167] TOTAL 116,785 [55–192] 259,650 [129–625] 1,875 [1–6] 229 [0–1] Eurasia Middle East 65,538 [19–138] 31 [0–0] 18,090 [3–40] 8,272 [4–44] Eastern Europe 8 [0–0] Central America 21,133 [16–28] 3 [0–0] 137 [0–0] Caribbean 34,255 [18–46] 40,950 [12–93] 274,041 [118–818] 258,341 [214–431] 119,707 [105–186] 195,237 [61–404] 108,322 [118–1428] 4,581,603 [1553– 9619] 8,318,530 [1514– 10993] 262,841 [179–917] 244,767 [144–1139] 766,842 [238–1313] Exposures Southern Cone 18 [0–0] 16 [0–0] Andean Brazil 6,005 [3–15] 1,971 [1–7] North Africa 6,330 [2–27] 197 [0–3] Indonesia West Africa 20,847 [7–55] India SADC 265 [0–1] 6,002 [1–11] 2,438 [1–13] Asia 3 Asia 4 5,423 [1–10] Asia 2 China Deaths [95% CI]* Cluster Table 2. Estimates of rabies deaths, exposures, PEP use, prevented deaths, DALYs (due to rabies and to NTVs), and average dog vaccination coverage, probability that a dog is rabid (RP), and probability of receiving PEP (PP) by cluster in canine rabies endemic regions. The Global Burden of Endemic Canine Rabies 10 / 20 The Global Burden of Endemic Canine Rabies Fig 3. The distribution of the global burden of rabies: A) human rabies deaths, B) per capita death rates (per 100,000 persons), and C) expenditure on dog vaccination (per 100,000 persons). Countries shaded in grey are free from canine rabies. doi:10.1371/journal.pntd.0003709.g003 substantial burden (518,000 DALYs in terms of YLD), although we do not include this component in our total estimate due to a lack of data to validate assumptions about disability weighting (0.108) and its application (60 days for true exposures). Fig 4. Model sensitivity to parameter uncertainty. PEP = post-exposure prophylaxis. doi:10.1371/journal.pntd.0003709.g004 PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 11 / 20 The Global Burden of Endemic Canine Rabies Economic Costs The overall economic costs of canine rabies were estimated as 8.6 billion USD (95% CIs: 2.9– 21.5 billion). These costs were mainly due to productivity losses from premature deaths (2.27 billion USD), direct expenditure on PEP (totalling 1.70 billion USD) and lost income whilst seeking PEP (1.31 billion USD). However, there was considerable variation in the breakdown of costs by region (Fig 5): the largest proportion of costs was due to premature death in Asia and Africa; much less was lost due to PEP (direct costs, and from travel and lost income whilst seeking PEP) in Africa compared to Asia and the Americas, and in the Americas a large proportion of costs were due to dog vaccination. Livestock deaths amounted to 512 million USD per year, with major losses in African countries with livestock-dependent economies (e.g. Ethiopia, Sudan, and Tanzania) and in more populous countries in Asia (China, India, Bangladesh and Pakistan). Table 3 provides a breakdown of estimated costs and a country breakdown is given in S1 Table. Globally, over 70% of the estimated economic burden was societal (from premature deaths and losses from seeking PEP); 20% fell to the medical sector or to bite victims (direct costs) and ~8% to the veterinary sector or directly to communities, from livestock losses and control interventions (dog vaccination and population management i.e. culling and/or sterilization/ birth control). Only around 0.01% of costs were from laboratory-based surveillance. The breakdown of costs by region varied dramatically (Figs 5 and S1 and S1 Table gives country breakdown). Dog vaccination accounted for less than 1.5% (~$130 million) of the economic burden. In the Americas, average per capita expenditure on dog vaccination was approximately $0.11, amounting to almost 20% of the economic burden. In most other endemic low-income countries, per capita expenditure on dog vaccination was negligible (less than $0.02 or <2% of the economic burden respectively, Figs 3C and 5 inset). Unit costs differed greatly between countries both for dog vaccination (for example costing $6-7/dose in some West African countries, $1/dose in Laos, $0.45/dose in the Philippines, $0.5/dose in Chad and $0.2–0.3 in Tanzania) and human PEP (range: $11–150 per dose), as well as the regimens and types of vaccines and RIG used (compressed file archive). Most Fig 5. Division of costs associated with rabies, prevention and control across sectors by region. Inset shows proportional expenditure in different regions. The breakdown of costs by cluster is given in S1 Fig and Table 3 and detailed by country in S1 Table. doi:10.1371/journal.pntd.0003709.g005 PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 12 / 20 The Global Burden of Endemic Canine Rabies Table 3. Breakdown of economic costs of rabies by cluster in thousands of USD. Cluster Direct Costs Travel costs Lost Income Productivity losses from premature death Livestock losses Dog vaccination Dog population management Surveillance Asia 2 75.15 15.96 20.90 252.298 2.073 0.074 24.163 0.000 Asia 3 16.67 4.42 7.08 104.774 0.564 0.214 65.706 0.042 Asia 4 41.81 2.32 33.26 45.658 11.248 0.123 1.728 0.021 China 648.27 49.25 807.32 1,642.646 4.235 0.195 40.777 0.046 India 491.23 42.60 138.03 1,646.650 9.050 0.417 62.348 0.002 Indonesia 21.18 0.95 10.03 37.123 6.384 0.811 1.717 0.000 North Africa 38.11 1.45 19.70 106.002 2.756 1.013 89.661 0.040 Congo Basin 14.47 1.28 2.33 154.424 0.481 0.003 19.670 0.001 West Africa 48.53 2.94 5.37 313.348 6.684 0.026 60.086 0.004 SADC 55.00 4.81 19.61 199.579 4.600 0.263 110.129 0.079 Andean 32.00 1.24 27.26 7.867 10.753 0.396 3.130 0.014 Brazil 45.37 1.82 63.36 11.070 16.620 0.342 0.007 0.288 Caribbean 11.66 0.24 2.43 7.702 2.575 0.113 0.296 0.006 Central America 34.13 0.49 11.74 1.873 31.308 0.809 0.001 0.020 Southern Cone 6.18 0.13 9.00 1.730 4.710 0.521 8.753 0.007 Eastern Europe 51.09 1.84 41.88 18.350 10.460 0.053 0.627 0.062 Eurasia 28.93 2.17 68.66 89.932 4.451 1.413 13.758 0.083 Middle East 42.58 0.35 26.06 35.907 0.592 0.128 9.543 0.025 TOTAL 1702.35 134.28 1314.01 4676.93 129.55 6.91 512.10 0.74 Estimates by country are in S1 Table including which cluster countries were assigned to. Asia 4 comprises the Philippines, Sri Lanka, Thailand (High PEP use); Asia 3 comprises Bhutan, Nepal, Bangladesh, Pakistan (Himalayan region); Asia 2 comprises Cambodia, Myanmar, Laos, Vietnam and Democratic People’s Republic of Korea; SADC comprises countries in the Southern African Development Community, Eurasia comprises Afghanistan, Kazakhstan, Kyrgyzstan, Mongolia, the Russian Federation, Turkmenistan, Tajikistan, and Uzbekistan. doi:10.1371/journal.pntd.0003709.t003 countries with a high disease burden reported negligible use of RIG. Countries with more substantial RIG use were in Eastern Europe, North Africa and a few Asian countries (Sri Lanka, Thailand and the Philippines). A few countries reported use of NTVs, notably Ethiopia, several countries in Latin America (Peru, Venezuela, Bolivia, Honduras and Argentina), Myanmar and Pakistan, with Bangladesh discontinuing use of NTVs in 2011 [1] (though were classified as using NTVs for the year 2010). Only a few countries reported widespread ID vaccination, including the Philippines, Sri Lanka, Thailand, and to a much lesser extent in India. Discussion This study highlights that the mortality risks and per capita rabies burden fall disproportionately upon the poorest regions of the world, with impacts on local communities, public health sector budgets and livestock economies. Mortality and loss of economic productivity due to premature death are the most serious effects of canine rabies. Highest mortality rates occur in areas with limited dog vaccination, where PEP is the only lifeline for at-risk populations, yet PEP supply and distribution systems are wholly inadequate in many of these places and often very costly. PEP costs, the second largest component of the economic burden, could be reduced in many areas through more judicious and cost-effective administration. The methods PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 13 / 20 The Global Burden of Endemic Canine Rabies developed here shed light on important gaps in knowledge, provide a preliminary picture of the distribution of the rabies burden by country and underline the lack of investment in rabies control and prevention measures. Improved surveillance and reporting of bites and rabies cases is needed, both for better burden estimates, and most importantly to monitor the impacts of control efforts. We estimate that annually canine rabies causes around 59,000 deaths and 3.7 million DALYs, which is considerably higher than previous estimates [18], however this is largely due to methods. In our study we did not apply time discounting or age-weighting to calculate DALYs for consistency with the Global Burden of Disease (GBD) study 2010 [32, 33]. Our DALY estimates are correspondingly higher than the GBD study as a result of our higher estimates of mortality. Applying 3% discounting and age-weighting we estimated slightly higher but more consistent DALY estimates to the Knobel et al. 2005 study [18]. S1 Table shows differences in DALY estimates due to the life table used and age-weighting and time discounting. We estimated a much higher economic burden of around 8.6 billion USD annually, compared to 583.5 million USD reported by Knobel et al. [18]. These differences are due to the inclusion of lost income from premature deaths and from all countries (not just those of Africa and Asia), in addition to a major increase in estimated direct PEP costs from 300 million USD [18] to 1.7 billion USD. Reasons contributing to this increase include higher prices and increased availability of cell-culture vaccines compared to NTVs, now discontinued in most endemic areas. Furthermore, we estimated higher bite incidence and PEP use (29.2 versus 4.3 million PEP delivered, respectively), based on published data, with over 10 million PEP delivered in China alone [34]. According to our estimates, India, the world’s second most populous country (with close to 18% of the global population) accounts for over 35% of the global rabies burden (approximately 20,800 deaths). This is broadly consistent with a recent verbal autopsy study, that estimated ~12,700 deaths from furious rabies alone in India [12]. Generally, our country estimates are in line with active surveillance studies [10, 13]. However, our estimates of deaths and DALYs are considerably higher than the GBD study, which attributed only 26,400 deaths and 1.46 million DALYs worldwide to rabies in 2010 [32, 33]. The GBD study drew upon officially recorded data for rabies, which grossly underestimate the disease burden due to extensive underreporting from countries where rabies is most prevalent. A large proportion of rabies deaths are misdiagnosed in areas with high general mortality [16]. Hospitalization provides little palliative care and death is inevitable, therefore many victims, particularly those from poor communities, either do not attend a facility or do not stay until death. Most victims (>75%) die at home [10,11,13], and these deaths are absent from official records. Furthermore, few clinical rabies diagnoses will be made from verbal autopsies unless interviewers probe for a history of a dog bite (as in [12]). Therefore even in countries where rabies is notifiable, many rabies deaths are not recorded, and burden estimates must therefore rely on predictive approaches. Our model provides a point of comparison for individual countries, but, due to the nature of extrapolating across large and heterogeneous populations, inaccuracies are likely and active surveillance studies are warranted. A major question is how to quantify anxiety associated with a life-threatening bite from a rabid animal. Previous burden estimates have ignored this component. We were also unable to find empirical evidence to validate a disability weighting. However, using assumptions agreed upon by the PRP group, we show that anxiety could be substantial (>10% of the total burden, ~518,000 DALYs), but research is needed to validate this weighting and its application to bite victims. In our study, AEs from NTV use account for a very small proportion of DALYs (0.8%, 30,400 DALYs), and this has declined (44,900 DALYs in 2005 [18]) due to the discontinuation of NTVs in most countries. PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 14 / 20 The Global Burden of Endemic Canine Rabies Lost productivity due to premature death ($4.7 billion 55.2%) was the largest component of the economic burden followed by direct costs of PEP ($1.7 billion, 20.1%). Investment in PEP has reduced rabies deaths in some countries, but at a high cost [35], whilst there has been very little investment in dog vaccination (Figs 3 and 5). More judicious administration of PEP could substantially reduce PEP costs (as indicated by the divergence of estimates of PEP use and prevented deaths detailed in Table 2). For example, >1 million PEP are delivered annually in the Americas [36], but most are precautionary for healthy animal bites. Most countries use IM delivery of PEP, but substantial savings (>60%) could be achieved globally by switching to the more cost-effective ID route as recommended by WHO [1]. Indirect costs of seeking PEP (1.84 billion, both travel and lost income) are a major cost to households of exposed individuals. This is a particular problem in rural areas, since PEP is typically only available in urban centres (sometimes only capital cities [10]). Recent improvements in PEP provision have reduced mortality in some countries [13, 37] and should be considered more widely. Vaccination of dogs, the proven way of preventing human exposures and eliminating the disease at source, comprised a very small proportion of the economic burden (<1.5%). Outside North America and Europe, a large investment in dog vaccination has only been sustained in one region (~$0.11/person/year in the Americas). The result is that the rabies burden in the Americas is very small (<200 deaths per year across the continent, mostly in Haiti), in contrast to other countries where rabies is endemic and expenditure on dog vaccination is negligible (<$0.02/person/year). Unlike the international government-coordinated control effort across Latin America, many developing countries have relegated rabies control and prevention to the private sector with no regulatory requirements or incentives (for example, as part of structural adjustment programmes in sub-Saharan Africa). As a result, rabies has been neglected in comparison to economically important livestock diseases. Generally, medical sector costs were much higher than veterinary costs (Fig 5), but investment in dog vaccination could bring down costs to the medical sector, demonstrating the need for intersectoral coordination [38, 39]. The World Bank has been supporting the strengthening of veterinary services (e.g. through the OIE Performance of Veterinary Services pathway) and considering zoonoses prevention and control as a ‘public good’, but resources are still lacking. Standardization of vaccine procurement is greatly needed to assist poorer countries to implement mass vaccinations, given the wide variation in vaccine prices shown by our study. Vaccine banks such as those administered by OIE could have a pivotal role to play. Livestock losses are a relatively small component of the global economic burden (6%), but represent an important cost to impoverished and livestock dependent communities, particularly in Africa [40]. Our estimates should be considered on the low side due to the limited data, drawing from only a few cross-sectional studies [14, 23, 25, 30], including laboratory confirmed cases, which underestimate the true burden [24]. Better reporting of livestock cases (suspect and confirmed) and further active surveillance studies are therefore necessary. There are number of limitations to our study. The most critical relate to uncertainty surrounding parameter estimates, particularly in relation to bite incidence (Fig 4). A Bayesian Hierarchical approach could better incorporate uncertainties, but estimates will be constrained by the data scarcity and quality. By fitting the relationship between vaccination coverage and dog rabies incidence phenomenologically, using longitudinal time series, the effect of differences in surveillance quality between locations was reduced. High turnover of dog populations means that single vaccination campaigns have short-lived impacts, whereas sustained campaigns progressively reduce disease incidence [3, 41]. Using lagged average coverage from consecutive campaigns improved the model fit and meant estimates were less subject to stochastic fluctuations (Fig 2). However, we had little power to estimate incidence at negligible coverage levels due to large multiannual epidemic fluctuations, compounded by limited surveillance. PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 15 / 20 The Global Burden of Endemic Canine Rabies Consistency in estimates of the basic reproductive ratio for canine rabies [3, 42] provides reassurance for extrapolations, but localized heterogeneities and landscape and demographic characteristics will influence vaccination impacts. Improved surveillance should enable future use of more mechanistic dynamical models [43]. Nonetheless, the direct relationship between dog vaccination and disease incidence provides a logical means of comparison that frames the problem in terms of investment in disease control and prevention. The proportion of reported bites due to suspect rabid animals contributes uncertainty to our results (Fig 4), and varies according to rabies incidence and treatment-seeking behaviour (Fig 2). But we were unable to find data to classify this systematically. By assuming that the maximum probability that a bite is by a rabid animal was 0.74 in countries with negligible vaccination coverage [23], we set an upper limit on the mortality burden, with reduced mortality rates in countries with higher coverage. More generally, major uncertainties in treatment-seeking behaviour, PEP availability and dog bite incidence limited the accuracy of our estimates. For example, little is known about the variation in PEP seeking in different socio-economic and cultural settings. However, cumulative evidence, confirmed by responses to our questionnaire, shows that PEP accessibility is very poor, often restricted to capital cities in the poorest countries [10], whilst in richer countries, or where PEP is provided free-of-charge, PEP is sought more readily [24, 44]. We ignored mortality and costs due to imported cases in rabies-free countries, which can be individually expensive [45], but are negligible compared to endemic rabies. Our use of average per capita GDP will mean that productive losses are overestimated, as rabies disproportionately affects impoverished communities. A further limitation of our study is that the burden is not broken down between urban and rural areas, due to a lack of data. However, dog vaccinations are implemented mostly in urban areas, which are easiest to access; dog:human ratios are typically higher in rural areas [46–48]; and PEP access is best in capital cities. Hence most rabies cases are expected to be from rural areas [18]. Finally, our estimates do not include the impacts of wildlife-transmitted rabies (from terrestrial wildlife and bats maintaining rabies virus transmission independently from domestic dogs). However, as canine rabies accounts for well over 95% of all human cases, our estimates are expected to be close to the overall rabies mortality burden globally. On the other hand, livestock losses due to wildlife rabies (for example, vampire bat rabies in the Americas [1]), will add substantially to the economic burden of rabies in certain parts of the world. This study demonstrates that the global burden of canine rabies is substantial, even though the disease is entirely preventable. Success in tackling the problem is contingent on investment in dog rabies control, which we show has been severely lacking. Long-term mass dog vaccination efforts could reduce medical sector and societal costs, and elimination is feasible with currently available methods [40, 49], however innovative financing models are required to overcome institutional barriers. Supporting Information S1 Text. Supporting bibliography. (DOCX) S1 Table. Estimates by country of rabies deaths, exposures, PEP use, prevented deaths, dog vaccination coverage, probability that a dog is rabid (RP), of bite victims receiving PEP (PP), DALYs, costs and 95% confidence intervals of estimates. Clusters to which countries are assigned are shown and inputs used for estimating parameters including the human development index and whetehr a country s rabies-free or endemic (RISK). Estimates of years of life lost (YLL) and DALYs (due to rabies and to adverse events from the use of nerve tissue PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 16 / 20 The Global Burden of Endemic Canine Rabies vaccines) are shown under different assumptions (estimates under the assumption of no time discounting or age-weighting should be directly comparable to the 2010 Global Burden of Disease study). (XLSX) S1 Fig. Division of costs associated with rabies, prevention and control across sectors by cluster. Inset shows proportional expenditure in different clusters. Full details of countries by cluster are given in S1 Table. Asia 4 comprises: Philippines, Sri Lanka, Thailand (High PEP use); Asia 3 comprises Bhutan, Nepal, Bangladesh, Pakistan (Himalayan region); Asia 2 comprises Cambodia, Myanmar, Laos, Vietnam and Democratic People’s Republic of Korea; SADC comprises countries in the Southern African Development Community, Eurasia comprises Afghanistan, Kazakhstan, Kyrgyzstan, Mongolia, the Russian Federation, Turkmenistan, Tajikistan, and Uzbekistan. (EPS) S1 Dataset. Model code and input data files including references, rationale and detail of Delphi process. The code folder contains seven R scripts: burden_model.R runs the model using data compiled in burden_1.R, after estimating parameters using: FitCovInc.R, FitPP.R, and creating Fig 2 (RabiesBurdenFig2.R). The script burden_results.R summarizes findings using the output of burden_model.R and burden_sensitivity.R runs the sensitivity analyses. The data folder contains 12 csv files called by the R code for the analyses, and one excel file (Vet. xlsx) with additional details about the data sources in vcountry2.csv and vcluster2.csv and with Delphi process estimates for dog vaccination coverage. Data sources are detailed in the relevant data sources and the details of the sources of data used in the analysis are in the supporting bibliography, S1 text. (ZIP) Acknowledgments We are very grateful to the extended networks of contributors and survey respondents vital to this study. These include the Asian Rabies Expert Bureau (AREB), the African Rabies Expert Bureau (AfroREB), the Southern and Eastern African Rabies Group (SEARG), the Middle East & Eastern Europe Rabies Expert Group (MEEREB), the Pan-American Health Organization and REDIPRA. We are grateful to the many contributors to this work including, B. Bolker, H. Bourhy, B. Dodet, J. Maki, N. Rotsztajn, J. Yoder, S. Acha, R. Alarcoo, M.M. Alqadasi, M. Barrios, O. Bharti, L.S. Cambara, M. Colman, A. Deressa, B.N. Diot, B.F. Diuedonne, M. Domitien, C Domingues, A.S. Domingues, B. Douangngeun, P.M. Dzimbiri, L.P. Escobar, P.L. Ferrer, I.A. C. Flores, J.W.K. Githinji, L.R. Guzman, L.F. Gwehure, O. Jaimes, E.A. Jimenez Suazo, P. Kabambingabu, E. Kayesa, S. Khaiseb, S. Khatiwada, B.M. Leite, L. Lind, J.R.P. Marins, W.I.N. Majeste, K. Monyame, A. Ngangnou, C. Ngeleja, S. Ntegeyibizaza, A. Nunez, A.A. Okurut, G. R. Pant, R.J.P. Parrales, E.S. Patrick, F.E. Peralta, D. Perez, G.M. Peters, L.H. Queiroz, M.Rahman, P. Romijn, V.H. de Paz Rozalez, C. Sabeta, N. Samson, A.R. Sanchez, H.F. Sanchez, V.H. Sancho, L.Q. Santi, T.O al Shekaili, W. Shumba, I.S. Soto, Q. Tang, E.P. Torres, I. Tiembre, A.E. Valery, F. Vargos-Pino, R.T. Velíçsquez, J.Villatoro, S. Vong, A.A. Wadair, Y.Wenwu, J. Weyer, G. Willat, and Y.Z. Zhang. Author Contributions Conceived and designed the experiments: KH LC JD LHN CER DJB. Performed the experiments: KH MS AK LC. Analyzed the data: KH JD. Contributed reagents/materials/analysis PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 17 / 20 The Global Burden of Endemic Canine Rabies tools: KH LC MA JB JDB SC PC CMF EH LK FL FXM AM MEM TM LHN SR CER CS LT MANV JZ JD. Wrote the paper: KH TL. References 1. WHO. Expert Consultation on Rabies. Second Report. Geneva: WHO, 2013. 2. Hemachudha T, Laothamatas J, Rupprecht CE. Human rabies: a disease of complex neuropathogenetic mechanisms and diagnostic challenges. Lancet Neurology. 2002; 1(2):101–9. 3. Hampson K, Dushoff J, Cleaveland S, Haydon DT, Kaare M, Packer C, et al. Transmission Dynamics and Prospects for the Elimination of Canine Rabies. Plos Biology. 2009; 7(3):e1000053. doi: 10.1371/ journal.pbio.1000053 PMID: 19278295 4. Coleman PG, Fevre EM, Cleaveland S. Estimating the public health impact of rabies. Emerging infectious diseases. 2004; 10(1):140–2. PMID: 15078611 5. Zulu GC, Sabeta CT, Nel LH. Molecular epidemiology of rabies: Focus on domestic dogs (Canis familiaris) and black-backed jackals (Canis mesomelas) from northern South Africa. Virus Research. 2009; 140(1–2):71–8. doi: 10.1016/j.virusres.2008.12.003 PMID: 19118588 6. Sabeta CT, Mansfield KL, McElhinney LM, Fooks AR, Nel LH. Molecular epidemiology of rabies in bateared foxes (Otocyon megalotis) in South Africa. Virus Research. 2007; 129(1):1–10. PMID: 17537536 7. Bingham J, Foggin CM, Wandeler AI, Hill FWG. The epidemiology of rabies in Zimbabwe. 2. Rabies in jackals (Canis adustus and Canis mesomelas). Onderstepoort Journal of Veterinary Research. 1999a; 66:11–23. PMID: 10396757 8. Cleaveland S, Fevre EM, Kaare M, Coleman PG. Estimating human rabies mortality in the United Republic of Tanzania from dog bite injuries. Bulletin of the World Health Organization. 2002; 80(4):304– 10. PMID: 12075367 9. Deressa A, Ali A, Beyene M, Selassie BN, Yimer E, Hussen K. The status of rabies in Ethiopia: A retrospective record review. Ethiopian Journal of Health Development. 2010; 24(2):127–32. 10. Ly S, Buchy P, Heng NY, Ong S, Chhor N, Bourhy H, et al. Rabies Situation in Cambodia. Plos Neglected Tropical Diseases. 2009; 3(9):e511. doi: 10.1371/journal.pntd.0000511 PMID: 19907631 11. Tenzin, Dhand NK, Gyeltshen T, Firestone S, Zangmo C, Dema C, et al. Dog Bites in Humans and Estimating Human Rabies Mortality in Rabies Endemic Areas of Bhutan. PLoS Neglected Tropical Diseases. 2011; 5(11):e1391. doi: 10.1371/journal.pntd.0001391 PMID: 22132247 12. Suraweera W, Morris SK, Kumar R, Warrell DA, Warrell MJ, Jha P, et al. Deaths from Symptomatically Identifiable Furious Rabies in India: A Nationally Representative Mortality Survey. PLoS Negl Trop Dis. 2012; 6(10):e1847. doi: 10.1371/journal.pntd.0001847 PMID: 23056661 13. Hossain M, Ahmed K, Bulbul T, Hossain S, Rahman A, Biswas MNU, et al. Human rabies in rural Bangladesh. Epidemiology and Infection. 2012; 140(11):1964–71. PMID: 22185694 14. Jemberu WT, Molla W, Almaw G, Alemu S. Incidence of Rabies in Humans and Domestic Animals and People's Awareness in North Gondar Zone, Ethiopia. PLoS Neglected Tropical Diseases. 2013; 7(5): E2216–E. doi: 10.1371/journal.pntd.0002216 PMID: 23675547 15. Hampson K, Dobson A, Kaare M, Dushoff J, Magoto M, Sindoya E, et al. Rabies Exposures, Post-Exposure Prophylaxis and Deaths in a Region of Endemic Canine Rabies. PLoS Neglected Tropical Diseases. 2008; 2(11):e339. doi: 10.1371/journal.pntd.0000339 PMID: 19030223 16. Mallewa M, Fooks AR, Banda D, Chikungwa P, Mankhambo L, Molyneux E, et al. Rabies Encephalitis in Malaria-Endemic Area, Malawi, Africa. Emerging Infectious Diseases. 2007; 13(1):136–9. PMID: 17370529 17. Maganga S, Cleaveland S, Ferguson H, Lembo T, Simon C, Urassa H, et al. The Burden of Rabies in Tanzania and its Impacts on Local Communities. PLoS Negl Trop Dis. 2013; 7(11):e2510. doi: 10. 1371/journal.pntd.0002510 PMID: 24244767 18. Knobel DL, Cleaveland S, Coleman PG, Fevre EM, Meltzer MI, Miranda MEG, et al. Re-evaluating the burden of rabies in Africa and Asia. Bulletin of the World Health Organization. 2005; 83(5):360–8. PMID: 15976877 19. Vigilato MAN, Cosivi O, Knobl T, Clavijo A, Silva HMT. Rabies update for Latin America and the Caribbean. Emerging infectious diseases. 2013; 19(4):678–9. doi: 10.3201/eid1904.121482 PMID: 23750499 20. Wu X, Hu R, Zhang Y, Dong G, Rupprecht CE. Reemerging rabies and lack of systemic surveillance in People’s Republic of China. Emerging Infectious Diseases. 2009; 15(8):1159–64. doi: 10.3201/ eid1508.081426 PMID: 19751575 PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 18 / 20 The Global Burden of Endemic Canine Rabies 21. Putra AAG, Hampson K, Girardi J, Hiby E, Knobel D, Mardiana IW, et al. Response to a rabies epidemic, Bali, Indonesia, 2008–2011. Emerging infectious diseases. 2013; 19(4):648–51. doi: 10.3201/ eid1904.120380 PMID: 23632033 22. Lembo T, Craig PS, Miles MA, Hampson KR, Meslin FX. Zoonoses prevention, control, and elimination in dogs. Macpherson CNL, Meslin FX, Wandeler AI, editors 2013. 205–58 p. 23. Lembo T, Hampson K, Haydon DT, Craft M, Dobson A, Dushoff J, et al. Exploring reservoir dynamics: a case study of rabies in the Serengeti ecosystem. Journal of Applied Ecology. 2008; 45(4):1246–57. PMID: 22427710 24. Thiptara A, Atwill ER, Kongkaew W, Chomel BB. Epidemiologic Trends of Rabies in Domestic Animals in Southern Thailand, 1994–2008. American Journal of Tropical Medicine and Hygiene. 2011; 85 (1):138–45. doi: 10.4269/ajtmh.2011.10-0535 PMID: 21734139 25. Tenzin, Dhand NK, Dorjee J, Ward MP. Re-emergence of rabies in dogs and other domestic animals in eastern Bhutan, 2005–2007. Epidemiology and Infection. 2011; 139(2):220–5. doi: 10.1017/ S0950268810001135 PMID: 20492745 26. Shim E, Hampson K, Cleaveland S, Galvani AP. Evaluating the cost-effectiveness of rabies post-exposure prophylaxis: a case study in Tanzania. Vaccine. 2009; 27:7167–72. doi: 10.1016/j.vaccine.2009. 09.027 PMID: 19925948 27. Bates D, Maechler M, Bolker B, Walker S. lme4: Linear mixed-effects models using Eigen and S4. R package version 1.1–7. http://CRAN.R-project.org/package = lme4. 2014. 28. Murray CJL, Ezzati M, Flaxman AD, Lim S, Lozano R, Michaud C, et al. GBD 2010: design, definitions, and metrics. The Lancet. 2012; 380(9859):2063–6. doi: 10.1016/S0140-6736(12)61899-6 PMID: 23245602 29. Undurraga EA, Halasa YA, Shepard DS. Use of Expansion Factors to Estimate the Burden of Dengue in Southeast Asia: A Systematic Analysis. PLoS Negl Trop Dis. 2013; 7(2):e2056. doi: 10.1371/journal. pntd.0002056 PMID: 23437407 30. Tenzin, Sharma B, Dhand NK, Timsina N, Ward MP. Reemergence of Rabies in Chhukha District, Bhutan, 2008. Emerging infectious diseases. 2010; 16(12):1925–30. doi: 10.3201/eid1612.100958 PMID: 21122223 31. Chilonda P, Otte J. Indicators to monitor trends in livestock production at national, regional and international levels. Livestock Research for Rural Development. 2006; 18(8):117. 32. Lozano R, Naghavi M, Foreman K, AlMazroa MA, Memish ZA. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012; 381: 2095–2128. 33. Murray CJL, Vos T, Lozano R, AlMazroa MA, Memish ZA. Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012; 380: 2197–2223 doi: 10.1016/S0140-6736(12)61689-4 PMID: 23245608 34. Yin C-p, Zhou H, Wu H, Tao X-y, Rayner S, Wang S-m, et al. Analysis on factors related to rabies epidemic in China from 2007–2011. Virologica Sinica. 2012; 27(2):132–43. doi: 10.1007/s12250-0123244-y PMID: 22492004 35. Denduangboripant J, Wacharapluesadee S, Lumlertdacha B, Ruankaew N, Hoonsuwan W, Puanghat A, et al. Transmission dynamics of rabies virus in Thailand: Implications for disease control. BMC Infectious Diseases. 2005; 5. 36. Schneider MC, Belotto A, Ade MP, Hendrickx S, Leanes LF, de Freitas Rodrigues MJ, et al. Current status of human rabies transmitted by dogs in Latin America. Cadernos De Saude Publica. 2007; 23 (9):2049–63. PMID: 17700940 37. Lapiz SMD, Miranda MEG, Garcia RG, Daguro LI, Paman MD, Madrinan FP, et al. Implementation of an Intersectoral Program to Eliminate Human and Canine Rabies: The Bohol Rabies Prevention and Elimination Project. PMID: PLoS Neglected Tropical Diseases. 2012; 6(12) e1891 doi: 10.1371/journal. pntd.0001891 38. Zinsstag J, Durr S, Penny MA, Mindekem R, Roth F, Gonzalez SM, et al. Transmission dynamics and economics of rabies control in dogs and humans in an African city. Proceedings of the National Academy of Sciences of the United States of America. 2009; 106(35):14996–5001. doi: 10.1073/pnas. 0904740106 PMID: 19706492 39. Fitzpatrick MC, Hampson K, Cleaveland S, Mzimbiri I, Lankester F, Lembo T, et al. Cost-effectiveness of Canine Vaccination to Prevent Human Rabies in Rural Tanzania. Annals of Internal Medicine. 2014; 160(2):91–100. PMID: 24592494 40. Lembo T, Hampson K, Kaare MT, Ernest E, Knobel D, Kazwala RR, et al. The Feasibility of Canine Rabies Elimination in Africa: Dispelling Doubts with Data. PLoS Neglected Tropical Diseases. 2010; 4(2). PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 19 / 20 The Global Burden of Endemic Canine Rabies 41. Hampson K, Dushoff J, Bingham J, Bruckner G, Ali YH, Dobson A. Synchronous cycles of domestic dog rabies in sub-Saharan Africa and the impact of control efforts. Proceedings of the National Academy of Sciences of the United States of America. 2007; 104(18):7717–22. PMID: 17452645 42. Townsend S, Sumantra IP, Pudjiatmoko, Bagus GN, E. B, Cleaveland S, et al. Designing Programs for Eliminating Canine Rabies from Islands: Bali, Indonesia as a Case Study. PLoS Negl Trop Dis. 2013; 7 (8):e2372. doi: 10.1371/journal.pntd.0002372 PMID: 23991233 43. Simons E, Ferrari M, Fricks J, Wannemuehler K, Anand A, Burton A, et al. Assessment of the 2010 global measles mortality reduction goal: results from a model of surveillance data. Lancet. 2012; 379 (9832): 2173–8. doi: 10.1016/S0140-6736(12)60522-4 PMID: 22534001 44. Kumarapeli V, Awerbuch-Friedlander T. Human rabies focusing on dog ecology-A challenge to public health in Sri Lanka. Acta Tropica. 2009; 112(1):33–7. doi: 10.1016/j.actatropica.2009.06.009 PMID: 19540826 45. Lardon Z, Watier L, Brunet A, Bernède C, Goudal M, Dacheux L, et al. Imported Episodic Rabies Increases Patient Demand for and Physician Delivery of Antirabies Prophylaxis. PLoS Negl Trop Dis. 2010; 4(6): e723. doi: 10.1371/journal.pntd.0000723 PMID: 20582307 46. Davlin SL, VonVille HM. Canine rabies vaccination and domestic dog population characteristics in the developing world: A systematic review. Vaccine. 2012; 30(24):3492–502. doi: 10.1016/j.vaccine.2012. 03.069 PMID: 22480924 47. Kayali U, Mindekem R, Hutton G, Ndoutamia AG, Zinsstag J. Cost-description of a pilot parenteral vaccination campaign against rabies in dogs in N'Djamena, Chad. Tropical Medicine & International Health. 2006; 11(9):1482–. 48. Kaare M, Lembo T, Hampson K, Ernest E, Estes A, Mentzel C, et al. Rabies control in rural Africa: Evaluating strategies for effective domestic dog vaccination. Vaccine. 2009; 27(1):152–60. doi: 10.1016/j. vaccine.2008.09.054 PMID: 18848595 49. Rupprecht CE, Barrett J, Briggs D, Cliquet F, Fooks AR, Lumlertdacha B, et al. Can rabies be eradicated? Developments in Biological Standardization (Basel). 2008; 131:95–121. PMID: 18634470 50. Belotto AJ. Organization Of Mass Vaccination For Dog Rabies In Brazil. Reviews of Infectious Diseases. 1988; 10(SUPPL. 4):S693–S6. PMID: 3206082 PLOS Neglected Tropical Diseases DOI:10.1371/journal.pntd.0003709 April 16, 2015 20 / 20