# Seasonal lipid dynamics of four Arctic bivalves: Implications for their physiological capacities to cope with future changes in coastal ecosystems

Guillaume Bridier, Frédéric Olivier, Jacques Grall, Laurent Chauvaud, Mikael Sejr, Réjean Tremblay

<figure><img src="image_1.png" /></figure>

## To cite this version:

Guillaume Bridier, Frédéric Olivier, Jacques Grall, Laurent Chauvaud, Mikael Sejr, et al.. Seasonal lipid dynamics of four Arctic bivalves: Implications for their physiological capacities to cope with future changes in coastal ecosystems. Ecology and Evolution, 2023, 13 (11), 10.1002/ece3.10691. hal-04302705

HAL Id: hal-04302705

https://hal.univ-brest.fr/hal-04302705

Submitted on 24 Nov 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Seasonal lipid dynamics of four Arctic bivalves: Implications for their physiological capacities to cope with future changes in coastal ecosystems

Guillaume Bridier<fnref n="1" /> ID |

Frédéric Olivier<fnref n="2" /> ID |

Jacques Grall<fnref n="3" />,<fnref n="4" /> |

Laurent Chauvaud<fnref n="3" /> |

Mikael K. Sejr<fnref n="5" /> | Réjean Tremblay<fnref n="1" />

## Correspondence

Guillaume Bridier, Institut des Sciences de la Mer, Université du Québec à Rimouski, 310 Allée des Ursulines, Rimouski, QC G5L 2Z9, Canada.

Email: guillaume.bridier@live.ie

## Funding information

Université de Bretagne Occidentale/ LIA BeBEST; Brittany Regional Council (Allocation de Recherche Doctorale); Université du Québec à Rimouski; Institut Polaire Français Paul-Emile Victor (IPEV - PRIVARC project); European H2020 INTAROS program; Agence Nationale de la Recherche (ANR GAAP); Fondation UBO (FAC 2018); Observatoire Marin de l'IUEM (UMS3113)

## Abstract

The Arctic is exposed to unprecedented warming, at least three times higher than the global average, which induces significant melting of the cryosphere. Freshwater inputs from melting glaciers will subsequently affect coastal primary production and organic matter quality. However, due to a lack of basic knowledge on the physiology of Arctic organisms, it remains difficult to understand how these future trophic changes will threaten the long-term survival of benthic species in coastal habitats. This study aimed to gain new insights into the seasonal lipid dynamics of four dominant benthic bivalves (Astarte moerchi, Hiatella arctica, Musculus discors, and Mya truncata) collected before and after sea ice break-up in a high-Arctic fjord (Young Sound, NE Greenland). Total lipid content and fatty acid composition of digestive gland neutral lipids were analyzed to assess bivalve energy reserves while the fatty acid composition of gill polar lipids was determined as a biochemical indicator of interspecies variations in metabolic activity and temperature acclimation. Results showed a decrease in lipid reserves between May and August, suggesting that bivalves have only limited access to fresh organic matter until sea ice break-up. The lack of seasonal variation in the fatty acid composition of neutral lipids, especially essential ω3 fatty acids, indicates that no fatty acid transfer from the digestive glands to the gonads occurs between May and August, and therefore, no reproductive investment takes place during this period. Large interspecies differences in gill fatty acid composition were observed, which appear to be related to differences in species life span and metabolic strategies. Such differences in gill fatty acid composition of polar lipids, which generally influence metabolic rates and energy needs, may imply that not all benthic species will be equally sensitive to future changes in primary production and organic matter quality in Arctic coastal habitats.

## KEYWORDS

Arctic coastal ecosystems, bivalve, climate change, energy reserves, fatty acids, homeoviscous adaptation, meltwater inputs, metabolic rate depression

<footnote n="1">¹Institut des Sciences de la mer de Rimouski, Université du Québec à Rimouski, Rimouski, Quebec, Canada</footnote>

<footnote n="2">²Biologie des Organismes et Ecosystèmes Aquatiques (BOREA) UMR 8067 MNHN, CNRS, SU, IRD 207, UCN, UA, Paris, France</footnote>

<footnote n="3">³Laboratoire des Sciences de l'Environnement Marin (LEMAR) UMR 6539 UBO, CNRS, IRD, Ifremer, Plouzané, France</footnote>

<footnote n="4">⁴Observatoire Marin de l'Institut Universitaire Européen de la Mer UMS 3113, Université de Bretagne Occidentale, Plouzané, France</footnote>

<footnote n="5">⁵Arctic Research Centre and Ecoscience, Aarhus University, Aarhus C, Denmark</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
© 2023 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.</footnote>
# TAXONOMY CLASSIFICATION

## Ecophysiology

## 1 | INTRODUCTION

Arctic ecosystems are undergoing unprecedented changes due to atmospheric warming at a rate of three to four times higher than the global average (Rantanen et al., 2022). Since the 1970s, rising temperatures have led to a 40% reduction in summer (September) sea ice extent and a 65% reduction in sea ice thickness (AMAP, 2017; Kwok, 2018). These environmental changes have caused a shift in phytoplankton physiology (Ardyna et al., 2014) and a 57% increase in phytoplankton productivity between 1998 and 2018 (Lewis et al., 2020). However, these trends mask strong disparities at smaller spatial scales during the same period, as primary production has increased very strongly on inflow shelves (e.g., Fram Strait, Bering Sea), while remaining stable or decreasing on outflow shelves (e.g., Canadian Arctic Archipelago, Ardyna & Arrigo, 2020; Tremblay et al., 2015). The effects of climate change are also heterogeneous in coastal areas, which are exposed to both specific local (e.g., freshwater inputs) and broader constraints acting at the entire Arctic scale (Bridier, 2020; Wassmann et al., 2020). Inputs of meltwater and terrestrial particles in coastal areas modify the composition of phytoplankton and bacterial communities through changes in seawater stratification and turbidity (Arendt et al., 2010; Delpech et al., 2021; Paulsen et al., 2017), which reduces phytoplankton production (Holding et al., 2019; Meire et al., 2017; Murray et al., 2015). In addition, meltwater inputs provide refractory terrestrial organic carbon from melting glaciers and thawing permafrost to coastal marine ecosystems (Hernes et al., 2021; Wadham et al., 2019). All these environmental changes subsequently affect the coastal food webs by changing the origin of food sources and reducing their quality, in terms of polyunsaturated fatty acids (PUFAs), available to benthic primary consumers (Bridier et al., 2019; Kędra et al., 2015). Such a decrease in PUFA from organic matter could have subsequent consequences on the physiology of benthic invertebrates, especially since some PUFAs are considered essential (essential fatty acids (EFAs), e.g., 20:5ω3, 22:6ω3). EFAs are recognized to be weakly synthesized de novo by marine invertebrates but are essential for their growth and survival through their crucial role in physiological function (Parrish, 2009).

The fatty acids stored in reserve tissues (the neutral lipid fraction) serve as the main energy reserve in benthic invertebrates and are mobilized to meet metabolic needs (Haider et al., 2020; Mayrand et al., 2017; Pazos et al., 1997), or to support growth and reproduction (e.g., egg production, Leroy et al., 2013). In contrast to fatty acids used as energy stores, those incorporated into cell membranes (the polar lipid fraction) play a crucial role in the functioning of the organism by maintaining the integrity of cell membranes, particularly in relation to changes in temperature. In poikilotherm animals, the fluidity of cell membranes decreases at lower temperatures (Hochachka & Somero, 2002). Reduced membrane fluidity leads to a decrease in the activity of membrane-bound proteins and ion

diffusion rates, ultimately resulting in a decrease in metabolism (Hochachka & Somero, 2002). To counteract these effects, poikilotherms can adjust their membrane lipid composition when they are exposed to cold temperature (a process called Homeoviscous adaptation) by increasing the proportion of monounsaturated fatty acids (MUFAs) and PUFAs relative to saturated fatty acids, lengthening the carbon chains of fatty acids, or decreasing the cholesterol content of cell membranes (Hazel, 1995; Pernet et al., 2006; Sinensky, 1974). According to the membrane-pacemaker theory of aging (Hulbert, 2005; Hulbert et al., 2007), the fluidity of cell membranes is also influenced by their fatty acid composition, which affects the metabolic level and longevity of the organism. A high metabolic level promotes the formation of reactive oxygen species that are prone to damage cell membranes, especially phospholipids composed of PUFA (Kraffe et al., 2004). To avoid the oxidation of lipid membranes, long-lived species reduce the level of unsaturation in their phospholipids (i.e., the number of double bonds of their fatty acids) with subsequent negative consequences on the metabolic rate (Abele et al., 2009; Blier et al., 2017; Hulbert et al., 2007).

In this context, the presence of sufficient amounts of EFA in organic matter sources ingested by primary consumers is crucial to maintain the integrity of their cell membranes and to ensure their optimal physiological function. Thus, the availability of EFA in organic matter might be expected to ultimately determine the survival of some species in their environment. However, because of the relatively limited understanding of the physiology of Arctic species and their requirements for EFA, it remains a challenge to determine their vulnerability to future trophic changes in coastal ecosystems. This lack of fundamental knowledge is particularly evident for benthic invertebrates (Renaud et al., 2015, 2019), which are less studied than other taxa such as fish or seabirds (Poloczanska et al., 2016), despite their key role in Arctic marine food webs (e.g., Grebmeier, 2012; Kędra et al., 2015). Therefore, the objectives of this study are to: (1) provide new basic information on the physiology of bivalves with different longevity (e.g., cold acclimation and metabolism), (2) identify some basis of physiological requirements of benthic invertebrates in terms of EFAs, and (3) propose new hypotheses on the physiological capacity of these species to cope with future changes in trophic environment quality in Arctic ecosystems.

## 2 | MATERIALS AND METHODS

### 2.1 | Study site and choice of model species

Samples were collected in Young Sound, a High-Arctic fjord that is part of the Young Sound-Tyrolerfjord system located in the northeast Greenland (74.2-74.3°N, 19.7-21.9°W, Figure 1). Due
to the long duration of sea ice cover (9–10 months per year), this fjord is considered one of the least productive coastal ecosystems in the world, with pelagic primary production estimated at only $10 \text{ g C m}^{-2} \text{ year}^{-1}$ (Rysgaard et al., 1999). The Young Sound-Tyrolerfjord system is 90 km long, varies in width between 2 and 7 km, and covers an area of $390 \text{ km}^2$ (Rysgaard et al., 2003). The maximum depth is 330 m while the average depth is around 100 m (Rysgaard et al., 2003). The Young Sound-Tyrolerfjord system is characterized by two shallow sills, the outermost of which is 45 m deep and separates the water masses of the fjord from those of the Greenland Sea (Bendtsen et al., 2007). The fjord is exposed to meltwater inputs from June to September, to total volume of which is estimated to be about $0.9–1.4 \text{ km}^3$ per year (Bendtsen et al., 2014). The land-terminating glaciers contribute significantly to this discharge, accounting on average for 50%–80% of the annual freshwater inputs, especially in the inner region of the fjord (Citterio et al., 2017). These summer meltwater plumes induce strong spatial gradients in salinity and temperature in surface waters above the pycnocline (i.e., 5–10 m), ranging from 8 to 30 PSU and from to 7 to 2°C, respectively, between the inner and outer fjords (Bendtsen et al., 2014; Rysgaard et al., 2003; Sejr et al., 2017). However, subsurface and bottom waters (i.e., below the pycnocline) remain stable throughout the year (temperature < to 0.5°C and salinity > 27 PSU throughout the year at 17 m, Sejr et al., 2022).

Such thermal stability of bottom water masses thus provides an opportunity to study the impact of seasonality on bivalve physiology

from a trophic (organic matter quality and availability) and reproductive perspective. Four Arctic bivalve species characterized by different longevity were selected: *Astarte moerchi* (maximum longevity [ML] ≈ 110 years, Olivier et al., 2020), *Hiatella arctica* (ML ≈ 120 years, Sejr et al., 2002), *Musculus discors* (ML ≈ 10 years, Selin, 2010), and *Mya truncata* (ML ≈ 50 years, Sejr & Christensen, 2007). These four bivalve species are all suspension feeders and are among the most abundant mollusks in Young Sound, with densities reaching up to $47 \text{ inds} \cdot \text{m}^{-2}$ for *A. moerchi*, $803 \text{ inds} \cdot \text{m}^{-2}$ for *H. arctica*, $1882 \text{ inds} \cdot \text{m}^{-2}$ for *M. discors*, and $66 \text{ inds} \cdot \text{m}^{-2}$ for *M. truncata* (Bridier, 2020; Sejr et al., 2000).

Little information is available on the reproduction of Arctic bivalve in general, except that spawning periods are usually strongly synchronized with the onset of phytoplankton blooms (Brandner et al., 2017). However, these species can be divided into two distinct groups according to the mode of larval development. Larvae produced by *H. arctica* and *M. truncata* have a pelagic development (with wide larval dispersal), whereas those produced by *A. moerchi* and *M. discors* have a benthic development (i.e., without larval dispersal, Brandner et al., 2017; Larsen et al., 2007). A specific study carried out on *H. arctica* shows that this species spawns between June and October in Young Sound (Veillard et al., 2023). The reproductive cycles of *M. truncata* and *H. arctica* are generally synchronized in the Arctic (Brandner et al., 2017), which could suggest that these two species also spawn at the same time in Young Sound (i.e., from June to the end of October). However, to our knowledge, there is no

<figure><img src="image_3.png" /><figcaption>FIGURE 1 Map of the study site showing the location of the Young Sound-Tyrolerfjord system and the Basalt Island sampling station.</figcaption></figure>
information on the reproductive strategies of *M. discors* and *A. moerchi* in the literature.

## 2.2 | Sampling

Bivalves were collected at the Basalt Island station (depth = 21.5 m, 74.33°N, 20.36°W, see Figure 1) in May 2018 under approximately 2 m of sea ice cover (Laurent Chauvaud, personal observations) by SCUBA divers and in August 2018 (only 2–3 days after sea ice breakup) using a 0.04 m² Van Veen grab. Although winter and summer bivalves were collected using two different methods (divers vs. Van Veen grab), both were conducted randomly (i.e., without selection of a specific area or a bivalve size). As the aim of this study was not to compare species abundances or densities (which would obviously be influenced by the sampling methods), but only to gain new insights into the physiology of Arctic bivalves, we believe that these two sampling methods will not influence our seasonal comparisons of bivalve lipid profiles. Unlike other species, the grab does not allow sufficient penetration into the sediment to collect *M. truncata*. Instead, *M. truncata* specimens were collected by scuba divers during a field campaign conducted in August 2015 (to avoid methodological bias, seasonal variations are not discussed for this species). After collection, samples were stored in the field (Daneborg Station) in a –80°C freezer, then transferred to Canada in a –20°C cooler case and stored in a –80°C freezer at UQAR/ISMER until analysis.

## 2.3 | Lipid analyses

A total of 76 individuals (23 *A. moerchi*, 20 *H. arctica*, 18 *M. discors*, and 15 *M. truncata*) were dissected in the laboratory to collect digestive glands and gills and estimate their wet weight prior to lipid extraction. Digestive gland weights and sizes of dissected individuals are provided in Table S1. Individuals were not sexed prior to lipid analyses. Although sex can control the lipid composition of specific tissues (e.g., gonads, mantle), it generally does not influence the lipid profiles of gills and digestive glands (e.g., Birkely et al., 2003; Fernández-Reiriz et al., 2015). To limit potential life stage bias affecting the interindividual variability in lipid profiles, we decided to analyze only individuals that had reached adult size (i.e., assumed to be sexually mature). In addition, all individuals collected for each species were generally evenly distributed across several size classes to avoid an effect of size and/or age on the interindividual variability in lipid profiles. This effect is moreover known to be negligible compared to interspecies physiological differences in mollusk bivalves (Munro & Blier, 2015).

In the digestive gland, the neutral lipid fraction was used as a proxy for food availability and quality because dietary fatty acids are incorporated into triacylglycerol storage lipids with almost no modification (Couturier et al., 2020; Dalsgaard et al., 2003). For gills, only the polar lipid fraction was recovered and used as a proxy for metabolic rate

and acclimation to cold temperature since gill polar lipids are mostly composed of fatty acids from phospholipids involved in the good physiological functioning of organisms (Couturier et al., 2020; Hazel & Williams, 1990; Pernet et al., 2007). Lipid extractions were performed according to the method of Folch et al. (1957) as modified in Chen et al. (1981). Lipids were extracted by grinding digestive glands and gills in tissue grinders with a dichloromethane/methanol solution (2:1, v:v). Total lipid extracts from digestive glands were weighed prior to separation of lipid fractions for use as a proxy for the nutritional status of the bivalves. Neutral and polar fractions were separated on silica gel columns hydrated with 6% deionized water. The columns were first preconditioned with 10mL of methanol and 10mL of dichloromethane. Dry lipid extracts were retrieved three times in 0.5 mL of dichloromethane (total volume = 1.5mL) to ensure that no lipids remained in the tubes and then loaded onto the upper part of silica columns. Silica columns containing lipid extracts (i.e., neutral and polar fractions) were first washed with 10mL dichloromethane/methanol (98:2, v:v) to elute neutral lipids and then washed with 10mL methanol to elute polar lipids (Marty et al., 1992). Eluted fractions were evaporated under a N2 flow, diluted in a methylation solution (Butylhydroxytoluene/sulfuric acid/toluene), heated at 100°C for 10 min to methylate fatty acids and purified on silica gel columns to remove free sterols (Mejri et al., 2014). Fatty acids were finally separated on an Omegawax 250 capillary column (Supelco) and determined by gas chromatography (Trace GC Ultra, Thermo Scientific) coupled to a mass detector (model ITQ900, Thermo Scientific). Fatty acids were identified using commercial standards (Supelco® 37 component FAME mix, CRM47885, Sigma-Aldrich) or mass spectra for peaks not assigned to a known standard using Xcalibur v 3.1 software (Thermo Scientific).

Two indicators based on the fatty acid composition of the gills were calculated to assess the thermal acclimation of bivalves to cold temperatures and their physiological requirements of EFAs. First, thermal acclimation was assessed by the unsaturation index (UI) based on the equation proposed by Logue et al. (2000):

$$ UI = (1 \times \text{Monoenoics}) + (2 \times \text{Dienoics}) + (3 \times \text{Trienoics}) + (4 \times \text{Tetraenoics}) + (5 \times \text{Pentaenoics}) + (6 \times \text{Hexaenoics}) $$

UI represents the average number of double bonds (i.e., unsaturation) between the carbons of a single fatty acid multiplied by 100 (Hazel & Williams, 1990). An UI of 100 would mean that each fatty acid has on average a single unsaturation, while an UI of 300 would mean that each fatty acid has three double bonds. The UI does not provide a direct measure of membrane fluidity (i.e., permeability of membrane proteins), but it is a useful proxy for comparing membrane fluidity between different species and relating these biochemical differences to differences in metabolism or cold acclimation (Hulbert et al., 2007; Munro & Blier, 2015; Pernet et al., 2007). A species with a higher UI (i.e., higher membrane fluidity) is expected to have a higher metabolism and to be more acclimated to low temperatures. Conversely, a species with a low UI (i.e., low membrane fluidity) is considered to have a lower metabolism and to be less acclimated to low temperatures (Munro & Blier, 2015, Pernet et al., 2007). Second, the potential physiological
requirements of bivalves for EFAs were estimated by calculating an index based on the ratio of fatty acids present in gill polar lipids (structural lipids) to the fatty acids from total lipids of food sources (usually particulate organic matter [POM]), as described by the following formula (example with fatty acid 20:ω3):

$$ \text{Gill: Diet ratio} = \frac{\% \text{ of 20: } \omega 3 \text{ in POM total lipids}}{\% \text{ of 20: } 5\omega 3 \text{ in gill polar fatty acids}} $$

This comparison allows us to assess whether specific fatty acids are selectively retained by the organism to ensure the integrity of cell membranes (Copeman et al., 2002; Gendron et al., 2013; Mejri et al., 2021). If the relative proportion of a specific fatty acid in gill polar lipids is equal to or greater than the dietary level (i.e., Gill:Diet ratio equal to or less than 1), then the bivalve requirement for that specific fatty acid is presumably met. Conversely, if the Gill:Diet ratio of a fatty acid is greater than 1, it means that this specific fatty acid is selectively incorporated by bivalve and thus reveals a potential dietary deficiency. To calculate this ratio, we used data from Bridier et al. (2019) on the fatty acid composition of particulate organic matter filtered by bivalves in Young Sound as a proxy for the lipid profile of food sources consumed by the bivalves sampled in this study. We chose to use only summer fatty acid profiles of gill polar lipids and total lipids of particulate organic matter to calculate Gill:Diet ratios, as it has already been shown in Young Sound that bivalves feed only marginally in winter due to depletion of the trophic environment (Bridier et al., 2019).

## 2.4 | Statistical analyses

Data for total lipid content of digestive glands were analyzed using a univariate two-way PERMANOVA (species and season as factors). The PERMANOVAs were preferred to classical ANOVA tests because they do not require assumptions of normality and homoscedasticity for balanced designs (Anderson & Walsh, 2013). When significant, we used a pairwise PERMANOVA to determine whether the observed differences were significant across all groups (i.e., all

species and/or seasons). Seasonal and interspecies variation in the fatty acid composition of both neutral and polar lipids was visualized using two principal component analyses (PCAs). The PCAs were based on transformed fatty acid data using the Hellinger distance to reduce the weight of rare fatty acids on the ordination (Legendre & Gallagher, 2001). The effects of season, species, and their interactions on the fatty acid composition of neutral and polar lipids were assessed by separate multivariate two-way PERMANOVAs. In addition, univariate two-way PERMANOVAs and pairwise PERMANOVAs were applied for lipid content in digestive glands, lipid indicators (UI; sum of saturated fatty acids (SFAs), MUFA (monounsaturated fatty acids), PUFA, branched fatty acids (BrFAs), EFA, 20:2 & 22:2 NMI (non-methylene-interrupted fatty acids) and unknown PUFA; MUFA/PUFA, and EPA/DHA (i.e., 20:5ω3/22:6ω3) ratios) as well as each fatty acid to test their specific variability between the two seasons and four species. Univariate one-way PERMANOVAs with species as factor were applied for Gill:Diet ratios of specific fatty acids (i.e., 20:5ω3 and 22:6ω3) and lipid indicators (i.e., sum of SFA, MUFA, and PUFA), since Gill:Diet ratios were calculated only for the summer season. All statistical analyses were performed using R (R Core Team, 2019) and the “vegan” package (Oksanen et al., 2019).

# 3 | RESULTS

## 3.1 | Nutritional level status indicator

The winter lipid contents in the digestive glands were 40.0±11.3 mg g⁻¹ for *Mya truncata* (mean±standard error of the mean), 43.5±3.2 mg g⁻¹ for *Astarte moerchi*, 47.5±11.7 mg g⁻¹ for *Hiatella arctica*, and 57.6±10.1 mg g⁻¹ for *Musculus discors* (Figure 2). The summer lipid contents in the digestive glands were 26.1±2.8 mg g⁻¹ for *H. arctica*, 30.9±2.8 mg g⁻¹ for *A. moerchi*, 43.3±7.7 mg g⁻¹ for *M. discors*, and 49.0±8.3 mg g⁻¹ for *M. truncata* (Figure 2).

<figure><img src="image_4.png" /><figcaption>FIGURE 2 Seasonal variation in digestive gland lipid contents (mg per g of wet weight) of four bivalve species. Data are mean±Standard error of the mean (SEM). Letters indicate significant interspecific variation in lipid content (winter and summer seasons combined), asterisks indicate significant seasonal variation between lipid content of a same species, and N is the number of replicates per species and season.</figcaption></figure>
The total lipid content of the digestive glands showed variability between the two seasons (degrees of freedom for treatment and error [DF]=1 and 47, pseudo-F=1931.92, p-perm=.007) and between the four species (DF=3 and 47, pseudo-F=902.68, p-perm=.022) without interactions between these two factors (DF=3 and 47, pseudo-F=384.61, p-perm=.216). The levels were significantly higher in the digestive glands of *M. discors* than in the glands of *A. moerchi* and *H. arctica* (*p* < .05) while *M. truncata* showed intermediate values (Figure 2, Table S1). A seasonal effect was observed in the digestive gland from *A. moerchi* and *H. arctica* decreasing by more than 35% between winter and summer (*p* < .05) while no seasonal effect was observed for *M. discors* and *M. truncata* (Figure 2).

## 3.2 | Neutral lipid profiles in digestive glands

Overall, neutral lipid profiles showed high levels of PUFAs, ranging from 40.5±3.0% (*A. moerchi*) to 46.6±0.9% (*M. discors*) in winter, and from 43.0±3.0% (*H. arctica*) to 47.6±1.3% (*M. truncata*) in summer (Table S2). These high levels of PUFAs are mainly explained by the high levels of 20:5ω3 and 22:6ω3. Levels of 20:5ω3 varied from 2.7±0.5% (*M. truncata*) to 4.8±0.3% (*H. arctica*) in winter and from 3.4±0.3% (*M. truncata* and *M. discors*) to 5.6±0.5% (*H. arctica*) in summer (Table S2). Levels of 22:6ω3 varied from 24.6±3.0% (*A. moerchi*) to 33.8±2.7% (*M. truncata*) in winter and from 24.4±1.1% (*H. arctica*) to 37.0±0.5% (*M. truncata*) in summer (Table S2).

The two-way PERMANOVA showed no variation between the two sampled seasons (DF=1 and 62, pseudo-F=1.3410, p-perm=.2367) and no interactions between species and season

factors (DF=3 and 62, pseudo-F=1.6413, p-perm=.1184). However, significant differences were observed between the four species (DF=3 and 62, pseudo-F=12.1583, p-perm=.0001). According to pairwise comparisons, lipid profiles in the neutral fraction of digestive glands were significantly different between all species (*p* < .05). The PCA showed that 36.7% of the total fatty acid variation was explained by the first principal component (Figure 3a,b). The first principal component discriminated the lipid profiles of *A. morchi* and *H. arctica* from those of *M. discors* and *M. truncata* while the second principal component (explaining 16.2% of the total variation) distinguished *A. moerchi* and *M. discors* from *H. artica* and *M. truncata* (Figure 3a). Overall, the fatty acid composition of the neutral fraction of digestive glands from all species was first dominated by a high relative contribution of 20:5ω3 (ranging from 24.4±1.1% [*H. arctica*] to 37.0±0.5% [*M. truncata*]), followed by 16:1ω7 (ranging from 13.7±1.0% [*H. arctica*] to 24.6±1.0% [*M. discors*]), then 16:0 (ranging from 10.0±1.6% [*M. truncata*] to 14.1±0.6% [*H. arctica*]) (Table S2).

## 3.3 | Polar lipid profiles in gills

Overall, the polar lipid profiles showed high levels of PUFAs, ranging from 40.1±3.2% (*A. moerchi*) to 61±1.6% (*H. arctica*) in winter and from 39.0±3.0% (*A. moerchi*) to 63.7±1.1% (*M. truncata*) in summer (Table S3). These high levels of PUFAs are mainly explained by the high levels of 20:2 NMI, 20:5ω3, and 22:6ω3. Levels of 20:2 NMI ranged from 8.0±0.6% (*A. moerchi*) to 12.0±0.8% (*H. arctica*) in winter and from 7.8±0.7% (*A. moerchi*) to 13.0±0.9% (*M. truncata*) in summer (Table S3). Levels of 20:5ω3 varied from 3.7±1.2% (*A. moerchi*) to 14.9±1.9% (*M. truncata*) in winter and from 3.9±0.8% (*A. moerchi*) to 11.2±0.7% (*M. discors*) in summer

<figure><img src="image_5.png" /><figcaption>FIGURE 3 Principal component analysis (PCA) based on Hellinger-transformed neutral fatty acid percentages of bivalve digestive glands from winter (W) and summer (S) seasons (a: individual factor map, b: variable factor map). Variable factor map includes only the 15 most discriminant fatty acids.</figcaption></figure>
(Table S3). Finally, levels of 22:6ω3 varied from 3.6±1.2% (*A. moerchi*) to 15.2±0.8% (*H. arctica*) in winter and from 3.0±0.8% (*A. moerchi*) to 16.2±0.5% (*M. truncata*) in summer (Table S3).

The two-way PERMANOVA revealed the presence of interactions between species and season factors (DF=3 and 68, pseudo-F=2.717, *p*-perm=.0041). A pairwise PERMANOVA showed that all lipid profiles of each factor combination were different from each other, except for winter and summer lipid profiles of *A. moerchi* and *M. discors*, respectively. Most of the variation in lipid profiles was observed on the first principal component, which discriminated *A. moerchi* from the other three species (Figure 4a). When the two principal components were combined, the lipid profiles of *M. truncata* were positioned between the lipid profiles of *A. moerchi* and those of *H. arctica* and *M. discors* (Figure 4a). Overall, *A. moerchi* was characterized by lower levels of PUFAs (especially 22:6ω3 and 20:5ω3) and higher levels of MUFAs (e.g., 20:1ω11; Figure 4b).

These differences in fatty acid profiles were confirmed by univariate two-way PERMANOVAs realized separately for each specific fatty acid (Table S3). The UI and the sum of PUFAs displayed lower values for *A. moerchi* than for *H. arctica*, *M. discors*, and *M. truncata* (*p* < .05, Figure 5a, Table S3). These differences were mainly attributed to 22:6ω3, which displayed a twofold lower relative contribution in polar fatty acid profiles of *A. moerchi* than those of *M. discors* and a fourfold lower relative contribution than those observed for *H. arctica* and *M. truncata* (Figure 5b, Table S3). An opposite trend was observed for 20:1ω11, with *A. moerchi* showing the highest relative contribution (Figure 5c, Table S3). An intermediate relative contribution was observed for *M. truncata* while *H. arctica* and *M. discors* showed the lowest levels, approximately 10-fold lower than those of *A. moerchi* (Figure 5c, Table S3).

Finally, the sum of the 20:2 and 22:2 NMI showed higher values for *A. moerchi*, intermediate values for *M. discors*, and lower values for *H. arctica* and *M. truncata* (the latter characterized by the absence of 22:2) (Figure 5d, Table S3).

## 3.4 | Specific fatty acid ratio from polar gill lipids to diet

For each specific fatty acid tested, the polar lipid fraction from gill to diet ratios showed variable values among the bivalves (Figure 6). Overall, *Astarte moerchi* showed lower Gill:Diet ratios for 20:5ω3 (Gill: Diet ratio=1.5, DF=3 and 37, pseudo-F=23.409, *p*-perm=.0001) and 22:6ω3 (Gill: Diet ratio=3.4, DF=3 and 37, pseudo-F=42.27, *p*-perm=.0001) than the other three species (Gill: Diet ratios >3 and 11 for 20:5ω3 and 22:6ω3, respectively). The same trend was observed for the sum of PUFAs as *A. moerchi* showed the lowest Gill:Diet ratios (DF=3 and 37, pseudo-F=16.922, *p*-perm=.0001). However, the variations among all species were less important, with 25% interspecies variation. The Gill:Diet ratio for MUFA remained around 1 for all species except *A. moerchi*, which had twice as much MUFA in its gills as in particulate organic matter (DF=3 and 37, pseudo-F=49.049, *p*-perm=.0001).

## 4 | DISCUSSION

With climate change, Arctic coastal ecosystems will experience unprecedented environmental changes in sea ice extent and thickness, temperature, hydrological cycle, and stratification (Wassmann et al., 2020). All of these changes are likely to affect

<figure><img src="image_6.png" /><figcaption>**FIGURE 4** Principal component analysis (PCA) based on Hellinger-transformed polar fatty acid percentages of bivalve gills from winter (W) and summer (S) seasons (a: individual factor map, b: variable factor map). Variable factor map includes only the 15 most discriminant fatty acids.</figcaption></figure>
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

FIGURE 5 Seasonal and interspecific variations in unsaturation index (a), relative contributions of 22:6ω3 (b), 20:1ω11 (c) and sum of 20:2 and 22:2 NMI fatty acids (d) in polar lipids of four bivalve species. Data are mean ± Standard error of the mean (SEM). Letters indicate significant interspecific variations (winter and summer seasons combined) while asterisks indicate significant seasonal variations between lipid contents of a same species.
<figure><img src="image_11.png" /><figcaption>**FIGURE 6** Ratio of polar fatty acids of gills to total fatty acids of particulate organic matter. The dashed line (ratio = 1) indicates equal amounts of fatty acids in gills and particulate organic matter. A Gill:Diet ratio for a specific fatty acid greater than 1 indicates selective retention of that specific fatty acid by the organism while a ratio lower than 1 reflects the absence of a physiological role for that fatty acid for the organism. Σ SFA, Σ MUGA, and Σ PUFA correspond to the sum of saturated, monounsaturated, and polyunsaturated fatty acids, respectively. The 20:5ω3 and 22:6ω3 represent the two fatty acids most involved in the metabolism and regulation of membrane fluidity in bivalves. Error bars represent the standard error of the mean (SEM).</figcaption></figure>

primary production and the quality of organic matter sources (e.g., PUFA content) in these low-productivity ecosystems (Ardyna & Arrigo, 2020; Bridier et al., 2019; Michel et al., 2015). However, the lack of fundamental knowledge on the physiology of benthic invertebrates, which are essential links in Arctic marine food webs (Grebmeier, 2012), precludes a proper assessment of their vulnerability to these future qualitative and quantitative changes in the trophic environment. Considering the importance of PUFAs for the physiology of benthic invertebrates, especially for cold acclimation and maintenance of membrane fluidity (Hazel, 1995; Pernet et al., 2006; Sinensky, 1974), we believe it is essential to assess their fatty acid requirements to better evaluate their physiological vulnerability to future trophic changes. To address this knowledge gap, we investigated in this study the role of longevity on the fatty acid requirements of four bivalve species.

## 4.1 | Feeding indicators

On both sampling occasions, *Musculus discors* showed at least 20% more lipid content in their digestive glands than the other species, suggesting a higher accumulation of food possibly related to increased feeding activity. Surprisingly, this species did not show significant seasonal differences suggesting a possible feeding behavior during the winter season. However, while such a behavior may have an energetic advantage (the species can assimilate lipid, protein, or carbohydrate compounds), this is not the case from a physiological

point of view. EFA levels are probably too low at this time of year to be retained by the organisms to meet physiological requirements and allow the optimal control of membrane fluidity (and its effects on metabolism or cold acclimation). This was observed by Bridier et al. (2019), as the levels of EFAs (20:5ω3 and 22:6ω3) in particulate organic matter were so low in winter that they were below detection or quantification thresholds.

Seasonal variability was observed in *Astarte moerchi* and *Hiatella arctica*, which displayed an ~35% decrease in lipid content between winter and summer. The seasonal decrease in lipid content observed for *A.moerchi* and *H.arctica* may reflect (1) a strong allocation of energy reserves (i.e., egg production and spawning) between these two sampling periods toward reproduction (Amiraux et al., 2021; Birkely et al., 2003; Fokina et al., 2018; Pazos et al., 1997) and/or (2) high mobilization of energy reserves to meet the energetic needs of bivalves in the face of low food availability and/or interruption of bivalve feeding during winter (Ahn et al., 2003; Mayrand et al., 2017; McMeans et al., 2015). The mobilization of energy reserves for egg production usually involves a transfer of EFAs (especially 20:5ω3 and 22:6ω3) from the digestive glands to the gonads, resulting in a drastic decrease of these fatty acids in the digestive glands after reproductive investment (Leroy et al., 2013; Martínez-Pita et al., 2012). The absence of such seasonal variations in our digestive gland lipid profiles suggests that no reproductive investment has occurred between our two sampling seasons. Although there is no information to our knowledge on the reproductive cycles of *A.moerchi* and *M.discors* in the literature, observations on *Hiatella arctica* and *Mya truncata* suggest that these two species spawn between June and October in Young Sound (Brandner et al., 2017; Veillard et al., 2023). Considering the absence of seasonal variations in digestive gland lipid profiles, it is possible that spawning occurred later (e.g., September or October), or was even completely suppressed this year if the bivalves had too low energy reserves to sustain this reproductive effort (Amiraux et al., 2021). Future changes in primary production (both in quantity and in quality) could therefore have dramatic consequences for the population dynamics of some bivalve species if the organisms do not have sufficient energy reserves to support the reproductive effort each year.

Seasonal differences seem more likely to reflect a mobilization of energy reserves during winter to meet the energy needs of bivalves. Considering the short period of primary production in Young Sound (restricted to June–September, Holding et al., 2019; Rysgaard et al., 1999) it is likely that *A.moerchi* and *H.arctica*, unlike *M.dis-cors*, are unable to feed for most of the year. Although fjord-shelf exchange may bring organic matter into the fjord from primary production occurring in the open ocean (which is ice-free in May, Boone et al., 2017), the low quality of organic matter (i.e., high dominance of SFA) during this period suggests that this food source is primarily refractory (Bridier et al., 2019). Sea ice breakup occurred very late for Young Sound during our study (August 2, 2018, i.e., the latest observed breakup in a century, Greenland Ecosystem Monitoring database), as sea ice retreat usually occurs around mid-July (Holding et al., 2019). Considering that primary production in Young Sound
peaks 1-2 weeks prior to sea ice breakup (Rysgaard et al., 1999), the lipid content measured in our study is likely to be more depleted than would typically be observed at the same time in other years (e.g., Amiraux et al., 2021). However, these potential interannual differences remain small compared to many other Arctic sites with much more regular access to fresh organic matter throughout the year (e.g., Renaud et al., 2011). Our results thus highlight the unique trophic conditions in Young Sound where primary consumers have only 2-3 months to replenish their energy reserves for the winter period, as previously suggested by Sejr and Christensen (2007). In all bivalves collected, the dominant fatty acids in the neutral lipid fraction of the digestive gland were the 16:1ω7 and 20:5ω3 during the two seasons confirming their primary feeding on pelagic and/or benthic diatoms (Bridier et al., 2021). *M. discors* seems to be better adapted to ingest the available food, as its lipid content was systematically higher than that of other species without a significant decrease during winter conditions.

## 4.2 | Physiological indicators

The fatty acid composition of polar lipids did not vary seasonally for *Astarte moerchi*, *Hiatella arctica*, and *Musculus discors*. The absence of seasonal variation in the fatty acid composition of gills supports the homeoviscous adaptation theory (Hazel, 1995; Sinensky, 1974), which has been observed in marine bivalves (e.g., Hall et al., 2002; Parent et al., 2008; Pernet et al., 2007). To counteract the effect of seasonal temperature variations on cell membrane fluidity, marine bivalves undergo a significant lipid remodeling between winter and summer by increasing the proportion of MUFA to SFA (i.e., decreasing the SFA/ MUFA ratio, Munro & Blier, 2012) and/or by increasing the proportion of PUFAs (especially in 20:5ω3 and 22:6ω3, Hall et al., 2002; Hazel, 1995; Pernet et al., 2007). However, all these studies were mainly conducted in subarctic or temperate or tropical systems, characterized by moderate or strong seasonal temperature variations (i.e., thermal amplitudes >6°C, Mathieu-Resuge et al., 2020). The relatively stable temperatures in Young Sound (i.e., thermal amplitude <2.5°C at 17 m, Sejr et al., 2022) suggest that this temperature range is insufficient to induce a seasonal remodeling of lipid membranes. However, the fatty acid composition of polar lipids observed throughout the year likely reflects the acclimation of organisms to the low temperatures observed in the fjord (temperature <to 0.5°C throughout the year at 17 m, Sejr et al., 2022). The strong relative contributions of PUFAs observed in *H. arctica*, *M. discors*, and *Mya truncata* (i.e., >47%) are similar to those observed in other studies that also found an effect of cold temperature on the lipid composition of cell membranes of Arctic/sub-Arctic bivalve species (Gaillard et al., 2015; Hacker Teper et al., 2022; Thyrring et al., 2017).

In contrast to the other three species, *A. moerchi* is characterized by much lower levels of PUFAs, especially 22:6ω3 (i.e., two to four times lower than the other species). Such low levels probably do not reflect a failure of *A. moerchi* to acclimate to cold temperatures but rather highlight the physiological characteristics of this

species. According to the *membrane-pacemaker theory of aging*, the differences in UI and relative contribution of 22:6ω3 observed between *A. moerchi*, *M. discors*, and *M. truncata* suggest that *A. moerchi* may have a lower metabolic rate than the other two species (Hulbert, 2005; Hulbert et al., 2007). Such a difference could reflect the longer longevity of *A. moerchi*, as long-lived species typically have lower metabolic rates (associated with a lower production of reactive oxygen species) and modify the lipid composition of cell membranes by preferentially selecting fatty acids with the lowest number of double bonds to avoid peroxidation of membrane lipids (Blier et al., 2017). The higher levels of 20:2 and 22:2 NMI observed in gill lipid profiles of *A. moerchi* also support this hypothesis, as these fatty acids are known to be more resistant to oxidative processes than other PUFAs (Barnathan, 2009). A previous study highlighted the ability of *Astarte borealis* to significantly reduce its metabolism under anoxic conditions to less than 0.5% of its aerobic rate (Oeschger, 1990). One might expect that such a metabolic strategy, known as metabolic rate depression (Abele et al., 2009), might also occur in *A. moerchi* to cope with other stressful environmental conditions such as low temperature and/or poor trophic conditions, as is the case in Young Sound.

Surprisingly, *H. arctica* considered with long lifespan (i.e., ML ≈ 120 years, Sejr et al., 2002) showed similar lipid dynamics in their cell membranes of gills than short-lived species (i.e., similar UI and PUFA level). Thus, *H. arctica* could have metabolic rate at similar levels than *M. truncata* and *M. discors* and probably different from the potentially reduced metabolic rate of long-lived *A. moerchi* showing the lowest levels of PUFAs and 22:6ω3 in their gill membranes. The digestive glands of *H. arctica* also show the greatest seasonal variation in lipid content suggesting that this species may mobilize its energy reserves during winter to potentially cover the metabolic costs. Such results may indicate that *H. arctica* does not enter a metabolically depressed state like *A. moerchi* might when environmental conditions are stressful.

One hypothesis can be suggested to explain why two long-lived species (i.e., *A. moerchi* and *H. arctica*) adopt different fatty acid dynamics in gill membranes during winter. It may be that cold acclimation and metabolism act in opposite directions in regulating the fatty acid composition of cell membranes. If bivalves regulate the cell membrane fluidity by increasing the proportion of PUFAs and especially 22:6ω3 to adapt to cold (i.e., homeoviscous adaptation, e.g., Pernet et al., 2006), such a lipid remodeling would also inevitably lead to an increase in metabolic rate (according to the membrane pacemaker theory of metabolism, Hulbert, 2005). We can therefore speculate that *A. moerchi* and *H. arctica* use different strategies to acclimate to cold temperature. Homeoviscous adaptation does not necessarily involve a mere increase in the proportion of PUFAs but may also involve additional mechanisms such as an increase in the proportion of MUFAs, a decrease in cholesterol content, or an elongation of fatty acid carbon chains (Hazel, 1995; Pernet et al., 2006; Sinensky, 1974). *H. arctica* may mainly regulate the fluidity of cell membranes by increasing the proportion of PUFAs while *A. moerchi* may favor some of these additional mechanisms to acclimate
to cold temperature without affecting the metabolic rate. The low relative contribution of PUFAs in *A.merchi* cell membranes combined with the high relative contribution of de novo synthetized MUFAs with long carbon chains (e.g., 20:1ω11) or NMI fatty acids (which also ensure good membrane fluidity with less unsaturation than PUFAs, Barnathan, 2009; Kraffe et al., 2004) seems to confirm this hypothesis. Therefore, maintaining high levels of PUFAs in cell membranes for cold acclimation would inevitably lead to an increase in metabolism and could ultimately result in a mismatch between the high energy needs of *H.arctica* and the limited lipid reserves of long-lived species. This hypothesis seems to be confirmed by the seasonal dynamics of digestive gland lipid content, since *H.arctica* has both the lowest total lipid content and the highest seasonal variation among the four species. On the other hand, *A.merchi* might maintain a reduced metabolism and thus better match the size of its energy reserves with the rate of their depletion by favoring other cold acclimation mechanisms such as a decrease in the SFA/MUFA ratio (Munro & Blier, 2012) or decrease in the polar lipid/sterol ratio (i.e., an increase in the proportion of sterol compared to the proportion of polar lipids, Parent et al., 2008). Such contrasting physiological strategies for regulating membrane fluidity between two long-lived species may reflect differences in their geographic distribution. While *A.merchi* is restricted to Arctic waters (De Cesare et al., 2017; Olivier et al., 2020; Petersen, 2001), *H.arctica* has a broad geographic distribution from temperate latitudes (e.g., Azores, Mediterranean Sea) to the Arctic Ocean (Cardigos et al., 2006; Hunter, 1949; Manousis, 2021). The biogeography of *H.arctica* may therefore imply that this species is not as well adapted as *A.merchi* to cold Arctic environments. It is possible that biogeography also explains some of the differences observed between the lipid profiles of *A.merchi* and those of *M.discors* and *M.truncata*, since the latter two species are also characterized by a wider geographic distribution than *A.merchi* (Günther & Fedyakov, 2000; Petersen, 1978). We believe that the role of biogeography in these interspecific physiological differences is probably small compared to the role of longevity, since longevity is known to be the main factor controlling the fatty acid composition of cell membranes (Munro & Blier, 2015). However, it would be interesting for future studies to attempt to separate the confounding effects of biogeography and longevity by selecting circumpolar, boreal and/or cosmopolitan species with different longevities within each geographic distribution.

## 4.3 | Selective retention of EFAs

The variation in Gill:Diet ratios among the four species reflects the trends observed in the fatty acid composition of polar lipids (see previous section). All bivalves have a Gill:Diet ratio for total PFAs (Σ MUFA) higher than 1, reflecting the need for these species to selectively retain PUFAs from their diet to maintain the integrity of their cell membranes (Gendron et al., 2013; Mejri et al., 2021). Such incorporation of PUFAs into membrane lipids may be partly explained by the physiological requirements of bivalves for highly

unsaturated fatty acids with long carbon chains (e.g., 20:5ω3 and 22:6ω3), which are crucial for maintaining cell membrane fluidity in the face of cold (i.e., homeoviscous adaptation, Pernet et al., 2006). However, as feeding regime was dominated by 20:5ω3 from diatoms, their selective retention is largely lower than that observed for 22:6ω3, which is poorly represented in the diet. Some variation in Gill:Diet ratios is also observed among the four bivalve species. *A.merchi* has significantly lower ratios for PUFAs, especially for EFAs such as 20:5ω3 and 22:6ω3 which are two to six times less selectively retained in *A.merchi* cell membranes than in the other three species. As discussed above, these differences probably highlight the physiological characteristics of *A.merchi* and in particular its long lifespan, which may be related to the lower proportion of PUFAs in lipid membranes to reduce oxidative stress (Blier et al., 2017) and maintain a slower metabolism (Hulbert, 2005; Hulbert et al., 2007). In contrast to the other species, the Σ MUFA Gill:Diet ratio is twice as high in *A.merchi* (Gill: Diet ratio=2.0), mainly due to 20:1ω11 (representing about one third of the total relative contribution of MUFAs). However, this probably does not reflect a selective retention of MUFAs by *A.merchi* but rather reflects de novo synthesis of 20:1ω11 to adjust membrane fluidity (Kraffe et al., 2004; Mathieu-Resuge et al., 2020).

## 4.4 | Need for new knowledge of benthic invertebrate physiology on a pan-Arctic scale

The results discussed above provide some insight on the vulnerability of benthic invertebrates to change in Arctic ecosystems. However, it should be kept in mind that our results come from the study of a single Arctic region and a single study site (i.e., Young Sound, North-East Greenland). Therefore, it is important that a larger number of future studies focus on acquiring new fundamental information on the physiology of benthic invertebrates in order to accurately assess their vulnerability to climate change (Renaud et al., 2015, 2019). The comparison of several stations distributed along a trophic gradient (e.g., a gradient extending from the inner to the outer part of a fjord, Bridier et al., 2019) could help to better assess the degree of importance of interspecies physiological differences in the resilience/vulnerability of invertebrates to future trophic changes. The acquisition of new data on the physiology of benthic invertebrates in different Arctic regions would also allow a better understanding of how these physiological acclimations may vary locally according to the local environmental characteristics (latitude, depth, distance from shore, etc.). Only the synthesis of these future studies will provide a comprehensive understanding of the vulnerability of benthic invertebrates to future trophic changes on a pan-Arctic scale.

# 5 | CONCLUSIONS AND OUTLOOK

Arctic coasts will experience multiple environmental changes related to cryosphere loss that will subsequently affect overall ecosystem
FUTURES TROPHIC CHANGES IN ARCTIC COASTAL ECOSYSTEMS :
YOUNG SOUND STUDY CASE

<figure><img src="image_12.png" /><figcaption>FIGURE 7 A schematic attempting to summarize future trophic constraints in Young Sound (in terms of primary producer phenology and organic matter quantity/quality) with the expected physiological capabilities of the four studied bivalve species to cope with these future trophic changes.</figcaption></figure>

functioning. Increased glacier melt and meltwater inputs to coastal
ecosystems are likely to change the phenology of primary producers,
modify the ecosystem productivity, and alter the nutritional qual-
ity of organic matter sources available to primary consumers. New
fundamental knowledge about the physiology and nutritional needs
of Arctic benthic invertebrates is therefore essential if we want to
properly assess their vulnerability to future environmental changes
in coastal ecosystems. Fatty acids are a powerful tool that can pro-
vide information on the spatiotemporal changes in the source, qual-
ity, and quantity of food sources, as well as reflect physiological
responses to environmental stress. The strong disparities observed
among the four bivalve species in their fatty acid dynamics and re-
modeling in relation to their nutritional needs suggest that not all
benthic invertebrates will be equally vulnerable to future trophic
changes (Figure 7). Species characterized by long life expectancy
and evidence of low metabolism (e.g., *A. moerchi*) are likely to be
best adapted to future trophic conditions, as their low EFA require-
ments (e.g., 20:5ω3 and 22:6ω3) may help them to cope with a de-
crease in organic matter quality. In contrast, species with a shorter
life expectancy (e.g., *Musculus discors* and *Mya truncata*) may have a
higher energetic turnover and thus a higher demand for EFAs. On
the other hand, more regular inputs of fresh organic matter (due to
longer ice-free period and growing season for phytoplankton) may
help these species to meet their high metabolic needs during winter
(due to the absence of winter dormancy) without depleting their en-
ergy reserves. Finally, species characterized by long life expectancy
and evidence of high metabolic rate (e.g., *Hiatella arctica*) may be the

least adapted to these future trophic conditions because they have
both high EFA and energy requirements (due to their high metabolic
rates) but limited energy reserves.

These interspecific differences in vulnerability to future envi-
ronmental change may have broader consequences for the func-
tioning of Arctic coastal ecosystems. Invertebrates with short
life spans and high metabolic rates also have higher growth and
productivity, and therefore represent optimal carbon transfer
pathways within Arctic food webs. This is the case for *M. truncata*,
which is an essential link in the transfer of organic matter pro-
duced by phytoplankton or microphytobenthos to higher trophic
levels (e.g., walruses and eiders, Born et al., 2003; Frimer, 1997).
Conversely, species with long life spans and low metabolic rates
(e.g., *A. moerchi*) have lower growth and productivity, and there-
fore contribute less to carbon transfers within trophic webs.
Consequently, the increased vulnerability of high metabolic rate
species to future trophic changes may have further consequences
for coastal Arctic food webs by reducing carbon flows to higher
trophic levels.

AUTHOR CONTRIBUTIONS

**Guillaume Bridier:** Conceptualization (equal); formal analysis (equal); investigation (equal); methodology (equal); software (equal); visualization (equal); writing – original draft (equal).

**Frédéric Olivier:** Conceptualization (equal); funding acquisition (equal); project administration (equal); supervision (equal); writing – review and editing (equal). **Jacques Grall:** Conceptualization
(equal); funding acquisition (equal); project administration (equal); supervision (equal); writing - review and editing (equal). **Laurent Chauvaud:** Conceptualization (equal); funding acquisition (equal); project administration (equal); supervision (equal); writing - review and editing (equal). **Mikael K. Sejr:** Conceptualization (equal); funding acquisition (equal); project administration (equal); writing - review and editing (equal). **Réjean Tremblay:** Conceptualization (equal); funding acquisition (equal); project administration (equal); resources (equal); supervision (equal); writing - review and editing (equal).

## ACKNOWLEDGMENTS

We thank Egon Randa Frandsen, Jakob Thyrring, Mie Sichlau Winding, Carl Isaken, and Erwan Amice for their help in the field collecting samples. We are very grateful to the MarineBasis program and the staff at the Zackenberg and Daneborg stations for their logistical support in Greenland. Thanks to Jean-Bruno Nadalini for his invaluable help in the UQAR-ISMER laboratory. We also thank Mathieu Babin for carrying out the analyses by gas chromatography-mass spectrometry and Thibault Androuain for FA identification. This work was a part of the PhD project of Guillaume Bridier, co-funded by the "Allocation de Recherche Doctorale" from the Brittany Regional Council and the Université de Bretagne Occidentale/LIA BeBEST, and continued during the postdoctoral fellowship of Guillaume Bridier funded by the Université du Québec à Rimouski. Additional fundings were provided by the Institut Polaire Français Paul-Emile Victor (IPEV-PRIVARC project), the European H2020 INTAROS program, the Agence Nationale de la Recherche (ANR GAAP), the Fondation UBO (FAC 2018), and the Observatoire Marin de l'IUEM (UMS3113). We would also like to thank LabexMer for funding a three-month mobility grant to UQAR-ISMER to carry out laboratory analyses. We would like to warmly thank the two anonymous reviewers and the associate editor for their valuable comments, which greatly improved this manuscript. This study is an international collaboration between the BeBEST International Laboratory (UBO, UQAR, CNRS, and MNHN) and the Aarhus University.

## DATA AVAILABILITY STATEMENT

The data that support the findings of this study are available in the Supporting Information. Authors are available upon request for additional information on data analyses.

## ORCID

Guillaume Bridier ID https://orcid.org/0000-0002-6746-8818

Frédéric Olivier ID https://orcid.org/0000-0002-5776-8761

## REFERENCES

Abele, D., Brey, T., & Philipp, E. (2009). Bivalve models of aging and the determination of molluscan lifespans. *Experimental Gerontology, 44*, 307-315.

Ahn, I. Y., Surh, J., Park, Y. G., Kwon, H., Choi, K. S., Kang, S. H., Choi, H. J., Kim, K. W., & Chung, H. (2003). Growth and seasonal energetics of the Antarctic bivalve *Laternula elliptica* from King George Island, Antarctica. *Marine Ecology Progress Series*, 257, 99-110.

AMAP (Arctic Monitoring and Assessment Programme). (2017). *Snow, water, ice and permafrost in the Arctic (SWIPA)*. Arctic Monitoring and Assessment Programme (AMAP).

Amiraux, R., Archambault, P., Moriceau, B., Lemire, M., Babin, M., Memery, L., Massé, G., & Tremblay, J. E. (2021). Efficiency of symmetric-benthic coupling revealed by analyses of n-3 fatty acids, IP25 and other highly branched isoprenoids in two filter-feeding Arctic benthic molluscs: *Mya truncata* and *Serripes groenlandicus*. *Organic Geochemistry*, 151, 104160.

Anderson, M. J., & Walsh, D. C. (2013). PERMANOVA, ANOSIM, and the mantel test in the face of heterogeneous dispersions: What null hypothesis are you testing? *Ecological Monographs*, 83, 557-574.

Ardyna, M., & Arrigo, K. R. (2020). Phytoplankton dynamics in a changing Arctic Ocean. *Nature Climate Change*, 10, 892-903.

Ardyna, M., Babin, M., Gosselin, M., Devred, E., Rainville, L., Rainville, L., & Tremblay, J. É. (2014). Recent Arctic Ocean sea ice loss triggers novel fall phytoplankton blooms. *Geophysical Research Letters*, 41, 6207-6212.

Arendt, K. E., Nielsen, T. G., Rysgaard, S., & Tönnesson, K. (2010). Differences in plankton community structure along the Godthåbsfjord, from the Greenland ice sheet to offshore waters. *Marine Ecology Progress Series*, 401, 49-62.

Barnathan, G. (2009). Non-methylene-interrupted fatty acids from marine invertebrates: Occurrence, characterization and biological properties. *Biochimie*, 91, 671-678.

Bendtsen, J., Gustafsson, K. E., Rysgaard, S., & Vang, T. (2007). Physical conditions, dynamics and model simulations during the ice-free period of the young sound/Tyrolerfjord system. In S. Rysgaard & R. N. Glud (Eds.), *Carbon cycling in Arctic marine ecosystems: Case study young sound* (Vol. 58, pp. 46-59). Meddr. Grønland. Bioscience.

Bendtsen, J., Mortensen, J., & Rysgaard, S. (2014). Seasonal surface layer dynamics and sensitivity to runoff in a high Arctic fjord (young sound/Tyrolerfjord, 74°N). *Journal of Geophysical Research, Oceans*, 119, 6461-6478.

Birkely, S. R., Grahl-Nielsen, O., & Gulliksen, B. (2003). Temporal variations and anatomical distributions of fatty acids in the bivalve *Mya truncata*, L. 1758, from Isfjorden, Spitsbergen. *Polar Biology*, 26, 83-92.

Blier, P. U., Abele, D., Munro, D., Degletagne, C., Rodriguez, E., & Hagen, T. (2017). What modulates animal longevity? Fast and slow aging in bivalves as a model for the study of lifespan. *Seminars in Cell & Developmental Biology*, 70, 130-140.

Boone, W., Rysgaard, S., Kirillov, S., Dmitrenko, I., Bendtsen, J., Mortensen, J., Meire, L., Petrusevich, V., & Barber, D. G. (2017). Circulation and fjord-shelf exchange during the ice-covered period in young sound-Tyrolerfjord, Northeast Greenland (74°N). *Estuarine, Coastal and Shelf Science*, 194, 205-216.

Born, E. W., Rysgaard, S., Ehlmé, G., Sejr, M., Acquarone, M., & Levermann, N. (2003). Underwater observations of foraging free-living Atlantic walruses (*Odobenus rosmarus rosmarus*) and estimates of their food consumption. *Polar Biology*, 26, 348-357.

Brandner, M. M., Stübner, E., Reed, A. J., Gabrielsen, T. M., & Thatje, S. (2017). Seasonality of bivalve larvae within a high Arctic fjord. *Polar Biology*, 40, 263-276.

Bridier, G. (2020). Functioning of Arctic and sub-Arctic shallow benthic ecosystems in highly-stratified coastal systems [PhD dissertation, Université de Bretagne Occidentale]. https://tel.archives-ouvertes.fr/tel-03588863/

Bridier, G., Meziane, T., Grall, J., Chauvaud, L., Sejr, M. K., Menneteau, S., & Olivier, F. (2019). Coastal waters freshening and extreme seasonality affect organic matter sources, quality, and transfers in a high Arctic fjord (young sound, Greenland). *Marine Ecology Progress Series*, 610, 15-31.

Bridier, G., Olivier, F., Chauvaud, L., Sejr, M. K., & Grall, J. (2021). Food source diversity, trophic plasticity, and omnivory enhance the stability of a shallow benthic food web from a high-Arctic fjord
exposed to freshwater inputs. *Limnology and Oceanography, 66*, S259-S272.

Cardigos, F., Tempera, F., Ávila, S., Gonçalves, J., Colaço, A., & Santos, R. S. (2006). Non-indigenous marine species of the Azores. *Helgoland Marine Research, 60*, 160–169.

Chen, I. S., Shen, C. S. J., & Sheppard, A. J. (1981). Comparison of methylene chloride and chloroform for the extraction of fats from food products. *Journal of the American Oil Chemists' Society, 58*, 599–601.

Citterio, M., Sejr, M. K., Langen, P. L., Mottram, R. H., Abermann, J., Larsen, J. H., Skov, K., & Lund, M. (2017). Towards quantifying the glacial runoff signal in the freshwater input to Tyrolerfjord-Young Sound, NE Greenland. *Ambio, 46*, S146-S159.

Copeman, L. A., Parrish, C. C., Brown, J. A., & Harel, M. (2002). Effect of docosahexanoic, eicosapentanoic, and arachidonic acids on the early growth, survival, lipid composition and pigmentation of yellowtail flounder (*Limanda ferruginea*): A live food enrichment experiment. *Aquaculture, 210*, 285–304.

Couturier, L. I. E., Michel, L. N., Amaro, T., Budge, S. M., da Costa, E., De Troch, M., Di Dato, V., Fink, P., Giraldo, C., Le Grand, F., Loaiza, I., Mathieu-Resuge, M., Nichols, P. D., Parrish, C. C., Sardenne, F., Vagner, M., Pernet, F., & Soudant, P. (2020). State of art and best practices for fatty acid analysis in aquatic sciences. *ICES Journal of Marine Science, 77*, 2375–2395.

Dalsgaard, J., John, M. S., Kattner, G., Müller-Navarra, D., & Hagen, W. (2003). Fatty acid trophic markers in the pelagic marine environment. *Advances in Marine Biology, 46*, 225–340.

De Cesare, S., Meziane, T., Chauvaud, L., Richard, J., Sejr, M. K., Thébault, J., Winkler, G., & Olivier, F. (2017). Dietary plasticity in the bivalve *Astarte moerchi* revealed by a multimarker study in two Arctic fjords. *Marine Ecology Progress Series, 567*, 157–172.

Delpech, L. M., Vonnahme, T. R., McGovern, M., Gradinger, R., Przemel, K., & Poste, A. E. (2021). Terrestrial inputs shape coastal bacterial and archaeal communities in a high Arctic Fjord (Isfjorden, Svalbard). *Frontiers in Microbiology, 12*, 614634.

Fernández-Reiriz, M. J., Garrido, J. L., & Irisarri, J. (2015). Fatty acid composition in *Mytilus galloprovincialis* organs: Trophic interactions, sexual differences and differential anatomical distribution. *Marine Ecology Progress Series, 528*, 221–234.

Fokina, N., Storhaug, E., Bakhmet, I., Maximovich, N., Frantzen, M., & Nahrgang, J. (2018). Seasonal changes in lipid class content in mussels *Mytilus* spp. from Rakkfjorden in the Norwegian Sea and Kandalaksha Bay of the White Sea. *Polar Biology, 41*, 2103–2117.

Folch, J., Lees, M., & Stanley, G. H. S. (1957). A simple method for the isolation and purification of total lipids from animal tissues. *The Journal of Biological Chemistry, 226*, 497–509.

Frimer, O. (1997). Diet of moulting king eiders *Somateria spectabilis* at Disko Island, West Greenland. *Ornis Fennica, 74*, 187–194.

Gaillard, B., Meziane, T., Tremblay, R., Archambault, P., Layton, K. K. S., Martel, A. L., & Olivier, F. (2015). Dietary tracers in *Bathyarca glacialis* from contrasting trophic regions in the Canadian Arctic. *Marine Ecology Progress Series, 536*, 175–186.

Gendron, L., Tremblay, R., Belvin, S., Génard, B., Motnikar, S., & Côté, J. (2013). Condition, survival and growth in situ of hatchery-reared stage IV lobster (*Homarus americanus*) fed Artemia and lipid-rich wild zooplankton. *Aquaculture, 416*, 380–389.

Grebeier, J. M. (2012). Shifting patterns of life in the Pacific Arctic and sub-Arctic seas. *Annual Review of Marine Science, 4*, 63–78.

Günther, C. P., & Fedyakov, V. V. (2000). Seasonal changes in the bivalve larval plankton of the White Sea. *Senckenbergiana Maritima, 30*, 141–151.

Hacker Teper, S., Parrish, C. C., & Gagnon, P. (2022). Multiple trophic tracer analyses of subarctic Rhodolith (*Lithothamnion glaciale*) bed Trophodynamics uncover bottom-up forcing and benthic-pelagic coupling. *Frontiers in Marine Science, 9*, 899812.

Haider, F., Timm, S., Bruhns, T., Noor, M. N., & Sokolova, I. M. (2020). Effects of prolonged food limitation on energy metabolism and burrowing activity of an infaunal marine bivalve, *Mya arenaria*. *Comparative Biochemistry and Physiology. Part A, Molecular & Integrative Physiology, 250*, 110780.

Hall, J. M., Parrish, C. C., & Thompson, R. J. (2002). Eicosapentaenoic acid regulates scallop (*Placopecten magellanicus*) membrane fluidity in response to cold. *The Biological Bulletin, 202*, 201–203.

Hazel, J. R. (1995). Thermal adaptation in biological membranes: Is homeoviscous adaptation the explanation? *Annual Review of Physiology, 57*, 19–42.

Hazel, J. R., & Williams, E. E. (1990). The role of alterations in membrane lipid composition in enabling physiological adaptation of organisms to their physical environment. *Progress in Lipid Research, 29*, 167–227.

Hernes, P. J., Tank, S. E., Sejr, M. K., & Glud, R. N. (2021). Element cycling and aquatic function in a changing Arctic. *Limnology and Oceanography, 66*, S1–S16.

Hochachka, P. W., & Somero, G. N. (2002). *Biochemical adaptation*. Princeton University Press.

Holding, J. M., Markager, S., Juul-Pedersen, T., Paulsen, M. L., Møller, E. F., Meire, L., & Sejr, M. K. (2019). Seasonal and spatial patterns of primary production in a high-latitude fjord affected by Greenland ice sheet run-off. *Biogeosciences, 16*, 3777–3792.

Hulbert, A. J. (2005). On the importance of fatty acid composition of membranes for aging. *Journal of Theoretical Biology, 234*, 277–288.

Hulbert, A. J., Pamplona, R., Buffenstein, R., & Buttemer, W. A. (2007). Life and death: Metabolic rate, membrane composition, and life span of animals. *Physiological Reviews, 87*, 1175–1213.

Hunter, W. R. (1949). The structure and behaviour of *Hiatella gallicana* (lam.) and *H. Arctica* (L.) with special reference to the boring habit. *Proceedings of the Royal Society of Edinburgh B, 63*, 271–289.

Kędra, M., Moritz, C., Choy, E. S., David, C., Degen, R., Duerksen, S., Ellingsen, I., Górska, B., Grebmeier, J. M., Kirievskaya, D., van Oevelen, D., Piwosz, K., Samuelsen, A., & Węsławski, J. M. (2015). Status and trends in the structure of Arctic benthic food webs. *Polar Research, 34*, 23775.

Kraffe, E., Soudant, P., & Marty, Y. (2004). Fatty acids of serine, ethanolamine, and choline plasmalogens in some marine bivalves. *Lipids, 39*, 59–66.

Kwok, R. (2018). Arctic Sea ice thickness, volume, and multiyear ice coverage: Losses and coupled variability (1958–2018). *Environmental Research Letters, 13*, 105005.

Larsen, J. B., Frischer, M. E., Ockelmann, K. W., Rasmussen, L. J., & Hansen, B. W. (2007). Temporal occurrence of planktotrophic bivalve larvae identified morphologically and by single step nested multiplex PCR. *Journal of Plankton Research, 29*, 423–436.

Legendre, P., & Gallagher, E. D. (2001). Ecologically meaningful transformations for ordination of species data. *Oecologia, 129*, 271–280.

Leroy, F., Meziane, T., Riera, P., & Comtet, T. (2013). Seasonal variations in maternal provisioning of *Crepidula fornicata* (Gastropoda): Fatty acid composition of females, embryos and larvae. *PLoS ONE, 8*, e75316.

Lewis, K. M., van Dijken, G. L., & Arrigo, K. R. (2020). Changes in phytoplankton concentration now drive increased Arctic Ocean primary production. *Science, 369*, 198–202.

Logue, J., de Vries, A., Fodor, E., & Cossins, A. (2000). Lipid compositional correlates of temperature-adaptive interspecific differences in membrane physical structure. *The Journal of Experimental Biology, 203*, 2105–2115.

Manousis, T. (2021). The marine Mollusca of Greece: An up-to-date, systematic catalogue, documented with bibliographic and pictorial references. *Xenophora Taxon, 34*, 26–47.

Martínez-Pita, I., Sánchez-Lazo, C., Ruíz-Jarabo, I., Herrera, M., & Mancera, J. M. (2012). Biochemical composition, lipid classes, fatty
acids and sexual hormones in the mussel *Mytilus galloprovincialis* from cultivated populations in South Spain. *Aquaculture*, 358-359, 274-283.

Marty, Y., Delaunay, F., Moal, J., & Samain, J. F. (1992). Changes in the fatty acid composition of *Pecten maximus* during larval development. *Journal of Experimental Marine Biology and Ecology*, *163*, 221-234.

Mathieu-Resuge, M., Le Grand, F., Schaal, G., Lluch-Cota, S. E., Racotta, I. S., & Kraffe, E. (2020). Specific regulations of gill membrane fatty acids in response to environmental variability reveal fitness differences between two suspension-feeding bivalves (*Nodipecten sub nodosus* and *Spondylus crassisquama*). *Conservation Physiology*, *8*, coaa079.

Mayrand, E., Comeau, L. A., & Mallet, A. (2017). Physiological changes during overwintering of the eastern oyster Crassostrea virginica (Gmelin, 1791). *Journal of Molluscan Studies*, *83*, 333-339.

McMeans, B. C., McCann, K. S., Humphries, M., Rooney, N., & Fisk, A. T. (2015). Food web structure in temporally-forced ecosystems. *Trends in Ecology & Evolution*, *30*, 662-672.

Meire, L., Mortensen, J., Meire, P., Juul-Pedersen, T., Sejr, M. K., Rysgaard, S., Nygaard, R., Huybrechts, P., & Meysman, F. J. R. (2017). Marine-terminating glaciers sustain high productivity in Greenland fjords. *Global Change Biology*, *23*, 5344-5357.

Mejri, S. C., Audet, C., Vandenberg, G. W., Parrish, C. C., & Tremblay, R. (2014). Biochemical egg quality in a captive walleye (*Sander vitreus*) broodstock population relative to ovulation timing following hormonal treatment. *Aquaculture*, *431*, 99-106.

Mejri, S. C., Tremblay, R., Audet, C., Wills, P. S., & Riche, M. (2021). Essential fatty acid requirements in tropical and cold-water marine fish larvae and juveniles. *Frontiers in Marine Science*, *8*, 680003.

Michel, C., Hamilton, J., Hansen, E., Barber, D., Reigstad, M., Iacozza, J., Seuthe, L., & Niemi, A. (2015). Arctic Ocean outflow shelves in the changing Arctic: A review and perspectives. *Progress in Oceanography*, *139*, 66-88.

Munro, D., & Blier, P. U. (2012). The extreme longevity of *Arctica islandica* is associated with increased peroxidation resistance in mitochondrial membranes. *Aging Cell*, *11*, 845-855.

Munro, D., & Blier, U. (2015). Age, diet, and season do not affect longevity-related differences in peroxidation index between *Spisula solidissima* and *Arctica islandica*. *The Journals of Gerontology. Series A, Biological Sciences and Medical Sciences*, *70*, 434-443.

Murray, C., Markager, S., Stedmon, C. A., Juul-Pedersen, T., Sejr, M. K., & Bruhn, A. (2015). The influence of glacial melt water on bio-optical properties in two contrasting Greenlandic fjords. *Estuarine, Coastal and Shelf Science*, *163*, 72-83.

Oeschger, R. (1990). Long-term anaerobiosis in sublittoral marine invertebrates from the Western Baltic Sea: *Halicryptus spinulosus* (Priapulida), *Astarte borealis* and *Arctica islandica* (Bivalvia). *Marine Ecology Progress Series*, *59*, 133-143.

Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., Minchin, P. R., O'Hara, R. B., Simpson, G. L., Solymos, P., Stevens, M. H. H., Szoecs, E., & Wagner, H. (2019). *Vegan: Community ecology package*. R package version 2.5-6. Retrieved from https://CRAN.R-project.org/package=vegan

Olivier, F., Gaillard, B., Thébault, J., Meziane, T., Tremblay, R., Dumont, D., Bélanger, S., Gosselin, M., Jolivet, A., Chauvaud, L., Martel, A. L., Rysgaard, S., Olivier, A.-H., Pettré, J., Mars, J., Gerber, S., & Archambault, P. (2020). Shells of the bivalve *Astarte moerchi* give new evidence of a strong pelagic-benthic coupling shift occurring since the late 1970s in the north water polynya. *Philosophical Transactions of the Royal Society A*, 378, 20190353.

Parent, G. J., Pernet, F., Tremblay, R., Sévigny, J. M., & Ouellette, M. (2008). Modeling of membrane lipids in gills of adult hard clam *Mercenaria mercenaria* during declining temperature. *Aquatic Biology*, *3*, 101-109.

Parrish, C. C. (2009). Essential fatty acids in aquatic food webs. In M. T. Arts, M. T. Brett, & M. J. Kainz (Eds.), *Lipids in aquatic ecosystems* (pp. 309-326). Springer.

Paulsen, M. L., Nielsen, S. E. B., Müller, O., Møller, E. F., Stedmon, C. A., Juul-Pedersen, T., Markager, S., Sejr, M. K., Huertas, A. D., Larsen, A., & Middelboe, M. (2017). Carbon bioavailability in a high Arctic fjord influenced by glacial meltwater, NE Greenland. *Frontiers in Marine Science*, *4*, 176.

Pazos, A. J., Román, G., Acosta, C. P., Abad, M., & Sánchez, J. L. (1997). Seasonal changes in condition and biochemical composition of the scallop *Pecten maximus* L. from suspended culture in the ria de Arousa (Galicia, NW Spain) in relation to environmental conditions. *Journal of Experimental Marine Biology and Ecology*, *211*, 169-193.

Pernet, F., Tremblay, R., Comeau, L., & Guderley, H. (2007). Temperature adaptation in two bivalve species from different thermal habitat: Energetic and remodeling of membrane lipids. *The Journal of Experimental Biology*, *210*, 2999-3014.

Pernet, F., Tremblay, R., Gionet, C., & Landry, T. (2006). Lipid remodeling in wild and selectively bred hard clams at low temperatures in relation to genetic and physiological parameters. *The Journal of Experimental Biology*, *209*, 4663-4675.

Petersen, G. H. (1978). Life cycles and population dynamics of marine benthic bivalves from the Disko Bugt area in West Greenland. *Ophelia*, *17*, 95-120.

Petersen, G. H. (2001). *Studies on some Arctic and Baltic Astarte species (Bivalvia, Mollusca). Meddelelser om Grønland, bioscience 52.* Danish Polar Center.

Polozczanska, E. S., Burrows, M. T., Brown, C. J., Molinos, J. G., Halpern, B. S., Hoegh-Guldberg, O., Kappel, C. V., Moore, P. J., Richardson, A. J., Schoeman, D. S., & Sydeman, W. J. (2016). Responses of marine organisms to climate change across oceans. *Frontiers in Marine Science*, *3*, 62.

R Core Team. (2019). *R: A language and environment for statistical computing*. R Foundation for Statistical Computing.

Rantanen, M., Karpechko, A. Y., Lipponen, A., Nordling, K., Hyvärinen, O., Ruosteenoja, K., Vihma, T., & Laaksonen, A. (2022). The Arctic has warmed nearly four times faster than the globe since 1979. *Communications Earth & Environment*, *3*, 168.

Renaud, P. E., Sejr, M. K., Bluhm, B. A., Sirenko, B., & Ellingsen, I. H. (2015). The future of Arctic benthos: Expansion, invasion, and biodiversity. *Progress in Oceanography*, *139*, 244-257.

Renaud, P. E., Tessmann, M., Evenset, A., & Christensen, G. N. (2011). Benthic food-web structure of an Arctic fjord (Kongsfjorden, Svalbard). *Marine Biology Research*, *7*, 13-26.

Renaud, P. E., Wallhead, P., Kotta, J., Włodarska-Kowalczuk, M., Bellerby, R. G., Rätsep, M., Slagstad, D., & Kukliński, P. (2019). Arctic sensitivity? Suitable habitat for benthic taxa is surprisingly robust to climate change. *Frontiers in Marine Science*, *6*, 538.

Rysgaard, S., Nielsen, T. G., & Hansen, B. W. (1999). Seasonal variation in nutrients, pelagic primary production and grazing in a high-Arctic coastal marine ecosystem, young sound, Northeast Greenland. *Marine Ecology Progress Series*, *179*, 13-25.

Rysgaard, S., Vang, T., Stjernholm, M., Rasmussen, B., Windelin, A., & Kiilsholm, S. (2003). Physical conditions, carbon transport, and climate change impacts in a Northeast Greenland fjord. *Arctic, Antarctic, and Alpine Research*, *35*, 301-312.

Sejr, M. K., Bruhn, A., Dalsgaard, T., Juul-Pedersen, T., Stedmon, C. A., Blicher, M., Meire, L., Mankoff, K. D., & Thyrring, J. (2022). Glacial meltwater determines the balance between autotrophic and heterotrophic processes in a Greenland fjord. *Proceedings of the National Academy of Sciences of the United States of America*, *119*, e2207024119.

Sejr, M. K., & Christensen, P. B. (2007). Growth, production and carbon demand of macrofauna in young sound, with special emphasis on the bivalves *Hiatella arctica* and *Mya truncata*. In S. Rysgaard & R. N. Glud (Eds.), *Carbon cycling in Arctic marine ecosystems*:
Case study young sound (Vol. 58, pp. 122-135). Meddr Grønland. Bioscience.

Sejr, M. K., Jensen, K. T., & Rysgaard, S. (2000). Macrozoobenthic structure in a high-Arctic East Greenland fjord. *Polar Biology*, **23**, 792-801.

Sejr, M. K., Sand, M. K., Jensen, K. T., Petersen, J. K., Christensen, P. B., & Rysgaard, S. (2002). Growth and production of *Hiatella arctica* (Bivalvia) in a high-Arctic fjord (young sound, Northeast Greenland). *Marine Ecology Progress Series*, **244**, 163-169.

Sejr, M. K., Stedmon, C. A., Bendtsen, J., Abermann, J., Juul-Pedersen, T., Mortensen, J., & Rysgaard, S. (2017). Evidence of local and regional freshening of Northeast Greenland coastal waters. *Scientific Reports*, **7**, 13183.

Selin, N. I. (2010). The growth and life span of bivalve mollusks at the northeastern coast of Sakhalin Island. *Russian Journal of Marine Biology*, **36**, 258-269.

Sinensky, M. (1974). Homeoviscous adaptation – A homeostatic process that regulates the viscosity of membrane lipids in *Escherichia coli*. *Proceedings of the National Academy of Sciences of the United States of America*, **71**, 522-525.

Thyrring, J., Tremblay, R., & Sejr, M. K. (2017). Importance of ice algae and pelagic phytoplankton as food sources revealed by fatty acid trophic markers in a keystone species (*Mytilus trossulus*) from the high Arctic. *Marine Ecology Progress Series*, **572**, 155-164.

Tremblay, J.-É., Anderson, L. G., Matrai, P., Coupel, P., Bélanger, S., Michel, C., & Reigstad, M. (2015). Global and regional drivers of nutrient supply, primary production and CO₂ drawdown in the changing Arctic Ocean. *Progress in Oceanography*, **139**, 171-196.

Veillard, D., Tremblay, R., Sejr, M. K., Chauvaud, L., Cam, E., & Olivier, F. (2023). Recruitment dynamics of *Hiatella arctica* within a high Arctic

site (young Sound Fjord, NE Greenland). *Polar Biology*. https://doi.org/10.1007/s00300-023-03201-0

Wadham, J. L., Hawkings, J. R., Tarasov, L., Gregoire, L. J., Spencer, R. G. M., Gutjahr, M., Ridgwell, A., & Kohfeld, K. E. (2019). Ice sheets matter for the global carbon cycle. *Nature Communications*, **10**, 1-17.

Wassmann, P., Carmack, E. C., Bluhm, B. A., Duarte, C. M., Berge, J., Brown, K., Grebmeier, M., Holding, J., Kosobokova, K., Kwok, R., Matrai, P., Agusti, S., Babin, M., Bhatt, U., Eicken, H., Polyakov, I., Rysgaard, S., & Huntington, H. P. (2020). Towards a unifying pan-arctic perspective: A conceptual modelling toolkit. *Progress in Oceanography*, **189**, 102455.

## SUPPORTING INFORMATION

Additional supporting information can be found online in the Supporting Information section at the end of this article.

**How to cite this article:** Bridier, G., Olivier, F., Grall, J., Chauvaud, L., Sejr, M. K., & Tremblay, R. (2023). Seasonal lipid dynamics of four Arctic bivalves: Implications for their physiological capacities to cope with future changes in coastal ecosystems. *Ecology and Evolution*, **13**, e10691.
https://doi.org/10.1002/ece3.10691