Behavioral plasticity mitigates the effect of warming on white‐tailed deer

Climate change is expected to create novel environments in which extant species cannot persist, therefore leading to the loss of them and their associated ecological functions within the ecosystem. However, animals may employ behavioral mechanisms in response to warming that could allow them to maintain their functional roles in an ecosystem despite changed temperatures. Specifically, animals may shift their activity in space or time to make use of thermal heterogeneity on the landscape. However, few studies consider the role of behavioral plasticity and spatial or temporal heterogeneity in mitigating the effects of climate change. We conducted experiments to evaluate the potential importance of behavior in mediating the net effects of warming on white‐tailed deer (Odocoileus virginianus). We used shade structures to manipulate the thermal environment around feeding stations to monitor deer feeding activity and measure total consumption. In individual experiments where deer only had access to unshaded feeders, deer fed less during the day but compensated by increasing feeding during times when temperature was lower. In group experiments where deer had access to both shaded and unshaded feeders, deer often fed during the day but disproportionally preferred the cooler, shaded feeders. Our results suggest that deer can capitalize on temporal and spatial heterogeneity in the thermal environment to meet nutritional and thermal requirements, demonstrating the importance of behavioral plasticity when predicting the net effects of climate change.     

Impact of white‐tailed deer on the spread of Borrelia burgdorferi

There is a public perception that the white‐tailed deer Odocoileus virginianus (Artiodactyla: Cervidae) is the main reservoir supporting the maintenance and spread of the causative agent of Lyme disease, Borrelia burgdorferi. This study examines the pathogen prevalence rate of Borrelia in adult Ixodes scapularis (Ixodida: Ixodidae), the black‐legged tick, collected from white‐tailed deer and compares it with pathogen prevalence rates in adult ticks gathered by dragging vegetation in two contiguous counties west of the Hudson Valley in upstate New York. In both Broome and Chenango Counties, attached and unattached ticks harvested from white‐tailed deer had significantly lower prevalences of B. burgdorferi than those collected from vegetation. No attached ticks on deer (n = 148) in either county, and only 2.4 and 7.3% of unattached ticks (n = 389) in Broome and Chenango Counties, respectively, were harbouring the pathogen. This contrasts with the finding that 40.8% of ticks in Broome County and 46.8% of ticks in Chenango County collected from vegetation harboured the pathogen. These data suggest that a mechanism in white‐tailed deer may aid in clearing the pathogen from attached deer ticks, although white‐tailed deer do contribute to the spatial distribution of deer tick populations and also serve as deadend host breeding sites for ticks.     

Active immunization against gonadotropin‐releasing hormone in female white‐tailed deer

The ability of gonadotropin‐releasing hormone (GnRH) immunization to disrupt estrous cycles in captive white‐tailed deer (Odocoileus virginianus) was tested. Four does were each injected subcutaneously with 1 mg of a GnRH analog‐ovalbumin conjugate, using a diethylaminoethyl (DEAE)‐dextran solution as the adjuvant. Control deer (n = 4) received ovalbumin alone in DEAE‐dextran. The immunization schedule consisted of a primary immunization on 30 December 1994, followed by three booster immunizations at 4‐week intervals. In addition, a booster was administered the following year (3 Nov 1995) before the start of the breeding season. Serum from control females did not contain GnRH antibodies, and estrous cycles of these deer were not disrupted (as determined by serum progesterone concentrations). In contrast, all GnRH‐immunized deer had detectable GnRH antibody titers by 1 week after the first booster. Estrous cycles were disrupted in two of these four deer. The booster given the following year failed to stimulate an increase in mean antibody titers, and all deer began to cycle, including one that had maintained high titers from the previous year. These data suggest that a GnRH analog conjugated to ovalbumin, using DEAE‐dextran as adjuvant, is immunogenic in female white‐tailed deer and may be able to prevent ovulation in some individuals. However, pending further work to increase the immunogenic and biological responses and to decrease the variation among animals, this vaccine does not appear to be an effective means of contraception for deer. Zoo Biol 18:385–396, 1999. © 1999 Wiley‐Liss, Inc.     

Beyond the browse line: complex cascade effects mediated by white‐tailed deer

In many ecosystems, browsing of large mammals can affect plant species compositions. However, much less is known about potential above‐ and below‐ground trophic interactions of large browsing mammals. This study focused on the direct and indirect effects of browsing on trophic and abiotic interactions within forest ecosystems. To quantify above‐ and below‐ground cascade effects, white‐tailed deer have been excluded for over 18 years from three 4‐ha plots, which were paired with same sized deer access plots. Our results demonstrate complex direct and indirect cascade effects on forest food webs. Deer exclusion directly altered woody species composition and significantly increased shrub and sapling density. Above‐ground cascade effects include greater leaf litter accumulation and higher arthropod density and biomass within the exclosures. Below‐ground indirect effects include significant decrease in soil nutrients, and higher arbuscular mycorrhizal fungal inoculum potential in the exclosures. Because ecosystems have finite resource availability, high deer density may imbalance the system by redirecting resources toward maintaining deer biomass at the expense of multiple trophic levels throughout the forest community. Both complex bottom up and top down trophic cascade effects demonstrated largely unidirectional negative responses suggesting that high deer density has reduced the biodiversity of the forest community.     

Breeding behavior of female white‐tailed deer relative to conception: Evidence for female mate choice

Female white‐tailed deer (Odocoileus virginianus) are thought to choose between two behavioral strategies to maximize the quality of potential mates: sit and wait, characterized by concentrating activity within a restricted area, and excursive behavior, characterized by increased activity and excursions outside the home range. As movement patterns may influence conception, our goal was to examine the patterns of female white‐tailed deer movements to evaluate which breeding strategy was employed. We equipped 36 female white‐tailed deer with GPS collars from August 2013 to December 2015. We found that movement rate and probability of activity were greatest near the peak of the breeding season, and we observed increases in both metrics during the 40 days prior to estimated conception. Peak size of home range and core area occurred in the days surrounding conception. We found that 11 deer performed an excursion, ranging from 43 days before until 36 days after conception, with the peak probability of being outside of an individual home range occurring 1 day prior to conception. Our results suggest that female white‐tailed deer may attempt to maximize the quality of their mates by advertising availability for breeding through excursive behaviors just prior to entering estrus.     

The effect of terrain and female density on survival of neonatal white‐tailed deer and mule deer fawns

Juvenile survival is a highly variable life‐history trait that is critical to population growth. Antipredator tactics, including an animal's use of its physical and social environment, are critical to juvenile survival. Here, we tested the hypothesis that habitat and social characteristics influence coyote (Canis latrans) predation on white‐tailed deer (Odocoileus virginianus) and mule deer (O. hemionus) fawns in similar ways during the neonatal period. This would contrast to winter when the habitat and social characteristics that provide the most safety for each species differ. We monitored seven cohorts of white‐tailed deer and mule deer fawns at a grassland study site in Alberta, Canada. We used logistic regression and a model selection procedure to determine how habitat characteristics, climatic conditions, and female density influenced fawn survival during the first 8 weeks of life. Fawn survival improved after springs with productive vegetation (high integrated Normalized Difference Vegetation Index values). Fawns that used steeper terrain were more likely to survive. Fawns of both species had improved survival in years with higher densities of mule deer females, but not with higher densities of white‐tailed deer females, as predicted if they benefit from protection by mule deer. Our results suggest that topographical variation is a critical resource for neonates of many ungulate species, even species like white‐tailed deer that use more gentle terrain when older. Further, our results raise the possibility that neonatal white‐tailed fawns may benefit from associating with mule deer females, which may contribute to the expansion of white‐tailed deer into areas occupied by mule deer.     

Addressing regional relationships between white‐tailed deer densities and land classes

White‐tailed deer (Odocoileus virginianus) populations have recovered to about 30 million animals in the United States, but land cover has changed during the interval of recovery. To address the relationship between deer densities and current land cover at regional scales, I applied random forests and extreme gradient boosting classifiers to model low and high deer density classes, at two different thresholds (5.8 and 11.6 deer/km²), and land classes in three regions during approximately 2003. For low and high deer density classes divided at 5.8 deer/km², deciduous broadleaf forest overall was the most influential and positive variable in the central east and central regions and crop and pasture were the most influential and negative variables in the southeast region. Deer density increased with area of deciduous and mixed forests, woody wetlands, and shrub in all regions. Deer density decreased with area of crop, developed open space, and developed low and medium residential density in all regions. For density classes divided at 11.6 deer/km², deer density had the strongest relationship with woody wetlands in the central east region, mixed and deciduous forest in the southeast region, and woody wetlands and herbaceous vegetation in the central region. Deer density increased with deciduous and mixed forests, woody wetlands, and shrub in all regions. Conversely, deer density decreased with herbaceous vegetation, crop, and developed low residential densities in all regions. Therefore, at regional scales, deer overall occurred at greater densities in forests and woody wetlands and lower densities in agricultural and residential development, which did not appear to support more deer. Deer preference for forests does result in damage to forest products, but alternatively, some may consider that deer provide important socioeconomic and ecological services by reducing number of small trees, particularly in the absence of other disturbances that historically controlled tree biomass.     

Red oak seedlings as indicators of deer browse pressure: Gauging the outcome of different white‐tailed deer management approaches

After decades of high deer populations, North American forests have lost much of their previous biodiversity. Any landscape‐level recovery requires substantial reductions in deer herds, but modern societies and wildlife management agencies appear unable to devise appropriate solutions to this chronic ecological and human health crisis. We evaluated the effectiveness of fertility control and hunting in reducing deer impacts at Cornell University. We estimated spring deer populations and planted Quercus rubra seedlings to assess deer browse pressure, rodent attack, and other factors compromising seedling performance. Oak seedlings protected in cages grew well, but deer annually browsed ≥60% of unprotected seedlings. Despite female sterilization rates of >90%, the deer population remained stable. Neither sterilization nor recreational hunting reduced deer browse rates and neither appears able to achieve reductions in deer populations or their impacts. We eliminated deer sterilization and recreational hunting in a core management area in favor of allowing volunteer archers to shoot deer over bait, including at night. This resulted in a substantial reduction in the deer population and a linear decline in browse rates as a function of spring deer abundance. Public trust stewardship of North American landscapes will require a fundamental overhaul in deer management to provide for a brighter future, and oak seedlings may be a promising metric to assess success. These changes will require intense public debate and may require new approaches such as regulated commercial hunting, natural dispersal, or intentional release of important deer predators (e.g., wolves and mountain lions). Such drastic changes in deer management will be highly controversial, and at present, likely difficult to implement in North America. However, the future of our forest ecosystems and their associated biodiversity will depend on evidence to guide change in landscape management and stewardship.     

Estimating preharvest density,adult sex ratio,and fecundity of white‐tailed deer using noninvasive sampling techniques

Adult sex ratio and fecundity (juveniles per female) are key population parameters in sustainable wildlife management, but inferring these requires abundance estimates of at least three age/sex classes of the population (male and female adults and juveniles). Prior to harvest, we used an array of 36 wildlife camera traps during 2 and 3 weeks in the early autumn of 2016 and 2017, respectively. We recorded white‐tailed deer adult males, adult females, and fawns from the pictures. Simultaneously, we collected fecal DNA (fDNA) from 92 20 m × 20 m plots placed in 23 clusters of four plots between the camera traps. We identified individuals from fDNA samples with microsatellite markers and estimated the total sex ratio and population density using spatial capture–recapture (SCR). The fDNA‐SCR analysis concluded equal sex ratio in the first year and female bias in the second year, and no difference in space use between sexes (fawns and adults combined). Camera information was analyzed in a spatial capture (SC) framework assuming an informative prior for animals’ space use, either (a) as estimated by fDNA‐SCR (same for all age/sex classes), (b) as assumed from the literature (space use of adult males larger than adult females and fawns), or (c) by inferring adult male space use from individually identified males from the camera pictures. These various SC approaches produced plausible inferences on fecundity, but also inferred total density to be lower than the estimate provided by fDNA‐SCR in one of the study years. SC approaches where adult male and female were allowed to differ in their space use suggested the population had a female‐biased adult sex ratio. In conclusion, SC approaches allowed estimating the preharvest population parameters of interest and provided conservative density estimates.     

Spatial correlations between browsing on balsam fir by white‐tailed deer and the nutritional value of neighboring winter forage

Associational effects, that is, the influence of neighboring plants on herbivory suffered by a plant, are an outcome of forage selection. Although forage selection is a hierarchical process, few studies have investigated associational effects at multiple spatial scales. Because the nutritional quality of plants can be spatially structured, it might differently influence associational effects across multiple scales. Our objective was to determine the radius of influence of neighbor density and nutritional quality on balsam fir (Abies balsamea) herbivory by white‐tailed deer (Odocoileus virginianus) in winter. We quantified browsing rates on fir and the density and quality of neighboring trees in a series of 10‐year‐old cutovers on Anticosti Island (Canada). We used cross‐correlations to investigate relationships between browsing rates and the density and nutritional quality of neighboring trees at distances up to 1,000 m. Balsam fir and white spruce (Picea glauca) fiber content and dry matter in vitro true digestibility were correlated with fir browsing rate at the finest extra‐patch scale (across distance of up to 50 m) and between cutover areas (300–400 m). These correlations suggest associational effects, that is, low nutritional quality of neighbors reduces the likelihood of fir herbivory (associational defense). Our results may indicate associational effects mediated by intraspecific variation in plant quality and suggest that these effects could occur at scales from tens to hundreds of meters. Understanding associational effects could inform strategies for restoration or conservation; for example, planting of fir among existing natural regeneration could be concentrated in areas of low nutritional quality.     

Patterns of testosterone in male white‐tailed deer (Odocoileus virginianus): Seasonal and lifetime variation

Testosterone is strongly associated with the annual development of antlers in cervids, but endocrine research on wild, freely breeding ungulates is often done without repeated capture of known‐aged individuals. As a result, our knowledge on how testosterone fluctuates over the course of a lifetime and variation in lifetime patterns among individuals is limited. We investigated patterns of testosterone in a freely breeding population of white‐tailed deer (Odocoileus virginianus) in Alabama, USA, that breeds in January. Testosterone peaked during the height of the breeding season, despite this period occurring approximately 2 months later than in most temperate‐region white‐tailed deer populations. Age‐related differences in testosterone were only prevalent during the breeding season, with bucks ≥3.5 years old having greater testosterone (853 ng/dl ± 96 SE; p = 0.012) than bucks 1.5–2.5 years old (364 ng/dl ± 100 SE). Additionally, an individual''s testosterone level as a yearling was not positively associated with their lifetime maximum testosterone level (p = 0.583), and an individual''s mean testosterone level was positively associated with lifetime testosterone variation (p < 0.001). To our knowledge, our study is one of the first to assess how testosterone early in life might relate to individual testosterone later in life. We believe these data provide insight into lifetime hormonal patterns in cervids, and that these patterns may indicate intraspecific variation of lifetime reproductive strategies.     

Modelling geographic patterns of population density of the white‐tailed deer in central Mexico by implementing ecological niche theory

Conservation and management of species require basic knowledge on their geographic distribution and abundance. Here, we propose a novel approach, based on the theory of the ecological niche, to model the spatial patterns of the white‐tailed deer Odocoileus virginianus population density in two regions of central Mexico (Balsas Basin and Tehuacán‐Cuicatlán Valley). We used an ecological niche model to generate binary geographic distribution maps of the white‐tailed deer in each region based on occurrence data and a set of environmental variables. Then, the centroid of the distributions was calculated in ecological space (niche centroid) and the multidimensional Euclidian ecological distance of each pixel to the niche centroid was estimated. Finally, for each region the distance to the niche centroid (DNC) was regressed against 14 independent occurrence points in each site containing white‐tailed deer density information to determine the function describing the DNC‐density relationship, which was used to generate maps describing the distribution of white‐tailed deer density. Our results indicated an inverse DNC‐density relationship in both regions (Balsas Basin: r² = 0.90 and Tehuacán‐Cuicatlán: r² = 0.76) that was validated via bootstrapping resulting in a predicting capacity of near 62% for Balsas Basin and 65% for Tehuacán‐Cuicatlán Valley. Our results suggest that the distance to the niche centroid method is a robust, science‐based correlative approach that resulted useful to predict the population density of the white‐tailed deer in a spatially explicit fashion. The proposed approach is suitable for predicting the distribution of density for white‐tailed deer for which occurrence data with accompanying density information exists, but relative abundance can also be estimated when no abundance data are available.     

Absence of founder effect and evidence for adaptive divergence in a recently introduced insular population of white‐tailed deer (Odocoileus virginianus)

Islands are generally colonized by few individuals which could lead to a founder effect causing loss of genetic diversity and rapid divergence by strong genetic drift. Insular conditions can also induce new selective pressures on populations. Here, we investigated the extent of genetic differentiation within a white‐tailed deer (Odocoileus virginianus) population introduced on an island and its differentiation with its source mainland population. In response to their novel environmental conditions, introduced deer changed phenotypically from mainland individuals, therefore we investigated the genetic bases of the morphological differentiation. The study was conducted on Anticosti Island (Québec, Canada) where 220 individuals were introduced 120 years ago, resulting in a population size over 160,000 individuals. We used genotyping‐by‐sequencing (GBS) to generate 8,518 filtered high‐quality SNPs and compared patterns of genetic diversity and differentiation between the continental and Anticosti Island populations. Clustering analyses indicated a single panmictic island population and no sign of isolation by distance. Our results revealed a weak, albeit highly significant, genetic differentiation between the Anticosti Island population and its source population (mean F_ST = 0.005), which allowed a population assignment success of 93%. Also, the high genetic diversity maintained in the introduced population supports the absence of a strong founder effect due to the large number of founders followed by rapid population growth. We further used a polygenic approach to assess the genetic bases of the divergent phenotypical traits between insular and continental populations. We found loci related to muscular function and lipid metabolism, which suggested that these could be involved in local adaptation on Anticosti Island. We discuss these results in a harvest management context.     

Using heterozygosity–fitness correlations to study inbreeding depression in an isolated population of white‐tailed deer founded by few individuals

A heterozygosity–fitness correlations (HFCs) may reflect inbreeding depression, but the extent to which they do so is debated. HFCs are particularly likely to occur after demographic disturbances such as population bottleneck or admixture. We here study HFC in an introduced and isolated ungulate population of white‐tailed deer Odocoileus virginianus in Finland founded in 1934 by four individuals. A total of 422 ≥ 1‐year‐old white‐tailed deer were collected in the 2012 hunting season in southern Finland and genotyped for 14 microsatellite loci. We find significant identity disequilibrium as estimated by g₂. Heterozygosity was positively associated with size‐ and age‐corrected body mass, but not with jaw size or (in males) antler score. Because of the relatively high identity disequilibrium, heterozygosity of the marker panel explained 51% of variation in inbreeding. Inbreeding explained approximately 4% of the variation in body mass and is thus a minor, although significant source of variation in body mass in this population. The study of HFC is attractive for game‐ and conservation‐oriented wildlife management because it presents an affordable and readily used approach for genetic monitoring that allowing identification of fitness costs associated with genetic substructuring in what may seem like a homogeneous population.     

Cyclic voles, prey switching in red fox, and roe deer dynamics – a test of the alternative prey hypothesis

Medium-sized predators sometimes switch to alternative prey species as their main prey declines. Our objective of this study was to test the alternative prey hypothesis for a medium sized predator (red fox, Vulpes vulpes ), a small cyclically fluctuating main prey (microtine voles) and larger alternative prey (roe deer fawns, Capreolus capreolus ). We used long-term time series (28 years) on voles, red fox and roe deer from the Grimsö Wildlife Research Area (59°40'N, 15°25'E) in south-central Sweden to investigate interspecific relationships in the annual fluctuations in numbers of the studied species. Annual variation in number of roe deer fawns in autumn was significantly and positively related to vole density and significantly and negatively related to the number of fox litters in the previous year. In years of high vole density, predation on roe deer fawns was small, but in years of low vole density predation was more severe. The time lag between number of fox litters and predation on fawns was due to the time lag in functional response of red fox in relation to voles. This study demonstrates for the first time that the alternative prey hypothesis is applicable to the system red fox, voles and roe deer fawns.     

Comparing habitat quality within and between environments using giving up densities: an example based on the winter habitat of white‐tailed deer Odocoileus virginianus

We previously developed a model, based on the precepts of optimal patch use, to compare habitat quality both within and between environments. Here we illustrate the use of this model in a study estimating quality of winter habitats (deer yards) of white‐tailed deer Odocoileus virginianus near the northern limit of their range by following their foraging behaviour. We compare giving up densities (GUDs), the amount of food remaining in a patch when a forager ceases foraging there, with and without the presence of supplemental food in order to draw inferences about the relative quality either of habitats within an environment or of distinct environments. We use our model to evaluate the impact of alterations to the winter habitat of deer at two distinct sites that differed in their level of predation risk and food availability. The first site, the Mont Rigaud deer yard, was surrounded by farm land and gradually‐expanding suburbs. Predators were rare and food was available in winter either in farm fields or around private homes but deer browsing in the past had left little food in the forest. At the second site, the Calumet deer yard, deer experienced a higher predation risk and did not have access to supplemental food from farm fields or private homes. However, past browsing by deer had not drastically reduced food in the forest. We offered food to deer in four habitats per site (forest, forest edge, clearing, clearing edge) with four to six replicates per site and measured the GUDs after 24 h. Analysis of these data, interpreted according to our model, suggests that deer are more sensitive to metabolic costs at Mont Rigaud and food availability at Calumet; predation risk does not alter deer behaviour between the two sites. Within habitats, deer at Mont Rigaud reacted to clearings as though they imposed higher metabolic costs than the forest. They also reacted to an interaction in which missed food costs influenced GUD only when metabolic costs were not too high. Thus our model appears to provide a convenient tool for comparing habitat quality both within and between environments.     

Genetic structure and phylogeography of a European flagship species,the white‐tailed sea eagle Haliaeetus albicilla

We analysed 120 white‐tailed sea eagles Haliaeetus albicilla from eastern (Poland and Estonia) and southeastern (Serbian Danube population) Europe for genetic variability and structuring at the mitochondrial control region and seven nuclear microsatellite loci. We combined this new dataset with sequence and genotype data from previous analyses covering Greenland and Eurasia (total sample sizes of 420 and 186 individuals for mtDNA and microsatellites, respectively) to address the following questions: 1) does the large eastern population in Europe add significantly to the species‘ overall genetic diversity? 2) Do the new sequence data match the clinal distribution pattern (west to east) of the two major mtDNA lineages? 3) Does the preliminary hypothesis of two nuclear genetic clusters recently found in this species hold for the whole of Europe, and do these clusters show a geographic pattern? Our results confirmed Europe as a stronghold of genetic diversity in white‐tailed sea eagles, and the east of the continent contributed disproportionately to this, the reason being the admixture of eagles with different genetic background. As hypothesised, both mitochondrial lineages were recovered also in eastern Europe, but the globally more eastern lineage was dominant. The presence of two microsatellite clusters was also confirmed, and these groups, too, show a non‐random geographic distribution, with, except for Poland, a high proportion of ‘eastern‐type’ eagles in the populations of east–central and eastern Europe.