Large-scale assessment of genetic structure to assess risk of populations of a large herbivore to disease.

Chronic wasting disease (CWD) can spread among cervids by direct and indirect transmission, the former being more likely in emerging areas. Identifying subpopulations allows the delineation of focal areas to target for intervention. We aimed to assess the population structure of white-tailed deer (Odocoileus virginianus) in the northeastern United States at a regional scale to inform managers regarding gene flow throughout the region. We genotyped 10 microsatellites in 5701 wild deer samples from Maryland, New York, Ohio, Pennsylvania, and Virginia. We evaluated the distribution of genetic variability through spatial principal component analysis and inferred genetic structure using non-spatial and spatial Bayesian clustering algorithms (BCAs). We simulated populations representing each inferred wild cluster, wild deer in each state and each physiographic province, total wild population, and a captive population. We conducted genetic assignment tests using these potential sources, calculating the probability of samples being correctly assigned to their origin. Non-spatial BCA identified two clusters across the region, while spatial BCA suggested a maximum of nine clusters. Assignment tests correctly placed deer into captive or wild origin in most cases (94%), as previously reported, but performance varied when assigning wild deer to more specific origins. Assignments to clusters inferred via non-spatial BCA performed well, but efficiency was greatly reduced when assigning samples to clusters inferred via spatial BCA. Differences between spatial BCA clusters are not strong enough to make assignment tests a reliable method for inferring the geographic origin of deer using 10 microsatellites. However, the genetic distinction between clusters may indicate natural and anthropogenic barriers of interest for management.

Species and physiographic factors drive Indian cucumber root and Canada mayflower plant chemistry: Implications for white-tailed deer forage quality.

Nutrition is fundamental to white-tailed deer (Odocoileus virginianus) management given its relationship to habitat carrying capacity and population productivity. Ecological Sites (ESs) are a United States federal landscape management unit of specific land potential due to unique soils, topography, climate, parent material, and perhaps deer forage nutritional value. We present results of a study that extends the use of ESs to inform white-tailed deer management by evaluating indicator plant chemistry in two spring forb species, Indian cucumber root (Medeola virginiana) and Canada mayflower (Maianthemum canadense), across the northcentral Appalachians. We sampled spring forbs and underlying soils across two ESs: Dry, upland, oak-maple-hemlock hardwood forest (OMH) and Deep soil, high slope, northern hardwood forests (NHF). Plant elemental content, soil pH, and site aspect, slope and elevation were measured. Our results show that forb chemistry differs between species and within a species geographically. Indian cucumber root, as compared to Canada mayflower, has significantly higher Mg, Na, Cu, Fe, and Zn, and lower Mn. Canada mayflower in the NHF ES, versus OMH ES, was found to have significantly higher K, Mn, and B. Indian cucumber root in the NHF ES, versus the OMH ES, was found to have significantly higher Mg, Al, Fe, and Ca:P ratio but lower K. Linear discriminant analysis shows that plant tissue Mn was the best discriminator between species, and between ESs, Canada mayflower plant tissue Mn and Indian cucumber plant tissue P, K, Ca, Mg and Mn were best discriminators. Given that nutrition determines habitat carrying capacity, differences in forage nutrition between ESs may have different potentials to support deer. Forage nutrition is an important aspect of deer habitat conditions and carrying capacity, thus ESs are likely to support deer populations with different growth potential, which means that even if the same plant species occur in different ESs their nutritional value to deer may differ.

Variability in prion protein genotypes by spatial unit to inform susceptibility to chronic wasting disease.

Chronic wasting disease (CWD) is a fatal encephalopathy affecting North American cervids. Certain alleles in a host's prion protein gene are responsible for reduced susceptibility to CWD. We assessed for the first time variability in the prion protein gene of elk (Cervus canadensis) present in Pennsylvania, United States of America, a reintroduced population for which CWD cases have never been reported. We sequenced the prion protein gene (PRNP) of 565 elk samples collected over 7 years (2014-2020) and found two polymorphic sites (codon 21 and codon 132). The allele associated with reduced susceptibility to CWD is present in the population, and there was no evidence of deviations from Hardy-Weinberg equilibrium in any of our sampling years (p-values between 0.14 and 1), consistent with the lack of selective pressure on the PRNP. The less susceptible genotypes were found in a frequency similar to the ones reported for elk populations in the states of Wyoming and South Dakota before CWD was detected. We calculated the proportion of less susceptible genotypes in each hunt zone in Pennsylvania as a proxy for their vulnerability to the establishment of CWD, and interpolated these results to obtain a surface representing expected proportion of the less susceptible genotypes across the area. Based on this analysis, hunt zones located in the southern part of our study area have a low proportion of less susceptible genotypes, which is discouraging for elk persistence in Pennsylvania given that these hunt zones are adjacent to the deer Disease Management Area 3, where CWD has been present since 2014.

Caution is warranted when using animal space-use and movement to infer behavioral states.

BACKGROUND: Identifying the behavioral state for wild animals that can't be directly observed is of growing interest to the ecological community. Advances in telemetry technology and statistical methodologies allow researchers to use space-use and movement metrics to infer the underlying, latent, behavioral state of an animal without direct observations. For example, researchers studying ungulate ecology have started using these methods to quantify behaviors related to mating strategies. However, little work has been done to determine if assumed behaviors inferred from movement and space-use patterns correspond to actual behaviors of individuals. METHODS: Using a dataset with male and female white-tailed deer location data, we evaluated the ability of these two methods to correctly identify male-female interaction events (MFIEs). We identified MFIEs using the proximity of their locations in space as indicators of when mating could have occurred. We then tested the ability of utilization distributions (UDs) and hidden Markov models (HMMs) rendered with single sex location data to identify these events. RESULTS: For white-tailed deer, male and female space-use and movement behavior did not vary consistently when with a potential mate. There was no evidence that a probability contour threshold based on UD volume applied to an individual's UD could be used to identify MFIEs. Additionally, HMMs were unable to identify MFIEs, as single MFIEs were often split across multiple states and the primary state of each MFIE was not consistent across events. CONCLUSIONS: Caution is warranted when interpreting behavioral insights rendered from statistical models applied to location data, particularly when there is no form of validation data. For these models to detect latent behaviors, the individual needs to exhibit a consistently different type of space-use and movement when engaged in the behavior. Unvalidated assumptions about that relationship may lead to incorrect inference about mating strategies or other behaviors.

Risky movements? Natal dispersal does not decrease survival of a large herbivore.

Natal dispersal is assumed to be a particularly risky movement behavior as individuals transfer, often long distances, from birth site to site of potential first reproduction. Though, because this behavior persists in populations, it is assumed that dispersal increases the fitness of individuals despite the potential for increased risk of mortality. The extent of dispersal risk, however, has rarely been tested, especially for large mammals. Therefore, we aimed to test the relationship between dispersal and survival for both males and females in a large herbivore. Using a radio-transmittered sample of 398 juvenile male and 276 juvenile female white-tailed deer (Odocoileus virginianus), we compared survival rates of dispersers and nondispersers. We predicted that dispersing deer would experience greater overall mortality than philopatric deer due to direct transfer-related risks (e.g., vehicular collision), indirect immigration-related mortality attributable to colonization of unfamiliar habitat, and increased overwinter mortality associated with energetic costs of movement and unfamiliarity with recently colonized habitat. For both male and female yearlings, survival rates of dispersers (male = 49.9%, female = 64.0%) did not differ from nondispersers (male = 51.6%, female = 70.7%). Only two individuals (both female) were killed by vehicular collision during transfer, and overwinter survival patterns were similar between the two groups. Although dispersal movement likely incurs energetic costs on dispersers, these costs do not necessarily translate to decreased survival. In many species, including white-tailed deer, dispersal is likely condition-dependent, such that larger and healthier individuals are more likely to disperse; therefore, costs associated with dispersal are more likely to be borne successfully by those individuals that do disperse. Whether low-risk dispersal of large mammals is the rule or the exception will require additional research. Further, future research is needed to evaluate nonsurvival fitness-related costs and benefits of dispersal (e.g., increased reproductive opportunities for dispersers).

Soil chemistry, and not short-term (1-2 year) deer exclusion, explains understory plant occupancy in forests affected by acid deposition.

The loss of species diversity and plant community structure throughout the temperate deciduous forests of North America have often been attributed to overbrowsing by white-tailed deer (Odocoileus virginanus). Slow species recovery following removal from browsing, or reduction in deer density, has been termed a legacy effect of past deer herbivory. However, vegetation legacy effects have also coincided with changes to soil chemistry throughout the north-eastern USA. In this paper, we assess the viability of soil chemistry (i.e. pH, extractable nutrients and extractable metals) and other factors (topography, light, overstory basal area and location) as alternative explanations for a lack of vegetation recovery. We compared the relative effects of soil chemistry, site conditions and short-term (1-2 year) deer exclusion on single-species occupancy probabilities of 10 plant taxa common to oak-hickory forests in central Pennsylvania. We found detection for all modelled species was constant and high ( p ^ > 0.65), and occupancy probability of most taxa was best explained by at least one soil chemistry parameter. Specifically, ericaceous competing vegetation was more likely to occupy acidic (pH < 3.5), base cation-poor (K < 0.20 cmolc kg-1) sites, while deer-preferred plants were less likely to occur when soil manganese exceeded 0.1 cmolc kg-1. Short-term deer exclusion did not explain occupancy of any plant taxon, and site conditions were of nominal importance. This study demonstrates the importance of soil chemistry in shaping plant community composition in the north-central Appalachians, and suggests soil as an alternative, or additional, explanation for deer vegetation legacy effects. We suggest that the reliance on phyto-indicators of deer browsing effects may overestimate the effects of browsing if those species are also limited by unfavourable soil conditions. Future research should consider study designs that address the complexity of deer forest interactions, especially in areas with complex site-vegetation histories.

Heterogeneity of a landscape influences size of home range in a North American cervid.

In the northeastern United States, chronic wasting disease has recently been detected in white-tailed deer (Odocoileus virginianus) populations, and understanding the relationship between landscape configuration and home range may improve disease surveillance and containment efforts. The objectives of our study were to compare size of home range for deer occupying a continuum of forested landscapes and to investigate relationships between size of home range and measures of landscape configuration. We used a movement-based kernel density estimator to estimate home range at five spatial scales among deer across study areas. We developed 7 linear regression models that used measures of the configuration of the forested landscape to explain size of home range. We observed differences in size of home range between sexes among areas that differed based on landscape configuration. We documented size of home range changed with various metrics that identifying connectivity of forested patches. Generally, size of home range increased with an increasing proportion of homogenous forest. Our results suggest that deer in our region occupy a landscape at hierarchically-nested scales that is controlled by the connectivity of the forested landscape across local or broad geographical regions.

Accounting for tagging-to-harvest mortality in a Brownie tag-recovery model by incorporating radio-telemetry data.

The Brownie tag-recovery model is useful for estimating harvest rates but assumes all tagged individuals survive to the first hunting season; otherwise, mortality between time of tagging and the hunting season will cause the Brownie estimator to be negatively biased. Alternatively, fitting animals with radio transmitters can be used to accurately estimate harvest rate but may be more costly. We developed a joint model to estimate harvest and annual survival rates that combines known-fate data from animals fitted with transmitters to estimate the probability of surviving the period from capture to the first hunting season, and data from reward-tagged animals in a Brownie tag-recovery model. We evaluated bias and precision of the joint estimator, and how to optimally allocate effort between animals fitted with radio transmitters and inexpensive ear tags or leg bands. Tagging-to-harvest survival rates from >20 individuals with radio transmitters combined with 50-100 reward tags resulted in an unbiased and precise estimator of harvest rates. In addition, the joint model can test whether transmitters affect an individual's probability of being harvested. We illustrate application of the model using data from wild turkey, Meleagris gallapavo, to estimate harvest rates, and data from white-tailed deer, Odocoileus virginianus, to evaluate whether the presence of a visible radio transmitter is related to the probability of a deer being harvested. The joint known-fate tag-recovery model eliminates the requirement to capture and mark animals immediately prior to the hunting season to obtain accurate and precise estimates of harvest rate. In addition, the joint model can assess whether marking animals with radio transmitters affects the individual's probability of being harvested, caused by hunter selectivity or changes in a marked animal's behavior.

The relationship between species detection probability and local extinction probability.

In community-level ecological studies, generally not all species present in sampled areas are detected. Many authors have proposed the use of estimation methods that allow detection probabilities that are <1 and that are heterogeneous among species. These methods can also be used to estimate community-dynamic parameters such as species local extinction probability and turnover rates (Nichols et al. Ecol Appl 8:1213-1225; Conserv Biol 12:1390-1398). Here, we present an ad hoc approach to estimating community-level vital rates in the presence of joint heterogeneity of detection probabilities and vital rates. The method consists of partitioning the number of species into two groups using the detection frequencies and then estimating vital rates (e.g., local extinction probabilities) for each group. Estimators from each group are combined in a weighted estimator of vital rates that accounts for the effect of heterogeneity. Using data from the North American Breeding Bird Survey, we computed such estimates and tested the hypothesis that detection probabilities and local extinction probabilities were negatively related. Our analyses support the hypothesis that species detection probability covaries negatively with local probability of extinction and turnover rates. A simulation study was conducted to assess the performance of vital parameter estimators as well as other estimators relevant to questions about heterogeneity, such as coefficient of variation of detection probabilities and proportion of species in each group. Both the weighted estimator suggested in this paper and the original unweighted estimator for local extinction probability performed fairly well and provided no basis for preferring one to the other.