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

HAL
open science

# Metalloenzyme signatures in authigenic carbonates from the Chukchi Borderlands in the western Arctic Ocean

Dong-Hun Lee, Jung-Hyun Kim, Yung Mi Lee, Germain Bayon, Dahae Kim, Young Jin Joe, Xudong Wang, Kyung-Hoon Shin, Young Keun Jin

► To cite this version:

Dong-Hun Lee, Jung-Hyun Kim, Yung Mi Lee, Germain Bayon, Dahae Kim, et al., Metalloenzyme signatures in authigenic carbonates from the Chukchi Borderlands in the western Arctic Ocean. Scientific Reports, 2022, 12 (1), 16597 (15p.). 10.1038/s41598-022-21184-6. hal-04203895

HAL Id: hal-04203895

https://hal.science/hal-04203895v1

Submitted on 6 Feb 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

OPEN

# Metalloenzyme signatures in authigenic carbonates from the Chukchi Borderlands in the western Arctic Ocean

Dong-Hun Lee<fnref n="1" />,⁵, Jung-Hyun Kim²✉, Yung Mi Lee², Germain Bayon³, Dahae Kim², Young Jin Joe², Xudong Wang⁴, Kyung-Hoon Shin<fnref n="1" /> & Young Keun Jin²

Migration of methane-rich fluids at submarine cold seeps drives intense microbial activity and precipitation of authigenic carbonates. In this study, we analyzed microbially derived authigenic carbonate samples recently recovered from active gas hydrate mounds on the southwestern slope of the Chukchi Borderlands (CB), western Arctic Ocean. Our main aim was to characterize the distribution patterns of trace elements in carbonate-hosted lipid fractions to assess metalloenzyme requirements of microbes involved in anaerobic oxidation of methane (AOM). We measured stable isotopes, trace elements, lipid biomarkers, and genomic DNA, and results indicate the dominance of AOM-related lipid biomarkers in studied carbonate samples, as well as a predominant occurrence of the anaerobic methanotrophic archaea (ANME)-1. We also report evidence for significant preferential enrichments of various trace elements (Li, Ni, Co, Cu, Zn, and Mo) in the total lipid fractions of CB carbonates, relative to elemental compositions determined for corresponding carbonate fractions, which differ from those previously reported for other seep sites. We hypothesize that trace element enrichments in carbonate-hosted lipid fractions could vary depending on the type of AOM microbial assemblage. Additional work is required to further investigate the mechanisms of lipid-bound trace elements in cold seep carbonates as potential metalloenzymes in AOM.

Cold seeps occur at ocean margins worldwide<fnref n="1" />,², corresponding to methane-rich fluid emanations from the seabed into the water column³,⁴. The seepage rates of reduced fluids can vary substantially at the seafloor⁵, leading to a large spectrum of redox conditions in subsurface sediments at cold seeps⁶. A key biogeochemical process in submarine methane seeps is the anaerobic oxidation of methane (AOM), which efficiently consumes a substantial fraction of methane released into the overlying water column⁷,⁸. During AOM, alkalinity levels strongly increase in the surrounding pore waters because of the production of biocarbonate ($HCO_3^-$), promoting carbonate supersaturation and precipitation²,⁹. Accordingly, authigenic carbonate precipitation is typically encountered in cold seeps¹⁰,¹¹, which can serve as an archive for investigating methane release events in both modern and ancient seeps¹²,¹³.

AOM is typically mediated by a consortium of anaerobic methanotrophic archaea (ANME) and sulfate-reducing bacteria (SRB)¹⁴,¹⁵. However, recent studies have also reported that ANMEs can also use other electron acceptors, such as manganese (Mn)- and iron (Fe)-rich oxyhydroxides, to promote methane oxidation in cold seeps¹⁶. Notably, other trace metals such as nickel (Ni), cobalt (Co), molybdenum (Mo), tungsten (W), and zinc (Zn) can also be involved in AOM, acting as enzymatic co-factors¹⁷,¹⁸. Most of these findings have been obtained from culture experiments¹⁷,¹⁹. A recent study characterized, for the first time, the trace element geochemistry of the total lipid fractions preserved in authigenic carbonates from various seeps worldwide, such as the Congo, Nile deep-sea, and Niger fans, and the Gulf of Mexico²⁰, suggesting that such an approach could be used to identify trace metals essential to microbial activity in cold seeps. A previous study reported marked enrichment of Ni, Co, Mo, and W in carbonate-hosted lipid fractions, which were linked to previously identified enzymatic

<footnote n="1">¹Department of Marine Sciences and Convergent Technology, Hanyang University ERICA Campus, Ansan 15588, Republic of Korea. ²Korea Polar Research Institute, Incheon 21990, Republic of Korea. ³CNRS, Ifremer, Brest University, Geo-Ocean, 29280 Plouzané, France. ⁴Shanghai Engineering Research Center of Hadal Science and Technology, College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China. ⁵Present address: Marine Environment Research, Division, National Institute of Fisheries Science, Busan, Republic of Korea. ⁶email: jhkim123@kopri.re.kr</footnote>
pathways involved in AOM²⁰. Nonetheless, information about how AOM may be influenced by changes in trace metal bioavailability in cold seeps remains scarce.

Recently, gas hydrate mound structures have been discovered in the Chukchi Borderlands (CB) of the western Arctic Ocean²¹. Geochemical information on the origin of the emitted gas and pore fluid properties was obtained within this mound structure during the R/V ARAON expedition ARA09C in 2018²². However, thus far, biogeochemical signatures preserved in authigenic carbonates have not been investigated in these newly discovered gas hydrate mounds. In this study, we provide a detailed geochemical characterization of a series of authigenic carbonate samples using oxygen and carbon stable isotopes and mineralogy, as well as trace element abundances, lipid biomarkers, and nucleic acids. The main objectives of this study were to i) evaluate key environmental factors controlling carbonate precipitation in the CB mounds, ii) characterize geochemical and microbial community characteristics in the CB mounds, and iii) assess the relation of trace metals in microbial activity associated with AOM by comparing our new data for CB geochemical and microbial signatures with those previously obtained at other seeps worldwide²⁰. We aimed to test the hypothesis that specific trace element enrichments in the total lipid fractions extracted from authigenic carbonates could be indicative of preferential trace metal utilization by different ANME communities. Our study sheds new light on the metalloenzyme requirements of microbes involved in AOM.

## Results

**Mineralogy and stable isotopic compositions.** The authigenic carbonate samples collected in this study (Fig. 1, see also Fig. S1 in the Supplementary Information) were dominated by calcite, with a minor contribution from dolomite (Table 1). The δ¹³C values of carbonates (δ¹³C<sub>carbonate</sub>) were depleted, varying between -37.0 and -32.8‰. The δ¹⁸O values of carbonates (δ¹⁸O<sub>carbonate</sub>) ranged from 3.9 to 5.6‰, with more enriched values occurring in samples from the deeper sediment layers (Table 1). The measured δ¹³C<sub>carbonate</sub> and δ¹⁸O<sub>carbonate</sub> values of bulk CB carbonates partially overlapped with those determined for a series of seep carbonate samples previously investigated²⁰ (Table 1 and Fig. S2<sup>10,23–27</sup>).

**Trace element compositions.** Trace element data for carbonate (1 M AA leachates), sulfide (3 M HNO₃ leachates), detrital silicate (HF-HCl digestion), and lipid fractions are reported in Tables S1–S4 in the Supplementary Information. Compared to JLs-1 (carbonate reference material; Triassic limestone), the carbonate fractions investigated in this study were generally characterized by much higher trace element abundances (Fig. 2 and Fig. S3), many of which (lithium (Li), transition metals, strontium (Sr), REE, lead (Pb), and thorium (Th)) were enriched up to a few hundred times. In contrast, three elements (titanium (Ti), barium (Ba), and W) showed relatively lower concentrations (up to tenfold depletion compared to JLs-1).

Trace element data for the extracted sulfide-rich solid phases (i.e., 3 M HNO₃ leachates) were normalized to the data for corresponding carbonate fractions, indicating that some elements (Ti, copper (Cu), zirconium (Zr), and Th) in the carbonate and sulfide fractions were enriched as much as 10 to 100 times (Fig. S3). Trace element abundances in the detrital silicate fractions were also normalized to the Post Archean Australian Shale (PAAS) reference values²⁸. Some elements, such as Li, Rb (rubidium), Ba, W, Pb, Th, and U, were relatively enriched in the detrital silicate fractions relative to PAAS (Fig. S3). For the lipid fractions, several trace elements displayed a large range of abundance values (Li: 2.8–29.9 ppm, transition metals: 0.5–300 ppm, Rb: 2.2–4.1 ppm, Sr: 12.0–45.0 ppm, Ba: 13.0–22.0 ppm, Tables S1–S4). In contrast, REEs displayed very low abundances (<1.34 ppm) in the lipid fractions. Notably, the lipid fraction data normalized to JLs-1 and the carbonate fractions of the CB carbonates clearly showed marked enrichment of Ti, Ni, Cu, Zn, and Pb (as much as 10–1000 times; Fig. 2, see also Fig. S4).

**Lipid biomarker characteristics.** Concentrations of archaeal and bacterial lipids showed large variations in the CB authigenic carbonates (Fig. 3, see also Table S5). Among isoprenoid hydrocarbons, the irregular tail-to-tail isoprenoids, such as crocetane, are often coeluted with phytane. Considering that the isotopic compositions of the crocetane/phytane peaks in CB carbonates were with -28.5 to -25.4‰ (Table S5), i.e., not ¹³C-depleted, it appears that crocetane was nearly absent in our samples, but phytane derived from the degradation of chlorophyll and other pigments was present with the values of 0.01–0.05 µg g⁻¹ dry weight (dw). The C₂₅ compound 2,6,10,15,19-pentamethylicosane (PMI) was detected in the range of 0.02–0.04 µg g⁻¹ dw. Moreover, isoprenoid dialkyl glycerol diethers (isoprenoid DGDs), such as archaeol and *sn*-2-hydroxyarchaeol, were the most predominant archaeal lipids in the samples, displaying concentration ranges of 0.02–0.18 µg g⁻¹ dw and 0.02–0.32 µg g⁻¹ dw, respectively. For bacterial lipids, non-isoprenoid DGDs (i.e., If, IIa, and IIId) with anteiso pentadecyl moieties or cyclopropyl groups attached at both the *sn*-1 and *sn*-2 positions were identified, ranging between 0.01 and 0.11 µg g⁻¹ dw. With respect to other DGDs (series I, II, and III), we could not fully determine their presence due to the detection limit. Among fatty acids (FAs) detected, saturated FAs (e.g., C16:0 and C18:0) were the most predominant, ranging between 0.04 and 0.62 µg g⁻¹ dw (Table S5). Other FAs (e.g., i-C15:0, *ai*-C15:0, C16:1ω7, C18:1ω9, C18:1ω7, C20:0, and C22:0) were approximately 2 to 10 times lower in concentration (0.01–0.23 µg g⁻¹ dw) compared with the predominant saturated FAs (Table S5). The δ¹³C values of PMI ranged from -88.2 to -61.5‰ (Fig. 3, see also Table S5). The δ¹³C values of archaeol and *sn*-2-hydroxyarchaeol varied between -102.9 and -82.9‰ and between -108.0 and -87.5‰, respectively (Fig. 3, see also Table S5). Finally, the δ¹³C values of non-isoprenoid DGDs ranged from -88.6 to -50.6‰ (Fig. 3, see also Table S5).

**Microbial compositions.** The archaeal phyla detected in the studied samples were Euryarchaeota, Lok archaeota, Bathyarchaeota, and Thaumarchaeota (Table S6), with Euryarchaeota having the highest abundance, accounting for > 85% of the total archaeal sequences. Among them, the relative abundance of archaea ASV032,
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Map showing (a) the mound structures discovered in the Chukchi Borderlands (CB) with authigenic carbonates investigated in this study and (b) four comparison study sites: Congo fan (CF), Nile deep-sea fan (NDSF), Niger fan (NF), and Gulf of Mexico (GoM). The bathymetric map was generated with the International Bathymetric Chart of the Arctic Ocean Version 3.0 (https://www.ngdc.noaa.gov/mgg/bathymetry/arctic). The sub-bottom profile (SBP) image was acquired using a SBP120 equipment (frequency; 12 kHz, wide-beam angle; 90–120°, Kongsberg). The global ocean map was generated with the Ocean Data View version 5.1.4 (GlobalHR, https://odv.awi.de).</figcaption></figure>

which belongs to ANME-1, was dominant (up to 99.5%; Fig. 4). Another archaea, ASV095, which belongs to *Methanosarcinales*, comprised 0.5% to 64.4% of archaeal sequences. In the bacterial communities, Firmicutes, Proteobacteria, Atribacteria, Actinobacteria, Bacteroidetes, Planctomycetes, and WM88 were predominant (Table S6).

## Discussion

**Constraints on fluid sources and environmental conditions during carbonate precipitation.** Magnesium (Mg) calcite was the most dominant mineral (23–46 wt. %) in all carbonate samples retrieved from the CB mounds (Table 1). Generally, the environmental conditions during which carbonate precipitation occurs at cold seeps can be inferred from the mineralogical composition of carbonates²³,²⁹. For instance, aragonite precipitation is believed to be favored at sub-oxic conditions characterized by high dissolved sulfate concentrations in pore waters (i.e., near the seafloor), while high-Mg calcite and dolomite form preferentially under more reducing conditions and low sulfate levels³⁰. Mg-bearing carbonate minerals such as dolomite are often encountered in relatively deeply buried sulfate-depleted sediment layers³¹. By analogy, the Mg-calcite carbonates encountered at the CB mounds most likely indicate precipitation under reducing conditions within the sediment column³². Moreover, considering that dissolved sulfide can catalyze dehydration of Mg²⁺ ions in fluids, enabling more rapid precipitation of Mg-calcites as opposed to aragonite³³, the CB mounds may be sub-
<table><thead><tr><th rowspan="2">Sites</th><th rowspan="2">Sample names</th><th rowspan="2">Sample ID</th><th colspan="2">Location</th><th rowspan="2">Water depth (m)</th><th rowspan="2">Core depth (cm)</th><th colspan="2">Carbon and oxygen isotopes</th><th colspan="3">Mineral compositions</th><th rowspan="2">Reference</th></tr><tr><th>Longitude</th><th>Latitude</th><th>δ<sup>13</sup>C (VPDB)</th><th>δ<sup>18</sup>O (VPDB)</th><th>Arag. (wt %)</th><th>Cal. (wt %)</th><th>Dolo. (wt %)</th></tr></thead><tbody><tr><td rowspan="10">Chukchi borderland (Western Arctic)</td><td>ARA09C-ST07_GC03 45</td><td>CB-1</td><td>169° 47.6843' W</td><td>75° 42.7173' N</td><td>699</td><td>45</td><td>-35.7</td><td>5.2</td><td>n.d</td><td>25</td><td>n.d</td><td>This study</td></tr><tr><td>ARA09C-ST07_GC03 451</td><td>CB-2</td><td>169° 47.6843' W</td><td>75° 42.7173' N</td><td></td><td>451</td><td>-33.3</td><td>3.9</td><td>n.d</td><td>23</td><td>n.d</td><td>This study</td></tr><tr><td>ARA09C-ST14_GC03 71-74</td><td>CB-3</td><td>169° 44.3467' W</td><td>75° 40.3835' N</td><td>646</td><td>71</td><td>-37.0</td><td>5.2</td><td>n.d</td><td>24</td><td>3</td><td>This study</td></tr><tr><td>ARA09C-ST14_GC03 73</td><td>CB-4</td><td>169° 44.3467' W</td><td>75° 40.3835' N</td><td></td><td>73</td><td>-33.3</td><td>5.4</td><td>n.d</td><td>23</td><td>3</td><td>This study</td></tr><tr><td>ARA09C-ST04_GC04 75-80</td><td>CB-5</td><td>169° 44.2104' W</td><td>75° 40.7934' N</td><td>605</td><td>75</td><td>-34.8</td><td>5.6</td><td>n.d</td><td>32</td><td>6</td><td>This study</td></tr><tr><td>ARA09C-ST04_GC03 116-123</td><td>CB-6</td><td>169° 44.2104' W</td><td>75° 40.7934' N</td><td>605</td><td>116</td><td>-32.8</td><td>5.0</td><td>n.d</td><td>22</td><td>6</td><td>This study</td></tr><tr><td>ARA09C-ST06_GC01 38</td><td>CB-7</td><td>169° 44.1945' W</td><td>75° 40.8402' N</td><td>611</td><td>38</td><td>-35.1</td><td>5.5</td><td>n.d</td><td>46</td><td>n.d</td><td>This study</td></tr><tr><td>ZR2-PL13-P04</td><td>CF-1</td><td>-</td><td>-</td><td>2830</td><td>0</td><td>-58.5</td><td>5.3</td><td>n.d</td><td>&gt;50</td><td>n.d</td><td>Pierre and Fouquet<sup>25</sup></td></tr><tr><td>ZR2-PL14-P05</td><td>CF-2</td><td>-</td><td>-</td><td>3150</td><td>0</td><td>-53.6</td><td>5.6</td><td>n.d</td><td>&gt;50</td><td>n.d</td><td>Pierre and Fouquet<sup>25</sup></td></tr><tr><td>BZ1-GBT3-PL7-83</td><td>CF-3</td><td>-</td><td>-</td><td>3150</td><td>0</td><td>-49.5</td><td>2.9</td><td>n.d</td><td>&gt;50</td><td>n.d</td><td>Pierre and Fouquet<sup>25</sup></td></tr><tr><td rowspan="4">Nile deep-sea fan (NDSF)</td><td>NL4-CC1</td><td>NDSF-1</td><td>-</td><td>-</td><td>3032</td><td>0</td><td>-38.9</td><td>2.8</td><td>72</td><td>4</td><td>1</td><td>Gontharet et al.<sup>24</sup></td></tr><tr><td>NL7-CC1</td><td>NDSF-2</td><td>-</td><td>-</td><td>1691</td><td>0</td><td>-41.8</td><td>4.2</td><td>n.d</td><td>51</td><td>27</td><td>Gontharet et al.<sup>24</sup></td></tr><tr><td>NL14-CC5</td><td>NDSF-3</td><td>-</td><td>-</td><td>2129</td><td>0</td><td>-29.2</td><td>3.2</td><td>70</td><td>20</td><td>n.d</td><td>Gontharet et al.<sup>24</sup></td></tr><tr><td>NL20-CC1</td><td>NDSF-4</td><td>-</td><td>-</td><td>3018</td><td>0</td><td>-38.4</td><td>3.1</td><td>82</td><td>5</td><td>n.d</td><td>Gontharet et al.<sup>24</sup></td></tr><tr><td rowspan="4">Niger fan (NF)</td><td>N1-KS-07</td><td>NF-1</td><td>-</td><td>-</td><td>1633</td><td>0</td><td>-45.9</td><td>4.0</td><td>87</td><td>12</td><td>n.d</td><td>Rongemaille et al.<sup>27</sup></td></tr><tr><td>N1-KS-22</td><td>NF-2</td><td>-</td><td>-</td><td>1150</td><td>300</td><td>-27.8</td><td>4.9</td><td>90</td><td>10</td><td>n.d</td><td>Rongemaille et al.<sup>27</sup></td></tr><tr><td>N1-KSF-45</td><td>NF-3</td><td>-</td><td>-</td><td>1546</td><td>185</td><td>-47.1</td><td>6.1</td><td>n.d</td><td>6</td><td>87</td><td>Rongemaille et al.<sup>27</sup></td></tr><tr><td>N1-KI-47</td><td>NF-4</td><td>-</td><td>-</td><td>1540</td><td>12</td><td>-48.0</td><td>5.7</td><td>n.d</td><td>8</td><td>86</td><td>Rongemaille et al.<sup>27</sup></td></tr></tbody><tfoot><tr><td colspan="14">Continued</td></tr></tfoot></table>
<table><caption>**Table 1.** Site information and mineral and stable isotopic compositions of authigenic carbonate samples considered in this study. n.d. denotes 'not determined'. '—' indicates 'not available'.</caption><thead><tr><th rowspan="2">Sites</th><th rowspan="2">Sample names</th><th rowspan="2">Sample ID</th><th colspan="2">Location</th><th rowspan="2">Water depth (m)</th><th rowspan="2">Core depth (cm)</th><th colspan="2">Carbon and oxygen isotopes</th><th colspan="3">Mineral compositions</th><th rowspan="2">Reference</th></tr><tr><th>Longitude</th><th>Latitude</th><th>δ<sup>13</sup>C (VPDB)</th><th>δ<sup>18</sup>O (VPDB)</th><th>Arag. (wt %)</th><th>Cal. (wt %)</th><th>Dolo. (wt %)</th></tr></thead><tbody><tr><td rowspan="7">Gulf of Mexico (GoM)</td><td>#2 4173-2 (AT340)</td><td>GoM-1</td><td>-</td><td>-</td><td>2216</td><td>0</td><td>-54.2</td><td>4.0</td><td>63</td><td>n.d</td><td>n.d</td><td>Feng et al.<sup>10</sup>, Roberts et al.<sup>26</sup>, Feng and Roberts<sup>23</sup></td></tr><tr><td>#7 4174-2 (GC600)</td><td>GoM-2</td><td>-</td><td>-</td><td>1250</td><td>0</td><td>-21.4</td><td>4.8</td><td>n.d</td><td>52</td><td>n.d</td><td>Feng et al.<sup>10</sup>, Roberts et al.<sup>26</sup>, Feng and Roberts<sup>23</sup></td></tr><tr><td>#65 271-1 (MC462)</td><td>GoM-3</td><td>-</td><td>-</td><td>973</td><td>0</td><td>-40.5</td><td>4.6</td><td>80</td><td>n.d</td><td>n.d</td><td>Feng et al.<sup>10</sup>, Roberts et al.<sup>26</sup>, Feng and Roberts<sup>23</sup></td></tr><tr><td>#67 272-1 (GC415)</td><td>GoM-4</td><td>-</td><td>-</td><td>1110</td><td>0</td><td>-38.2</td><td>3.4</td><td>n.d</td><td>68</td><td>1</td><td>Feng et al.<sup>10</sup>, Roberts et al.<sup>26</sup>, Feng and Roberts<sup>23</sup></td></tr><tr><td>#68 273-1 (GC852)</td><td>GoM-5</td><td>-</td><td>-</td><td>1633</td><td>0</td><td>-40.2</td><td>3.7</td><td>n.d</td><td>70</td><td>1</td><td>Feng et al.<sup>10</sup>, Roberts et al.<sup>26</sup>, Feng and Roberts<sup>23</sup></td></tr><tr><td>#69 273-2 (GC852)</td><td>GoM-6</td><td>-</td><td>-</td><td>1633</td><td>0</td><td>-40.1</td><td>3.7</td><td>n.d</td><td>62</td><td>2</td><td>Feng et al.<sup>10</sup>, Roberts et al.<sup>26</sup>, Feng and Roberts<sup>23</sup></td></tr><tr><td>#70 273-3 (GC852)</td><td>GoM-7</td><td>-</td><td>-</td><td>1633</td><td>0</td><td>-48.1</td><td>4.5</td><td>n.d</td><td>54</td><td>2</td><td>Feng et al.<sup>10</sup>, Roberts et al.<sup>26</sup>, Feng and Roberts<sup>23</sup></td></tr></tbody></table>

jected to anoxic conditions over long periods of time, since sulfide-enriched microenvironments could have been shaped by AOM processes under anoxic conditions. Furthermore, preferential precipitation of calcite over aragonite may occur at relatively low temperatures³⁴, given that the temperature in the bottom seawaters of the Chukchi Sea is typically below 1 °C²¹. Thus, these properties could potentially explain the predominance of calcite over aragonite in the carbonate samples.

The measured δ¹³Ccarb values (−37.0 to −32.8‰; Table 1 and Fig. S2) were significantly lower than the typical values (~0‰) of marine bicarbonate (i.e., dissolved inorganic carbon in seawater)³⁵. In general, the δ¹³Ccarb values reflect the source of dissolved carbon incorporated during carbonate precipitation³⁶. The potential contributions of these different carbon isotopic signatures include (1) methane (δ¹³CCH₄ for biogenic origin; <−60‰, and δ¹³CCH₄ for thermogenic origin; −50 to −30‰) and (2) organic matter (ca. −25‰)³⁷,³⁸. Furthermore, the extent of ¹³C-depletion in carbonates can be an indicator of specific microbial processes³⁹,⁴⁰. Considering that the carbon isotopic composition of ascending methane (δ¹³CCH₄; −95.7 to −44.1‰) investigated at the CB mounds reflects a mixture of biogenic and thermogenic sources²², the microbial oxidation using this ascending methane can produce ¹³C-depleted bicarbonate (HCO₃⁻) in sediments. We could not determine the precise age of carbonate samples using the U-Th method, due to the large amounts of detrital Th in our samples. Nonetheless, considering the sampling depths of our carbonate samples (see Fig. S1), they might be weakly related to the current sulfate-methane transition zone (SMTZ). Thus, the isotopic depletion of all carbonates (δ¹³Ccarb < −30‰) at the CB mounds suggests that a substantial fraction of the carbon involved in carbonate precipitation was derived from the ascending methane-rich fluids, likely related to ancient methane releases similar to those that occurred in the Barents Sea⁴¹. Pore fluids become ¹⁸O-enriched (δ¹⁸Ofluid) during clay mineral dehydration or gas hydrate dissociation, which results in high δ¹⁸Ocarb values³⁰,⁴²,⁴³. At the studied site, gas hydrate layers have been identified at sediment depths below 3–4 m²². In this regard, the calcite-dominated carbonates (δ¹⁸Ocarb; 5.1 ± 0.6‰ VPDB (n = 7), Table 1 and Fig. S2) investigated at the CB mounds display ¹⁸O-enriched isotopic signatures. Therefore, these results may indicate that they precipitated from similarly enriched pore waters, resulting from gas hydrate dissociation.

**Trace element enrichments in authigenic carbonates.** Compared to marine bioskeletal carbonates (i.e., JLs-1 limestone), most trace elements in the carbonate fractions are enriched in the studied cold seep carbonates, except for Ti, Ba, and W (Fig. 2a). These enrichments are well explained as reflecting the geochemical composition of surrounding pore waters, inherited from various early diagenetic processes that typically result in
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Enrichment factors of (a) carbonate fraction (1 M AA leachates) data normalized to values for the JLs-1 carbonate reference material, (b) total lipid fraction data normalized to values for the JLs-1 carbonate reference material, and (c) total lipid fraction data normalized to values for the corresponding carbonate fractions (1 M AA leachates).</figcaption></figure>

selective trace element enrichment⁴⁴. In cold seeps, authigenic carbonates dominated by calcite typically incorporate higher amounts of trace elements (including REEs) compared to aragonite-rich samples, which form near the seafloor²⁷. Moreover, sustained anoxic conditions in cold seeps generally explain why many redox-sensitive trace elements (Mo, U, Ni, V, Co, and Zn) are enriched in corresponding authigenic carbonates, such as in the CB carbonate samples⁴⁵. Enrichment of Mo typically relates to Fe-Mn oxyhydroxide cycling in subsurface sediments and/or sulfide mineral formation, resulting from the presence of hydrogen sulfide in pore waters released from AOM⁴⁶. Similarly, the observed enrichments identified for Ni, Co, and Zn reflect the remineralization of organic matter in sub-surface sediments, which also results in enriched sedimentary pore waters⁴⁷,⁴⁸.

The total lipid extracts of CB carbonates displayed particular trace element enrichments (Li, Ti, Ni, Co, Cu, Zn, and Mo) compared to the corresponding carbonate fractions (Fig. 2c). In contrast to the enrichments observed in leached carbonate fractions, which are directly inherited from the chemical composition of the
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Lipid biomarker signatures of authigenic carbonates from the CB: (a) concentrations of major lipid biomarkers and ratios of *sn*-2-hydroxyarchaeol relative to archaeol, and (b) carbon isotopic compositions of major lipid biomarkers. Note that the dashed line indicates the average value of δ¹³C of methane (δ¹³C~CH₄~) while the shaded area shows the full range of δ¹³C~CH₄~ values.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Relative abundances of archaea taxa at the phylum and order level of Euryarchaeota in authigenic carbonates of the Chukchi Borderland (CB). ‘n.d.’ denotes ‘not determined’.</figcaption></figure>
pore waters from which they have precipitated, the enrichment associated with the lipid fractions hosted by authigenic carbonates can provide information on the bio-essential elements related to AOM-driven microbial activity at seeps²⁰. Thus far, Ni and Zn have been identified in microbial organism cells, especially in seep-related microbial assemblages, where they are bound to specific sites of proteins and enzymes⁴⁹,⁵⁰. Co is present in cobamides involved in methyl group transfer during methanogenesis and methanotrophic processes¹⁷. Mo is bound to a pterin cofactor to form molybdopterin, which catalyzes electron redox reactions⁵⁰, while Zn occurs as a single structural atom in both enzymes⁵⁰ (e.g., methanol-coenzyme M methyltransferase enzyme (Mta) and Heterodisulfide reductase (Hdr)). Considering that microbial metalloenzymes play a key role in catalyzing major biogeochemical reactions, in particular AOM⁵¹, these enzymes may require some trace elements (e.g., Co, Ni, Mo, and Zn) as co-factors for electron transport or as catalytic centers at active sites⁵⁰. Although trace elements (especially Cu and Ni) may be enriched in the organic (i.e., lipid) fractions due to their chemicophysical affinities to organic matter⁴⁵, trace elements, such as Co, Ni, Mo, and Zn, were enriched, not only in the lipid fractions (Fig. 2b) but also in the carbonate fractions (Fig. 2a). Hence, the enrichment of Co, Ni, Mo, and Zn identified in the total lipid fractions associated with CB carbonates is most likely to reflect their implication in specific metal-rich enzymatic pathways related to AOM-driven methanotrophic activity. Fe is primarily present as Fe-S clusters used for electron transport and/or catalysis, while Ni is either bound to Fe-S clusters or in the center of a porphyrin unique to methanogen/methanotrophs containing cofactor F430⁵². Thus, their enrichment in the total lipid extracts separated from CB carbonates could similarly reflect their potential role as a limiting factor for methanotrophic activities because of their high bioavailability as essential metalloenzymes during AOM⁵³.

Concerning Cu, which was also significantly enriched in the analyzed lipid fractions, it is worth mentioning that Cu-dependent anaerobic methanotrophic enzymes have not yet been identified, although Cu availability is important for aerobic methanotrophs⁵⁴. For instance, an increase in Cu concentration can cause up to a 55-fold expression of particulate methane monooxygenase (pMMO) as a membrane protein found in aerobic methanotrophs⁵⁰. Similar enrichments have been recently documented in cold seeps for light rare earth elements (e.g., La, Ce), also reflecting their preferential utility in aerobic methanotroph activity⁵⁵. In the CB carbonates, the total lipid fractions extracted did not show any particular enrichment in light REE, indicating a decoupling between Cu and light REE. Given that Mg-calcite carbonates were predominantly precipitated under anoxic conditions in the CB mounds, the observed Cu enrichments in the CB carbonates would not be associated with the aerobic methanotrophy. Hence, additional studies are required to better assess the potential utility of Cu in AOM treatment. Similarly, it still needs to be elucidated why Li and Ti were enriched in the studied lipid fractions in the CB carbonates, and whether the potential bioavailability of Li and Ti as metalloenzymes could represent a limiting factor in methanogenesis and methanotrophic activity under anoxic conditions. Notably, these elements have previously been identified in microorganism cells, especially in seep-related microbial assemblages, where they were known to be bound to specific sites of proteins and enzymes⁴⁹,⁵⁰. Thus, to some degree, we speculate that the observed Cu, Li, and Ti enrichments in the total lipids of CB carbonates might be indicative of alternative microbial processes for AOM as metalloenzymes. Hence, in many physiological and biogeochemical aspects, further research is required to investigate the potential availability of Cu, Li, and Ti in anaerobic methanogens/methanotrophs in pure cultures.

**Microbial signatures in authigenic carbonates.** In this study, irregular hydrocarbons (e.g., PMI) and isoprenoid DGDs (e.g., archaeol and *sn*-2-hydroxyarchaeol) were identified as AOM-related lipids, showing the predominance of isoprenoid DGDs in most CB carbonate samples (Fig. 3). The δ¹³C composition of these lipids indicated relatively ¹³C-depleted signatures, compared to those reported for methane in the CB mounds (Fig. 3). In general, methane-derived microbial lipids are depleted in ¹³C in comparison to the methane source as a result of isotopic fractionation during methane assimilation⁵⁶. The more ¹³C-depleted isoprenoid DGDs (i.e., archaeol and *sn*-2-hydroxyarchaeol) may indicate the occurrence of methanotrophs involved in AOM. Notably, CB-6 showed the most predominant ¹³C-depleted isoprenoid DGDs. This ¹³C-depleted pattern could be closely linked to active AOM taking place in sediment layers where the SMTZ currently occurs²² (90–120 cm sediment depths). On the other hand, because the depth of the SMTZ in methane seeps can fluctuate as a result of the upward methane fluxes⁵⁷,⁵⁸, low amounts of ¹³C-depleted lipids investigated in other CB carbonate samples could alternatively be explained as a result of the diffusive migration of methane-rich fluids⁵⁹. In this context, the low amount of ¹³C-depleted microbial lipids in the CB carbonates below and above the present SMTZ may potentially reflect the weak AOM activities as a fossil signature during the change in methane fluxes⁶⁰. Non-isoprenoid DGDs (If, Ila, and IId) derived from thermophilic and halophilic SRB⁶¹,⁶² were also detected in the studied CB carbonates (Fig. 3). Overall, the abundance of bacterial lipids in the CB carbonates was lower than that of AOM-related archaeal lipids. Considering previous studies in other cold seepages, such as the North Sea⁶³ and the Eastern Mediterranean Sea⁶⁴, the observed pattern between both lipid groups appear to reflect the process of AOM-coupled sulfate reduction. Furthermore, these bacterial lipids are more enriched in ¹³C compared to archaeal lipids (Δ¹³C = 10–15‰), which reflects their involvement in AOM as syntrophic partners⁷,⁶⁵. Based on these isotopic fractionations, the ¹³C-depleted non-isoprenoid DGDs in the CB carbonates (CB-1, -2, -6, and -7) may serve as evidence of syntrophic organisms involved in AOM. However, the δ¹³C values of the non-isoprenoid DGDs (–53.9 ± 3.6‰, n = 5) in other CB carbonates (CB-3, -4, and -5) were slightly enriched in ¹³C relative to the ascending methane reported in the CB mounds. Therefore, these ¹³C-enriched lipids may indicate that these compounds originate from a mixed source mediating AOM as well as other processes. Our results should be further confirmed by determining quantitative and qualitative properties of source organisms obtained from different analytical approaches, such as potassium hydroxide (KOH) hydrolysis.

Based on the microbial lipids and nucleic acids, the chemotaxonomy of key microbes involved in AOM can be determined⁶⁶. Three groups of ANME (ANME-1, ANME-2, and ANME-3) have often been reported in various
cold seep environments, which are related to methanogens from the order *Methanosarcinales* and *Methanomi-
crobiales*⁵⁸,⁶⁷. The microbial communities dominated by ANME-1 biosynthesize archaeol and PMI with relatively
low amounts of pentamethyllicosenes (PMEs), whereas ANME-2 is characterized by the predominance of *sn-
2*-hydroxyarchaeol (relative to archaeol) and abundant PMEs and crocetane/crocetenes⁶⁶. Together with the
presence of PMI in CB carbonates, the ratio of *sn*-2-hydroxyarchaeol to archaeol varied within the range 0.9–1.8
(Fig. 3), suggesting the dominant occurrence of ANME-1, which lived under anaerobically sulfate-depleted
conditions⁶⁸. These distinctive lipid patterns in most CB carbonates were also supported by higher propor-
tions of ASV032, which belongs to ANME-1 (Fig. 4). Furthermore, the co-occurrence of syntrophic bacterial
partners associated with ANME-1 can be indirectly inferred from the abundance of bacterial lipids, such as
non-isoprenoid DGDs and fatty acids⁶⁶. To date, ¹³C-depleted non-isoprenoid DGDs have been mainly found in
prominent ANME-1 systems⁶⁶. Thus, the less ¹³C-depleted non-isoprenoid DGDs in the CB carbonates appear
to be associated with an unknown SRB linked to AOM⁶⁹. Furthermore, the low abundance of branched fatty
acids (i.e., i-C15:0 and *ai*-C15:0) in the CB carbonates may also indicate a negligible presence of typical SRBs
(i.e., *Desulfosarcina/Desulfococcus*) associated with ANME-1. This inference may be partially supported by the
incomplete sequence of *Deltaproteobacteria*. Hence, these lipid characteristics in CB carbonates, in combination
with 16S rRNA gene pyrosequencing data, suggest that SRBs, not belonging to *Deltaproteobacteria* but possibly
to other bacterial communities, would have been involved in AOM in the CB mounds. Notably, ANME-1 has
been microscopically shown to occur as a single cell in cold seep environments⁷⁰. Accordingly, ANME-1-related
syntrophic relationships within the CB mounds appear to be more linked to other bacterial communities than
typical SRB consortia.

**Potential metalloenzymes involved in methanotrophic activity.** To further investigate the distri-
bution patterns of trace elements (W, Li, Ti, Ni, Co, Cu, Zn, and Mo) presumably acting as potential metal-
loenzymes involved in AOM at cold seeps, principal component analysis (PCA) was performed using the data
for total lipid extracts normalized to JLS-1 (Fig. 5a). We used these data to compare authigenic carbonate sam-
ples collected from different regions worldwide instead of those normalized to associated carbonate fractions
(Fig. S5). Although W was not detected in our total lipid samples of the CB, this element is included because pro-
nounced lipid-bound W anomalies are reported at other seep sites worldwide²⁰. PCA was also performed using
the archaeal and bacterial lipid data to compare the CB signatures with those from other cold seeps²⁰ (Fig. 5b–c).
The PCA results show that the combined distribution of trace elements and bacterial lipids in CB carbonates
differs from those reported for other methane seeps worldwide (e.g., Congo fan; CF, Nile deep-sea fan; NDSF,
Niger fan; NF, and Gulf of Mexico; GoM). This finding suggests that the enrichment of Ni, Co, Zn, Cu, Ti, and
Mo in carbonate-hosted lipid fractions could potentially be related to different methanotrophic communities
and their metalloenzymes. This hypothesis is plausible because ANME-1 is predominant in CB carbonates,
while ANME-1 and *Methanosarcinales* (potential ANME-2) co-exist at other sites²⁰ (Figs. 6, 7). Ni-containing
methyl-coenzyme MCR is the key enzyme in AOM¹⁹. Two different MCR homologues, designated as proteins
I and II⁷¹,⁷² have been isolated from cold seep areas: protein I contains (17S)-17²-methylthio-F430, whereas
protein II contains F430¹⁸,⁷³. Considering that protein I is predominantly detected in ANME-1⁷⁴, we may infer
that particular environmental conditions (e.g., temperature and partial pressure of methane) in the CB mound
area specifically promote the synthesis of protein I for AOM⁷⁵ compared to what can be observed at other seep
sites. Future studies targeting deep gas-fluid sources and subsurface lithology are needed to unravel the crucial
factors for determining discriminative ANME activities involved in the stimulated response under permafrost
dissociation. The impact of lipid/enzyme degradation on the correlations between AOM-related lipids and trace
elements should also be tested by comparing their relationships using data from modern and ancient carbonate
samples. Previous studies have shown that there is a significant trace element enrichment difference between
intra- and inter-crystalline lipid extractions in seep carbonates⁴⁸. Hence, future studies should also explore the
effect of the two-step extraction procedure on trace element distributions.

Conclusions

We analyzed a series of authigenic carbonate samples collected at methane seeps from the southwestern slope of the CB using various inorganic and organic proxies of fluid sources and biogeochemical processes (stable isotopes, mineralogy, trace elements, lipid biomarkers, and nucleic acids). The absence of crocetane, the abundance of PMI, and the relatively low ratio of *sn*-2-hydroxyarchaeol to archaeol in the CB carbonates indicated the dominant occurrence of ANME-1, as confirmed by 16S rRNA gene pyrosequencing. Specific trace elements (Ni, Co, Zn, Cu, Ti, and Mo) appeared to be enriched in the total lipid fractions extracted from CB carbonates, confirming recent findings of their utility as metalloenzymes in microbial processes related to AOM at cold seeps. The patterns of trace element enrichment and microbial lipid assemblages in CB carbonates reveal distinct specific characteristics compared to those from other seeps worldwide (e.g., Congo fan, Nile deep-sea fan, Niger fan, and Gulf of Mexico). This finding suggests that trace element enrichment in carbonate-hosted lipid fractions at cold seeps could vary depending on the type of AOM microbial assemblages among the ANME groups. Further work is required to elucidate the role of trace elements, such as Co, Zn, Cu, Ti, and Mo, as potential metalloenzymes in AOM during N₂ fixation or aerobic methane oxidation coupled with denitrification. In this regard, pure-culture experiments are necessary to understand their bioavailability for metalloenzymes involved in AOM.

**Material and methods**

**Sampling.** Previous geophysical surveys (e.g., multibeam and sub-bottom profiling) conducted at the south-
western slope of the CB²¹ led to the discovery of gas hydrate mound structures characterized by acoustic blank-
ing below the crest at certain depth intervals (Fig. 1). During the ARA09C expedition in 2018, sediment cores
<figure><img src="image_8.png" /></figure>

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

<figure><img src="image_10.png" /><figcaption>**Figure 5.** PCA results for (a) the trace metal data for the total lipids normalized to the values of the JLs-1 carbonate reference material, and for the concentration data of (b) archaeal lipids and (c) bacterial lipids.</figcaption></figure>

(ARA09C-ST07, -ST14, and -ST04/06) were retrieved from three mound structures, together with various authigenic carbonate samples recovered at different sediment depths between 38 and 415 cm (Fig. 1). The recovered carbonate samples were cleaned with ultrapure Milli-Q (MQ) water and air-dried prior to being crushed into powder using an agate and mortar.

**Mineral and stable isotope analysis.** The bulk authigenic carbonate samples were ground to a fine powder (<20 µm particle size) for mineral analysis using a Philips X-ray diffraction (XRD) system at the laboratories of the Crystallography research group (Central Laboratory for Crystallography and Applied Material Sciences, University of Bremen). The X-ray source was a Cu anode operated at 40 kV and 45 mA using CuKα radiation equipped with a diffracted beam graphite monochromator. Step scans were run from 5 to 65° (2θ) with a step size of 0.02° and a count time of 2 s per step (0.6°/min). Mineral identification was performed using the Philips X' Pert HighScore software, which provides a semi-quantitative value for each identified mineral on the basis of relative intensity ratio values. The oxygen and carbon isotopic compositions (δ¹⁸O and δ¹³C) of the carbonate samples were measured at Beta Analytic (service sample number: 536510 to 536,516). In brief, the sample was placed into a clean vial and dissolved in H₃PO₄ to release CO₂ into the headspace of the vial. The sample was placed in a temperature-controlled tray (72 °C), where CO₂ equilibration took place for one hour. Finally,
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Comparison of lipid biomarker signatures of authigenic carbonates from the CB, Congo fan (CF), Nile deep-sea fan (NDSF), Niger fan (NF), and Gulf of Mexico (GoM): (a) concentrations of major lipid biomarkers and ratios of *sn*-2-hydroxyarchaeol relative to archaeol, and (b) carbon isotopic compositions of major lipid biomarkers with those of carbonates.</figcaption></figure>

helium (He) was added to the sample and a mixture of He and CO₂ was injected into a Finnegan MAT 251 mass spectrometer. All isotopic values were expressed using the δ-notation relative to the Vienna Pee Dee Belemnite (VPDB) standard and were reported in per mil (‰) with a standard deviation of less than 0.3‰ for both δ¹⁸O and δ¹³C values.

**Trace element analysis.** Homogenized carbonate samples (~5 g) were treated using the sequential extraction procedure previously described²⁰ to selectively separate carbonate, sulfide, and detrital silicate fractions. We used the lipid fractions obtained from the solvent extraction method using a mixture of organic solvents (dichloromethane (DCM):methanol (MeOH), 2:1 v/v) to compare trace elements in lipid fractions with lipid biomarker data. The total lipid extracts (TLEs) were split into two aliquots for trace element analysis at IFREMER (Institut Français de Recherche pour l'Exploitation de la Mer; i.e., "French Research Institute for Exploitation of the Sea") and lipid biomarker analysis at the Korea Polar Research Institute (KOPRI). All analyses were performed using an Element XR™ inductively coupled plasma-mass spectrometer (ICP-MS) at the Pôle Spectrométrie Océan (Brest, France). Polyatomic oxide and hydroxide interferences on the rare earth element (REE) were corrected using oxide formation rates determined by analyzing solutions of MQ–H₂O, Ba + Ce, Pr + Nd, and Sm + Tb at the beginning of each measurement session and applied to all samples. Elemental concentrations were calculated using the Tm addition method⁷,⁷,⁷. The internal precision of all measurements was generally better than 5%. Repeated analyses of the JLs-1 reference material (Triassic limestone) were also performed, with a precision of < 10% for most trace elements, except for Li (11.8% relative standard deviation; RSD), Ti (61% RSD), and Zr (17% RSD). Note that the relatively poor RSD determined for Ti mostly reflects the very low abundance of this element in JLs-1 (4.1 µg/g). Significant Ti enrichments were measured in the lipid fractions, up to three orders of magnitude (Supplementary Tables); hence, much higher than the RSD determined for JLs-1. Therefore, we are confident that the measured Ti enrichments in the studied lipid extracts are meaningful, and thus, Ti data are still reported in this study. Because of the high Ba/REE ratios, several carbonate samples analyzed (including the carbonate reference material JLs-1) displayed anomalously high Eu (and to a lesser extent Gd) concentrations as a result of under-corrected interferences; hence, these two elements were not reported here.

**Lipid biomarker analysis.** Detailed procedures for lipid biomarker analyses have been previously described⁶⁰. In short, one-half of the TLE was dried over anhydrous Na₂SO₄ and treated with tetrabutylam-
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Phylogenetic tree of archaeal rDNA sequences obtained from authigenic carbonates.</figcaption></figure>

monium sulfite reagent to remove elemental sulfur. The TLE was chromatographically separated into apolar and polar fractions over an Al₂O₃ column (activated for 2 h at 150 °C). The apolar fraction was eluted using hexane:DCM (9:1), and 40 µL of 5α-androstane (10 µg mL⁻¹) was added as an internal standard. The polar fraction was recovered with DCM:MeOH (1:1 v/v) as an eluent and divided into two aliquots, to which either C₂₂, 7,16-diol (10 µg mL⁻¹) or C₁₉ nonadecanoic acid (10 µg mL⁻¹) was added as an internal standard. Each aliquot was derivatized through silylation and methylation, prior to quantification by gas chromatography (GC) and identification by gas chromatography-mass spectrometry (GC-MS). GC and GC-MS conditions were as previously described⁶⁰. Molecular compounds were determined by comparing their mass spectral fragmentation patterns and retention times with previously published data⁶⁰. The δ¹³C values of lipid compounds were determined using an isotope ratio mass spectrometer (IRMS) connected to a GC via a combustion interface (glass tube packed with copper oxide, CuO, operated at 850 °C), as previously described⁶⁰. Isotopic values were expressed as δ¹³C values per mil relative to VPDB. The δ¹³C values were corrected for the introduction of additional carbon atoms during silylation and methylation. The analytical errors were less than ±0.4‰ for all lipid compounds.

**Nucleic acid analysis.** Genomic DNA was extracted using the FastDNA Spin Kit for Soil (MP Biomedicals, USA). DNA samples were subjected to polymerase chain reaction (PCR) amplification, library preparation, and paired-end Illumina MiSeq sequencing (2×300 bp) at the Integrated Microbiome Resource (IMR), Dalhousie University, Canada (http://cgeb-imr.ca). The primer pairs 515F/ 926R and 956F/ 1401R targeting the V4–V5 and V6–V8 regions were used to amplify bacterial and archaeal 16S rRNA genes, respectively⁷⁸,⁷⁹. Amplicons were sequenced using the paired-end (2×300 bp) Illumina MiSeq system (Illumina, USA) at the Integrated Microbiome Resource (IMR), Dalhousie University, Canada (http://cgeb-imr.ca). Samples that were not amplified were excluded from further analyses. The adapter and primer sequences were removed using Cutadapt version 2.10⁸⁰, and the resultant sequences were processed using DADA2 version 0.9.5⁸¹ to infer amplicon sequence variants (ASVs). For the quality trimming process, we applied a filtering option of maxN=0, maxEE=c (2,2), and truncQ=2. The low-quality tails of both reads were removed with truncLen=c (270,210). Denoising was performed after trimming based on the DADA2 error model. Sequences were dereplicated, and the core sample inference algorithm was applied to the dereplicated data. Paired reads were merged, and the chi-
meric sequences were removed. The following processes were performed after constructing a sequence table of ASVs to assign taxonomy using the software mothur⁸². Taxonomic assignments of representative ASV sequences were determined against the EzBiocloud database using sequence similarity searches⁸³. After taxonomic assignment, archaea and unknown ASVs were removed for bacterial analysis, and bacteria and unknown ASVs were removed for archaeal analysis. All sequence data used in this study were deposited in the Sequence Read Archive at the National Center for Biotechnology Information under accession number PRJNA825649. Phylogenetic trees of major archaeal ASVs of Methanomicrobia with relative abundances greater than 3% were constructed using the maximum-likelihood algorithm⁸⁴ with the general time reversible evolutionary model using MEGA X⁸⁵. The robustness of the tree topologies was assessed using bootstrap analyses based on 1,000 replications of the sequences.

**Statistical analysis.** To provide a general view of the variability of trace element abundances and microbial lipid distributions in the carbonate samples, gaps in the data set were first filled as previously described⁸⁶. When some components were not determined, a value of one-half of the minimum value detected for that variable in the whole data set was set as the limit of detection. These values were then replaced by a random number between zero and the limit of detection. Finally, the fractional abundances of trace elements and microbial lipids were obtained by normalizing each concentration to the summed concentrations of all trace elements and microbial lipids, respectively. Based on the fractional abundances, PCA was performed using R software version 3.4.2 (package information; FactoMineR, Vienna, Austria).

**Data availability**

All DNA sequence data generated for this study were deposited in the Sequence Read Archive at the National Center for Biotechnology Information under accession number PRJNA825649.

Received: 29 March 2022; Accepted: 23 September 2022

Published online: 05 October 2022

**References**

1. Campbell, K. A. Hydrocarbon seep and hydrothermal vent paleoenvironments and paleontology: Past developments and future research directions. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **232**, 362–407 (2006).

2. Suess, E. Marine cold seeps and their manifestations: Geological control, biogeochemical criteria and environmental conditions. *Int. J. Earth Sci.* **103**, 1889–1916 (2014).

3. Judd, A. & Hovland, M. *Seabed Fluid Flow* (Cambridge University Press, 2007).

4. Torres, M. E. et al. Fluid and chemical flux in and out of sediments hosting methane hydrate deposits on Hydrate Ridge, OR, I: Hydrological processes. *Earth Planet. Sci. Lett.* **201**, 525–540 (2002).

5. Gründger, F. et al. Methane-fuelled biofilms predominantly composed of methanotrophic ANME-1 in Arctic gas hydrate-related sediments. *Sci. Rep.* **9**, 9725 (2019).

6. Solomon, E. A., Kastner, M., Jannasch, H., Robertson, G. & Weinstein, Y. Dynamic fluid flow and chemical fluxes associated with a seafloor gas hydrate deposit on the northern gulf of mexico slope. *Earth Planet. Sci. Lett.* **270**, 95–105 (2008).

7. Hinrichs, K. U. & Boetius, A. The anaerobic oxidation of methane: new insights in microbial ecology and biogeochemistry. in *Ocean Margin Systems* 457–477 (Berlin, Germany: Springer-Verlag, 2002).

8. Reeburgh, W. Oceanic methane biogeochemistry. *Am. Chem. Soc.* **107**, 486–513 (2007).

9. Ussler, W. & Paull, C. K. Rates of anaerobic oxidation of methane and authigenic carbonate mineralization in methane-rich deep-sea sediments inferred from models and geochemical profiles. *Earth Planet. Sci. Lett.* **266**, 271–287 (2008).

10. Feng, D. & Chen, D. Authigenic carbonates from methane seeps of the northern Congo fan: Microbial formation mechanism. *Mar. Pet. Geol.* **27**, 748–756 (2010).

11. Mazzini, A., Ivanov, M. K., Parnell, J., Stadnitskaia, A. & Cronin, B. T. Methane-related authigenic carbonates from the black Sea: Geochemical characterisation and relation to seeping fluids. *Mar. Geol.* **212**, 153–181 (2004).

12. Bayon, G. et al. Multi-disciplinary investigation of fluid seepage on an unstable margin: The case of the Central Nile deep sea fan. *Mar. Geol.* **261**, 92–104 (2009).

13. Yao, H., Panieri, G., Lehmann, M. F., Himmler, T. & Niemann, H. Biomarker and isotopic composition of seep carbonates record environmental conditions in two arctic methane seeps. *Front. Earth Sci.* **8**, 1–12 (2021).

14. Boetius, A. et al. A marine microbial consortium apparently mediating anaerobic oxidation of methane. *Nature* **407**, 623–626 (2000).

15. Hinrichs, K. U., Hayes, J. M., Sylva, S. P., Brewer, P. G. & DeLong, E. F. Methane-consuming archaebacteria in marine sediments. *Nature* **398**, 802–805 (1999).

16. Beal, E. J., House, C. H. & Orphan, V. J. Manganese- and iron-dependent marine methane oxidation. *Science (80-)* **325**, 184–187 (2009).

17. Glass, J. B. et al. Geochemical, metagenomic and metaproteomic insights into trace metal utilization by methane-oxidizing microbial consortia in sulphidic marine sediments. *Environ. Microbiol.* **16**, 1592–1611 (2014).

18. Krüger, M. et al. A conspicuous nickel protein in microbial mats that oxidize methane anaerobically. *Nature* **426**, 878–881 (2003).

19. Scheller, S., Goenrich, M., Boecher, R., Thauer, R. K. & Jaun, B. The key nickel enzyme of methanogenesis catalyses the anaerobic oxidation of methane. *Nature* **465**, 606–608 (2010).

20. Wang, X. et al. Trace element systematics in cold seep carbonates and associated lipid compounds. *Chem. Geol.* **528**, 119277 (2019).

21. Kim, Y. G. et al. Occurrence of active gas hydrate mounds in the southwestern slope of the Chukchi plateau. *Arctic Ocean. Epis.* **43**, 811–823 (2020).

22. Kim, J. H. et al. Upwarding gas source and postgenetic processes in the shallow sediments from the ARAON mounds, Chukchi sea. *J. Nat. Gas Sci. Eng.* **76**, 103223 (2020).

23. Feng, D. & Roberts, H. H. Geochemical characteristics of the barite deposits at cold seeps from the northern Gulf of Mexico continental slope. *Earth Planet. Sci. Lett.* **309**, 89–99 (2011).

24. Gontharet, S. et al. Nature and origin of diagenetic carbonate crusts and concretions from mud volcanoes and pockmarks of the nile deep-sea fan (eastern Mediterranean Sea). *Deep. Res. Part II Top. Stud. Oceanogr.* **54**, 1292–1311 (2007).

25. Pierre, C. & Fouquet, Y. Authigenic carbonates from methane seeps of the Congo deep-sea fan. *Geo-Marine Lett.* **27**, 249–257 (2007).
26. Roberts, H. H., Feng, D. & Joye, S. B. Cold-seep carbonates of the middle and lower continental slope, northern Gulf of Mexico. *Deep. Res. Part II Top. Stud. Oceanogr.* **57**, 2040–2054 (2010).

27. Rongemaille, E. et al. Rare earth elements in cold seep carbonates from the Niger delta. *Chem. Geol.* **286**, 196–206 (2011).

28. Taylor, S. R. & McIennan, S. M. *The Continental Crust: Its Composition and Evolution. An Examination of the Geochemical Record Preserved in Sedimentary Rocks*. (Blackwell Scientific Publications, 1985).

29. Peckmann, J., Goedert, J. L., Thiel, V., Michaelis, W. & Reitner, J. A comprehensive approach to the study of methane-seep deposits from the lincoln creek formation, western Washington State, USA. *Sedimentology* **49**, 855–873 (2002).

30. Aloisi, G. et al. CH4-consuming microorganisms and the formation of carbonate crusts at cold seeps. *Earth Planet. Sci. Lett.* **203**, 195–203 (2002).

31. Campbell, K. A., Farmer, J. D. & Des Marais, D. Ancient hydrocarbon seeps from the Mesozoic convergent margin of California: Carbonate geochemistry, fluids and paleoenvironments. *Geofluids* **2**, 63–94 (2002).

32. Baker, P. A. & Kastner, M. Constraints on the formation of sedimentary dolomite. *Science*(80-) **213**, 214–216 (1981).

33. Smrzka, D. et al. Trace element distribution in methane-seep carbonates: The role of mineralogy and dissolved sulfide. *Chem. Geol.* **580**, 120357 (2021).

34. Burton, E. A. & Walter, L. M. Relative precipitation rates of aragonite and Mg calcite from seawater: Temperature or carbonate ion control?. *Geology* **15**, 111–114 (1987).

35. Aharon, P., Schwarcz, H. P. & Roberts, H. H. Radiometric dating of submarine hydrocarbon seeps in the Gulf of Mexico. *Geol. Soc. Am. Bull.* **109**, 568–579 (1997).

36. Pierre, C. & Rouchy, J. M. Carbonate replacements after sulfate evaporites in the middle Miocene of Egypt. *J. Sediment. Res.* **58**, 446–456 (1988).

37. Killops, S. D. & Killops, V. J. *Introduction to Organic Geochemistry* (Blackwell publishing, 2005).

38. Whiticar, M. J. Carbon and hydrogen isotope systematics of bacterial formation and oxidation of methane. *Chem. Geol.* **161**, 291–314 (1999).

39. Chen, Z. et al. Discovery of seep carbonate nodules as new evidence for gas venting on the northern continental slope of South China Sea. *Chinese Sci. Bull.* **51**, 1228–1237 (2006).

40. Ritger, S., Carson, B. & Suess, E. Methane-derived authigenic carbonates formed by subduction-induced pore-water expulsion along the Oregon/Washington margin: Reply. *Bull. Geol. Soc. Am.* **98**, 147–156 (1987).

41. Andreassen, K. et al. Massive blow-out craters formed by hydrate-controlled methane expulsion from the Arctic seafloor. *Science*(80-) **356**, 948–953 (2017).

42. Feng, D. et al. Time integrated variation of sources of fluids and seepage dynamics archived in authigenic carbonates from Gulf of Mexico Gas Hydrate Sea floor observatory. *Chem. Geol.* **385**, 129–139 (2014).

43. Naehr, T. H. et al. Authigenic carbonate formation at hydrocarbon seeps in continental margin sediments: A comparative study. *Deep. Res. Part II Top. Stud. Oceanogr.* **54**, 1268–1291 (2007).

44. Haley, B. A., Klinkhammer, G. P. & McManus, J. Rare earth elements in pore waters of marine sediments. *Geochim. Cosmochim. Acta* **68**, 1265–1279 (2004).

45. Smrzka, D. et al. *The behavior of trace elements in seawater, sedimentary pore water, and their incorporation into carbonate minerals: a review*. *Facies* vol. 65 (Springer Berlin Heidelberg, 2019).

46. Peketi, A. et al. Tracing the Paleo sulfate-methane transition zones and H2S seepage events in marine sediments: An application of C-S-Mo systematics. *Geochem. Geophys. Geosystems* **13**, 1–11 (2012).

47. Algeo, T. J. & Maynard, J. B. Trace-element behavior and redox facies in core shales of upper pennsylvanian kansas-type cyclothems. *Chem. Geol.* **206**, 289–318 (2004).

48. Smrzka, D. et al. Trace elements in methane-seep carbonates: Potentials, limitations, and perspectives. *Earth-Sci. Rev.* **208**, 103263 (2020).

49. Barton, L. L. et al. The bacterial metallome: Composition and stability with specific reference to the anaerobic bacterium desulfovibrio desulfuricans. *Biometals* **20**, 291–302 (2007).

50. Glass, J. B. & Orphan, V. J. Trace metal requirements for microbial enzymes involved in the production and consumption of methane and nitrous oxide. *Front. Microbiol.* **3**, 1–20 (2012).

51. Zerkle, A. L., House, C. H. & Brantley, S. L. Biogeochemical signatures through time as inferred from whole microbial genomes. *Am. J. Sci.* **305**, 467–502 (2005).

52. Ferry, J. G. *Methanogenesis: Ecology, physiology, biochemistry & genetics* (Springer Science & Business Media, 2012).

53. Glass, J. B. et al. Trace metal imaging of sulfate-reducing bacteria and methanogenic archaea at single-cell resolution by synchrotron X-Ray fluorescence imaging. *Geomicrobiol. J.* **35**, 81–89 (2018).

54. Sushkevich, V. L., Palagin, D., Ranocchiari, M. & Van Bokhoven, J. A. Selective anaerobic oxidation of methane enables direct synthesis of methanol. *Science*(80-) **356**, 523–527 (2017).

55. Bayon, G. et al. Microbial utilization of rare earth elements at cold seeps related to aerobic methane oxidation. *Chem. Geol.* **555**, 119832 (2020).

56. Niemann, H. & Elvert, M. Diagnostic lipid biomarker and stable carbon isotope signatures of microbial communities mediating the anaerobic oxidation of methane with sulphate. *Org. Geochem.* **39**, 1668–1677 (2008).

57. Lee, D.-H. et al. Chemosynthetic bacterial signatures in frenulata tubeworm (Oligobrachia sp.) in an active mud volcano of the Canadian Beaufort Sea. *Mar. Ecol. Prog. Ser.* **628**, 95–104 (2019).

58. Niemann, H. et al. Novel microbial communities of the Häkon Mosby mud volcano and their role as a methane sink. *Nature* **443**, 854–858 (2006).

59. Peckmann, J. & Thiel, V. Carbon cycling at ancient methane-seeps. *Chem. Geol.* **205**, 443–467 (2004).

60. Lee, D.-H. et al. Biogeochemical evidence of anaerobic methane oxidation on active submarine mud volcanoes on the continental slope of the Canadian Beaufort Sea. *Biogeosciences* **15**, 7419–7433 (2018).

61. Langworthy, T. A., Holzer, G., Zeikus, J. G. & Tornabene, T. G. Iso-and anteiso-branched glycerol diethers of the thermophilic anaerobe Thermodesulfobacterium commune. *Syst. Appl. Microbiol.* **4**, 1–17 (1983).

62. Sturt, H. F., Summons, R. E., Smith, K., Elvert, M. & Hinrichs, K.-U. Intact polar membrane lipids in prokaryotes and sediments deciphered by high-performance liquid chromatography/electrospray ionization multistage mass spectrometry—new biomarkers for biogeochemistry and microbial ecology. *Rapid Commun. Mass Spectrom.* **18**, 617–628 (2004).

63. Niemann, H. et al. Methane emission and consumption at a north sea gas seep (Tommeliten area) methane emission and consumption at a North Sea gas seep (Tommeliten area). *Biogeosciences* **2**, 1197–1241 (2005).

64. Gontharet, S., Stadnitskaia, A., Bouloubassi, I., Pierre, C. & Damsté, J. S. S. Palaeo methane-seepage history traced by biomarker patterns in a carbonate crust, Nile deep-sea fan (Eastern mediterranean Sea). *Mar. Geol.* **261**, 105–113 (2009).

65. Wegener, G., Niemann, H., Elvert, M., Hinrichs, K. & Boetius, A. Assimilation of methane and inorganic carbon by microbial communities mediating the anaerobic oxidation of methane. *Environ. Microbiol.* **10**, 2287–2298 (2008).

66. Blumenberg, M., Seifert, R., Reitner, J., Pape, T. & Michaelis, W. Membrane lipid patterns typify distinct anaerobic methanotrophic consortia. *Proc. Natl. Acad. Sci.* **101**, 11111–11116 (2004).

67. Knittel, K. & Boetius, A. Anaerobic oxidation of methane: progress with an unknown process. *Annu. Rev. Microbiol.* **63**, 311–334 (2009).
68. Lazar, C. S., Dinasquet, J., L'Haridon, S., Pignet, P. & Toffin, L. Distribution of anaerobic methane-oxidizing and sulfate-reducing communities in the G11 Nyegga pockmark, Norwegian Sea. *Antonie van Leeuwenhoek. Int. J. Gen. Mol. Microbiol.* **100**, 639–653 (2011).

69. Pancost, R. D., Hopmans, E. C., Sinninghe Damsté, J. S., The M. S. Scientific Party. Archaeal lipids in mediterranean cold seeps: Molecular proxies for anaerobic methane oxidation. *Geochim Cosmochim Acta* **65**, 1611–1627 (2001).

70. Orphan, V. J., House, C. H., Hinrichs, K.-U., McKeegan, K. D. & DeLong, E. F. Multiple archaeal groups mediate methane oxidation in anoxic cold seep sediments. *Proc. Natl. Acad. Sci.* **99**, 7663–7668 (2002).

71. Nauhaus, K., Albrecht, M., Elvert, M., Boetius, A. & Widdel, F. In vitro cell growth of marine archaeal-bacterial consortia during anaerobic oxidation of methane with sulfate. *Environ. Microbiol.* **9**, 187–196 (2007).

72. Thauer, R. K., Kaster, A. K., Seedorf, H., Buckel, W. & Hedderich, R. Methanogenic archaea: Ecologically relevant differences in energy conservation. *Nat. Rev. Microbiol.* **6**, 579–591 (2008).

73. Mayr, S. et al. Structure of an F430 variant from archaea associated with anaerobic oxidation of methane. *J. Am. Chem. Soc.* **130**, 10758–10767 (2008).

74. Kaneko, M. et al. Quantitative analysis of coenzyme F430 in environmental samples: A new diagnostic tool for methanogenesis and anaerobic methane oxidation. *Anal. Chem.* **86**, 3633–3638 (2014).

75. Knittel, K. et al. Diversity and distribution of methanotrophic archaea at cold seeps. *Appl. Environ. Microbiol.* **71**, 467–479 (2005).

76. Barrat, J. A. et al. Determination of rare earth elements in sixteen silicate reference samples by ICP-MS after TM addition and ion exchange separation. *Geostand. Newsl.* **20**, 133–139 (1996).

77. Bayon, G. et al. Multi-disciplinary investigation of fluid seepage on an unstable margin: The case of the central Nile deep sea fan. *Mar. Geol.* **261**, 92–104 (2009).

78. Comeau, A. M., Li, W. K. W., Tremblay, J. É., Carmack, E. C. & Lovejoy, C. Arctic ocean microbial community structure before and after the 2007 record sea ice minimum. *PLoS One* **6**(11), e27492 (2011).

79. Walters, W. et al. Improved bacterial 16s rrna gene (v4 and v4–5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. *MSystems* **1**(1), e00009–e00015 (2016).

80. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. *EMBNet. J.* **17**(1), 10–12 (2011).

81. Callahan, B. J. et al. DADA2: High-resolution sample inference from illumina amplicon data. *Nat. Methods* **13**, 581–583 (2016).

82. Schloss, P. D. et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. *Appl. Environ. Microbiol.* **75**, 7537–7541 (2009).

83. Yoon, S. H. et al. Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. *Int. J. Syst. Evol. Microbiol.* **67**, 1613–1617 (2017).

84. Felsenstein, J. Evolutionary trees from DNA sequences: A maximum likelihood approach. *J. Mol. Evol.* **17**, 368–376 (1981).

85. Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. *Mol. Biol. Evol.* **35**, 1547–1549 (2018).

86. Yunker, M. B., Belicka, L. L., Harvey, H. R. & Macdonald, R. W. Tracing the inputs and fate of marine and terrigenous organic matter in Arctic ocean sediments: A multivariate analysis of lipid biomarkers. *Deep. Res. Part II Top. Stud. Oceanogr.* **52**, 3478–3508 (2005).

## Acknowledgements

We would like to thank the captain and crew of the R/V ARAON for their help at sea. This study was supported by the Korea Institute of Marine Science and Technology Promotion (KIMST) (KIMST Grant 20210632; KOPRI-PM22050) and a National Research Foundation of Korea (NRF) grant funded by the Ministry of Science and ICT (NRF-2021M1A5A1075512; KOPRI-PN22013).

## Author contributions

D.H.L. and J.H.K. designed and coordinated the study. D.H.L., J.H.K., Y.M.L., Y.J.J., and Y.K.J. conducted the field work during the ARAON expedition. G.B. and X.W. conducted the trace element analysis. D.H.L. and D.K. carried out the lipid biomarker analysis. Y.M.L. conducted the nucleic acid analysis. D.H.L. and J.H.K. led the writing and interpretation of the data. All authors commented on the manuscript and contributed to the writing.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-21184-6.

**Correspondence** and requests for materials should be addressed to J.-H.K.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022