Impact of metabolism and temperature on 2 H / 1 H fractionation in lipids of the marine bacterium Shewanella piezotolerans WP3

. Compound-speciﬁc hydrogen isotopes have increasingly been used as a powerful proxy for investigating biogeochemical cycles and climate change over the past 2 decades. Understanding the hydrogen isotope in extant organisms is fundamental for us to interpret such isotope signals preserved in natural environmental samples. Here, we studied the controls on hydrogen isotope fractionation between fatty acids and growth water by an Fe-reducing heterotrophic marine bacterium Shewanella piezo-tolerans WP3 growing on different organic substrates, including N-acetyl-D-glucosamine (GlcNac), glucose, acetate, pyruvate, L-alanine, and L-glutamate. Meanwhile, we also evaluated the impact of growth temperature on the hydrogen isotope composition of fatty acids using GlcNac as the sole organic substrate. Our results show that the abundance-weighted mean fatty-acid


Introduction
The hydrogen isotope composition of lipids preserved in sediments is widely applied to trace biogeochemical processes and study regional paleo-hydrological cycles (Hayes, 2001;Huang et al., 2004;Sachse et al., 2012;Sessions, 2016).During the past 2 decades, compound-specific hydrogen isotopes have increasingly been utilized as a valid proxy for reconstructing precipitation δ 2 H because relatively constant 2 H/ 1 H fractionations in terrestrial plants and phototrophic algae are assumed (Sachse et al., 2012;Sessions, 2016;Huang and Meyers, 2019).However, increasing numbers of studies have found that there are large ranges of hydrogen isotope ratios in lipids (up to 700 ‰) from various environmental samples, including microbial mats from hot springs (Naraoka et al., 2010;Osburn et al., 2011), marine particulate organic matter (Jones et al., 2008;Heinzelmann et al., 2016), X. Chen et al.: Impact of metabolism and temperature on 2 H / 1 H fractionation and marine and lacustrine sediments (Li et al., 2009;Chen et al., 2021).This raises questions as to what factors cause such isotope signals in natural environmental samples.Clues are found in a breakthrough culture study by X. Zhang et al. (2009), who reported that hydrogen isotope compositions of lipids are correlated with microbial metabolisms, with strong depletion in 2 H for chemoautotrophs but slight depletion for phototrophs and enrichment for heterotrophs.These features provide potential insights for interpreting hydrogen isotope signals in natural environments, and we can further infer that the metabolism of microbes thrives in modern and ancient environments.
Among these microbes, heterotrophs are of particular concern due to their essential role in the organic carbon cycle on earth (Moran et al., 2022).Understanding the hydrogen isotope composition of lipids in heterotrophs with different central metabolic pathways is important for decoding such isotope signals in geologic archives.Recently, in a systematic laboratory culture, Wijker et al. (2019) showed that aerobic heterotrophic bacteria display the largest range of variations, with ε lipids/water values ranging from ∼ −150 ‰ to +400 ‰, which are quantitatively related to the central metabolic pathways.In particular, heterotrophic bacteria exhibit the strongest 2 H enrichment growing on tricarboxylic acid (TCA) cycle substrates (e.g., acetate, succinate and citrate), but they show 2 H depletion on sugars (Wijker et al., 2019).The most plausible explanations of hydrogen isotopic differences are related to central metabolic pathways in the production of nicotinamide adenine dinucleotide phosphate (NADPH), which is the major source of H (> 50 %) in fatty acids (Jackson et al., 1999;X. Zhang et al., 2009;Heinzelmann et al., 2015a, b;Wijker et al., 2019).Combined with laboratory experiments and metabolic flux model analysis, ε lipids/water values show positive correlations with the relative flux (% of glucose uptake) of Entner-Doudoroff (ED) and TCA cycle pathways, while they show negative correlations with Embden-Meyerhof-Parnas (EMP) and pentose phosphate (PP) pathways (Wijker et al., 2019).However, the most investigated species by previous studies are restricted to model bacteria (e.g., Escherichia coli and Rhizobium radiobacter) and are unlikely to be widely distributed in natural environments.More importantly, understanding the hydrogen isotopes generated in a bacterium with wide distributions and ecological functions is essential for providing further constraints on the isotope signals in natural environments.
Compared with metabolisms, several other factors such as the lipid biosynthetic pathway, growth phase and rate, salinity, and pressure have been thought to potentially influence the hydrogen isotope composition of lipids (Heinzelmann et al., 2015a, b;Zhao et al., 2020;Dirghangi and Pagani, 2013a).For example, long-chain polyunsaturated fatty acids biosynthesized through the polyketide pathway have much lower δ 2 H values than those of short-chain lengths via the fatty-acid synthase (FAS) pathway (Fang et al., 2014).The ε lipids/water values of the C 16 fatty acid in Pseudomonas strain LFY10 are increased by only about 10 ‰-15 ‰ from exponential to stationary phases (Heinzelmann et al., 2015b).For environmental factors, Heinzelmann et al. (2015a) showed that the hydrogen isotope ratios of fatty acids exhibited no obvious correlations with salinity in cultures of heterotrophic microbe Pseudomonas strain LFY10, while they were enriched in 2 H with increasing salinity in the phototrophic algae Isochrysis galbana.Hydrogen isotope fractionations exhibited an inverse correlation with growth pressure in a deepsea gram-positive piezotolerant bacterium Sporosarcina sp.DSK25 (Zhao et al., 2020).Still, the variations associated with these factors are rather small (∼ 10 ‰-20 ‰) compared with metabolisms.Temperature is one of the most important environmental factors, and culture studies have demonstrated its impacts on the hydrogen isotope composition of lipids in phototrophs and archaea (Schouten et al., 2006;Dirghangi and Pagani, 2013b).Until now, however, there have been no cultures focusing on whether or not temperature can affect lipid hydrogen isotopes in heterotrophic bacteria and, if so, by what relationship between fractionations and growth temperatures.
The Shewanella genus is ubiquitously distributed around the globe, such as in surface freshwater and the deepest marine trenches, and can use various electron acceptors for metabolic activities (Hau and Gralnick, 2007;Fredrickson et al., 2008).Shewanella piezotolerans WP3 is a gram-negative, moderately halophilic bacterium within Shewanella genus and has an important role in the biochemical cycle of organic matter in the deep sea (Xiao et al., 2007;Lemaire et al., 2020).In our previous studies, the metabolisms of S. piezotolerans WP3 have been systematically studied through pure cultures and genomics (Wang et al., 2008;Dufault-Thompson et al., 2017).Specifically, S. piezotolerans WP3 can use a variety of organic substrates such as GlcNac, acetate, glucose, pyruvate, and amino acid as sole carbon sources (Xiao et al., 2007).Furthermore, it has relatively wide ranges of growth temperatures from 0-28 • C, with an optimal temperature at 15-20 • C (Wang et al., 2004).Here, we investigate the compound-specific hydrogen isotopes in S. piezotolerans WP3 growing on different substrates, and also discuss the effect of temperatures on hydrogen isotope fractionations using GlcNac as the sole carbon source.Then, we compare the fractionations with a metabolic flux model to decode the detail fractionation mechanisms in heterotrophic bacteria.Our study can provide the characterization of hydrogen isotope fractionations in S. piezotolerans WP3 lipids and may give insights into the biochemical cycle and ecological functions in Fe-reducing oceanic sediments.
2 Materials and methods 2.1 Strains, media, and growth conditions S. piezotolerans WP3 was first isolated from west Pacific deep-ocean sediment at a depth of 1914 m (Xiao et al., 2007).S. piezotolerans WP3 cultures have been reported by our previous study (Chen et al., 2022).Briefly, S. piezotolerans WP3 was grown in batch culture in a growth medium containing the following reagents (g L −1 ): NaCl (26.0), MgCl 2 • 6H 2 O (5.0), CaCl 2 • 2H 2 O (1.4), Na 2 SO 4 (4.0),NH 4 Cl (1.5), KH 2 PO 4 (0.1), and KCl (0.5) as well as trace minerals and vitamins.To examine the hydrogen fractionations of metabolic type, GlcNac, glucose, sodium acetate, sodium pyruvate, and amino acids were supplemented to each culture as the sole carbon source, with a final concentration of 4 g L −1 and culturing at 20 • C and shaking at 200 rpm.WP3 was grown at 4, 10, 15, 20, and 25 • C using GlcNac as the sole carbon source under static conditions for investigating the relationship between growth temperature and hydrogen isotope fractionation.All media were passed through a 0.22 µm filter.Final pH was adjusted to about 7.0.For water isotope analysis, 2 mL samples were taken from each culture before incubation.Duplicate cultures were maintained for each condition, and cells were harvested in early stationary phase.Samples of 20 mL were taken, washed with culture salt solutions, and centrifuged for 20 min at 4500 g; cell pellets were stored at −20 • C before analysis.

Extraction of cellular lipids
Cellular lipid extraction and analysis follows the procedure described in Rodríguez-Ruiz et al. (1998).Briefly, about 1-20 mg freeze-dried biomass was transesterified using 1 mL hexane and 2 mL 20 : 1 anhydrous methanol : acetyl chloride; then it was heated at 100 • C for 10 min.After cooling at room temperature, 2 mL water was added into the mixture, followed three times by extraction using 3 mL hexane.The extracts were dried under nitrogen at room temperature.The fatty-acid methyl ester (FAME) fraction was purified via a silica gel flash column using dichloromethane (DCM) as the eluent.After purification, the FAMEs were analyzed using a gas chromatography/mass spectrometry (GC-MS) system with an Agilent DB-1 column (30 m × 320 µm × 0.17 µm).Compounds were identified by comparison of mass spectra with the published and library data (Chen et al., 2022).

Hydrogen isotope analysis
The hydrogen isotope compositions of FAMEs were determined using an HP 6890 gas chromatograph interfaced to a MAT 253 isotope ratio mass spectrometer using a pyrolysis interface.The GC was fit with a wax column (30 m × 320 µm × 0.25 µm) and the temperature program was as follows: the oven temperature was held at 100 • C for 1 min, then increased by 20 • C min −1 to 200 • C and then followed by a 10 • C min −1 increase up to 310 • C and was isothermal for 7 min.The C 20 FAME standard with known δ 2 H values from Indiana University was measured after every sixth sample injection.The low and high relative-abundance FAMEs were analyzed separately to ensure proper signal amplitudes and to obtain accurate hydrogen isotope ratios.The δ 2 H values of FAMEs were corrected for the isotopic contribution of the hydrogen in the methyl group added during methylation using the following formula: δ 2 H fatty acids = [(2n + 2) × δ 2 H measured + 55.8 ‰ × 3]/(2n − 1), where n is the number of carbons in the fatty acids and −55.8 ‰ is the δ 2 H value of the added methyl group.The average standard deviation of the triplicate analysis of each compound was smaller than 5 ‰.Hydrogen isotope fractionations between fatty acids and growth water were calculated as ε FA/water = 1000(δ 2 H FA − δ 2 H water )/(δ 2 H water + 1000).
Culture water samples were analyzed using a MAT 253 isotope ratio mass spectrometer.The measurement was normalized to the Vienna Standard Mean Ocean Water and the Standard Light Antarctic Precipitation (VSMOW-SLAP) scale using a two-point linear calibration generated from reference waters supplied by IAEA, with an analytical precision and accuracy of ±0.6 ‰ for δ 2 H.

Fatty-acid distribution
The composition of fatty acids in S. piezotolerans WP3 growing on different organic substrates has been reported by our previous study (Chen et al., 2022).The lipid compositions are similar, growing on different organic substrates and temperatures.There are no systematic relationships between fatty-acid abundance and growth substrate.Briefly, the carbon chain lengths of n-alkanoic acids vary from C 12 to C 20 , peaking at C 16 and C 16 : 1 .Many kinds of low-melting-point fatty acids including branched-chain fatty acids (iC 14 , iC 15 , iC 16 , and iC 17 ), monounsaturated fatty acids (C 14 : 1 , C 16 : 1 , and C 18 : 1 ) and C 20 : 5 are found.The relative abundance of these three low-melting-point fatty acids displays a trend towards increasing by decreasing the culture temperature.

Hydrogen isotope ratio of fatty acids growing on different organic substrates
The hydrogen isotope composition of fatty acids is listed in Table 1.The δ 2 H values of fatty acids vary considerably between cultures on different substrates.On the whole, the hydrogen isotope composition of fatty acids is most depleted in 2 H growing on GlcNac and glucose and slightly 2 H-depleted on L-alanine (  1).The abundance-weighted mean δ 2 H values growing on GlcNac are −94 ± 1 ‰ (n = 2), and they are a little lower than that on glucose (−86 ± 35 ‰, n = 2).Fatty acids have a similar range on cultures when pyruvate (24 ± 4 ‰, n = 2) and acetate (3 ± 3 ‰, n = 2) are used as the sole organic carbon.The abundance-weighted mean δ 2 H values are −24 ± 9 ‰ (n = 2) and 32 ± 14 ‰ (n = 2) growing on Lalanine and L-glutamate, respectively.

Hydrogen isotope ratios of fatty acids at different growth temperatures
We investigate the hydrogen isotope composition of fatty acids in S. piezotolerans WP3 using GlcNac as the sole carbon growing at different temperatures (4, 10, 15, 20, and 25 • C).The abundance-weighted mean δ 2 H values exhibit considerable ranges, with lower values at optimal growth temperatures (15 and 20 • C) and higher values at nonoptimal temperatures (Table 2 and Fig. 2).Specifically, the mean δ 2 H values are −23 ± 2 ‰ (n = 2) and −23 ‰ growing at 15 and 20

Variations in fractionation growing on different organic substrates
The range of δ 2 H values for different fatty-acid structures in heterotrophic bacteria is often large (> 100 ‰) (X.Zhang et al., 2009;Osburn et al., 2016).In this study, S. piezotolerans WP3 grown on glucose exhibits a range of δ 2 H value varying from −176 ‰ to −40 ‰, while the range is relatively smaller on acetate (from −41 ‰-19 ‰) (Table 1 and Fig. 1).To facilitate quantitative comparison between 2 H/ 1 H fractionation and different organic substrates, we calculate an abundance-weighted mean hydrogen isotope fractionation between fatty acids and growth water (ε FA/water ) and standard deviation for each duplicate cultures.Our results show that ε FA/water values are substantially influenced by the type of organic carbon assimilated during growth (Fig. 1).Specifically, fatty acids are mostly enriched in 2 H growing on the direct precursor of the TCA cycle such as amino acids (Lglutamate and L-alanine), pyruvate, and acetate with average ε FA/water values from −4 ‰-52 ‰, while fatty acids exhibited 2 H depletion using GlcNac and glucose as the sole carbon source with ε FA/water values of −76 ‰ and −67 ‰, respectively.As demonstrated by many previous culture experiments (X.Zhang et al., 2009;Wijker et al., 2019), the variation in lipid δ 2 H values is not related to organic substrate isotopic composition because relatively little H from organic substrate is incorporate into fatty acids (Sessions et al., 1999;White et al., 2005).In contrast, the 2 H/ 1 H fractionation is substantially influenced by the central metabolic pathways, which are driven by the type of organic carbon source assimilated in all analyzed strains (X.Zhang et al., 2009;Dirghangi and Pagani, 2013a;Heinzelmann et al., 2015a, b;Wijker et al., 2019).
Comparison of our results with other aerobic heterotrophic bacterial cultures exhibits consistent 2 H/ 1 H fractionation patterns.For example, heterotrophs grown on glucose produced fatty acids were 2 H-depleted relative to source water (X.Zhang et al., 2009;Wijker et al., 2019).The metabolism of glucose by S. piezotolerans WP3 involves the Embden-Meyerhof pathway to pyruvate, and then pyruvate is transformed to acetyl-CoA via decarboxylation, and most of the acetyl-CoA enters into TCA cycle (Fig. 3).When grown on pyruvate and a direct precursor (acetate) of the TCA cycle, S. piezotolerans WP3 fatty acids were 2 H-enriched compared to water, similar to that observed in Cupriavidus necator strain H16 and Bacillus subtilis cultures (Wijker et al.,   2019; X. Zhang et al., 2009).Finally, when amino acids (Lalanine and L-glutamate) were used as organic substrates for heterotrophic bacteria, fatty acids were 2 H-enriched relative to water or there was no substantial 2 H/ 1 H fractionation.
According to the metabolism of these two amino acids in S. piezotolerans WP3 (Fig. 3), L-alanine was transformed to pyruvate that subsequently participates in the TCA cy-   2019), as NADPH is the dominant source (> 50 %) of hydrogen in fatty acids (Sessions, 2016).The dehydrogenase reactions and transhydrogenases interconversion between nicotinamide adenine dinucleotide (NADH) and NADPH are all associated with a large hydrogen isotopic fractionation (Wijker et al., 2019).Specifically, the NADP + is reduced through the pentose phosphate and glycolysis pathway resulting in 2 H depletion of NADPH; then fatty acids were 2 H-depleted relative to source water (Wijker et al., 2019;X. Zhang et al., 2009).In contrast, heterotrophic bacteria growing on TCA cycle organic substrates produce strongly 2 H-enriched NADPH via the associated NADPH-reducing reactions and the conversion of NADPH to NADH catalyzed by transhydrogenases.
For S. piezotolerans WP3 growing on GlcNac and glucose, NADPH is mainly produced via the pentose phosphate and glycolysis pathway (Fig. 3).Thus, relatively lower δ 2 H values of fatty acids are observed growing on GlcNac and glucose; the same order has also been observed previously (X.Zhang et al., 2009;Wijker et al., 2019).However, growth on glucose yields different 2 H/ 1 H fractionations in this study, with five heterotrophic bacteria investigated by  4; Wijker et al., 2019).Fatty acids in this study are −67 ‰ 2 H-depleted relative to growth water.It is surprising that the glucose metabolism in the above aerobic heterotrophic bacteria produces significantly different isotope ratios of NADPH.A systematic experiment by Wijker et al. (2019) demonstrated that changes in metabolic fluxes have an essential role in controlling the hydrogen isotope composition of lipids, with high ED and TCA fluxes relative to PP and EMP being associated with more 2 H-enriched lipids and vice versa (Wijker et al., 2019).Therefore, com-pared with R. radiobacter, E. meliloti, and P. fluorescens strains, relatively higher PP and EMP fluxes in S. piezotolerans WP3 may cause more 2 H-depleted lipids.
Another is that the conversion of pyruvate to amino acids also consumes NADPH (Fig. 3); this process may affect the hydrogen isotope composition of the NADPH pool.Further, different 2 H/ 1 H fractionations are also observed in different aerobic heterotrophic bacteria growing on a TCA cycle intermediate (acetate).Acetate can be metabolized directly to form acetyl-CoA, which is an important intermediate substrate for various biochemical reactions (Fig. 3).As a result, the pool of precursors or intermediates of the TCA cycle and δ 2 H values of NADPH utilized for lipid synthesis could be different, as reflected by the lipid 2 H/ 1 H fractionation. Up to now, however, there are no flux measurements for these investigated species growing on TCA cycle substrates.Therefore, isotope labeling experiments should be used to determine the associated metabolic flux growing on pyruvate and TCA cycle substrates for decoding the variations in 2 H/ 1 H fractionation in different strains.

Influence of growth temperature on the hydrogen isotope composition of fatty acids
In this study, the abundance-weighted mean δ 2 H values of fatty acids in S. piezotolerans WP3 exhibit considerable ranges, with lower values at optimal growth temperatures (15 and 20 • C) and higher values at non-optimal temperatures (4, 10, and 25 • C; Fig. S2).Many previous culture experiments have shown that growth temperature exerts significant effects on the hydrogen isotope composition of lipids in photoautotrophic algae (Z.Zhang et al., 2009;Schouten et al., 2006).The lipid δ 2 H values are decreased by increasing growth temperatures below the optimal growth temper- ature, while values increase by further increasing temperatures above the optimal growth temperature (Schouten et al., 2006;Z. Zhang et al., 2009) These similar results are also observed in a halophilic archaean Haloarcula marismortui studied by Dirghangi and Pagani (2013b).In contrast, one culture study (Tetrahymena thermophila) indicates fatty acids are more 2 H-enriched by increasing growth tempera-ture (Dirghangi and Pagani, 2013b).The mechanisms may be associated with growth rates and enzyme activities in organisms controlled by temperature (Dirghangi and Pagani, 2013b;X. Zhang et al., 2009).The possible mechanisms are that temperature affects lipid hydrogen isotope composition through growth rates and enzyme activities in organisms (Dirghangi and Pagani, 2013b;X. Zhang et al., 2009).
https://doi.org/10.5194/bg-20-1491-2023Biogeosciences, 20, 1491-1504, 2023 The effect of growth rates on hydrogen isotope composition is mainly determined by the proportion of intracellular water from the metabolism (Kreuzer-Martin et al., 2006).Enzyme activities are also strongly related to temperature because organisms produce NADPH through different metabolisms in response to temperature change (X.Zhang et al., 2009).We will discuss these potential factors separately below.S. piezotolerans WP3 is psychrotolerant and can grow at 0-28 • C, with an optimal growth temperature at 15 to 20 • C (Wang et al., 2004).The growth rate of S. piezotolerans WP3 is temperature dependent on the highest rates of optimal growth temperature (Xiao et al., 2007).Our results show that relatively larger fractionation is present at a higher growth rate (15 and 20 • C), while smaller fractionation is observed at relatively lower rates (4, 10, and 25 • C; Fig. 2).This is similar to previous results according to which halophilic archaea yield more 2 H-depleted lipids at an optimal growth temperature, while lipids become 2 Henriched (11 ‰) with the temperature increasing above 40 • C (Dirghangi and Pagani, 2013a).These observations point to growth rate being the dominant factor causing the observed variations in lipid 2 H/ 1 H fractionations.A probable mechanism is that the growth rate affects lipid δ 2 H via hydrogen atom exchange between NADPH and intracellular water (Wijker et al., 2019).A culture experiment by X. Zhang et al. (2009) reported that lipids showed markedly different δ 2 H values with two heterotrophic bacteria (C.necator and E. coli) grown on the same organic substrate, with greater fractionation for faster-growing microbes.Hydrogen isotope exchange between NADPH and water in heterotrophs is lower under higher growth rates compared to lower growth rates.However, the strength of hydrogen isotope exchange between NADPH and water depends on the turnover time of NADPH, with a longer cellular residence time for more significant hydrogen isotope exchange.Recently, an important experiment with six heterotrophic bacteria provided evidence that no significant correlation between NADPH turnover time and 2 H/ 1 H fractionation was observed, suggesting such exchange effects were minor (Wijker et al., 2019).By contrast, variation in 2 H/ 1 H fractionations in different strains growing on glucose is mainly driven by fluxes through NADP + -reducing and NADPH-NADHbalancing reactions.Thus, the observed correlations between 2 H/ 1 H fractionation and growth temperature in this study are likely caused by other mechanisms associated with growth rate.
The different central metabolic pathways produce NADPH with a characteristic isotope composition due to differences in the kinetic isotope effects (KIEs) accompanying NADP + reduction by dehydrogenases and transhydrogenases (X.Zhang et al., 2009).The relative carbon flux associated with G6PDH and 6PGDH enzymes shows a negative relationship with ε FA/water values, while positive correlations are shown for ME (ED pathway) and ICDH (TCA cycle; Wijker et al., 2019).Only increasing carbon flux through G6PDH and 6PGDH can result in decreasing δ 2 H values of fatty acids.In the S. piezotolerans WP3 strain, the production of NADPH growing on GlcNac is mainly associated with three dehydrogenase enzymes including glucose-6-phosphate dehydrogenase (G6PDH), 6phosphogluconate dehydrogenase (6PGDH), and isocitrate dehydrogenase (ICDH; Fig. 3).Under the optimal growth temperature, relatively larger isotope fractionations are most likely related to the higher fluxes of NADPH generated by the enzymes G6PDH and 6PGDH.According to the culture experiment and transcriptome data, the expression level of G6PDH and 6PGDH encoding genes at an optimal growth temperature is higher than at non-optimal temperatures (Meng et al., 2020).Given the above constraints, changes in the fluxes of NADP + reduction by the associated enzymes G6PDH and 6PGDH can plausibly explain the variation in lipid 2 H/ 1 H fractionation grown at different temperatures.
The interconversion between NADPH and NADH by transhydrogenase enzymes is also regarded as a potential factor for determining the lipid hydrogen isotope composition (Wijker et al., 2019).The ε FA/water values displayed positive correlations with the NADPH transhydrogenase fluxes in aerobic heterotrophic bacteria, with NADPH underproduction corresponding to lower ε FA/water values (Wijker et al., 2019).The supply of NADPH by enzymatic reactions should be nearly balanced with the demand for NADPH in fattyacid biosynthesis reactions (Fuhrer and Sauer, 2009).Under high growth rates, the NADPH is overproduced; thus a portion of NADPH should be converted to NADH for balancing between catabolic and anabolic fluxes (Spaans et al., 2015).This process will lead to 2 H-enrichment in residual NADPH, indirectly causing 2 H-enriched lipids.Therefore, if the predominant interconversion from NADPH to NADH occurs in S. piezotolerans WP3 under high growth rates (optimal growth temperature), the positive relationship between fattyacid δ 2 H values and growth rates should be observed.However, the estimates are not in accordance with our measurements, indicating transhydrogenase-enzyme-associated reactions are negligible in this study.According to the above discussions, therefore, we assumed that the effect of growth temperature on the fatty-acid δ 2 H value is primarily related to the activity of G6PDH and 6PGDH enzymes in S. piezotolerans WP3.

Geochemical implications
We summarized previously published data (including this study) on the 2 H/ 1 H fractionations in aerobic heterotrophic bacteria growing on a variety of organic substrates, including sugars, pyruvate, and TCA cycle substrates (Figs. 4  and 5).On the whole, variations in 2 H/ 1 H fractionation are strongly correlated with the central metabolic pathways because different metabolism-associated enzymes produce NADPH with different hydrogen isotope compositions (Wi- jker et al., 2019).Specifically, fatty acids have average ε lipids/water values ranging from −100 ‰ to +100 ‰ using sugars as substrates and from ∼ −50 ‰ to +200 ‰ and +100 ‰ to +400 ‰ using pyruvate and TCA cycle substrates, respectively (Figs. 4 and 5).However, considerable variations in fractionation are observed among these investigated strains growing on the same organic substrates, and they are predominantly driven by metabolic fluxes (Wijker et al., 2019).Previous laboratory experiments showed that the average ε FA/water values in P. fluorescens and Rhizobia growing on fructose are about 25 % lower than those growing on glucose (Wijker et al., 2019).The ED pathway is predominant in the glucose metabolism for P. fluorescens, while the flux is only about 54 % through the ED pathway growing on fructose.The metabolic flux analysis can also excellently interpret the 2 H/ 1 H fractionations in E. coli growing on gluconate and glucose (X.Zhang et al., 2009).For heterotrophs growing on sugars, B. subtilis and E. coli exhibit the most 2 H-depleted fatty acids, while two species (P.fluorescens and Pseudomonas str.LFY10) are slightly enriched in 2 H composition (Fig. 4).This is similar to the flux data, according to which B. subtilis and E. coli strongly rely on the EMP pathway for glucose catabolism, while the TCA cycle is predominant in P. fluorescens and Pseudomonas str.LFY10 (Heinzelmann et al., 2015b).For S. piezotolerans WP3 growing on glucose and GlcNac, the average ε lipids/water values are −67 ‰ and −76 ‰, respectively, and they are higher than B. subtilis and E. coli but lower than P. fluorescens and Pseudomonas str.LFY10 (Fig. 4).We speculate that the metabolic flux of the TCA cycle in S. piezotolerans WP3 may vary between E. coli and P. fluorescens.
It is interesting to note that aerobic heterotrophs growing on TCA cycle substrates (acetate, succinate, citrate, and benzoic acid) exhibit the largest variations in ε lipids/water values (from ∼ −50 ‰-400 ‰).Heterotrophs including P. fluorescens and Pseudomonas also have relatively higher ε lipids/water values compared with other investigated strains (Fig. 5).The NADPH is predominately produced through the TCA cycle, which results in strong 2 H enrichment of NADPH.Moreover, overproduction of NADPH is generally present in heterotrophic bacteria (e.g., E. coli) growing on the direct precursors or intermediates of the TCA cycle (Gerosa et al., 2015;Wijker et al., 2019).The most important thing is that a proportion of NADPH is firstly used as an H − donor for amino acid and protein synthesis, causing residual NADPH enriched in 2 H (Fig. 3).Both of these two processes can yield fatty acids with strong 2 H enrichment, and the metabolism and range of NADPH overproduction may be responsible for the different ε lipids/water values in these species.
Compound-specific isotope measurements can provide further insights into the biogeochemical cycle and ecological change in deep oceans.Many previous studies reported that C 20 : 5 fatty acid can be used as the specific lipids for phytoplankton (Heinzelmann et al., 2016).However, our data showed that S. piezotolerans WP3 can also produce a certain amount of C 20 : 5 fatty acid (Wang et al., 2009;Chen et al., 2022).Fatty acids produced from phototrophs have relative larger hydrogen and carbon isotope fractionations with ε FA/water and δ 13 C FA/Substrate values ranging from about −200 ‰ to −150 ‰ and about −28 ‰, respectively (Jones et al., 2008;Li et al., 2009;Guan et al., 2019).By contrast, heterotrophs obviously have different fractionations, with ε FA/water and δ 13 C FA/Substrate values from about −100 ‰-100 ‰ and −10 ‰-4 ‰.Among heterotrophs, the ε FA/water and δ 13 C FA/Substrate are characterized by different values assimilating different types of organic matter (Fig. 6).The carbon and hydrogen isotope compositions of these specific biomarker lipids are a promising tool to infer the community metabolisms and to then decode carbon cycling processes in Fe-reducing oceanic sediments.

Conclusions
Our results showed that the hydrogen isotope composition of individual fatty acids in S. piezotolerans WP3 varies considerably in the same culture.The branched fatty acids and polyunsaturated fatty acid (C 20 : 5 ) have the lowest δ 2 H values among these compounds, indicating that different biosynthesis pathways have an obvious influence on the δ 2 H values of different-chain-length fatty acids.The hydrogen isotope fractionations between fatty acids and growth water vary systematically by the type of central carbon metabolic pathway.Specifically, WP3 growing on GlcNac and glucose produce NADPH through pentose phosphate and glycolysis pathways; fatty acids are most depleted in 2 H. Fatty acids are slightly enriched in 2 H using pyruvate and L-alanine as the sole carbon source and strongly 2 H-enriched using TCA cycle substrates (acetate and L-glutamate) because NADP + is produced through the TCA cycle pathway.The most plausible reasons are that hydrogen isotope fractionations in aerobic heterotrophs are highly correlated to the type of central metabolic pathways due to different enzymes accompanying NADP + reduction.
We also observed large variations in fractionation in WP3 using GlcNac as the sole carbon source at different temperatures.Fatty acids have relatively lower ε FA/water values at optimal temperatures (−3 ± 3 ‰ and 9 ‰ at 15 and 20 • C, respectively), while there are higher values at 4 • C (24 ± 5 ‰), 10 • C (15 ± 12 ‰), and 25 • C (35 ± 41 ‰).We assumed that the influence of growth temperature on the δ 2 H values of fatty acids is most likely related to the activity of G6PDH and 6PGDH enzymes.

Figure 1 .
Figure 1.The hydrogen isotope fractionation between fatty acids and growth water in S. piezotolerans WP3 growing on different substrates including GlcNac, acetate, D-glucose, pyruvate, and amino acids at 20 • C and shaking at 200 rpm.The error bars represent the 1σ standard deviations of duplicate cultures.
cle, while L-glutamate was directly metabolized via the TCA cycle pathway to acetyl-CoA.All these results indicate that variation in the 2 H/ 1 H fractionation in S. piezotolerans WP3 grown on different types of organic substrates is most likely correlated with metabolisms.Processes associated with NADPH production are crucial in controlling the lipid 2 H/ 1 H fractionation(Wijker et al., https://doi.org/10.5194/bg-20-1491-2023Biogeosciences, 20, 1491-1504, 2023

Figure 2 .
Figure 2. The hydrogen isotope composition of fatty acids in S. piezotolerans WP3 growing at different temperatures.Panel (f) shows the abundance-weighted mean ε FA/water values at different temperatures.The error bars represent the 1σ standard deviations of duplicate cultures.None of the cultures have been shaken.

Figure 5 .
Figure 5. Compilation of measured hydrogen isotope fractionations growing on pyruvate and TCA cycle substrates.For each organism, bars represent the range of values for all individual fatty acids.The vertical line is the average value.Data sources: 1 this study; 2 Wijker et al. (2019); 3 X.Zhang et al. (2009); 4,5 Heinzelmann et al. (2015a).

Figure 6 .
Figure 6.The values of carbon and hydrogen fractionation of fatty acids in S. piezotolerans WP3 growing on different organic substrates.The δ 13 C FA/Substrate values are cited from Chen et al. (2022).
Fang et al. (2014)et al. (2014), polyunsaturated fatty acids (C 22:6 ) synthesized through the polyketide pathway have much lower δ 2 H values than C 16 via the fatty-acid synthase (FAS) pathway in a piezophilic bacterium Moritella japonica DSK1.All these results suggest that fatty-acid biosynthesis processes including desaturation, different precursors, and biosynthetic pathway have considerable influence on the hydrogen isotope composition in S. piezotolerans WP3.
• C, respectively.WP3 growing at 4, 10, and 25 • C has mean δ 2 H values of 4 ± 5 ‰, −4 ± 12 ‰, and to 15 ± 41 ‰.4 Discussions4.1 Hydrogen isotope composition of individual fatty acidsThe δ 2 H values of individual fatty acids vary considerably in the same culture with values ranging from 60 ‰-138 ‰.This indicates biosynthesis pathways have an obvious influence on the δ 2 H values of individual fatty acids.Specifically, the different chain length of saturated fatty acids have similar δ 2 H values, suggesting chain elongation processes by enzymes have no obvious influence on hydrogen isotope values.The δ 2 H values of monounsaturated fatty acids (C 16:1 ) are slightly lower than those of saturated fatty acids (C jker et al. (2019) also found that branched fatty acids (iC 15 and aC 15 ) have relatively lower δ 2 H values compared with C 16 fatty acid in five species of heterotrophic bacterium.It is interesting to note that C 13 and C 20:5 in S. piezotolerans WP3 have the lowest δ 2 H values among these fatty acids.

Table 1 .
The hydrogen isotopic values of fatty acids in S. piezotolerans WP3 growing on different organic substrates at 20 • C and shaking at 200 rpm; "-" indicates that the concentrations are too low for δ 2 H measurements.Standard deviations of two duplicate cultures are given in parentheses.

Table 2 .
The hydrogen isotopic values of fatty acids in S. piezotolerans WP3 growing at different temperatures using GlcNac as the sole carbon source.For better comparison, none of the cultures have been shaken; "-" indicates that the concentrations are too low for δ 2 H measurements.Standard deviations of two duplicate cultures are given in parentheses.