==== Front PeerJ PeerJ PeerJ PeerJ 2167-8359 PeerJ Inc. San Diego, USA 15595 10.7717/peerj.15595 Ecology Marine Biology Zoology Opportunistic consumption of marine pelagic, terrestrial, and chemosynthetic organic matter by macrofauna on the Arctic shelf: a stable isotope approach http://orcid.org/0000-0002-9644-0179 Kokarev Valentin 1 http://orcid.org/0000-0001-5590-5673 Zalota Anna K. 1 Zuev Andrey 2 Tiunov Alexei 2 Kuznetsov Petr 3 Konovalova Olga 45 Rimskaya-Korsakova Nadezhda 3nadezdarkorsakova@gmail.com 1 Laboratory of Ecology of Coastal Bottom Communities, Shirshov Institute of Oceanology RAS, Moscow, Russia 2 Laboratory of Soil Zoology and General Entomology, A.N. Severtsov Institute of Ecology and Evolution RAS, Moscow, Russia 3 Department of Invertebrate Zoology, Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia 4 Centre of Marine Research, Lomonosov Moscow State University, Moscow, Russia 5 National Research Tomsk State University, Tomsk, Russia de Rezende Carlos 29 6 2023 2023 11 e1559517 1 2023 29 5 2023 © 2023 Kokarev et al. 2023 Kokarev et al. https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited. Macrofauna can contribute substantially to the organic matter cycling on the seafloor, yet the role of terrestrial and chemosynthetic organic matter in the diets of microphagous (deposit and suspension) feeders is poorly understood. In the present study, we used stable isotopes of carbon and nitrogen to test the hypothesis that the terrestrial organic matter supplied with river runoff and local chemosynthetic production at methane seeps might be important organic matter sources for macrofaunal consumers on the Laptev Sea shelf. We sampled locations from three habitats with the presumed differences in organic matter supply: “Delta” with terrestrial inputs from the Lena River, “Background” on the northern part of the shelf with pelagic production as the main organic matter source, and “Seep” in the areas with detected methane seepage, where chemosynthetic production might be available. Macrobenthic communities inhabiting each of the habitats were characterized by a distinct isotopic niche, mostly in terms of δ13C values, directly reflecting differences in the origin of organic matter supply, while δ15N values mostly reflected the feeding group (surface deposit/suspension feeders, subsurface deposit feeders, and carnivores). We conclude that both terrestrial and chemosynthetic organic matter sources might be substitutes for pelagic primary production in the benthic food webs on the largely oligotrophic Laptev Sea shelf. Furthermore, species-specific differences in the isotopic niches of species belonging to the same feeding group are discussed, as well as the isotopic niches of the symbiotrophic tubeworm Oligobrachia sp. and the rissoid gastropod Frigidoalvania sp., which are exclusively associated with methane seeps. Isotopic niche Laptev Sea Methane seeps Macrobenthic communities Deposit feeders Thyasiridae Bivalvia Siboglinidae Annelida Mollusca Russian Science Foundation20-74-10011 Lomonosov Moscow State University121032300121-0 Tomsk State University Development ProgramPriority-2030 This research was supported by the Russian Science Foundation (grant # 20-74-10011). This research was carried out as part of the Scientific Project of the State Order of the Government of Russian Federation to Lomonosov Moscow State University No. 121032300121-0. Contribution of Olga Konovalova was supported by the Tomsk State University Development Program (Priority-2030). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. ==== Body pmcIntroduction It is a widely accepted paradigm that environment and biodiversity jointly drive ecosystem functioning (van der Plas, 2019). Specifically, accounting for the functional diversity of benthic biota might improve understanding of various ecosystem functions on the seafloor including organic matter cycling (Snelgrove et al., 2018). Within marine sediments, organic matter mineralization in part depends on the processing of organic matter by macrofauna (Woulds et al., 2016; Middelburg, 2018). Different macrobenthic communities respond differently to inputs of organic matter depending on species and functional composition (Josefson, Forbes & Rosenberg, 2002; Norling et al., 2007; Karlson et al., 2010; Belley & Snelgrove, 2017). Moreover, organic matter differs in its bioavailability for macrobenthic consumers as well as its nutritional value depending on its origin and quality (Mayor et al., 2012; Hunter et al., 2013; Campanyà-Llovet, Snelgrove & Parrish, 2017). Therefore, identifying the organic matter sources for benthic food webs and functional differences within macrobenthic consumers are prerequisites for an improved understanding of benthic organic matter cycling. There are several major sources of organic matter on the Arctic shelf: phytoplankton and ice algae are the main sources of marine primary production (Sakshaug, 2004; Søreide et al., 2006), while river runoff and coastal erosion supply organic matter of terrestrial origin (Goñi et al., 2000; Knies & Martinez, 2009; Vonk et al., 2012; Xiao, Fahl & Stein, 2013). Macrobenthic communities mainly rely on the export flux of organic matter from the euphotic zone, both of sympagic (ice algae) and pelagic (phytoplankton) origin (Søreide et al., 2013; Mäkelä, Witte & Archambault, 2017). There is also evidence of terrestrial organic matter being incorporated into benthic food webs, although it is considered more refractory compared to organic matter of marine origin (Dunton, Schonberg & Cooper, 2012; Bell, Bluhm & Iken, 2016; Harris et al., 2018). The exact quality and quantity of organic matter reaching the seafloor might depend on many factors and processes, both biotic and abiotic, including depth, water mass properties and currents, overall primary production, terrestrial inputs, grazing by zooplankton, and microbial degradation of sinking organic matter (Iken, Bluhm & Gradinger, 2005; Roy et al., 2015; Stasko et al., 2018). Locally, benthic primary production is important for benthic food webs, including microphytobenthos and macroalgae in the coastal habitats (Woelfel et al., 2010, Renaud et al., 2015) and chemosynthetic production in the areas with seabed emission of reduced compounds, such as methane and hydrogen sulfide, also known as “methane seeps” and “hydrothermal vents” (Sweetman et al., 2013; Åström, Bluhm & Rasmussen, 2022). Consequently, regional and local ecosystem features might drive differences in organic matter availability and consumption by macrofauna. In this respect, the Laptev Sea appears to be among the least studied Arctic seas (Kędra et al., 2015). This part of the Arctic is considered a region of high biogeochemical activity due to considerable river runoff, coastal erosion, and methane seepage (Dittmar & Kattner, 2003, Vonk et al., 2012, Savvichev et al., 2018, Baranov et al., 2020). Macrofauna can be generally divided into microphages (e.g., deposit and suspension feeders) and macrophages (e.g., carnivores and herbivores), but some taxa, such as siboglinid tubeworms, are symbiotrophic and host chemosynthetic bacteria (Jumars, Dorgan & Lindsay, 2015; Åström et al., 2020). Soft-bottom macrobenthic communities are typically dominated by deposit and suspension feeding fauna (Kokarev et al., 2017; Kokarev et al., 2021a). The food items consumed by microphages include suspended and freshly deposited plankton, such as diatoms and flagellates, but also organic matter of various origins in the form of particulate matter, or detritus, and associated microorganisms (Levinton, 1972; North et al., 2014; Kędra et al., 2021). Since the visual analysis of food resources utilized by microphagous feeders can be problematic due to the small size of food items and high amount of consumed inorganic particles, indirect methods, such as stable isotopes of carbon and nitrogen, proved to be useful in identifying sources of organic matter consumed by macrofauna (Yokoyama & Ishihi, 2007; North et al., 2014). The ratio of 15N to 14N (δ15N) is reflective of trophic level due to enrichment at each trophic level by 2.5–5‰ while the ratio of 13C to 12C (δ13C) changes little in the food chains and is therefore used to trace carbon sources within food webs (Bearhop et al., 2004). Stable isotopes can be used both for comparing resource use between species (i.e., isotopic niche overlap) or community-wide aspects of trophic structure or resource use (Karlson, Gorokhova & Elmgren, 2015; Swanson et al., 2015; Jackson et al., 2011; Włodarska-Kowalczuk et al., 2019). However, in order to differentiate the consumption of different carbon sources, these carbon sources must differ in their isotopic signature. For instance, on the Laptev Sea shelf, macrobenthic consumers can be expected to exploit terrestrial, marine, and chemosynthetic organic matter (Kokarev et al., 2017; Vedenin et al., 2020). Carbon isotopic values of terrestrial organic matter are in the range of −28 to −26‰, while phytoplankton-derived organic matter is typically in the range of −22 to −20‰, although marine signatures can be more variable in the Arctic, especially when ice algae are considered with carbon isotopic signature in the range of −15 to −8‰ (Stein & Macdonald, 2004). However, there are seasonal and geographic variation in the baseline isotopic signatures of marine and terrestrial organic matter (Tamelander et al., 2009; Søreide et al., 2013; Bell, Bluhm & Iken, 2016; McClelland et al., 2016). In the areas of methane seepage, distinct microbial communities in the sediments are formed that include chemosynthetic microorganisms capable of utilizing either methane (methanotrophs) or hydrogen sulfide (thiotrophs) as an energy source (Åström et al., 2020). Carbon in the biomass of sulfur-oxidizing bacteria and methanotrophic microorganisms is usually characterized by very low δ13C values (<−30‰) (Levin & Michener, 2002; Decker & Olu, 2012; Åström, Bluhm & Rasmussen, 2022). In the present study, we hypothesized that macrofauna on the Laptev Sea shelf might utilize organic matter of three origins: marine pelagic, terrestrial, and chemosynthetic. To test this hypothesis, we have obtained data on the carbon and nitrogen stable isotope composition of macrofauna from three soft-bottom habitats with the presumed differences in organic matter sources. The “Background” habitat on the northern part of the shelf was chosen to represent the baseline for the consumption of pelagic primary production. In the “Delta” habitat reliance on the terrestrial organic matter supplied with Lena River runoff was expected in addition to the marine primary production. Finally, the “Seep” habitat in the areas with detected methane seepage where additional bacterial chemosynthetic production might be available for consumption, which is absent in the “Background” and “Delta”. Specifically, we aimed to answer the following questions: (1) How different are the isotopic niches of the studied communities, and (2) whether pronounced isotopic differences exist between species with similar feeding habits? Materials and Methods Study area and definition of habitats The Laptev Sea is a Siberian epicontinental sea with an average water depth of around 50 meters (Fig. 1). Prolonged ice coverage and considerable freshwater runoff, mainly supplied by the Lena River, largely shape ecological patterns on the eastern part of the Laptev Sea shelf (Schmid et al., 2006). The primary production and the fluxes of phytoplankton-derived organic matter to the seafloor generally decrease with increasing latitude (Stein & Fahl, 2004; Demidov, Sheberstov & Gagarin, 2020). The highest organic carbon content in the sediments is recorded east of the Lena Delta and is the consequence of the accumulation of terrestrial organic matter supplied with the river runoff (Stein & Fahl, 2004; Xiao, Fahl & Stein, 2013). There is a pronounced latitudinal gradient in the composition of macrobenthic communities associated with the riverine discharge and regime of sedimentation (Kokarev et al., 2017). Close to the Lena Delta, species-poor communities are mainly dominated by the bivalve Portlandia arctica, while on the northern part of the shelf, the number of species increases, and tube-dwelling annelid worms are abundant, including Myriochele heeri, Owenia polaris, and Maldane sarsi. In addition, on the northern part of the shelf, specific benthic habitats are associated with the methane seeps (Baranov et al., 2020; Vedenin et al., 2020). These habitats are characterized by the presence of carbonate crusts and white microbial mats on the sediment surface. Few species are found exclusively at the methane seeps, most notably the siboglinid tubeworm Oligobrachia sp. and the rissoid gastropod Frigidoalvania sp. 10.7717/peerj.15595/fig-1 Figure 1 Map of the study area. For this study, the sampling stations, nine in total, were divided into three types of soft-bottom benthic habitats based on the presumed differences in organic matter supply to the sediments (Fig. 1): the habitat in the vicinity of the Lena Delta with high sedimentation rates of riverine organic matter (“Delta”), the northern shelf habitat in the area with no detected methane seepage (“Background”) and methane seeps (“Seep”). Two stations from the “Delta” habitat, 6976 and 6977, were treated as one location as these stations are located close to each other (<500 m) and no small-scale heterogeneity in terms of organic matter sources was expected in this habitat based on the prior knowledge (Kokarev et al., 2017). Contrary, small-scale heterogeneity was expected at the “Seep” habitat (stations 5625, 6939, 6947, 6952, 6953, 6992) as there is evidence that chemosynthetic production is distributed very locally, e.g., patches of white microbial mats (Baranov et al., 2020; Vedenin et al., 2020). The “Background” station (6950) was located outside of the influence of the Lena River runoff but close to the “Seep” stations so that they would be similar in terms of bentho-pelagic coupling and fluxes of organic matter to the seafloor but differ in the sediment biogeochemical processes associated with microbial activities (Savvichev et al., 2018; Savvichev et al., 2023). Although the “Background” and “Delta” had only one location sampled, we consider the number of animals sampled from each of these habitats to be sufficient to represent the isotopic niche of respective communities (Table 1; Table S1). 10.7717/peerj.15595/table-1 Table 1 List of the macrobenthic taxa sampled. Numbers indicate the number of samples used for stable isotope analyses. Surface Deposit Feeders, SDF; Suspension Feeders, SF; Subsurface Deposit Feeders, SSDF (for details see section “Classification of the macrofauna based on feeding habits”). An asterisk (*) indicates taxa unassigned to a particular feeding habit Station-wise data and all the measured δ13C, δ15N, and C:N values for each sample are available in the Table S1. Taxonomic group Species Feeding habit Habitat Background Delta Seep Bivalvia Ennucula tenuis (Montagu, 1808) SDF/SF 1 1 Bivalvia Macoma calcarea (Gmelin, 1791) SDF/SF 2 2 Bivalvia Macoma moesta (Deshayes, 1855) SDF/SF 1 2 Bivalvia Macoma sp. SDF/SF 1 2 Bivalvia Macoma torelli (A. S. Jensen, 1905) SDF/SF 1 Bivalvia Nuculana pernula (O. F. Müller, 1779) SDF/SF 1 12 Bivalvia Periploma aleuticum (A. Krause, 1885) SDF/SF 1 1 Bivalvia Portlandia arctica (Gray, 1824) SDF/SF 2 5 6 Bivalvia Thracia myopsis Møller, 1842 SDF/SF 2 Bivalvia Yoldia hyperborea (A. Gould, 1841) SDF/SF 1 Bivalvia Yoldiella lenticula (Møller, 1842) SDF/SF 3 7 Bivalvia Yoldiella solidula Warén, 1989 SDF/SF 3 Bivalvia Parathyasira dunbari (Lubinsky, 1976)* n/a 1 1 Bivalvia Thyasira cf. gouldii* n/a 1 6 5 Bivalvia Montacuta spitzbergensis Knipowitsch, 1901* n/a 1 Gastropoda Frigidoalvania sp.* n/a 6 Annelida Nephtyidae gen. sp. Carnivore 1 3 14 Annelida Maldane sarsi Malmgren, 1865 SSDF 5 19 Annelida Pectinaria (Cistenides) hyperborea (Malmgren, 1866) SSDF 2 3 8 Annelida Sternapsis sp. SSDF 2 3 12 Annelida Myriochele heeri Malmgren, 1867 SDF/SF 1 28 Annelida Owenia polaris Koh, Bhaud & Jirkov, 2003 SDF/SF 1 9 Annelida Spiochaetopterus typicus M Sars, 1856 SDF/SF 2 Annelida Oligobrachia sp. Symbiotroph 27 Collection and processing of macrofauna For the present study, most of the material was collected during the 82nd cruise of the R/V “Akademik Mstislav Keldysh” in October 2020 (data on the sampling stations are available in the Table S2). The sediment samples were obtained at eight stations using a grab sampler “Okean” (0.25 m2). Sediments from the grab samples were sieved on the 0.5 mm mesh, and the residues were sorted under a stereomicroscope onboard. Specimens were cleared of debris (and polychaetes removed from their tubes), identified to the species level wherever possible, and frozen at −20 °C before further processing. In addition, the gastropods Frigidoalvania sp., collected in September 2017 during the 69th cruise of the R/V “Akademik Mstislav Keldysh” with a Sigsbee trawl and kept frozen at −20 °C, were used in this study (station 5625, for more details on 69th cruise see Vedenin et al., 2020 and Baranov et al., 2020). These gastropods were not found during the 82nd cruise, most probably due to a very patchy distribution based on grab samples (Vedenin et al., 2020). Overall, 24 taxa were sampled in this study (Table 1; Figs. 2 and 3). In the lab, the shells of mollusks were removed with a razor blade and tweezers, and for larger specimens the foot muscle was dissected for stable isotope analyses. For the smaller bivalves, annelids, and Frigidoalvania sp. all the soft tissues were used. For Thyasira cf. gouldii, Parathyasira dunbari, Yoldiella lenticula, Yoldiella solidula, and one sample of Ennucula tenuis, several individuals of the same species were pooled in one sample due to the very low amount of soft tissues in a single specimen. Subsequently, samples were dried at 60 °C and homogenized with mortar and pestle. Field samplings and experiments were approved by the Ministry of Education and Science of the Russian Federation (the field study approval number DN-09-54/52). 10.7717/peerj.15595/fig-2 Figure 2 Macrophotographs of the mollusk species studied. (A) Macoma calcarea, (B) Ennucula tenuis, (C) Yoldia hyperborea, (D) Thyasira cf. gouldii, (E) Periploma aleuticum, (F) Portlandia arctica, (G) Yoldiella lenticula, (H) Nuculana pernula, (I) Montacuta spitzbergensis, (G) Frigidoalvania sp., (K) Thracia myopsis, (L) Yoldiella solidula, (M) Parathyasira dunbari. All scale bars equal 1 mm. Photo credit: Valentin Kokarev. 10.7717/peerj.15595/fig-3 Figure 3 Macrophotographs of the studied annelids. All alive except (A) and (G). (A) Owenia polaris, (B) Sternapsis sp., (C) Aglaophamus malmgreni (Nephtyidae), (D) Spiochaetopterus typicus, (E) Pectinaria hyperborea, (F) Oligobrachia sp. (with metatrochophore larvae), (G) Maldane sarsi, (H) Myriochele heeri. All scale bars equal 1 mm. Photo credit: (A, G) Valentin Kokarev, (B, D–F, H) Petr Kuznetsov, (C) Tina Molodtsova. Collection and processing of sediments In addition to animal samples, at each station during the 82nd cruise (except for 6976) surface sediments were sampled from the grab prior to sieving with a spoon and immediately frozen at −20 °C. Subsequently, sediment samples were dried at 60 °C. For the stable isotope analyses, dried samples were divided into two replicates. One replicate was acidified with 2M HCl to remove inorganic carbon and was used to determine the δ13C value. Since acidification alters the δ15N value, non-acidified samples were used to determine δ15N (Silberberger et al., 2021). Stable isotope analyses Dried samples (200–500 μg for animal tissues and 3 to 10 mg for sediment samples) were wrapped into tin capsules and analyzed for stable carbon and nitrogen isotope composition using Flash 1112 Elemental Analyzer (Thermo Fisher Scientific, Waltham, MA, USA) and a Thermo Delta V Plus isotope ratio mass spectrometer connected via a Conflo IV peripheral at the A.N. Severtsov Institute of Ecology and Evolution RAS. Results are reported in permille (‰) using the delta (δ) notation and were normalized to VPDB and air N2 using USGS40 and USGS41 reference materials (US Geological Survey, Reston, VA, USA). The drift was corrected using an internal laboratory standard (casein). The standard deviation of δ13C and δ15N values in the laboratory standard (n = 8) was <0.2‰. Although C:N ratios were variable, particularly between foot vs. whole clam samples, which is clearly a consequence of differences in lipid content (Weems et al., 2012), the data were not corrected for lipids using mathematical corrections as available equations might not be suitable for benthic invertebrates such as bivalves and polychaetes (Silberberger et al., 2021). All the measured δ13C, δ15N, and C:N values for each animal sample are available in the Table S1. The values in the text are presented as mean ± SD. Classification of the macrofauna based on feeding habits Prior to the data analyses, species were assigned to trophic groups based on Jumars, Dorgan & Lindsay (2015) and Kokarev et al. (2017) (Table 1). Surface deposit feeders and suspension feeders were combined for two reasons: (1) some species can switch from deposit feeding to suspension feeding, e.g., Owenia, Spiochaetopterus, Macoma (Jumars, Dorgan & Lindsay, 2015; Reid & Reid, 1969) and (2) both groups utilize similar food sources: fresh phytoplankton material, food resources deposited in the upper layer of the sediments, which are available to suspension feeders through resuspension (Christensen & Kanneworff, 1985; Sokołowski et al., 2014). Moreover, only two species could be classified as strict suspension feeders: bivalves Thracia myopsis and Periploma aleuticum (Sartori & Domaneschi, 2005; Morton, 1981). Subsurface deposit feeders are represented by annelids that feed in the deeper layers of the sediment and rely on more microbially degraded organic matter, which is reflected in higher δ15N values (North et al., 2014; Jumars, Dorgan & Lindsay, 2015). The symbiotrophic tubeworm Oligobrachia sp. and carnivores nephtyid polychaetes were the only taxa in their respective feeding group. In the section “Seep-associated fauna” we used additional data on Oligobrachia and δ13C values of methane from North Atlantic and Arctic seeps for comparison on a larger scale (Gebruk et al., 2003; Decker & Olu, 2012; Lee et al., 2019a; Åström, Bluhm & Rasmussen, 2022; Kravchishina et al., 2021). Several taxa were not assigned to a particular feeding habit: the gastropod Frigidoalvania sp., the bivalve Montacuta spitzbergensis, and thyasirid bivalves. Frigidoalvania sp. is currently known only from the seeps in the Laptev Sea, and little is known about its feeding habits (Nekhaev & Krol, 2020). M. spitzbergensis belongs to the family Montacutidae, which comprises commensal bivalves, but its host is unknown with all the specimens (including this study) found freely in sediments (Kamenev, 2008). Feeding strategies in the thyasirid bivalves are diverse, particularly in the cryptic Thyasira cf. gouldii, which might derive its nutrition from sulfide-oxidizing bacteria in the gills, subsurface pedal feeding, and suspension feeding (Zanzerl & Dufour, 2017; Zanzerl et al., 2019). Data analysis Stable isotope data were visualized as biplots using R version 4.2.2 (R Core Team, 2022) and package ggplot2 version 3.4.1 (Wickham, 2016). Since the data partially were not normally distributed in bivariate isotopic space (tested with Shapiro-Wilk multivariate normality test in the package RVAideMemoire version 0.9-81-2; Hervé, 2022), the data analysis was performed with non-parametric approaches that do not assume normal distribution or homogeneity of variance. Differences between pre-defined groups (“Habitat”, “Feeding habit” and “Station”) were tested with a two-way non-parametric Scheirer Ray Hare test in the package rcompanion version 2.4.21 (Mangiafico, 2023) or one-way Kruskall-Wallis test and subsequent pairwise Dunn’s test with Holm’s procedure to adjust p-value for multiple comparisons in the package dunn.test version 1.3.5. (Dinno, 2017). The multivariate differences on the species level were investigated using non-parametric rank-based Wilks’ Lambda type statistic, for which p-value was calculated using F approximation, and the empirical nonparametric relative treatment effects for each variable in the package npmv version 2.4.0 (Ellis et al., 2017). The relative effects quantify the probability that a randomly chosen individual from one species exhibits a higher value for the variable (δ13C or δ15N) than a randomly chosen individual across all species. Results Habitat-associated patterns The sediment stable isotope data showed some differences between “Background”, “Seep” and “Delta” habitats: lower δ13C values were observed for the “Seep” and the “Delta” as well as a considerably lower δ15N value for “Delta” sediments (Table 2). “Background” and “Seep” habitats had similar δ15N values. 10.7717/peerj.15595/table-2 Table 2 δ13C and δ15N measured for the sediment samples. An asterisk (*) indicates average values for two samples. Background Delta Seep 6950 6977 6939 6947* 6952 6953 6992 δ13C, ‰ −23.5 −25.8 −26.0 −24.2 −24.4 −24.4 −23.6 δ15N, ‰ 5.9 3.0 5.3 5.6 6.0 5.9 5.7 The patterns within macrobenthic consumers were more prominent. There was a considerable difference between the symbiotrophic annelid Oligobrachia sp. from “Seep” stations and all the other macrobenthic taxa due to non-overlapping δ13C values (−48.4 ± 3.1‰ for Oligobrachia sp. and −24.3 ± 2.7‰ for all the other taxa). The dispersion of the Oligobrachia sp. samples in the isotopic space was of comparable magnitude to all the other taxa combined (Fig. S1). Therefore, subsequent comparisons of different habitats were performed excluding Oligobrachia sp. The exclusion of Oligobrachia sp. revealed the significant differences between macrobenthic consumers from the three habitats (Fig. 4A). The “Delta” community was characterized by lower values of δ13C (−27.9 ± 1.6‰) compared to the “Background” community (−22.7 ± 1.2‰). The “Seep” community (−24.0 ± 2.6‰) with the widest isotopic niche overlapped both with the “Delta” and “Background” communities. The two-way Scheirer Ray Hare test (excluding taxa unassigned to feeding group and carnivores due to their low representation (n < 5) in the “Delta” and “Background”) revealed that the effect of the interaction of “Habitat” and “Feeding habit” on the δ13C is not significant (H = 5.2, df = 2, p = 0.075), as well as the main effect of “Feeding habit” (H = 3.6, df = 1, p = 0.056), while the main effect of “Habitat” was significant (H = 31.7, df = 2, p < 0.001). Further, differences in the range of δ13C values were revealed between different “Seep” stations (Fig. 5): while stations 6939 and 6992 were completely within the “Background” range, considerably lower δ13C values were recorded for the stations 6952, 6953, and 6947, with stations 6953 and 6947 being significantly different from “Background” (results of Kruskal-Wallis test and subsequent Dunn’s test are presented in the Table S3). 10.7717/peerj.15595/fig-4 Figure 4 Stable isotope biplot of δ13C and δ15N values for the macrobenthic taxa excluding Oligobrachia sp. Habitats are locations on the Laptev Sea shelf with the presumed differences in organic matter sources for macrofauna (marine vs. terrestrial vs. chemosynthetic). Group centroids (calculated as mean δ13C and δ15N values) are shown. (A) and (B) are the same plot but differently color-coded. Surface Deposit Feeders, SDF; Suspension Feeders, SF; Subsurface Deposit Feeders, SSDF. An asterisk (*) indicates taxa unassigned to a particular feeding habit. Thyasridae was represented by two species: Thyasira gouldi and Parathyasira dunbari. 10.7717/peerj.15595/fig-5 Figure 5 Boxplot of δ13C values for the macrobenthic taxa (excluding Oligobrachia sp) grouped by sampling station. The station 5625, where only Frigidoalvania sp. was sampled, is not shown. Blue dashed lines correspond to δ13C range (min-max) of the “Background” habitat. The red dashed line corresponds to the carbon isotopic signature (−35) of sulphur-oxidizing bacteria at high-latitude seeps (Åström, Bluhm & Rasmussen, 2022). The green dashed line corresponds to the isotopic signature of Lena River (−31.1) particulate organic matter (McClelland et al., 2016). Macrobenthic consumers with different feeding habits could be differentiated based on their δ15N values (Fig. 4B): surface deposit feeders and suspension feeders had the lowest and carnivores the highest δ15N values (7.1 ± 0.9‰ and 11.6 ± 0.6‰, respectively) with subsurface deposit feeders in between (10.3 ± 1.2‰). In addition, centroid of the “Delta” community was clearly skewed towards lower δ15N (7.1 ± 2.2‰) compared to “Background” (8.6 ± 2.3‰) and “Seep” (8.5 ± 2.0‰). The two-way Scheirer Ray Hare test (excluding taxa unassigned to feeding group and carnivores due to their low representation (n < 5) in the “Delta” and “Background”) revealed that the effect of interaction of “Habitat” and “Feeding habit” on the δ15N is not significant (H = 1.4, df = 2, p = 0.500), while the main effects of “Habitat” (H = 13.5, df = 2, p = 0.001) and “Feeding habit” (H = 93.7, df = 1, p < 0.001) are significant. Among taxa with unassigned feeding habits, thyasirid bivalves and the gastropod Frigidoalvania sp. clearly grouped with surface deposit and suspension feeders (Fig. 4B). The single sample of the bivalve Montacuta spitzbergensis grouped with a few outliers from the “Seep” community with the lowest δ13C values. In the following section “Background” and “Seep” samples for surface deposit/suspension feeders and subsurface deposit feeders are discussed together since samples from these two habitats clearly overlap (Figs. 4 and 5). Species isotopic niches Surface deposit/suspension feeders “Delta” surface deposit/suspension feeders were represented only by Portlandia arctica (δ13C = −29.0 ± 1.0‰; δ15N = 5.6 ± 0.2‰) and Macoma spp. (δ13C = −29.0 ± 0.4‰; δ15N = 6.0 ± 0.3 ‰). This group occupied a narrow and distinct isotopic niche that did not overlap with the isotopic niches of the same taxa from “Seep” and “Background” (P. arctica: δ13C = −22.5 ± 1.7‰; δ15N = 6.5 ± 0.4‰; Macoma spp.: δ13C = −23.5 ± 1.4‰; δ15N = 6.9 ± 0.4‰). Surface deposit/suspension feeders from “Background” and “Seep” habitats were represented by 12 taxa, and there was a general trend in lower δ13C values being associated with lower δ15N values (Fig. 6). Two “Seep” samples of the bivalve Thracia myopsis had the lowest δ13C values (≈−30‰). Oweniid annelids Myriochele heeri and Owenia polaris along with annelid Spiochaetopterus typicus and bivalve Periploma aleuticum had values of the δ13C higher than −24‰. Contrary, the majority of the other bivalve species (Yoldiella lenticula, Nuculana pernula, Portlandia arctica, Macoma spp.) showed a wide range of δ13C values with some of the “Seep” samples depleted in 13C. The bivalve Y. solidula had higher δ15N values compared to Y. lenticula, N. pernula, P. arctica, and Macoma spp. The lowest δ13C value (−25.3‰) from the “Background” was observed for the bivalve Ennucula tenuis, which appeared to be an outlier compared to the rest of the “Background” species, but a similar value for this species was recorded from the “Seep” habitat. However, the samples of E. tenuis had a high C:N ratio (>6), which might have resulted in decreased δ13C values due to high lipid content. 10.7717/peerj.15595/fig-6 Figure 6 Stable isotope biplot of δ13C and δ15N values for surface deposit and suspension feeders from the northern part of the Laptev Sea shelf. Group centroids (calculated as mean δ13C and δ15N values) are shown for groups with 7 or more samples.The following samples are not shown to improve readability: “Background” E. tenuis (δ13C = −25.3‰, δ15N = 6.1‰), “Seep” E. tenuis (δ13C = −25.2‰, δ15N = 5.8‰), “Background” P. aleuticum (δ13C = −22.6‰, δ15N = 6.8‰), “Seep” P. aleuticum (δ13C = −21.9‰, δ15N = 6.5‰), “Seep” S. typicus (δ13C = −22.8‰, δ15N = 7.7‰ and δ13C = −21.3‰, δ15N = 8.1‰), “Seep” Y. hyperborea (δ13C = −24.9‰, δ15N = 6.4‰). The multivariate differences among surface deposit and suspension feeding species from “Background” and “Seep” habitats with sample size ≥7 were significant (Wilks Lambda = 10.3, df1 = 10, df2 = 140, p < 0.001). The relative effects revealed that these differences were mainly due to the difference between surface deposit feeding bivalves and oweniid polychaetes, both in terms of δ13C and δ15N (Table 3). Moreover, O. polaris had the most distinct niche. 10.7717/peerj.15595/table-3 Table 3 Nonparametric relative effects for surface deposit and suspension feeders. The relative effects quantify probability that a randomly chosen individual from one species exhibits a higher value for the variable (δ13C or δ15N) than a randomly chosen individual across all species. Species δ13C δ15N Macoma spp. 0.35714 0.41466 Myriochele heeri 0.60099 0.66906 Nuculana pernula 0.28921 0.27123 Owenia polaris 0.79805 0.80195 Portlandia arctica 0.59578 0.25731 Yoldiella lenticula 0.20649 0.25909 Subsurface deposit feeders The polychaetes Sternapsis sp. and Pectinaria hyperborea from the “Delta” habitat had very different δ15N values (6.8–7.3‰ for Sternapsis sp. and 9.7–10.7‰ for P. hyperborea), and therefore their isotopic niches did not overlap. Within “Background” and “Seep” habitats three sampled species had three overlapping but significantly different isotopic niches (Fig. 7; Wilks Lambda = 10.2, df1 = 4, df2 = 88, p < 0.001). The relative effects revealed that M.sarsi had generally higher δ15N, while Sternapsis sp. had the highest δ13C value (Table 4). 10.7717/peerj.15595/fig-7 Figure 7 Stable isotope biplot of δ13C and δ15N values for subsurface deposit feeders from the northern part of the Laptev Sea shelf. Group centroids (calculated as mean δ13C and δ15N values) are shown. 10.7717/peerj.15595/table-4 Table 4 Nonparametric relative effects for subsurface deposit. The relative effects quantify probability that a randomly chosen individual from one species exhibits a higher value for the variable (δ13C or δ15N) than a randomly chosen individual across all species. Species δ13C δ15N Maldane sarsi 0.41667 0.62153 Pectinaria hyperborea 0.37708 0.36771 Sternapsis sp. 0.73065 0.38616 Carnivores Represented by the single taxon (nephtyid polychaetes), this group showed the highest δ15N (11.6 ± 0.6‰) values but had a similar δ13C (24.0 ± 2.7‰) range to the surface deposit/suspension feeders and subsurface deposit feeders (Fig. 4B). Thyasiridae Thyasirid bivalves grouped with suspension/surface deposit feeders (Fig. 4B). However, two samples of Parathyasira dunbari could be distinguished from Thyasira cf. gouldii on the northern part of the shelf by lower δ15N values (4.9‰ and 5.9‰ vs. 7.5 ± 0.5‰). Moreover, P. dunbari had the lowest recorded δ15N value for “Background” and “Seep” habitats among all the sampled taxa, excluding Oligobrachia sp. Seep-associated fauna Two species were associated exclusively with the “Seep” sites: the rissoid gastropod Frigidoalvania sp. and the tubeworm Oligobrachia sp. Frigidoalvania sp. had relatively low values of δ13C (−27.2 ± 0.5‰) and a distinct niche among surface deposit/suspension feeders from the seep sites (Fig. 4B). The samples of Oligobrachia sp. had a very wide dispersion in the isotopic space and had the lowest values both for δ13C and δ15N (down to −53.9‰ and to −4.0‰, respectively; Fig. S1). The analysis of Oligobrachia isotopic signatures from different seep locations in the North Atlantic and the Arctic revealed no clear relationship between δ13C of methane and δ13C of Oligobrachia (Fig. S2). Discussion Habitat-associated patterns Macrobenthic communities from the three studied habitats differed in terms of the δ13C range, reflecting variability in organic matter sources for macrobenthic communities. However, there were no differences in the δ13C range between different feeding groups, including primary consumers (deposit and suspension feeders) and secondary consumers (carnivores). The differences in the δ13C range between habitats were much more pronounced for macrofauna than for sediments, thus, the isotopic signatures of the sediments cannot be used as a food source baseline. “Background” macrofaunal samples had δ13C values larger than −24‰ (except for two samples with values ≈−25‰), indicating the predominant reliance of the community on marine organic matter (Vonk et al., 2012; Bell, Bluhm & Iken, 2016). However, the role of marine organic matter decreased for “Delta” and partially for “Seep” habitats. Delta On the eastern Laptev Sea and East Siberian Sea shelf macrobenthic communities are influenced by riverine input and regime of sedimentation, which results in low species and functional diversity close to the deltas due to environmental filtering (Kokarev et al., 2017, 2021a). The species found close to the Lena Delta, are wide-spread shelf species, which can successfully establish a population in waters with warmer temperatures, fluctuating salinity, and enhanced sedimentation (e.g., Portlandia arctica, Macoma spp., Thyasira cf. gouldii, Sternapsis sp., Pectinaria hyperborea, nephtyid annelid Aglaophamus malmgreni). In the present study, we showed that these species might consume organic matter of different origins depending on availability as inferred from their variable δ13C signatures. Macrofauna from the “Delta” habitat, especially surface deposit feeding bivalves, was considerably depleted in 13C and slightly depleted in 15N relative to “Background” (down to −30.5‰ and 5.3‰, respectively) directly reflecting the isotopic signature of Lena River particulate organic matter (δ13C = −31.1 ± 0.8‰; McClelland et al., 2016) and terrestrial end-member values used for the Arctic shelf (δ13C = −28.8 ± 3.2‰, δ15N = 0.8 ± 1.0‰; Bell, Bluhm & Iken, 2016). This fact suggests that terrestrial organic matter is a primary food source for the “Delta” community at least during the season the samples were taken. Such low values for δ13C and δ15N were not previously recorded for the Arctic macrofauna in the areas influenced by river runoff (Dunton, Schonberg & Cooper, 2012; Bell, Bluhm & Iken, 2016; Harris et al., 2018; Stasko et al., 2018), which either reflects the higher contribution of terrestrial organic matter to the benthic food web close to the Lena Delta, or seasonal/geographical variations in the isotopic signatures of fluvial particulate organic matter (McClelland et al., 2016). Seep Little is known about the nutrition of macrobenthos, particularly deposit feeders, from the arctic cold seeps (Åström et al., 2020). Our data revealed a wider δ13C range for the “Seep” community relative to “Background” indicating that carbon of chemosynthetic origin is actively consumed by the non-symbiotrophic deposit feeding macrofauna in addition to photosynthetic carbon. Individuals of many species were more depleted in 13C at seep sites compared to “Background”, which is consistent with patterns observed at some other methane seeps, usually located below shelf depths (Levin & Michener, 2002; Levin, 2005; Thurber et al., 2010; Sellanes et al., 2011), and indicates partial reliance on chemosynthetic production which is 13C-depleted (δ13C ≈−35‰ for sulfur-oxidising bacteria from white microbial mats at high-latitude seeps, Åström et al., 2020). It should be noted, however, that for many species from “Seep” sites along with more depleted isotopic signatures, isotopic signatures within the “Background” range were observed. Moreover, samples from stations 6939 and 6992 were completely within “Background” δ13C range. Both of these facts suggest the patchy distribution of chemosynthetic organic matter and its opportunistic consumption. Indeed, video observations from the seeps in the Laptev Sea documented patchy distributions of white microbial mats at scales <1 m (Baranov et al., 2020). However, these results differ from the shelf of the Barents Sea, where differences between isotopic niches of seep and non-seep communities were less pronounced with no apparent consumption of chemosynthetic organic matter at shallow (<200 m) seeps (Åström, Bluhm & Rasmussen, 2022). This supports the previous hypothesis that in the less productive Laptev Sea chemosynthetic organic matter is more readily consumed by macrofauna supporting higher abundances at the seep sites (Vedenin et al., 2020), which likely makes the Laptev Sea seeps more of an exception among the shallow-water seeps (Dando, 2010). Species isotopic niches Surface deposit and suspension feeders Surface deposit feeders from the “Delta” habitat, represented by bivalves Macoma spp. and Portlandia arctica, had very similar isotopic signatures. On the northern part of the shelf at “Seep” and “Background” habitats, similarly, deposit feeding bivalves Macoma spp., P. arctica, Yoldiella lenticula, and Nuculana pernula had high niche overlap and appeared to have a contribution from chemosynthetic carbon to their diet at “Seep” sites. This finding corroborates previous results that different species of surface deposit feeding bivalves generally consume similar food sources from the sediment pool (North et al., 2014; Oxtoby et al., 2016). Interestingly, Yoldiella solidula had slightly higher δ15N values compared to the above-mentioned bivalves, possibly indicating a different resource use pattern for this species. Oweniids Myriochele heeri and Owenia polaris had niches distinct from each other and limited overlap with most abundant deposit feeding bivalves. Oweniids might selectively consume algal material recently deposited on the seafloor or from the near-bottom suspension (Jumars, Dorgan & Lindsay, 2015). The high densities of oweniid polychaetes are also sometimes linked to the input of fresh pelagic organic matter to the seafloor (Fiege, Kröncke & Barnich, 2000; Kokarev et al., 2021b). In our study, M. heeri and O. polaris had δ13C >−24‰, indicative of phytoplankton consumption (Vonk et al., 2012; Bell, Bluhm & Iken, 2016). Probably, dense populations of these species on the northern Laptev Sea shelf (Kokarev et al., 2017) rely less on bacterially reworked sediment pool of organic matter compared to deposit feeding bivalves and selectively feed on fresh phytodetritus, avoiding consumption of bacterially derived chemosynthetic organic matter at seep sites. The differences in the isotopic niches between M. heeri and O. polaris are most probably related to morpho-functional differences between the two genera: Owenia uses a tentacular crown to feed, while Myriochele lacks any anterior appendages (Fig. 3; Jumars, Dorgan & Lindsay, 2015). Consumption of chemosynthetic carbon was not evident also for the polychaete Spiochaetopterus typicus and the bivalves Periploma aleuticum, Yoldiella solidula, and Ennucula tenuis, but only a few samples of these species were analyzed. However, two samples of the bivalve Thracia myopsis and a single sample of the bivalve Montacuta spitzbergensis were depleted in 13C characteristic of chemosynthetic carbon consumption. T. myopsis is a strict suspension feeder and has common for the genus mucous lining of the siphon passage in the sediment that restricts the deposited sediment from entering the inhalant siphon aperture (Sartori & Domaneschi, 2005). Thus, T. myopsis feeds on suspended organic matter, which was shown to be depleted in 13C at the seep sites (δ13C = −34.8 to –35.9‰; Savvichev et al., 2018). Subsurface deposit feeders Subsurface deposit feeders, e.g., the polychaete Pectinaria hyperborea, consume mainly bacterially reworked detritus and are characterized by higher δ15N relative to surface deposit feeding bivalves (North et al., 2014). In the present study, we observed a similar pattern. However, three of the studied species had distinct isotopic niches. In the “Delta” community, Sternapsis sp. and P. hyperborea consumed different resources, while on the northern part of the shelf, at “Background” and “Seep” sites, their isotopic niches partially overlapped with each other and Maldane sarsi. It is not clear how these species partition resources, but possible mechanisms include feeding at different sediment levels or selection for different particle sizes (Hughes, 1979; Whitlatch, 1980). For all three species consumption of chemosynthetic organic matter was evident at the seep sites. Thyasiridae Thyasirid bivalves grouped with surface deposit feeders, although they form pedal tracts that penetrate deep into the sediment (Zanzerl & Dufour, 2017), while specimens identified as Thyasira cf. gouldii can be also symbiotrophic (Batstone et al., 2014). Isotopic signatures of T. cf. gouldii in the present study varied widely among habitats suggesting direct consumption of deposited organic matter in the surface sediment layer rather than symbiont-based nutrition. Interestingly, samples of Parathyasira dunbari with similar δ13C values to T. cf. gouldii had lower δ15N values. Such a pattern was also observed in Bonne Bay, Newfoundland, where symbiotic T. cf. gouldii, asymbiotic T. cf. gouldii, and Parathyasira sp. could be distinguished by δ15N but not δ13C values, and in conjunction with fatty acid analysis it was hypothesized that Parathyasira sp. might partially rely on free-living sulfur-oxidizing bacteria in the sediment (Zanzerl et al., 2019). We observed numerous bacteria on the inner fold of the mantle margin for Parathyasira dunbari during SEM investigation of this species, but not for T. cf. gouldii, supporting this hypothesis (Kokarev et al., in prep.). Seep-associated fauna Rissoid gastropods are characteristic of seep and vent fauna in the Arctic and North Atlantic, where they are associated with bacterial mats of sulfur-oxidizing bacteria, particularly of the genus Sulfurovum (Gebruk et al., 2003, Decker & Olu, 2012; Sweetman et al., 2013; Sen et al., 2019). Sulfurovum is abundant also in seep bacterial communities in the Laptev Sea (Savvichev et al., 2018). Although in our study rissoid gastropod Frigidoalvania sp. showed quite low δ13C values (−27.2% ± 0.5‰) relative to the “Background” surface deposit feeders, these values were quite high relative to the ones measured for Alvania sp. from the seeps in the North Atlantic (Decker & Olu, 2012). This indicates either a feeding habit on a mixture of different bacteria (Sweetman et al., 2013) or probably a mixed diet with significant contributions of photosynthetic organic matter. The siboglinid worm Oligobrachia sp. showed very variable isotopic signatures in the Laptev Sea and across a wider geographical range, which is not directly related to the isotopic signature of methane (Fig. S2). Although there are several cryptic Oligobrachia species in the Arctic and North Atlantic, similar symbiotic sulfur-oxidizing bacteria were detected in all the species (Lösekann et al., 2008; Sen et al., 2018; Lee et al., 2019b; Karaseva et al., 2021). A variable pattern in the isotopic signatures might be a result of autotrophic fixation of inorganic carbon from a mixture of sources: the inorganic carbon from the overlaying seawater and a more 13C-depleted pore-water inorganic carbon resulting from anaerobic oxidation of methane in the sediment (Lösekann et al., 2008; Lee et al., 2019a). This contrasts with the pattern observed for the siboglinid Siboglinum poseidoni hosting methanotrophic bacteria, which carbon isotopic signature directly reflects the isotopic signature of methane (Schmaljohann et al., 1990). High variation in δ15N values might also indicate a wide spectrum of nitrogen sources available for Oligobrachia, possibly including local assimilation of inorganic nitrogen (Levin, 2005) or amino acids from dissolved organic matter (Southward & Southward, 1981). Therefore, Oligobrachia might benefit from a wide spectrum of carbon and nitrogen sources available for the worm simultaneously, resulting in high variability of isotopic signatures even on small scales. More studies are needed on the biology of this species and its symbiotic bacteria in order to assess its role in carbon and nutrient cycling at methane seeps. Conclusion Our study revealed that the terrestrial and chemosynthetic organic matter might be more substantial food sources for macrofauna than previously demonstrated for the Arctic shelf. Distribution of many species, e.g., surface deposit feeding bivalves and subsurface deposit feeding polychaetes found in all three habitats, is not constrained by the origin of organic matter as demonstrated by the high variability of their isotopic signatures and suggests opportunistic consumption of available resources. However, differences in resource use among some species with the same feeding habits were observed, particularly between oweniid polychaetes and bivalves, suggesting that species-specific traits might be important for ecosystem functioning (Norling et al., 2007; Godbold, Rosenberg & Solan, 2009). These findings might have important implications for assessing macrobenthic carbon demand and the carbon budget of the Arctic shelf ecosystem (Schmid et al., 2006; Zaborska et al., 2018). It must be noted, however, that both consumptions of terrestrial and chemosynthetic organic matter results in lower δ13C values of macrofaunal consumers, making the interpretation of data contextual. Therefore, a combination of a stable isotope approach with supplemental methods, such as fatty acid analysis, might be a more beneficial yet more time-consuming and less cost-effective strategy to further assess the consumption of terrestrial/chemosynthetic organic matter by macrofauna. Supplemental Information 10.7717/peerj.15595/supp-1 Supplemental Information 1 Stable isotope biplot of δ13C and δ15N values for all macrobenthic taxa sampled showing the difference between symbiotrophic Oligobrachia sp. and heterotrophic taxa. Group centroids (calculated as mean δ13C and δ15N values) are shown. Click here for additional data file. 10.7717/peerj.15595/supp-2 Supplemental Information 2 Relationship between δ13C of methane and δ13C of Oligobrachia (mean ± sd) at several seep locations in the Arctic and North Atlantic. Based on the data from the present study and Gebruk et al., 2003; Decker & Olu, 2012; Lee et al., 2019a; Åström, Bluhm & Rasmussen, 2022; Kravchishina et al., 2021. HMMV, Haakon Mosby Mud Volcano. Click here for additional data file. 10.7717/peerj.15595/supp-3 Supplemental Information 3 δ13C, δ15N and C/N values for macrofaunal samples. Click here for additional data file. 10.7717/peerj.15595/supp-4 Supplemental Information 4 Data on the sampling stations. Click here for additional data file. 10.7717/peerj.15595/supp-5 Supplemental Information 5 Results of Kruskal-Wallis test and Dunn’s test for differences in the range of δ13C values at different sampling stations (excluding Oligobrachia sp.). Significant pair-wise comparisons are in bold. Click here for additional data file. We thank the crew and participants of the 82nd cruise of R/V “Akademik Mstislav Keldysh” for successful sampling in the Laptev Sea. Additional Information and Declarations Competing Interests Author Contributions Field Study Permissions Data Availability The authors declare that they have no competing interests. Valentin Kokarev conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft. Anna K. Zalota performed the experiments, authored or reviewed drafts of the article, and approved the final draft. Andrey Zuev performed the experiments, authored or reviewed drafts of the article, and approved the final draft. Alexei Tiunov performed the experiments, authored or reviewed drafts of the article, and approved the final draft. Petr Kuznetsov performed the experiments, authored or reviewed drafts of the article, and approved the final draft. Olga Konovalova performed the experiments, authored or reviewed drafts of the article, and approved the final draft. Nadezhda Rimskaya-Korsakova conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the article, and approved the final draft. The following information was supplied relating to field study approvals (i.e., approving body and any reference numbers): Field samplings and experiments were approved by the Ministry of Education and Science of the Russian Federation. The following information was supplied regarding data availability: The raw data are available in the Supplemental Files. ==== Refs References Åström, Bluhm & Rasmussen (2022) Åström EKL Bluhm BA Rasmussen TL Chemosynthetic and photosynthetic trophic support from cold seeps in Arctic benthic communities Frontiers in Marine Science 2022 9 August 4411 10.3389/fmars.2022.910558 Åström et al. (2020) Åström EKL Sen A Carroll ML Carroll JL Cold seeps in a warming arctic: insights for benthic ecology Frontiers in Marine Science 2020 7 587 10.3389/fmars.2020.00244 Baranov et al. (2020) Baranov B Galkin S Vedenin A Dozorova K Gebruk A Flint M Methane seeps on the outer shelf of the Laptev Sea: characteristic features, structural control, and benthic fauna Geo-Marine Letters 2020 40 4 541 557 10.1007/s00367-020-00655-7 Batstone et al. (2014) Batstone RT Laurich JR Salvo F Dufour SC Divergent chemosymbiosis-related characters in thyasira cf. gouldi (Bivalvia: Thyasiridae) PLOS ONE 2014 9 3 e92856 10.1371/journal.pone.0092856 24658402 Bearhop et al. (2004) Bearhop S Adams CE Waldron S Fuller RA Macleod H Determining trophic niche width: a novel approach using stable isotope analysis Journal of Animal Ecology 2004 73 5 1007 1012 10.1111/j.0021-8790.2004.00861.x Bell, Bluhm & Iken (2016) Bell LE Bluhm BA Iken K Influence of terrestrial organic matter in marine food webs of the Beaufort Sea shelf and slope Marine Ecology Progress Series 2016 550 1 24 10.3354/meps11725 Belley & Snelgrove (2017) Belley R Snelgrove PVR The role of infaunal functional and species diversity in short-term response of contrasting benthic communities to an experimental food pulse Journal of Experimental Marine Biology and Ecology 2017 491 38 50 10.1016/j.jembe.2017.03.005 Campanyà-Llovet, Snelgrove & Parrish (2017) Campanyà-Llovet N Snelgrove PVR Parrish CC Rethinking the importance of food quality in marine benthic food webs Progress in Oceanography 2017 156 July 240 251 10.1016/j.pocean.2017.07.006 Christensen & Kanneworff (1985) Christensen H Kanneworff E Sedimenting phytoplankton as major food source for suspension and deposit feeders in the Øresund Ophelia 1985 24 3 223 244 10.1080/00785326.1985.10429730 Dando (2010) Dando PR Kiel S Biological communities at marine shallow-water vent and seep sites paul R The Vent and Seep Biota Aspects from Microbes to Ecosystems 2010 Berlin Springer 333 378 Decker & Olu (2012) Decker C Olu K Habitat heterogeneity influences cold-seep macrofaunal communities within and among seeps along the Norwegian margin-Part 2: contribution of chemosynthesis and nutritional patterns Marine Ecology 2012 33 2 231 245 10.1111/j.1439-0485.2011.00486.x Demidov, Sheberstov & Gagarin (2020) Demidov AB Sheberstov SV Gagarin VI Seasonal variability and annual primary production of phytoplankton in the laptev sea assessed by MODIS-aqua data Izvestiya, Atmospheric and Oceanic Physics 2020 56 9 950 962 10.1134/S000143382009008X Dinno (2017) Dinno A dunn.test: Dunn’s test of multiple comparisons using rank sums. R package version 1.3.5 2017 https://cran.r-project.org/web/packages/dunn.test/ Dittmar & Kattner (2003) Dittmar T Kattner G The biogeochemistry of the river and shelf ecosystem of the Arctic Ocean: a review Marine Chemistry 2003 83 3–4 103 120 10.1016/S0304-4203(03)00105-1 Dunton, Schonberg & Cooper (2012) Dunton KH Schonberg SV Cooper LW Food web structure of the alaskan nearshore shelf and estuarine lagoons of the beaufort sea Estuaries and Coasts 2012 35 2 416 435 10.1007/s12237-012-9475-1 Ellis et al. (2017) Ellis AR Burchett WW Harrar SW Bathke AC Nonparametric inference for multivariate data: the R package npmv Journal of Statistical Software 2017 76 1 1 18 10.18637/jss.v076.i04 36568334 Fiege, Kröncke & Barnich (2000) Fiege D Kröncke I Barnich R High abundance of Myriochele fragilis Nilsen & Holthe, 1985 (Polychaeta: Oweniidae) in the deep sea of the Eastern Mediterranean Hydrobiologia 2000 426 1 97 103 10.1023/A:1003947318925 Gebruk et al. (2003) Gebruk AV Krylova EM Lein AY Georgy M Anderson E Pimenov NV Georgy A Crane K Methane seep community of the Håkon Mosby mud volcano (the Norwegian Sea): composition and trophic aspects Sarsia 2003 88 6 37 41 10.1080/00364820310003190 Godbold, Rosenberg & Solan (2009) Godbold JA Rosenberg R Solan M Species-specific traits rather than resource partitioning mediate diversity effects on resource use PLOS ONE 2009 4 10 e7423 10.1371/journal.pone.0007423 19829703 Goñi et al. (2000) Goñi MA Yunker MB MacDonald RW Eglinton TI Distribution and sources of organic biomarkers in arctic sediments from the Mackenzie River and Beaufort Shelf Marine Chemistry 2000 71 1–2 23 51 10.1016/S0304-4203(00)00037-2 Harris et al. (2018) Harris CM McTigue ND McClelland JW Dunton KH Do high Arctic coastal food webs rely on a terrestrial carbon subsidy? Food Webs 2018 15 e00081 10.1016/j.fooweb.2018.e00081 Hervé (2022) Hervé M RVAideMemoire: testing and plotting procedures for biostatistics. R package version 0.9-81-2 2022 https://cran.r-project.org/package=RVAideMemoire Hughes (1979) Hughes TG Mode of life and feeding in maldanid polychaetes from St. Margaret’s Bay, Nova Scotia Journal of the Fisheries Research Board of Canada 1979 36 12 1503 1507 10.1139/f79-218 Hunter et al. (2013) Hunter WR Jamieson A Huvenne VAI Witte U Sediment community responses to marine vs. terrigenous organic matter in a submarine canyon Biogeosciences 2013 10 1 67 80 10.5194/bg-10-67-2013 Iken, Bluhm & Gradinger (2005) Iken K Bluhm BA Gradinger R Food web structure in the high Arctic Canada Basin: evidence from δ13C and δ15N analysis Polar Biology 2005 28 3 238 249 10.1007/s00300-004-0669-2 Jackson et al. (2011) Jackson AL Inger R Parnell AC Bearhop S Comparing isotopic niche widths among and within communities: SIBER-Stable Isotope Bayesian Ellipses in R Journal of Animal Ecology 2011 80 3 595 602 10.1111/j.1365-2656.2011.01806.x 21401589 Josefson, Forbes & Rosenberg (2002) Josefson AB Forbes TL Rosenberg R Fate of phytodetritus in marine sediments: functional importance of macrofaunal community Marine Ecology Progress Series 2002 230 71 85 10.3354/meps230071 Jumars, Dorgan & Lindsay (2015) Jumars PA Dorgan KM Lindsay SM Diet of worms emended: an update of polychaete feeding guilds Annual Review of Marine Science 2015 7 1 497 520 10.1146/annurev-marine-010814-020007 Kamenev (2008) Kamenev GM Little-known arctic species Montacuta spitzbergensis (Bivalvia: Montacutidae) from the north-western Pacific with notes on Montacuta substriata and Tellimya ferruginosa Journal of the Marine Biological Association of the United Kingdom 2008 88 2 347 356 10.1017/S0025315408000568 Karaseva et al. (2021) Karaseva NP Rimskaya-Korsakova NN Ekimova IA Kokarev VN Simakov MI Gantsevich MM Malakhov VV The first discovery of pogonophora (Annelida, Siboglinidae) in the east siberian sea coincides with the areas of methane seeps Doklady Biological Sciences 2021 501 1 187 191 10.1134/S0012496621060016 34962604 Karlson, Gorokhova & Elmgren (2015) Karlson AML Gorokhova E Elmgren R Do deposit-feeders compete? Isotopic niche analysis of an invasion in a species-poor system Scientific Reports 2015 5 1 1 8 10.1038/srep09715 Karlson et al. (2010) Karlson AML Nascimento FJA Näslund J Elmgren R Higher diversity of deposit-feeding macrofauna enhances phytodetritus processing Ecology 2010 91 5 1414 1423 10.1890/09-0660.1 20503873 Kędra et al. (2021) Kędra M Cooper LW Silberberger MJ Zhang M Biasatti D Grebmeier JM Organic carbon source variability in Arctic bivalves as deduced from the compound specific carbon isotopic composition of amino acids Journal of Marine Systems 2021 219 April 103547 10.1016/j.jmarsys.2021.103547 Kędra et al. (2015) Kędra M Moritz C Choy ES David C Degen R Duerksen S Ellingsen I Górska B Grebmeier JM Kirievskaya D van Oevelen D Piwosz K Samuelsen A Węsławski JM Status and trends in the structure of Arctic benthic food webs Polar Research 2015 34 1 23775 10.3402/polar.v34.23775 Knies & Martinez (2009) Knies J Martinez P Organic matter sedimentation in the western Barents Sea region: terrestrial and marine contribution based on isotopic composition and organic nitrogen content Norsk Geologisk Tidsskrift 2009 89 1–2 79 89 Kokarev et al. (2021b) Kokarev V Tachon M Austad M McGovern M Reiss H Strong macrobenthic community differentiation among sub-Arctic deep fjords on small spatial scales Estuarine, Coastal and Shelf Science 2021b 252 February 107271 10.1016/j.ecss.2021.107271 Kokarev et al. (2017) Kokarev VNN Vedenin AAA Basin ABB Azovsky AII Taxonomic and functional patterns of macrobenthic communities on a high-Arctic shelf: a case study from the Laptev Sea Journal of Sea Research 2017 129 3 61 69 10.1016/j.seares.2017.08.011 Kokarev et al. (2021a) Kokarev VN Vedenin AA Polukhin AA Basin AB Taxonomic and functional patterns of macrobenthic communities on a high Arctic shelf: a case study from the East Siberian Sea Journal of Sea Research 2021a 174 April 102078 10.1016/j.seares.2021.102078 Kravchishina et al. (2021) Kravchishina MD Lein AY Flint MV Baranov BV Miroshnikov AY Dubinina EO Dara OM Boev AG Savvichev AS Methane-derived authigenic carbonates on the seafloor of the laptev sea shelf Frontiers in Marine Science 2021 8 July 195 10.3389/fmars.2021.690304 Lee et al. (2019a) Lee D Kim J Lee Y Jin Y Paull C Kim D Shin K Chemosynthetic bacterial signatures in Frenulata tubeworm Oligobrachia sp. in an active mud volcano of the Canadian Beaufort Sea Marine Ecology Progress Series 2019a 628 1 95 104 10.3354/meps13084 Lee et al. (2019b) Lee YM Noh HJ Lee DH Kim JH Jin YK Paull C Bacterial endosymbiont of Oligobrachia sp. (Frenulata) from an active mud volcano in the Canadian Beaufort Sea Polar Biology 2019b 42 12 2305 2312 10.1007/s00300-019-02599-w Levin (2005) Levin LA Ecology of cold seep sediments: interactions of fauna with flow, chemistry and microbes Oceanography and Marine Biology 2005 Boca Raton CRC Press 11 56 Levin & Michener (2002) Levin LA Michener RH Isotopic evidence for chemosynthesis-based nutrition of macrobenthos: the lightness of being at Pacific methane seeps Limnology and Oceanography 2002 47 5 1336 1345 10.4319/lo.2002.47.5.1336 Levinton (1972) Levinton J Stability and trophic structure in deposit-feeding and suspension-feeding communities The American Naturalist 1972 106 950 472 486 10.1086/282788 Lösekann et al. (2008) Lösekann T Robador A Niemann H Knittel K Boetius A Dubilier N Endosymbioses between bacteria and deep-sea siboglinid tubeworms from an Arctic Cold Seep (Haakon Mosby Mud Volcano, Barents Sea) Environmental Microbiology 2008 10 12 3237 3254 10.1111/j.1462-2920.2008.01712.x 18707616 Mangiafico (2023) Mangiafico S rcompanion: functions to support extension education program evaluation 2023 R package version 2.4.21 https://CRAN.R-project.org/package=rcompanion Mayor et al. (2012) Mayor DJ Thornton B Hay S Zuur AF Nicol GW McWilliam JM Witte UFM Resource quality affects carbon cycling in deep-sea sediments ISME Journal 2012 6 9 1740 1748 10.1038/ismej.2012.14 22378534 McClelland et al. (2016) McClelland JW Holmes RM Peterson BJ Raymond PA Striegl RG Zhulidov AV Zimov SA Zimov N Tank SE Spencer RGM Staples R Gurtovaya TY Griffin CG Particulate organic carbon and nitrogen export from major Arctic rivers Global Biogeochem. Cycles 2016 30 5 629 643 10.1002/2015GB005351 Middelburg (2018) Middelburg JJ Reviews and syntheses: to the bottom of carbon processing at the seafloor Biogeosciences 2018 15 2 413 427 10.5194/bg-15-413-2018 Morton (1981) Morton B The biology and functional morphology of Periploma (Offadesma) angasai (Bivalvia: Anomaiodesmata: Periplomatidae) Journal of Zoology 1981 193 1 39 70 10.1111/j.1469-7998.1981.tb01490.x Mäkelä, Witte & Archambault (2017) Mäkelä A Witte U Archambault P Ice algae versus phytoplankton: resource utilization by Arctic deep sea macroinfauna revealed through isotope labelling experiments Marine Ecology Progress Series 2017 572 1 18 10.3354/meps12157 Nekhaev & Krol (2020) Nekhaev IO Krol EN Hidden under ice and mud: diversity of shell-bearing microgastropods in the eastern Arctic seas Systematics and Biodiversity 2020 18 8 794 809 10.1080/14772000.2020.1785577 Norling et al. (2007) Norling K Rosenberg R Hulth S Grémare A Bonsdorff E Importance of functional biodiversity and species-specific traits of benthic fauna for ecosystem functions in marine sediment Marine Ecology Progress Series 2007 332 2004 11 23 10.3354/meps332011 North et al. (2014) North CA Lovvorn JR Kolts JM Brooks ML Cooper LW Grebmeier JM Deposit-feeder diets in the Bering Sea: potential effects of climatic loss of sea ice-related microalgal blooms Ecological Applications 2014 24 6 1525 1542 10.1890/13-0486.1 29160671 Oxtoby et al. (2016) Oxtoby LE Budge SM Iken K O’Brien DM Wooller MJ Feeding ecologies of key bivalve and polychaete species in the Bering Sea as elucidated by fatty acid and compound-specific stable isotope analyses Marine Ecology Progress Series 2016 557 161 175 10.3354/meps11863 R Core Team (2022) R Core Team R: a language and environment for statistical computing 2022 Vienna R Foundation for Statistical Computing https://www.r-project.org Reid & Reid (1969) Reid RGB Reid A Feeding processes of members of the genus Macoma (Mollusca: Bivalvia) Canadian Journal of Zoology 1969 47 4 649 657 10.1139/z69-111 Renaud et al. (2015) Renaud PE Løkken TS Jørgensen LL Berge J Johnson BJ Macroalgal detritus and food-web subsidies along an Arctic fjord depth-gradient Frontiers in Marine Science 2015 2 31 10.3389/fmars.2015.00031 Roy et al. (2015) Roy V Iken K Gosselin M Tremblay JÉ Bélanger S Archambault P Benthic faunal assimilation pathways and depth-related changes in food-web structure across the Canadian Arctic Deep Sea Research Part I: Oceanographic Research Papers 2015 102 4 55 71 10.1016/j.dsr.2015.04.009 Sakshaug (2004) Sakshaug E Stein R Macdonald RW Primary and secondary production in the arctic seas The Organic Carbon Cycle in the Arctic Ocean 2004 Berlin Springer 57 82 Sartori & Domaneschi (2005) Sartori AF Domaneschi O The functional morphology of the antarctic bivalve Thracia meridionalis Smith, 1885 (Anomalodesmata: Thraciidae) Journal of Molluscan Studies 2005 71 3 199 210 10.1093/mollus/eyi028 Savvichev et al. (2018) Savvichev AS Kadnikov VV Kravchishina MD Galkin SV Novigatskii AN Sigalevich PA Merkel AY Ravin NV Pimenov NV Flint MV Methane as an organic matter source and the trophic basis of a laptev sea cold seep microbial community Geomicrobiology Journal 2018 35 5 411 423 10.1080/01490451.2017.1382612 Savvichev et al. (2023) Savvichev AS Rusanov II Kadnikov VV Beletsky AV Zakcharova EE Samylina OS Sigalevich PA Semiletov IP Ravin NV Pimenov NV Biogeochemical activity of methane-related microbial communities in bottom sediments of cold seeps of the laptev sea Microorganisms 2023 11 2 250 10.3390/microorganisms11020250 36838215 Schmaljohann et al. (1990) Schmaljohann R Faber E Whiticar M Dando P Co-existence of methane- and sulphur-based endosymbioses between bacteria and invertebrates at a site in the Skagerrak Marine Ecology Progress Series 1990 61 1 119 124 10.3354/meps061119 Schmid et al. (2006) Schmid MK Piepenburg D Golikov AA Juterzenka K Petryashov VV Spindler M Trophic pathways and carbon flux patterns in the Laptev Sea Progress in Oceanography 2006 71 2–4 314 330 10.1016/j.pocean.2006.09.002 Sellanes et al. (2011) Sellanes J Zapata-Hernández G Pantoja S Jessen GL Chemosynthetic trophic support for the benthic community at an intertidal cold seep site at Mocha Island off central Chile Estuarine, Coastal and Shelf Science 2011 95 4 431 439 10.1016/j.ecss.2011.10.016 Sen et al. (2018) Sen A Duperron S Hourdez S Piquet B Léger N Gebruk A Le Port AS Svenning MM Andersen AC Cryptic frenulates are the dominant chemosymbiotrophic fauna at Arctic and high latitude Atlantic cold seeps PLOS ONE 2018 13 12 e0209273 10.1371/journal.pone.0209273 30592732 Sen et al. (2019) Sen A Himmler T Hong WL Chitkara C Lee RW Ferré B Lepland A Knies J Atypical biological features of a new cold seep site on the Lofoten-Vesterålen continental margin (northern Norway) Scientific Reports 2019 9 1762 10.1038/s41598-018-38070-9 30741962 Silberberger et al. (2021) Silberberger MJ Koziorowska-Makuch K Kuliński K Kędra M Stable isotope mixing models are biased by the choice of sample preservation and pre-treatment: implications for studies of aquatic food webs Frontiers in Marine Science 2021 7 1 18 10.3389/fmars.2020.621978 Snelgrove et al. (2018) Snelgrove PVR Soetaert K Solan M Thrush S Wei CL Danovaro R Fulweiler RW Kitazato H Ingole B Norkko A Parkes RJ Volkenborn N Global carbon cycling on a heterogeneous seafloor Trends in Ecology & Evolution 2018 33 2 96 105 10.1016/j.tree.2017.11.004 29248328 Sokołowski et al. (2014) Sokołowski A Szczepańska A Richard P Kędra M Wołowicz M Węsławski JM Trophic structure of the macrobenthic community of Hornsund, Spitsbergen, based on the determination of stable carbon and nitrogen isotopic signatures Polar Biology 2014 37 9 1247 1260 10.1007/s00300-014-1517-7 Southward & Southward (1981) Southward AJ Southward EC Dissolved organic matter and the nutrition of the Pogonophora: a reassessment based on recent studies of their morphology and biology Kieler Meeresforsch 1981 5 445 453 Stasko et al. (2018) Stasko AD Bluhm BA Michel C Archambault P Majewski A Reist JD Swanson H Power M Benthic-pelagic trophic coupling in an Arctic marine food web along vertical water mass and organic matter gradients Marine Ecology Progress Series 2018 594 1 19 10.3354/meps12582 Stein & Fahl (2004) Stein R Fahl K The laptev sea: distribution, sources, variability and burial of organic carbon The Organic Carbon Cycle in the Arctic Ocean 2004 Berlin Springer 213 237 Stein & Macdonald (2004) Stein R Macdonald RW Stein R Macdonald RW Geochemical proxies used for organic carbon source identification in arctic ocean sediments The Organic Carbon Cycle in the Arctic Ocean 2004 Berlin Springer 24 32 Swanson et al. (2015) Swanson HK Lysy M Power M Stasko AD Johnson JD Reist J A new probabilistic method for quantifying n-dimensional ecological niches and niche overlap HEIDI Ecology 2015 96 2 318 324 10.1890/14-0235.1 26240852 Sweetman et al. (2013) Sweetman AK Levin LA Rapp HT Schander C Faunal trophic structure at hydrothermal vents on the southern mohn’s ridge, arctic ocean Marine Ecology Progress Series 2013 473 115 131 10.3354/meps10050 Søreide et al. (2013) Søreide JE Carroll ML Hop H Ambrose WG Hegseth EN Falk-Petersen S Sympagic-pelagic-benthic coupling in Arctic and Atlantic waters around Svalbard revealed by stable isotopic and fatty acid tracers Marine Biology Research 2013 9 9 831 850 10.1080/17451000.2013.775457 Søreide et al. (2006) Søreide JE Hop H Carroll ML Falk-Petersen S Hegseth EN Seasonal food web structures and sympagic-pelagic coupling in the European Arctic revealed by stable isotopes and a two-source food web model Progress in Oceanography 2006 71 1 59 87 10.1016/j.pocean.2006.06.001 Tamelander et al. (2009) Tamelander T Kivimäe C Bellerby RGJ Renaud PE Kristiansen S Base-line variations in stable isotope values in an Arctic marine ecosystem: effects of carbon and nitrogen uptake by phytoplankton Hydrobiologia 2009 630 1 63 73 10.1007/s10750-009-9780-2 Thurber et al. (2010) Thurber AR Kröger K Neira C Wiklund H Levin LA Stable isotope signatures and methane use by New Zealand cold seep benthos Marine Geology 2010 272 1–4 260 269 10.1016/j.margeo.2009.06.001 van der Plas (2019) van der Plas F Biodiversity and ecosystem functioning in naturally assembled communities Biological Reviews 2019 94 4 1220 1245 10.1111/brv.12499 30724447 Vedenin et al. (2020) Vedenin AA Kokarev VN Chikina MV Basin AB Galkin SV Gebruk AV Fauna associated with shallow-water methane seeps in the Laptev Sea PeerJ 2020 8 S1 e9018 10.7717/peerj.9018 32411521 Vonk et al. (2012) Vonk JE Sanchez-Garca L Van Dongen BE Alling V Kosmach D Charkin A Semiletov IP Dudarev OV Shakhova N Roos P Eglinton TI Andersson A Gustafsson A Activation of old carbon by erosion of coastal and subsea permafrost in Arctic Siberia Nature 2012 489 7414 137 140 10.1038/nature11392 22932271 Weems et al. (2012) Weems J Iken K Gradinger R Wooller MJ Journal of experimental marine biology and ecology carbon and nitrogen assimilation in the bering sea clams nuculana radiata and macoma moesta Journal of Experimental Marine Biology and Ecology 2012 430 32 42 10.1016/j.jembe.2012.06.015 Whitlatch (1980) Whitlatch RB Patterns of resource utilization and coexistence in marine intertidal deposit-feeding communities Journal of Marine Research 1980 38 4 743 765 Wickham (2016) Wickham H ggplot2 2016 Berlin Springer International Publishing Włodarska-Kowalczuk et al. (2019) Włodarska-Kowalczuk M Aune M Michel LN Zaborska A Legeżyńska J Is the trophic diversity of marine benthic consumers decoupled from taxonomic and functional trait diversity? Isotopic niches of Arctic communities Limnology and Oceanography 2019 64 5 2140 2151 10.1002/lno.11174 Woelfel et al. (2010) Woelfel J Schumann R Peine F Flohr A Kruss A Tegowski J Blondel P Wiencke C Karsten U Microphytobenthos of arctic kongsfjorden (Svalbard, Norway): biomass and potential primary production along the shore line Polar Biology 2010 33 9 1239 1253 10.1007/s00300-010-0813-0 Woulds et al. (2016) Woulds C Bouillon S Cowie GL Drake E Middelburg JJ Witte U Patterns of carbon processing at the seafloor: the role of faunal and microbial communities in moderating carbon flows Biogeosciences 2016 13 15 4343 4357 10.5194/bg-13-4343-2016 Xiao, Fahl & Stein (2013) Xiao X Fahl K Stein R Biomarker distributions in surface sediments from the Kara and Laptev seas (Arctic Ocean): indicators for organic-carbon sources and sea-ice coverage Quaternary Science Reviews 2013 79 1979 40 52 10.1016/j.quascirev.2012.11.028 Yokoyama & Ishihi (2007) Yokoyama H Ishihi Y Variation in food sources of the macrobenthos along a land-sea transect: a stable isotope study Marine Ecology Progress Series 2007 346 127 141 10.3354/meps07010 Zaborska et al. (2018) Zaborska A Włodarska-Kowalczuk M Legeżyńska J Jankowska E Winogradow A Deja K Sedimentary organic matter sources, benthic consumption and burial in west Spitsbergen fjords–Signs of maturing of Arctic fjordic systems? Journal of Marine Systems 2018 180 1 112 123 10.1016/j.jmarsys.2016.11.005 Zanzerl & Dufour (2017) Zanzerl H Dufour SC The burrowing behaviour of symbiotic and asymbiotic thyasirid bivalves Journal of Conchology 2017 42 5 299 308 Zanzerl et al. (2019) Zanzerl H Salvo F Jones SW Dufour SC Feeding strategies in symbiotic and asymbiotic thyasirid bivalves Journal of Sea Research 2019 145 16 23 10.1016/j.seares.2018.12.005