Use of Stress Signals of Their Attached Bacteria to Monitor Sympagic Algae Preservation in Canadian Arctic Sediments

Based on the strong aggregation of sympagic (ice-associated) algae and the high mortality or inactivity of bacteria attached to them, it was previously hypothesized that sympagic algae should be significant contributors to the export of carbon to Arctic sediments. In the present work, the lipid content of 30 sediment samples collected in the Canadian Arctic was investigated to test this hypothesis. The detection of high proportions of trans vaccenic fatty acid (resulting from cis-trans isomerase (CTI) activity of bacteria under hypersaline conditions) and 10S-hydroxyhexadec-8(trans)-enoic acid (resulting from 10S-DOX bacterial detoxification activity in the presence of deleterious free palmitoleic acid) confirmed: (i) the strong contribution of sympagic material to some Arctic sediments, and (ii) the impaired physiological status of its associated bacterial communities. Unlike terrestrial material, sympagic algae that had escaped zooplanktonic grazing appeared relatively preserved from biotic degradation in Arctic sediments. The expected reduction in sea ice cover resulting from global warming should cause a shift in the relative contributions of ice-associated vs. pelagic algae to the seafloor, and thus to a strong modification of the carbon cycle.


Introduction
In the Canadian Arctic, primary production is supported by sympagic (ice-associated) algae dominated by pennate diatoms and more specifically by the Naviculaceae [1] or by the centric diatom Melosira arctica [2] during the ice-covered period, and then by phytoplankton in open waters [3,4]. The contribution of sympagic algae to total primary production varies widely depending on the season and the region [5,6]. Because of global warming (causing a decrease in sea ice extent and duration), we are currently witnessing a reduced primary production of the sympagic algae and an increase in that of pelagic phytoplankton. Sympagic algae are assumed to be a main source of organic matter reaching the seafloor [2,[7][8][9]. It is therefore feared that in the future, the biological pump would have a positive feedback effect on global warming. The high contribution of sympagic algae to Arctic sediments was previously attributed to: (i) their strong aggregation resulting from the high concentrations of extracellular polymeric substances (EPSs) produced by these organisms in the ice, which protects biogenic silica of diatom frustules from dissolution [10] and induces quick settling in the water column (100-500 m d −1 ; [11]); (ii) the fact that aggregated microalgae are not favored prey particles for grazers [12]; and (iii) the poor physiological status of their associated bacterial communities (inducing weak mineralization [7] and biogenic silica dissolution [13] in the water column and surface sediments). These bacterial communities are damaged by: (i) intense osmotic stress induced by changes portions of highly branched isoprenoids (HBIs; IP 25 and tri-unsaturated HBIs) to estimate relative proportions of sympagic vs. pelagic production [40][41][42][43], some uncertainties remain about the nature of POM reaching the seafloor. In particular, processes affecting the sinking organic matter and its lipid content through the water column down to the seafloor (such as degradation and remineralization) need further investigation. The aim of this work was thus to use new indicators of the osmotic and chemical stress status of bacteria associated with sympagic material (trans-vaccenic acid and 10S-hydroxyhexadec-8(trans)-enoic acid, respectively) to estimate: (i) the contribution of the sympagic material to the investigated sediments and (ii) the status of the bacteria associated with it. Ratios of stanols (widespread bacterial biohydrogenation products of sterols [44,45]) to their parent sterols were also used to confirm the impaired ability of stressed bacteria to degrade sympagic material.
Microorganisms 2021, 9, x FOR PEER REVIEW 3 of 16 To confirm this expected strong contribution of sympagic material to the seafloor, we examined the lipid content of 30 sediment samples collected in the Canadian Arctic (Figure 1). Despite the recent development of several indices based on the relative proportions of highly branched isoprenoids (HBIs; IP25 and tri-unsaturated HBIs) to estimate relative proportions of sympagic vs. pelagic production [40][41][42][43], some uncertainties remain about the nature of POM reaching the seafloor. In particular, processes affecting the sinking organic matter and its lipid content through the water column down to the seafloor (such as degradation and remineralization) need further investigation. The aim of this work was thus to use new indicators of the osmotic and chemical stress status of bacteria associated with sympagic material (trans-vaccenic acid and 10S-hydroxyhexadec-8(trans)-enoic acid, respectively) to estimate: (i) the contribution of the sympagic material to the investigated sediments and (ii) the status of the bacteria associated with it. Ratios of stanols (widespread bacterial biohydrogenation products of sterols [44,45]) to their parent sterols were also used to confirm the impaired ability of stressed bacteria to degrade sympagic material.

Sediment Sampling
Two contrasting regions were selected, namely the Canadian Beaufort Shelf and Baffin Bay. The Beaufort Shelf is a perennially stratified interior shelf influenced by Pacific-

Sediment Sampling
Two contrasting regions were selected, namely the Canadian Beaufort Shelf and Baffin Bay. The Beaufort Shelf is a perennially stratified interior shelf influenced by Pacific-derived waters supplied via the Beaufort Gyre and the Alaskan coastal current [46]. Sediments of this seasonally ice-covered shelf are strongly influenced by the Mackenzie River [47], the largest river draining into the Arctic in terms of sediment and POM [48]. The Mackenzie shelf functions as a vast estuary receiving inputs of both terrestrial and marine sources of organic matter [49]. By contrast, Baffin Bay is characterized by a strong spatial sea ice variability resulting from the inflow of the West Greenland Current (WGC) composed of the relatively warm Atlantic Irminger Current (IC), which restricts sea ice extent in its eastern part [50]. The sedimentation rates appeared to be relatively similar in the two zones investigated ranging from 0.04 to 0.20 cm yr −1 in the Beaufort Shelf [51,52] and from 0.06 to 0.11 cm yr −1 in Baffin Bay [52]. Depths of the different sampling stations (ranging from 7 to 2017 m) are given in Supplementary Table S2.
Sediment samples were collected from a broad range of locations ( Figure 1) within the Canadian Arctic using an USNEL box corer (50 [54][55][56]) seemed relatively unaffected. Moreover, similar values were measured in the Amundsen Gulf (ranging from 2.4 and 6.8 g C m −2 yr −1 [56,57]) and in Baffin Bay (ranging from 1.1 to 6.7 g C m −2 yr −1 [56,58]) justifying a comparative study of the different sediment samples investigated. From each box core, one sample of ca. 50 cm 2 was collected from intact sediment surface (0-1 cm) and frozen immediately at −80 • C for later analysis.

Lipid Analysis
Samples were reduced with excess NaBH 4 after addition of MeOH (25 mL, 30 min) to reduce labile hydroperoxides to alcohols, which are more amenable to analysis using gas chromatography-mass spectrometry (GC-MS). Water (25 mL) and KOH (2.8 g) were then added and the resulting mixture saponified (to break complex lipids down into their constituent fatty acids) by refluxing (2 h). After cooling, the mixture was acidified (HCl, 2 N) to pH 1 and extracted with dichloromethane (DCM; 3 × 20 mL). The combined DCM extracts were dried over anhydrous Na 2 SO 4 , filtered and concentrated by rotary evaporation at 40 • C to give total lipid extracts (TLEs). Aliquots of TLEs were either silylated and analyzed by gas chromatography-electron impact quadrupole time-of-flight mass spectrometry (GC-QTOF) for monounsaturated fatty acid oxidation product quantification, or methylated, then treated with dimethyldisulfide (DMDS) and analyzed by GC-MS/MS to determine double-bond stereochemistry. Cis and trans isomers of monounsaturated fatty acid methyl esters react with DMDS stereospecifically to form threo and erythro adducts, which exhibit similar mass spectra, but are well-separated by gas chromatography, allowing unambiguous double-bond stereochemistry determination [59].

Derivatization
The derivatization method used was silylation allowing the replacement of the active hydrogen atom of carboxylic and alcoholic groups of acids and hydroxyacids by a trimethylsilyl group in one step. TLEs were silylated by dissolving them in 300 µL of a mixture of pyridine and BSTFA (N,O-Bis(triméthylsilyl)trifluoroacétamide; Supelco; 2:1, v/v) and heating to 50 • C (1 h). After evaporation to dryness under a stream of N 2 , the derivatized residue was dissolved in a mixture of hexane and BSTFA (to avoid desilylation) and analyzed by GC-MS/MS or GC-QTOF.

Determination of Double-Bond Stereochemistry
TLEs were dissolved in 2 mL of BF 3 /methanol (10%) (Sigma-Aldrich, St. Louis, MO, USA) and heated at 80 • C (1 h) in a screw-cap flask to obtain fatty acid methyl esters (FAMEs). After cooling, an excess of water was added and FAMEs were extracted three times with hexane, dried over anhydrous Na 2 SO 4 , filtered on Whatman cellulose filters (diameter 90 mm, porosity 11 µm), concentrated using rotary evaporation, and transferred to screw-cap flasks. After solvent removal (N 2 ), 200 µL of DMDS (Sigma-Aldrich, St. Louis, MO, USA) and 50 µL of iodine solution (60 µg µL −1 in diethyl ether) were added. The mixtures were shaken and heated at 50 • C (48 h), excess iodine was removed by addition of 2 mL of a 5% Na 2 S 2 O 3 solution, and lipids were extracted three times with 1 mL of hexane. The extracts were dried with anhydrous Na 2 SO 4 , filtered, and concentrated prior to analysis by GC-MS/MS.

Gas Chromatography-EI Quadrupole Time of Flight Mass Spectrometry
Accurate mass measurements were carried out in full scan mode using an Agilent 7890B/7200 GC-QTOF System (Agilent Technologies, Parc Technopolis-ZA Courtaboeuf, Les Ulis, France). A cross-linked 5% phenyl-methylpolysiloxane (Agilent; HP-5MS ultra inert) (30 m × 0.25 mm, 0.25 µm film thickness) capillary column was used. Analyses were performed with an injector operating in pulsed splitless mode set at 270 • C and the oven temperature was programmed from 70 • C to 130 • C at 20 • C min −1 and then to 300 • C at 5 • C min −1 . The pressure of the carrier gas (He) was maintained at 0.69 × 10 5 Pa until the end of the temperature program. Instrument temperatures were 300 • C for the transfer line and 230 • C for the ion source. Nitrogen (1.5 mL min −1 ) was used as collision gas. Accurate mass spectra were recorded across the range m/z 50-700 at 4 GHz with the collision gas opened. The QTOF-MS instrument provided a typical resolution ranging from 8009 to 12,252 from m/z 68.9955 to 501.9706. Perfluorotributylamine (PFTBA) was used for daily MS calibration. Compounds were identified by comparing their TOF mass spectra, accurate masses, and retention times with those of standards. Each compound was quantified by extraction of specific accurate fragment ions, peak integration, and determination of individual response factors using external standards and Mass Hunter software (Agilent Technologies, Parc Technopolis-ZA Courtaboeuf, Les Ulis, France).

Statistical Analysis
The variable investigated being non-parametric, Spearman correlations were performed to determine the correlation between the depth and the trans/cis ratio or 10S-DOX degradation percentage. Mann-Whitney-Wilcoxon tests were performed to identify any significant differences in (i) trans/cis ratio and 10S-DOX degradation percentage between Baffin Bay and the Beaufort Sea and (ii) sterol:stanol ratio percentage between sitosterol and brassicasterol.
Trans/cis ratios of vaccenic acid were measured after DMDS treatment in the different surface sediment samples investigated ( Figure 2). The results obtained ( Figure 1, Table 1,  and Supplementary Table S2) showed significantly higher values in the Beaufort Sea than in Baffin Bay (mean ± SE = 0.20 ± 0.03 and 0.06 ± 0.02 (g:g), respectively; W = 210, p < 0.01). The efficiency of sympagic-benthic coupling should logically be better in shallower zones. This assumption is well supported by the good anticorrelation observed between trans/cis ratio and depth of the different stations (Spearman's rho = −0.67, p < 0.01).

Statistical Analysis
The variable investigated being non-parametric, Spearman correlations were performed to determine the correlation between the depth and the trans/cis ratio or 10S-DOX degradation percentage. Mann-Whitney-Wilcoxon tests were performed to identify any significant differences in (i) trans/cis ratio and 10S-DOX degradation percentage between Baffin Bay and the Beaufort Sea and (ii) sterol:stanol ratio percentage between sitosterol and brassicasterol.
Trans/cis ratios of vaccenic acid were measured after DMDS treatment in the different surface sediment samples investigated ( Figure 2). The results obtained ( Figure 1, Table 1,  and Supplementary Table S2) showed significantly higher values in the Beaufort Sea than in Baffin Bay (mean ± SE = 0.20 ± 0.03 and 0.06 ± 0.02 (g:g), respectively; W = 210, p < 0.01). The efficiency of sympagic-benthic coupling should logically be better in shallower zones. This assumption is well supported by the good anticorrelation observed between trans/cis ratio and depth of the different stations (Spearman's rho = −0.67, p < 0.01). A strong contribution sympagic material to the sediments of the Beaufort Sea is well supported by the high ratios previously measured in surface sediments of other stations in this zone (0.65 ± 0.15 g:g, n = 6) [39]. It was previously observed that bacterial osmotic stress resulting from hypersaline conditions in brine channels occurred only during the early stages of ice melting and not at the end of the melting season [7,14]. The sympagic material reaching the seafloor of the Beaufort Sea thus seems to be discharged from sea ice at the beginning of melting. The high trans/cis ratios observed ( Figure 1; Table 1) ruled out a significant contribution of sympagic algae released in the water column at the advanced stages of ice melting and open water phytoplankton to the sediments. A strong contribution of such material (exhibiting very low and no CTI activity, respectively [14,69]) to the sediment should strongly lower the trans/cis ratios (by adding bacteria unstressed by salinity with trans/cis ratios < 0.1). This assumption is firmly supported by the results of Morata et al. [70], who observed that inputs of sympagic algae to the sediment started to increase from January/February in the Beaufort Sea, while sympagic algal bloom occurred between mid-March and late May [71]. Juul-Pedersen et al. [72] also observed a sinking export of sympagic algal material in the same zone during winter and early spring (i.e., well before the onset of spring melt). The amount of diatoms exported to the seafloor strongly depends on a match or mismatch between algal production and zooplankton grazing [73]. It is generally considered that copepods do not feed at chlorophyll concentrations < 1 µg L −1 [74]. The very low chlorophyll concentrations measured in POM between February and March in the Beaufort Sea (0.002-0.003 µg L −1 [70]) should thus not favor zooplankton grazing on sympagic algae during this period. We note that an increase in zooplankton fecal pellet production was observed during this same period [75]; however, those authors attributed it to the grazing of non-pigmented sources of food, such as microzooplankton.
Chlorophyll concentration increases strongly in surface waters in late spring and summer at the end of ice melting and during open water diatom bloom (up to 5 µg L −1 [76]) favoring grazing of sympagic [77] and pelagic algae by copepods. Coprophagy can then reduce losses by sinking of the resulting fecal pellets at only a few percent of the pellet production rate [78]. Bacterial decomposition by non-stressed internal bacteria is another mechanism that could prevent pellets sinking to the seafloor [79,80]. Forest et al. [81] previously estimated that 97% of the primary-produced C was grazed or degraded in the water column of the Amundsen Gulf. The remaining 3% reaching the seafloor thus seem mainly composed of sympagic algae associated with bacteria stressed by salinity and released during the early stages of ice melting. This view is well supported by earlier observations of Morata and Renaud [82]. These authors analyzed sedimentary pigments in the Beaufort Sea and concluded that in the spring, sympagic-algal production largely influenced organic matter inputs to the benthos, while in the summer grazing was responsible for inputs of degraded material.
Whereas a high CTI activity was previously observed in sea ice and sinking particles during the early stages of ice melting and in sediments from a landfast ice station located at Qikiqtarjuaq near Broughton Island (western Baffin Bay, [7]), trans/cis ratios were found to be particularly weak in middle and eastern Baffin Bay (Figure 1). Interestingly, Yunda-Guarin et al. [8] found more sea ice-derived particulate organic carbon in surface sediments in the western side than in the eastern side of Baffin Bay. These authors attributed this difference to the timing of sea ice retreat. Ice cover is greater in the western part of Baffin Bay than in its eastern side, which is influenced by warmer waters from the western streams of Greenland [83]. Moreover, recently calculated ecological network analysis indices reveal that the complex eastern Baffin Bay food web favors the classical grazing chain, while the shorter western food web induces a higher carbon export [84]. Based on these results, the contrasting high trans/cis ratio observed in sediments from Qikiqtarjuaq (Figure 1) may be Microorganisms 2021, 9, 2626 8 of 15 due to a large contribution of sympagic material containing bacteria stressed by salinity. This hypothesis is supported by the similarity of trans/cis ratios previously observed at this station in sinking particles and in surface sediments [7]. The very low ratios observed in the other samples from eastern Baffin Bay (Figure 1) are in close agreement with the higher contributions of pelagic carbon observed by Yunda-Guarin et al. [8] in these zones.
10S-Hydroxyhexadec-8(trans)-enoic acid arising from NaBH 4 reduction in the corresponding hydroperoxide was detected in most of the sediments investigated (Figures 1 and 3; Supplementary Table S2). The percentages of this hydroxyacid (relative to the parent palmitoleic acid and its degradation products) observed in Baffin Bay and the Beaufort Sea were not significatively different (mean ± SE = 15.20 ± 2.80 and 13.64 ± 1.80%, respectively; W = 130.5, p = 0.69). The presence of this compound is indicative of bacterial 10S-DOX oxidation of palmitoleic acid [85]. This enzymatic activity seems to be a detoxification strategy enabling sympagic bacteria to survive the production of deleterious free palmitoleic acid [33,86] by sympagic algae [14]. We note that despite this detoxification strategy a high mortality of attached bacteria could previously be observed in sea ice during the sympagic bloom [14]. It was recently observed that 10S-DOX could also be employed by some pelagic bacteria to detoxify free palmitoleic acid released by wounded diatoms in the presence of copepods [38]. The lack of correlation observed between the percentage of 10Shydroxyhexadec-8(trans)-enoic acid and depth (Spearman's rho = −0.22, p = 0.25) suggests a mixed contribution of sympagic and pelagic bacteria to the sedimentary 10S-DOX signal. Although 10S-hydroxyhexadec-8(trans)-enoic acid is thus not sufficiently specific for use as a marker of the contribution of sympagic material to Arctic sediments, it is indicative of the presence of bacteria whose membranes are altered by FFAs. The bacterial communities present in Arctic sediments are thus in part non-growing (i.e., stressed by salinity in brine channels of ice) or dead (i.e., stressed by FFAs in ice or in the water column).  It is generally considered that the deposited algal aggregates are overgrown wi a few weeks to months by specific bacterial groups of the surrounding sediment [ which are likely better adapted to deep-sea environment than surface-derived bact [88]. In the Arctic, these benthic bacteria are dominated by members of the Roseoba clade (i.e., the genera Neptunomonas, Arcobacter, and Sedimentitalea [87]), which all con large proportions of cis-vaccenic acid [89][90][91]. The lack of dilution of trans/cis ratios of acid observed in surface sediments from the Beaufort Sea and northern and western Ba Bay (Figure 1; Table 1) is thus very surprising. In these zones, sympagic-derived ag It is generally considered that the deposited algal aggregates are overgrown within a few weeks to months by specific bacterial groups of the surrounding sediment [87], which are likely better adapted to deep-sea environment than surface-derived bacteria [88]. In the Arctic, these benthic bacteria are dominated by members of the Roseobacter clade (i.e., the genera Neptunomonas, Arcobacter, and Sedimentitalea [87]), which all contain large proportions of cis-vaccenic acid [89][90][91]. The lack of dilution of trans/cis ratios of this acid observed in surface sediments from the Beaufort Sea and northern and western Baffin Bay ( Figure 1; Table 1) is thus very surprising. In these zones, sympagic-derived aggregates reaching the seafloor do not appear to be intensively colonized by benthic bacteria. This absence of colonization may be attributed to the presence of bactericidal FFAs (notably palmitoleic acid) in these aggregates attested by the 10S-DOX activity observed in most of them (Figure 1; Supplementary Table S2).

Impact of Stress State of Bacteria on the Degradation of Sympagic Material in Sediments
It is generally considered that the lifetime of organic compounds in marine sediments depends on environmental conditions such as bioturbation, physical mixing, and the presence or absence of oxygen and other electron acceptors [92]. The benthic fauna inhabiting the seafloor is known to: (i) modify the vertical zonation of respiration reactions, (ii) control microbial community assembly, (iii) change the distribution of organic matter (OM), and (iv) influence rates of microbial OM remineralization in surface sediments [93]. Preservation/degradation of OM also depends on its sources [94]; terrestrial vascular debris generally contain greater proportions of refractory organic carbon than algal material, owing to their geochemical composition [95]. The non-growing or dead state of bacteria associated with sympagic algae should strongly impact the degradation conditions of this material in Arctic sediments. Since the conversion of ∆ 5 -stenol to stanol is generally considered as indicating an intense bacterial degradation [96,97], we compared the stanol/∆ 5 -stenol ratios of 24-methylcholesta-5,22-dien-3β-ol (brassicasterol) (common sterol of sympagic diatoms; [98]) and 24-ethylcholest-5-en-3β-ol (sitosterol) (sterol present in some diatoms but also widely distributed in terrestrial vascular plants [99]) in sediment samples from the Beaufort Sea where the highest bacterial stresses were observed (Table 1). We note that in these sediments sitosterol is known to result mainly from terrigenous material [64,100]. The stanol/∆ 5 -stenol ratios obtained for the two sterols were significantly different (mean ± SE = 18.2 ± 1.9% and 52.8 ± 5.9%, for brassicasterol and sitosterol respectively; W = 36, p < 0.01, Table 2). In the case of brassicasterol, the values observed were close to those found in healthy phototrophic organisms (5-20% [101,102]), attesting to the good biotic preservation of the algal material. By contrast, the measured sitostanol/sitosterol ratio reflects a strong biodegradation of this terrigenous sterol and supports the results of Gõni et al. [103], who estimated that 65% of the old and fossil carbon is respired or buried in the Mackenzie delta. On entering the sea, terrestrial organic matter thus seems to undergo strong bacterial remineralization. This could be attributed to the involvement of "priming effects" (enhanced remineralization of terrestrial OM in the presence of fresh substrates from pelagic algal sources; [104][105][106]). By contrast, aggregated sympagic material sinking quickly through the water column with bacteria of poor health status appears well-preserved in surface sediments. Under warmer conditions, the Arctic carbon cycle, which is mostly driven by the primary production (sympagic and pelagic algae) and the riverine inputs (including permafrost [107]), will be impacted. Although the contributions of these sources to the exported carbon have recently been estimated [2,[108][109][110], their burial efficiency, impacting Arctic CO 2 storage remains uncertain. Here, we monitored their biotic behavior with specific lipid tracers in several surface sediment samples collected in the Canadian Arctic. We show that in the Beaufort Sea and northern and western Baffin Bay, unlike phytoplankton and permafrost, sympagic algae organic matter released during the early stages of ice melt by brine drainage contributes significantly to the exported and buried carbon, owing to: (i) the weak zooplanktonic grazing activity during this period, (ii) its strong aggregation enhancing rates of sinking to the seafloor, and (iii) its recalcitrance to demineralizing processes. This particular resistance to biodegradation results from the non-growing (inactive) and dead state of bacteria associated with this sympagic material (Figure 4).  Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Table S1: Depth, coordinates, and concentrations (µg g −1 ) of the main lipids of interest in the different surface sediments investigated, Table S2  We note that Koch et al. [111] previously estimated fluxes of sympagic diatoms in northern Bering and Chukchi Seas with HBI tracers and observed the highest fluxes in July (i.e., at the end and not in the early stages of ice melting). These observations contrasting with ours may be attributed to the particularly strong pelagic-benthic coupling of northern Bering and Chukchi Seas resulting from low grazing pressure due to temporal mismatch between zooplankton and sympagic algae production [112] and the shallowness of these zones [113]. However, an effect of the climate of specific years when the research was carried out cannot be totally excluded. Match or mismatch of zooplanktonic grazing with sympagic algae fluxes seems to be a true key factor controlling the export of sympagic material to the seafloor (Figure 4). In the future, it will be of interest to analyze several Arctic sediments to determine whether the conclusions of this work can be extended to other Arctic Zones. During these comparative studies, it will be necessary to use the same research protocol (i.e., the same tracers).