Data-Driven Management-A Dynamic Occupancy Approach to Enhanced Rabies Surveillance Prioritization

Evaluation of contingency actions to control the spread of raccoon rabies in Ohio and Virginia

The raccoon (Procyon lotor) variant of the rabies virus (RRV) is enzootic in the eastern United States and oral rabies vaccination (ORV) is the primary strategy to prevent and control landscape spread. Breaches of ORV management zones occasionally occur, and emergency "contingency" actions may be implemented to enhance local control. Contingency actions are an integral part of landscape-scale wildlife rabies management but can be very costly and routinely involve enhanced rabies surveillance (ERS) around the index case. We investigated two contingency actions in Ohio (2017-2019 and 2018-2021) and one in Virginia (2017-2019) using a dynamic, multi-method occupancy approach to examine relationships between specific management actions and RRV occurrence, including whether ERS was sufficient around the index case. The RRV occupancy was assessed seasonally at 100-km2 grids and we examined relationships across three spatial scales (regional management zone, RRV free regions, and local contingency areas). The location of a grid relative to the ORV management zone was the strongest predictor of RRV occupancy at the regional scale. In RRV free regions, the neighbor effect and temporal variability were most important in influencing RRV occupancy. Parenteral (hand) vaccination of raccoons was important across all three contingency action areas, but more influential in the Ohio contingency action areas where more raccoons were hand vaccinated. In the Virginia contingency action area, ORV strategies were as important in reducing RRV occupancy as a hand vaccination strategy. The management action to trap, euthanize, and test (TET) raccoons was an important method to increase ERS, yet the impacts of TET on RRV occupancy are not clear. The probability of detecting additional cases of RRV was exceptionally high (>0.95) during the season the index case occurred. The probability of detecting RRV through ERS declined in the seasons following initial TET efforts but remained higher after the contingency action compared to the ERS detection probabilities prior to index case incidence. Local RRV cases were contained within one year and eliminated within 2-3 years of each contingency action.

Assessing the Efficiency of Local Rabies Vaccination Strategies for Raccoons (Procyon lotor) in an Urban Setting

Raccoon rabies virus (RRV) has been managed using multiple vaccination strategies, including oral rabies vaccination and trap-vaccinate-release (TVR). Identifying a rabies vaccination strategy for an area is a nontrivial task. Vaccination strategies differ in the amount of effort and monetary costs required to achieve a particular level of vaccine seroprevalence (efficiency). Simulating host movement relative to different vaccination strategies in silico can provide a useful tool for exploring the efficiency of different vaccination strategies. We refined a previously developed individual-based model of raccoon movement to evaluate vaccination strategies for urban Hamilton, Ontario, Canada. We combined different oral rabies vaccination baiting (hand baiting, helicopter, and bait stations) with TVR strategies and used GPS data to parameterize and simulate raccoon movement in Hamilton. We developed a total of 560 vaccination strategies, in consultation with the Ontario Ministry of Natural Resources and Forestry, for RRV control in Hamilton. We documented the monetary costs of each vaccination strategy and estimated the population seroprevalence. Intervention costs and seroprevalence estimates were used to calculate the efficiency of each strategy to meet targets set for the purpose of RRV control. Estimated seroprevalence across different strategies varied widely, ranging from less than 5% to more than 70%. Increasing bait densities (distributed using by hand or helicopter) led to negligible increase in seroprevalence. Helicopter baiting was the most efficient and TVR was the least efficient, but helicopter-based strategies led to lower levels of seroprevalence (6-12%) than did TVR-based strategies (17-70%). Our simulations indicated that a mixed strategy including at least some TVR may be the most efficient strategy for a local urban RRV control program when seroprevalence levels >30% may be required. Our simulations provide information regarding the efficiency of different vaccination strategies for raccoon populations, to guide local RRV control in urban settings.

Are rabid raccoons (Procyon lotor) ready for the rapture? Determining the geographic origin of rabies virus-infected raccoons using RADcapture and microhaplotypes

North America is recognized for the exceptional richness of rabies virus (RV) wildlife reservoir species. Management of RV is accomplished through vaccination targeting mesocarnivore reservoir populations, such as the raccoon (Procyon lotor) in Eastern North America. Raccoons are a common generalist species, and populations may reach high densities in developed areas, which can result in contact with humans and pets with potential exposures to the raccoon variant of RV throughout the eastern United States. Understanding the spatial movement of RV by raccoon populations is important for monitoring and refining strategies supporting the landscape-level control and local elimination of this lethal zoonosis. We developed a high-throughput genotyping panel for raccoons based on hundreds of microhaplotypes to identify population structure and genetic diversity relevant to rabies management programs. Throughout the eastern United States, we identified hierarchical population genetic structure with clusters that were connected through isolation-by-distance. We also illustrate that this genotyping approach can be used to support real-time management priorities by identifying the geographic origin of a rabid raccoon that was collected in an area of the United States that had been raccoon RV-free for 8 years. The results from this study and the utility of the microhaplotype panel and genotyping method will provide managers with information on raccoon ecology that can be incorporated into future management decisions.

Biodiversity indices and Random Forests reveal the potential for striped skunk (Mephitis mephitis) fecal microbial communities to function as a biomarker for oral rabies vaccination

Wildlife disease surveillance and monitoring poses unique challenges when assessing rates of population vaccination, immunity, or infection prevalence. Non-invasively detected biomarkers can help reduce risk to both animal and field personnel during wildlife disease management activities. In this study, we investigated the utility of fecal microbiome data collected from captive striped skunks (Mephitis mephitis) in predicting rabies virus vaccination and infection status. We sequenced the hypervariable region 4 (V4) of the bacterial 16S gene and estimated alpha and beta diversity across timepoints in three groups of skunks: vaccination then rabies virus infection, sham vaccination then rabies virus infection, and rabies virus infected without vaccination. Alpha diversity did not differ among treatment groups but beta diversity between treatments was statistically significant. The phyla Firmicutes and Proteobacteria were dominant among all samples. Using Random Forests, we identified operational taxonomic units (OTUs) that greatly influenced classification of fecal samples into treatment groups. Each of these OTUs was correlated with fecal volatile organic compounds detected from the samples for companion treatment groups in another study. This research is the first to highlight striped skunk microbiome biodiversity as a vaccination biomarker which pushes the frontier on alternative methods for surveillance and monitoring of vaccination and disease in wildlife populations.

Spatial ecology of translocated raccoons

Raccoons (Procyon lotor) are routinely translocated both legally and illegally to mitigate conflicts with humans, which has contributed to the spread of rabies virus across eastern North America. The movement behavior of translocated raccoons has important ramifications for disease transmission yet remains understudied and poorly quantified. To examine the spatial ecology of raccoons following experimental translocation, we performed reciprocal 16 km-distance translocations of 30 raccoons between habitats of high and low raccoon density (bottomland hardwood and upland pine, respectively) across the Savannah River Site (SRS) in Aiken, South Carolina, USA (2018-2019). Translocation influenced patterns of raccoon space use, with translocated animals exhibiting a 13-fold increase in 95% utilization distributions (UDs) post- compared to pre-translocation (mean 95% UD 35.8 ± 36.1 km² vs 1.96 ± 1.17 km²). Raccoons originating from upland pine habitats consistently had greater space use and larger nightly movement distances post-translocation compared to raccoons moved from bottomland hardwood habitats, whereas these differences were generally not observed prior to translocation. Estimated home ranges of male raccoons were twice the area as estimated for female raccoons, on average, and this pattern was not affected by translocation. After a transient period lasting on average 36.5 days (SD = 30.0, range = 3.25-92.8), raccoons often resumed pre-experiment movement behavior, with 95% UD sizes not different from those prior to translocation (mean = 2.27 ± 1.63km²). Most animals established new home ranges after translocation, whereas three raccoons moved > 16 km from their release point back to the original capture location. Four animals crossed a 100-m wide river within the SRS post-translocation, but this behavior was not documented among collared raccoons prior to translocation. Large increases in space use combined with the crossing of geographic barriers such as rivers may lead to elevated contact rates with conspecifics, which can heighten disease transmission risks following translocation. These results provide additional insights regarding the potential impacts of raccoon translocation towards population level risks of rabies outbreaks and underscore the need to discourage mesocarnivore translocations to prevent further spread of wildlife rabies.

Spatial and temporal overlap of domestic cats (Felis catus) and native urban wildlife

Free-roaming domestic cats (Felis catus) are known to pose threats to ecosystem health via transmission of zoonotic diseases and predation of native wildlife. Likewise, free-roaming cats are also susceptible to predation or disease transmission from native wildlife. Physical interactions are required for many of these risks to be manifested, necessitating spatial and temporal overlap between cats and wildlife species. Therefore, knowledge of the location and extent of shared habitat and activity periods would benefit management programs. We used data from a 3-year camera trap survey to model species-specific occupancy and identify landscape variables that contribute to the distribution of free-roaming domestic cats and eight native mammal species in Washington, DC. (USA). Our analysis includes five species that are common prey items of domestic cats, and three species that are potential disease vectors or are otherwise known to be a risk to cats. We then predicted the probability of occupancy and estimated the probability of spatial overlap between cats and each native wildlife species at multiple scales. We also used kernel density estimations to calculate temporal overlap between cats and each native wildlife species. Across spatial scales, occupancy for potential disease vector species was generally positively correlated with canopy cover and open water. Prey species were also generally positively correlated with canopy cover, but displayed negative associations with human population density and inconsistent associations with average per capita income. Domestic cat occupancy was negatively correlated with natural habitat characteristics and positively correlated with human population density. Predicted spatial overlap between domestic cats and native wildlife was greatest for potential disease vector species. Temporal overlap was high (>0.50) between cats and all but two native wildlife species, indicating that temporal overlap is probable wherever species overlap spatially. Our findings indicate that the risk to and from domestic cats varies across urban landscapes, but primarily arises from human activities. As such, humans are implicated in the negative outcomes that result from cats interacting with wildlife. Data-driven management to reduce such interactions can aid in cat population management, biodiversity conservation, and public health campaigns.