White-tailed Deer Populations Research Articles

Modeling platform to assess the effectiveness of single and integrated Ixodes scapularis tick control methods

BackgroundLyme disease continues to expand in Canada and the USA and no single intervention is likely to curb the epidemic.MethodsWe propose a platform to quantitatively assess the effectiveness of a subset of Ixodes scapularis tick management approaches. The platform allows us to assess the impact of different control treatments, conducted either individually (single interventions) or in combination (combined efforts), with varying timings and durations. Interventions include three low environmental toxicity measures in differing combinations, namely reductions in white-tailed deer (Odocoileus virginianus) populations, broadcast area-application of the entomopathogenic fungus Metarhizium anisopliae, and fipronil-based rodent-targeted bait boxes. To assess the impact of these control efforts, we calibrated a process-based mathematical model to data collected from residential properties in the town of Redding, southwestern Connecticut, where an integrated tick management program to reduce I.xodes scapularis nymphs was conducted from 2013 through 2016. We estimated parameters mechanistically for each of the three treatments, simulated multiple combinations and timings of interventions, and computed the resulting percent reduction of the nymphal peak and of the area under the phenology curve.ResultsSimulation outputs suggest that the three-treatment combination and the bait boxes–deer reduction combination had the overall highest impacts on suppressing I. scapularis nymphs. All (single or combined) interventions were more efficacious when implemented for a higher number of years. When implemented for at least 4 years, most interventions (except the single application of the entomopathogenic fungus) were predicted to strongly reduce the nymphal peak compared with the no intervention scenario. Finally, we determined the optimal period to apply the entomopathogenic fungus in residential yards, depending on the number of applications.ConclusionsComputer simulation is a powerful tool to identify the optimal deployment of individual and combined tick management approaches, which can synergistically contribute to short-to-long-term, costeffective, and sustainable control of tick-borne diseases in integrated tick management (ITM) interventions.Graphical

Maternal behaviors influence survival of ungulate neonates under heavy predation risk.

Predators impose top-down forces on prey populations, with the strength of those effects often varying over space and time and among demographic groups. In ungulates, predation risk is typically greatest for neonatal offspring, with some suggesting that predators can key in on adult activity to locate hidden neonates. However, few field studies to date have been able to directly assess the influence of maternal care on ungulate neonate survival. Using a population of white-tailed deer under heavy coyote predation pressure, we tested the maternal dispersion hypothesis, which suggests the dispersion of maternal activity temporally and spatially attenuates risk of predation for ungulate neonates during this vulnerable altricial phase. We compared support for this hypothesis with more commonly tested hypotheses regarding the influence of habitat conditions and intrinsic factors on neonatal survival. Fawn survival to 16 weeks was 27.7%, with coyotes accounting for 59% of fawn mortalities. In support of our maternal temporal diffusion hypothesis, we found that neonatal survival decreased as more maternal visits (proportionally) occurred at night. The only other significant (p < .1) predictor of fawn survival was birth timing, with fawn survival decreasing as the season progressed. Given that fawn survival declined as the proportion of nighttime visits increased, and that wild pig presence and human disturbance can push doe and fawn activity toward nocturnal hours, additional research is needed to determine whether managing pig and human disturbance can decrease fawn mortality. More broadly, given the importance of recruitment in ungulate population dynamics, our finding opens a potentially important new line of inquiry on how maternal behaviors influence predation risk in large animal predator-prey ecology.

Using UAVs Equipped With Thermal Cameras To Survey A Known Population Of White-tailed Deer

Using UAVs Equipped With Thermal Cameras To Survey A Known Population Of White-tailed Deer

Surveillance and detection of Haemaphysalis longicornis (Acari: Ixodidae) in protected areas from Georgia, USA.

The invasion of the Asian longhorned tick, Haemaphysalis longicornis Neumann, into the United States poses a significant ecological, veterinarian, and public health threat. We conducted a comprehensive tick survey using standard tick flagging protocol for collection over 3 field seasons, March-August, and 56 surveyed sites to identify the presence of H. longicornis in the native ecosystem and prove its establishment according to USDA criteria. Of the total 56 state parks and wildlife management areas (WMA) surveyed, only one was found to be invaded by H. longicornis; detection of H. longicornis occurred at Buck Shoals Wildlife Management area in White County, GA. This site is maintained by the state of Georgia, has no agricultural animals present, and hosts a large white-tailed deer population. After the initial detection of H. longicornis in 2022, an additional field season occurred in 2023, where H. longicornis was confirmed as established based on USDA criteria. The increase in H. longicornis populations from 2021 to 2023 at Buck Shoals WMA points to the rapid spread of this tick within the environment. Our findings provide evidence of the rapid establishment of H. longicornis in the southern edge of suitability for this tick and within the native ecosystem beyond farmlands and private land.

Evaluating the Interaction of Emerging Diseases on White-Tailed Deer Populations Using an Agent-Based Modeling Approach.

Disease co-occurrence in wildlife populations is common yet understudied. In the case of disease-caused mortality, the mortality attributed to one disease has the potential to buffer populations against subsequent alternative disease outbreaks by reducing populations and thus contacts needed to sustain disease transmission. However, substantial disease-driven population declines may also prevent populations from recovering, leading to localized extinctions. Hemorrhagic disease (HD), a vector-transmitted, viral disease in white-tailed deer (WTD), similar to chronic wasting disease (CWD), a prion disease, has increased in frequency and distribution in the United States. However, unlike CWD, which progresses slowly, HD can cause mortality only days after infection. Hemorrhagic disease outbreaks can result in substantial localized mortality events in WTD near vector habitats such as wetlands and may reduce local deer densities and consequent CWD transmission. The objective of our study was to evaluate the potential for HD outbreaks to buffer CWD risk where the diseases co-occur. Using an agent-based modeling approach, we found that frequent, intense HD outbreaks have the potential to mitigate CWD risk, especially if those outbreaks occur shortly after CWD introduction. However, HD outbreaks that do not result in substantial WTD mortality are unlikely to impact CWD or WTD population dynamics. Severe HD outbreaks may reduce CWD cases and could present an opportunity for managers to boost CWD control initiatives in a post-HD outbreak year.

White-tailed deer population dynamics in a multipredator landscape shaped by humans.

Large terrestrial mammals increasingly rely on human-modified landscapes as anthropogenic footprints expand. Land management activities such as timber harvest, agriculture, and roads can influence prey population dynamics by altering forage resources and predation risk via changes in habitat, but these effects are not well understood in regions with diverse and changing predator guilds. In northeastern Washington state, USA, white-tailed deer (Odocoileus virginianus) are vulnerable to multiple carnivores, including recently returned gray wolves (Canis lupus), within a highly human-modified landscape. To understand the factors governing predator-prey dynamics in a human context, we radio-collared 280 white-tailed deer, 33 bobcats (Lynx rufus), 50 cougars (Puma concolor), 28 coyotes (C. latrans), and 14 wolves between 2016 and 2021. We first estimated deer vital rates and used a stage-structured matrix model to estimate their population growth rate. During the study, we observed a stable to declining deer population (lambda = 0.97, 95% confidence interval: 0.88, 1.05), with 74% of Monte Carlo simulations indicating population decrease and 26% of simulations indicating population increase. We then fit Cox proportional hazard models to evaluate how predator exposure, use of human-modified landscapes, and winter severity influenced deer survival and used these relationships to evaluate impacts on overall population growth. We found that the population growth rate was dually influenced by a negative direct effect of apex predators and a positive effect of timber harvest and agricultural areas. Cougars had a stronger effect on deer population dynamics than wolves, and mesopredators had little influence on the deer population growth rate. Areas of recent timber harvest had 55% more forage biomass than older forests, but horizontal visibility did not differ, suggesting that timber harvest did not influence predation risk. Although proximity to roads did not affect the overall population growth rate, vehicle collisions caused a substantial proportion of deer mortalities, and reducing these collisions could be a win-win for deer and humans. The influence of apex predators and forage indicates a dual limitation by top-down and bottom-up factors in this highly human-modified system, suggesting that a reduction in apex predators would intensify density-dependent regulation of the deer population owing to limited forage availability.

Demographic Risk Factors Vary in the Invasion Front of Chronic Wasting Disease in West Virginia, USA.

After detecting chronic wasting disease (CWD) in white-tailed deer (Odocoileus virginianus) in Hampshire County, West Virginia, USA, in 2005, we investigated the change of CWD apparent prevalence and potential factors influencing infection risk during the invasion front. Over eight sampling years (2006-2012 and 2017) during a 12-yr period within a 101-km2-area monitoring zone, we sampled and tested a total of 853 deer for CWD by ELISA and immunohistochemistry. Bayesian logistic regression of risk factors included collection year, age class, sex, and adjusted body weight (weight after accounting for sex, age, kidney fat index, and number of fetuses). In the whole-herd model (n=634), collection year, age, and adjusted body weight were associated with increased odds of CWD, whereas an age-weight interaction had a negative relationship. We found that males drove the positive associations with age and adjusted body weight, whereas females were responsible for the negative interaction effect. These findings suggest potential behavioral and physiological mechanisms related to sex that may influence CWD exposure. Older males exhibited higher CWD prevalence, aligning with previous studies. Notably, the novel finding of adjusted body weight as a risk factor in males warrants further investigation, and this study highlights the need for future research on social behavior and its role in CWD transmission within white-tailed deer populations.

Island demographics and trait associations in white-tailed deer.

When a population is isolated and composed of few individuals, genetic drift is the paramount evolutionary force and results in the loss of genetic diversity. Inbreeding might also occur, resulting in genomic regions that are identical by descent, manifesting as runs of homozygosity (ROHs) and the expression of recessive traits. Likewise, the genes underlying traits of interest can be revealed by comparing fixed SNPs and divergent haplotypes between affected and unaffected individuals. Populations of white-tailed deer (Odocoileus virginianus) on islands of Saint Pierre and Miquelon (SPM, France) have high incidences of leucism and malocclusions, both considered genetic defects; on the Florida Keys islands (USA) deer exhibit smaller body sizes, a polygenic trait. Here we aimed to reconstruct island demography and identify the genes associated with these traits in a pseudo case-control design. The two island populations showed reduced levels of genomic diversity and a build-up of deleterious mutations compared to mainland deer; there was also significant genome-wide divergence in Key deer. Key deer showed higher inbreeding levels, but not longer ROHs, consistent with long-term isolation. We identified multiple trait-related genes in ROHs including LAMTOR2 which has links to pigmentation changes, and NPVF which is linked to craniofacial abnormalities. Our mixed approach of linking ROHs, fixed SNPs and haplotypes matched a high number (~50) of a-priori body size candidate genes in Key deer. This suite of biomarkers and candidate genes should prove useful for population monitoring, noting all three phenotypes show patterns consistent with a complex trait and non-Mendelian inheritance.

Persistence of viral RNA in North American elk experimentally infected with an ancestral strain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)

White-tailed deer (Odocoileus virginianus) have emerged as a reservoir host for SARS-CoV-2 given their susceptibility to infection and demonstrated high rates of seroprevalence and infection across the United States. As SARS-CoV-2 circulates within free-ranging white-tailed deer populations, there is the risk of transmission to other wildlife species and even back to the human population. The goal of this study was to determine the susceptibility, shedding, and immune response of North American elk (Cervus elaphus canadensis) to experimental infection with SARS-CoV-2, to determine if another wide-ranging cervid species could potentially serve as a reservoir host for the virus. Here we demonstrate that while North American elk do not develop clinical signs of disease, they do develop a neutralizing antibody response to infection, suggesting the virus is capable of replicating in this mammalian host. Additionally, we demonstrate SARS-CoV-2 RNA presence in the medial retropharyngeal lymph nodes of infected elk three weeks after experimental infection. Consistent with previous observations in humans, these data may highlight a mechanism of viral persistence for SARS-CoV-2 in elk.

Habitat alteration or climate: What drives the densities of an invading ungulate?

Anthropogenic habitat alteration and climate change are two well-known contributors to biodiversity loss through changes to species distribution and abundance; yet, disentangling the effects of these two factors is often hindered by their inherent confound across both space and time. We leveraged a contrast in habitat alteration associated with the jurisdictional boundary between two Canadian provinces to evaluate the relative effects of spatial variation in habitat alteration and climate on white-tailed deer (Odocoileus virginianus) densities. White-tailed deer are an invading ungulate across much of North America, whose expansion into Canada's boreal forest is implicated in the decline of boreal caribou (Rangifer tarandus caribou), a species listed as Threatened in Canada. We estimated white-tailed deer densities using 300 remote cameras across 12 replicated 50 km2 landscapes over 5 years. White-tailed deer densities were significantly lower in areas where winter severity was higher. For example, predicted deer densities declined from 1.83 to 0.35 deer/km2 when winter severity increased from the lowest value to the median value. There was a tendency for densities to increase with increasing habitat alteration; however, the magnitude of this effect was approximately half that of climate. Our findings suggest that climate is the primary driver of white-tailed deer populations; however, understanding the mechanisms underpinning this relationship requires further study of over-winter survival and fecundity. Long-term monitoring at the invasion front is needed to evaluate the drivers of abundance over time, particularly given the unpredictability of climate change and increasing prevalence of extreme weather events.

An individual-based model for direct and indirect transmission of chronic wasting disease in free-ranging white-tailed deer

Chronic wasting disease (CWD) is an infectious prion disease that infects members of the Cervidae family (i.e., deer) resulting in widespread ecological, economic, and recreational ramifications. We introduce a spatially explicit individual-based model (IBM) that integrates individual deer movement and behavior with population and disease dynamics to forecast CWD in populations of free-ranging white-tailed deer (Odocoileus virginianus). We use a Susceptible-Exposed-Infectious-Dead (S-E-I-D) epidemiological framework to explore spatiotemporal dynamics of CWD within an agriculturally dominated area in Michigan, USA. The IBM results closely mimicked documented short- and long-term dynamics of white-tailed deer populations and CWD in the Midwestern, USA. We applied pattern-oriented modeling using annual apparent CWD prevalence rates reported by Midwestern state wildlife agencies to validate the disease model. The introduction of a single infected deer to the modeled landscape (93 km2) led to an outbreak of CWD in 100 out of 350 model simulations (29 %); prevalence never exceeded 1.47 % for repetitions where the outbreak ended. For the 100 simulations where disease persisted, the deer population declined by 87 % by year 50 following initial introduction of CWD. Mean (±SD) prevalence after 5, 10, 25, and 50 years was 1.1 % (±1.0 %), 3.4 % (±3.3 %), 46.5 % (±18.8 %), and 51.8 % (±18.1 %), respectively, which highly correlated (r = 0.99) with annual CWD prevalence reported in Wisconsin white-tailed deer populations for years 1–21 post initial detection. Combined with a global sensitivity analysis, the IBM indicated that prevalence of CWD at year 20 was most sensitive to harvest rate of yearling and adult female deer and least sensitive to prion shedding rate, prion half-life, and deer group numbers, indicating that deer population parameters were more influential than disease parameters on CWD dynamics. Our IBM serves as a tool to explore and better understand indirect and direct transmission of CWD within free-ranging cervid populations. Users of this model can adjust parameter values to explore how interactions among individual deer and between deer and their environment affect CWD dynamics. This IBM also serves as a framework for applying and assessing spatially and temporally explicit management scenarios.

Spatial Variation of Population Density and Age Structure of White-Tailed Deer (Odocoileus virginianus) in a Colombian Moor

Wildlife population density is affected by resource availability, predators and competitors, and anthropogenic stressors. Its analysis and the factors affecting it are important for conservation, use, or population control. Based on estimates of density and population parameters, together with niche modeling, previous studies considered that the white-tailed deer population at the Chingaza National Natural Park (NNP) is undergoing exponential growth due to the marginal influence of factors regulating its abundance, such as predators. The aim of this study was to evaluate the population density and age structure of the white-tailed deer in the Chingaza NNP, Colombia. To evaluate the population density and age structure of white-tailed deer in the Monterredondo and La Paila sectors of the Chingaza NNP, 10 transects were installed with circular plots in each. These were visited twice; in each visit, we collected all deer pellet groups, measuring each pellet to estimate its volume. Population density was calculated using the Eberhardt and Van Etten method and the semi-automated PELLET procedure using FSC and FAR techniques; the age structure was estimated by k-means and fuzzy clustering analyses. With the FAR technique, the estimated deer density was 2.09 ind/km2 in Monterredondo and 0.94 ind/km2 in La Paila. With PELLET, these values were 0.37 ± 0.4 ind/km2 for Monterredondo and 0.16 ± 0.2 ind/km2 for La Paila. The proportion of age classes (adult, juvenile, and fawn) was 1:0.29:0 for Monterredondo and 1:0.57:0.04 for La Paila. The population density of white-tailed deer in two sectors of the Chingaza NNP is not homogeneous, but shows temporal and spatial variations. Population density has decreased versus reports from previous studies in this same area, with the highest deer abundance recorded in 2004, mostly in the Monterredondo sector. The proportion of age classes suggests that the population is stable, as the largest proportion of individuals corresponds to the adult stage

Wildlife vaccination strategies for eliminating bovine tuberculosis in white-tailed deer populations.

Many pathogens of humans and livestock also infect wildlife that can act as a reservoir and challenge disease control or elimination. Efficient and effective prioritization of research and management actions requires an understanding of the potential for new tools to improve elimination probability with feasible deployment strategies that can be implemented at scale. Wildlife vaccination is gaining interest as a tool for managing several wildlife diseases. To evaluate the effect of vaccinating white-tailed deer (Odocoileus virginianus), in combination with harvest, in reducing and eliminating bovine tuberculosis from deer populations in Michigan, we developed a mechanistic age-structured disease transmission model for bovine tuberculosis with integrated disease management. We evaluated the impact of pulse vaccination across a range of vaccine properties. Pulse vaccination was effective for reducing disease prevalence rapidly with even low (30%) to moderate (60%) vaccine coverage of the susceptible and exposed deer population and was further improved when combined with increased harvest. The impact of increased harvest depended on the relative strength of transmission modes, i.e., direct vs indirect transmission. Vaccine coverage and efficacy were the most important vaccine properties for reducing and eliminating disease from the local population. By fitting the model to the core endemic area of bovine tuberculosis in Michigan, USA, we identified feasible integrated management strategies involving vaccination and increased harvest that reduced disease prevalence in free-ranging deer. Few scenarios led to disease elimination due to the chronic nature of bovine tuberculosis. A long-term commitment to regular vaccination campaigns, and further research on increasing vaccines efficacy and uptake rate in free-ranging deer are important for disease management.

Moving beyond the mean: an analysis of faecal corticosterone metabolites shows substantial variability both within and across white-tailed deer populations.

Glucocorticoid (GC) levels have significant impacts on the health and behaviour of wildlife populations and are involved in many essential body functions including circadian rhythm, stress physiology and metabolism. However, studies of GCs in wildlife often focus on estimating mean hormone levels in populations, or a subset of a population, rather than on assessing the entire distribution of hormone levels within populations. Additionally, explorations of population GC data are limited due to the tradeoff between the number of individuals included in studies and the amount of data per individual that can be collected. In this study, we explore patterns of GC level distributions in three white-tailed deer (Odocoileus virginianus) populations using a non-invasive, opportunistic sampling approach. GC levels were assessed by measuring faecal corticosterone metabolite levels ('fCMs') from deer faecal samples throughout the year. We found both population and seasonal differences in fCMs but observed similarly shaped fCM distributions in all populations. Specifically, all population fCM cumulative distributions were found to be very heavy-tailed. We developed two toy models of acute corticosterone elevation in an effort to recreate the observed heavy-tailed distributions. We found that, in all three populations, cumulative fCM distributions were better described by an assumption of large, periodic spikes in corticosterone levels every few days, as opposed to an assumption of random spikes in corticosterone levels. The analyses presented in this study demonstrate the potential for exploring population-level patterns of GC levels from random, opportunistically sampled data. When taken together with individual-focused studies of GC levels, such analyses can improve our understanding of how individual hormone production scales up to population-level patterns.

Over-winter body mass and conceptions of white-tailed deer in central Texas

At lower latitudes ungulate reproduction is often seasonal, but the duration of mating and parturition seasons can be long. To further understand the adaptive significance of synchrony in reproductive events, we examined two hypotheses to link over-winter body mass with lower synchrony of conceptions in a wild population of white-tailed deer ( Odocoileus virginianus (Zimmermann, 1780)) in central Texas. If females maintain over-winter body mass, then the window of time that conception occurs should be relaxed and longer relative to females where over-winter body mass declines. Across 9 years, we measured body mass in 1.5-year-old and older females in November ( n = 170) and in March (118). In March, we also measured forehead–rump lengths in offspring of gravid females and calculated conception dates. Over-winter body mass did not decline but was maintained, on average, across the years of study. Conception dates ranged over 50 days and 52% of conceptions occurred in a 2-week period centered on the median conception date. Relative to northern populations of deer, there was lower synchrony in conception dates. White-tailed deer across the geographic range appear to be flexible in over-winter body mass dynamics that affects duration of reproductive events such as conception.

Feeding habits of a white-tailed deer (Odocoileus virginianus cariacou) population as inferred by dental wear at a protected area of the Colombian Orinoquia

Feeding habits of a white-tailed deer (Odocoileus virginianus cariacou) population as inferred by dental wear at a protected area of the Colombian Orinoquia

Population Dynamics of a Declining White-Tailed Deer Population in the Southern Appalachian Region of the United States

Although generally abundant, white-tailed deer (Odocoileus virginianus) populations in the southeastern United States have recently experienced several localized declines attributed to reduced fawn recruitment following the establishment of coyotes (Canis latrans). The Southern Appalachians is a mountainous region suggested to be experiencing white-tailed deer declines, as harvest numbers and hunter success rates have substantially decreased in northern Georgia since 1979. Low fawn survival (16%) was also recently documented in the Chattahoochee National Forest (CNF) in northern Georgia, necessitating further examination. We radio-collared 14 yearling and 45 adult female white-tailed deer along with 71 fawns during 2018–2020 in the CNF to estimate field-based vital rates (i.e., survival and fecundity) and parameterize stage-structured population models. We projected population growth rates (λ) over 10 years to evaluate the current rate of decline and various other management scenarios. Our results indicated that the observed population would decline by an average of 4.0% annually (λ = 0.960) under current conditions. Only scenarios including antlerless harvest restrictions in addition to improved fawn survival resulted in positive growth (λ = 1.019, 1.085), suggesting these measures are likely necessary for population recovery in the region. This approach can be applied by wildlife managers to inform site-specific management strategies.

Impacts of a Severe Hailstorm on a White-Tailed Deer Population in East-Central Alabama

Impacts of a Severe Hailstorm on a White-Tailed Deer Population in East-Central Alabama

White-tailed deer (Odocoileus virginianus) fawn survival and the influence of landscape characteristics on fawn predation risk in the Southern Appalachian Mountains, USA.

In the Southern Appalachian region of the United States, harvest data has indicated the occurrence of low deer densities while exposing a trend of declining white-tailed deer (Odocoileus virginianus) populations over the past several decades in northern Georgia. A triumvirate of increasing fawn predator populations reside in the Southern Appalachian Mountains including coyotes (Canis latrans), black bears (Ursus americanus) and bobcats (Lynx rufus). This region is also characterized by a homogenous landscape composed of mature forests and sparse understory vegetation, likely lacking adequate cover to offer fawns refugia from predators. Our objectives were to estimate survival and cause-specific mortality rates of fawns while assessing a possible link between mortality risk, intrinsic fawn characteristics (i.e., birth mass, Julian birth date, sibling status), and landscape features within fawn usage areas. During 2018-2020, we radio-collared 71 fawns within the Chattahoochee National Forest of northern Georgia, USA and monitored survival to 12 weeks of age. We observed low fawn survival (cumulative = 0.157, 95% CI = 0.091-0.273; vaginal implant transmitter = 0.196, 95% CI = 0.096-0.403) with predation as the leading cause of all known mortalities (45 of 55 mortalities; 82%) due primarily to coyotes (n = 22), black bears (n = 12), and bobcats (n = 7). Relationships between landscape features and fawn predation risk were minimal with only one informative covariate. Increasing amounts of early successional land cover within fawn usage areas decreased fawn mortality risk within the first 20 days of life, but elevated mortality risk thereafter. All fawns with any amount of early successional land cover in their usage areas died of predation (n = 13) at various time intervals, suggesting limited areas of potential fawning cover may be targeted by predators. However, fawn predation risk seemed to be high regardless of landscape covariates due to the limited number of surviving fawns. Coyote-caused mortality occurred over a longer period at a consistently higher magnitude than all other forms of mortality, indicating possible delayed prey-switching behavior and coyote predation as an important factor of fawn survival. The low recruitment of fawns influenced by high predation rates and homogenous habitat conditions is likely the cause of deer population declines in the region.

Accelerated evolution of SARS-CoV-2 in free-ranging white-tailed deer

The zoonotic origin of the COVID-19 pandemic virus highlights the need to fill the vast gaps in our knowledge of SARS-CoV-2 ecology and evolution in non-human hosts. Here, we detected that SARS-CoV-2 was introduced from humans into white-tailed deer more than 30 times in Ohio, USA during November 2021-March 2022. Subsequently, deer-to-deer transmission persisted for 2–8 months, disseminating across hundreds of kilometers. Newly developed Bayesian phylogenetic methods quantified how SARS-CoV-2 evolution is not only three-times faster in white-tailed deer compared to the rate observed in humans but also driven by different mutational biases and selection pressures. The long-term effect of this accelerated evolutionary rate remains to be seen as no critical phenotypic changes were observed in our animal models using white-tailed deer origin viruses. Still, SARS-CoV-2 has transmitted in white-tailed deer populations for a relatively short duration, and the risk of future changes may have serious consequences for humans and livestock.