Improving dog bite victim survey and estimation of annual human deaths due to suspected rabies cases in three selected Liberian cities and environs, 2008–2017

Rabies remains a public health challenge of unknown magnitude in Liberia in spite of the goal of ensuring that no human in the country dies of rabies by 2030. The annual prevalence of Dog Bite Victims (DBVs) and true load of Annual Human Deaths (AHDs) due to rabies were not known. We investigated three selected cities of Liberia for annual prevalence of DBVs and true load of AHD due to suspected rabies, using 10–year retrospective record, 2008–2017 obtained from Buchanan, Gbarnga, and Voinjama, three socio-economically important cities in post-conflict Liberia. Data were sourced at County Reference Hospitals and at the Liberia National Institute of Health for these cities and their local environs. In addition, household questionnaire survey was used to identify and audit data quality for unreported DBVs, and treatment received from traditional caregivers. The proportion was used to audit the 10-year data on unreported DBVs in the cities. Descriptive statistics was used to summarize annual DBVs over the 10-year period in the three cities, respectively. A standardized clinical decision tree model was used to estimate AHDs due to suspected rabies. Based on questionnaire survey, 140/365, 148/375 and 146/350 DBVs did not visit any orthodox health facility in Buchanan, Gbarnga and Voinjama cities, respectively in 2014. An estimated total of 559 DBVs died of suspected rabies in the three cities and their environs during the 10-year period. Mean yearly prevalence of DBVs was 179±106.82, 393±257.85 and 76.9±38.11 per 100,000 population, while mean AHDs due to suspected rabies was 14.3±8.47, 35.5±23.25, and 6.1±3.21 per 100,000 population in Buchanan, Gbarnga, and Voinjama cities, respectively. The present findings provide annual prevalence of suspected rabies cases, corrected for under-reporting in three selected cities of Liberia. The findings would be useful in planning for stepwise actions towards rabies elimination, ensuring that no human dies of rabies in Liberia by 2030.

Introduction An in-country assessment of Stepwise Action for Rabies Elimination (SARE) protocol was conducted in Liberia in 2018 with a score of 1.5 out of 5 ascribed to the country, indicating that the foundational requirements for developing a national rabies control program has been put in place [1]. Despite this development, the true load of dog bite victims (DBVs) and annual human deaths due to rabies were not known for all, but Monrovia, the capital city of Liberia [2]. The lack of cold facts about the true load of rabies deaths in the country makes the goal of its elimination by 2030 less certain. Efforts to restore public health surveillance in Liberia have been intensified for over a decade commencing five years after the Second Liberian Civil War, which lasted from 1999 to 2003. One Health approach to improving rabies surveillance and control in canine and human populations [3,4] has gradually improved. Nonetheless, the impact of these efforts might have been small in view of the fact that in Liberia, rabies remains a neglected public health challenge decades after the first report of its clinical diagnosis in DBVs in the country [5].
Based on a recent report of co-circulation of African and Asian lineages of rabies virus in Liberia [6], the epidemiology of rabies in the country might also be undergoing significant changes. The current situation of rabies in Liberia therefore demands more effective One Health action geared towards improving surveillance of DBVs and control of the disease at the human-animal interface, with smart allocation of scarce resources [2,6,7]. In order to address the critical aspect of under-reporting of DBVs associated with care-seeking preferences [8,9], this study was designed to determine the true prevalence of dog bite victims and annual human deaths due to suspected rabies in three selected cities in Liberia, namely Buchanan, Gbarnga and Voinjama.

Study areas
This study was conducted in three cities selected from three of fifteen counties of Liberia, being three of four counties that submitted DBV reports to the Ministry of Health and Social Welfare during 2008-2013. Findings on the fourth county, Montserrado county, are not included in this study because they have been reported elsewhere [2]. The three cities included in the present study are Buchanan, the capital city of Grand Bassa County; Gbarnga, the capital city of Bong County; and Voinjama, the capital city of Lofa County, Liberia (Fig 1). The first study area, Buchanan and environs (latitude 5.8872˚north, longitude 10.0304˚west) is a coastal city that shares boundary with the Atlantic Ocean on the southwest, with a population of 34,270 persons in 2008 [12,13]. Liberia has conducted only four censuses since 1960 till date. They were conducted in 1962, 1974, 1984 and 2008, respectively. Although no other more recent census was conducted in Liberia since 2008, the annual population growth rate in the city from 2009 to 2017 was estimated at 1.4%, 1.5%, and 1.2% in 2008-2011, 2012-2014, and 2015-2017, respectively, in Buchanan and environs. In 2008-2013, reports of cases of DBVs from Buchanan city to the Ministry of Health and Social Welfare (MoHSW) came only from one human hospital, the Liberian Government Hospital (LGH). By 2017, there were nine other clinics in the city that reported DBVs to the Liberian National Institute of Health (LNIH) because of integrated and enhanced disease surveillance and response deployed after Ebola virus disease outbreak in the country [14].
The second city, Gbarnga (latitude 7.0024˚north, longitude 9.4728˚west) in Bong County, is situated in the central region of Liberia. Gbarnga had a population of 34,046 persons in 2008 [12,13]. There were no other census figures produced since 2008, but the annual population growth, 2009-2017 for the city was estimated at 1.0% in Gbarnga and environs. In 2008 to 2013, Phebe Hospital and Nursing School (PHNS) alone reported DBV cases to the MoHSW. By 2017, DBVs were reported by 22 other clinics/hospitals, in addition to PHNS [14]. The third study area, Voinjama city and environs (latitude 8.4202 o N, longitude 9.7539 o W), had

Study design
The study was based on an explanatory sequential mixed-method design, which included retrospective clinical case review of dog bite victims, household questionnaire survey of DBVs, their care-seeking preferences, and multi-criteria analysis of PEP compliance for estimation of annual human deaths due to suspected rabies cases in three selected cities of Liberia. The three cities were purposively selected because they were among the only four cities that submitted monthly reports on DBVs to the Ministry of Health and Social Welfare (MOHSW) in Liberia from 2008 to 2013. The submission of reports was also partly because they were designated reference hospitals in their respective county. All data were anonymized to protect the privacy of patients. In most cases, residential addresses of DBVs were available but anonymised. The addresses indicated site names that were geospatially referenced to produce map points for a Geospatial Observation-Linked Data on Suspected Rabies Exposure at the Human-Dog interface (GOLD-SREHD) profile.

Case definition and classification of rabies exposure
During this study, hospital records of DBVs were classified, based on WHO guidelines for case classification of human exposure to rabies, as no exposure, suspected (or possible) exposure, probable exposure, and confirmed exposure cases [15,16]. Cases of no exposure included records of humans touching or feeding animals or having animal licks on intact skin. A case of suspected exposure was defined as a person who has had close contact (usually a bite or scratch) with a rabies-susceptible animal in (or originating from) a rabies-infected area. A case of probable exposure was defined as a person who has had close contact (usually a bite or scratch) with an animal displaying clinical signs consistent with rabies at the time of exposure, or within 10 days following exposure in a rabies-infected area [17]. A case of confirmed exposure is a person who has had close contact (usually a bite or scratch) with a laboratory-confirmed rabid animal [15][16][17].

Data sources
Review of retrospective hospital case records. The initial dataset comprised all clinical cases of DBVs (n = 618; 460 in Buchanan, 85 in Gbarnga, and 73 in Voinjama) that were presented for treatment at the three selected reference hospitals during the first 6-year period, between January 1, 2008 and December 31, 2013, in the three cities ( Table 1). The 2014-2016 DBVs datasets were not collected in Liberia due to Ebola virus disease (EVD) outbreak during that period. Non-dog bite cases and case records of dog bite incidents that occurred outside each of the cities and their immediate environs, not exceeding the county boundaries were excluded from the final dataset used for this study. However, we utilized predicted cases of DBVs from deterministic model of projected data for 2014-2017 dog bite cases. Details about the deterministic model procedure is presented below.
For each DBV, information on the date of bite incident, gender, age, residential location, injured body part(s), the length of time that elapsed before victim showed up for anti-rabies prophylaxis, medications given, and treatment outcome were extracted from the records. The data generated were sorted on monthly basis for the 10-year period, from January 2008 to December 2017 and used to assess treatment compliance and time of presentation of DBVs for treatment, especially when the bite was assumed to have come from a rabid dog.
Household questionnaire survey on DBVs and care-seeking preference. A pre-tested questionnaire on care-seeking preferences of DBVs and their PEP compliance was used for household survey in the three cities, between January and August 2014. A total of 1,090 households (365, 375 and 350, respectively, in Buchanan, Gbarnga and Voinjama) were administered the questionnaire, just at the onset of the West African Ebola virus epidemic in Liberia. A respondent was provided by each household. Respondents were selected based on a three-stage population-representative stratified random sampling method. In each city, human dwellings (that is, stand-alone buildings) were randomly selected based on geographic size (that is, probability proportional to size). Dwellings were considered to contain multiple households if meals or living spaces were not shared. Hence, number of households in each dwelling was determined, and in-house identification number was assigned to each household. A simple random sampling method was then used to select a total of 1090 households enlisted and used for the study. These however were not integrated into the site names that were geospatially referenced to develop the Geospatial Observation-Linked Data Profile of Human Exposure to Suspected Rabid Dogs (GOLD-PHESRD).
Deterministic model for forecasting 2014-2016 missed DBV data. Due to the prevailing challenges of EVD epidemic in West African that resulted in severe humanitarian crisis of international concern [14,18], no human hospital records of dog bite victims in 2014-2016 were available to the authors at the county reference hospitals and at the Liberian NIH. Thus, we utilized predicted case of DBVs from deterministic model of projected data for 2014-2017 expected for Buchanan, Gbarnga and Voinjama. The data were previously developed and reported by co-authors of the present study [19]. A full description of this deterministic model is given in the full text of the paper online (https://www.researchgate.net/publication/ 302575115). At the resumption of reporting of DBV cases in January 2017, the Liberia Ministry of Health and Social Welfare provided dog bite data for 2017, as collated by the Liberian National Institute of Health (LNIH) which was established in 2016 after the EVD epidemic [14]. We then compared the model projected data for 2017 [19] with the actual cases reported to LNIH in 2017. The findings were evaluated and integrated into a continuous 10-year dataprofile of DBVs in the three cities with their respective county environs. The LNIH-enhanced reported data profile for 2017 DBVs comprised hospital-records of DBVs from 10 clinics/hospitals in Buchanan, including LGH; some 23 clinics/hospitals in Gbarnga, including PHNS; and 21 clinics/hospitals in Voinjama, including TMH. One health data audit and rabies laboratory confirmation. Evidence of confirmation of rabies infection among DBVs was sourced from the Liberian National Diagnostic Unit. However, this was limited to identification of clinical cases that were consistent with WHO case definition of rabies. The evidence included neurological signs, hydrophobia in dogs, and/or human death with history of dog bite. By 2017, the Central Veterinary Laboratory installed a fluorescent microscope in collaboration with Swiss Tropical and Public Health Institute and funding by the Global Alliance for Vaccines and Immunizations (GAVI). However, routine confirmatory diagnosis of rabies cases were not performed at that time as both the reagents and manpower for use of the equipment were just being deployed in that year [20,21].
A field survey and laboratory diagnosis of rabies from suspected human and dog specimens in Liberia, 2014-2016 had utilized RNA extraction technique, PCR detection and gene sequencing protocol [6] to confirm rabies virus infection in dogs and human exposures. The present authors have reported in an earlier paper [6], three rabies virus (RABV) strains isolated during the 2014-2016 survey in Liberia. The three RABVs sequences were accessioned in NCBI GenBank (MF765758, MH507336 and MH507337). However, the great majority of offending dogs and their bite victims were not tested. They remained rabies suspect as recommended by WHO classification for Liberia, being a recognized rabies endemic West African state [15]. Likewise, humans bitten by suspected rabid dogs were rabies exposure suspects

Conversion of site name to map point of DBVs presented to hospitals
The residential addresses of DBVs were available in most cases and were used to generate map points of street locations of DBVs. A combination of hand-held Global Positioning System (GPS, Garmin eTrex) and Google Earth Pro (Google, California, USA) were used to convert site names of DBVs to map points for all 618 victims enlisted in the study. Table 1 shows the distribution and proportion of 618 DBVs in the cities, including LGH, Buchanan, PHNS, Gbarnga, and TMH, Voinjama having 460, 85, and 73 DBVs, respectively, over the first 6-year period. The map points were saved in Microsoft Excel, version 2007 (Microsoft, Redmond, Washington), and then used for map visualization and spatio-temporal analyses, using Kulldorff's space-time scan statistics as previously described [22][23][24][25][26][27].

Statistical analyses and modelling
Descriptive and categorical analysis of DBVs. A descriptive map (Fig 2) of retrospective annual DBVs that were presented for treatment in Liberian cities and environs was drawn to visualize the geographic distribution of reported cases by county health authorities during the 10-year period, 2008-2017, using ArcGIS 10.1 (Environmental Systems Research Institute, ESRI, Redlands, CA) platform. Epidemiological profile of the victims was then compiled to summarize their ages, gender, time of care after bite (categorized into early, when less than 24 hours; or late, when more than 24 hours), site of dog-bite mark on the human body, initial treatment received, number of doses of rabies vaccine and/or number of doses of Equine Rabies immunoglobulin administered. The compliance of DBV to PEP treatment appointments was aggregated into: (i) completed; or (ii) defaulted. In addition, a descriptive statistics of treatment outcomes was computed under three categories, namely: (i) prevented; (ii) dead; and (iii) unknown (Tables 2-4). A multivariate matrix was computed to display the frequency distribution of age, gender, and seasonal frequency of dog bite cases in the three cities. Chi Square test was used to analyze the frequencies of each of the categorical variables. Statistical significance was determined at α 0.05 . DBVs were categorized based on site of bite mark, as bite to head, hands, legs, trunk (Table 5), and multiple sites.
In addition, a descriptive statistical summary of non-hospital treated DBVs was audited and presented in simple tables (Tables 6, 7 and 8), based on data generated from household questionnaire survey of DBVs. The tables summed victims that did and those that did not visit an orthodox health facility in Buchanan, Gbarnga and Voinjama cities, respectively within the 10-year period. Total DBVs presented at the hospitals were estimated based on the ratio of 2017 cases presented at each of the designated county level referral hospital per city, compared to 2017 enhanced reporting from several health facilities in each city/county [14].
Spatio-temporal scan statistics. Individual DBV was presented as map point in SaTScan 9.0 environment after they were first converted by site name to map point. A case file of DBVs with time (date) of occurrence in each city was created. Using Kulldorff permutation model [24], a space-time scan was used to analyze and identify the most likely space-time cluster in each city. The Kulldorff's space time scan statistic was preferred for each city based on its cylindrical window, in which the circular base represents the geographic area of the potential cluster of dog bite and the vertical arm represents the time span. The null hypothesis of space-time permutation scan assumes that the relative risk of occurrence of DBV is the same within the cylindrical window compared with outside it. The model controls for co-variants and eliminates preselection bias by not specifying a priori the observed set of cases within a cluster [22,23].
Spatial size of the scan window was limited to 50% of the total population per city. This was based on reported cases of dog bite from 2008 to 2013 (n = 460, 85 and 73 in Buchanan, Gbarnga and Voinjama, respectively), estimation for 2014-2017 [19], and actual reported cases in 2017. In each analysis, cluster radius of 1.0 km with time scale of one month was specified. Statistical significance of a cluster was determined by comparing the expected against observed number of cases of DBVs based on the null distribution obtained through Monte Carlo simulation with 999 replications and significance level was set as 0.05. The window with the maximum test statistic was considered the cluster least likely to be due to chance (primary  [28], and applied by Knobel et al. [29] and Tenzin et al. [30]. The model comprises a set of ten probability steps that incorporate the distribution of bite injuries on different parts of the body of victims, and the probability of a victim developing rabies, based on verified clinical records. The decision tree for determining the probability of rabies following the bite of a suspect rabid dog has been described by Tenzin et al. [30]. The probability of human deaths resulting from the bite of suspected rabid dog, PR, was estimated using the following model equation: where, P1 is the probability of suspected rabid dog being rabid; P2 is the probability of bite   Table 7. https://doi.org/10.1371/journal.pntd.0008957.t003 injury to the head; P3 is the probability of bite injury to the hands; P4 is the probability of bite injury to the trunk; P5 is the probability of bite injury to the legs; P6 is the probability of developing rabies following bite injury to the head; P7 is the probability of developing rabies following bite injury to the hands; P8 is the probability of developing rabies following bite injury to the trunk; P9 is the probability of developing rabies following bite injury to the legs; and P10 is the probability of receiving post-exposure treatment if bitten by a suspected rabid dog. The probability of developing rabies following the bite of a rabid dog to the head (P6), hands (P7), trunk (P8) and legs (P9) were taken to be 45.0%, 28.0%, 5.0% and 5.0%, respectively, as previously reported [28][29][30].
The total number of suspected human rabies cases per year, N, was estimated by: where, i is the incidence of suspected rabid dog bites per 100,000 per year; Q is the population at risk; and PR is the probability of a human death following a bite from a suspected rabid dog [28][29][30].   Table 8. https://doi.org/10.1371/journal.pntd.0008957.t004 The population distribution of the three selected cities and environs was calculated based on Liberia 2008 census (which remains the current, official report for the country), and estimated annual growth rate of 1.4% in Buchanan and environs, Grand Bassa County, 1.0% in Gbarnga and environs, and 1.3% in Voinjama and environs in 2008-2011 [13].  Head (Y t = 1.867+0.086t); Hands (Y t = 1.2+0.229t); Leg (Y t = 2.933+0.4t); Trunk (Y t = 3.067+0.171t) � Where multiple bites were sustained by a victim, the victim was assigned to one of the anatomical locations bitten that portends highest probability of developing rabies. https://doi.org/10.1371/journal.pntd.0008957.t005

Table 6. Proportion of hospital presented dog bite victims compared with non-hospital visiting victims used to compute normalized annual cases of dog bite victims in Buchanan City and environs, Grand Bassa County, Liberia, 2008-2017.
Year

Normalized annual count of dog bite victims, 2008-2017
Normalized annual cases of DBVs, 2008-2017 were computed from the proportion of cases presented at each county designated referral hospital compared to community survey of nonhospital visiting DBVs in each year. Tables 6, 7, 8 and 9 provide city-specific computation of  (Table 9). https://doi.org/10.1371/journal.pntd.0008957.g003 Table 9. Estimation of annual prevalence of human deaths due to rabies in three selected cities and environs, Liberia, 2008-2017.

Estimated annual human deaths due to suspected rabies, 2008-2017
An estimated total of 559 DBVs died of suspected rabies in the three cities and their environs during 2008-2017. Majority of DBVs in Buchanan, Gbarnga and Voinjama (47.3%, 58.8% and 46.6%, respectively) that sought medical care were presented more than 24 hours after sustaining dog bite injury (Tables 2, 3 and 4). The estimated annual DBVs for the three cities, from 2008 to 2017, are presented in Table 5. The corresponding parameters for estimating suspected rabies deaths in each of the three cities, based on probability decision tree and anatomical locations of dog-bite injuries are presented in Table 10. Annual human deaths estimation parameters were abbreviated as P1, P2, P3, P4, P5, P6, P7, P8, P9, and P10 (Table 10) Tables 6, 7 and 8, respectively. The estimated Mean Annual Human Deaths due to suspected rabies cases from 2008 to 2017 was 14.3±8.47, 35.5±23.25, and 6.1±3.21 per 100,000 population in Buchanan, Gbarnga and Voinjama cities and environs, respectively (Table 9).

Space-time clusters of dog bite victims
There were two significant space-time clusters in Buchanan. The most significant cluster occurred from May to December 2013. The second cluster occurred earlier, in June 2010, and lasted for only 30 days (p<0.013), with a radius of 0.32 km (Table 11; Figs 4 and 5). Four significant space-time clusters of DBVs were detected in Gbarnga City. One of the clusters had a duration of 31 days. In Voinjama, dog bite events were diffuse in space and time without a significant cluster. Notable about these clusters was that they were located within city environment, or in peri-urban areas, relative to other locations outside the clusters. This spatial pattern appear to indicate that persons living closer to county-specific reference hospital were more prone to present victims of dog bite to orthodox hospital nearby compared to people that lived farther away from the hospitals but still within the same County.

Discussion
The study was conducted to improve dog bite victim survey in three selected cities of Liberia using an explanatory sequential mixed-method study. This included retrospective clinical case review of dog bite victims, household questionnaire survey of victims' care-seeking preferences and multi-criteria analysis of their post-exposure prophylactic treatment compliance. All these were done to estimate both the absolute number of annual human deaths due to suspected rabies cases and the respective prevalence of cases in the three selected cities of Liberia. Thus, we conducted city-specific data audit of dog-bite victims in Buchanan, Gbarnga and Voinjama cities and used a standardized clinical decision tree model, earlier proposed [30] to compute the annual human deaths in each of the three cities. By accounting for previously unreported dog-bite victims and suspected rabies cases, a more accurate estimate of DBVs and suspected rabies cases were determined as the true burdens of rabies in each of Buchanan, Gbarnga and Voinjama cities, along with their respective county environs in Liberia, 2008-2017. We here reported the mean yearly prevalence of dog bite victims as 179±106.82, 393±257.85 and 76.9  Thus, the two critical government agencies that were needed to enlist their employees who were willing, able, and available for participation in human-animal disease survey, laboratory detection and reporting in Liberia supported the project. Two personnel, one each from MoHSW and MoA enrolled with the CCPZ to embark on a postgraduate program in humananimal disease surveillance at the University of Ibadan, Nigeria. The two sectors thereby worked with the academia to pool together essential data from animal and human health sectors. Questionnaire survey were jointly designed and conducted to gather additional field and archival specimens of suspected dog and human rabies cases for confirmation of rabies at the grassroots level in the three Liberian cities from 2012 to 2017. All secondary data used in this study were provided by the Liberia Ministry of Health and Social Welfare. All the data were shared with personnel of the Central Agricultural Research Institute (CARI) and academia in one health partnership.
Gbarnga city had the highest (approximately 36 deaths per 100,000 population) of annual human deaths due to suspected rabies cases out of the three cities. The 36 cases were much lower than what was reported [2] for Monrovia, the capital city of Liberia, being 155 AHD due to rabies, about the same period. The Gbarnga rabies burden was comparable to reported estimate for Plateau State, Nigeria, and Abidjan in Cote d'Ivoire, being an average annual human deaths of 40 per 100,000 population, due to rabies [1,31,32]. Gbarnga contrast with Monrovia is to be expected, since there is a high human population density difference between the two Liberian cities, with Monrovia, having about ten times the population density of Gbarnga [12,13]. Yet, the situation in Gbarnga portends a relatively high prevalence of suspected rabies deaths when compared to Plateau State of Nigeria. In each of the 10-year chronological prevalence, the annual human deaths due to suspected rabies cases in Buchanan was higher than Voinjama city, but comparable to other municipalities in West Africa, including Techiman Municipality in the Middle Belt of Ghana [33].
The two space-time clusters (Table 11) of dog bite victims identified in Buchanan and the four in Gbarnga city (based on hospital-presented cases only) were proximal to the locations of each reporting hospital facility during 2008-2013, whereas rabies-susceptible dogs were present across the entire city and county landscapes as reported by Jomah et al. [34] based on rabies neutralizing antibody survey, 2015. This finding corroborates the tendency of dog bite victims to self-medicate or present to traditional caregivers in rural areas than to conventional hospitals among low income populace in Liberia [2]. A hospital-based surveillance system would usually miss this segment of the public, except where additional survey is done to capture non-hospital visiting sub-group. Poverty and related low-or non-formal education have been associated with non-hospital care-seeking tendencies among people in Liberia [8,9] in general and dog bite victims in particular [2]. Often the victims were unable to afford the cost of rabies vaccines for post-exposure prophylaxis, and would naturally seek care from traditional caregivers. This is similar to the findings in Ilorin, a city of Nigeria with notable traditional sector [35], as well as in the slum areas of Monrovia in Liberia [2,36], and slums of Abidjan in Cote d'Ivoire [32].
The detection of cluster radii lower than 1.0 km and time scale larger than a month, which were different from those earlier specified at statistical analysis were due to a priori guide to be analysis. It was notable that shortly after Ebola virus disease epidemic ended in 2017, the number of hospitals included in the integrated disease surveillance system had remarkably increased from three hospitals only (one each from the three selected counties, of Grand Bassa, Bong and Lofa, respectively) as obtained during 2008 to 2014, to a minimum of ten hospitals/ clinics per county in each of the 15 counties of Liberia by 2017 [14]. This transformation was facilitation by the WHO following Ebola virus disease outbreak that attracted global attention and support to Liberia [14,20,21].
We have exercised high caution in our claims about preliminary laboratory confirmation backing of the secondary data sourced form hospital records and the MoHSW data because of the retrospective nature, 2008-2017 that did not have accompanying laboratory confirmation of each of the rabies suspected cases. However, rabies suspect status and the deterministic model of dog bite victims here derived conformed to WHO 2004 case categorization. We have reported that one of the RABV strains (MF765758) isolated in Liberia clustered with China Lineage 2 of dogs. This strain co-circulated with two other strains of Africa lineages [6]. Five other RABV strains, two each first isolated in France (GU992321, GQ918139), and India (KF535200, KF535201) and one from China (DQ875050) have since been added by Zaho et al. [37], to the list of RABV strains that were 99% similar to Liberia MF 765758 strain [6], thus, being phylogenetically of the same lineage, China Lineage 2 (cosmopolitan). These indices have added value to closing the critical data gap in the 10-year prevalence of suspected rabies cases [38] in the three cities. The sequential mixed-method approach here used could be adaptable to other counties of Liberia and West African countries [10,11,25,[38][39][40].
These authors have learned that a horizontal [41] and a vertical line of interaction with stakeholders in rabies surveillance are essential to achieve the one health action in Liberia here reported. The achievement would be sustained through the same lines of engagement and research coordination. As a result, a forum for Rabies in West Africa (RIWA), which was inaugurated in Ibadan, Nigeria in December 2012 [42] had administratively incorporated a board of trustees for a Society for Rabies in West Africa in November 2017. The new incorporation will enhance a sustainable coordination of the improvement now achieved in rabies surveillance, by further facilitating partnership network with the Central Veterinary Laboratory in Liberia and with regional and global veterinary laboratory diagnosticians, such as from the University of Ibadan CCPZ, the Global Alliance for Rabies Control, the OIE-FAO-WHO tripartite and the Swiss Tropical and Public Health Institute. In addition, fostering partnership with the World Association of Veterinary Laboratory Diagnosticians, would enhance quality in veterinary diagnostic services as required by OIE [43], and supported by the Pan-African Association of Veterinary Laboratory Diagnosticians.
The present study presented a 10-year chronological estimates of annual prevalence of human deaths due to suspected rabies cases in Buchanan, Gbarnga and Voinjama cities of Liberia with their local environs, 2008-2017. This information will be useful for planning more effective control and stepwise actions to achieving the goal that no human should die of rabies by 2030.