skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Insights into bear evolution from a Pleistocene polar bear genome
The polar bear ( Ursus maritimus ) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears ( Ursus arctos ). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear’s lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.  more » « less
Award ID(s):
1854550
PAR ID:
10411812
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
119
Issue:
24
ISSN:
0027-8424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Molecular Biology and Evolution)
    Yoder, Anne (Ed.)
    Abstract Over the past 15 years, the D-statistic, a four-taxon test for organismal admixture (hybridization, or introgression) which incorporates single nucleotide polymorphism data with allelic patterns ABBA and BABA, has seen considerable use. This statistic seeks to discern significant deviation from either a given species tree assumption, or from the balanced incomplete lineage sorting that could otherwise defy this species tree. However, while the D-statistic can successfully discriminate admixture from incomplete lineage sorting, it is not a simple matter to determine the directionality of admixture using only four-leaf tree models. As such, methods have been developed that use 5 leaves to evaluate admixture. Among these, the DFOIL method, which tests allelic patterns on the “symmetric” tree S = (((1,2),(3,4)),5), succeeds in finding admixture direction for many five-taxon examples. However, DFOIL does not make full use of all symmetry, nor can DFOIL function properly when ancient samples are included because of the reliance on singleton patterns (such as BAAAA and ABAAA). Here, we take inspiration from DFOIL to develop a new and completely general family of five-leaf admixture tests, dubbed Δ-statistics, that can either incorporate or exclude the singleton allelic patterns depending on individual taxon and age sampling choices. We describe two new shapes that are also fully testable, namely the “asymmetric” tree A = ((((1,2),3),4),5) and the “quasisymmetric” tree Q = (((1,2),3),(4,5)), which can considerably supplement the “symmetric“ S = (((1,2),(3,4)),5) model used by DFOIL. We demonstrate the consistency of Δ-statistics under various simulated scenarios, and provide empirical examples using data from black, brown and polar bears, the latter also including two ancient polar bear samples from previous studies. Recently DFOIL and one of these ancient samples was used to argue for a dominant polar bear → brown bear introgression direction. However, we find, using both this ancient polar bear and our own, that by far the strongest signal using both DFOIL and Δ-statistics on tree S is actually bidirectional gene flow of indistinguishable direction. Further experiments on trees A and Q instead highlight what were likely two phases of admixture: one with stronger brown bear → polar bear introgression in ancient times, and a more recent phase with predominant polar bear → brown bear directionality. Code and documentation available at https://github.com/KalleLeppala/Delta-statistics. 
    more » « less
  2. (, NSF Arctic Data Center)
    This dataset contains information about bear specimens (Ursus spp.) identified from two archaeological sites in Unalaska, Alaska: Margaret Bay (UNL-048, 4700 Before Present (BP)) and Amaknak Bridge (UNL-050, 2500 BP). The material is accessioned at the Museum of the Aleutians in Unalaska. Zooarchaeological data included are three dimensional scans (OBJ, MTL, and JPG files for rendering in MeshLab software) and digital photographs. 3D scans were obtained by using a DAVID SLS-3 HD Structured White Light 3D Scanner at the University of Oklahoma Laboratories of Molecular Anthropology and Microbiome Research (LMAMR). Photographs were taken with a Nikon D3400 24.2-megapixel DX format DSLR camera. Zooarchaeological analyses were performed to investigate how bears came to be in Unalaska Bay, a region where bears do not currently live. It has been suggested that this assemblage is evidence of Neoglacial expansion of sea ice in the region and subsequent range expansion of polar bears. Our goals here were to assess whether the bears can be distinguished to species and determine whether these animals were harvested locally and represent a range expansion. The results suggest that bears are very difficult to distinguish using morphological characteristics alone: we argue there is likely a mix of brown and polar bear in this small assemblage, but that morphological analyses alone are inadequate for reconstructing bear distribution in this context. However, the age profiles and butchery patterns do suggest that bears were harvested locally, and we contend that expanding Neoglacial sea ice facilitated their presence around Unalaska Island. 
    more » « less
  3. ; ; ; ; ; ; ; (, Genome Biology and Evolution)
    Fraser, Bonnie (Ed.)
    Abstract The brown bear (Ursus arctos) is the second largest and most widespread extant terrestrial carnivore on Earth and has recently emerged as a medical model for human metabolic diseases. Here, we report a fully phased chromosome-level assembly of a male North American brown bear built by combining Pacific Biosciences (PacBio) HiFi data and publicly available Hi-C data. The final genome size is 2.47 Gigabases (Gb) with a scaffold and contig N50 length of 70.08 and 43.94 Megabases (Mb), respectively. Benchmarking Universal Single-Copy Ortholog (BUSCO) analysis revealed that 94.5% of single copy orthologs from Mammalia were present in the genome (the highest of any ursid genome to date). Repetitive elements accounted for 44.48% of the genome and a total of 20,480 protein coding genes were identified. Based on whole genome alignment to the polar bear, the brown bear is highly syntenic with the polar bear, and our phylogenetic analysis of 7,246 single-copy orthologs supports the currently proposed species tree for Ursidae. This highly contiguous genome assembly will support future research on both the evolutionary history of the bear family and the physiological mechanisms behind hibernation, the latter of which has broad medical implications. 
    more » « less
  4. ; ; ; (, Quaternary Science Reviews)
    Polar bear (Ursus maritimus) is the apex predator of the Arctic, largely dependent on sea-ice. The expected disappearance of the ice cover of the Arctic seas by the mid 21st century is predicted to cause a dramatic decrease in the global range and population size of the species. To place this scenario against the backdrop of past distribution changes and their causes, we use a fossil dataset to investigate the polar bear’s past distribution dynamics during the Late Glacial and the Holocene. Fossil results indicate that during the last deglaciation, polar bears were present at the southwestern margin of the Scandinavian Ice Sheet, surviving until the earliest Holocene. There are no Arctic polar bear findings from 8,000-6,000 years ago (8-6 ka), the Holocene’s warmest period. However, fossils that date from 8-9 ka and 5-6 ka suggest that the species likely survived this period in cold refugia located near the East Siberian Sea, northern Greenland and the Canadian Archipelago. Polar bear range expansion is documented by an increase in fossils during the last 4,000 years in tandem with cooling climate and expanding Arctic sea ice. The results document changes in polar bear’s distribution in response to Late Glacial and Holocene Arctic temperature and sea ice trends. 
    more » « less
  5. ; ; ; ; ; ; ; (, Ecosphere)
    Abstract Wildlife ecologists throughout the world strive to monitor trends in population abundance to help manage wildlife populations and conserve species at risk. Spatial capture–recapture studies are the gold standard for monitoring density, yet they can be difficult to apply because researchers must be able to distinguish all detected individuals. Spatial mark–resight (SMR) models only require a subset of the population to be marked and identifiable. Recent advances in SMR models with radio‐collared animals required a two‐staged analysis. We developed a one‐stage generalized SMR (gSMR) model that used detection histories of marked and unmarked animals in a single analysis. We used simulations to assess the performance of one‐ and two‐stage gSMR models. We then applied the one‐stage gSMR with telemetry and remote camera data to estimate grizzly bear (Ursus arctos) abundance from 2012 to 2023 within the Canadian Rocky Mountains. We estimated abundance trends for the population and reproductive females (females with cubs of the year). Simulations suggest that one‐ and two‐stage models performed equally well. One‐stage models are more dependable as they use exact likelihoods, whereas two‐stage models have shorter computation times for large data sets. Both methods had >95% credible interval coverage and minimal bias. Increasing the number of marked animals increased the accuracy and precision of abundance estimates, and ≥10 marked animals were required to obtain coefficients of variation <20% in most scenarios. The grizzly bear population increased slightly (growth rate λmean = 1.02) to a 2023 density of 10.4 grizzly bears/1000 km2. Reproductive female abundance had high interannual variability and increased to 1.0 bears/1000 km2. Population density was highest within protected areas, within high‐quality habitat and far from paved roads. The density of activity centers declined near paved roads over time. Mechanisms of decline may have included direct mortality and shifting activity centers to avoid human activity. Our study demonstrates the influence of human activity on localized density and the importance of protected areas for carnivore conservation. Finally, our study highlights the widespread utility of remote camera and telemetry‐based SMR models for monitoring spatiotemporal trends in abundance. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email