PLoS OnePLoS ONEplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, USA252957504189788PONE-D-14-1054410.1371/journal.pone.0107539Research ArticleBiology and Life SciencesEcologyVeterinary ScienceEcology and Environmental SciencesHabitatsNature-Society InteractionsMedicine and Health SciencesEpidemiologyInfectious DiseasesTropical DiseasesRight Place, Wrong Species: A 20-Year Review of Rabies Virus Cross Species Transmission among Terrestrial Mammals in the United StatesRight Place, Wrong Species: Cross Species Transmission of Rabies VirusWallaceRyan M.1*GilbertAmy2SlateDennis3ChipmanRichard3SinghAmber1Cassie Wedd1BlantonJesse D.1Division of High-Consequence Pathogens and Pathology, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of AmericaNational Wildlife Research Center, United States Department of Agriculture, Fort Collins, Colorado, United States of AmericaWildlife Services, United States Department of Agriculture, Concord, New Hampshire, United States of AmericaMarkotterWandaEditorUniversity of Pretoria, South Africa* E-mail: euk5@cdc.gov

Competing Interests: The authors have declared that no competing interests exist.

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.

20148102014910e107539283201420820142014This 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.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.

These authors have no support or funding to report.Data AvailabilityThe 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.
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.

MethodsData 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].

ResultsAnnual Trends in Rabies SpilloverRaccoon 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).

10.1371/journal.pone.0107539.g001Annual Rates of Raccoon Rabies Virus Variant Cross Species Transmission, 1990–2011.10.1371/journal.pone.0107539.t001Trends in rabies cross species transmission (CST) by species or species grouping, 1990–2011.
Cross Species TransmissionRabid Animals 1990–20111990 CST Rate2011 CST RateOverall CST ChangeLatest Trend Time PeriodModel APC*95% CIP-Value
RaccoonCST by AnimalSkunk16,6000.160.35118.0%1999–2011−0.9%−3.2, 1.40.4
Fox6,5370.050.18251.0%1996–20114.0%3.0, 5.0<0.01
Domestic5,0280.060.14139.7%1990–20115.2%4.5, 5.9<0.01
Cattle6440.010.0277.8%1990–20114.0%2.5, 5.7<0.01
Other2,0670.030.0457.1%1990–20112.7%1.7, 3.8<0.01
Total30,8760.310.73138.2%199820111.2%−0.1, 2.60.06
CST by RegionNortheast12,5420.130.62376.2%1999–2011−3.9%−7.2, −0.60.03
Mid-Atlantic11,1290.370.79110.5%1990–20112.7%2.0, 3.5<0.01
Southeast7,2050.390.7694.9%1994–20115.8%5.0, 6.6<0.01
SkunkCST by AnimalRaccoon1480.010.02183.3%1990–20116.3%2.4, 10.3<0.01
Fox3090.020.0355.6%1998–20119.5%3.5, 15.8<0.01
Domestic1,6010.100.07−35.9%1990–2011−1.0%−2.5, 0.50.19
Cattle1,1400.110.02−83.0%1990–2011−5.7%−7.8, −3.4<0.01
Other5800.050.04−13.3%1990–20110.50%−2.3, 3.30.73
Total3,7780.280.17−40.3%19902011−2.0%−3.6, −0.40.02
CST by RegionMidwest1,3050.480.20−59.0%1990–2011−1.5%−3.6, 0.70.16
South1,2120.160.178.8%2001–20117.0%3.0, 11.1<0.01
Southeast1680.280.18−36.4%2001–201113.6%1.9, 26.7<0.01
West1,0930.220.11−48.8%1990–20111.1%−0.8, 3.10.25
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.

10.1371/journal.pone.0107539.g002Annual Rates of Skunk Rabies Virus Variant Cross Species Transmission, 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).

10.1371/journal.pone.0107539.t002Animal Submission for Rabies Diagnosis and Cross Species Transmission of Rabies Virus by Non-Reservoir Animal, 2007–2011.
Raccoon Variant RegionSkunk Variant Region
Testing SubmissionsNumber Rabid (Positive Test Rate)CST RateTesting SubmissionsNumber Rabid (Positive Test Rate)CST RateRaccoon: Skunk Rate Ratio
Raccoon53,87711,134 (20.7%)-4,13278 (1.9%)0.0218.7 (14.2, 24.5)
Skunk10,6053,379 (31.9%)0.307,4853,329 (44.5%)-
Fox5,4041,720 (31.8%)0.1569866 (9.5%)0.028.1 (6.3, 10.5)
Domestic Animals100,1631,254 (1.3%)0.1136,234251 (0.7%)0.081.5 (1.3, 1.7)
Cattle1,923154 (8.0%)0.011,05593 (8.8%)0.030.5 (0.4, 0.6)
Other Animals19,422436 (2.2%)0.043,411101 (3.0%)0.031.3 (1.2, 1.4)
Total191,39418,077 (9.4%)0.6253,0153,918 (7.4%)0.173.6 (3.4, 3.9)
Spatial AnalysisRaccoon 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).

10.1371/journal.pone.0107539.g003Five-year aggregate spillover rates of raccoon rabies virus variant, by County: 1992–2011.
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).

10.1371/journal.pone.0107539.g004Five-year aggregate spillover rates of skunk rabies virus variants by County: 1992–2011.

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%.

10.1371/journal.pone.0107539.g005Clusters of High Rates of Raccoon Variant Cross Species Transmission over Varying Spatial Cluster Values, 2007–2011.
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.

10.1371/journal.pone.0107539.g006Comparison of Confirmed Rabies Virus Variant and Predicted Rabies Virus Variant, 2007–2011.
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).

10.1371/journal.pone.0107539.t003Linear Regression of Reservoir Submission Rate and Cross Species Transmission Rate: Skunk and Raccoon Variant Regions, 2007–2011.
VariableDFParameter EstimateStandard Errort ValuePr> |t|
Raccoon Variant RegionRaccoon Intercept182.563.0460827.10<.0001
Raccoon Submission Rate1−0.0180.00707−2.520.0118
Skunk Variant RegionsSkunk Intercept127.832.2804612.20<.0001
Skunk Submission Rate1−0.147280.05520−2.670.0080
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.

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.

ReferencesKnobelDL, CleavelandS, ColemanPG, FevreEM, MeltzerMI, et al (2005) Re-evaluating the burden of rabies in Africa and Asia. Bull World Health Organ83: 36036815976877NettlesVF, ShaddockJH, SikesRK, ReyesCR (1979) Rabies in translocated raccoons. Am J Public Health69: 601602443502Velasco-VillaA, ReederSA, OrciariLA, YagerPA, FrankaR, et al (2008) Enzootic rabies elimination from dogs and reemergence in wild terrestrial carnivores, United States. Emerg Infect Dis14: 1849185419046506BlantonJD, DyerJ, McBrayerJ, RupprechtCE (2012) Rabies surveillance in the United States during 2011. J Am Vet Med Assoc241: 71272222947154Freer L (2004) North American wild mammalian injuries. Emerg Med Clin North Am 22: : 445–473, ix.KizerKW (1979) Epidemiologic and clinical aspects of animal bite injuries. JACEP8: 134141430939BlantonJD, PalmerD, DyerJ, RupprechtCE (2011) Rabies Surveillance in the United States during 2010. JAVMA239: 1121718191Constantine 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.StreickerDG, TurmelleAS, VonhofMJ, KuzminIV, McCrackenGF, et al (2010) Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science329: 67667920689015McQuistonJH, YagerPA, SmithJS, RupprechtCE (2001) Epidemiologic characteristics of rabies virus variants in dogs and cats in the United States, 1999. J Am Vet Med Assoc218: 1939194211417737KuzminIV, ShiM, OrciariLA, YagerPA, Velasco-VillaA, et al (2012) Molecular Inferences Suggest Multiple Host Shifts of Rabies Viruses from Bats to Mesocarnivores in Arizona during 2001–2009. PLoS Pathog8: e100278622737076GuerraMA, CurnsAT, RupprechtCE, HanlonCA, KrebsJW, et al (2003) Skunk and raccoon rabies in the eastern United States: temporal and spatial analysis. Emerg Infect Dis9: 1143115014519253OertliEH, WilsonPJ, HuntPR, SidwaTJ, RohdeRE (2009) Epidemiology of rabies in skunks in Texas. J Am Vet Med Assoc234: 61662019250039DyerJL, WallaceR, OrciariL, HightowerD, YagerP, et al (2013) Rabies surveillance in the United States during 2012. J Am Vet Med Assoc243: 80581524004227Kim 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.GordonER, CurnsAT, KrebsJW, RupprechtCE, RealLA, et al (2004) Temporal dynamics of rabies in a wildlife host and the risk of cross-species transmission. Epidemiol Infect132: 51552415188720BoruckiMK, Chen-HarrisH, LaoV, VanierG, WadfordDA, et al (2013) Ultra-Deep Sequencing of Intra-host Rabies Virus Populations during Cross-species Transmission. PLoS Negl Trop Dis7: e255524278493Bureau USC (2011) 2010 Census Summary File 1— United States.SmithJS (1996) New aspects of rabies with emphasis on epidemiology, diagnosis, and prevention of the disease in the United States. Clin Microbiol Rev9: 1661768964034SidwaTJ, WilsonPJ, MooreGM, OertliEH, HicksBN, 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 Assoc227: 78579216178403ClarkKA, NeillSU, SmithJS, WilsonPJ, WhadfordVW, et al (1994) Epizootic canine rabies transmitted by coyotes in south Texas. J Am Vet Med Assoc204: 5365408163414KimHJ, FayMP, FeuerEJ, MidthuneDN (2000) Permutation tests for joinpoint regression with applications to cancer rates. Stat Med19: 33535110649300ESRI (2011) ArcGIS Desktop. Release 10. ed. Redlands, CA: Environmental Systems Research Institute.ChenJ, RothRE, NaitoAT, LengerichEJ, MaceachrenAM (2008) Geovisual analytics to enhance spatial scan statistic interpretation: an analysis of U.S. cervical cancer mortality. Int J Health Geogr7: 5718992163Centers for Disease Control and Prevention (2000) Update: raccoon rabies epizootic–United States and Canada, 1999. MMWR Morb Mortal Wkly Rep49: 313510680602KappusKD, BiglerWJ, McLeanRG, TrevinoHA (1970) The raccoon an emerging rabies host. J Wildl Dis6: 50750916512166KuzminaNA, LemeyP, KuzminIV, MayesBC, EllisonJA, et al (2013) The phylogeography and spatiotemporal spread of South-central skunk rabies virus. PLoS One8: e8234824312657AndersonRM, JacksonHC, MayRM, SmithAM (1981) Population dynamics of fox rabies in Europe. Nature289: 7657717464941SmithHT, EngemanRM (2002) An extraordinary Raccoon, Procyon lotor, density at an Urban Park. Canadian Field-Naturalist116: 636639Rosatte 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.RileySPD, HadidianJ, ManskiDA (1998) Population density, survival, and rabies in raccoons in an urban national park. Canadian Journal of Zoology-Revue Canadienne De Zoologie76: 11531164PrangeS, GehrtSD, WiggersEP (2003) Demographic factors contributing to high raccoon densities in urban landscapes. Journal of Wildlife Management67: 324333BroadfootJD, RosatteRC, O'LearyDT (2001) Raccoon and skunk population models for urban disease control planning in Ontario, Canada. Ecological Applications11: 295303Fascione N, Delach A, Smith ME, Defenders of Wildlife. (2004) People and predators: from conflict to coexistence. Washington: Island Press. xvi: , 285 p. p.HirschBT, PrangeS, HauverSA, GehrtSD (2013) Raccoon social networks and the potential for disease transmission. PLoS One8: e7583024130746Feldhamer 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.WeissingerMD, TheimerTC, BergmanDL, DelibertoTJ (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 Dis45: 38839719395748Bowen-Davies J, Lowings P (2000) Current perspectives on rabies 1. The biology of rabies and rabies-related viruses. In Practice 22: : 118-+.RivestP, BergeronJM (1981) Density, Food-Habits, and Economic Importance of Raccoons (Procyon-Lotor) in Quebec Agrosystems. Canadian Journal of Zoology-Revue Canadienne De Zoologie59: 17551762PageLK, SwihartRK, KazacosKR (1999) Implications of raccoon latrines in the epizootiology of baylisascariasis. J Wildl Dis35: 47448010479081TottonSC, TinlineRR, RosatteRC, BiglerLL (2002) Contact rates of raccoons (Procyon lotor) at a communal feeding site in rural eastern Ontario. J Wildl Dis38: 31331912038131Breed MD, Moore J (2012) Animal behavior. Amsterdam; Boston: Academic Press. xiii: , 475 p. p.Shirer HWFH (1970) Comparison from radiotracking of movements and denning habits of the raccoon, striped skunk, and opossum in northeastern Kansas. Journal of Mammalogy51: 13HillREJr, BeranGW (1992) Experimental inoculation of raccoons (Procyon lotor) with rabies virus of skunk origin. J Wildl Dis28: 51561548802HillRE, SmithKE, BeranGW, BeardPD (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 Diseases29: 4754778355351WinklerWG, SchneiderNJ, JenningsWL (1972) Experimental rabies infection in wild rodents. J Wildl Dis8: 991035007906SlateD, RupprechtCE, DonovanD, BadcockJ, MessierA, et al (2008) Attaining raccoon rabies management goals: history and challenges. Dev Biol (Basel)131: 43944718634506Houle MFD, MainguyJ, Canac-MarquisP (2011) Landscape composition and structure influence the abundance of mesopredators: implications for the control of the raccoon (Procyon lotor) variant of rabies. Can J Zool89: 11071116HennerCM, ChamberlainMJ, LeopoldBD, BurgerLW (2004) A multi-resolution assessment of raccoon den selection. Journal of Wildlife Management68: 179187