Glacial ice supports a distinct and undocumented polar bear subpopulation persisting in late 21st-century sea-ice conditions

Patterns and temporal trends in canine breakage and scarring in Western Hudson Bay polar bears (Ursus maritimus)

Canines are used by carnivores for prey capture and social interactions but are often damaged. The highly carnivorous polar bear (Ursus maritimus) has a female defence polygyny mating system where males compete for access to females and injuries to males, such as broken canines and cuts, are common. The Western Hudson Bay polar bear subpopulation has declined in abundance in recent decades and shifted from a female-biased to a male-biased adult sex ratio, which may have affected their mating system. We hypothesize that if changes in subpopulation structure have affected the mating system, then canine breakage and scarring may have changed over time. We assessed age- and sex-specific occurrences of canine breakage and scarring in 3493 individuals between 1981-2023 using non-parametric statistical analyses and linear mixed effect models. We found age- and sex-related differences in mean values of breakage and scarring. These injury occurrences increased with age in both sexes and males had greater amounts of both breakage and scarring compared to females. As the only main effect, sampling year was significant and indicated increasing breakage in both sexes over time; however, the top breakage model did not include year, indicating it was not as strong a predictor of breakage and scarring as age or sex. Age, sex, and year were all in the top model for predicting maximum scarring probabilities. We found some evidence that injuries changed over time, suggesting there could be changes to intraspecific interactions, but additional monitoring is needed.

Heterochronous mitogenomes shed light on the Holocene history of the Scandinavian brown bear

Following glacial retreat after the last ice age, brown bears (Ursus arctos) recolonised Scandinavia. Previous research based on mitochondrial markers suggests that bears recolonised from both the north and the south, with a contact zone in central Scandinavia. More recently, the Scandinavian brown bear was subjected to a strong population decline with only ca. 130 remaining individuals, due to intense human persecution approximately 100 years ago. Here, we analyse 41 ancient, historical, and modern mitochondrial genomes, to examine the number of female lineages involved in the postglacial recolonisation event and temporal changes in the Scandinavian brown bears' mitochondrial genetic diversity. Our results support the bi-directional recolonisation hypothesis, indicating multiple mitochondrial lineages from clade 1a possibly followed a southern route, while only a single lineage from clade 3a appears to have followed a northern route. Furthermore, we found that the recent bottleneck had a strong impact on the southern subpopulation, resulting in only one remaining haplotype in the contemporary brown bears. For the northern subpopulation, the impact was moderate, and most haplotypes were retained throughout the bottleneck. By exploring the postglacial recolonisation and recent population pressures, our study enhances understanding of how these factors have influenced the genetic diversity of Scandinavian brown bears.

Characterization of Nuclear Mitochondrial Insertions in Canine Genome Assemblies

BACKGROUND

The presence of mitochondrial sequences in the nuclear genome (Numts) confounds analyses of mitochondrial sequence variation, and is a potential source of false positives in disease studies. To improve the analysis of mitochondrial variation in canines, we completed a systematic assessment of Numt content across genome assemblies, canine populations and the carnivore lineage.

RESULTS

Centering our analysis on the UU_Cfam_GSD_1.0/canFam4/Mischka assembly, a commonly used reference in dog genetic variation studies, we found a total of 321 Numts located throughout the nuclear genome and encompassing the entire sequence of the mitochondria. A comparison with 14 canine genome assemblies identified 63 Numts with presence-absence dimorphism among dogs, wolves, and a coyote. Furthermore, a subset of Numts were maintained across carnivore evolutionary time (arctic fox, polar bear, cat), with eight sequences likely more than 10 million years old, and shared with the domestic cat. On a population level, using structural variant data from the Dog10K Consortium for 1879 dogs and wolves, we identified 11 Numts that are absent in at least one sample, as well as 53 Numts that are absent from the Mischka assembly.

CONCLUSIONS

We highlight scenarios where the presence of Numts is a potentially confounding factor and provide an annotation of these sequences in canine genome assemblies. This resource will aid the identification and interpretation of polymorphisms in both somatic and germline mitochondrial studies in canines.

Late Pleistocene polar bear genomes reveal the timing of allele fixation in key genes associated with Arctic adaptation

The polar bear (Ursus maritimus) occupies a relatively narrow ecological niche, with many traits adapted for cold temperatures, movement across snow, ice and open water, and for consuming highly lipid-dense prey species. The divergence of polar bears from brown bears (Ursus arctos) and their adaptation to their Arctic lifestyle is a well-known example of rapid evolution. Previous research investigating whole genomes uncovered twelve key genes that are highly differentiated between polar and brown bears, show signatures of selection in the polar bear lineage, and are associated with polar bear adaptation to the Arctic environment. Further research suggested fixed derived alleles in these genes arose from selection on both standing variation and de novo mutations in the evolution of polar bears. Here, we reevaluate these findings based on a larger and geographically more representative dataset of 119 polar bears and 135 brown bears, and assess the timing of derived allele fixation in polar bears by incorporating the genomes of two Late Pleistocene individuals (aged 130-100,000 years old and 100-70,000 years old). In contrast with previous results, we found no evidence of derived alleles fixed in present-day polar bears within the key genes arising from de novo mutation. Most derived alleles fixed in present-day polar bears were also fixed in the Late Pleistocene polar bears, suggesting selection occurred prior to 70,000 years ago. However, some derived alleles fixed in present-day polar bears were not fixed in the two Late Pleistocene polar bears, including at sites within APOB, LYST, and TTN. These three genes are associated with cardiovascular function, metabolism, and pigmentation, suggesting selection may have acted on different loci at different times.

Impact of Holocene environmental change on the evolutionary ecology of an Arctic top predator

The Arctic is among the most climatically sensitive environments on Earth, and the disappearance of multiyear sea ice in the Arctic Ocean is predicted within decades. As apex predators, polar bears are sentinel species for addressing the impact of environmental variability on Arctic marine ecosystems. By integrating genomics, isotopic analysis, morphometrics, and ecological modeling, we investigate how Holocene environmental changes affected polar bears around Greenland. We uncover reductions in effective population size coinciding with increases in annual mean sea surface temperature, reduction in sea ice cover, declines in suitable habitat, and shifts in suitable habitat northward. Furthermore, we show that west and east Greenlandic polar bears are morphologically, and ecologically distinct, putatively driven by regional biotic and genetic differences. Together, we provide insights into the vulnerability of polar bears to environmental change and how the Arctic marine ecosystem plays a vital role in shaping the evolutionary and ecological trajectories of its inhabitants.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

Defining biologically relevant and hierarchically nested population units to inform wildlife management

Wildlife populations are increasingly affected by natural and anthropogenic changes that negatively alter biotic and abiotic processes at multiple spatiotemporal scales and therefore require increased wildlife management and conservation efforts. However, wildlife management boundaries frequently lack biological context and mechanisms to assess demographic data across the multiple spatiotemporal scales influencing populations. To address these limitations, we developed a novel approach to define biologically relevant subpopulations of hierarchically nested population levels that could facilitate managing and conserving wildlife populations and habitats. Our approach relied on the Spatial "K"luster Analysis by Tree Edge Removal clustering algorithm, which we applied in an agglomerative manner (bottom-to-top). We modified the clustering algorithm using a workflow and population structure tiers from least-cost paths, which captured biological inferences of habitat conditions (functional connectivity), dispersal capabilities (potential connectivity), genetic information, and functional processes affecting movements. The approach uniquely included context of habitat resources (biotic and abiotic) summarized at multiple spatial scales surrounding locations with breeding site fidelity and constraint-based rules (number of sites grouped and population structure tiers). We applied our approach to greater sage-grouse (Centrocercus urophasianus), a species of conservation concern, across their range within the western United States. This case study produced 13 hierarchically nested population levels (akin to cluster levels, each representing a collection of subpopulations of an increasing number of breeding sites). These closely approximated population closure at finer ecological scales (smaller subpopulation extents with fewer breeding sites; cluster levels ≥2), where >92% of individual sage-grouse's time occurred within their home cluster. With available population monitoring data, our approaches can support the investigation of factors affecting population dynamics at multiple scales and assist managers with making informed, targeted, and cost-effective decisions within an adaptive management framework. Importantly, our approach provides the flexibility of including species-relevant context, thereby supporting other wildlife characterized by site fidelity.

Drivers of polar bear behavior and the possible effects of prey availability on foraging strategy

BACKGROUND

Change in behavior is one of the earliest responses to variation in habitat suitability. It is therefore important to understand the conditions that promote different behaviors, particularly in areas undergoing environmental change. Animal movement is tightly linked to behavior and remote tracking can be used to study ethology when direct observation is not possible.

METHODS

We used movement data from 14 polar bears (Ursus maritimus) in Hudson Bay, Canada, during the foraging season (January-June), when bears inhabit the sea ice. We developed an error-tolerant method to correct for sea ice drift in tracking data. Next, we used hidden Markov models with movement and orientation relative to wind to study three behaviors (stationary, area-restricted search, and olfactory search) and examine effects of 11 covariates on behavior.

RESULTS

Polar bears spent approximately 47% of their time in the stationary drift state, 29% in olfactory search, and 24% in area-restricted search. High energy behaviors occurred later in the day (around 20:00) compared to other populations. Second, olfactory search increased as the season progressed, which may reflect a shift in foraging strategy from still-hunting to active search linked to a shift in seal availability (i.e., increase in haul-outs from winter to the spring pupping and molting seasons). Last, we found spatial patterns of distribution linked to season, ice concentration, and bear age that may be tied to habitat quality and competitive exclusion.

CONCLUSIONS

Our observations were generally consistent with predictions of the marginal value theorem, and differences between our findings and other populations could be explained by regional or temporal variation in resource availability. Our novel movement analyses and finding can help identify periods, regions, and conditions of critical habitat.