2H/1H variation in microbial lipids is controlled by NADPH metabolism

Hydrogen isotope fractionation is controlled by CO2 in coccolithophore lipids

Hydrogen isotope ratios (δ2H) represent an important natural tracer of metabolic processes, but quantitative models of processes controlling H-fractionation in aquatic photosynthetic organisms are lacking. Here, we elucidate the underlying physiological controls of 2H/1H fractionation in algal lipids by systematically manipulating temperature, light, and CO2(aq) in continuous cultures of the haptophyte Gephyrocapsa oceanica. We analyze the hydrogen isotope fractionation in alkenones (αalkenone), a class of acyl lipids specific to this species and other haptophyte algae. We find a strong decrease in the αalkenone with increasing CO2(aq) and confirm αalkenone correlates with temperature and light. Based on the known biosynthesis pathways, we develop a cellular model of the δ2H of algal acyl lipids to evaluate processes contributing to these controls on fractionation. Simulations show that longer residence times of NADPH in the chloroplast favor a greater exchange of NADPH with 2H-richer intracellular water, increasing αalkenone. Higher chloroplast CO2(aq) and temperature shorten NADPH residence time by enhancing the carbon fixation and lipid synthesis rates. The inverse correlation of αalkenone to CO2(aq) in our cultures suggests that carbon concentrating mechanisms (CCM) do not achieve a constant saturation of CO2 at the Rubisco site, but rather that chloroplast CO2 varies with external CO2(aq). The pervasive inverse correlation of αalkenone with CO2(aq) in the modern and preindustrial ocean also suggests that natural populations may not attain a constant saturation of Rubisco with the CCM. Rather than reconstructing growth water, αalkenone may be a powerful tool to elucidate the carbon limitation of photosynthesis.

Large enrichments in fatty acid 2H/1H ratios distinguish respiration from aerobic fermentation in yeast Saccharomyces cerevisiae

Shifts in the hydrogen stable isotopic composition (2H/1H ratio) of lipids relative to water (lipid/water 2H-fractionation) at natural abundances reflect different sources of the central cellular reductant, NADPH, in bacteria. Here, we demonstrate that lipid/water 2H-fractionation (2εfattyacid/water) can also constrain the relative importance of key NADPH pathways in eukaryotes. We used the metabolically flexible yeast Saccharomyces cerevisiae, a microbial model for respiratory and fermentative metabolism in industry and medicine, to investigate 2εfattyacid/water. In chemostats, fatty acids from glycerol-respiring cells were >550‰ 2H-enriched compared to those from cells aerobically fermenting sugars via overflow metabolism, a hallmark feature in cancer. Faster growth decreased 2H/1H ratios, particularly in glycerol-respiring cells by 200‰. Variations in the activities and kinetic isotope effects among NADP+-reducing enzymes indicate cytosolic NADPH supply as the primary control on 2εfattyacid/water. Contributions of cytosolic isocitrate dehydrogenase (cIDH) to NAPDH production drive large 2H-enrichments with substrate metabolism (cIDH is absent during fermentation but contributes up to 20 percent NAPDH during respiration) and slower growth on glycerol (11 percent more NADPH from cIDH). Shifts in NADPH demand associated with cellular lipid abundance explain smaller 2εfattyacid/water variations (<30‰) with growth rate during fermentation. Consistent with these results, tests of murine liver cells had 2H-enriched lipids from slower-growing, healthy respiring cells relative to fast-growing, fermenting hepatocellular carcinoma. Our findings point to the broad potential of lipid 2H/1H ratios as a passive natural tracker of eukaryotic metabolism with applications to distinguish health and disease, complementing studies that rely on complex isotope-tracer addition methods.

Energy flux couples sulfur isotope fractionation to proteomic and metabolite profiles in Desulfovibrio vulgaris

Microbial sulfate reduction is central to the global carbon cycle and the redox evolution of Earth's surface. Tracking the activity of sulfate reducing microorganisms over space and time relies on a nuanced understanding of stable sulfur isotope fractionation in the context of the biochemical machinery of the metabolism. Here, we link the magnitude of stable sulfur isotopic fractionation to proteomic and metabolite profiles under different cellular energetic regimes. When energy availability is limited, cell-specific sulfate respiration rates and net sulfur isotope fractionation inversely covary. Beyond net S isotope fractionation values, we also quantified shifts in protein expression, abundances and isotopic composition of intracellular S metabolites, and lipid structures and lipid/water H isotope fractionation values. These coupled approaches reveal which protein abundances shift directly as a function of energy flux, those that vary minimally, and those that may vary independent of energy flux and likely do not contribute to shifts in S-isotope fractionation. By coupling the bulk S-isotope observations with quantitative proteomics, we provide novel constraints for metabolic isotope models. Together, these results lay the foundation for more predictive metabolic fractionation models, alongside interpretations of environmental sulfur and sulfate reducer lipid-H isotope data.

M. Rangel W, Kopáček M, Jílková V, Jansa J, Meador TB. Constraining activity and growth substrate of fungal decomposers via assimilation patterns of inorganic carbon and water into lipid biomarkers

Fungi are among the few organisms on the planet that can metabolize recalcitrant carbon (C) but are also known to access recently produced plant photosynthate. Therefore, improved quantification of growth and substrate utilization by different fungal ecotypes will help to define the rates and controls of fungal production, the cycling of soil organic matter, and thus the C storage and CO2 buffering capacity in soil ecosystems. This pure-culture study of fungal isolates combined a dual stable isotope probing (SIP) approach, together with rapid analysis by tandem pyrolysis-gas chromatography-isotope ratio mass spectrometry to determine the patterns of water-derived hydrogen (H) and inorganic C assimilated into lipid biomarkers of heterotrophic fungi as a function of C substrate. The water H assimilation factor (αW) and the inorganic C assimilation into C18:2 fatty acid isolated from five fungal species growing on glucose was lower (0.62% ± 0.01% and 4.7% ± 1.6%, respectively) than for species grown on glutamic acid (0.90% ± 0.02% and 7.4% ± 3.7%, respectively). Furthermore, the assimilation ratio (RIC/αW) for growth on glucose and glutamic acid can distinguish between these two metabolic modes. This dual-SIP assay thus delivers estimates of fungal activity and may help to delineate the predominant substrates that are respired among a matrix of compounds found in natural environments.IMPORTANCEFungal decomposers play important roles in food webs and nutrient cycling because they can feed on both labile and more recalcitrant forms of carbon. This study developed and applied a dual stable isotope assay (13C-dissolved inorganic carbon/2H) to improve the investigation of fungal activity in the environment. By determining the incorporation patterns of hydrogen and carbon into fungal lipids, this assay delivers estimates of fungal activity and the different metabolic pathways that they employ in ecological and environmental systems.

Biosynthetic and catabolic pathways control amino acid δ2H values in aerobic heterotrophs

The hydrogen isotope ratios (δ2HAA values) of amino acids in all organisms are substantially fractionated relative to growth water. In addition, they exhibit large variations within microbial biomass, animals, and human tissues, hinting at rich biochemical information encoded in such signals. In lipids, such δ2H variations are thought to primarily reflect NADPH metabolism. Analogous biochemical controls for amino acids remain largely unknown, but must be elucidated to inform the interpretation of these measurements. Here, we measured the δ2H values of amino acids from five aerobic, heterotrophic microbes grown on different carbon substrates, as well as five Escherichia coli mutant organisms with perturbed NADPH metabolisms. We observed similar δ2HAA patterns across all organisms and growth conditions, which-consistent with previous hypotheses-suggests a first-order control by biosynthetic pathways. Moreover, δ2HAA values varied systematically with the catabolic pathways activated for substrate degradation, with variations explainable by the isotopic compositions of important cellular metabolites, including pyruvate and NADPH, during growth on each substrate. As such, amino acid δ2H values may be useful for interrogating organismal physiology and metabolism in the environment, provided we can further elucidate the mechanisms underpinning these signals.

Time-resolved, deuterium-based fluxomics uncovers the hierarchy and dynamics of sugar processing by Pseudomonas putida

Pseudomonas putida, a microbial host widely adopted for metabolic engineering, processes glucose through convergent peripheral pathways that ultimately yield 6-phosphogluconate. The periplasmic gluconate shunt (PGS), composed by glucose and gluconate dehydrogenases, sequentially transforms glucose into gluconate and 2-ketogluconate. Although the secretion of these organic acids by P. putida has been extensively recognized, the mechanism and spatiotemporal regulation of the PGS remained elusive thus far. To address this challenge, we adopted a dynamic 13C- and 2H-metabolic flux analysis strategy, termed D-fluxomics. D-fluxomics demonstrated that the PGS underscores a highly dynamic metabolic architecture in glucose-dependent batch cultures of P. putida, characterized by hierarchical carbon uptake by the PGS throughout the cultivation. Additionally, we show that gluconate and 2-ketogluconate accumulation and consumption can be solely explained as a result of the interplay between growth rate-coupled and decoupled metabolic fluxes. As a consequence, the formation of these acids in the PGS is inversely correlated to the bacterial growth rate-unlike the widely studied overflow metabolism of Escherichia coli and yeast. Our findings, which underline survival strategies of soil bacteria thriving in their natural environments, open new avenues for engineering P. putida towards efficient, sugar-based bioprocesses.

Mixoplankton and mixotrophy: future research priorities

Phago-mixotrophy, the combination of photoautotrophy and phagotrophy in mixoplankton, organisms that can combine both trophic strategies, have gained increasing attention over the past decade. It is now recognized that a substantial number of protistan plankton species engage in phago-mixotrophy to obtain nutrients for growth and reproduction under a range of environmental conditions. Unfortunately, our current understanding of mixoplankton in aquatic systems significantly lags behind our understanding of zooplankton and phytoplankton, limiting our ability to fully comprehend the role of mixoplankton (and phago-mixotrophy) in the plankton food web and biogeochemical cycling. Here, we put forward five research directions that we believe will lead to major advancement in the field: (i) evolution: understanding mixotrophy in the context of the evolutionary transition from phagotrophy to photoautotrophy; (ii) traits and trade-offs: identifying the key traits and trade-offs constraining mixotrophic metabolisms; (iii) biogeography: large-scale patterns of mixoplankton distribution; (iv) biogeochemistry and trophic transfer: understanding mixoplankton as conduits of nutrients and energy; and (v) in situ methods: improving the identification of in situ mixoplankton and their phago-mixotrophic activity.

The metabolic sensitivity of hydrogen isotope fractionation differs between plant compounds

Hydrogen stable isotope analyses (δ²H) of plant derived organic compounds are a useful tool for ecological, environmental, and palaeoclimatological research. However, during organic compound synthesis, variable biosynthetic ²H-fractionation has been suggested to occur as a result of changes in plant carbon fluxes. So far, inference has been based on examining the δ²H patterns of plant compounds along environmental gradients, among plant species, and between plant organs. In an alternative approach, we used four plant species with four different types of mutations that cause impaired starch synthesis to assess whether variability in carbon metabolism affects the biosynthetic ²H-fractionation during cellulose, phytol, and acetogenic lipid synthesis. We found that mutants with impaired starch synthesis always had higher cellulose and phytol δ²H values compared to the wild type. By contrast, ²H-fractionation during acetogenic lipid biosynthesis generally did not show strong metabolic sensitivity. We rationalise these differences by considering the biosynthetic pathway of each compound and the likely source of the variable isotope fractionation. In different organic compounds, the sensitivity of variable biosynthetic ²H-fractionation to changes in C-metabolism depends on incorporation of specific H atoms from precursor molecules. As such, we determined that the similar increase in cellulose and phytol δ²H values as an effect of impaired starch synthesis most likely originates in triose-phosphates.

In situ diversity of metabolism and carbon use efficiency among soil bacteria

The central carbon (C) metabolic network harvests energy to power the cell and feed biosynthesis for growth. In pure cultures, bacteria use some but not all of the network's major pathways, such as glycolysis and pentose phosphate and Entner-Doudoroff pathways. However, how these pathways are used in microorganisms in intact soil communities is unknown. Here, we analyzed the incorporation of ¹³C from glucose isotopomers into phospholipid fatty acids. We showed that groups of Gram-positive and Gram-negative bacteria in an intact agricultural soil used different pathways to metabolize glucose. They also differed in C use efficiency (CUE), the efficiency with which a substrate is used for biosynthesis. Our results provide experimental evidence for diversity among microbes in the organization of their central carbon metabolic network and CUE under in situ conditions. These results have important implications for our understanding of how community composition affects soil C cycling and organic matter formation.

Species variation in the hydrogen isotope composition of leaf cellulose is mostly driven by isotopic variation in leaf sucrose

Experimental approaches to isolate drivers of variation in the carbon-bound hydrogen isotope composition (δ² H) of plant cellulose are rare and current models are limited in their application. This is in part due to a lack in understanding of how ² H-fractionations in carbohydrates differ between species. We analysed, for the first time, the δ² H of leaf sucrose along with the δ² H and δ¹⁸ O of leaf cellulose and leaf and xylem water across seven herbaceous species and a starchless mutant of tobacco. The δ² H of sucrose explained 66% of the δ² H variation in cellulose (R²  = 0.66), which was associated with species differences in the ² H enrichment of sucrose above leaf water ( ε sucrose <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0001" wiley:location="equation/pce14362-math-0001.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><msub><mtext>\unicode{x003B5}</mtext><mtext>sucrose</mtext></msub></mrow></math> : -126% to -192‰) rather than by variation in leaf water δ² H itself. ε sucrose <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0002" wiley:location="equation/pce14362-math-0002.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><msub><mtext>\unicode{x003B5}</mtext><mtext>sucrose</mtext></msub></mrow></math> was positively related to dark respiration (R²  = 0.27), and isotopic exchange of hydrogen in sugars was positively related to the turnover time of carbohydrates (R²  = 0.38), but only when ε sucrose <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0003" wiley:location="equation/pce14362-math-0003.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mrow><msub><mi mathvariant="normal">\unicode{x003B5}</mi><mtext>sucrose</mtext></msub></mrow></mrow></math> was fixed to the literature accepted value of - 171 <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0004" wiley:location="equation/pce14362-math-0004.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mrow><mo>\unicode{x02212}</mo><mn>171</mn></mrow></mrow></math> ‰. No relation was found between isotopic exchange of hydrogen and oxygen, suggesting large differences in the processes shaping post-photosynthetic fractionation between elements. Our results strongly advocate that for robust applications of the leaf cellulose hydrogen isotope model, parameterization utilizing δ² H of sugars is needed.

Central Metabolism and Growth Rate Impacts on Hydrogen and Carbon Isotope Fractionation During Amino Acid Synthesis in E. coli

Compound specific stable isotope analysis (CSIA) of amino acids from bacterial biomass is a newly emerging powerful tool for exploring central carbon metabolism pathways and fluxes. By comparing isotopic values and fractionations relative to water and growth substrate, the impact of variable flow path for metabolites through different central metabolic pathways, perturbations of these paths, and their resultant consequences on intracellular pools and resultant biomass may be elucidated. Here, we explore the effects that central carbon metabolism and growth rate can have on stable hydrogen (δ²H) and carbon (δ¹³C) compound specific isotopic values of amino acids, and whether diagnostic isotopic fingerprints are revealed by these paired analyses. We measured δ²H and δ¹³C in amino acids in the wild type Escherichia coli (MG1655) across a range of growth rates in chemostat cultures to address the unknown isotopic consequences as metabolic fluxes are shuffled between catabolic and anabolic metabolisms. Additionally, two E. coli knockout mutants, one with deficiency in glycolysis –pgi (LC1888) and another inhibiting the oxidative pentose phosphate pathway (OPPP) –zwf (LC1889), were grown on glucose and used as a comparison against the wild type E. coli (MG1655) to address the isotopic changes of amino acids produced in these perturbed metabolic pathways. Amino acid δ²H values, which collectively vary in composition by more than 400‰, are altered along with δ¹³C values demonstrating fundamental shifts in central metabolic pathways and/or fluxes. Within our linear discriminant analysis with a simple model organism to examine potential amino acid fingerprinting, our knockout strains and variable growth rate samples plot across a wider array of organism classification than merely within the boundaries of other bacterial data.

Deuterium in marine organic biomarkers: toward a new tool for quantifying aquatic mixotrophy

The traditional separation between primary producers (autotrophs) and consumers (heterotrophs) at the base of the marine food web is being increasingly replaced by the paradigm that mixoplankton, planktonic protists with the nutritional ability to use both phago(hetero)trophy and photo(auto)trophy to access energy are widespread globally. Thus, many 'phytoplankton' eat, while 50% of 'protozooplankton' also perform photosynthesis. Mixotrophy may enhance primary production, biomass transfer to higher trophic levels and the efficiency of the biological pump to sequester atmospheric CO₂ into the deep ocean. Although this view is gaining traction, science lacks a tool to quantify the relative contributions of autotrophy and heterotrophy in planktonic protists. This hinders our understanding of their impacts on carbon cycling within marine pelagic ecosystems. It has been shown that the hydrogen (H) isotopic signature of lipids is uniquely sensitive to heterotrophy relative to autotrophy in plants and bacteria. Here, we explored whether it is also sensitive to the trophic status in protists. The new understanding of H isotope signature of lipid biomarkers suggests it offers great potential as a novel tool for quantifying the prevalence of mixotrophy in diverse marine microorganisms and thus for investigating the implications of the 'mixoplankton' paradigm.

Inverse hydrogen isotope fractionation indicates heterotrophic microbial production of long-chain n-alkyl lipids in desolate Antarctic ponds

Long-chain (>C₂₅ ) n-alkyl lipids have long been considered biomarkers for higher plant leaf waxes and widely applied for paleoclimate and paleoenvironmental reconstructions. However, recent experimental and lacustrine sediment studies suggest long-chain n-alkanes and n-alkanoic acids can also be produced by aerobic microbes, probably heterotrophic microbes based on carbon isotope data. Here we show that sedimentary long-chain n-alkanes and n-alkanoic acids in two desolate Antarctic ponds where vascular plants are absent in the surroundings display hydrogen isotopic values up to 300 per mil higher than those of lake water. It is the first time that such strongly inverse or reduced hydrogen isotopic fractionation of lipid biomarkers is observed in natural sediment samples. Based on recent extensive experimental data on microbial hydrogen isotopic fractionation, our data can only be explained by the predominant production of long-chain n-alkanes and n-alkanoic acids from heterotrophic micro-organisms. Together with preliminary 16S rRNA gene sequencing data, our results represent the first unambiguous example of predominant heterotrophic microbial production of long-chain n-alkyl waxes in a natural environment.

Reconfiguration of metabolic fluxes in Pseudomonas putida as a response to sub-lethal oxidative stress

As a frequent inhabitant of sites polluted with toxic chemicals, the soil bacterium and plant-root colonizer Pseudomonas putida can tolerate high levels of endogenous and exogenous oxidative stress. Yet, the ultimate reason of such phenotypic property remains largely unknown. To shed light on this question, metabolic network-wide routes for NADPH generation-the metabolic currency that fuels redox-stress quenching mechanisms-were inspected when P. putida KT2440 was challenged with a sub-lethal H₂O₂ dose as a proxy of oxidative conditions. ¹³C-tracer experiments, metabolomics, and flux analysis, together with the assessment of physiological parameters and measurement of enzymatic activities, revealed a substantial flux reconfiguration in oxidative environments. In particular, periplasmic glucose processing was rerouted to cytoplasmic oxidation, and the cyclic operation of the pentose phosphate pathway led to significant NADPH-forming fluxes, exceeding biosynthetic demands by ~50%. The resulting NADPH surplus, in turn, fueled the glutathione system for H₂O₂ reduction. These properties not only account for the tolerance of P. putida to environmental insults-some of which end up in the formation of reactive oxygen species-but they also highlight the value of this bacterial host as a platform for environmental bioremediation and metabolic engineering.

Determination of the δ2 H values of high molecular weight lipids by high-temperature gas chromatography coupled to isotope ratio mass spectrometry

RATIONALE

The hydrogen isotopic composition of lipids (δ² H_lipid ) is widely used in food science and as a proxy for past hydrological conditions. Determining the δ² H values of large, well-preserved triacylglycerides and other microbial lipids, such as glycerol dialkyl glycerol tetraether (GDGT) lipids, is thus of widespread interest but has so far not been possible due to their low volatility which prohibits analysis by traditional gas chromatography/pyrolysis/isotope ratio mass spectrometry (GC/P/IRMS).

METHODS

We determined the δ² H values of large, polar molecules and applied high-temperature gas chromatography (HTGC) methods on a modified GC/P/IRMS system. The system used a high-temperature 7-m GC column, and a glass Y-splitter for low thermal mass. Methods were validated using authentic standards of large, functionalised molecules (triacylglycerides, TGs), purified standards of GDGTs. The results were compared with δ² H values determined by high-temperature elemental analyser/pyrolysis/isotope ratio mass spectrometry (HTEA/P/IRMS), and subsequently applied to the analysis of GDGTs in a sample from a methane seep and a Welsh peat.

RESULTS

The δ² H values of TGs agreed within error between HTGC/P/IRMS and HTEA/IRMS, with HTGC/P/IRMS showing larger errors. Archaeal lipid GDGTs with up to three cyclisations could be analysed: the δ² H values were not significantly different between methods with standard deviations of 5 to 6 ‰. When environmental samples were analysed, the δ² H values of isoGDGTs were 50 ‰ more negative than those of terrestrial brGDGTs.

CONCLUSIONS

Our results indicate that the HTGC/P/IRMS method developed here is appropriate to determine the δ² H values of TGs, GDGTs with up to two cyclisations, and potentially other high molecular weight compounds. The methodology will widen the current analytical window for biomarker and food light stable isotope analyses. Moreover, our initial measurements suggest that bacterial and archaeal GDGT δ² H values can record environmental and ecological conditions.

Substrate-dependent incorporation of carbon and hydrogen for lipid biosynthesis by Methanosarcina barkeri

Dual stable isotope probing has been used to infer rates of microbial biomass production and modes of carbon fixation. In order to validate this approach for assessing archaeal production, the methanogenic archaeon Methanosarcina barkeri was grown either with H₂ , acetate or methanol with D₂ O and ¹³ C-dissolved inorganic carbon (DIC). Our results revealed unexpectedly low D incorporation into lipids, with the net fraction of water-derived hydrogen amounting to 0.357 ± 0.042, 0.226 ± 0.003 and 0.393 ± 0.029 for growth on H₂ /CO₂ , acetate and methanol respectively. The variability in net water H assimilation into lipids during the growth of M. barkeri on different substrates is possibly attributed to different Gibbs free energy yields, such that higher energy yield promoted the exchange of hydrogen between medium water and lipids. Because NADPH likely serves as the portal for H transfer, increased NADPH production and/or turnover associated with high energy yield may explain the apparent differences in net water H assimilation into lipids. The variable DIC and water H incorporation into M. barkeri lipids imply systematic, metabolic patterns of isotope incorporation and suggest that the ratio of ¹³ C-DIC versus D₂ O assimilation in environmental samples may serve as a proxy for microbial energetics in addition to microbial production and carbon assimilation pathways.