Intramolecular 13C signals inform about carbon 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 treering 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 photosynthesising 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 of 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. Highlight Intramolecular 13C/12C analysis resolves information about carbon uptake and allocation (and associated environmental controls) enabling more comprehensive assessments of carbon metabolism, plant-environment interactions, and environmental variability than wholemolecule 13C/12C analysis.


Introduction
Plant carbon metabolism is a central component of the global carbon cycle. It both depends on and affects environmental properties. Improved understanding of long-term plant-environment interactions relies on information from plant archives (such as tree rings) because manipulation and monitoring experiments can cover short to medium timescales only. Stable carbon isotope ( 13 C/ 12 C) analysis is among the most advanced tools to extract physiological and environmental information from plant archives. Conventionally, average 13 C/ 12 C ratios of whole plant metabolites are analysed. However, this approach neglects 13 C/ 12 C differences known to occur among individual carbon positions of plant metabolites (Abelson and Hoering, 1961). By contrast, we recently analysed intramolecular 13 C/ 12 C ratios in glucose extracted across an annually-resolved Pinus nigra tree-ring timeseries  and reported intramolecular 13 C signals (i.e., systematic 13 C/ 12 C variation confined to individual glucose carbon positions; Wieloch et al., 2018). Only after their ecophysiological origins have been elucidated, may these archived signals become useful for applications within the plant and Earth sciences.
Based on our previous dataset (Wieloch et al., 2018), we have already pinpointed a 13 C signal at tree-ring glucose C-4 and proposed it informs about carbon flux around leaf-cytosolic glyceraldehyde-3-phosphate dehydrogenases and associated energy metabolism (Wieloch, 2021;Wieloch et al., 2021). Here, we utilise the same dataset to isolate a 13 C signal at tree-ring glucose C-5 and C-6. Since intramolecular 13 C variation is governed (inter alia) by enzyme isotope effects and metabolite partitioning (Hayes, 2001), we hypothesise the signal can be linked to shifts in carbon allocation and underlying environmental controls. Thus, we develop experimentally testable theories on ecophysiological mechanisms that can introduce the signal at glucose C-5 and C-6. To this end, we consider all enzyme reactions within central carbon metabolism of C3 plants. This includes the Calvin-Benson cycle (CBC), the photosynthetic carbon oxidation (PCO) cycle, starch and sucrose synthesis and degradation, cellulose synthesis, the pentose phosphate pathway, glycolysis, and carbon metabolism downstream of phosphoenolpyruvate (PEP). Carbon exchange between other biochemical pathways and the pathway leading to the formation of tree-ring glucose are presumably small, particularly when integrated over the course of growing seasons, the timeframe of tree-ring formation. Thus, these processes cannot introduce 13 C signals of substantial size into tree-ring archives. Furthermore, we only consider primary isotope effects (which occur at atoms with altered binding after chemical reactions). Sizes of secondary isotope effects (which occur at atoms with unaltered 5 binding after chemical reactions due to indirect involvement in reaction mechanisms) are usually small and therefore unlikely to introduce detectable 13 C signals into tree-ring archives.
Finally, we present evidence supporting our theory. For this part, we reanalyse our own treering dataset (Wieloch et al., 2018) in combination with publicly accessible climate data and 13 C/ 12 C data from an ozone-treatment experiment published by Saurer et al. (1995).
We distinguish two major types of 13 C fractionation; diffusion-Rubisco fractionation, and post-Rubisco fractionation (Wieloch et al., 2018). Diffusion-Rubisco fractionation accompanies CO2 diffusion from ambient air into plant chloroplasts and subsequent carbon fixation by Rubisco (Figs. 1-2;Farquhar et al., 1982). It affects all carbon positions of plant glucose equally (Wieloch et al., 2018). By contrast, post-Rubisco fractionation results from metabolic processes downstream of Rubisco and is position-specific . Deconvolution of the two fractionation types requires the intramolecular approach.
Our work makes several conceptual advances. (i) We show how constraints on signal propagation posed by metabolic networks can be used to narrow down signal origins. (ii) A conceptual model describes how the signal propagates from its origin to other glucose carbon positions and metabolite pools. Due to space restrictions, this model is presented in SI 1. (iii) We revise current theory on plant isotope fractionation by ozone exposure. (iv) The present paper and a companion paper on the C-4 signal  develop theories that consider all relevant parts of metabolism and link intramolecular 13 C signals with specific shifts in carbon allocation and their environmental causes. For isotope signals generated within complex metabolic networks, such comprehensive theories are required as starting point for subsequent tailored experimental tests.
6 Figure 1 Central carbon metabolism in trees. Solid and dashed arrows represent substantial and negligible metabolite fluxes, respectively. Red arrows: Reactions introducing the Δ5-6' signal.

Material and methods
Intramolecular 13 C/ 12 C ratios in tree-ring glucose of Pinus nigra from Vienna (Austria) were reported in Wieloch et al. (2018). They are expressed in terms of intramolecular 13 C discrimination, Δi', where i denotes individual carbon positions in tree-ring glucose (Wieloch et al., 2018; Abbreviations and symbols in Table 1). In this notation, positive values denote discrimination against 13 C. The prime denotes measurements subjected to a procedure that removes the 13 C redistribution effect by triose phosphate cycling (SI 2; Wieloch et al., 2018). This correction restores leaf-level 13 C signals. The dataset comprises six annually resolved timeseries (one per glucose carbon) each covering the period 1961 to 1995 and containing 31 timepoints (n=6*31=186).
These authors used two different methods to determine ΔCi. Here, we calculated ΔCi averages.
For the regression analysis and linear modelling, we used 13 C/ 12 C data from Wieloch et al. (2018)

Results and discussion
Tree-ring glucose exhibits a post-Rubisco signal at C-5 and C-6 Figure 4 shows results of a Hierarchical Cluster Analysis which groups Δi' timeseries according to co-variability (Wieloch et al., 2018). That is, Δi' timeseries carrying common 13 C signals form clusters. Primary separation occurs between the Δ1' to Δ3' cluster and the Δ4' to Δ6' cluster.
Thus, these clusters convey entirely different ecophysiological information. Wieloch et al. (2018) justified the use of air vapour pressure deficits, VPD, as proxy of diffusion-Rubisco fractionation. While the average timeseries of the Δ1' to Δ3' cluster correlates highly significantly with VPD (r=-0.70, p=0.00001, n=31), the average timeseries of the Δ4' to Δ6' cluster is not significantly correlated (r=-0.30, p>0.05, n=31). This indicates that the diffusion-Rubisco signal is preserved at glucose C-1 to C-3 but not at C-4 to C-6. Among all Δi', Δ5' and Δ6' exhibit the most significant correlation (r=0.61, p≤0.001, n=31). Since the diffusion-Rubisco signal is confined to glucose C-1 to C-3, we argue C-5 and C-6 exhibit a strong post-Rubisco signal denoted Δ5-6' signal.

Exclusion of metabolic locations as origin of the Δ5-6' signal
Much is known about plant carbon metabolism. Based on this knowledge, we can exclude several metabolic locations as origin of the Δ5-6' signal as a first step in the theory development.
Note that the Δ5-6' signal is introduced at the level of three-carbon compounds because reactions at other levels do not modify carbon bonds that become glucose C-5 and C-6.
The raw dataset of intramolecular 13 C discrimination, Δi, in tree-ring glucose (Wieloch et al., 2018) exhibits significant correlations among all pairs of symmetry-related timeseries (Table   3). That is, significant correlations occur between Δ1 and Δ6, Δ2 and Δ5, and Δ3 and Δ4. These correlations likely result from carbon redistribution by triose phosphate cycling (TPC) which involves triose phosphate equilibration. Wieloch et al. (2018) describe this process mathematically and used the model to remove the TPC effect from Δi yielding a TPC-free dataset, Δi'. In this latter dataset, significant correlations among pairs of symmetry-related timeseries are absent (Table 4). This provides strong evidence for the occurrence of substantial triose phosphate equilibration in tree-ring cells of the samples discussed here.
If a process in tree-ring cells had introduced a signal at carbon positions corresponding to glucose C-5 and C-6, triose phosphate equilibration would have transmitted it to carbon positions corresponding to glucose C-2 and C-1. The signal at C-5 would have had the same size as the signal at C-2 and the signal at C-6 would have had the same size as the signal at C-1. Please note that equally-sized signals at symmetry-related glucose carbon positions are not removed by the method removing TPC effects (Wieloch et al., 2018). Since the Δ5-6' signal is absent in Δ1' and Δ2' (Fig. 4), it must have been introduced at the leaf level. Δi denotes timeseries of intramolecular 13 C discrimination (Wieloch et al., 2018). Bold numbers refer to pairs of timeseries at symmetry-related glucose carbon positions. Data were measured on tree-ring glucose of Pinus nigra laid down from 1961 to 1995 at a dry site in the Vienna basin. This table was originally published as Table 1 in Wieloch et al. (2018). Δi' denotes timeseries of intramolecular 13 C discrimination corrected for triose phosphate cycling (Wieloch et al., 2018). Bold numbers refer to pairs of timeseries at symmetry-related glucose carbon positions. Data were measured on tree-ring glucose of Pinus nigra laid down from 1961 to 1995 at a dry site in the Vienna basin. This table was originally published as Table   2 in Wieloch et al. (2018).

Exclusion of the CBC and PCO cycle as origin of the Δ5-6' signal
Introduction of the Δ5-6' signal within the CBC or PCO cycle can be excluded because hexose phosphate synthesis includes conversion of photosynthetic/photorespiratory glyceraldehyde 3- Figs. 1-2). This would transmit any 13 C signal present at GAP carbon positions corresponding to glucose C-5 and C-6 to DHAP carbon positions corresponding to glucose C-2 and C-1. More generally, metabolites feeding into the stromal GAP pool can be excluded as origin of the Δ5-6' signal based on the same reasoning.

Exclusion of reactions downstream of OAA, pyruvate, and DAHP as Δ5-6' signal origin
Pyruvate kinase (PK) and pyruvate orthophosphate dikinase (PPDK) interconvert PEP and pyruvate (Figs. 1, and 3). The PK reaction is strongly on the side of pyruvate and considered nearly irreversible (Nageswara Rao et al., 1979;Tcherkez et al., 2011). In illuminated leaves of C3 plants, PPDK activity is either very low or undetectable, except for orchids, and grasses (Hocking and Anderson, 1986). In illuminated leaves of Xanthium strumarium, flux from pyruvate to PEP is very small at ≈0.05% of net CO2 assimilation (Tcherkez et al., 2011). In Arabidopsis thaliana, PPDK is upregulated during leaf senescence which is believed to facilitate nitrogen remobilisation (Taylor et al., 2010). This is of minor importance here because leaf senescence occurs during a short period relative to the multi-year lifespan of conifer needles. In Nicotiana tabacum, PPDK activity is increased up to 2.7-fold under strong drought (Doubnerová Hýsková et al., 2014). However, this should not result in significant flux in relation to fluxes in carbohydrate metabolism since basal PPDK activities in C3 plants are generally low (Hocking and Anderson, 1986;Tcherkez et al., 2011). Thus, flux from pyruvate to PEP should be small, and transmission of 13 C signals in pyruvate to cytosolic carbohydrates by the PK/PPDK interface should be negligible.
Phosphoenolpyruvate carboxylase (PEPC) and phosphoenolpyruvate carboxykinase (PEPCK) interconvert PEP and oxaloacetate (OAA). Conversion of PEP to OAA by PEPC is irreversible (Chollet et al., 1996). To our knowledge, there are no reports of PEPCK activity in mesophyll cells where bulk carbohydrate synthesis takes place (Pyke, 2001). PEPCK RNA and protein were not detected in leaves of Solanum lycopersicum irrespective of their developmental stage (Bahrami et al., 2001;Famiani et al., 2016). PEPCK protein or activity were not detected in leaves of Hordeum vulgare (Chen et al., 2000). In mature leaves of Arabidopsis thaliana, PEPCK protein amount and activity are low and likely confined to specific cell types (Malone et al., 2007). In leaves of Nicotiana tabacum, PEPCK occurs in trichomes and stomata (Leegood et al., 1999;Malone et al., 2007). In leaves of Cucumis sativus, PEPCK occurs in trichomes and phloem cells (Leegood et al., 1999;Chen et al., 2004). In leaves of Oryza sativa, PEPCK occurs in hydathodes, stomata, and the vascular parenchyma (Bailey and Leegood, 2016

Exclusion of starch, sucrose, and cellulose metabolism, and the pentose phosphate pathway as origin of the Δ5-6' signal
Reactions leading directly from stromal GAP to the formation of starch, sucrose, and cellulose, reactions remobilising starch, and reactions of the pentose phosphate pathway do not simultaneously modify carbon bonds that become glucose C-5 and C-6. This excludes these pathways for Δ5-6' signal introduction.

Origin of the Δ5-6' signal
After excluding several metabolic locations as the origin of the Δ5-6' signal (see 'Exclusion of metabolic locations as origin of the Δ5-6' signal'), the glycolytic pathway and PEP metabolism in leaves are left for consideration. Within this system, cytosolic GAP is used for PEP metabolism and sucrose synthesis (Fig. 1). Thus, GAP constitutes a central branch point in leaf carbon metabolism enabling isotope fractionation.

Leaf-level enolase, PEPC, PK, and/or DAHPS introduce the Δ5-6' signal
Enolase, PEPC, PK, and DAHPS are the only enzymes which simultaneously modify carbon bonds that become glucose C-5 and C-6 ( Figs. 1-3). Reactions catalysed by these enzymes may be accompanied by 13 C effects of substantial size and may thus introduce the Δ5-6' signal.
Enolase interconverts 2-phosphoglycerate (2PGA) and PEP (Figs. 1, and 3). In vivo, the reaction operates close to equilibrium (Kubota and Ashihara, 1990) and might thus be accompanied by an equilibrium isotope effect. Formation of the C=C double bond in PEP likely favours turnover of 13 C isotopologues of 2PGA leading to 13 C enrichment in PEP. A 13 C signal might then arise from varying allocation of PEP to downstream processes (Fig. 1). Increased downstream consumption would remove more 13 C-enriched PEP and leave behind more 12 Cenriched 2PGA for glucose synthesis.
Kinetic isotope effects may accompany the unidirectional conversions of PEP to OAA by PEPC, PEP to pyruvate by PK, and PEP and E4P to DAHP by DAHPS (Figs. 1, and 3). These reactions break the C=C double bond in PEP and can therefore be expected to favour turnover of 12 C-isotopologues of PEP leaving behind 13 C-enriched PEP for glucose synthesis. Due to the usually larger size of kinetic isotope effects compared to equilibrium isotope effects, effects by PEPC, PK and DAHPS can be expected to outweigh any reciprocal effect by enolase. Thus, considering all four enzymes together, increasing turnover of PEP by PEPC, PK, and DAHPS can be expected to result in 13 C-enriched tree-ring glucose, i.e., Δ5-6' decreases.

Signal transmission to tree-ring glucose
Isotope signals generated at the level of leaf-cytosolic PEP or 2PGA need to be transmitted to GAP to then enter hexose phosphates and tree-ring glucose (Fig. 1). Transmission of a signal introduced by cytosolic enzymes is straightforward since the cytosolic glycolytic reactions between PEP and GAP are at equilibrium (Kubota and Ashihara, 1990). Transmission of a signal introduced by stromal enzymes is more intricate. First, it requires an incomplete or lowactivity glycolytic pathway in leaf chloroplasts because signal equilibration with stromal triose phosphates would result in even signal distribution over all glucose carbon positions (SI 1.6).
An incomplete glycolytic pathway in leaf chloroplasts is supported by a reported lack of enolase in Arabidopsis thaliana and Oryza sativa (van der Straeten et al., 1991;Prabhakar et al., 2009;Fukayama et al., 2015). Second, signal transmission from stromal PEP to C-5 and C-6 of cytosolic hexose phosphate requires chloroplast export of PEP. Transport of PEP across the chloroplasts' inner membrane is mediated by the PEP/Pi translocator as counter-exchange with Pi, PEP, or 2PGA and the putative in vivo preference for the transport of Pi, and PEP (Fischer et al., 1997;Flügge et al., 2011). Numerous stromal processes, such as the shikimate pathway and fatty acid biosynthesis, rely on PEP import from the cytosol (Streatfield et al., 1999;Flügge et al., 2011). Therefore, a net flux of PEP from the cytosol to chloroplasts can be expected.
However, members of the phosphate translocator family are believed to be highly inefficient.
For instance, merely 10% of the activity of the triose phosphate translocator is used for net export of triose phosphate from chloroplasts; 90% is wasted on futile counter-exchanges (Flügge, 1987(Flügge, , 1999. In addition, low stromal and high cytosolic Pi levels (Sharkey and Vanderveer, 1989) can be expected to promote chloroplast export of PEP. Consequently, efficient equilibration of cytosolic and stromal PEP pools, and 13 C signals therein, can be expected. Thus, both cytosolic and stromal enzymes may contribute to the Δ5-6' signal.

Signal introduction requires substantial carbon fluxes and flux variability
For the introduction of a 13 C signal, a substantial share of the photosynthetically fixed carbon must be directed towards enolase, PEPC, PK, and/or DAHPS and their downstream derivatives.
This share must vary substantially; in the present case, on the interannual timescale. Therefore, we will now discuss carbon fluxes through enolase, PEPC, PK, and DAHPS.  (Dizengremel, 2001). In forest trees, the Rubisco/PEPC activity ratio can change from up to 25 in ozone-free air to about 2 under realistic levels of ambient ozone redirecting carbon flux to maintenance and repair processes (Dizengremel, 2001).
Isoforms of PK are localised in both the cytosol and chloroplasts (Figs. 1, and 3; Ambasht and Kayastha, 2002). They provide pyruvate (inter alia) for mitochondrial respiration, fatty acid 20 biosynthesis, and the non-mevalonate pathway (MEP). To our knowledge, estimates of the respiratory flux via PK in actively photosynthesising leaves are unavailable. However, this flux may be substantial when photorespiration is low and thus co-vary with photorespiration and its environmental controls (SI 3).
In illuminated photosynthetic tissue of Arabidopsis thaliana, fatty acid biosynthesis can occur at a rate of 2.3 µmol C mg chl -1 h -1 (Bao et al., 2000). Based on this, we estimate a ≈2% carbon flux relative to net CO2 assimilation into fatty acid biosynthesis (SI 4). In leaves, this flux is predominantly controlled at the level of acetyl-CoA carboxylase (Page et al., 1994;Harwood, 2005;Ohlrogge et al., 2015). It responds to the stromal redox state (energy status) and associated environmental controls (Rawsthorne, 2002;Harwood, 2005;Geigenberger and Fernie, 2014).
The plastid-localised MEP pathway is yet another metabolic route carrying substantial flux.
Since PEPC, PK, and DAHPS are located downstream of enolase (Figs. 1, and 3), all four enzymes may contribute to the Δ5-6' signal. Based on arguments given above, associated carbon fluxes and their variability can be expected to be substantial. Other leaf-level pathways consuming PEP, such as the cytosolic mevalonate pathway, may exert additional control over the Δ5-6' signal.

22
The study site is ≈10 km away from the city centre of Vienna and frequently exposed to substantial levels of tropospheric ozone (Oltmans et al., 1998;Ainsworth et al., 2012). Lefohn (1992) classified Pinus nigra as an ozone-sensitive tree species. Radiation stimulates the photochemical reactions of ozone formation (Ainsworth et al., 2012). Ozone triggers relative flux increases through the anaplerotic and shikimate pathways via PEPC and DAHPS, respectively (see 'Signal introduction requires substantial carbon fluxes and flux variability') and may thus cause 13 C increases at PEP carbon positions that become glucose C-5 and C-6 (see 'Leaf-level enolase, PEPC, PK, and/or DAHPS introduce the Δ5-6' signal'). This may introduce an isotope signal independent of the diffusion-Rubisco signal due to independence at the level of environmental controls.
By contrast, a process mitigating ozone entry into plant leaves may explain the component of the Δ5-6' signal which is inversely correlated with the diffusion-Rubisco signal. Dizengremel (2001) proposed that drought and ozone combined is a main recurring stress factor in forest ecosystems. Isohydric plant species, such as Pinus nigra, respond to drought by closing their stomata (Sade et al., 2012). Reduced stomatal conductance impedes ozone uptake (Tingey and Hogsett, 1985;Dobson et al., 1990;Dizengremel, 2001). In needles of Pinus halepensis, PEPC activities in control plants and plants exposed to mild drought stress were similar, strongly increased under ozone stress, but significantly less so under combined ozone and drought stress (Fontaine et al., 2003). Thus, anaplerotic flux rates can be expected to be highest under ozone stress but lower when ozone stress is accompanied by drought. While drought causes 13 C enrichments at all glucose carbon positions due to diffusion-Rubisco fractionation (Wieloch et al., 2018), it can be expected to reduce ozone-induced 13 C enrichments at glucose C-5 and C-6. This drought component of the ozone response may remove the diffusion-Rubisco signal from glucose C-5 and C-6. In SI 3, we discuss how changes in substrate supply to mitochondrial oxidative phosphorylation (glycolytic pyruvate versus photorespiratory glycine) may additionally contribute to the component of the Δ5-6' signal that is inversely correlated with diffusion-Rubisco fractionation.

Effects of tropospheric ozone on whole-molecule 13 C/ 12 C composition of plant cellulose
Growing Betula pendula at increased ozone levels, several authors reported decreased 13 C discrimination, Δ, in leaf and stem cellulose (Matyssek et al., 1992;Saurer et al., 1995).
Intriguingly, these Δ decreases coincided with increased ratios of intercellular to ambient CO2 concentrations, Ci/Ca. As pointed out by Matyssek et al. (1992) and Saurer et al. (1995), this cannot be explained by the standard model of diffusion-Rubisco fractionation which predicts a positive correlation between Ci/Ca and Δ (Farquhar et al., 1982). Thus, post-Rubisco fractionation can be expected to cause these ozone-related isotope effects. Matyssek et al. (1992) and Saurer et al. (1995) proposed increased relative carbon fixation by PEPC due to ozone explains the Δ decreases because carbon fixed by PEPC is strongly 13 C enriched compared to carbon fixed by Rubisco (Melzer and O'Leary, 1987). While this proposal is in line with significantly increased relative PEPC activities observed under ozone (Table 2), it conflicts with the setup of carbon metabolism. PEPC-fixed carbon supplies downstream metabolism, yet no pathway carrying substantial flux exists that could transfer it into carbohydrate metabolism (see 'Exclusion of reactions downstream of OAA, pyruvate, and DAHP as Δ5-6' signal origin'). Above, we propose an ozone-dependent mechanism for the introduction of the Δ5-6' signal which reconciles observations of Matyssek et al. (1992) and Saurer et al. (1995) with the setup of carbon metabolism (see 'Ecophysiological effects').
In Betula pendula, post-Rubisco fractionation causes average whole-molecule ΔΔ decreases of ≈-1.98 ± 0.58SE ‰ (Fig. 5). Below, we propose a fraction of the Δ5-6' signal enters glucose C-1 to C-4 through indirect signal propagation via chloroplast metabolism (see 'Signal propagation to all glucose carbons via chloroplast metabolism'). We estimate the signal at C-5 and C-6 is 6.625-fold larger than at C-1 to C-4 (SI 1.8). Thus, a ≈-1.98 ± 0.58SE ‰ effect at the whole-molecule level scales to ≈-4.56 ± 1.34SE ‰ effects at cellulose glucose C-5 and C-24 6 and to ≈-0.69 ± 0.20SE ‰ effects at C-1 to C-4 (SI 6). In Pinus nigra, measured Δ5-6' values fall within a 5.80 ± 1.55SE ‰ range (maximum=22.71 ± 0.99SE ‰, minimum=16.91 ± 0.56SE ‰). Wieloch et al. (2018) estimated the Δ5-6' timeseries contains 79% systematic and 21% error variance. Assuming the error is fully expressed in both the maximum and minimum value, we estimate a systematic timeseries range of ≈4.58 ± 1.22SE ‰ (5.80 ± 1.55SE ‰ * 0.79). This largely equals the estimated effect at glucose C-5 and C-6 in ozone-treated Betula pendula corroborating the theory proposed above. Notably, occurrence of the post-Rubisco fractionation effect in leaf cellulose of Betula pendula corroborates the proposed leaf-level origin of the Δ5-6' signal. values were estimated using a model by Farquhar et al. (1982). This model describes 13 C discrimination associated with plant carbon uptake including CO2 diffusion into plant leaves and assimilation by Rubisco. Numbers inside bars denote differences between measured and expected ΔΔ values. Statistically significant differences are marked by asterisks (one-tailed ttest: *, p<0.05; **, p<0.01). The difference of the Leaf/HF treatment is close to being statistically significant (p<0.08). This analysis is based on data published by Saurer et al. (1995).

Signal propagation at the level of TPI in the cytosol of leaves
The post-Rubisco signal at glucose C-5 and C-6 is independent of a signal at glucose C-1 and C-2 ( Fig. 4; Wieloch et al., 2018). That is, substantial signal propagation from C-5 and C-6 to C-2 and C-1 is not supported by the data. This is surprising for the following reason.
Transmission of the Δ5-6' signal from its origin, the lower end of the glycolytic pathway, into carbohydrate metabolism occurs via GAP (Fig. 1). Hence, signal independency requires negligible conversion of GAP (a precursor of glucose C-4 to C-6) to DHAP (a precursor of glucose C-1 to C-3) via leaf-cytosolic TPI (Figs. 1-2). Since TPI is often referred to as the prime example for the efficiency of enzyme catalysis, one would expect full equilibration of GAP and DHAP and inherent isotope signals. This view, however, is based on in vitro measurements of TPI kinetics. The following mechanisms may explain the apparent lack of equilibration and signal propagation in vivo.
In the light, chloroplast export of DHAP is favoured by the equilibrium position of stromal TPI, which is strongly on the side of DHAP (Walker, 1976;Knowles and Albery, 1977). Sharkey and Weise (2012) calculated there should be over 20 times more DHAP than GAP at equilibrium (Meyerhof and Junowicz-Kocholaty, 1943;Bassham and Krause, 1969). The substrate affinities of the triose phosphate translocator for DHAP and GAP are similar at Km=0.13 mM and Km=0.08 mM, respectively (Fliege et al., 1978). Thus, DHAP and GAP will be transported according to their concentrations. That is, 20 times more DHAP will be exported from chloroplasts to the cytosol. However, synthesis of fructose 1,6-bisphosphate uses DHAP and GAP at a 1:1 ratio. This may keep the concentration of leaf-cytosolic GAP low. Flux of GAP into glycolysis and processes consuming glycolytic intermediates will additionally contribute to low cytosolic GAP concentrations. This may restrict the GAP to DHAP back conversions. Furthermore, numerous common metabolites inhibit TPI competitively (Anderson, 1971;Grüning et al., 2014;Flügel et al., 2017;Li et al., 2019). In addition, the activity of cytosolic TPI decreases significantly upon treatment with reactive oxygen species, especially H2O2 (Lopez-Castillo et al., 2016). Thus, during active photosynthesis, a lack of isomeric and isotopic equilibrium between leaf-cytosolic GAP and DHAP is conceivable. This would block the propagation of 13 C signals in GAP to DHAP and enable independent 13 C signals in Δ5-6' and Δ1-2' as observed.

Signal propagation to other plant metabolites
Carbon flux changes around leaf-cytosolic enolase, PEPC, PK, and DAHPS proposedly introduce the Δ5-6' signal. Hence, upstream derivatives of 2PGA carbons corresponding to glucose C-5 and C-6 will inherit the signal (Fig. 1). Compared to plant cellulose, the signal will be distinctly smaller in chloroplast derivatives, such as starch, and distinctly larger in leaf sucrose synthesised during the photoperiod (SI 1). This latter metabolite may be used to follow the Δ5-6' signal on an hourly basis. Downstream derivatives of PEP carbons corresponding to glucose C-5 and C-6 will obtain an inverse Δ5-6' signal. These differences may help to test the theory.

Implications for whole-molecule 13 C/ 12 C analysis
The Δ5-6' signal has two components of variance (see 'Ecophysiological effects'). One is inversely correlated with diffusion-Rubisco fractionation, the other is independent. Both components have implications for studies of plant carbon uptake and associated properties by whole-molecule 13 C/ 12 C analysis. The inversely correlated component removes the diffusion-Rubisco signal from glucose C-5 and C-6. In addition, this signal is absent at glucose C-4 (Wieloch et al., 2018). That is, three out of six glucose carbon positions lack the diffusion-29 Rubisco signal. Thus, whole-molecule 13 C/ 12 C analysis captures an attenuated diffusion-Rubisco signal and underestimates the variability of the original signal and associated physiological properties, such as Ci/Ca and photosynthetic water-use efficiency.
The independent component of the Δ5-6' signal weakens signal extractions from whole-molecule 13 C/ 12 C measurements because it constitutes pseudorandom noise with respect to diffusion-Rubisco fractionation. This may explain why models of whole-molecule diffusion-Rubisco fractionation as functions of environmental properties often suffer from low explanatory powers, R 2 ≤0.5 (Barbour and Song, 2014). By contrast, intramolecular 13 C/ 12 C analysis resolves information about distinct ecophysiological processes; a fundamental conceptual advancement enabling more adequate modelling of the variability of plant carbon uptake and associated environmental/developmental controls.

Tracking carbon allocation in other biological organisms
Whole-molecule 13 C/ 12 C analysis enables assessments of plant carbon uptake (Farquhar et al., 1982). According to theory reported here, intramolecular 13 C/ 12 C analysis enables additional assessment of downstream carbon allocation in actively photosynthesising leaves. This includes carbon flux into the anaplerotic, shikimate, MEP, and fatty acid biosynthesis pathways, and mitochondrial respiration (Fig. 1). Intramolecular 13 C signals are governed by a small set of physicochemical principles that apply generally (Schmidt et al., 2015). Thus, intramolecular 13 C/ 12 C analysis can be expected to enable retrospective assessment of carbon allocation in any biological organism including, for instance, disease-related shifts.

Utility of the Δ5-6' signal
Laboratory experiments are limited to short timescales and in their capabilities to reproduce complex natural systems. Manipulation experiments on natural systems are limited to timescales of years and may suffer from spurious effects due to unnatural step changes in ambient conditions. By contrast, tree-ring analysis offers extensive temporal, spatial, species, and genotype coverage of natural systems that have not been subjected to unnatural step changes.
Thus, signal analysis may enable a better understanding of plant and ecosystem carbon balances including the so-called CO2 fertilisation effect.
Furthermore, intramolecular 13 C/ 12 C analysis not only enables analysis of carbon uptakeenvironment relationships but also of carbon allocation-environment relationships  and thus more comprehensive assessments of flux-level plant performance. For instance, atmospheric CO2 and ozone concentrations have increased over recent years (Fig. 6b). Under business-as-usual scenarios, this will continue over the next decades (Turnock et al., 2018).
While CO2 promotes leaf photosynthesis and net primary productivity, ozone has the very opposite effect (Ainsworth et al., 2012;IPCC, 2014). When this highly reactive chemical enters plant leaves through stomata, it causes harm to structure and function and leads to major rearrangements in carbon metabolism. While ozone decreases carbon fixation, it increases carbon allocation to costly maintenance and repair processes (Dizengremel, 2001;Ainsworth et al., 2012). This includes increased carbon flux into the anaplerotic and shikimate pathway, and this resource investment is likely recorded in the Δ5-6' signal. Ozone tolerance varies among species with metabolic changes depending on the duration of ozone exposure (Fontaine et al., 2003;Ainsworth et al., 2012). Thus, the Δ5-6' signal may support flux-level screenings for species/genotypes with the capacity to optimally adjust to prolonged ozone exposure (requires further investigation).
While glucose positions C-1 to C-3 preserve the VPD-dependent carbon uptake signal (Wieloch et al., 2018), this signal was removed from glucose C-5 and C-6 and replaced by an independent ozone-sensitive carbon allocation signal (Figs. 5-6). Thus, intramolecular 13 C/ 12 C analysis yields information about several environmental variables and may enable more powerful paleoenvironment reconstructions than whole-molecule analysis.
Lastly, sampling glucose at different developmental stages may enable the detection of shifts in carbon uptake and allocation over the lifespan of plants to better understand basic physiological processes such as plant senescence. In conclusion, intramolecular 13 C/ 12 C analysis opens up numerous new avenues of research within the plant and Earth sciences.