Population genomic structure of a widespread, urban-dwelling mammal: The eastern grey squirrel (Sciurus carolinensis).

Urbanization is a persistent and widespread driver of global environmental change, potentially shaping evolutionary processes due to genetic drift and reduced gene flow in cities induced by habitat fragmentation and small population sizes. We tested this prediction for the eastern grey squirrel (Sciurus carolinensis), a common and conspicuous forest-dwelling rodent, by obtaining 44K SNPs using reduced representation sequencing (ddRAD) for 403 individuals sampled across the species' native range in eastern North America. We observed moderate levels of genetic diversity, low levels of inbreeding, and only a modest signal of isolation-by-distance. Clustering and migration analyses show that estimated levels of migration and genetic connectivity were higher than expected across cities and forested areas, specifically within the eastern portion of the species' range dominated by urbanization, and genetic connectivity was less than expected within the western range where the landscape is fragmented by agriculture. Landscape genetic methods revealed greater gene flow among individual squirrels in forested regions, which likely provide abundant food and shelter for squirrels. Although gene flow appears to be higher in areas with more tree cover, only slight discontinuities in gene flow suggest eastern grey squirrels have maintained connected populations across urban areas in all but the most heavily fragmented agricultural landscapes. Our results suggest urbanization shapes biological evolution in wildlife species depending strongly on the composition and habitability of the landscape matrix surrounding urban areas.

Human recreation impacts seasonal activity and occupancy of American black bears (Ursus americanus) across the anthropogenic-wildland interface.

Protected areas serve an important role in wildlife conservation, yet most wildlife occur outside these areas, subject to varying degrees of human disturbance. In the Upper Peninsula of Michigan, American black bears (Ursus americanus), a highly mobile, opportunistic species, are common despite an extensive outdoor recreation industry with the potential to affect black bear spatial and temporal activity. We investigated how environmental and anthropogenic factors influence black bear occupancy, detection, and diel activity patterns across the anthropogenic-wildland interface before and after hibernation. Using 30 camera traps deployed across a rural-wildland interface, we captured black bears at 23 camera sites (~ 77%), which exhibited co-occurrence with humans at 10 sites (~ 33%), revealing that human presence and human population density exert negative effects on black bear seasonal occupancy. Bears were more nocturnal during the hunting season, before hibernation. Human recreational activity increased ~ 38% after hibernation, but bear diurnal activity also increased ~ 36%, except when cubs were present. Our results suggest bears prioritize avoiding humans spatially, rather than temporally, except during the hunting season and when cubs are present. Understanding black bear responses to human recreation patterns and environmental variation is essential for minimizing human-mediated disturbance, and fueling conservation efforts of large, charismatic carnivores.

The gut microbiome of wild American marten in the Upper Peninsula of Michigan.

Carnivores are ecologically important and sensitive to habitat loss and anthropogenic disruption. Here we measured trophic level and gut bacterial composition as proxies of carnivore ecological status across the Upper Peninsula, Michigan, for wild American marten (Martes americana; hereafter marten). In contrast to studies that have focused on omnivorous and herbivorous species, we find that marten, like other carnivore species without a cecum, are dominated by Firmicutes (52.35%) and Proteobacteria (45.31%) but lack Bacteroidetes. Additionally, a majority of the 12 major bacterial genera (occurring at ≥1%) are known hydrogen producers, suggesting these taxa may contribute to host energy requirements through fermentative production of acetate. Our study suggests that live trapping and harvest methods yield similar marten gut microbiome data. In addition, preserving undisturbed forest likely impacts marten ecology by measurably increasing marten trophic level and altering the gut microbiome. Our study underscores the utility of the gut microbiome as a tool to monitor the ecological status of wild carnivore populations.

Correlating gut microbial membership to brown bear health metrics.

The internal mechanisms responsible for modulating physiological condition, particularly those performed by the gut microbiome (GMB), remain under-explored in wildlife. However, as latitudinal and seasonal shifts in resource availability occur, the myriad micro-ecosystem services facilitated by the GMB may be especially important to wildlife health and resilience. Here, we use brown bears (Ursus arctos) as an ecological model to quantify the relationship between wildlife body condition metrics that are commonly used to assess individual and population-level health and GMB community composition and structure. To achieve these aims, we subsampled brown bear fecal samples collected during United States National Park Service research activities at three National Parks and Preserves (Katmai, Lake Clark, and Gates of the Arctic) and extracted microbial DNA for 16S rRNA amplicon sequencing and microbial taxonomic classification. We analyzed GMB communities using alpha diversity indices, subsequently using Spearman's correlation analysis to examine relationships between alpha diversity and brown bear health metrics. We found no differences in GMB composition among bears with differing body conditions, nor any correlations between alpha diversity and body condition. Our results indicate that GMB composition reflects diverse foraging strategies while allowing brown bears to achieve similar body condition outcomes.

Intrinsic and extrinsic factors influence on an omnivore's gut microbiome.

Gut microbiomes (GMBs), complex communities of microorganisms inhabiting the gastrointestinal tracts of their hosts, perform countless micro-ecosystem services such as facilitating energy uptake and modulating immune responses. While scientists increasingly recognize the role GMBs play in host health, the role of GMBs in wildlife ecology and conservation has yet to be realized fully. Here, we use brown bears (Ursus arctos) as an ecological model to (1) characterize GMB community composition associated with location, season, and reproductive condition of a large omnivore; (2) investigate how both extrinsic and intrinsic factors influence GMB community membership and structure; and (3) quantify differences in GMB communities among different locations, seasons, sex, and reproductive conditions. To achieve these aims, we subsampled brown bear fecal samples collected during United States National Park Service research activities at three National Parks and Preserves (Katmai, Lake Clark, and Gates of the Arctic) and extracted microbial DNA for 16S rRNA amplicon sequencing and microbial taxonomic classification. We analyzed GMB communities using alpha and beta diversity indices, subsequently using linear mixed models to examine relationships between alpha diversity and extrinsic and intrinsic factors. Katmai brown bears hosted the greatest alpha diversity, whereas Gates brown bears hosted the least alpha diversity. Our results indicate that location and diet drive GMB variation, with bears hosting less phylogenetic diversity as park distance inland increases. Monitoring brown bear GMBs could enable managers to quickly detect and assess the impact of environmental perturbations on brown bear health. By integrating macro and micro-ecological perspectives we aim to inform local and landscape-level management decisions to promote long-term brown bear conservation and management.

SNAPSHOT USA 2020: A second coordinated national camera trap survey of the United States during the COVID-19 pandemic.

Managing wildlife populations in the face of global change requires regular data on the abundance and distribution of wild animals, but acquiring these over appropriate spatial scales in a sustainable way has proven challenging. Here we present the data from Snapshot USA 2020, a second annual national mammal survey of the USA. This project involved 152 scientists setting camera traps in a standardized protocol at 1485 locations across 103 arrays in 43 states for a total of 52,710 trap-nights of survey effort. Most (58) of these arrays were also sampled during the same months (September and October) in 2019, providing a direct comparison of animal populations in 2 years that includes data from both during and before the COVID-19 pandemic. All data were managed by the eMammal system, with all species identifications checked by at least two reviewers. In total, we recorded 117,415 detections of 78 species of wild mammals, 9236 detections of at least 43 species of birds, 15,851 detections of six domestic animals and 23,825 detections of humans or their vehicles. Spatial differences across arrays explained more variation in the relative abundance than temporal variation across years for all 38 species modeled, although there are examples of significant site-level differences among years for many species. Temporal results show how species allocate their time and can be used to study species interactions, including between humans and wildlife. These data provide a snapshot of the mammal community of the USA for 2020 and will be useful for exploring the drivers of spatial and temporal changes in relative abundance and distribution, and the impacts of species interactions on daily activity patterns. There are no copyright restrictions, and please cite this paper when using these data, or a subset of these data, for publication.

SNAPSHOT USA 2019: a coordinated national camera trap survey of the United States.

With the accelerating pace of global change, it is imperative that we obtain rapid inventories of the status and distribution of wildlife for ecological inferences and conservation planning. To address this challenge, we launched the SNAPSHOT USA project, a collaborative survey of terrestrial wildlife populations using camera traps across the United States. For our first annual survey, we compiled data across all 50 states during a 14-week period (17 August-24 November of 2019). We sampled wildlife at 1,509 camera trap sites from 110 camera trap arrays covering 12 different ecoregions across four development zones. This effort resulted in 166,036 unique detections of 83 species of mammals and 17 species of birds. All images were processed through the Smithsonian's eMammal camera trap data repository and included an expert review phase to ensure taxonomic accuracy of data, resulting in each picture being reviewed at least twice. The results represent a timely and standardized camera trap survey of the United States. All of the 2019 survey data are made available herein. We are currently repeating surveys in fall 2020, opening up the opportunity to other institutions and cooperators to expand coverage of all the urban-wild gradients and ecophysiographic regions of the country. Future data will be available as the database is updated at eMammal.si.edu/snapshot-usa, as will future data paper submissions. These data will be useful for local and macroecological research including the examination of community assembly, effects of environmental and anthropogenic landscape variables, effects of fragmentation and extinction debt dynamics, as well as species-specific population dynamics and conservation action plans. There are no copyright restrictions; please cite this paper when using the data for publication.

Wild black bears harbor simple gut microbial communities with little difference between the jejunum and colon.

The gut microbiome (GMB), comprising the commensal microbial communities located in the gastrointestinal tract, has co-evolved in mammals to perform countless micro-ecosystem services to facilitate physiological functions. Because of the complex inter-relationship between mammals and their gut microbes, the number of studies addressing the role of the GMB on mammalian health is almost exclusively limited to human studies and model organisms. Furthermore, much of our knowledge of wildlife-GMB relationships is based on studies of colonic GMB communities derived from the feces of captive specimens, leaving our understanding of the GMB in wildlife limited. To better understand wildlife-GMB relationships, we engaged hunters as citizen scientists to collect biological samples from legally harvested black bears (Ursus americanus) and used 16S rRNA gene amplicon sequencing to characterize wild black bear GMB communities in the colon and jejunum, two functionally distinct regions of the gastrointestinal tract. We determined that the jejunum and colon of black bears do not harbor significantly different GMB communities: both gastrointestinal sites were dominated by Firmicutes and Proteobacteria. However, a number of bacteria were differentially enriched in each site, with the colon harboring twice as many enriched taxa, primarily from closely related lineages.

Noninvasive measures of physiological stress are confounded by exposure.

Glucocorticoids and glucocorticoid metabolites are increasingly used to index physiological stress in wildlife. Although feces is often abundant and can be collected noninvasively, exposure to biotic and abiotic elements may influence fecal glucocorticoid metabolite (FGM) concentrations, leading to inaccurate conclusions regarding wildlife physiological stress. Using captive snowshoe hares (Lepus americanus) and simulated environmental conditions, we evaluated how different realistic field conditions and temporal sampling constraints might influence FGM concentrations using an 11-oxoetiocholanolone-enzyme immunoassay. We quantified how fecal pellet age (i.e., 0-6 days), variable summer temperatures, and precipitation affected FGM concentrations. Fecal pellet age had a strong effect on FGM concentrations (ßAge = 0.395, s.d. = 0.085; ß2Age = -0.061, s.d. = 0.012), which were lowest at the beginning and end of our exposure period (e.g., meanday6 = 37.7 ng/mg) and typically highest in the middle (meanday3 = 51.8 ng/mg). The effect of fecal pellet age on FGM concentrations varied across treatments with warm-dry and cool-wet conditions resulting in more variable FGM concentrations relative to control samples. Given the confounding effects of exposure and environmental conditions, if fresh fecal pellet collection is not an option, we encourage researchers to develop a temporally consistent sampling protocol to ensure all samples are exposed to similar environmental conditions.

Adaptive introgression underlies polymorphic seasonal camouflage in snowshoe hares.

Snowshoe hares (Lepus americanus) maintain seasonal camouflage by molting to a white winter coat, but some hares remain brown during the winter in regions with low snow cover. We show that cis-regulatory variation controlling seasonal expression of the Agouti gene underlies this adaptive winter camouflage polymorphism. Genetic variation at Agouti clustered by winter coat color across multiple hare and jackrabbit species, revealing a history of recurrent interspecific gene flow. Brown winter coats in snowshoe hares likely originated from an introgressed black-tailed jackrabbit allele that has swept to high frequency in mild winter environments. These discoveries show that introgression of genetic variants that underlie key ecological traits can seed past and ongoing adaptation to rapidly changing environments.

Winter color polymorphisms identify global hot spots for evolutionary rescue from climate change.

Maintenance of biodiversity in a rapidly changing climate will depend on the efficacy of evolutionary rescue, whereby population declines due to abrupt environmental change are reversed by shifts in genetically driven adaptive traits. However, a lack of traits known to be under direct selection by anthropogenic climate change has limited the incorporation of evolutionary processes into global conservation efforts. In 21 vertebrate species, some individuals undergo a seasonal color molt from summer brown to winter white as camouflage against snow, whereas other individuals remain brown. Seasonal snow duration is decreasing globally, and fitness is lower for winter white animals on snowless backgrounds. Based on 2713 georeferenced samples of known winter coat color-from eight species across trophic levels-we identify environmentally driven clinal gradients in winter coat color, including polymorphic zones where winter brown and white morphs co-occur. These polymorphic zones, underrepresented by existing global protected area networks, indicate hot spots for evolutionary rescue in a changing climate.

The influence of sulfur and hair growth on stable isotope diet estimates for grizzly bears.

Stable isotope ratios of grizzly bear (Ursus arctos) guard hair collected from bears on the lower Stikine River, British Columbia (BC) were analyzed to: 1) test whether measuring δ34S values improved the precision of the salmon (Oncorhynchus spp.) diet fraction estimate relative to δ15N as is conventionally done, 2) investigate whether measuring δ34S values improves the separation of diet contributions of moose (Alces alces), marmot (Marmota caligata), and mountain goat (Oreamnos americanus) and, 3) examine the relationship between collection date and length of hair and stable isotope values. Variation in isotope signatures among hair samples from the same bear and year were not trivial. The addition of δ34S values to mixing models used to estimate diet fractions generated small improvement in the precision of salmon and terrestrial prey diet fractions. Although the δ34S value for salmon is precise and appears general among species and areas, sulfur ratios were strongly correlated with nitrogen ratios and therefore added little new information to the mixing model regarding the consumption of salmon. Mean δ34S values for the three terrestrial herbivores of interest were similar and imprecise, so these data also added little new information to the mixing model. The addition of sulfur data did confirm that at least some bears in this system ate marmots during summer and fall. We show that there are bears with short hair that assimilate >20% salmon in their diet and bears with longer hair that eat no salmon living within a few kilometers of one another in a coastal ecosystem. Grizzly bears are thought to re-grow hair between June and October however our analysis of sectioned hair suggested at least some hairs begin growing in July or August, not June and, that hair of wild bears may grow faster than observed in captive bears. Our hair samples may have been from the year of sampling or the previous year because samples were collected in summer when bears were growing new hair. The salmon diet fraction increased with later hair collection dates, as expected if samples were from the year of sampling because salmon began to arrive in mid-summer. Bears that ate salmon had shorter hair and δ15N and δ34S values declined with hair length, also suggesting some hair samples were grown the year of sampling. To be sure to capture an entire hair growth period, samples must be collected in late fall. Early spring samples are also likely to be from the previous year but the date when hair begins to grow appears to vary. Choosing the longest hair available should increase the chance the hair was grown during the previous year and, maximize the period for which diet is measured.

Within-sample variation in snowshoe hare faecal glucocorticoid metabolite measurements.

Faecal glucocorticoid metabolite (FGM) concentrations are used increasingly as a non-invasive measure to index physiological stress experienced by diverse taxa. However, FGM may not be evenly distributed throughout a faecal mass or faecal pellet group. Moreover, within-sample variation in FGM measurements associated with different sampling and/or processing techniques is rarely reported despite potentially having important implications for inferring stress levels in free-ranging wildlife. Using a captive collection of snowshoe hares (Lepus americanus), we (i) assessed repeatability of FGM measurements (i.e. precision) from two processing techniques (measurements derived from dividing whole pellet groups into equal proportions prior to processing [G1], measurements from subsamples derived from thoroughly homogenized whole pellet groups [G2]) and (ii) conducted a power analysis to estimate sample-size requirements for detecting statistically significant differences in FGM concentrations at a population level. Our results indicate that the mean percent coefficient of variation (%CV) for within-sample FGM variation was slightly higher for G1 (%CV = 35, range 13.45-65.37) than for G2 (%CV = 23, range 7.26-47.94), though not statistically significant (two sample t-test, n = 8, t = 1.57, P = 0.16). Thus, FGM is relatively evenly distributed within snowshoe hare faecal pellet groups. However, subsampling from homogenized whole pellet groups may be more appropriate when the sampling time frame is less controlled (e.g. multiple defecation events) because a subsample derived from a homogenized whole pellet group might be more representative of the animal's 'average' physiological state compared to FGM concentrations derived from a few haphazardly selected faecal pellets. Power analysis results demonstrated the importance of a priori consideration of sample sizes. Relatively small effect sizes (e.g. <20%) may require sampling that is logistically and/or cost prohibitive. Yet for many situations of ecological or conservation interest, treatment effects may be substantial (>25%) and thus moderate sample sizes may be sufficient for testing research hypotheses regarding changes FGM concentrations.

Sex, Diet, and the Social Environment: Factors Influencing Hair Cortisol Concentration in Free-Ranging Black Bears (Ursus americanus).

Increasingly, measures of glucocorticoid levels (e.g., cortisol), key components of the neuroendocrine stress axis, are being used to measure past hypothalamic-pituitary-adrenal (HPA) activity to index psychological and physiological stress exhibited by wildlife for assessing individual and population-level well-being. However, many intrinsic and extrinsic factors affect HPA activity in animals. Using American black bears (Ursus americanus; n = 116) as an ecological model and hair cortisol concentration (HCC) as an integrative measure of past HPA activity, we evaluated the influence of diet, sex and the social environment on black bear HCC in a free-ranging population that spanned adjoining ecoregions with differing densities of potential conspecific and heterospecific competitors. HCC varied by sex, with female HCC ranging from 0.6 to 10.7 pg/mg (median = 4.5 ± 1.2 mean absolute deviation [MAD]) and male HCC ranging from 0.5 to 35.1 pg/mg (median = 6.2 ± 2.6 MAD). We also observed a three-way interaction among sex, δ14C and ecoregion, which may indicate that some differences in HCC between female and male black bears results from variability in the nutritional needs of larger-bodied males relative to smaller-bodied females, slight differences in food resources use between ecoregions as well as sex-based differences regarding the social environment. Once we understand what drives sex-specific differences in HCC, HCC may aid our understanding of the physiological responses by bears and other wildlife to diverse environmental challenges.