Physiology of thermoregulation in the pika, Ochotona princeps

Life on the edge-a changing genetic landscape within an iconic American pika metapopulation over the last half century

Declines and extirpations of American pika (Ochotona princeps) populations at historically occupied sites started being documented in the literature during the early 2000s. Commensurate with global climate change, many of these losses at peripheral and lower elevation sites have been associated with changes in ambient air temperature and precipitation regimes. Here, we report on a decline in available genetic resources for an iconic American pika metapopulation, located at the southwestern edge of the species distribution in the Bodie Hills of eastern California, USA. Composed of highly fragmented habitat created by hard rock mining, the ore dumps at this site were likely colonized by pikas around the end of the 19th century from nearby natural talus outcrops. Genetic data extracted from both contemporary samples and archived natural history collections allowed us to track population and patch-level genetic diversity for Bodie pikas across three distinct sampling points during the last half- century (1948-1949, 1988-1991, 2013-2015). Reductions in within-population allelic diversity and expected heterozygosity were observed across the full time period. More extensive sampling of extant patches during the 1988-1991 and 2013-2015 periods revealed an increase in population structure and a reduction in effective population size. Furthermore, census records from the last 51 years as well as archived museum samples collected in 1947 from a nearby pika population in the Wassuk range (Nevada, USA) provide further support of the increasing isolation and genetic coalescence occurring in this region. This study highlights the importance of museum samples and long-term monitoring in contextualizing our understanding of population viability.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Return of the pika: American pikas re-occupy long-extirpated, warm locations

American pikas (Ochotona princeps), small mammals related to rabbits, occur in mountainous regions of western North America, where they live in shattered‐rock habitats (talus). Aspects of their physiology and life history create situations that appear to put pikas at risk from warming climates. Some low‐elevation, warm sites that historically harbored pikas have become extirpated, and the assumption is that these will not be re‐colonized under current climate trends. Unexpectedly, in 2021, we found that pikas had re‐colonized two very warm, low‐elevation, dry sites in eastern California, USA, in the Bodie Mountains and Mono Craters. Resident pikas appear to have been absent at both sites for ≥10 years. These findings suggest that pikas, which are normally diurnally active, are able to overcome thermal dispersal barriers and re‐colonize long‐extirpated sites, perhaps by moving during cool nights. Our data also highlight the often unrecognized suitability of pika habitat in warm regions where the interiors of taluses can remain stably cool even when external air temperatures are hot.

Factors influencing distributional shifts and abundance at the range core of a climate-sensitive mammal

Species are frequently responding to contemporary climate change by shifting to higher elevations and poleward to track suitable climate space. However, depending on local conditions and species' sensitivity, the nature of these shifts can be highly variable and difficult to predict. Here, we examine how the American pika (Ochotona princeps), a philopatric, montane lagomorph, responds to climatic gradients at three spatial scales. Using mixed-effects modeling in an information-theoretic approach, we evaluated a priori model suites regarding predictors of site occupancy, relative abundance, and elevational-range retraction across 760 talus patches, nested within 64 watersheds across the Northern Rocky Mountains of North America, during 2017-2020. The top environmental predictors differed across these response metrics. Warmer temperatures in summer and winter were associated with lower occupancy, lower relative abundances, and greater elevational retraction across watersheds. Occupancy was also strongly influenced by habitat patch size, but only when combined with climate metrics such as actual evapotranspiration. Using a second analytical approach, acute heat stress and summer precipitation best explained retraction residuals (i.e., the relative extent of retraction given the original elevational range of occupancy). Despite the study domain occurring near the species' geographic-range center, where populations might have higher abundances and be at lower risk of climate-related stress, 33.9% of patches showed evidence of recent extirpations. Pika-extirpated sites averaged 1.44℃ warmer in summer than did occupied sites. Additionally, the minimum elevation of pika occupancy has retracted upslope in 69% of watersheds (mean: 281 m). Our results emphasize the nuance associated with evaluating species' range dynamics in response to climate gradients, variability, and temperature exceedances, especially in regions where species occupy gradients of conditions that may constitute multiple range edges. Furthermore, this study highlights the importance of evaluating diverse drivers across response metrics to improve the predictive accuracy of widely used, correlative models.

Surviving winter on the Qinghai-Tibetan Plateau: Pikas suppress energy demands and exploit yak feces to survive winter

The Qinghai-Tibetan Plateau, with low precipitation, low oxygen partial pressure, and temperatures routinely dropping below -30 °C in winter, presents several physiological challenges to its fauna. Yet it is home to many endemic mammalian species, including the plateau pika (Ochotona curzoniae). How these small animals that are incapable of hibernation survive the winter is an enigma. Measurements of daily energy expenditure (DEE) using the doubly labeled water method show that pikas suppress their DEE during winter. At the same body weight, pikas in winter expend 29.7% less than in summer, despite ambient temperatures being approximately 25 °C lower. Combined with resting metabolic rates (RMRs), this gives them an exceptionally low metabolic scope in winter (DEE/RMRt = 1.60 ± 0.30; RMRt is resting metabolic rate at thermoneutrality). Using implanted body temperature loggers and filming in the wild, we show that this is achieved by reducing body temperature and physical activity. Thyroid hormone (T₃ and T₄) measurements indicate this metabolic suppression is probably mediated via the thyroid axis. Winter activity was lower at sites where domestic yak (Bos grunniens) densities were higher. Pikas supplement their food intake at these sites by eating yak feces, demonstrated by direct observation, identification of yak DNA in pika stomach contents, and greater convergence in the yak/pika microbiotas in winter. This interspecific coprophagy allows pikas to thrive where yak are abundant and partially explains why pika densities are higher where domestic yak, their supposed direct competitors for food, are more abundant.

Field Physiology: Studying Organismal Function in the Natural Environment

Continuous physiological measurements collected in field settings are essential to understand baseline, free-ranging physiology, physiological range and variability, and the physiological responses of organisms to disturbances. This article presents a current summary of the available technologies to continuously measure the direct physiological parameters in the field at high-resolution/instantaneous timescales from freely behaving animals. There is a particular focus on advantages versus disadvantages of available methods as well as emerging technologies "on the horizon" that may have been validated in captive or laboratory-based scenarios but have yet to be applied in the wild. Systems to record physiological variables from free-ranging animals are reviewed, including radio (VHF/UFH) telemetry, acoustic telemetry, and dataloggers. Physiological parameters that have been continuously measured in the field are addressed in seven sections including heart rate and electrocardiography (ECG); electromyography (EMG); electroencephalography (EEG); body temperature; respiratory, blood, and muscle oxygen; gastric pH and motility; and blood pressure and flow. The primary focal sections are heart rate and temperature as these can be, and have been, extensively studied in free-ranging organisms. Predicted aspects of future innovation in physiological monitoring are also discussed. The article concludes with an overview of best practices and points to consider regarding experimental designs, cautions, and effects on animals. © 2021 American Physiological Society. Compr Physiol 11:1979-2015, 2021.

Freezing in a warming climate: Marked declines of a subnivean hibernator after a snow drought

Recent snow droughts associated with unusually warm winters are predicted to increase in frequency and affect species dependent upon snowpack for winter survival. Changes in populations of some cold-adapted species have been attributed to heat stress or indirect effects on habitat from unusually warm summers, but little is known about the importance of winter weather to population dynamics and how responses to snow drought vary among sympatric species. We evaluated changes in abundance of hoary marmots (Marmota caligata) over a period that included a year of record-low snowpack to identify mechanisms associated with weather and snowpack. To consider interspecies comparisons, our analysis used the same a priori model set as a concurrent study that evaluated responses of American pikas (Ochotona princeps) to weather and snowpack in the same study area of North Cascades National Park, Washington, USA. We hypothesized that marmot abundance reflected mechanisms related to heat stress, cold stress, cold exposure without an insulating snowpack, snowpack duration, atmospheric moisture, growing-season precipitation, or select combinations of these mechanisms. Changes in marmot abundances included a 74% decline from 2007 to 2016 and were best explained by an interaction of chronic dryness with exposure to acute cold without snowpack in winter. Physiological stress during hibernation from exposure to cold, dry air appeared to be the most likely mechanism of change in marmot abundance. Alternative mechanisms associated with changes to winter weather, including early emergence from hibernation or altered vegetation dynamics, had less support. A post hoc assessment of vegetative phenology and productivity did not support vegetation dynamics as a primary driver of marmot abundance across years. Although marmot and pika abundances were explained by strikingly similar models over periods of many years, details of the mechanisms involved likely differ between species because pika abundances increased in areas where marmots declined. Such differences may lead to diverging geographic distributions of these species as global change continues.

Genomic variation in the American pika: signatures of geographic isolation and implications for conservation

Background

Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada.

Results

Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θ_W = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θ_W = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θ_W = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’s D was positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide.

Conclusions

Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.

Conservation status of American pikas (Ochotona princeps)

The American pika (Ochotona princeps) is commonly perceived as a species that is at high risk of extinction due to climate change. The purpose of this review is two-fold: to evaluate the claim that climate change is threatening pikas with extinction, and to summarize the conservation status of the American pika. Most American pikas inhabit major cordilleras, such as the Rocky Mountain, Sierra Nevada, and Cascade ranges. Occupancy of potential pika habitat in these ranges is uniformly high and no discernible climate signal has been found that discriminates between the many occupied and relatively few unoccupied sites that have been recently surveyed. Pikas therefore are thriving across most of their range. The story differs in more marginal parts of the species range, primarily across the Great Basin, where a higher percentage of available habitat is unoccupied. A comprehensive review of Great Basin pikas revealed that occupied sites, sites of recent extirpation, and old sites, were regularly found within the same geographic and climatic space as extant sites, and suggested that pikas in the Great Basin tolerated a broader set of habitat and climatic conditions than previously understood. Studies of a small subset of extirpated sites in the Great Basin and in California found that climate variables (most notably measures of hot temperature) were associated more often with extirpated sites than occupied sites. Importantly, upward contraction of the lower elevation boundary also was found at some sites. However, models that incorporated variables other than climate (such as availability of upslope talus habitat) often were better predictors of site persistence. Many extirpations occurred on small habitat patches, which were subject to stochastic extinction, as informed by a long-term pika metapopulation study in Bodie, California. In addition, several sites may have been compromised by cattle grazing or other anthropogenic factors. In contrast, several low, hot sites (Bodie, Mono Craters, Craters of the Moon National Monument and Preserve, Lava Beds National Monument, Columbia River Gorge) retain active pika populations, demonstrating the adaptive capacity and resilience of pikas in response to adverse environmental conditions. Pikas cope with warm temperatures by retreating into cool interstices of their talus habitat and augment their restricted daytime foraging with nocturnal activity. Pikas exhibit significant flexibility in their foraging tactics and are highly selective in their choice of available vegetation. The trait that places pikas at greatest risk from climate change is their poor dispersal capability. Dispersal is more restricted in hotter environments, and isolated low-elevation sites that become extirpated are unlikely to be recolonized in a warming climate. The narrative that American pikas are going extinct appears to be an overreach. Pikas are doing well across most of their range, but there are limited, low-elevation losses that are likely to be permanent in what is currently marginal pika habitat. The resilience of pikas in the face of climate change, and their ability or inability to persist in marginal, hot environments, will continue to contribute to our understanding of the impact of climate change on individual species.

Vector-Borne Pathogens in Ectoparasites Collected from High-Elevation Pika Populations

The American pika, Ochotona princeps, is projected to decline throughout North America as climate change reduces its range, and pikas have already disappeared from several locations. In addition to climate, disease spillover from lower elevation mammalian species might affect pikas. We sampled pika fleas in Colorado and Montana across elevations ranging from 2896 to 3612 m and screened them for the presence of DNA from rodent-associated bacterial pathogens (Bartonella species and Yersinia pestis) to test the hypothesis that flea exchange between pikas and rodents may lead to occurrence of rodent-associated pathogens in pika ectoparasites. We collected 275 fleas from 74 individual pikas at 5 sites in Colorado and one site in Montana. We found that 5.5% of 275 pika fleas in this study tested positive for rodent-associated Bartonella DNA but that variation in Bartonella infection prevalence in fleas among sites was not driven by elevation. Specifically, we detected DNA sequences from two loci (gltA and rpoB) that are most similar to Bartonella grahamii isolates collected from rodents in Canada. We did not detect Y. pestis DNA in our survey. Our results demonstrate evidence of rodent-associated flea-borne bacteria in pika fleas. These findings are also consistent with the hypothesis that rodent-associated pathogens could be acquired by pikas. Flea-borne pathogen spillover from rodents to pikas has the potential to exacerbate the more direct effects of climate that have been suggested to drive pika declines.

Investigating the role of environment in pika (Ochotona) body size patterns across taxonomic levels, space, and time

Body size is an important trait in animals because it influences a multitude of additional life history traits. The causal mechanisms underlying body size patterns across spatial, temporal, and taxonomic hierarchies are debated, and of renewed interest in this era of climate change. Here, we tested multiple hypotheses regarding body mass patterns at the intraspecific and interspecific levels. We investigated body size patterns within a climate-sensitive small mammal species, Ochotona princeps (n = 2,873 individuals), across their range with local environmental variables. We also examined body mass of populations over time to determine if body size has evolved in situ in response to environmental change. At the interspecific level we compared the mean mass of 26 pika species (genus Ochotona) to determine if environmental temperatures, food availability, habitat variability, or range area influence body size. We found correlations between temperature, vegetation, and particularly precipitation variables, with body mass within O. princeps, but no linear relationship between body size and any climate or habitat variable for Ochotona species. Body size trends in relation to climate were stronger at the intraspecific than the interspecific level. Our results suggest that body size within O. princeps likely is related to food availability, and that body size evolution is not always a viable response to temperature change. Different mechanisms may be driving body size at the interspecific and intraspecific levels and factors other than environment, such as biotic interactions, may also be influential in determining body size over space and time.

Ecological consequences of anomalies in atmospheric moisture and snowpack

Although increased frequency of extreme-weather events is one of the most secure predictions associated with contemporary climate change, effects of such events on distribution and abundance of climate-sensitive species remain poorly understood. Montane ecosystems may be especially sensitive to extreme weather because of complex abiotic and biotic interactions that propagate from climate-driven reductions in snowpack. Snowpack not only protects subnivean biotas from extreme cold, but also influences forage availability through timing of melt-off and water availability. We related relative abundances of an alpine mammal, the American pika (Ochotona princeps), to measures of weather and snowpack dynamics over an 8-yr period that included before and after a year of record-low snowpack in Washington, USA. We sought to (1) quantify any change in pika abundance associated with the snowpack anomaly and (2) identify aspects of weather and snowpack that influenced abundance of pikas. Pikas showed a 1-yr lag response to the snowpack anomaly and exhibited marked declines in abundance at elevations below 1,400 m simultaneous with increased abundances at higher elevations. Atmospheric moisture, indexed by vapor pressure deficit (VPD), was especially important, evidenced by strong support for the top-ranked model that included the interaction of VPD with snowpack duration. Notably, our novel application of VPD from gridded climate data for analyses of animal abundances shows strong potential for improving species distribution models because VPD represents an important aspect of weather that influences the physiology and habitat of biota. Pikas were apparently affected by cold stress without snowpack at mid elevations, whereas changes to forage associated with snowpack and VPD were influential at high and low elevations. Our results reveal context dependency in pika responses to weather and illustrate how snow drought can lead to rapid change in the abundance of subnivean animals.

Behavioural plasticity modulates temperature-related constraints on foraging time for a montane mammal

Contemporary climate change is altering temperature profiles across the globe. Increasing temperatures can reduce the amount of time during which conditions are suitable for animals to engage in essential activities, such as securing food. Behavioural plasticity, the ability to alter behaviour in response to the environment, may provide animals with a tool to adjust to changes in the availability of suitable thermal conditions. The extent to which individuals can alter fitness-enhancing behaviours, such as food collection, to proximately buffer variation in temperature, however, remains unclear. Even less well understood are the potential performance advantages of flexible strategies among endotherms. We examined the degree to which individuals altered rates of food collection in response to temperature, and two potential benefits, using the American pika (Ochotona princeps), a temperature-sensitive, food-hoarding mammal, as a model. From July-September 2013-2015, we used motion-activated cameras and in situ temperature loggers to examine pika food-caching activity for 72 individuals across 10 sites in the central Rocky Mountains, USA. We quantified % nitrogen by cache volume as a metric of cache quality, and the number of events during which pikas were active in temperatures ≥25°C as a measure of potential thermoregulatory stress. We found a strong negative effect of temperature on the rate at which pikas cached food. Individual responses to temperature varied substantially in both the level of food-collecting activity and in the degree to which individuals shifted activity with warming temperature. After accounting for available foraging time, individuals that exhibited greater plasticity collected a comparable amount of nitrogen, while simultaneously experiencing fewer occasions in which temperatures eclipsed estimated thermal tolerances. By varying food-collection norms of reaction, individuals were able to plastically respond to temperature-driven reductions in foraging time. Through this increased flexibility, individuals amassed food caches of comparable quality, while minimizing exposure to potentially stressful thermal conditions. Our results suggest that, given sufficient resource quality and availability, plasticity in foraging activity may help temperature-limited endotherms adjust to climate-related constraints on foraging time.

Altai pika (Ochotona alpina) alarm calls: individual acoustic variation and the phenomenon of call-synchronous ear folding behavior

Non-hibernating pikas collect winter food reserves and store them in hay piles. Individualization of alarm calls might allow discrimination between colony members and conspecifics trying to steal food items from a colony pile. We investigated vocal posture, vocal tract length, and individual acoustic variation of alarm calls, emitted by wild-living Altai pikas Ochotona alpina toward a researcher. Recording started when a pika started calling and lasted as long as possible. The alarm call series of 442 individual callers from different colonies consisted of discrete short (0.073-0.157 s), high-frequency (7.31-15.46 kHz), and frequency-modulated calls separated by irregular intervals. Analysis of 442 discrete calls, the second of each series, revealed that 44.34% calls lacked nonlinear phenomena, in 7.02% nonlinear phenomena covered less than half of call duration, and in 48.64% nonlinear phenomena covered more than half of call duration. Peak frequencies varied among individuals but always fitted one of three maxima corresponding to the vocal tract resonance frequencies (formants) calculated for an estimated 45-mm oral vocal tract. Discriminant analysis using variables of 8 calls per series of 36 different callers, each from a different colony, correctly assigned over 90% of the calls to individuals. Consequently, Altai pika alarm calls are individualistic and nonlinear phenomena might further increase this acoustic individualization. Additionally, video analysis revealed a call-synchronous, very fast (0.13-0.23 s) folding, depression, and subsequent re-expansion of the pinna confirming an earlier report of this behavior that apparently contributes to protecting the hearing apparatus from damage by the self-generated high-intensity alarm calls.

Seasonal temperature acclimatization in a semi-fossorial mammal and the role of burrows as thermal refuges

Small mammals in habitats with strong seasonal variation in the thermal environment often exhibit physiological and behavioral adaptations for coping with thermal extremes and reducing thermoregulatory costs. Burrows are especially important for providing thermal refuge when above-ground temperatures require high regulatory costs (e.g., water or energy) or exceed the physiological tolerances of an organism. Our objective was to explore the role of burrows as thermal refuges for a small endotherm, the pygmy rabbit (Brachylagus idahoensis), during the summer and winter by quantifying energetic costs associated with resting above and below ground. We used indirect calorimetry to determine the relationship between energy expenditure and ambient temperature over a range of temperatures that pygmy rabbits experience in their natural habitat. We also measured the temperature of above- and below-ground rest sites used by pygmy rabbits in eastern Idaho, USA, during summer and winter and estimated the seasonal thermoregulatory costs of resting in the two microsites. Although pygmy rabbits demonstrated seasonal physiological acclimatization, the burrow was an important thermal refuge, especially in winter. Thermoregulatory costs were lower inside the burrow than in above-ground rest sites for more than 50% of the winter season. In contrast, thermal heterogeneity provided by above-ground rest sites during summer reduced the role of burrows as a thermal refuge during all but the hottest periods of the afternoon. Our findings contribute to an understanding of the ecology of small mammals in seasonal environments and demonstrate the importance of burrows as thermal refuge for pygmy rabbits.

Apparent climate-mediated loss and fragmentation of core habitat of the American pika in the Northern Sierra Nevada, California, USA

Contemporary climate change has been widely documented as the apparent cause of range contraction at the edge of many species distributions but documentation of climate change as a cause of extirpation and fragmentation of the interior of a species’ core habitat has been lacking. Here, we report the extirpation of the American pika (Ochotona princeps), a temperature-sensitive small mammal, from a 165-km² area located within its core habitat in California’s Sierra Nevada mountains. While sites surrounding the area still maintain pikas, radiocarbon analyses of pika fecal pellets recovered within this area indicate that former patch occupancy ranges from before 1955, the beginning of the atmospheric spike in radiocarbon associated with above ground atomic bomb testing, to c. 1991. Despite an abundance of suitable rocky habitat climate warming appears to have precipitated their demise. Weather station data reveal a 1.9°C rise in local temperature and a significant decline in snowpack over the period of record, 1910–2015, pushing pika habitat into increasingly tenuous climate conditions during the period of extirpation. This is among the first accounts of an apparently climate-mediated, modern extirpation of a species from an interior portion of its geographic distribution, resulting in habitat fragmentation, and is the largest area yet reported for a modern-era pika extirpation. Our finding provides empirical support to model projections, indicating that even core areas of species habitat are vulnerable to climate change within a timeframe of decades.

Comparative analyses of basal rate of metabolism in mammals: data selection does matter

Basal rate of metabolism (BMR) is a physiological parameter that should be measured under strictly defined experimental conditions. In comparative analyses among mammals BMR is widely used as an index of the intensity of the metabolic machinery or as a proxy for energy expenditure. Many databases with BMR values for mammals are available, but the criteria used to select metabolic data as BMR estimates have often varied and the potential effect of this variability has rarely been questioned. We provide a new, expanded BMR database reflecting compliance with standard criteria (resting, postabsorptive state; thermal neutrality; adult, non-reproductive status for females) and examine potential effects of differential selectivity on the results of comparative analyses. The database includes 1739 different entries for 817 species of mammals, compiled from the original sources. It provides information permitting assessment of the validity of each estimate and presents the value closest to a proper BMR for each entry. Using different selection criteria, several alternative data sets were extracted and used in comparative analyses of (i) the scaling of BMR to body mass and (ii) the relationship between brain mass and BMR. It was expected that results would be especially dependent on selection criteria with small sample sizes and with relatively weak relationships. Phylogenetically informed regression (phylogenetic generalized least squares, PGLS) was applied to the alternative data sets for several different clades (Mammalia, Eutheria, Metatheria, or individual orders). For Mammalia, a 'subsampling procedure' was also applied, in which random subsamples of different sample sizes were taken from each original data set and successively analysed. In each case, two data sets with identical sample size and species, but comprising BMR data with different degrees of reliability, were compared. Selection criteria had minor effects on scaling equations computed for large clades (Mammalia, Eutheria, Metatheria), although less-reliable estimates of BMR were generally about 12-20% larger than more-reliable ones. Larger effects were found with more-limited clades, such as sciuromorph rodents. For the relationship between BMR and brain mass the results of comparative analyses were found to depend strongly on the data set used, especially with more-limited, order-level clades. In fact, with small sample sizes (e.g. <100) results often appeared erratic. Subsampling revealed that sample size has a non-linear effect on the probability of a zero slope for a given relationship. Depending on the species included, results could differ dramatically, especially with small sample sizes. Overall, our findings indicate a need for due diligence when selecting BMR estimates and caution regarding results (even if seemingly significant) with small sample sizes.

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.

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.

Predictors of Current and Longer-Term Patterns of Abundance of American Pikas (Ochotona princeps) across a Leading-Edge Protected Area

American pikas (Ochotona princeps) have been heralded as indicators of montane-mammal response to contemporary climate change. Pikas no longer occupy the driest and lowest-elevation sites in numerous parts of their geographic range. Conversely, pikas have exhibited higher rates of occupancy and persistence in Rocky Mountain and Sierra Nevada montane 'mainlands'. Research and monitoring efforts on pikas across the western USA have collectively shown the nuance and complexity with which climate will often act on species in diverse topographic and climatic contexts. However, to date no studies have investigated habitat, distribution, and abundance of pikas across hundreds of sites within a remote wilderness area. Additionally, relatively little is known about whether climate acts most strongly on pikas through direct or indirect (e.g., vegetation-mediated) mechanisms. During 2007-2009, we collectively hiked >16,000 km throughout the 410,077-ha Glacier National Park, Montana, USA, in an effort to identify topographic, microrefugial, and vegetative characteristics predictive of pika abundance. We identified 411 apparently pika-suitable habitat patches with binoculars (in situ), and surveyed 314 of them for pika signs. Ranking of alternative logistic-regression models based on AICc scores revealed that short-term pika abundances were positively associated with intermediate elevations, greater cover of mosses, and taller forbs, and decreased each year, for a total decline of 68% during the three-year study; whereas longer-term abundances were associated only with static variables (longitude, elevation, gradient) and were lower on north-facing slopes. Earlier Julian date and time of day of the survey (i.e., midday vs. not) were associated with lower observed pika abundance. We recommend that wildlife monitoring account for this seasonal and diel variation when surveying pikas. Broad-scale information on status and abundance determinants of montane mammals, especially for remote protected areas, is crucial for land and wildlife-resource managers trying to anticipate mammalian responses to climate change.

LaGomiCs-Lagomorph Genomics Consortium: An International Collaborative Effort for Sequencing the Genomes of an Entire Mammalian Order

The order Lagomorpha comprises about 90 living species, divided in 2 families: the pikas (Family Ochotonidae), and the rabbits, hares, and jackrabbits (Family Leporidae). Lagomorphs are important economically and scientifically as major human food resources, valued game species, pests of agricultural significance, model laboratory animals, and key elements in food webs. A quarter of the lagomorph species are listed as threatened. They are native to all continents except Antarctica, and occur up to 5000 m above sea level, from the equator to the Arctic, spanning a wide range of environmental conditions. The order has notable taxonomic problems presenting significant difficulties for defining a species due to broad phenotypic variation, overlap of morphological characteristics, and relatively recent speciation events. At present, only the genomes of 2 species, the European rabbit (Oryctolagus cuniculus) and American pika (Ochotona princeps) have been sequenced and assembled. Starting from a paucity of genome information, the main scientific aim of the Lagomorph Genomics Consortium (LaGomiCs), born from a cooperative initiative of the European COST Action "A Collaborative European Network on Rabbit Genome Biology-RGB-Net" and the World Lagomorph Society (WLS), is to provide an international framework for the sequencing of the genome of all extant and selected extinct lagomorphs. Sequencing the genomes of an entire order will provide a large amount of information to address biological problems not only related to lagomorphs but also to all mammals. We present current and planned sequencing programs and outline the final objective of LaGomiCs possible through broad international collaboration.

Low genetic diversity, restricted dispersal, and elevation-specific patterns of population decline in American pikas in an atypical environment

In the face of climate change, there is a growing need for research into the ability of organisms to persist at the limits of their bioclimatic envelope. American pikas ( Ochotona princeps ) have emerged as a focal mammalian species for investigating extinction risk related to climate change; however, most studies have occurred in characteristic alpine talus habitat within the range core. In the Columbia River Gorge (CRG), Oregon, American pikas inhabit low-elevation talus slopes previously considered outside the species’ bioclimatic range. We used microsatellite genotypic data to reconstruct levels of genetic variation, population connectivity, and demographic history at 11 CRG sites spanning an elevational gradient (104–1,292 m). Sampled sites separated into 2 genetic clusters largely explained by elevation, topography, and geographic proximity, with pairwise estimates of differentiation and migration rates suggesting little gene flow may be occurring. Sites were characterized by levels of allelic richness and heterozygosity substantially lower than values reported at characteristic alpine sites from the range core. Evidence of recent demographic contraction was found almost exclusively at high-elevation sites despite these areas being considered refuges from climate warming in more typical habitat. Given their unique genetic characteristics and persistence in an atypical environment, the CRG pika populations likely constitute a significant component of intraspecific biodiversity with high conservation value.

Modeling behavioral thermoregulation in a climate change sentinel

When possible, many species will shift in elevation or latitude in response to rising temperatures. However, before such shifts occur, individuals will first tolerate environmental change and then modify their behavior to maintain heat balance. Behavioral thermoregulation allows animals a range of climatic tolerances and makes predicting geographic responses under future warming scenarios challenging. Because behavioral modification may reduce an individual's fecundity by, for example, limiting foraging time and thus caloric intake, we must consider the range of behavioral options available for thermoregulation to accurately predict climate change impacts on individual species. To date, few studies have identified mechanistic links between an organism's daily activities and the need to thermoregulate. We used a biophysical model, Niche Mapper, to mechanistically model microclimate conditions and thermoregulatory behavior for a temperature‐sensitive mammal, the American pika (Ochotona princeps). Niche Mapper accurately simulated microclimate conditions, as well as empirical metabolic chamber data for a range of fur properties, animal sizes, and environmental parameters. Niche Mapper predicted pikas would be behaviorally constrained because of the need to thermoregulate during the hottest times of the day. We also showed that pikas at low elevations could receive energetic benefits by being smaller in size and maintaining summer pelage during longer stretches of the active season under a future warming scenario. We observed pika behavior for 288 h in Glacier National Park, Montana, and thermally characterized their rocky, montane environment. We found that pikas were most active when temperatures were cooler, and at sites characterized by high elevations and north‐facing slopes. Pikas became significantly less active across a suite of behaviors in the field when temperatures surpassed 20°C, which supported a metabolic threshold predicted by Niche Mapper. In general, mechanistic predictions and empirical observations were congruent. This research is unique in providing both an empirical and mechanistic description of the effects of temperature on a mammalian sentinel of climate change, the American pika. Our results suggest that previously underinvestigated characteristics, specifically fur properties and body size, may play critical roles in pika populations' response to climate change. We also demonstrate the potential importance of considering behavioral thermoregulation and microclimate variability when predicting animal responses to climate change.

The role of thermal physiology in recent declines of birds in a biodiversity hotspot

We investigated whether observed avian range contractions and population declines in the Fynbos biome of South Africa were mechanistically linked to recent climate warming. We aimed to determine whether there were correlations between preferred temperature envelope, or changes in temperature within species' ranges, and recent changes in range and population size, for 12 Fynbos-resident bird species, including six that are endemic to the biome. We then measured the physiological responses of each species at air temperatures ranging from 24 to 42°C to determine whether physiological thermal thresholds could provide a mechanistic explanation for observed population trends. Our data show that Fynbos-endemic species occupying the coolest regions experienced the greatest recent reductions in range and population size (>30% range reduction between 1991 and the present). In addition, species experiencing the largest increases in air temperature within their ranges showed the greatest declines. However, evidence for a physiological mechanistic link between warming and population declines was equivocal, with only the larger species showing low thermal thresholds for their body mass, compared with other birds globally. In addition, some species appear more vulnerable than others to air temperatures in their ranges above physiological thermal thresholds. Of these, the high-altitude specialist Cape rockjumper (Chaetops frenatus) seems most at risk from climate warming. This species showed: (i) the lowest threshold for increasing evaporative water loss at high temperatures; and (ii) population declines specifically in those regions of its range recording significant warming trends. Our findings suggest that caution must be taken when attributing causality explicitly to thermal stress, even when population trends are clearly correlated with rates of warming. Studies explicitly investigating the mechanisms underlying such correlations will be key to appropriate conservation planning.

American pika in a low-elevation lava landscape: expanding the known distribution of a temperature-sensitive species

In 2010, the American pika (Ochotona princeps fenisex) was denied federal protection based on limited evidence of persistence in low-elevation environments. Studies in nonalpine areas have been limited to relatively few environments, and it is unclear whether patterns observed elsewhere (e.g., Bodie, CA) represent other nonalpine habitats. This study was designed to establish pika presence in a new location, determine distribution within the surveyed area, and evaluate influences of elevation, vegetation, lava complexity, and distance to habitat edge on pika site occupancy. In 2011 and 2012, we conducted surveys for American pika on four distinct subalpine lava flows of Newberry National Volcanic Monument, Oregon, USA. Field surveys were conducted at predetermined locations within lava flows via silent observation and active searching for pika sign. Site habitat characteristics were included as predictors of occupancy in multinomial regression models. Above and belowground temperatures were recorded at a subsample of pika detection sites. Pika were detected in 26% (2011) and 19% (2012) of survey plots. Seventy-four pika were detected outside survey plot boundaries. Lava complexity was the strongest predictor of pika occurrence, where pika were up to seven times more likely to occur in the most complicated lava formations. Pika were two times more likely to occur with increasing elevation, although they were found at all elevations in the study area. This study expands the known distribution of the species and provides additional evidence for persistence in nonalpine habitats. Results partially support the predictive occupancy model developed for pika at Craters of the Moon National Monument, another lava environment. Characteristics of the lava environment clearly influence pika site occupancy, but habitat variables reported as important in other studies were inconclusive here. Further work is needed to gain a better understanding of the species’ current distribution and ability to persist under future climate conditions.

Climate Tolerances and Habitat Requirements Jointly Shape the Elevational Distribution of the American Pika (Ochotona princeps), with Implications for Climate Change Effects

Some of the most compelling examples of ecological responses to climate change are elevational range shifts of individual species, which have been observed throughout the world. A growing body of evidence, however, suggests substantial mediation of simple range shifts due to climate change by other limiting factors. Understanding limiting factors for a species within different contexts, therefore, is critical for predicting responses to climate change. The American pika (Ochotona princeps) is an ideal species for investigating distributions in relation to climate because of their unusual and well-understood natural history as well as observed shifts to higher elevation in parts of their range. We tested three hypotheses for the climatic or habitat characteristics that may limit pika presence and abundance: summer heat, winter snowpack, and forage availability. We performed these tests using an index of pika abundance gathered in a region where environmental influences on pika distribution have not been well-characterized. We estimated relative pika abundance via scat surveys and quantified climatic and habitat characteristics across two North-Central Rocky Mountain Ranges, the Wind River and Bighorn ranges in Wyoming, USA. Pika scat density was highest at mid-elevations and increased linearly with forage availability in both ranges. Scat density also increased with temperatures conducive to forage plant growth, and showed a unimodal relationship with the number of days below -5°C, which is modulated by insulating snowpack. Our results provide support for both the forage availability and winter snowpack hypotheses. Especially in montane systems, considering the context-dependent nature of climate effects across regions and elevations as well as interactions between climatic and other critical habitat characteristics, will be essential for predicting future species distributions.

Relating sub-surface ice features to physiological stress in a climate sensitive mammal, the American pika (Ochotona princeps)

The American pika (Ochotona princeps) is considered a sentinel species for detecting ecological effects of climate change. Pikas are declining within a large portion of their range, and ongoing research suggests loss of sub-surface ice as a mechanism. However, no studies have demonstrated physiological responses of pikas to sub-surface ice features. Here we present the first analysis of physiological stress in pikas living in and adjacent to habitats underlain by ice. Fresh fecal samples were collected non-invasively from two adjacent sites in the Rocky Mountains (one with sub-surface ice and one without) and analyzed for glucocorticoid metabolites (GCM). We also measured sub-surface microclimates in each habitat. Results indicate lower GCM concentration in sites with sub-surface ice, suggesting that pikas are less stressed in favorable microclimates resulting from sub-surface ice features. GCM response was well predicted by habitat characteristics associated with sub-surface ice features, such as lower mean summer temperatures. These results suggest that pikas inhabiting areas without sub-surface ice features are experiencing higher levels of physiological stress and may be more susceptible to changing climates. Although post-deposition environmental effects can confound analyses based on fecal GCM, we found no evidence for such effects in this study. Sub-surface ice features are key to water cycling and storage and will likely represent an increasingly important component of water resources in a warming climate. Fecal samples collected from additional watersheds as part of current pika monitoring programs could be used to further characterize relationships between pika stress and sub-surface ice features.

Radiocarbon dating of American pika fecal pellets provides insights into population extirpations and climate refugia

The American pika (Ochotona princeps) has become a species of concern for its sensitivity to warm temperatures and potential vulnerability to global warming. We explored the value of radiocarbon dating of fecal pellets to address questions of population persistence and timing of site extirpation. Carbon was extracted from pellets collected at 43 locations in the western Great Basin, USA, including three known occupied sites and 40 sites of uncertain status at range margins or where previous studies indicated the species is vulnerable. We resolved calibrated dates with high precision (within several years), most of which fell in the period of the mid-late 20th century bomb curve. The two-sided nature of the bomb curve renders far- and near-side dates of equal probability, which are separated by one to four decades. We document methods for narrowing resolution to one age range, including stratigraphic analysis of vegetation collected from pika haypiles. No evidence was found for biases in atmospheric 14C levels due to fossil-derived or industrial CO2 contamination. Radiocarbon dating indicated that pellets can persist for >59 years; known occupied sites resolved contemporary dates. Using combined evidence from field observations and radiocarbon dating, and the Bodie Mountains as an example, we propose a historical biogeographic scenario for pikas in minor Great Basin mountain ranges adjacent to major cordillera, wherein historical climate variability led to cycles of extirpation and recolonization during alternating cool and warm centuries. Using this model to inform future dynamics for small ranges in biogeographic settings similar to the Bodie Mountains in California, extirpation of pikas appears highly likely under directional warming trends projected for the next century, even while populations in extensive cordillera (e.g., Sierra Nevada, Rocky Mountains, Cascade Range) are likely to remain viable due to extensive, diverse habitat and high connectivity.

Stress hormone concentration in Rocky Mountain populations of the American pika (Ochotona princeps)

The American pika (Ochotona princeps) is considered a sentinel species for detecting ecological effects of climate change. Pikas are declining within a large portion of their range, but previous studies have focused only on local pika extinction as a metric of change. We designed a procedure which can provide an earlier warning signal, based on non-invasive sampling and analysis of physiological stress in living pikas. Pikas were sampled at several locations in the Rocky Mountains for the measurement of glucocorticoid metabolites (GCMs) in faeces. Using a time series of faecal pellets from 12 individuals, we detected a significant increase in faecal GCM level in response to capture, thus biologically validating the use of a corticosterone enzyme immunoassay. We also established baseline, peak, and post-peak GCM concentrations for pikas in the Rocky Mountains, which varied according to gender and individual. This is the first study to measure stress hormone metabolites in any species of pika. The methods developed and validated in this study can be used to add non-invasive measurements of physiological stress to pika monitoring programmes and other research designed to assess pika vulnerability to predicted changes in climate. Pika monitoring programmes currently in place use a protocol that relates current site use by pikas with data on local habitat characteristics, such as elevation, to infer potential effects of climate change. Data generated by these monitoring studies can be used to identify the trends in site use by pikas in relationship to habitat covariates. However, this approach does not take into account the role of behavioural thermoregulation and the pika's use of microhabitats to ameliorate variations in climate. Incorporating a stress metric, such as GCM concentration, will provide relatively direct evidence for or against the hypothesis that pikas can be stressed by climate regardless of behavioural adaptations.