Figures
Abstract
Introduction
In the continental US, four terrestrial mammalian species are reservoirs for seven antigenic rabies virus variants. Cross species transmission (CST) occurs when a rabies virus variant causes disease in non-reservoir species.
Methods
This study analyzed national surveillance data for rabies in terrestrial mammals. The CST rate was defined as: number of rabid non-reservoir animals/number of rabid reservoir animals. CST rates were analyzed for trend. Clusters of high CST rate counties were evaluated using space-time scanning statistics.
Results
The number of counties reporting a raccoon variant CST rate >1.0 increased from 75 in 1992 to 187 in 2011; counties with skunk variant CST rates >1.0 remained unchanged during the same period. As of 2011, for every rabid raccoon reported within the raccoon variant region, there were 0.73 cases of this variant reported in non-reservoir animals. Skunks were the most common non-reservoir animal reported with the raccoon rabies variant. Domestic animals were the most common non-reservoir animal diagnosed with a skunk rabies virus variant (n = 1,601). Cross species transmission rates increased fastest among domestic animals.
Conclusions
Cross species transmission of rabies virus variants into non-reservoir animals increases the risk of human exposures and threatens current advances toward rabies control. Cross species transmission in raccoon rabies enzootic regions increased dramatically during the study period. Pet owners should vaccinate their dogs and cats to ensure against CST, particularly in regions with active foci of rabies circulation. Clusters of high CST activity represent areas for further study to better understand interspecies disease transmission dynamics. Each CST event has the potential to result in a rabies virus adapted for sustained transmission in a new species; therefore further understanding of the dynamics of CST may help in early detection or prevention of the emergence of new terrestrial rabies virus variants.
Citation: Wallace RM, Gilbert A, Slate D, Chipman R, Singh A, Cassie Wedd, et al. (2014) Right Place, Wrong Species: A 20-Year Review of Rabies Virus Cross Species Transmission among Terrestrial Mammals in the United States. PLoS ONE 9(10): e107539. https://doi.org/10.1371/journal.pone.0107539
Editor: Wanda Markotter, University of Pretoria, South Africa
Received: March 28, 2014; Accepted: August 20, 2014; Published: October 8, 2014
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.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. Data were collected through routine, government, animal rabies surveillance. Data are available through written request to the author and after signature of approved data release agreement. Requests may be sent to the corresponding author.
Funding: These authors have no support or funding to report.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Rabies has one of the highest known infectious disease case fatality rates, with an estimated 55,000 deaths per year occurring primarily in developing countries where canine rabies has not been controlled [1]. Historically, rabies was responsible for hundreds of human deaths each year in the US, prior to the development of laws promoting responsible pet ownership, canine rabies vaccination campaigns, and increasing availability of effective rabies post-exposure prophylaxis (PEP) [2]. By 1970 the majority of human rabies cases resulted from exposures to rabid wildlife rather than dogs. In paradoxical fashion, by the mid -1970's, when the rabies burden in the US was at presumptive post-European colonization lows, one of the largest recorded wildlife disease epizootics was beginning in rural West Virginia [2], [3].
In 1977, translocation of raccoons from Florida to West Virginia for hunting purposes ignited a rabies epizootic in the local rabies-naïve raccoon populations [4]. Over the next 20 years raccoon variant rabies spread to all mid-Atlantic and northeastern states [2]. In stark comparison to the rabies epidemiology prior to 1970, over 90% of rabies cases are now reported from wildlife species and over 60% of these cases are reported from states that have enzootic circulation of the raccoon rabies virus variant [3], [4]. As evidence from this epizootic event, the introduction of novel rabies virus variants into naïve populations can have longstanding consequences on wildlife populations and human disease prevention. While domestically acquired human rabies deaths are at an historic low and canine rabies virus has been eliminated in the US, the risk of rabies transmission from wild animals remains substantial. More than four million animal bites are estimated to occur in the US each year [7]. Only a fraction of these are reported to public health officials [5], [6]. Despite poor reporting of animal bites, an estimated 35,000 persons (11/100,000 persons) still receive PEP each year due to potential rabies virus exposures [2], [10], [11].
In the continental US, wildlife species in the orders Carnivora and Chiroptera are responsible for enzootic circulation of rabies virus. These include: foxes (Vulpes lagopus and Urocyon cinereoargenteus), raccoons (Procyon lotor), skunks (Mephitis spp), and bats (multiple species) [7]. Seven distinct antigenic rabies virus variants are associated with these four terrestrial species [4]. In addition, more than 10 rabies virus variants are associated with as many species of bats [8]. Each distinct rabies virus variant is maintained within a specific animal reservoir species [9]–[12]. The geographic boundaries of these terrestrial rabies virus variants are dynamic and well described. Nearly all rabid terrestrial mammals are infected with the predicted terrestrial rabies virus variant associated with the geographic region where infection occurred, and very few cases of terrestrial rabies occur in geographic regions with only bat rabies virus variants [10], [13], [14].
Despite the apparent host adaptation and affinity displayed by rabies virus, all mammals are susceptible to this disease. Cross species transmission (CST) occurs when a rabies virus variant adapted to transmission in a specific reservoir animal is transmitted to a non-reservoir animal. Cross species transmission rarely results in successful propagation within a non-reservoir population, referred to as a host shift event [12], [15], [16]. However, the diversity of the seven terrestrial rabies virus variants in the US are apparently due to sustained host shifts from canine rabies virus variants introduced during colonization of the Americas as well as extant chiropteran variants [3]. More recently, repeated CST from bats into skunks and gray foxes in Northern Arizona has resulted in a focal region of sustained bat variant transmission within these terrestrial mammals [11]. Another recent study revealed possible molecular evidence of independent circulation of the California skunk variant in gray foxes following repeated CST events [17]. The impact of these and future potential host shift events are difficult to predict; however, it is possible that such events could impede or compromise current rabies management efforts and jeopardize the health and safety of people and their companion animals. Each instance of CST is an opportunity for viral adaptation resulting in sustained propagation in a non-reservoir species, and CST events are reported in the thousands each year [14]. Currently, public health surveillance is the primary approach for detecting potential host shift events. This study uses routinely collected surveillance data to describe spatiotemporal dynamics of rabies transmission among terrestrial mammals and serve as motivation for new evaluations and the development of techniques to better detect and predict areas at high risk of CST and subsequently where rabies virus host shift events may be more likely to occur.
Methods
Data Source
Rabies in animals is a nationally notifiable condition in the U.S. Nearly all animals tested for rabies are submitted due to a bite or other concerning interaction with a human or domestic animal. A small proportion of animals (∼5% annually) are reported through enhanced rabies surveillance involving wildlife that have not exposed humans or domestic animals [14]. National reporting of rabid animals has been conducted since 1938. The CDC program maintains a national database of aggregate reported animal rabies cases by county since 1990. Individual level data including reports for all animals submitted for rabies diagnosis was available since 2006 as described previously [14]. Regional analysis was performed using Health and Human Services regional designations. County-level human population-based rates were calculated from U.S. Census Bureau population data, 2010 census [18].
Data Selection and Categorization
Reports of all rabid terrestrial mammals within the continental US were compiled from 1990–2011. Rabid bats and humans were excluded from the study, as were data from Alaska, Hawaii, and Puerto Rico. States where only bat rabies virus variants have been reported were also excluded from the analysis (Washington, Oregon, Idaho, Nevada, and Mississippi). Rabies virus variant typing was not available for the majority of reported rabid animals. However, the variant associated with a terrestrial rabies case was assumed based on the geographic location where the animal was reported (e.g. a rabid raccoon in Texas was assumed to be infected with a skunk rabies virus variant). This assumption is based on previously published reports that have found that nearly all rabid terrestrial mammals were infected with the terrestrial rabies variant associated with the geographic region in which the case was reported [10], [13], [14]. This assumption was confirmed by analysis of species and rabies virus variant for reported cases where rabies virus variant information was available (2007–2011). All rabid terrestrial mammals within the defined raccoon variant region were assumed to have the raccoon rabies virus variant [4], [10], [14], [19]. Terrestrial mammals diagnosed with rabies in Arizona and a small region of west Texas were assumed to have either a skunk or fox rabies virus variant, as both variants are present in these regions. Rabid terrestrial mammals in all other counties were assumed to have had one of three skunk rabies virus variants. The three skunk virus variants were analyzed as one group to accommodate instability in predicted geographic boundaries among two adjacent variants. Non-reservoir animals were grouped as raccoons (P. lotor), skunks (all species), foxes (all species), domestic animals (cats and dogs), cattle, and other terrestrial mammals (e.g. bobcats, non-cattle ungulates, lagomorphs, rodents, and others).
Counties that reported fewer than three rabid reservoir animals during a 5-year period were excluded from analysis to control for jurisdictions with low surveillance activity. An epizootic of a Mexican canine rabies virus variant occurred in coyotes in southern Texas counties from 1990–1997, therefore rabid coyotes in select southern Texas counties in this time frame were excluded from analysis [20], [21]. In Tennessee, 17 counties were located in a region where the skunk and raccoon variant boundaries overlap. For these counties annual counts of rabid skunks and raccoons were compared and the reservoir was defined by the more frequent of the two potential reservoir species.
Rate Calculations
The reported rabies CST rate was defined as the rabid non-reservoir animals/rabid reservoir animals. These rates were calculated at the county and HHS regional level, annually and in 5-year aggregate time periods. Annual rates of animal submissions for rabies testing were calculated per 100,000 human population, 2007–2011; the only years when reliable submission data were available. To evaluate the impact of bias, a linear regression model was used to evaluate potential relationship between total submissions of reservoir animals for rabies diagnosis on the CST rate. All comparisons were evaluated at α = 0.05.
Temporal and Spatial Analysis
Trends in CST rates from 1990–2011 were evaluated using a Monte Carlo permutation model, which detected up to three time intervals where the trend in CST changed significantly [22]. The annual percent change in CST and accompanying 95% confidence interval were reported for the most recent time interval trend identified in the model. Cross species transmission rates were mapped by county by 5-year intervals [23]. A spatial-temporal scan statistic was used to identify county level clusters with a higher than expected CST rate in the raccoon rabies enzootic region from 2007–2011.The scan statistic was run and presented at multiple levels (2, 4, 8, 12, 25, and 50 percent) for the maximum spatial cluster size parameter (MSCS; i.e. the maximum percentage of the population at risk included in a scanned cluster) to increase sensitivity in relation to this parameter setting [24].
Results
Annual Trends in Rabies Spillover
Raccoon Variant Region.
A total of 67,058 rabid raccoons and 30,876 other animals were reported with raccoon rabies virus variant from 1990–2011. The CST rate ranged from 0.31 in 1990 to a high of 0.73 in 2011, representing a 138% overall increase (Table 1). Three time intervals with distinct CST trends were identified. From 1990–1993, there was no change in the annual CST rate. From 1993–1998, the CST rate increased at 15.9% each year (95% CI: 7.5–25.0), but plateaued again from 1998–2011 (Figure 1). Skunks were the most common non-reservoir animal reported with the raccoon rabies virus variant (n = 16,600). Overall, the annual CST rate of raccoon rabies into skunks increased 118.0% from 1990–2011 (0.16 to 0.35 respectively). However, CST of raccoon rabies into skunks did not change from 1999–2011 (95% CI: −3.2–1.4) (Table 1).
Skunk Variant Region.
A total of 19,247 rabid skunks and 3,778 other animals were reported from skunk rabies enzootic regions from 1990–2011. Overall, there was a 40.3% decrease in the CST rate of skunk rabies from 1990–2011 (0.28 to 0.17, respectively), with an average decrease of 2.0% (95% CI: −3.6–−0.4) per year (Figure 2).Domestic animals were the most commonly reported non-reservoir animal infected with skunk rabies (n = 1,601) (Table 1). Cross species transmission from skunk to raccoon was a rare event with only 148 occurrences reported from 1990–2011.
Raccoon and Skunk Variant CST Rates, 2007–2011
Reported CST of raccoon rabies into skunks was 18.7 times more frequent in comparison to skunk rabies transmitted into raccoons (95% CI 14.2–24.5) (Table 2). Transmission of the raccoon variant into foxes was 8.1 times more frequent when compared to skunk rabies virus variants (95% CI 6.3, 10.5). Transmission of the raccoon variant into domestic animals was 1.5 times more frequent than skunk variant transmission to domestic animals (95% CI 1.3, 1.7). In cattle, CST of the raccoon variant occurred less frequently than skunk variants (RR = 0.5, 95% CI 0.4, 0.6). Overall, CST of the raccoon variant was 3.7 times more frequent than skunk variants during the study period (95% CI 3.4, 4.0).
Spatial Analysis
Raccoon variant region.
During 1992–2011, 741 raccoon variant counties met the criteria for inclusion and reported an average of 6.2 rabid terrestrial animals each year. The CST rate doubled during the study period to 0.62 non-reservoir animals diagnosed for every rabid raccoon for time period 2007–2011.During 2007–2011, 187 of 592 (31.6%) counties reported a rate ≥1.0, and increase from 75 counties during 1992–1996 (Figure 3). From 1990–1999, the CST rate trend in the Northeast region increased at 23.8% per year (95% CI 17.1–38.0). This trend changed from 1999–2011, when the CST rate in the Northeast decreased 3.9% per year (95% CI −7.2–−0.6) (Table 1). In its most recent trend time period, 1994–2011, the spillover rate in the Southeast region increased 5.8% per year (95% CI 5.0–6.6). From 1990–2011, the CST rate increased at 2.7% per year in the mid-Atlantic region (95% CI 2.0–3.5).
Skunk Variant Region.
During 1992–2011, 754 skunk variant counties met the criteria for inclusion and reported an average of 1.4 rabid terrestrial animals each year. During 1992–1996, 0.22 other terrestrial animals were diagnosed rabid for every reported rabid skunk, and 35 of 392 (8.9%) counties reported a CST rate ≥1.0. During 2007–2011, the CST rate fell to 0.18 and only 20 of 324 (6.2%) counties reported a CST rate ≥1.0 (Figure 4).
In the Southeastern United States, the skunk rabies CST rate trend increased 13.6% per year from 2001–2011 (95% CI 1.9–26.7) (Table 1). In the most recent trend period in the southern United States from 2001–2011, the skunk rabies CST rate increased at 7.0% per year (95% CI 3.0–11.1). No significant changes in the trend or overall CST rate were identified for the Midwest and West regions.
Space-Time Cluster Analysis.
Three significant clusters of rabid non-reservoir animals were identified controlling for the number of reported rabid raccoons in the space-time scan statistic at a MSCS of 50% (Figure 5). The two clusters centered on New Hampshire and New York appear to remain fairly stable at lower MSCS levels. However, the risk levels appeared to be heterogenous across the large cluster centered on North Carolina. Additional smaller clusters were identified within this region starting at a MSCS level of 12%.
Validation of the Variant Assumption
Presence of variant typing data were examined for all 15,359 animals analyzed in this study between 2007–2011 (Figure 6). The study-defined raccoon variant region reported 13,044 rabid animals, of which 1,212 (9.3%) were variant typed. Of these 1,212 variant typed animals, 1,203 (99.3%) were infected with the predicted reservoir virus variant. The nine animals which did not have the predicted rabies virus variant were all found to have been infected with a skunk variant. These nine animals were from Tennessee, a state in which both skunk and raccoon rabies virus variants are present. The study-defined skunk variant region reported 2,315 rabid animals, of which 651 (28.1%) were variant typed. Of these 651 variant typed animals, 634 (97.4%) were infected with the predicted reservoir virus variant. Twelve of the 17 animals predicted to have skunk variant rabies virus were actually infected with a raccoon rabies virus variant. All 12 were from the state of Tennessee. The remaining five animals (0.8%) were infected with a bat variant. Of all variant typed animals, 98.6% had the study-defined, predicted rabies virus variant, 1.1% did not have the predicted variant but were from a state in which multiple variants are present, and 0.3% of animals were infected with a non-terrestrial variant.
Submission of Samples for Rabies Testing
Data on submission of rabies suspect animals for testing were available for 2007–2011. During this time period, 53,877 raccoons and 137,517 other terrestrial animals were submitted for rabies diagnosis in the raccoon variant region (20.7% and 5.0% rabid, respectively). The average annual number of submissions for rabies diagnosis was 57.9 per 100,000 persons residing in the raccoon variant region. From 2007–2011, 7,485 skunks and 45,530 other terrestrial animals were submitted for diagnosis from the skunk variant regions (44.5% and 1.3% rabid, respectively). The average annual number of submissions for rabies diagnosis was 67.9 per 100,000 persons residing in the skunk variant region.
A linear regression model revealed a relationship between CST and submission of raccoons for rabies diagnosis in the raccoon variant region (P = 0.012) (Table 3). However, changes in the number of raccoons submitted had little impact on the CST rate (β = −0.00018). A similar trend was identified in the relationship in the skunk variant region, but with a slightly greater impact on CST rates when the number of skunk submissions is increased or decreased (β = −0.0015, P = 0.008).
Discussion
Interspecies rabies transmission can be important in the process of virus adaptation and perpetuation in new hosts [3]. However, from a public health perspective, the more urgent concern is the role CST plays in potentially increasing the threat of rabies virus exposure to humans and domestic animals. The burden and frequency of CST appears to have increased in recent decades and was dramatically higher where the raccoon rabies virus variant is enzootic. This study provides epidemiologic evidence that the raccoon rabies virus variant may be more adept at infecting skunks than previously expected. Whether this is due to viral characteristics or host-host interactions requires further study. Furthermore, this study provides a unique analytic method for the prediction of rabies virus variant host shift events. Routine implementation of this analytic method to annual or real time data could be considered for incorporation into rabies surveillance programs.
Changing CST Trends
The rate of rabies virus transmission to non-reservoir animals was dynamic over the study period. Perhaps the most dramatic trend observed was the increase in CST of the raccoon rabies variant in the northeastern United States. This increase is most likely explained by the epizootic spread of this rabies virus variant into a previously rabies naïve region. The raccoon rabies variant did not expand into the Northeastern-Atlantic states until the mid-1980's, largely reaching its full northern distribution by 2000 [25]. Prior to the movement of rabies into this region, the raccoon population was rabies-naïve and populations had not experienced decimation due to the epizootic. During the epizootic, interspecies contact with the aberrantly high number of rabid raccoons likely accounts for a sizeable portion of the increasing CST rates observed in these northeastern states. In comparison, the second largest increasing trend in the CST rate was observed in Mid-Atlantic states where the raccoon rabies variant had been established for nearly two decades, and the smallest increase was found in Southeastern states where the variant has been established since the 1940's [26]. In contrast, the three skunk rabies virus variants have been established in the United States since at least the early 1800's; all three variants showed a decline in CST rates [27].
The apparent relationship between the CST rate and length of enzootic status, suggests that the increases in CST in the mid-Atlantic and Northeastern states might be expected to stabilize over time. However, interpretation of this apparent association is difficult, as many potential confounding factors are not collected during routine rabies surveillance efforts. Regional differences in animal population densities and species diversity have the potential to alter disease transmission dynamics. Therefore, ecological studies focusing on inter and intra-species contact rates may help explain the differences in the observed CST trends. Climate change, specifically drought, may also affect CST rates; regions experiencing drought may drive animal populations together as water sources become scarce. Likewise, changes in land use characteristics, such as deforestation or human encroachment into wildlife habitats may also result in changes in animal behaviors. Studies focusing on changing environmental features may help to explain regional differences in rabies transmission dynamics. This analysis describes general trends observed over the past 20 years and is meant to serve as a guide for more focused studies, as described above, to aid in the control of rabies in wildlife reservoir species.
Raccoons: Super Spreaders?
Raccoons were four times more likely to transmit rabies to other species when compared to skunks. All animal groups, with the exception of cattle, were more likely to be infected with the raccoon variant. Host and environmental factors which facilitate CST are largely unexplored. However, certain behavioral and anthropogenic factors are well described in skunks and raccoons which may explain some of the differences observed in the CST rates of these two species.
Transmission of diseases such as rabies that require direct contact are heavily influenced by the density of susceptible animals and their contact rate with infected animals [28]. Raccoon studies have recorded population densities of 1–250/km2, with higher densities occurring in urban settings [29]–[35]. These studies have also found frequent interspecies interactions. Fewer studies have examined skunk populations, but generally have reported lower population densities of only 3–6/km2 and discrepant reports of the impact of urban environments on population dynamics [36].
Factors other than population density, such as predation or interspecies interactions at communal resources, play integral roles in disease transmission dynamics [37]. Inquisitive species, such as raccoons, and territorial species, such as foxes, may be more likely to physically confront a rabid animal, resulting in increased disease transmission [36], [38]. Raccoons are opportunistic omnivores and prefer den sites near bodies of water and often prosper in high food resource areas such as suburban environment where human refuse, pet food, and other sources of reliable sustenance are likely to be abundant. In these settings, congregations of multiple species at feeding sites are commonplace [39]. Raccoons also commonly utilize latrine sites, which may facilitate interspecies disease transmission, as studies have found that up to 14 different mammalian species may frequent these sites [40]. Skunks, in contrast, typically have lower population densities and are considered an aposematic species; pelage and chemical spray are used as a warning to other species to stay away, which may account for reduced interspecies interactions and thereby lower levels of CST [41], [42]. Aposematic species may be accustomed to being avoided, rather than relying on fight or flight, and therefore less likely to flee when a rabid animal approaches. Skunks also frequently share den sites with raccoons [43]. These ecological and behavioral traits may play a role in their high rate of raccoon variant CST. However, interpretation of these interactions are complicated when the animal in question has rabies, as these animals no longer conform to the behavioral customs attributed to the species [12].
Cross species transmission of rabies among terrestrial animals is a complex interaction that likely depends on animal susceptibility to the virus, animal population densities, animal behavior, niche overlap, landscape characteristics, human population distribution, environmental conditions, and other factors. These multivariate associations can be difficult to evaluate over large geographic areas, even within a single virus variant region. Ecological studies of interspecies rabies transmission should focus on the geographic clusters of high CST rates to maximize detection of CST-related factors. Spatial models should be developed to evaluate the association of these factors with CST rates and identify significant factors that might provide insights which could be used to better understand local epidemiology of rabies and possibly targeted for interventions.
One Species, One Variant?
Classically, rabies virus variants have been associated with a single species which maintains the variant in nature by conspecific transmission. Cross species transmission of the raccoon variant was unexpectedly high, and the most commonly reported non-reservoir species were skunks. There are multiple theories for increased transmission of raccoon variant into skunks including aforementioned behavioral and ecological patterns that would cause increased interactions between the species [36], [37], [39]–[44]. However, it is not well described as to why these differences would result in increased transmission from raccoon to skunk, but not the inverse. There are potential virus virulence and host susceptibility differences that could affect transmission and result in increased transmission in only one direction, although these are not well studied either [38], [44]–[46].
One potential theory for the unusually high rate of raccoon variant rabies diagnosed in skunks is that this virus has already undergone a host shift into the skunk population. This would artificially inflate the raccoon variant CST rates observed in this study [12], [16]. Propagation of new rabies virus variants occur, presumably on a rare basis, when a virus adapts to a new species resulting in a host-shift and independent circulation in a new species. Over time the virus can be recognized as an independent variant through molecular or antigenic characterization. Molecular analysis of recent host shift events from bats to local mesocarnivores (skunks and foxes) suggested that the rabies virus in those events may have been genetically competent for a host shift prior to transmission (pre-shift adaptation) rather than undergoing genetic adaptation in the new host [11]. Contact rates between bats and mesocarnivores are relatively low, suggesting that transmission of the virus was likely restricted in one direction (i.e. bat to mesocarnivore). However, where contact rates in both directions between two species are relatively high (e.g. skunks and raccoons), it would seem plausible that a pre-shift adapted virus may allow for routine bi-directional transmission between species as well as conspecific circulation. Maintenance of the rabies virus variant between two hosts may prevent selective pressure on the virus complicating detection of a distinct rabies virus variant associated with a specific species. This theory would need to be tested, likely using both mathematical modeling and deep sequencing of viruses from raccoons and skunks. A comparison of high CST regions to low CST regions within the raccoon rabies virus territory might be an effective comparison for evaluating molecular differences in the rabies virus.
Limitations: Calculating CST Rates from Passive Surveillance
Nationally, rabies virus variant typing is prioritized among non-reservoir animals or locations where there is increased concern for new variant introduction. Variant typing is an added cost and burden for public health programs, and has no direct treatment impact for the human or domestic animal that was exposed. Therefore, rabies virus variant typing is performed on only a subset of samples, nationally. A critical assumption of this study was that all terrestrial animals were infected with the predicted rabies virus variant based on geographic epidemiology. The validity of this assumption is supported by the finding that nearly all variant typed animals were infected with the predicted study-defined variant, as has been reported in the literature [10], [13], [14] (Figure 6).
Public health testing of animals for rabies is directly related to the animal's probability of interacting with people or their companion animals. Differing state and local health department policies may also impact and bias animal submissions. Raccoons are often found at high population densities close to human habitats, which could potentially account for higher rates of submission for rabies testing [29]–[33], [47]–[49]. Contrary to this reasoning, our study found that on a per-capita basis, animal submission rates were similar between raccoon and skunk regions. Simple linear regression analysis indicated that the number of submissions of reservoir animals (skunk or raccoon) has little impact on CST rate. While there does appear to be some bias introduced due to surveillance activity it did not appear to be differential based on rabies virus variant and was of low magnitude, suggesting that submission and testing bias may not have a significant impact CST rates.
Conclusion
Rabies virus host shift events could threaten the rabies prevention success achieved in the past several decades in the United States. Traditionally, host shift events have been detected through astute local health programs. Given the consequences of delayed host shift recognition, new surveillance methods must be developed to rapidly identify potential high risk areas. This study provides two such methods for further exploration; monitoring of CST trends and evaluation of high CST areas. The development of the unique CST analysis has provided epidemiologic support to the theory of independent skunk to skunk transmission of the raccoon rabies virus variant. The development of algorithms to analyze areas at high risk for CST could be incorporated into standard reporting systems to raise awareness for host shift events among relevant health departments, thereby improving early detection of such events. Dogs and cats may have frequent opportunities for encounters with rabid animals and remain a critical barrier to human rabies exposure. The finding that CST rates increased in cats and dogs is a reminder of the importance of maintaining current vaccination status to protect animal health and prevent human exposure. Currently available animal rabies vaccines are effective against all known rabies virus variants in North America. The higher CST rate in the raccoon variant region highlights the burden of rabies in this region and increased risk even when encountering non-reservoir species. Each CST event has the potential to result in the establishment of a new reservoir species; therefore advancing our understanding of the dynamics of CST may help to prevent the emergence of new terrestrial rabies reservoir species.
Acknowledgments
Disclaimer: The findings and conclusions in this report are those of the authors and do not necessarily represent the views of their institutions. The authors thank the working group of Research and Policy for Infectious Disease Dynamics program of the Science and Technology Directorate, Department of Homeland Security.
Author Contributions
Conceived and designed the experiments: RMW AG JB. Analyzed the data: RMW AG JB AS CW RC DS. Wrote the paper: RMW AG JB AS CW RC DS.
References
- 1. Knobel DL, Cleaveland S, Coleman PG, Fevre EM, Meltzer MI, et al. (2005) Re-evaluating the burden of rabies in Africa and Asia. Bull World Health Organ 83: 360–368.
- 2. Nettles VF, Shaddock JH, Sikes RK, Reyes CR (1979) Rabies in translocated raccoons. Am J Public Health 69: 601–602.
- 3. Velasco-Villa A, Reeder SA, Orciari LA, Yager PA, Franka R, et al. (2008) Enzootic rabies elimination from dogs and reemergence in wild terrestrial carnivores, United States. Emerg Infect Dis 14: 1849–1854.
- 4. Blanton JD, Dyer J, McBrayer J, Rupprecht CE (2012) Rabies surveillance in the United States during 2011. J Am Vet Med Assoc 241: 712–722.
- 5.
Freer L (2004) North American wild mammalian injuries. Emerg Med Clin North Am 22: : 445–473, ix.
- 6. Kizer KW (1979) Epidemiologic and clinical aspects of animal bite injuries. JACEP 8: 134–141.
- 7. Blanton JD, Palmer D, Dyer J, Rupprecht CE (2011) Rabies Surveillance in the United States during 2010. JAVMA 239: 11.
- 8.
Constantine DG, Blehert DS, National Wildlife Health Center (U.S.) (2009) Bat rabies and other lyssavirus infections. Reston, Va.: U.S. Geological Survey. xi: , 68 p. p.
- 9. Streicker DG, Turmelle AS, Vonhof MJ, Kuzmin IV, McCracken GF, et al. (2010) Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science 329: 676–679.
- 10. McQuiston JH, Yager PA, Smith JS, Rupprecht CE (2001) Epidemiologic characteristics of rabies virus variants in dogs and cats in the United States, 1999. J Am Vet Med Assoc 218: 1939–1942.
- 11. Kuzmin IV, Shi M, Orciari LA, Yager PA, Velasco-Villa A, et al. (2012) Molecular Inferences Suggest Multiple Host Shifts of Rabies Viruses from Bats to Mesocarnivores in Arizona during 2001–2009. PLoS Pathog 8: e1002786.
- 12. Guerra MA, Curns AT, Rupprecht CE, Hanlon CA, Krebs JW, et al. (2003) Skunk and raccoon rabies in the eastern United States: temporal and spatial analysis. Emerg Infect Dis 9: 1143–1150.
- 13. Oertli EH, Wilson PJ, Hunt PR, Sidwa TJ, Rohde RE (2009) Epidemiology of rabies in skunks in Texas. J Am Vet Med Assoc 234: 616–620.
- 14. Dyer JL, Wallace R, Orciari L, Hightower D, Yager P, et al. (2013) Rabies surveillance in the United States during 2012. J Am Vet Med Assoc 243: 805–815.
- 15.
Kim BI, Blanton JD, Gilbert A, Castrodale L, Hueffer K, et al.. (2013) A Conceptual Model for the Impact of Climate Change on Fox Rabies in Alaska, 1980–2010. Zoonoses Public Health.
- 16. Gordon ER, Curns AT, Krebs JW, Rupprecht CE, Real LA, et al. (2004) Temporal dynamics of rabies in a wildlife host and the risk of cross-species transmission. Epidemiol Infect 132: 515–524.
- 17. Borucki MK, Chen-Harris H, Lao V, Vanier G, Wadford DA, et al. (2013) Ultra-Deep Sequencing of Intra-host Rabies Virus Populations during Cross-species Transmission. PLoS Negl Trop Dis 7: e2555.
- 18.
Bureau USC (2011) 2010 Census Summary File 1— United States.
- 19. Smith JS (1996) New aspects of rabies with emphasis on epidemiology, diagnosis, and prevention of the disease in the United States. Clin Microbiol Rev 9: 166–176.
- 20. Sidwa TJ, Wilson PJ, Moore GM, Oertli EH, Hicks BN, et al. (2005) Evaluation of oral rabies vaccination programs for control of rabies epizootics in coyotes and gray foxes: 1995–2003. J Am Vet Med Assoc 227: 785–792.
- 21. Clark KA, Neill SU, Smith JS, Wilson PJ, Whadford VW, et al. (1994) Epizootic canine rabies transmitted by coyotes in south Texas. J Am Vet Med Assoc 204: 536–540.
- 22. Kim HJ, Fay MP, Feuer EJ, Midthune DN (2000) Permutation tests for joinpoint regression with applications to cancer rates. Stat Med 19: 335–351.
- 23.
ESRI (2011) ArcGIS Desktop. Release 10. ed. Redlands, CA: Environmental Systems Research Institute.
- 24. Chen J, Roth RE, Naito AT, Lengerich EJ, Maceachren AM (2008) Geovisual analytics to enhance spatial scan statistic interpretation: an analysis of U.S. cervical cancer mortality. Int J Health Geogr 7: 57.
- 25. Centers for Disease Control and Prevention (2000) Update: raccoon rabies epizootic–United States and Canada, 1999. MMWR Morb Mortal Wkly Rep 49: 31–35.
- 26. Kappus KD, Bigler WJ, McLean RG, Trevino HA (1970) The raccoon an emerging rabies host. J Wildl Dis 6: 507–509.
- 27. Kuzmina NA, Lemey P, Kuzmin IV, Mayes BC, Ellison JA, et al. (2013) The phylogeography and spatiotemporal spread of South-central skunk rabies virus. PLoS One 8: e82348.
- 28. Anderson RM, Jackson HC, May RM, Smith AM (1981) Population dynamics of fox rabies in Europe. Nature 289: 765–771.
- 29. Smith HT, Engeman RM (2002) An extraordinary Raccoon, Procyon lotor, density at an Urban Park. Canadian Field-Naturalist 116: 636–639.
- 30.
Rosatte RC, Power MJ, Macinnes CD (1992) Density, Dispersion, Movements and Habitat of Skunks (Mephitis-Mephitis) and Raccoons (Procyon-Lotor) in Metropolitan Toronto. Wildlife 2001: Populations: 932–944.
- 31. Riley SPD, Hadidian J, Manski DA (1998) Population density, survival, and rabies in raccoons in an urban national park. Canadian Journal of Zoology-Revue Canadienne De Zoologie 76: 1153–1164.
- 32. Prange S, Gehrt SD, Wiggers EP (2003) Demographic factors contributing to high raccoon densities in urban landscapes. Journal of Wildlife Management 67: 324–333.
- 33. Broadfoot JD, Rosatte RC, O'Leary DT (2001) Raccoon and skunk population models for urban disease control planning in Ontario, Canada. Ecological Applications 11: 295–303.
- 34.
Fascione N, Delach A, Smith ME, Defenders of Wildlife. (2004) People and predators: from conflict to coexistence. Washington: Island Press. xvi: , 285 p. p.
- 35. Hirsch BT, Prange S, Hauver SA, Gehrt SD (2013) Raccoon social networks and the potential for disease transmission. PLoS One 8: e75830.
- 36.
Feldhamer GA, Thompson BC, Chapman JA (2003) Wild mammals of North America: biology, management, and conservation. Baltimore, Md.: Johns Hopkins University Press. xiii: , 1216 p. p.
- 37. Weissinger MD, Theimer TC, Bergman DL, Deliberto TJ (2009) Nightly and seasonal movements, seasonal home range, and focal location photo-monitoring of urban striped skunks (Mephitis mephitis): implications for rabies transmission. J Wildl Dis 45: 388–397.
- 38.
Bowen-Davies J, Lowings P (2000) Current perspectives on rabies 1. The biology of rabies and rabies-related viruses. In Practice 22: : 118-+.
- 39. Rivest P, Bergeron JM (1981) Density, Food-Habits, and Economic Importance of Raccoons (Procyon-Lotor) in Quebec Agrosystems. Canadian Journal of Zoology-Revue Canadienne De Zoologie 59: 1755–1762.
- 40. Page LK, Swihart RK, Kazacos KR (1999) Implications of raccoon latrines in the epizootiology of baylisascariasis. J Wildl Dis 35: 474–480.
- 41. Totton SC, Tinline RR, Rosatte RC, Bigler LL (2002) Contact rates of raccoons (Procyon lotor) at a communal feeding site in rural eastern Ontario. J Wildl Dis 38: 313–319.
- 42.
Breed MD, Moore J (2012) Animal behavior. Amsterdam; Boston: Academic Press. xiii: , 475 p. p.
- 43. Shirer HW FH (1970) Comparison from radiotracking of movements and denning habits of the raccoon, striped skunk, and opossum in northeastern Kansas. Journal of Mammalogy 51: 13.
- 44. Hill RE Jr, Beran GW (1992) Experimental inoculation of raccoons (Procyon lotor) with rabies virus of skunk origin. J Wildl Dis 28: 51–56.
- 45. Hill RE, Smith KE, Beran GW, Beard PD (1993) Further-Studies on the Susceptibility of Raccoons (Procyon-Lotor) to a Rabies Virus of Skunk Origin and Comparative Susceptibility of Striped Skunks (Mephitis-Mephitis). Journal of Wildlife Diseases 29: 475–477.
- 46. Winkler WG, Schneider NJ, Jennings WL (1972) Experimental rabies infection in wild rodents. J Wildl Dis 8: 99–103.
- 47. Slate D, Rupprecht CE, Donovan D, Badcock J, Messier A, et al. (2008) Attaining raccoon rabies management goals: history and challenges. Dev Biol (Basel) 131: 439–447.
- 48. Houle M FD, Mainguy J, Canac-Marquis P (2011) Landscape composition and structure influence the abundance of mesopredators: implications for the control of the raccoon (Procyon lotor) variant of rabies. Can J Zool 89: 1107–1116.
- 49. Henner CM, Chamberlain MJ, Leopold BD, Burger LW (2004) A multi-resolution assessment of raccoon den selection. Journal of Wildlife Management 68: 179–187.