A deep dive into fat: Investigating blubber lipidomic fingerprint of killer whales and humpback whales in northern Norway

Abstract In cetaceans, blubber is the primary and largest lipid body reservoir. Our current understanding about lipid stores and uses in cetaceans is still limited, and most studies only focused on a single narrow snapshot of the lipidome. We documented an extended lipidomic fingerprint in two cetacean species present in northern Norway during wintertime. We were able to detect 817 molecular lipid species in blubber of killer whales (Orcinus orca) and humpback whales (Megaptera novaeangliae). The profiles were largely dominated by triradylglycerols in both species and, to a lesser extent, by other constituents including glycerophosphocholines, phosphosphingolipids, glycerophosphoethanolamines, and diradylglycerols. Through a unique combination of traditional statistical approaches, together with a novel bioinformatic tool (LION/web), we showed contrasting fingerprint composition between species. The higher content of triradylglycerols in humpback whales is necessary to fuel their upcoming half a year fasting and energy‐demanding migration between feeding and breeding grounds. In adipocytes, we assume that the intense feeding rate of humpback whales prior to migration translates into an important accumulation of triacylglycerol content in lipid droplets. Upstream, the endoplasmic reticulum is operating at full capacity to supply acute lipid storage, consistent with the reported enrichment of glycerophosphocholines in humpback whales, major components of the endoplasmic reticulum. There was also an enrichment of membrane components, which translates into higher sphingolipid content in the lipidome of killer whales, potentially as a structural adaptation for their higher hydrodynamic performance. Finally, the presence of both lipid‐enriched and lipid‐depleted individuals within the killer whale population in Norway suggests dietary specialization, consistent with significant differences in δ15N and δ13C isotopic ratios in skin between the two groups, with higher values and a wider niche for the lipid‐enriched individuals. Results suggest the lipid‐depleted killer whales were herring specialists, while the lipid‐enriched individuals might feed on both herrings and seals.

Since 2011, several hundreds of killer whales (Orcinus orca) and humpback whales (Megapetara noveangliae) have gathered during the wintertime into some specific fjords of Troms in northern Norway to feed on the abundant schooling Norwegian springspawning (NSS) herring (Clupea harengus) that for unknown reasons have started to overwinter in these areas (Jourdain et al., 2020;Rikardsen, 2019). Humpback whales are long-range migratory species using the Norwegian Sea as a migration corridor. They travel between their summer feeding area characterized by high primary productivity mainly located in the northern Barents Sea east of the Svalbard archipelago (Leonard & Øien, 2020) and their breeding ground located in warm tropical waters of the West Indies or Cap Verde with low food availability during winter (Kettemer et al., in prep [https://whale track ing20 18.uit.no/]; Stevick et al., ,2003Stevick et al., , , 2006. As capital breeders, humpback whales accumulate large energy reserves in the feeding grounds to meet the costs of growth, maintenance, locomotion, and future reproduction (Jönsson, 1997;Lockyer & Brown, 1981;Oftedal, 2011;Stephens et al., 2009;Waugh et al., 2012). By contrast, as income breeders, killer whales meet their energy demand throughout the year by continual foraging and thus have less need to store large energy reserves (Jönsson, 1997;Stephens et al., 2009). Killer whales that feed off northern Norway largely rely upon the NSS herring, following their yearly migration along the Norwegian coast (Dietz et al., 2020;Jourdain et al., 2020;Similä et al., 1996;Vogel et al., 2021). Nonetheless, diverse individuals' diet preferences including both fish and marine mammal prey and alternative space-use strategy in the Barents Sea have also been indicated by stable isotopes and satellite tracking studies (Dietz et al., 2020;Foote et al., 2009;Jourdain et al., 2020).
Lipids function as a cellular energy reservoir and serve as major structural components of biological membranes such as plasma membranes and membranes of intracellular organelles (van Meer, 2005;van Meer & Vaz, 2005). They also participate in signaling as chemical messengers (Plakkal Ayyappan et al., 2016;Saliba et al., 2015). Glycerophospholipids and sphingolipids are major lipid constituents of cell membranes (Pomorski et al., 2001;van Meer et al., 2008), while triacylglycerols represent the bulk of lipid storage. The endoplasmic reticulum is a major biogenesis site of lipids, ultimately transported to other organelles as is the case for triacylglycerols in lipid droplets (Gregor & Hotamisligil, 2007;Lafontan & Langin, 2009;van Meer et al., 2008). Accordingly, lipid droplets can shrink and expand in a controlled manner in response to energetic demand and nutritional state (Kalantari et al., 2010). In short, when animals are feeding on a highly calorific diet, lipid storage increases and the lipid droplet expands in response to acute fat storage (Chitraju et al., 2012;Coelho et al., 2013).
To date, most studies investigating lipids in cetaceans only focused on a narrow snapshot of the lipidome, mainly on fatty acids and triacylglycerols (Bagge et al., 2012;Gross, 2020;Hansen & Cheah, 1969;Kawai et al., 1988;Krahn et al., 2004;Ruchonnet et al., 2006;Waugh et al., 2012Waugh et al., , 2014. With the advent of stateof-the-art analytical methods, together with the development of bioinformatics (e.g., Molenaar et al., 2019), the so-called "-omics" disciplines enable more comprehensive systemic profiling to understand the role of lipids in biological systems (Lydic & Goo, 2018;Miller et al., 2020). Thus, the objectives of this study were to document an extended lipidomic fingerprint in killer whales and humpback whales present in northern Norway during wintertime to better understand the determinants of lipid stores and uses in these cetaceans. (a) We described an extended lipid profile in blubber of killer whales and humpback whales. While triradylglycerols are expected to be the predominant lipid class in blubber for both species based on previous studies, we reasoned that the use of lipidomics in the present study would provide an in-depth profiling and therefore likely enable the identification and quantification of other minor lipid classes in cetacean blubber, which have so far been overlooked (e.g.,

K E Y W O R D S
capital breeder, cetacean, herring eater, income breeder, LION/web, lipid, seal eater, stable isotope Bernier-Graveline et al., 2021). (b) Through a unique combination of traditional statistical approaches, together with the use of a novel bioinformatic tool (Molenaar et al., 2019), we then compare the lipidomic fingerprint of both species. Although lipidomic composition similarities due to their phylogenetic affiliation were expected, the differences in response to their contrasted life-history strategies were also anticipated. Specifically, we hypothesized higher proportion of lipid storage in humpback whales in order to cope with their ongoing long fasting and energy-demanding migration at the time of sampling. (c) Finally, we investigated intraspecific variations of the lipidome within the killer whale population. Because of diverse individual dietary specialization within the killer whale population in Norway (Dietz et al., 2020;Jourdain et al., 2017Jourdain et al., , 2020, we sought to identify groups of individuals with closely related lipidomic fingerprint to search for specific feeding groups. To provide more information on potential dietary specialization in the killer whales, stable isotopes of nitrogen (δ 15 N) and carbon (δ 13 C) were determined in skin tissue of sampled individuals.

| Field sampling
Fieldwork was conducted in December 2017 (5 and 6) and January 2018 (13), in the coastal area of Skjervøy located in northern Norway (70°3′12.84″N, 20°57′23.07″E). We sampled 21 adult killer whales (18 males and 3 females) and 4 adult humpback whales (1 male and 3 females). All individuals were mostly biopsied on the main body in the region below or behind the dorsal fin with a biopsy dart with a 100-mm hollow stainless-steel tip (inner diameter: 6 mm, biopsy core maximum depth ~8.5 cm), shot from a crossbow. The dorsal region is enriched in lipids and considered as important fat storage, as compared to other body regions (Lockyer et al., 1984;Ruchonnet et al., 2006). The dart was attached to a string secured to the crossbow that enabled rapid recovery of the sample. Once collected, the biopsies were immediately removed from the dart to separate skin and blubber. Skin of humpback whales is generally thicker than skin of killer whales. The upper part of the skin (a few millimeters) was divided vertically into two parts, one part was used for stable isotope analyses (SIA) and the other was used for molecular sexing. The blubber closest to skin-blubber interface was used for lipid analysis (i.e., outer layer). Samples were snap-frozen and stored in a dry shipper during field trip duration. Back at the field station, samples for lipid analysis and SIA were transferred directly into a −80°C freezer.
The small piece of skin for molecular sexing was preserved in an ethanol-filled 2-ml cryotube at room temperature.

| Molecular sexing
Sex of the whales was identified from skin samples by molecular sexing tests using the methods described in Bérubé and Palsbøll (1996).
Each test consisted of one Y-specific fragment and one X-specific fragment in a multiplex PCR assay, as detailed in Appendix S1. Total genomic DNA was extracted in both killer whales and humpback whales using the DNeasy Blood & Tissue Kit (Qiagen), following the manufacturer's protocol ("Purification of total DNA from animal tissues" (spin-column protocol)), as detailed in the SI).

| Stable isotope analyses
Bulk nitrogen (δ 15 N) and carbon (δ 13 C) stable isotope ratios were analyzed from skin samples of killer whales (n = 19). The analyses were carried out following methods described in Marcoux et al. (2012) and Tartu et al. (2020). Briefly, skin samples were lyophilized at −48°C at a pressure of 133 × 103 mbar for 48 hr, homogenized into small pieces using a scalpel, lipid extracted using a 2:1 chloroform:methanol mixture to remove lipids (which can confound interpretation of δ 13 C data), and then the samples (400-600 μg) were placed into tin cups. Values of δ 13 C and δ 15 N were determined on a Delta V Advantage Thermo Scientific Continuous Flow Mass Spectrometer (Thermo Scientific, Bremen, Germany) coupled to a 4,010 Elemental Combustion System (Costech Instruments). Instrument accuracy met laboratory quality assurance criteria, as values were within 0.18‰ for NIST 8,547, NIST 8,573, and NIST 8,574 for δ 15 N, and 0.17‰ for NIST 8,573, 8,542, and 8,544 for δ 13 C. Precision was assessed with our laboratory standards (NIST 1577c, tilapia muscle, USGS 40, and Urea (n = 12 for each)), run every 15 samples, and was within 0.18. Bulk isotopes were used to investigate isotopic niche width as a proxy of the trophic niche (Newsome et al., 2007).
The δ 13 C of a predator reflects the origin of food sources, indicating the sources of primary production in the food consumed by the whales, for example, pelagic versus benthic, inshore versus offshore (Hobson, 1999;Hobson et al., 1996;Kelly, 2000). Typically, inshore/ benthic ecosystems are characterized by higher δ 13 C than pelagic/ offshore ecosystems. The δ 15 N is commonly used as an indicator of the trophic position of a consumer (Hobson et al., 1994(Hobson et al., , 1996Kelly, 2000;Newsome et al., 2007), owing to the stepwise enrichment from food source to consumer. Controlled diet experiment estimated the half-life rates in skin of bottlenose dolphin (Tursiops trucuntus) to be 48 ± 19 days for nitrogen and 24 ± 8 days for carbon (Giménez et al., 2016). Therefore, isotopic values measured in killer whales are expected to represent the individuals' feeding habits in the 4-7 weeks prior sampling.

| Statistics and data analysis
We first characterized the lipidomic fingerprint of blubber from killer whales and humpback whales. Molecular lipid species were grouped, and their concentrations summed according to which lipid class they belong to, in order to facilitate profile visualization. We then compared the composition of lipidomic fingerprint between killer whales and humpback whales (in terms of proportion, each lipid species represented in each sampled individual). Finally, we investigated intraspecific variations of lipidomic signature within the killer whale population. Proportions of each lipid species were normalized as the centered log-ratio concentrations (ln(lipid x /∑lipid)) to minimize any potential misleading results (e.g., due to body size difference between species). Testing the influence of sex was not possible due to the imbalanced sex ratio in each species (humpback whales: 1 male and 3 females, killer whales: 18 males and 3 females). Except for analyses performed with LION/web (Molenaar et al., 2019), all statistics have been carried out using R version 4.0.3 (R core Team 2020). A significance level of p-value = .05 was used for all tests.

| Interspecific comparison: killer whales versus humpback whales
We first performed a partial least-squares discriminant analysis (PLS-DA) on the lipidome using the "mix.omics" R package (Rohart et al., 2017) and tested for potential species clustering using the MANOVA Pillai's trace test from the "car" R package (Fox & Weisberg, 2019). The Pillai test is consistent with the use of an asymmetric data set. We further performed lipid ontology (LION) enrichment analysis to compare the lipidomic fingerprint of both species.
To do so, we used a novel bioinformatic tool developed for global and in-depth lipidomic data mining of mammalian system, the LION/ web-based interface (Molenaar et al., 2019). LION/web is an ontology database containing information related to lipid metabolism, which associates >50,000 lipid molecular species to biophysical, chemical, and cell biological features (referred as "LION terms"). In the present study, we used LION/web to connect the molecular lipid species with lipid classification, function, and distribution within the cellular components (i.e., the predominant cellular localization). As a general guideline, the lipid function branch comprises 3 subcategories: "membrane component" (i.e., associated with lipids that are primarily regarded as a structural component of lipid membranes), "lipid mediated signaling" (i.e., lipids that have been implicated in signaling), and "lipid storage" (i.e., lipids that are associated with storage). Enrichment analysis was performed using the ranking mode as described in Molenaar et al. (2019), on normalized lipid log-ratio concentrations to compare the lipidomic fingerprint of killer whales and humpback whales and report any potential lipid LION term enrichments in the condition of interest. The analysis has been parameterized with killer whales as the control condition and humpback whales as the condition of interest. Calculated p-values from LION/ web were adjusted for multiple testing (i.e., Benjamini-Hochberg;Benjamini & Hochberg, 1995

| Intraspecific variations in killer whales
We first performed a heat map to visualize the lipidomic fingerprint of each killer whale using the LION-PCA heat map module of LION/ web (Molenaar et al., 2019), which is a modified version from a gene ontology-PCA (Wagner, 2015) specifically developed for global and in-depth lipidomic data mining. Lipidomic signatures were sorted following a hierarchical clustering per individual in order to identify any potential groups of individuals with closely related lipidomic fingerprint. Differences between the two identified clusters (i.e., "lipid-enriched" and "lipid-depleted") were further confirmed with a PLS-DA and the MANOVA Pillai's trace test. We also compared the nitrogen (δ 15 N) and carbon (δ 13 C) isotopic ratios of the "lipidenriched" and "lipid-depleted" killer whales using linear models. The isotopic niche width (inferred from δ 13 C and δ 15 N) of the two identified groups was illustrated by standard ellipses (containing ∼95% of the data) on an isotopic biplot using the "SIBER" R package (Jackson et al., 2011). The areas of the resultant ellipses were then computed using both the maximum-likelihood (SEAc, adjusted for small sample size) and the Bayesian approaches (SEAb; parametrized as detailed in Jackson et al., 2011). SEA is a measure of the standard deviation for bivariate data. Estimated SEA values were directly compared in a probabilistic manner in terms of similarity between the two groups of killer whales.

| Interspecific comparison: killer whales versus humpback whales
PLS-DA indicated that blubber lipidomic fingerprint significantly differ between both species (Figure 2; MANOVA, Pillai's trace test, F = 20.98, p-value < .001). Out of the 817 molecular lipid species, 791 were matched to the LION/web database (i.e., 97%), meaning that 26 lipid species were not considered in the LION/web analysis.

| Intraspecific variations in killer whales
The lipidomic fingerprints of killer whales appeared to cluster into two distinct groups (Figure 4). The first group is composed of 13 individuals with enriched lipid content (referred as "lipid-enriched") compared with the second group composed of 8 individuals with depleted lipid content (referred as "lipid-depleted"). The suggested difference in lipidomic fingerprint between the two groups was further confirmed by the PLS-DA, which indicates significant differences ( Figure S2; MANOVA, Pillai's trace test, F = 7.90, p-value < .001).

| D ISCUSS I ON
To the best of our knowledge, this study provides the most comprehensive assessment of the lipidome in cetaceans. We were able to detect 817 molecular lipid species in blubber of killer whales and humpback whales. Through a unique combination of traditional statistical approaches, together with enrichment analyses using a novel bioinformatic tool specifically developed for global and in-depth F I G U R E 1 Proportion of lipid class (mean ± SD) in blubber lipidomic fingerprint of adult killer whales (n = 21) and humpback whales (n = 4) off northern Norway lipidomic data mining (Molenaar et al., 2019), we showed contrasting fingerprints between both species. Specifically, we found higher content of lipids involved in storage in humpback whales, which translates into higher triacylglycerol content stored in adipocyte lipid droplets. In contrast, killer whales showed higher content of membrane components, which translates into higher sphingolipid content. The reported differences in the lipidomic composition between both species are consistent with their contrasted life-history strategy (i.e., capital versus income breeders). Finally, the presence of both lipid-enriched and lipid-depleted individuals within the killer whale population in Norway suggests dietary specialization, consistent with significant differences in nitrogen (δ 15 N) and carbon (δ 13 C) stable isotope ratios in skin tissue between the two groups.  Kawai et al., 1988), killer whales (Krahn et al., 2004), fin whale (Balaenoptera physalus ;Lockyer et al., 1984;Ruchonnet et al., 2006), sei whale (Balaenoptera borealis; Bottino, 1978), and humpback whales (Waugh et al., 2012).

| Lipidomic fingerprint in blubber of killer whales and humpback whales
Besides, these former studies largely focused on triacylglycerols, and so far, an exceedingly limited number of studies have focused on minor constituents, despite their essential functions (Bernier-Graveline et al., 2021). With the advent of state-of-the-art analytical methods, lipidomics enable a better systemic profiling. Accordingly, our results underscored the substantial contribution of glycerophosphocholines (~15%) and phosphosphingolipids (~5%) in the blubber lipidomic fingerprint of humpback whale, while pioneering studies conducted on other baleen whale species, the gray whale (Krahn et al., 2001), the sei whale (Bottino, 1978), and the fin whale (Lockyer et al., 1984), reported significantly lower phospholipid content (<4%, 3%, and ~5%, respectively). Similarly, our study highlights the substantial contribution of phosphosphingolipids (~5%), glycerophosphocholines (~5%), and glycerophosphoethanolamines (~3%) in the blubber lipidomic fingerprint of killer whale, while a former study conducted on the same species and on white whales failed to detect the presence of any phospholipids (Krahn et al., 2004). Therefore, the present study provides new knowledge and baseline information for future investigations F I G U R E 2 PLS-DA of blubber lipidomic fingerprint in adult killer whales (n = 21; orange triangles) and humpback whales (n = 4; blue circles) off northern Norway. Ellipses represent 95% confidence interval. The colored background surface represents the prediction areas PC1: 66% expl. var PC2: 12% expl. var on lipids and energy stores in marine mammals. Lipid concentrations reported in Table S1 should be interpreted with caution as lipid loss may have occurred at sampling due to the force of the biopsy cutting tip (Ryan et al., 2013). As a result, the lipid concentrations of blubber measured from biopsy darting might have been underestimated as compared to the tissue in situ. Such lipid loss could occur without affecting the lipid class profiles (Krahn et al., 2004;Waugh et al., 2014).

| Interspecific comparison: Killer whales versus humpback whales
Blubber lipidomic fingerprint differs between killer whales and humpback whales. Yet, one might expect a higher degree of similarities between two phylogenetically affiliated species. Recently, a large-scale analysis showed that changes in the mammalian lipidomes do not accumulate proportionally to the phylogenetic distances, but rather represent phenotypic differences including species-specific adaptations (Khrameeva et al., 2018). In that respect, the reported difference in the lipidomic fingerprint composition between killer whales and humpback whales is consistent with their contrasted life-history strategy.
Among the cetaceans, baleen whales are typically capital breeders, while toothed whales are generally income breeders (Irvine et al., 2017;Oftedal, 1997). Consistent with our results, Irvine et al., (2017) showed differences in cetacean energy stores between capital and income breeders in Australian waters, reporting that humpback whales accumulate more lipids than sperm whales (Physeter macrocephalus) in order to fuel their annual migration.
According to Bengtson Nash et al., (2013), the utilization of humpback whale blubber lipids during the migration journey across the Southern Ocean was evidenced by an average of 23% reduction in lipids in outer blubber layer between early and late migrating cohorts. As capital breeders, humpback whales require to accumulate and store large energy reserves in their feeding grounds, in order to meet the costs of growth, maintenance, locomotion, and reproduction in the breeding grounds (Jönsson, 1997;Oftedal, 2011;Stephens et al., 2009;Waugh et al., 2012). Specifically, the humpback whales investigated in the present study belong to a long-range migratory population using the Norwegian Sea as a migration corridor between their main summer feeding areas characterized with high primary productivity located in the northern Barents Sea east of the Svalbard archipelago and their breeding ground in warm tropical waters of the West Indies or Cap Verde with low food availability during late winter and spring (Kettemer et al., in prep [https://whale F I G U R E 3 LION enrichment analysis of blubber lipidomic fingerprint composition in adult killer whales (n = 21) versus humpback whales (n = 4) off northern Norway. The gray vertical line indicates the cutoff value of significant enrichments. The size and color (from gray to red) of the horizontal bars (x-axis) are scaled with the enrichment. Killer whales were considered as the control condition and humpback whales as the condition of interest. "up" refers to an enrichment, and "down" refers to a depletion track ing20 18.uit.no/; Stevick et al., 2003;. The presence of humpback whales in northern Norway during wintertime represents a stop-over area to refuel at an early stage of their migration route, and thus a unique opportunity to replenish their energy stores before undertaking a 6-month fasting and energy-demanding journey across the north Atlantic Ocean diagonal (Rikardsen, 2019). In adipocytes, we assume that the intense feeding rate of humpback whales translates into an important accumulation of triacylglycerol content in lipid droplets. Upstream, the endoplasmic reticulum is certainly operating at full capacity to supply the acute storage of lipids (Gregor & Hotamisligil, 2007;van Meer et al., 2008), which is consistent with the reported enrichment of glycerophosphocholines in humpback whales, major components of the endoplasmic reticulum (van Meer & de Kroon, 2011). To gain a complete understanding of lipid storage and shed light on mobilization of lipid species (Waugh et al., 2012), further studies investigating lipidomic fingerprint of humpback whales when they return from the migration are required. As income breeders, killer whales meet their energy demand throughout the year by continuous foraging, and thus, storing large reserves of energy is less essential (Jönsson, 1997;Stephens et al., 2009). Specifically, the killer whales investigated in the present study largely rely upon the NSS herrings (Jourdain et al., 2020;Mul et al., 2020;Similä et al., 1996;Vogel et al., 2021), following their yearly migration along the Norwegian coast (Dietz et al., 2020;Similä et al., 1996), although alternative space-use strategy into the Barents Sea region toward Novaya Zemlya Island has been recently pointed out in three satellite-tagged individuals (Dietz et al., 2020).
As compared to the long-range migration of humpback whales, killer whales travel across shorter distances and spend most of their time performing area-restricted search behavior (Dietz et al., 2020;Mul et al., 2020;Vogel et al., 2021), indicating continual foraging and thus less lipid storage requirements.
In addition, the blubber has an important structural role and function as a streamlining tissue in cetaceans (Fish et al., 2008;Iverson & Koopman, 2018). Killer whales are efficient predators with specialized and elaborated feeding strategies, which require strong agility, high-speed swimming, precise body movements, and maneuverability. In northern Norway, they are known to perform carousel feeding, which consists of herding the herrings into a tight and dense fish ball close to the surface, tail slapping the edge of the
During the winter feeding on herring in the north Norwegian fjords, many of the killer whales are also documented to feed around the herring fishing vessels (Mul et al., 2020;Rikardsen, 2019). Moreover, previous field observations also reported killer whales chasing, harassing, and consuming harbor (Phoca vitulina) and gray (Halichoerus grypus) seals in Norwegian coastal waters (Jourdain et al., 2017;Vongraven & Bisther, 2014). In comparison, humpback whales are low-speed swimmers taking benefit from their large body size and using lunge feeding to swallow large amounts of prey (Goldbogen et al., 2008;Simon et al., 2012). In northern Norway, humpback whales are known to perform surface lunging through herring balls herded by killer whales (Jourdain et al., 2017;Rikardsen, 2019 Although the present study provides relevant insights about lipids stores and uses in cetaceans, the low sample size of humpback whales precludes to make any firm conclusions and further studies including a larger dataset are needed. The use of remote biopsy darting to collect blubber and skin has marked a step forward to F I G U R E 5 Comparison between δ 15 N and δ 13 C isotopic signatures between the "lipid-enriched" (n = 12, red) and "lipiddepleted" (n = 7, yellow) clusters in killer whales off northern Norway. Significant differences between the two groups are highlighted with an "*" δ N 15 δ C 13 * Lipid-enriched Lipid-depleted * F I G U R E 6 Isotopic niche width (inferred from δ 13 C and δ 15 N in skin) illustrated by standard ellipses (containing ∼95% of the data and computed with "SIBER" R package) of the "lipid-enriched" (red triangles) and "lipid-depleted" (yellow circles) killer whales off northern Norway (n = 19). Comparison of the standard ellipse area (SEA) between the two groups Lipid-enriched L ipid-depleted SEA(‰ 2 ) δ N 15 δ 13 C investigate physiological questions in free-living cetaceans (Hunt et al., 2013). However, inherent challenges of this sampling technique do not allow systematic collection of the full blubber core.
Thus, the findings presented here from the outer layer would gain benefits from complementary investigations conducted on healthy stranded animals, focusing on the full blubber core or the inner layer, which appears to be the most metabolically active blubber part (e.g., Lockyer et al., 1984; but see Waugh et al., 2014 andAguilar &Borrell, 1990).

| Intraspecific variations in killer whales
The overwintering in these areas (Dietz et al., 2020;Jourdain et al., 2017;Rikardsen, 2019;Similä et al., 1996). Moreover, the presence of gray seals pupping grounds in northern Norway between September and December might represent a valuable pinniped food source at this time of the year (Haug et al., 1994;Nilssen 2007). In that context, the lipid-depleted killer whales might fall into herring specialists (38% of all individuals), while the lipid-enriched individuals might feed on both herrings and seals (62% of all individuals).
Importantly, our conclusion about the diets of the two killer whale groups was supported by significantly higher values of δ 15 N and δ 13 C in lipid-enriched killer whales as compared to the lipid-depleted individuals, indicative of prey from higher trophic position and living in near-shore coastal habitats. Moreover, the wider isotopic niche of the lipid-enriched group is consistent with a mixed diet composed of both seals and herring, while the narrow and overlapping isotopic niche of the lipid-depleted killer whales is indicative of prey specialization exclusively composed of herring (Foote et al. 2012;Jourdain et al., 2020;Samarra et al., 2017;Vogel et al., 2021). As a matter of evidence, average isotopic values measured in our study for each group (lipid-enriched: δ 15 N = 12.7 ± 0.8‰ and δ 13 C = −18.9 ± 1.0‰; lipid-depleted: δ 15 N = 11.9 ± 0.4‰ and δ 13 C = −19.7 ± 0.3‰) were perfectly in line with those reported in Jourdain et al., (2020) (seal eaters: δ 15 N = 12.6 ± 0.3‰ and δ 13 C = −18.9 ± 0.6‰; herring eaters: δ 15 N = 11.7 ± 0.2‰ and δ 13 C = −19.3 ± 0.3‰). With a relatively high nutritional value (Beck et al., 2003), the occurrence of seals in the diet of killer whales may be considered, at first glance, as being beneficial. However, this also raises the issue of exposure to environmental pollutants as "preying on seals pushes killer whales from Norway above pollution effects thresholds" (Andvik et al., 2020).

This work was financed by the Fram Center Flagship Hazardous
Substances (WhaleHealth project; PI: Jenny Bytingsvik) and the WhaleFeast project supported by the Regional Research Council No. 282469. This study was approved by the Norwegian Animal Research Authority (FOTS ID 13639 and 14135). We thank Trond Johnsen and Evert Mul for their help collecting the samples in the field and Anders Goksøyr to communicate about LION/web. The authors thank the two reviewers and editor for their comments, which have greatly improved the manuscript.

CO N FLI C T O F I NTE R E S T
None declared.