Maintaining functional connectivity is critical for the long-term conservation of wildlife populations. Landscape genomics provides an opportunity to assess long-term functional connectivity by relating environmental variables to spatial patterns of genomic variation resulting from generations of movement, dispersal and mating behaviors. Identifying landscape features associated with gene flow at large geographic scales for highly mobile species is becoming increasingly possible due to more accessible genomic approaches, improved analytical methods and enhanced computational power. We characterized the genetic structure and diversity of migratory mule deer Odocoileus hemionus using 4051 single nucleotide polymorphisms in 406 individuals sampled across multiple habitats throughout Wyoming, USA. We then identified environmental variables associated with genomic variation within genetic groups and statewide using a stepwise approach to first evaluate nonlinear relationships of landscape resistance with genetic distances and then use mixed-effects modeling to choose top landscape genomic models. We identified three admixed genetic groups of mule deer and found that environmental variables associated with gene flow varied among genetic groups, revealing scale-dependent and regional variation in functional connectivity. At the statewide scale, more gene flow occurred in areas with low elevation and mixed habitat. In the southern genetic group, more gene flow occurred in areas with low elevation. In the northern genetic group, more gene flow occurred in grassland and forest habitats, while highways and energy infrastructure reduced gene flow. In the western genetic group, the null model of isolation by distance best represented genetic patterns. Overall, our findings highlight the role of different seasonal ranges on mule deer genetic connectivity, and show that anthropogenic features hinder connectivity. This study demonstrates the value of combining a large, genome-wide marker set with recent advances in landscape genomics to evaluate functional connectivity in a wide-ranging migratory species.

The green wave hypothesis (GWH) states that migrating animals should track or 'surf' high-quality forage at the leading edge of spring green-up. To index such high-quality forage, recent work proposed the instantaneous rate of green-up (IRG), i.e. rate of change in the normalized difference vegetation index over time. Despite this important advancement, no study has tested the assumption that herbivores select habitat patches at peak IRG. We evaluated this assumption using step selection functions parametrized with movement data during the green-up period from two populations each of bighorn sheep, mule deer, elk, moose and bison, totalling 463 individuals monitored 1-3 years from 2004 to 2014. Accounting for variables that typically influence habitat selection for each species, we found seven of 10 populations selected patches exhibiting high IRG-supporting the GWH. Nonetheless, large herbivores selected for the leading edge, trailing edge and crest of the IRG wave, indicating that other mechanisms (e.g. ruminant physiology) or measurement error inherent with satellite data affect selection for IRG. Our evaluation indicates that IRG is a useful tool for linking herbivore movement with plant phenology, paving the way for significant advancements in understanding how animals track resource quality that varies both spatially and temporally.

Accurately quantifying species' area requirements is a prerequisite for effective area-based conservation. This typically involves collecting tracking data on species of interest and then conducting home-range analyses. Problematically, autocorrelation in tracking data can result in space needs being severely underestimated. Based on the previous work, we hypothesized the magnitude of underestimation varies with body mass, a relationship that could have serious conservation implications. To evaluate this hypothesis for terrestrial mammals, we estimated home-range areas with global positioning system (GPS) locations from 757 individuals across 61 globally distributed mammalian species with body masses ranging from 0.4 to 4000 kg. We then applied block cross-validation to quantify bias in empirical home-range estimates. Area requirements of mammals <10 kg were underestimated by a mean approximately15%, and species weighing approximately100 kg were underestimated by approximately50% on average. Thus, we found area estimation was subject to autocorrelation-induced bias that was worse for large species. Combined with the fact that extinction risk increases as body mass increases, the allometric scaling of bias we observed suggests the most threatened species are also likely to be those with the least accurate home-range estimates. As a correction, we tested whether data thinning or autocorrelation-informed home-range estimation minimized the scaling effect of autocorrelation on area estimates. Data thinning required an approximately93% data loss to achieve statistical independence with 95% confidence and was, therefore, not a viable solution. In contrast, autocorrelation-informed home-range estimation resulted in consistently accurate estimates irrespective of mass. When relating body mass to home range size, we detected that correcting for autocorrelation resulted in a scaling exponent significantly >1, meaning the scaling of the relationship changed substantially at the upper end of the mass spectrum.