Intramolecular 13C analysis of tree rings provides multiple plant ecophysiology signals covering decades

NMR-Based Method for Intramolecular 13C Distribution at Natural Abundance Adapted to Small Amounts of Glucose

Quantitative nuclear magnetic resonance (NMR) for isotopic measurements, known as irm-NMR (isotope ratio measured by NMR), is well suited for the quantitation of ¹³C-isotopomers in position-specific isotope analysis and thus for measuring the carbon isotope composition (δ¹³C, mUr) in C-atom positions. Irm-NMR has already been used with glucose after derivatization to study sugar metabolism in plants. However, up to now, irm-NMR has exploited a "single-pulse" sequence and requires a relatively large amount of material and long experimental time, precluding many applications with biological tissues or extracts. To reduce the required amount of sample, we investigated the use of 2D-NMR analysis. We adapted and optimized the NMR sequence so as to be able to analyze a small amount (10 mg) of a glucose derivative (diacetonide glucofuranose, DAGF) with a precision better than 1 mUr at each C-atom position. We also set up a method to correct raw data and express ¹³C abundance on the usual δ¹³C scale (δ-scale). In fact, due to the distortion associated with polarization transfer and spin manipulation during 2D-NMR analyses, raw ¹³C abundance is found to be on an unusual scale. This was compensated for by a correction factor obtained via comparative analysis of a reference material (commercial DAGF) using both previous (single-pulse) and new (2D) sequences. Glucose from different biological origins (CO₂ assimilation metabolisms of plants, namely, C₃, C₄, and CAM) was analyzed with the two sequences and compared. Validation criteria such as selectivity, limit of quantification, precision, trueness, and robustness are discussed, including in the framework of green analytical chemistry.

Intramolecular carbon isotope signals reflect metabolite allocation in plants

Stable isotopes at natural abundance are key tools to study physiological processes occurring outside the temporal scope of manipulation and monitoring experiments. Whole-molecule carbon isotope ratios (13C/12C) enable assessments of plant carbon uptake yet conceal information about carbon allocation. Here, we identify an intramolecular 13C/12C signal at tree-ring glucose C-5 and C-6 and develop experimentally testable theories on its origin. More specifically, we assess the potential of processes within C3 metabolism for signal introduction based (inter alia) on constraints on signal propagation posed by metabolic networks. We propose that the intramolecular signal reports carbon allocation into major metabolic pathways in actively photosynthesizing leaf cells including the anaplerotic, shikimate, and non-mevalonate pathway. We support our theoretical framework by linking it to previously reported whole-molecule 13C/12C increases in cellulose of ozone-treated Betula pendula and a highly significant relationship between the intramolecular signal and tropospheric ozone concentration. Our theory postulates a pronounced preference for leaf cytosolic triose-phosphate isomerase to catalyse the forward reaction in vivo (dihydroxyacetone phosphate to glyceraldehyde 3-phosphate). In conclusion, intramolecular 13C/12C analysis resolves information about carbon uptake and allocation enabling more comprehensive assessments of carbon metabolism than whole-molecule 13C/12C analysis.

Metabolism is a major driver of hydrogen isotope fractionation recorded in tree-ring glucose of Pinus nigra

Stable isotope abundances convey valuable information about plant physiological processes and underlying environmental controls. Central gaps in our mechanistic understanding of hydrogen isotope abundances impede their widespread application within the plant and biogeosciences. To address these gaps, we analysed intramolecular deuterium abundances in glucose of Pinus nigra extracted from an annually resolved tree-ring series (1961-1995). We found fractionation signals (i.e. temporal variability in deuterium abundance) at glucose H¹ and H² introduced by closely related metabolic processes. Regression analysis indicates that these signals (and thus metabolism) respond to drought and atmospheric CO₂ concentration beyond a response change point. They explain ≈ 60% of the whole-molecule deuterium variability. Altered metabolism is associated with below-average yet not exceptionally low growth. We propose the signals are introduced at the leaf level by changes in sucrose-to-starch carbon partitioning and anaplerotic carbon flux into the Calvin-Benson cycle. In conclusion, metabolism can be the main driver of hydrogen isotope variation in plant glucose.

Anaplerotic flux into the Calvin-Benson cycle: hydrogen isotope evidence for in vivo occurrence in C3 metabolism

As the central carbon uptake pathway in photosynthetic cells, the Calvin-Benson cycle is among the most important biochemical cycles for life on Earth. A carbon flux of anaplerotic origin (i.e. through the chloroplast-localized oxidative branch of the pentose phosphate pathway) into the Calvin-Benson cycle was proposed recently. Here, we measured intramolecular deuterium abundances in leaf starch of Helianthus annuus grown at varying ambient CO₂ concentrations, C_a . Additionally, we modelled deuterium fractionations expected for the anaplerotic pathway and compared modelled with measured fractionations. We report deuterium fractionation signals at H¹ and H² of starch glucose. Below a C_a change point, these signals increase with decreasing C_a consistent with modelled fractionations by anaplerotic flux. Under standard conditions (C_a  = 450 ppm corresponding to intercellular CO₂ concentrations, C_i , of 328 ppm), we estimate negligible anaplerotic flux. At C_a  = 180 ppm (C_i  = 140 ppm), more than 10% of the glucose-6-phosphate entering the starch biosynthesis pathway is diverted into the anaplerotic pathway. In conclusion, we report evidence consistent with anaplerotic carbon flux into the Calvin-Benson cycle in vivo. We propose the flux may help to: maintain high levels of ribulose 1,5-bisphosphate under source-limited growth conditions to facilitate photorespiratory nitrogen assimilation required to build-up source strength; and counteract oxidative stress.

Establishment of a GC-MS-based 13 C-positional isotopomer approach suitable for investigating metabolic fluxes in plant primary metabolism

¹³ C-Metabolic flux analysis (¹³ C-MFA) has greatly contributed to our understanding of plant metabolic regulation. However, the generation of detailed in vivo flux maps remains a major challenge. Flux investigations based on nuclear magnetic resonance have resolved small networks with high accuracy. Mass spectrometry (MS) approaches have broader potential, but have hitherto been limited in their power to deduce flux information due to lack of atomic level position information. Herein we established a gas chromatography (GC) coupled to MS-based approach that provides ¹³ C-positional labelling information in glucose, malate and glutamate (Glu). A map of electron impact (EI)-mediated MS fragmentation was created and validated by ¹³ C-positionally labelled references via GC-EI-MS and GC-atmospheric pressure chemical ionization-MS technologies. The power of the approach was revealed by analysing previous ¹³ C-MFA data from leaves and guard cells, and ¹³ C-HCO₃ labelling of guard cells harvested in the dark and after the dark-to-light transition. We demonstrated that the approach is applicable to established GC-EI-MS-based ¹³ C-MFA without the need for experimental adjustment, but will benefit in the future from paired analyses by the two GC-MS platforms. We identified specific glucose carbon atoms that are preferentially labelled by photosynthesis and gluconeogenesis, and provide an approach to investigate the phosphoenolpyruvate carboxylase (PEPc)-derived ¹³ C-incorporation into malate and Glu. Our results suggest that gluconeogenesis and the PEPc-mediated CO₂ assimilation into malate are activated in a light-independent manner in guard cells. We further highlight that the fluxes from glycolysis and PEPc toward Glu are restricted by the mitochondrial thioredoxin system in illuminated leaves.

Carbon flux around leaf-cytosolic glyceraldehyde-3-phosphate dehydrogenase introduces a 13C signal in plant glucose

Within the plant and Earth sciences, stable isotope analysis is a versatile tool conveying information (inter alia) about plant physiological and paleoclimate variability across scales. Here, we identify a 13C signal (i.e. systematic 13C/12C variation) at tree-ring glucose C-4 and report an experimentally testable theory on its origin. We propose the signal is introduced by glyceraldehyde-3-phosphate dehydrogenases in the cytosol of leaves. It conveys two kinds of (potentially convoluted) information: (i) commitment of glyceraldehyde 3-phosphate to 3-phosphoglycerate versus fructose 1,6-bisphosphate metabolism; and (ii) the contribution of non-phosphorylating versus phosphorylating glyceraldehyde-3-phosphate dehydrogenase to catalysing the glyceraldehyde 3-phosphate to 3-phosphoglycerate forward reaction of glycolysis. The theory is supported by 13C fractionation modelling. Modelling results provide the first evidence in support of the cytosolic oxidation-reduction (COR) cycle, a carbon-neutral mechanism supplying NADPH at the expense of ATP and NADH, which may help to maintain leaf-cytosolic redox balances. In line with expectations related to COR cycling, we found a positive correlation between air vapour pressure deficit and 13C discrimination at glucose C-4. Overall, 13C-4 signal analysis may enable an improved understanding of leaf carbon and energy metabolism.

NMR-based isotopic and isotopomic analysis

Molecules exist in different isotopic compositions and most of the processes, physical or chemical, in living systems cause selection between heavy and light isotopes. Thus, knowing the isotopic fractionation of the common atoms, such as H, C, N, O or S, at each step during a metabolic pathway allows the construction of a unique isotope profile that reflects its past history. Having access to the isotope abundance gives valuable clues about the (bio)chemical origin of biological or synthetic molecules. Whereas the isotope ratio measured by mass spectrometry provides a global isotope composition, quantitative NMR measures isotope ratios at individual positions within a molecule. We present here the requirements and the corresponding experimental strategies to use quantitative NMR for measuring intramolecular isotope profiles. After an introduction showing the historical evolution of NMR for measuring isotope ratios, the vocabulary and symbols - for describing the isotope content and quantifying its change - are defined. Then, the theoretical framework of very accurate quantitative NMR is presented as the principle of Isotope Ratio Measurement by NMR spectroscopy, including the practical aspects with nuclei other than ²H, that have been developed and employed to date. Lastly, the most relevant applications covering three issues, tackling counterfeiting, authentication, and forensic investigation, are presented, before giving some perspectives combining technical improvements and methodological approaches.

Position-Specific Carbon Stable Isotope Ratios by Proton NMR Spectroscopy

Carbon stable isotopes provide insights into the origin and synthesis pathway of an organic molecule, and hence, contribute information that is fundamental to understanding chemical, physiological, and ecological processes. Organic carbon ¹³C/¹²C isotope ratios are commonly obtained as whole-molecule averages or as measurements of bulk samples. In contrast, position-specific isotope analysis (PSIA) provides isotope ratios for the individual carbons within a molecule, providing additional information that is masked by traditional analytical techniques. Here we introduce a ¹H NMR method for determining position-specific ¹³C/¹²C ratios within organic molecules. A peak shape superposition procedure is used to bypass the need for traditional peak integration, by exploiting relationships among the shapes of ¹H and ¹³C satellite peaks in ¹H NMR spectra. The method also has a significant sensitivity advantage over NMR methods that utilize direct detection of ¹³C. Furthermore, we demonstrate that isotope standard materials (such as those obtainable from U.S. Geological Survey) are indispensable in calibrating an NMR instrument, in order to obtain accurate isotope ratio results. Our analytical approach was applied to organic molecules of different complexity and origin, including ethanols, propionic acids, and thymidine. Results verify that chemically identical molecules from different sources can have different intramolecular isotope distributions; hence position-specific ¹³C/¹²C ratios provide an isotopic fingerprint of an organic molecule. Position-specific information for the nucleoside thymidine, where five of eight carbon positions were measured, is significant because its complexity would make it a difficult target for PSIA by mass spectrometry. The ¹H NMR method is complementary to other methods of PSIA, and will make ¹³C/¹²C PSIA employable to a wider range of organic molecules.

Tree-ring isotopes adjacent to Lake Superior reveal cold winter anomalies for the Great Lakes region of North America

Tree-ring carbon isotope discrimination (Δ¹³C) and oxygen isotopes (δ¹⁸O) collected from white pine (Pinus strobus) trees adjacent to Lake Superior show potential to produce the first winter-specific paleoclimate reconstruction with inter-annual resolution for this region. Isotopic signatures from 1976 to 2015 were strongly linked to antecedent winter minimum temperatures (T_min), Lake Superior peak ice cover, and regional to continental-scale atmospheric winter pressure variability including the North American Dipole. The immense thermal inertia of Lake Superior underlies the unique connection between winter conditions and tree-ring Δ¹³C and δ¹⁸O signals from the following growing season in trees located near the lake. By combining these signals, we demonstrate feasibility to reconstruct variability in T_min, ice cover, and continental-scale atmospheric circulation patterns (r ≥ 0.65, P < 0.001).