Spatiotemporal variation in the fecal microbiota of mule deer is associated with proximate and future measures of host health.

BACKGROUND: Mule deer rely on fat and protein stored prior to the winter season as an energy source during the winter months when other food sources are sparse. Since associated microorganisms ('microbiota') play a significant role in nutrient metabolism of their hosts, we predicted that variation in the microbiota might be associated with nutrient storage and overwintering in mule deer populations. To test this hypothesis we performed a 16S rRNA marker gene survey of fecal samples from two deer populations in the western United States before and after onset of winter. RESULTS: PERMANOVA analysis revealed the deer microbiota varied interactively with geography and season. Further, using metadata collected at the time of sampling, we were able to identify different fecal bacterial taxa that could potentially act as bioindicators of mule deer health outcomes. First, we identified the abundance of Collinsella (family: Coriobacteriaceae) reads as a possible predictor of poor overwintering outcomes for deer herds in multiple locations. Second, we showed that reads assigned to the Bacteroides and Mollicutes Order RF39 were both positively correlated with deer protein levels, leading to the idea that these sequences might be useful in predicting mule deer protein storage. CONCLUSIONS: These analyses confirm that variation in the microbiota is associated with season-dependent health outcomes in mule deer, which may have useful implications for herd management strategies.

Spatial processes decouple management from objectives in a heterogeneous landscape: predator control as a case study.

Predator control is often implemented with the intent of disrupting top-down regulation in sensitive prey populations. However, ambiguity surrounding the efficacy of predator management, as well as the strength of top-down effects of predators in general, is often exacerbated by the spatially implicit analytical approaches used in assessing data with explicit spatial structure. Here, we highlight the importance of considering spatial context in the case of a predator control study in south-central Utah. We assessed the spatial match between aerial removal risk in coyotes (Canis latrans) and mule deer (Odocoileus hemionus) resource selection during parturition using a spatially explicit, multi-level Bayesian model. With our model, we were able to evaluate spatial congruence between management action (i.e., coyote removal) and objective (i.e., parturient deer site selection) at two distinct scales: the level of the management unit and the individual coyote removal. In the case of the former, our results indicated substantial spatial heterogeneity in expected congruence between removal risk and parturient deer site selection across large areas, and is a reflection of logistical constraints acting on the management strategy and differences in space use between the two species. At the level of the individual removal, we demonstrated that the potential management benefits of a removed coyote were highly variable across all individuals removed and in many cases, spatially distinct from parturient deer resource selection. Our methods and results provide a means of evaluating where we might anticipate an impact of predator control, while emphasizing the need to weight individual removals based on spatial proximity to management objectives in any assessment of large-scale predator control. Although we highlight the importance of spatial context in assessments of predator control strategy, we believe our methods are readily generalizable in any management or large-scale experimental framework where spatial context is likely an important driver of outcomes.

Temporal patterns of ungulate herbivory and phenology of aspen regeneration and defense.

Ungulate herbivory can create strong top-down effects on forest recruitment, especially after fire. Defense strategies of tree species against ungulate herbivory include escape through vertical growth and resistance through the production of defense compounds. Using a four-way fence design and camera traps we characterized the differential impacts of ungulate herbivores (deer, elk, cattle) on aspen forest recruitment and plant defense responses and how they vary depending on the timing of herbivory. Aspen height growth was greatest between June and August and ungulate use of aspen was highest in July and August. Three years after fire, height of aspen differed among fence treatments with full ungulate exclusion > deer-only plots > native ungulate plots > fenceless plots: 108 ± 4 cm, 94 ± 4 cm, 89 ± 4 cm, and 65 ± 4 cm, respectively. Fenceless plots had the highest rates of removal of apical meristems by the end of 2014 and 2015 (61% and 53%, respectively). Native ungulate plots, and deer-only plots both had similar removal of apical meristems in 2014 (37% and 39%, respectively). The highest phenolic glycoside concentrations were associated with an 80% reduction in meristem removal and four-fold greater aspen height by the end of summer. Low nitrogen was associated with an 86% reduction in apical meristem removal and threefold greater aspen height. In conclusion, our study suggests that high ungulate abundance can have detrimental impacts on forest recruitment and that high aspen defense chemistry and lower leaf N deters ungulate herbivory, especially in the late summer.

Evaluating movement-based methods for estimating the frequency and timing of parturition in mule deer.

BACKGROUND: Information on reproduction of harvested species such as mule deer (Odocoileus hemionus) is vital for conservation and management. Furthermore, parturition in ungulates may be detected using patterns of movement logged by GPS transmitters. Several movement-based methods have been developed to detect parturition in ungulates including the Peterson method, behavioral change point analysis (BCPA), rolling minimum convex polygons (rMCP), individual-based method (IBM), and population-based method (PBM). Our objectives were to (1) test the accuracy and the precision of each previously described method and (2) develop an improved method optimized for mule deer that incorporated aspects of the other methods. METHODS: We determined parturition timing and status for female mule deer fitted with GPS collars and implanted with vaginal implant transmitters (VITs). We used movement patterns before and after parturition to set movement thresholds for each movement-based method. Following model training, we used location and birth date data from an external dataset to test the effectiveness of each movement-based method. Additionally, we developed a novel method for detecting parturition called the analysis of parturition indicators (API). We used two regression analyses to determine the accuracy and precision of estimates generated by each method. RESULTS: The six methods we employed varied in accuracy, with the API, rMCP, and BCPA being most accurate. Precision also varied among methods, with the API, rMCP, and PBM generating the most precise estimates of parturition dates. The API and the rMCP performed similarly and better overall than any of the other existing methods. CONCLUSIONS: We found that movement-based methods could be used to accurately and precisely detect parturition in mule deer. Further, we determined that the API and rMCP methods had the greatest overall success at detecting parturition in mule deer. The relative success of the API and rMCP may be attributed to the fact that both methods use home range size to detect parturition and are validated using known parturition dates of collared deer. We present the API as an efficient method of estimating birth status and timing of parturition of mule deer fitted with GPS transmitters, as well as affirm the effectiveness of a previously developed method, rMCP.

A fine-scale examination of parturition timing in temperate ungulates.

Parturition timing has long been a topic of interest in ungulate research. However, few studies have examined parturition timing at fine scale (e.g., <1 day). Predator activity and environmental conditions can vary considerably with diel timing, which may result in selective pressure for parturition to occur during diel times that maximize the likelihood of neonate survival. We monitored parturition events and early-life survival of elk (Cervus canadensis) and mule deer (Odocoileus hemionus) in Utah, USA to better understand diel timing of parturition in temperate ungulates. Diel timing of parturition was moderately synchronous among conspecifics and influenced by environmental variables on the date of parturition. For elk, parturition events were most common during the morning crepuscular period and generally occurred later (i.e., closer to 12:00) when a relatively large proportion of the moon was illuminated. For mule deer, parturition events were most common during the diurnal period and generally occurred later (i.e., closer to 15:00) on cold, wet dates. Diel timing of parturition did not influence neonate survival, but larger datasets may be required to verify the apparent lack of influence. Although additional work could evaluate alternative variables that might affect parturition timing, our data provide an improved and finer scale understanding of reproductive ecology and phenology in ungulates.

Estimating age of mule deer in the field: Can we move beyond broad age categories?

Age of individuals is an intrinsic demographic parameter used in the modeling and management of wildlife. Although analysis of cementum annuli from teeth is currently the most accurate method used to age ungulates, the age of live ungulates in the field can be estimated by examining tooth wear and tooth replacement patterns. However, there may be limitations to aging based on tooth wear as the rate of tooth wear likely varies among individuals due to factors such as age, diet, environment, and sex. Our objective was to determine the reliability of estimating age for mule deer based on tooth wear and tooth replacement patterns. We compared ages estimated by tooth wear (collected at time of capture for a statewide monitoring effort) to ages determined from cementum analysis (from teeth collected after mortalities of radio-tracked animals from the monitoring effort). Accuracy was high; ages estimated from tooth wear were within one year of cementum ages >75% of the time when aged by experienced observers. Bias in accuracy for estimates of age was low but slightly biased toward underestimation (i.e., 0.6 years on average)-especially as cementum age increased. Our results indicate that aging mule deer using patterns in tooth wear can be reliable if observers estimating age have experience using this method.

Summer watering patterns of mule deer in the Great Basin Desert, USA: implications of differential use by individuals and the sexes for management of water resources.

Changes in the abundance and distribution of free water can negatively influence wildlife in arid regions. Free water is considered a limiting factor for mule deer (Odocoileus hemionus) in the Great Basin Desert. Consequently, a better understanding of differential use of water by individuals and the sexes could influence the conservation and management of mule deer and water resources in their habitats. We deployed remote cameras at all known water sources (13 wildlife water developments and 4 springs) on one mountain range in western Utah, USA, during summer from 2007 to 2011 to document frequency and timing of water use, number of water sources used by males and females, and to estimate population size from individually identified mule deer. Male and female mule deer used different water sources but visited that resource at similar frequencies. Individual mule deer used few water sources and exhibited high fidelity to that resource. Wildlife water developments were frequently used by both sexes. Our results highlight the differing use of water sources by sexes and individual mule deer. This information will help guide managers when siting and reprovisioning wildlife water developments meant to benefit mule deer and will contribute to the conservation and management of this species.

Variation in movement patterns of mule deer: have we oversimplified migration?

BACKGROUND: Conservation and management of migratory animals has gained attention in recent years, but the majority of research has focused on stereotypical 'migrant' and 'resident' behaviors, often failing to incorporate any atypical behaviors or characterize migratory behaviors beyond distance and timing of the migration. With migration threatened by anthropogenic development and climate change, it is crucial that we understand the full range of migratory behaviors. Our objective was to demonstrate and characterize the variation in migration strategies, including typical and atypical migratory behaviors for mule deer (Odocoileus hemionus) in Utah, USA. METHODS: Because calculation of common metrics such as distance, timing, and use of stopovers during migration did not adequately describe the variation we observed in migratory behavior for this species-particularly when animals visited multiple (> 3) ranges for extended lengths of time-we developed additional methods and categories to describe observed variation in migratory behavior. We first categorized trajectories based on the number of discrete, separate ranges and range shifts between them. Then, we further characterized the variation in migration strategies by examining the timing, duration, and distance traveled within each of the categories. We also examined if and how frequently individual deer switched among categories from year to year. RESULTS: We classified 1218 movement trajectories from 722 adult female mule deer, and found that 54.4% were dual-range migrants, who made one round-trip to one distinct range. Multi-range migrants (23.6%) made one round-trip during which they stayed at multiple discrete ranges. Commuters (1.0%) traveled to the same range multiple times, and poly migrants (1.5%) made multiple round-trips to different ranges. Gradual movers (2.5%) did not show a discrete range shift but moved gradually between ranges, whereas residents (12.6%) never left their home ranges, and dispersers (4.4%) left but never returned. Of the deer that we monitored for multiple years, 51.2% switched among categories. CONCLUSION: We conclude that the substantial number of atypical migratory strategies, as well as the number of deer that switched categories, underlines the importance of studying these less-stereotyped behaviors that may be exhibited by large proportions of populations. Acknowledging and investigating the full complexity and diversity in migratory strategies might uncover unknowns with respect to underlying factors and drivers of migration, and can help shape effective conservation strategies.

De novo chromosome-length assembly of the mule deer (Odocoileus hemionus) genome.

The mule deer (Odocoileus hemionus) is an ungulate species that is distributed in a range from western Canada to central Mexico. Mule deer are an essential source of food for many predators, are relatively abundant, and commonly make broad migration movements. A clearer understanding of the mule deer genome can improve our knowledge of its population genetics, movements, and demographic history, aiding in conservation efforts. Their large population size, continuous distribution, and diversity of habitat make mule deer excellent candidates for population genomics studies; however, few genomic resources are currently available for this species. Here, we sequence and assemble the mule deer genome into a highly contiguous chromosome-length assembly for use in future research using long-read sequencing and Hi-C technologies. We also provide a genome annotation and compare demographic histories of the mule deer and white-tailed deer using the pairwise sequentially Markovian coalescent model. We expect this assembly to be a valuable resource in the continued study and conservation of mule deer.

The influence of habitat use on harvest vulnerability of cow elk (Cervus canadensis).

Pressure from hunting can alter the behavior and habitat selection of game species. During hunting periods, cervids such as elk (Cervus canadensis) typically select for areas further from roads and closer to tree cover, while altering the timing of their daily activities to avoid hunters. Our objective was to determine the habitat characteristics most influential in predicting harvest risk of elk. We captured 373 female elk between January 2015 and March 2017 in the Uinta-Wasatch-Cache National Forest and surrounding area of central Utah, USA. We determined habitat selection during the hunting season using a resource selection function (RSF) for 255 adult cow elk. Additionally, we used a generalized linear mixed model to evaluate risk of harvest based on habitat use within home ranges (3rd order selection) as well as the location of the home range on the landscape to evaluate vulnerability on a broader scale. Female elk selected for areas that reduced hunter access (rugged terrain, within tree cover, on private land). Age, elevation and distance to roads within a home range were most influential in predicting harvest risk (top model accounted for 36.2% of AIC weight). Elevation and distance to trees were most influential in predicting risk when evaluating the location of the home range (top model accounted for 42.1% of AIC weight). Vulnerability to harvest was associated with proximity to roads. Additionally, survival in our landscape decreased with age of femaleelk.

Using satellite-derived estimates of plant phenological rhythms to predict sage-grouse nesting chronology.

The "green wave" hypothesis posits that during spring consumers track spatial gradients in emergent vegetation and associated foraging opportunities. This idea has largely been invoked to explain animal migration patterns, yet the general phenomenon underlies trends in vertebrate reproductive chronology as well. We evaluated the utility of this hypothesis for predicting spatial variation in nest initiation of greater sage-grouse (Centrocerus urophasianus), a species of conservation concern in western North America. We used the Normalized Difference Vegetation Index (NDVI) to map the green wave across elevation and then compiled dates and locations of >450 sage-grouse nests from 20 study sites (2000-2014) to model nest initiation as a function of the start of the growing season (SOS), defined here as the maximum daily rate of increase in NDVI. Individual sites were drawn from three ecoregions, distributed over 4.5° latitude, and spanning 2,300 m in elevation, which captured the climatic, edaphic, and floristic diversity of sagebrush ecosystems in the southern half of current sage-grouse range. As predicted, SOS displayed a significant, positive relationship with elevation, occurring 1.3 days later for each 100 m increase in elevation. In turn, sage-grouse nest initiation followed SOS by 22 ± 10 days (r2  = .57), with hatch dates falling on or just prior to the peak of the growing season. By timing nesting to the green wave, sage-grouse chicks hatched when the abundance of protein-rich invertebrate biomass is hypothesized to be nearing a seasonal high. This adaptation likely represents a strategy for maximizing reproductive success in the arid, variable environments that define sagebrush ecosystems. Given projected changes in climate and land use, these results can be used to predict periods of relative sensitivity to habitat disturbance for sage-grouse. Moreover, the near real-time availability of satellite imagery offers a heretofore underutilized means of mapping the green wave, planning habitat restoration, and monitoring range conditions.

Behavioral differences at scent stations between two exploited species of desert canids.

Coyotes (Canis latrans) and kit foxes (Vulpes macrotis) are desert canids that share ecological similarities, but have disparate histories with anthropogenic pressure that may influence their responses towards novel stimuli. We used remote cameras to investigate response to novel stimuli for these two species. We predicted that coyotes (heavily pressured species) would be more wary towards novel stimuli on unprotected land (canid harvest activities are permitted) than in protected areas (canid harvest activities are not permitted), whereas kit foxes (less pressured species) would exhibit no difference. We examined differences in the investigative behaviors at 660 scent stations in both protected and unprotected areas. Coyotes showed no differences between protected and unprotected land and were generally more wary than kit foxes, supporting our prediction. Kit foxes were more investigative on protected land, contrary to our expectations. Our study provides evidence that anthropogenic pressure can alter the behaviors of wildlife species.

Correction: The effects of electric power lines on the breeding ecology of greater sage-grouse.

[This corrects the article DOI: 10.1371/journal.pone.0209968.].

The effects of electric power lines on the breeding ecology of greater sage-grouse.

Anthropogenic infrastructure can negatively affect wildlife through direct mortality and/or displacement behaviors. Some tetranoids (grouse spp.) species are particularly vulnerable to tall anthropogenic structures because they evolved in ecosystems void of vertical structures. In western North America, electric power transmission and distribution lines (power lines) occur in sagebrush (Artemisia spp.) landscapes within the range of the greater sage-grouse (Centrocercus urophasianus; sage-grouse). The U.S. Fish and Wildlife Service recommended using buffer zones near leks to mitigate the potential impacts of power lines on sage-grouse. However, recommended buffer distances are inconsistent across state and federal agencies because data are lacking. To address this, we evaluated the effects of power lines on sage-grouse breeding ecology within Utah, portions of southeastern Idaho, and southwestern Wyoming from 1998-2013. Overall, power lines negatively affected lek trends up to a distance of 2.7 and 2.8 km, respectively. Power lines died not affect lek persistence. Female sage-grouse avoided transmission lines during the nesting and brooding seasons at distances up to 1.1 and 0.8 km, respectively. Nest and brood success were negatively affected by transmission lines up to distances of 2.6 and 1.1 km, respectively. Distribution lines did not appear to affect sage-grouse habitat selection or reproductive fitness. Our analyses demonstrated the value of sagebrush cover in mitigating potential power line impacts. Managers can minimize the effects of new transmission power lines by placing them in existing anthropogenic corridors and/or incorporating buffers at least 2.8 km from active leks. Given the uncertainty we observed in our analyses regarding sage-grouse response to distribution lines coupled with their role in providing electric power service directly to individual consumers, we recommend that buffers for these power lines be considered on a case-by-case basis. Micrositing to avoid important habitats and habitat reclamation may reduce the potential impacts of new power line construction.

Sharp-Tailed Grouse Nest Survival and Nest Predator Habitat Use in North Dakota's Bakken Oil Field.

Recent advancements in extraction technologies have resulted in rapid increases of gas and oil development across the United States and specifically in western North Dakota. This expansion of energy development has unknown influences on local wildlife populations and the ecological interactions within and among species. Our objectives for this study were to evaluate nest success and nest predator dynamics of sharp-tailed grouse (Tympanuchus phasianellus) in two study sites that represented areas of high and low energy development intensities in North Dakota. During the summers of 2012 and 2013, we monitored 163 grouse nests using radio telemetry. Of these, 90 nests also were monitored using miniature cameras to accurately determine nest fates and identify nest predators. We simultaneously conducted predator surveys using camera scent stations and occupancy modeling to estimate nest predator occurrence at each site. American badgers (Taxidea taxus) and striped skunks (Mephitis mephitis) were the primary nest predators, accounting for 56.7% of all video recorded nest depredations. Nests in our high intensity gas and oil area were 1.95 times more likely to succeed compared to our minimal intensity area. Camera monitored nests were 2.03 times more likely to succeed than non-camera monitored nests. Occupancy of mammalian nest predators was 6.9 times more likely in our study area of minimal gas and oil intensity compared to the high intensity area. Although only a correlative study, our results suggest energy development may alter the predator community, thereby increasing nest success for sharp-tailed grouse in areas of intense development, while adjacent areas may have increased predator occurrence and reduced nest success. Our study illustrates the potential influences of energy development on the nest predator-prey dynamics of sharp-tailed grouse in western North Dakota and the complexity of evaluating such impacts on wildlife.

Negative effects of an exotic grass invasion on small-mammal communities.

Exotic invasive species can directly and indirectly influence natural ecological communities. Cheatgrass (Bromus tectorum) is non-native to the western United States and has invaded large areas of the Great Basin. Changes to the structure and composition of plant communities invaded by cheatgrass likely have effects at higher trophic levels. As a keystone guild in North American deserts, granivorous small mammals drive and maintain plant diversity. Our objective was to assess potential effects of invasion by cheatgrass on small-mammal communities. We sampled small-mammal and plant communities at 70 sites (Great Basin, Utah). We assessed abundance and diversity of the small-mammal community, diversity of the plant community, and the percentage of cheatgrass cover and shrub species. Abundance and diversity of the small-mammal community decreased with increasing abundance of cheatgrass. Similarly, cover of cheatgrass remained a significant predictor of small-mammal abundance even after accounting for the loss of the shrub layer and plant diversity, suggesting that there are direct and indirect effects of cheatgrass. The change in the small-mammal communities associated with invasion of cheatgrass likely has effects through higher and lower trophic levels and has the potential to cause major changes in ecosystem structure and function.

Long-lasting effects of maternal condition in free-ranging cervids.

Causes of phenotypic variation are fundamental to evolutionary ecology because they influence the traits acted upon by natural selection. One such cause of phenotypic variation is a maternal effect, which is the influence of the environment experienced by a female (and her corresponding phenotype) on the phenotype of her offspring (independent of the offspring's genotype). While maternal effects are well documented, the longevity and fitness impact of these effects remains unclear because it is difficult to follow free-living individuals through their reproductive lifetimes. For long-lived species, it has been suggested that maternal effects are masked by environmental variables acting on offspring in years following the period of dependence. Our objective was to use indirect measures of maternal condition to determine if maternal effects have long-lasting influences on male offspring in two species of cervid. Because antlers are sexually selected, we used measures of antler size at time of death, 1.5-21.5 years after gestation to investigate maternal effects. We quantified antler size of 11,000 male elk and mule deer born throughout the intermountain western US (6 states) over nearly 30 years. Maternal condition during development was estimated indirectly using a suite of abiotic variables known to influence condition of cervids (i.e., winter severity, spring and summer temperature, and spring and summer precipitation). Antler size of male cervids was significantly associated with our indirect measure of maternal condition during gestation and lactation. Assuming the correctness of our indirect measure, our findings demonstrate that antler size is a sexually selected trait that is influenced-into adulthood-by maternal condition. This link emphasizes the importance of considering inherited environmental effects when interpreting population dynamics or examining reproductive success of long-lived organisms.

Water developments and canids in two North American deserts: a test of the indirect effect of water hypothesis.

Anthropogenic modifications to landscapes intended to benefit wildlife may negatively influence wildlife communities. Anthropogenic provisioning of free water (water developments) to enhance abundance and distribution of wildlife is a common management practice in arid regions where water is limiting. Despite the long-term and widespread use of water developments, little is known about how they influence native species. Water developments may negatively influence arid-adapted species (e.g., kit fox, Vulpes macrotis) by enabling water-dependent competitors (e.g., coyote, Canis latrans) to expand distribution in arid landscapes (i.e., indirect effect of water hypothesis). We tested the two predictions of the indirect effect of water hypothesis (i.e., coyotes will visit areas with free water more frequently and kit foxes will spatially and temporally avoid coyotes) and evaluated relative use of free water by canids in the Great Basin and Mojave Deserts from 2010 to 2012. We established scent stations in areas with (wet) and without (dry) free water and monitored visitation by canids to these sites and visitation to water sources using infrared-triggered cameras. There was no difference in the proportions of visits to scent stations in wet or dry areas by coyotes or kit foxes at either study area. We did not detect spatial (no negative correlation between visits to scent stations) or temporal (no difference between times when stations were visited) segregation between coyotes and kit foxes. Visitation to water sources was not different for coyotes between study areas, but kit foxes visited water sources more in Mojave than Great Basin. Our results did not support the indirect effect of water hypothesis in the Great Basin or Mojave Deserts for these two canids.

Vigilance of kit foxes at water sources: a test of competing hypotheses for a solitary carnivore subject to predation.

Animals that are potential prey do not respond equally to direct and indirect cues related to risk of predation. Based on differential responses to cues, three hypotheses have been proposed to explain spatial variation in vigilance behavior. The predator-vigilance hypothesis proposes that prey increase vigilance where there is evidence of predators. The visibility-vigilance hypothesis suggests that prey increase vigilance where visibility is obstructed. Alternatively, the refuge-vigilance hypothesis proposes that prey may perceive areas with low visibility (greater cover) as refuges and decrease vigilance. We evaluated support for these hypotheses using the kit fox (Vulpes macrotis), a solitary carnivore subject to intraguild predation, as a model. From 2010 to 2012, we used infrared-triggered cameras to record video of kit fox behavior at water sources in Utah, USA. The refuge-vigilance hypothesis explained more variation in vigilance behavior of kit foxes than the other two hypotheses (AICc model weight=0.37). Kit foxes were less vigilant at water sources with low overhead cover (refuge) obstructing visibility. Based on our results, the predator-vigilance and visibility-vigilance hypotheses may not be applicable to all species of prey. Solitary prey, unlike gregarious prey, may use areas with concealing cover to maximize resource acquisition and minimize vigilance.