Each species occupies a distinct ecological niche, defined by a specific set of environmental conditions and resource requirements necessary for its survival and reproduction. However, climate change is altering species distributions, predator-prey relationships and resource partitioning between species with these changes being pronounced in the Arctic. Stable isotope analysis of carbon (δ¹³C) and nitrogen (δ¹⁵N) has been widely used to estimate isotopic niches and quantify niche overlap among species, a two-dimensional approach (2D). However, δ¹³C is not always sufficient to differentiate habitat and resource use among species due to minimal variation between end-members. Incorporating sulfur stable isotopes (δ³⁴S) can enhance resolution in such cases. Using an Arctic coastal food web as a model system, we used a three-dimensional isotopic niche approach (3D: δ¹³C-δ¹⁵N-δ³⁴S) with 717 individuals across 69 species spanning multiple taxonomic groups (invertebrates, fish, seabirds, and marine mammals) that utilize resources from benthic and pelagic habitats. We compared the traditional 2D isotopic niches with a 3D framework using nicheROVER to assess how the addition of a third dimension changes niche size estimates and probability of niche overlap between species. We found that benthic-associated species, such as common eider (Somateria mollissima) and various benthic invertebrates, exhibited greater changes in isotopic niche size with the addition of δ³⁴S than pelagic-associated species. In addition, niche overlap among benthic-associated taxa decreased with the 3D approach, indicating better resolution of habitat use and resource partitioning. This finding likely reflects the greater ecological diversity, foraging specialization and more complex food web structure characteristic of benthic ecosystems. We recommend incorporating δ³⁴S for aquatic studies that involve benthic habitats and emphasize the value of multidimensional approaches in ecological niche analysis.

Spatiotemporal environmental heterogeneity is a major evolutionary driver, which can cause profound phylogeographic complexity, particularly at the periphery of species ranges. Ringed seals display a highly disjoint distribution, occurring in high abundance throughout the circumpolar Arctic, as well as in the Baltic Sea, Lake Saimaa and Lake Ladoga. These relict Fennoscandian ringed seals were traditionally regarded as originating from a single colonisation event after the Last Glacial Maximum (LGM), however recent studies have challenged this perception. Here, we analyse 246 mitogenomes and 180 skulls to unravel the diversity and spatiotemporal pattern of diversification in Fennoscandian ringed seals. Contrary to previous assumptions, our results reveal a complex evolutionary history characterised by pre-LGM diversification from Arctic ringed seals and possibly several Fennoscandian colonisation events. We hypothesize that Saimaa seals originate from Arctic ringed seals, from which they diverged prior to their arrival in Lake Saimaa. Ladoga seals appear to also originate from the Arctic, with secondary colonisation events from paleo-Skagerrak-Kattegat-Baltic, while the Baltic ringed seal have mixed evolutionary origins. Lake Saimaa and to some extent Lake Ladoga ringed seals have experienced a loss of diversity and evolved divergent skull morphologies, likely as a result of colonisation bottlenecks, isolation and dietary specialisation, while Baltic Sea ringed seals have retained remarkably high levels of genetic and morphological diversity. Our study supports the classification of Saimaa, Ladoga and Baltic ringed seals as distinct taxa, and highlights the need for management and conservation efforts to mitigate cumulative impacts of human activities and climate change on Fennoscandian ringed seals. 

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 many 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.