The Arctic environment plays a critical role in the global climate system and marine biodiversity. The region’s ice-covered expanses provide essential breeding and feeding grounds for a diverse assemblage of marine species, who have adapted to thrive in these harsh conditions and consequently are under threat from global warming. The bearded seal (Erignathus barbatus) is an ice-obligate Arctic species using sea-ice for all aspects of its life-history, rendering it particularly vulnerable to sea-ice loss. It is one of the least studied and hence enigmatic of the Arctic marine mammals, with little knowledge regarding genetic structure, diversity, adaptations and demographic history, consequently hampering management and conservation efforts. Here, we sequenced 70 whole nuclear genomes from across most of the species’ circumpolar range, finding significant genetic structure between the Pacific and the Atlantic subspecies, which diverged during the Penultimate Glacial Period (~192 KYA). Remarkably, we found fine-scale genetic structure within both subspecies, with at least two distinct populations in the Pacific and three in the Atlantic. We hypothesize sea-ice dynamics and bathymetry had a prominent role of in shaping bearded seal genetic structure and diversity. Resulting genomic data can be used to complement the health, physiological, and behavioral research needed to conserve this species. In addition, we provide recommendations for management units that can be used to more specifically assess climatic and anthropogenic impacts on bearded seal populations.

Identification of phenotypic characteristics in reproductively successful individuals provides important insights into the evolutionary processes that cause range shifts due to environmental change. Female beluga whales (Delphinapterus leucas) from the Baffin Bay region (BB) of the Canadian Arctic in the core area of the species’ geographic range have larger body size than their conspecifics at the southern range periphery in Hudson Bay (HB). We investigated the mechanism for this north and south divergence as it relates to ovarian reproductive activity (ORA = total corpora) that combines morphometric data with ovarian corpora counted from female reproductive tracts. Based on the previous finding of reproductive senescence in older HB females, but not for BB whales, we compared ORA patterns of the two populations with age and body length. Female beluga whale ORA increased more quickly with age (63% partial variation explained) in BB than in HB (41%). In contrast, body length in HB female beluga whales accounted for considerably more of the total variation (12 vs 1%) in ORA compared to BB whales. We speculate that female HB beluga whale ORA was more strongly linked with body length due to higher population density resulting in food competition that favors the energetic advantages of larger body size during seasonal food limitations. Understanding the evolutionary mechanism of how ORA varies across a species’ range will assist conservation efforts in anticipating and mitigating future challenges associated with a warming planet.

Identifying phenotypic characteristics of evolutionarily fit individuals provides important insight into the evolutionary processes that cause range shifts with climate warming. Female beluga whales (Delphinapterus leucas) from the Canadian high Arctic (BB) residing in the core region of the species’ geographic range are 14% larger than their conspecifics at the southern periphery in Hudson Bay (HB). We investigated the causal mechanism for this north (core)-south (periphery) difference as it relates to fitness by combining morphometric data with ovarian corpora counted in female reproductive tracts. We found evidence for reproductive senescence in older HB females from the southern peripheral population but not for BB whales. Female beluga whale fitness in the more-northern BB increased faster with age (48% partial variation explained) versus a more gradual slope (25%) in HB. In contrast, body length in HB female beluga accounted for five times more of the total variation in fitness compared to BB whales. We speculate that female HB beluga fitness was more strongly linked with body length due to higher density, as larger body size provides survival advantages during seasonal food limitations. Understanding the evolutionary mechanism of how fitness changes will assist conservation efforts in anticipating and mitigating future challenges to peripheral populations.