Natural rewilding of the Great Basin: Genetic consequences of recolonization by black bears (Ursus americanus )

Spatial patterns of reproduction suggest marginal habitat limits continued range expansion of black bears at a forest-desert ecotone

Investigating spatial patterns of animal occupancy and reproduction in peripheral populations can provide insight into factors that form species range boundaries. Following historical extirpation, American black bears (Ursus americanus) recolonized the western Great Basin in Nevada from the Sierra Nevada during the late 1900s. This range expansion, however, has not continued further into the Great Basin despite the presence of additional habitat. We aimed to quantify whether reduced reproduction toward the range edge contributes to this range boundary. We analyzed black bear detections from 100 camera traps deployed across black bear distribution in western Nevada using a multistate occupancy model that quantified the probability of occupancy and reproduction (i.e., female bears with cubs occupancy) in relation to changes in habitat type and habitat amount toward the range boundary. We detected a strong effect of habitat amount and habitat type on the probability of black bear occupancy and reproduction. At similar levels of landscape-scale habitat amount (e.g., 50%), estimated probability of occupancy for adult bears in piñon-juniper woodlands near the range boundary was 0.39, compared to ~1.0 in Sierra Nevada mixed-conifer forest (i.e., core habitat). Furthermore, estimated probability of cub occupancy, conditional on adult bear occupancy, in landscapes with 50% habitat was 0.32 in Great Basin piñon-juniper woodlands, compared to 0.92 in Sierra Nevada mixed-conifer forest. Black bear range in the western Great Basin conforms to the center-periphery hypothesis, with piñon-juniper woodland at the range edge supporting ecologically marginal habitat for the species compared to habitat in the Sierra Nevada. Further geographic expansion of black bears in the Great Basin may be limited by lower occupancy of reproducing females in piñon-juniper woodland. Center-periphery range dynamics may be common in large carnivore species, as their dispersal ability allows them to colonize low-quality habitat near range edges.

Contrasting genomic consequences of anthropogenic reintroduction and natural recolonization in high-arctic wild reindeer

Anthropogenic reintroduction can supplement natural recolonization in reestablishing a species' distribution and abundance. However, both reintroductions and recolonizations can give rise to founder effects that reduce genetic diversity and increase inbreeding, potentially causing the accumulation of genetic load and reduced fitness. Most current populations of the endemic high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) originate from recent reintroductions or recolonizations following regional extirpations due to past overharvesting. We investigated and compared the genomic consequences of these two paths to reestablishment using whole-genome shotgun sequencing of 100 Svalbard reindeer across their range. We found little admixture between reintroduced and natural populations. Two reintroduced populations, each founded by 12 individuals around four decades (i.e. 8 reindeer generations) ago, formed two distinct genetic clusters. Compared to the source population, these populations showed only small decreases in genome-wide heterozygosity and increases in inbreeding and lengths of runs of homozygosity. In contrast, the two naturally recolonized populations without admixture possessed much lower heterozygosity, higher inbreeding and longer runs of homozygosity, possibly caused by serial population founder effects and/or fewer or more genetically related founders than in the reintroduction events. Naturally recolonized populations can thus be more vulnerable to the accumulation of genetic load than reintroduced populations. This suggests that in some organisms even small-scale reintroduction programs based on genetically diverse source populations can be more effective than natural recolonization in establishing genetically diverse populations. These findings warrant particular attention in the conservation and management of populations and species threatened by habitat fragmentation and loss.

Genetic architecture and evolution of color variation in American black bears

Color variation is a frequent evolutionary substrate for camouflage in small mammals, but the underlying genetics and evolutionary forces that drive color variation in natural populations of large mammals are mostly unexplained. The American black bear, Ursus americanus (U. americanus), exhibits a range of colors including the cinnamon morph, which has a similar color to the brown bear, U. arctos, and is found at high frequency in the American southwest. Reflectance and chemical melanin measurements showed little distinction between U. arctos and cinnamon U. americanus individuals. We used a genome-wide association for hair color as a quantitative trait in 151 U. americanus individuals and identified a single major locus (p < 10-13). Additional genomic and functional studies identified a missense alteration (R153C) in Tyrosinase-related protein 1 (TYRP1) that likely affects binding of the zinc cofactor, impairs protein localization, and results in decreased pigment production. Population genetic analyses and demographic modeling indicated that the R153C variant arose 9.36 kya in a southwestern population where it likely provided a selective advantage, spreading both northwards and eastwards by gene flow. A different TYRP1 allele, R114C, contributes to the characteristic brown color of U. arctos but is not fixed across the range.

Pleistocene-Holocene vicariance, not Anthropocene landscape change, explains the genetic structure of American black bear (Ursus americanus) populations in the American Southwest and northern Mexico

The phylogeography of the American black bear (Ursus americanus) is characterized by isolation into glacial refugia, followed by population expansion and genetic admixture. Anthropogenic activities, including overharvest, habitat loss, and transportation infrastructure, have also influenced their landscape genetic structure. We describe the genetic structure of the American black bear in the American Southwest and northern Mexico and investigate how prehistoric and contemporary forces shaped genetic structure and influenced gene flow. Using a suite of microsatellites and a sample of 550 bears, we identified 14 subpopulations organized hierarchically following the distribution of ecoregions and mountain ranges containing black bear habitat. The pattern of subdivision we observed is more likely a product of postglacial habitat fragmentation during the Pleistocene and Holocene, rather than a consequence of contemporary anthropogenic barriers to movement during the Anthropocene. We used linear mixed-effects models to quantify the relationship between landscape resistance and genetic distance among individuals, which indicated that both isolation by resistance and geographic distance govern gene flow. Gene flow was highest among subpopulations occupying large tracts of contiguous habitat, was reduced among subpopulations in the Madrean Sky Island Archipelago, where montane habitat exists within a lowland matrix of arid lands, and was essentially nonexistent between two isolated subpopulations. We found significant asymmetric gene flow supporting the hypothesis that bears expanded northward from a Pleistocene refugium located in the American Southwest and northern Mexico and that major highways were not yet affecting gene flow. The potential vulnerability of the species to climate change, transportation infrastructure, and the US-Mexico border wall highlights conservation challenges and opportunities for binational collaboration.

Habitat connectivity and resource selection in an expanding bobcat (Lynx rufus) population

Terrestrial carnivores are among the most imperiled species worldwide, yet some species are resilient and are recovering in human-dominated landscapes after decades or centuries of absence. Bobcat (Lynx rufus) populations were extirpated from much of Midwestern US in the mid-1800's, and are currently expanding and recolonizing their former range. In this study, we investigated multi-scale habitat selection for Ohio's expanding bobcat population, and examined habitat connectivity in order to evaluate the conduits for dispersal statewide. We used citizen observations collected between 1978 and 2019 and logistic regression to evaluate population-level habitat selection, and GPS telemetry data for 20 individuals collected between 2012 and 2014 and a distribution-weighted exponential Resource Selection Function to evaluate individual-level habitat selection within home ranges. At the population level, bobcats selected for higher amounts of forest and pasture (at a 50 km² scale) and herbaceous vegetation (at 15-50 50 km² scales), thus overall heterogeneous forested habitat. At individual (home range) level, bobcats selected for forested habitats with low road density and farther away from high traffic roads; they also showed weak selection for open habitat at the home range level. Male home ranges were significantly greater than female home ranges. Lastly, we used the population-level spatial outputs (i.e. habitat suitability map) to parameterize habitat connectivity models using circuit theory in the program Circuitscape. We tested three relationships between habitat suitability and resistance to movement and used a subset of data on potential dispersing individuals to evaluate which relationship performed best. All three relationships performed almost equally well, and we calculated a weighted averaged connectivity map as our final map. Habitat was highly permeable to movements between core areas of two genetically distinct subpopulations located in southeastern Ohio. We also identified potential dispersal corridors from the core areas to other regions of Ohio dominated by agriculture and suburban development via forested riparian corridors. Overall, our analysis offers new information on habitat selection and connectivity in a rebounding felid population and offers important ecological information for wildlife management strategies. We recommend that the suitability and connectivity models should be periodically updated until the population reaches an equilibrium, and be integrated with data from neighboring states for a comprehensive assessment of a conservation success story.

An apex carnivore's life history mediates a predator cascade

Apex predators can shape communities via cascading top-down effects, but the degree to which such effects depend on predator life history traits is largely unknown. Within carnivore guilds, complex hierarchies of dominance facilitate coexistence, whereby subordinate species avoid dominant counterparts by partitioning space, time, or both. We investigated whether a major life history trait (hibernation) in an apex carnivore (black bears Ursus americanus) mediated its top-down effects on the spatio-temporal dynamics of three sympatric mesocarnivore species (coyotes Canis latrans, bobcats Lynx rufus, and gray foxes Urocyon cinereoargenteus) across a 15,000 km² landscape in the western USA. We compared top-down, bottom-up, and environmental effects on these mesocarnivores using an integrated modeling approach. Black bears exerted top-down effects that varied as a function of hibernation and were stronger than bottom-up or environmental impacts. High black bear activity in summer and fall appeared to buffer the most subordinate mesocarnivore (gray foxes) from competition with dominant mesocarnivores (coyotes and bobcats), which were in turn released by black bear hibernation in winter and early spring. The mesocarnivore responses occurred in space (i.e., altered occupancy and site visitation intensity) rather than time (i.e., diel activity patterns unaffected). These results suggest that the spatio-temporal dynamics of mesocarnivores in this system were principally shaped by a spatial predator cascade of interference competition mediated by black bear hibernation. Thus, certain life history traits of apex predators might facilitate coexistence among competing species over broad time scales, with complex implications for lower trophic levels.

Recolonizing carnivores: Is cougar predation behaviorally mediated by bears?

Conservation and management efforts have resulted in population increases and range expansions for some apex predators, potentially changing trophic cascades and foraging behavior. Changes in sympatric carnivore and dominant scavenger populations provide opportunities to assess how carnivores affect one another. Cougars (Puma concolor) were the apex predator in the Great Basin of Nevada, USA, for over 80 years. Black bears (Ursus americanus) have recently recolonized the area and are known to heavily scavenge on cougar kills. To evaluate the impacts of sympatric, recolonizing bears on cougar foraging behavior in the Great Basin, we investigated kill sites of 31 cougars between 2009 and 2017 across a range of bear densities. We modeled the variation in feeding bout duration (number of nights spent feeding on a prey item) and the proportion of primary prey, mule deer (Odocoileus hemionus), in cougar diets using mixed-effects models. We found that feeding bout duration was driven primarily by the size of the prey item being consumed, local bear density, and the presence of dependent kittens. The proportion of mule deer in cougar diet across all study areas declined over time, was lower for male cougars, increased with the presence of dependent kittens, and increased with higher bear densities. In sites with feral horses (Equus ferus), a novel large prey, cougar consumption of feral horses increased over time. Our results suggest that higher bear densities over time may reduce cougar feeding bout durations and influence the prey selection trade-off for cougars when alternative, but more dangerous, large prey are available. Shifts in foraging behavior in multicarnivore systems can have cascading effects on prey selection. This study highlights the importance of measuring the impacts of sympatric apex predators and dominant scavengers on a shared resource base, providing a foundation for monitoring dynamic multipredator/scavenger systems.

Genetic differentiation and asymmetric gene flow among Carpathian brown bear (Ursus arctos) populations-Implications for conservation of transboundary populations

The abundance and distribution of large carnivores in Europe have been historically reduced. Their recovery requires multilevel coordination, especially regarding transboundary populations. Here, we apply nuclear and mitochondrial genetic markers to test for admixture level and its impact on population genetic structure of contemporary brown bears (Ursus arctos) from the Eastern, Southern, and Western Carpathians. Carpathian Mountains (Europe). Nearly 400 noninvasive brown bear DNA samples from the Western (Poland) and Eastern Carpathians (Bieszczady Mountains in Poland, Slovakia, Ukraine) were collected. Together with DNA isolates from Slovakia and Romania, they were analyzed using the set of eight microsatellite loci and two mtDNA regions (control region and cytochrome b). A set of 113 individuals with complete genotypes was used to investigate genetic differentiation across national boundaries, genetic structuring within and between populations, and movement between populations. Transboundary brown bear subpopulations (Slovakia and Poland) did not show significant internal genetic structure, and thus were treated as cohesive units. All brown bears from the Western Carpathians carried mitochondrial haplotypes from the Eastern lineage, while the Western lineage prevailed in the brown bears from the Bieszczady Mountains. Despite similar levels of microsatellite variability, we documented significant differentiation among the studied populations for nuclear markers and mtDNA. We also detected male-biased and asymmetrical movement into the Bieszczady Mountains population from the Western Carpathians. Our findings suggest initial colonization of the Western Carpathians by brown bears possessing mtDNA from the Eastern lineage. Genetic structuring among populations at microsatellite loci could be a result of human-mediated alterations. Detected asymmetric gene flow suggests ongoing expansion from more abundant populations into the Bieszczady Mountains and thus supports a metapopulation model. The knowledge concerning this complex pattern can be implemented in a joint Carpathian brown bear management plan that should allow population mixing by dispersing males.