When can we measure stress noninvasively? Postdeposition effects on a fecal stress metric confound a multiregional assessment

Measurement of stress hormone metabolites in fecal samples has become a common method to assess physiological stress in wildlife populations. Glucocorticoid metabolite (GCM) measurements can be collected noninvasively, and studies relating this stress metric to anthropogenic disturbance are increasing. However, environmental characteristics (e.g., temperature) can alter measured GCM concentration when fecal samples cannot be collected immediately after defecation. This effect can confound efforts to separate environmental factors causing predeposition physiological stress in an individual from those acting on a fecal sample postdeposition. We used fecal samples from American pikas (Ochotona princeps) to examine the influence of environmental conditions on GCM concentration by (1) comparing GCM concentration measured in freshly collected control samples to those placed in natural habitats for timed exposure, and (2) relating GCM concentration in samples collected noninvasively throughout the western United States to local environmental characteristics measured before and after deposition. Our timed‐exposure trials clarified the spatial scale at which exposure to environmental factors postdeposition influences GCM concentration in pika feces. Also, fecal samples collected from occupied pika habitats throughout the species' range revealed significant relationships between GCM and metrics of climate during the postdeposition period (maximum temperature, minimum temperature, and precipitation during the month of sample collection). Conversely, we found no such relationships between GCM and metrics of climate during the predeposition period (prior to the month of sample collection). Together, these results indicate that noninvasive measurement of physiological stress in pikas across the western US may be confounded by climatic conditions in the postdeposition environment when samples cannot be collected immediately after defecation. Our results reiterate the importance of considering postdeposition environmental influences on this stress metric, especially in multiregional comparisons. However, measurements of fecal GCM concentration should prove useful for population monitoring within an eco‐region or when postdeposition exposure can be minimized.     

Individual‐based analysis of hair corticosterone reveals factors influencing chronic stress in the American pika

Glucocorticoids are often measured in wildlife to assess physiological responses to environmental or ecological stress. Hair, blood, saliva, or fecal samples are generally used depending on the timescale of the stress response being investigated and species‐specific considerations. Here, we report the first use of hair samples to measure long‐term corticosterone levels in the climate‐sensitive American pika (Ochotona princeps). We validated an immunoassay‐based measurement of corticosterone extracted from hair samples and compared corticosterone estimates obtained from plasma, hair, and fecal samples of nine pikas. To demonstrate an ecological application of this technique, we characterized physiological stress in 49 pikas sampled and released at eight sites along two elevational transects. Microclimate variation was measured at each site using both ambient and subsurface temperature sensors. We used an information theoretic approach to compare support for linear, mixed‐effects models relating corticosterone estimates to microclimate, body size, and sex. Corticosterone was measured accurately in pika hair samples after correcting for the influence of sample mass on corticosterone extraction efficiency. Hair‐ and plasma‐based estimates of corticosterone were weakly correlated. The best‐supported model suggested that corticosterone was lower in larger, male pikas, and at locations with higher ambient temperatures in summer. Our results are consistent with a general negative relationship between body mass and glucocorticoid concentration observed across mammalian species, attributed to the higher mass‐specific metabolic rates of smaller bodied animals. The higher corticosterone levels in female pikas likely reflected the physiological demands of reproduction, as observed in a wide array of mammalian species. Additionally, we establish the first direct physiological evidence for thermal stress in the American pika through nonlethal sampling of corticosterone. Interestingly, our data suggest evidence for cold stress likely induced during the summer molting period. This technique should provide a useful tool to researchers wishing to assess chronic stress in climate‐sensitive mammals.     

Landscape effects on gene flow for a climate‐sensitive montane species,the American pika

Climate change is arguably the greatest challenge to conservation of our time. Most vulnerability assessments rely on past and current species distributions to predict future persistence but ignore species' abilities to disperse through landscapes, which may be particularly important in fragmented habitats and crucial for long‐term persistence in changing environments. Landscape genetic approaches explore the interactions between landscape features and gene flow and can clarify how organisms move among suitable habitats, but have suffered from methodological uncertainties. We used a landscape genetic approach to determine how landscape and climate‐related features influence gene flow for American pikas (Ochotona princeps) in Crater Lake National Park. Pikas are heat intolerant and restricted to cool microclimates; thus, range contractions have been predicted as climate changes. We evaluated the correlation between landscape variables and genetic distance using partial Mantel tests in a causal modelling framework, and used spatially explicit simulations to evaluate methods of model optimization including a novel approach based on relative support and reciprocal causal modelling. We found that gene flow was primarily restricted by topographic relief, water and west‐facing aspects, suggesting that physical restrictions related to small body size and mode of locomotion, as well as exposure to relatively high temperatures, limit pika dispersal in this alpine habitat. Our model optimization successfully identified landscape features influencing resistance in the simulated data for this landscape, but underestimated the magnitude of resistance. This is the first landscape genetic study to address the fundamental question of what limits dispersal and gene flow in the American pika.     

Variation in subsurface thermal characteristics of microrefuges used by range core and peripheral populations of the American pika (Ochotona princeps)

Microrefuges provide microclimates decoupled from inhospitable regional climate regimes that enable range‐peripheral populations to persist and are important to cold‐adapted species in an era of accelerated climate change. However, identifying and describing the thermal characteristics of microrefuge habitats is challenging, particularly for mobile organisms in cryptic, patchy habitats. We examined variation in subsurface thermal conditions of microrefuge habitats among different rock substrate types used by the American pika (Ochotona princeps ), a climate‐sensitive, rock‐dwelling Lagomorph. We compared subsurface temperatures in talus and lava substrates in pika survey sites in two US national park units; one park study area on the range periphery and the other in the range core. We deployed paired sensors to examine within‐site temperature variation. We hypothesized that subsurface temperatures within occupied sites and structurally complex substrates would be cooler in summer and warmer in winter than unoccupied and less complex sites. Although within‐site variability was high, with correlations between paired sensors as low as 47%, we found compelling evidence that pikas occupy microrefuge habitats where subsurface conditions provide more thermal stability than in unoccupied microhabitats. The percentage of days in which microhabitat temperatures were between ?2.5 and 25.5°C was significantly higher in occupied sites. Interestingly, thermal conditions were substantially more stable (p  < .05) in the lava substrate type identified to be preferentially used by pikas (pahoehoe vs. a'a) in a previous study. Our study and others suggest that thermal stability appears to be the defining characteristic of subsurface microrefuges used by American pikas and is a likely explanation for enigmatic population persistence at the range periphery. Our study exemplifies an integrated approach for studying complex microhabitat conditions, paired with site use surveys and contextualized with information about gene flow provided by complementary studies.     

Mechanistic variables can enhance predictive models of endotherm distributions: the American pika under current,past, and future climates

How climate constrains species’ distributions through time and space is an important question in the context of conservation planning for climate change. Despite increasing awareness of the need to incorporate mechanism into species distribution models (SDMs), mechanistic modeling of endotherm distributions remains limited in this literature. Using the American pika (Ochotona princeps) as an example, we present a framework whereby mechanism can be incorporated into endotherm SDMs. Pika distribution has repeatedly been found to be constrained by warm temperatures, so we used Niche Mapper, a mechanistic heat‐balance model, to convert macroclimate data to pika‐specific surface activity time in summer across the western United States. We then explored the difference between using a macroclimate predictor (summer temperature) and using a mechanistic predictor (predicted surface activity time) in SDMs. Both approaches accurately predicted pika presences in current and past climate regimes. However, the activity models predicted 8–19% less habitat loss in response to annual temperature increases of ~3–5 °C predicted in the region by 2070, suggesting that pikas may be able to buffer some climate change effects through behavioral thermoregulation that can be captured by mechanistic modeling. Incorporating mechanism added value to the modeling by providing increased confidence in areas where different modeling approaches agreed and providing a range of outcomes in areas of disagreement. It also provided a more proximate variable relating animal distribution to climate, allowing investigations into how unique habitat characteristics and intraspecific phenotypic variation may allow pikas to exist in areas outside those predicted by generic SDMs. Only a small number of easily obtainable data are required to parameterize this mechanistic model for any endotherm, and its use can improve SDM predictions by explicitly modeling a widely applicable direct physiological effect: climate‐imposed restrictions on activity. This more complete understanding is necessary to inform climate adaptation actions, management strategies, and conservation plans.     

达乌尔鼠兔扩散过程中的生境选择

通过调控牧压,研究了达乌尔鼠兔扩散过程中的生境选择。结果表明,放牧活动强了草原植被的斑块性,改变了植被为鼠兔提供的资源条件,导致鼠兔斑块间的扩散和扩散鼠兔对不同植被斑块的生境选择。     

Habitat availability and gene flow influence diverging local population trajectories under scenarios of climate change: a place‐based approach

Ecological niche theory holds that species distributions are shaped by a large and complex suite of interacting factors. Species distribution models (SDMs) are increasingly used to describe species’ niches and predict the effects of future environmental change, including climate change. Currently, SDMs often fail to capture the complexity of species’ niches, resulting in predictions that are generally limited to climate‐occupancy interactions. Here, we explore the potential impact of climate change on the American pika using a replicated place‐based approach that incorporates climate, gene flow, habitat configuration, and microhabitat complexity into SDMs. Using contemporary presence–absence data from occupancy surveys, genetic data to infer connectivity between habitat patches, and 21 environmental niche variables, we built separate SDMs for pika populations inhabiting eight US National Park Service units representing the habitat and climatic breadth of the species across the western United States. We then predicted occurrence probability under current (1981–2010) and three future time periods (out to 2100). Occurrence probabilities and the relative importance of predictor variables varied widely among study areas, revealing important local‐scale differences in the realized niche of the American pika. This variation resulted in diverse and – in some cases – highly divergent future potential occupancy patterns for pikas, ranging from complete extirpation in some study areas to stable occupancy patterns in others. Habitat composition and connectivity, which are rarely incorporated in SDM projections, were influential in predicting pika occupancy in all study areas and frequently outranked climate variables. Our findings illustrate the importance of a place‐based approach to species distribution modeling that includes fine‐scale factors when assessing current and future climate impacts on species’ distributions, especially when predictions are intended to manage and conserve species of concern within individual protected areas.     

Are pikas exposed to and affected by selenium deficiency?

Regional extirpations of pikas (Ochatona princeps) within the last few decades have been attributed to global warming. Other recent global alterations such as increased nitrogen (N) deposition and associated selenium (Se) deficiency may further stress pika populations. In 2003 and 2004, we live-trapped pikas from three populations in Wyoming and measured Se values in their hair. We also sampled hair and liver from museum specimens collected throughout the Northern Rocky Mountains in 1987 and 1988. Our results suggest that liver and hair values were related, and that pika hair reflected the Se concentrations of the geologic parent materials. We determined that animals residing in several remote areas in the Rocky Mountain region could be Se deficient and that increase in N deposition correlated with an increase rather than a decrease in Se values in pika hair. In addition, we found no relation between Se contents in hair and body condition index, suggesting that low Se levels may not have negative effects on individual pikas. Whether Se levels influence reproductive success of pikas is unknown and should be the focus of future studies.     

高原鼠兔血清无机元素含量的分析

测定分析了高原鼠兔血清电解质、无机元素、血清铁和总铁结合力, 并与人和其它实验动物的相应指标进行了比较。得出以下主要结论: 1.高原鼠兔血清电解质浓度不受生活环境和食物条件的影响, 其机体能自身调节血清离子浓度, 以维持电解质在血液中的相对衡定; 2.高原鼠兔血清铁含量与其栖息地海拔相关, 随海拔升高血清铁含量增加; 3.高原鼠兔通过增加运铁蛋白与铁的结合力作为适应高原低氧环境的方式之一。     

A diet rich in C3 plants reveals the sensitivity of an alpine mammal to climate change

Plant–herbivore interactions provide critical insights into the mechanisms that govern the spatiotemporal distributions of organisms. These interactions are crucial to understanding the impacts of climate change, which are likely to have an effect on the population dynamics of alpine herbivores. The Royle's pika (Ochotona roylei, hereafter pika) is a lagomorph found in the western Himalaya and is dependent on alpine plants that are at risk from climate change. As the main prey of many carnivores in the region, the pika plays a crucial role in trophic interactions. We examined topographical features, plant genera presence and seasonal dynamics as drivers of the plant richness in the pika's diet across an elevational gradient (2,600–4,450 m). We identified 79 plant genera in the faecal pellets of pikas, of which 89% were forbs, >60% were endemic to the Himalaya, and 97.5% of the diet plant genera identified followed the C₃ photosynthetic pathway. We found that, during the premonsoon season, the number of genera in the pika's diet decreased with increasing elevation. We demonstrate that a large area of talus supports greater plant diversity and, not surprisingly, results in higher species richness in the pika's diet. However, in talus habitat with deep crevices, pikas consumed fewer plant genera suggesting they may be foraging suboptimally due to predation risk. The continued increase in global temperature is expected to have an effect on the distribution dynamics of C₃ plants and consequently influence pika diet and distribution, resulting in a significant negative cascading effect on the Himalayan ecosystem.     

Influences of chronic and current season grazing by collared pikas on above-ground biomass and species richness in subarctic alpine meadows

We studied an alpine herbivory gradient established by collared pikas, a small central place foraging lagomorph, to examine the effects of multiple grazing levels on above-ground live biomass (AGB) and species richness (SR) in alpine tundra. The effects of within-season (four sampling periods), multi-season (across three summers) and longer-term dynamics (inferred from spatial location of vegetation with respect to pika haypiles) were examined. Along the grazing gradient, we found support for and against hypotheses that propose biphasic, increasing, or decreasing responses to herbivory, both in terms of AGB and SR. Our results suggest that plant–herbivore predictability is still weak. To further examine the impact of herbivory, we experimentally removed pikas using mesh exclosures placed at increasing distance from the edge of talus occupied by pikas. AGB after the second consecutive year of herbivore exclusion increased by 125% compared to control plots in highly grazed areas adjacent to talus (<1 m). In more lightly grazed sites at distances 1–6 m from talus, AGB increased by more than 40% after pikas were removed. No differences were observed in the ungrazed sites >6 m from talus. AGB was highest in meadow patches previously grazed by pikas compared to those with little grazing history, but this response was only observed after two seasons following release from herbivory. Grazed sites at distances of 1–6 m had the highest SR. These results indicate that multi-year measurements of growth are particularly relevant in ecosystems dominated by long-lived perennials in regions where productivity is low. Infrequent herbivore vacancies may provide local short-term release from pika grazing, thereby contributing to the persistence of productive, highly palatable vegetation.     

高原鼠兔对垂穗披碱草及豆科植物生长特性的影响

高原鼠兔通过取食对高寒草甸产生广泛的影响,从而影响高寒草甸植物的个体生长特征。2012年在青海省果洛州矮嵩草草甸,以高原鼠兔喜食植物(豆科和垂穗披碱草)为研究对象,通过小区控制研究高原鼠兔取食（对照组、低密度处理组、中密度处理组、高密度处理组）对植物个体生长特征的影响。研究显示：在实验前期（5月和6月）,喜食植物个体特征基本上无明显规律;在实验后期（7月和8月）,随着高原鼠兔取食强度增大,高原鼠兔喜食植物的地上生物量、平均高度、分蘖枝数、分蘖枝长度、丛幅均呈下降趋势;叶片长度和叶片数均呈先上升后下降态势。2.高原鼠兔喜食植物个体特征的增长率对照组显著高于高密度处理组（P<0.05）。以上结果表明(1)高寒草甸豆类植物在其生长后期与高原鼠兔的取食强度呈一定的负相关关系;高原鼠兔取食强度在阈值（中、低密度）内有利于垂穗披碱草的生长,但超过其阈值的取食强度不利于其生长。(2)高原鼠兔高强度的取食对其喜食植物个体的生长率有明显的抑制作用。     

Diversity of diazotrophic gut inhabitants of pikas (Ochotonidae) revealed by PCR-DGGE analysis

Diazotrophic gut symbionts are considered to act as nitrogen providers for their hosts, as was shown for various termite species. Although the diet of lagomorphs, like pikas or rabbits, is very poor in nitrogen and energy, their fecal matter contains 30–40% of protein. Since our hypothesis was that pikas maintained a diazotrophic consortium in their gastrointestinal tract, we conducted the first investigation of microbial diversity in pika guts. We obtained gut samples from animals of several Ochotona species, O. hyperborea (Northern pika), O. mantchurica (Manchurian pika), and O. dauurica (Daurian pika), in order to retrieve and compare the nitrogen-fixing communities of different pika species. The age and gender of the animals were taken into consideration. We amplified 320-bp long fragments of the nifH gene using the DNA extracted directly from the colon and cecum samples of pika’s gut, resolved them by DGGE, and performed phylogenetic reconstruction of 51 sequences obtained from excised bands. No significant difference was detected between the nitrogen-fixing gut inhabitants of different pika species. NifH sequences fell into two clusters. The first cluster contained the sequences affiliated with NifH Cluster I (Zehr et al., 2003) with similarity to Sphingomonas sp., Bradyrhizobium sp., and various uncultured bacteria from soil and rhizosphere. Sequences from the second group were related to Treponema sp., Fibrobacter succinogenes, and uncultured clones from the guts of various termites and belonged to NifH Cluster III. We suggest that diazotrophic organisms from the second cluster are genuine endosymbionts of pikas and provide nitrogen for further synthesis processes thus allowing these animals not to be short of protein.     

Keystone status of plateau pikas (Ochotona curzoniae): effect of control on biodiversity of native birds

The plateau pika (Ochotona curzoniae) of the Qinghai–Xizang (Tibetan) plateau, People's Republic of China, has been considered a pest because it putatively competes with native livestock for forage and contributes to rangeland degradation. As a result the plateau pika has been poisoned across vast areas of the high alpine meadows of the plateau. The plateau pika has also been considered a keystone species for biodiversity on the plateau. As one test of the keystone species hypothesis, we investigated the effects of poisoning plateau pikas on avian species richness and abundance. We conducted standardized censuses of birds on a number of sites across the alpine grassland of Qinghai province on which pikas either had or had not been poisoned. Avian species richness and abundance were higher on non-poisoned sites, in particular for species that nest in pika burrows such as Hume's groundpecker (Pseudopodoces humilis) and six species of snowfinch (Montifringilla spp., Pyrgilauda spp.), and species that prey on pikas (upland buzzard, Buteo hemilasius; black-eared kite, Milvus lineatus). The plateau pika thus appears to be both an allogenic engineer and a keystone species. Poisoning pikas reduces biodiversity of native species on the Qinghai–Xizang plateau, therefore management decisions concerning plateau pikas should reflect caution and careful assessment.     

高原鼠兔在不同植被覆盖的立地占据、局部灭绝和建群率研究

青藏高原的草地生态系统被野生动物和家畜广泛地利用。大型有蹄类动物的放牧和小型哺乳类动物的挖掘行为,会改变其栖息地的植物群落。基于拉萨市墨竹工卡县日多镇年村牦牛(Bos grunniens)放牧的实践,我们发现高原鼠兔(Ochotona curzoniae)立地占域率、建群和区域灭绝率与栖息地变量之间存在一定的格局关系。为了了解高原鼠兔立地占域率、建群和区域灭绝率与栖息地之间的关系,我们采用多季节性立地占据模型,利用Presence软件对相关数据进行分析。结果表明,高原鼠兔局部灭绝率与植被覆盖度之间存在正相关性,而与建群率之间则呈负相关性。     

Noninvasive sampling reveals population genetic structure in the Royle’s pika,Ochotona roylei,in the western Himalaya

Understanding population genetic structure of climate‐sensitive herbivore species is important as it provides useful insights on how shifts in environmental conditions can alter their distribution and abundance. Herbivore responses to the environment can have a strong indirect cascading effect on community structure. This is particularly important for Royle's pika (Lagomorpha: Ochotona roylei), a herbivorous talus‐dwelling species in alpine ecosystem, which forms a major prey base for many carnivores in the Himalayan arc. In this study, we used seven polymorphic microsatellite loci to detect evidence for recent changes in genetic diversity and population structure in Royle's pika across five locations sampled between 8 and 160 km apart in the western Himalaya. Using four clustering approaches, we found the presence of significant contemporary genetic structure in Royle's pika populations. The detected genetic structure could be primarily attributed to the landscape features in alpine habitat (e.g., wide lowland valleys, rivers) that may act as semipermeable barriers to gene flow and distribution of food plants, which are key determinants in spatial distribution of herbivores. Pika showed low inbreeding coefficients (F_IS) and a high level of pairwise relatedness for individuals within 1 km suggesting low dispersal abilities of talus‐dwelling pikas. We have found evidence of a recent population bottleneck, possibly due to effects of environmental disturbances (e.g., snow melting patterns or thermal stress). Our results reveal significant evidence of isolation by distance in genetic differentiation (F_ST range = 0.04–0.19). This is the first population genetics study on Royle's pika, which helps to address evolutionary consequences of climate change which are expected to significantly affect the distribution and population dynamics in this talus‐dwelling species.     

The Importance of Biologically Relevant Microclimates in Habitat Suitability Assessments

Predicting habitat suitability under climate change is vital to conserving biodiversity. However, current species distribution models rely on coarse scale climate data, whereas fine scale microclimate data may be necessary to assess habitat suitability and generate predictive models. Here, we evaluate disparities between temperature data at the coarse scale from weather stations versus fine-scale data measured in microhabitats required for a climate-sensitive mammal, the American pika (Ochotona princeps). We collected two years of temperature data in occupied talus habitats predicted to be suitable (high elevation) and unsuitable (low elevation) by the bioclimatic envelope approach. At low elevations, talus surface and interstitial microclimates drastically differed from ambient temperatures measured on-site and at a nearby weather station. Interstitial talus temperatures were frequently decoupled from high ambient temperatures, resulting in instantaneous disparities of over 30°C between these two measurements. Microhabitat temperatures were also highly heterogeneous, such that temperature measurements within the same patch of talus were not more correlated than measurements at distant patches. An experimental manipulation revealed that vegetation cover may cool the talus surface by up to 10°C during the summer, which may contribute to this spatial heterogeneity. Finally, low elevation microclimates were milder and less variable than typical alpine habitat, suggesting that, counter to species distribution model predictions, these seemingly unsuitable habitats may actually be better refugia for this species under climate change. These results highlight the importance of fine-scale microhabitat data in habitat assessments and underscore the notion that some critical refugia may be counterintuitive.     

Social and ecological influences on dispersal and philopatry in the plateau pika (Ochotona curzoniae)

Benefits and costs of dispersal and philopatry of the socialplateau pika (Ochotona curzoniae) were studied on the Tibetanplateau for 3 years. Although short-lived, plateau pikas livein cohesive family groups that occupy burrow systems in sedgemeadow habitat Most (57.8%) plateau pikas were philopatric,and dispersal movements were extremely restricted. No juvenilefemales or adult pikas moved more than two family ranges betweenyears; the greatest observed dispersal distances were by twojuvenile males that moved five family ranges from the familyof their birth. Traversing unfamiliar habitat was not a costof pika dispersal because most dispersers settled in familiesthat they could easily visit before dispersal. Dispersal movementsappeared to result in equalization of density among pika families,an expected result if competition for environmental resourcesinfluenced dispersal. Males did not disperse to gain advantagesin competition for mates, as evidenced by their moving to familieswith significantly fewer females. Females, however, moved tofamilies with significantly more males. Males provide abundantpaternal care, and significantly more offspring per female survivedto become adults from families with more adult males per adultfemale. Evidence concerning the influence of inbreeding avoidanceon natal dispersal was indirect. Some males exhibited natalphilopatry; thus some families had opportunity for dose inbreeding.Males and females that dispersed had no opposite-sex relativesin their new families. Philopatric pikas may have benefitedby remaining in families that exhibited low local densities,and philopatric females might have benefited from social cooperationwith relatives.     

Genetic variation and fine-scale population structure in American pikas across a human-modified landscape

Natural resource extraction can represent a major human modification to the landscape. Habitat reclamation is becoming an increasingly important strategy for abating the loss of biodiversity associated with these developments; however, the demographic and genetic consequences of colonizing artificial habitat remain unknown in many species. Here, we investigated the genetic consequences of landscape modifications for the American pika (Ochotona princeps) relative to two major developments in British Columbia, Canada: a large open-pit copper mine (Highland Valley Copper) under partial reclamation and a bisecting major highway (97C). We assessed microsatellite genotypic data for 109 individuals across 15 sites located either within the mine on artificial habitat or on adjacent natural habitat both north and south of the highway. There were no significant differences in levels of heterozygosity, allelic richness or inbreeding between natural (n?=?7) and artificial sites (n?=?8). However, pikas residing on artificial habitat exhibited significantly higher relatedness estimates. Bayesian clustering analyses revealed two distinct genetic units corresponding to north and south of the highway, with further substructure detected in the south. Likewise, high genetic friction was detected in the central region of the area, largely corresponding to the highway and modified landscape associated with the mine. At a finer scale, pairwise estimates of differentiation and migration rates suggest little gene flow may be occurring among sites across the sampling area, with some evidence for directional migration from artificial to natural sites. Overall, artificial habitat has been successful in promoting occupancy for American pikas, however, barriers to gene flow likely associated with resource extraction and road construction limit connectivity across the landscape.