skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, November 15 until 2:00 AM ET on Saturday, November 16 due to maintenance. We apologize for the inconvenience.


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
NSF-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. 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. Abstract

    Some animal species are responding to climate change by altering the timing of events like mating and migration. Such behavioral plasticity can be adaptive, but it is not always. Polar bears (Ursus maritimus) from the southern Beaufort Sea subpopulation have mostly remained on ice year‐round, but as the climate warms and summer sea ice declines, a growing proportion of the subpopulation is summering ashore. The triggers of this novel behavior are not well understood. Our study uses a parametric time‐to‐event model to test whether biological and/or time‐varying environmental variables thought to influence polar bear movement and habitat selection also drive decisions to swim ashore. We quantified the time polar bears spent occupying offshore sea ice of varying ice concentrations. We evaluated variations in the ordinal date bears moved to land with respect to local environmental conditions such as sea ice concentration and wind across 10 years (2005–2015). Results from our study suggest that storm events (i.e., sustained high wind speeds) may force polar bears from severely degraded ice habitat and catalyze seasonal movements to land. Unlike polar bears long adapted to complete summer ice melt, southern Beaufort Sea bears that summer ashore appear more tolerant of poor‐quality sea ice habitat and are less willing to abandon it. Our findings provide a window into emergent, climatically mediated behavior in an Arctic marine mammal vulnerable to rapid habitat decline.

     
    more » « less
  3. Abstract

    Arctic climate change poses serious threats to polar bears (Ursus maritimus) as reduced sea ice makes seal prey inaccessible and marine ecosystems undergo bottom‐up reorganization. Polar bears’ elongated skulls and reduced molar dentition, as compared to their sister species the grizzly bear (Ursus arctos), are adaptations associated with hunting seals on sea ice and a soft, lipid‐rich diet of blubber and meat. With significant declines in sea ice, it is unclear if and how polar bears may be altering their diets. Clarifying polar bear dietary responses to changing climates, both today and in the past, is critical to proper conservation and management of this apex predator. This is particularly important when a dietary strategy may be maladaptive. Here, we test the hypothesis that hard‐food consumption (i.e., less preferred foods including bone), inferred from dental microwear texture analysis, increased with Arctic warming. We find that polar bears demonstrate a conserved absence of hard‐object feeding in Alaska through time (including approximately 1000 years ago), until the 21st century, consistent with a highly conserved and specialized diet of soft blubber and flesh. Notably, our results also suggest that some 21st‐century polar bears may be consuming harder foods (e.g., increased carcass utilization, terrestrial foods including garbage), despite having skulls and metabolisms poorly suited for such a diet. Prior to the 21st century, only polar bears with larger mandibles demonstrated increased hard‐object feeding, though to a much lower degree than closely related grizzly bears which regularly consume mechanically challenging foods. Polar bears, being morphologically specialized, have biomechanical constraints which may limit their ability to consume mechanically challenging diets, with dietary shifts occurring only under the most extreme scenarios. Collectively, the highly specialized diets and cranial morphology of polar bears may severely limit their ability to adapt to a warming Arctic.

     
    more » « less
  4. 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
  5. Abstract

    Recent reductions in thickness and extent have increased drift rates of Arctic sea ice. Increased ice drift could significantly affect the movements and the energy balance of polar bears (Ursus maritimus) which forage, nearly exclusively, on this substrate. We used radio‐tracking and ice drift data to quantify the influence of increased drift on bear movements, and we modeled the consequences for energy demands of adult females in the Beaufort and Chukchi seas during two periods with different sea ice characteristics. Westward and northward drift of the sea ice used by polar bears in both regions increased between 1987–1998 and 1999–2013. To remain within their home ranges, polar bears responded to the higher westward ice drift with greater eastward movements, while their movements north in the spring and south in fall were frequently aided by ice motion. To compensate for more rapid westward ice drift in recent years, polar bears covered greater daily distances either by increasing their time spent active (7.6%–9.6%) or by increasing their travel speed (8.5%–8.9%). This increased their calculated annual energy expenditure by 1.8%–3.6% (depending on region and reproductive status), a cost that could be met by capturing an additional 1–3 seals/year. Polar bears selected similar habitats in both periods, indicating that faster drift did not alter habitat preferences. Compounding reduced foraging opportunities that result from habitat loss; changes in ice drift, and associated activity increases, likely exacerbate the physiological stress experienced by polar bears in a warming Arctic.

     
    more » « less