Isotopically nonstationary 13C flux analysis of changes in Arabidopsis thaliana leaf metabolism due to high light acclimation

Leaf day respiration: More than just catabolic CO2 production in the light

Day respiration is a net flux resulting from several CO2‐generating and CO2‐fixing reactions, not only related to catabolism but also to anabolism. We review pieces of evidence that decarboxylating reactions are partly fed by carbon sources disconnected from current photosynthesis and how they reflect various metabolic pathways.

Central transcriptional regulator controls photosynthetic growth and carbon storage in response to high light

Carbon capture and biochemical storage are some of the primary drivers of photosynthetic yield and productivity. To elucidate the mechanisms governing carbon allocation, we designed a photosynthetic light response test system for genetic and metabolic carbon assimilation tracking, using microalgae as simplified plant models. The systems biology mapping of high light-responsive photophysiology and carbon utilization dynamics between two variants of the same Picochlorum celeri species, TG1 and TG2 elucidated metabolic bottlenecks and transport rates of intermediates using instationary 13C-fluxomics. Simultaneous global gene expression dynamics showed 73% of the annotated genes responding within one hour, elucidating a singular, diel-responsive transcription factor, closely related to the CCA1/LHY clock genes in plants, with significantly altered expression in TG2. Transgenic P. celeri TG1 cells expressing the TG2 CCA1/LHY gene, showed 15% increase in growth rates and 25% increase in storage carbohydrate content, supporting a coordinating regulatory function for a single transcription factor.

The oxidative pentose phosphate pathway in photosynthesis: a tale of two shunts

CO2 release in the light (RL) and its presumed source, oxidative pentose phosphate pathways, were found to be insensitive to CO2 concentration. The oxidative pentose phosphate pathways form glucose 6-phosphate (G6P) shunts that bypass the nonoxidative pentose phosphate reactions of the Calvin-Benson cycle. Using adenosine diphosphate glucose and uridine diphosphate glucose as proxies for labeling of G6P in the stroma and cytosol respectively, it was found that only the cytosolic shunt was active. Uridine diphosphate glucose, a proxy for cytosolic G6P, and 6-phosphogluconate (6PG) were significantly less labeled than Calvin-Benson cycle intermediates in the light. But ADP glucose, a proxy for stromal G6P, is labeled to the same degree as Calvin-Benson cycle intermediates and much greater than 6PG. A metabolically inert pool of sedoheptulose bisphosphate can slowly equilibrate keeping the label in sedoheptulose lower than in other stromal metabolites. Finally, phosphorylation of fructose 6-phosphate (F6P) in the cytosol can allow some unlabeled carbon in cytosolic F6P to dilute label in phosphenolpyruvate. The results clearly show that there is oxidative pentose phosphate pathway activity in the cytosol that provides a shunt around the nonoxidative pentose phosphate pathway reactions of the Calvin-Benson cycle and is not strongly CO2-sensitive.

Resource allocation modeling for autonomous prediction of plant cell phenotypes

Predicting the plant cell response in complex environmental conditions is a challenge in plant biology. Here we developed a resource allocation model of cellular and molecular scale for the leaf photosynthetic cell of Arabidopsis thaliana, based on the Resource Balance Analysis (RBA) constraint-based modeling framework. The RBA model contains the metabolic network and the major macromolecular processes involved in the plant cell growth and survival and localized in cellular compartments. We simulated the model for varying environmental conditions of temperature, irradiance, partial pressure of CO2 and O2, and compared RBA predictions to known resource distributions and quantitative phenotypic traits such as the relative growth rate, the C:N ratio, and finally to the empirical characteristics of CO2 fixation given by the well-established Farquhar model. In comparison to other standard constraint-based modeling methods like Flux Balance Analysis, the RBA model makes accurate quantitative predictions without the need for empirical constraints. Altogether, we show that RBA significantly improves the autonomous prediction of plant cell phenotypes in complex environmental conditions, and provides mechanistic links between the genotype and the phenotype of the plant cell.

Modeling lignin biosynthesis: a pathway to renewable chemicals

Plant biomass contains lignin that can be converted into high-value-added chemicals, fuels, and materials. The precise genetic manipulation of lignin content and composition in plant cells offers substantial environmental and economic benefits. However, the intricate regulatory mechanisms governing lignin formation challenge the development of crops with specific lignin profiles. Mathematical models and computational simulations have recently been employed to gain fundamental understanding of the metabolism of lignin and related phenolic compounds. This review article discusses the strategies used for modeling plant metabolic networks, focusing on the application of mathematical modeling for flux network analysis in monolignol biosynthesis. Furthermore, we highlight how current challenges might be overcome to optimize the use of metabolic modeling approaches for developing lignin-engineered plants.

The end game(s) of photosynthetic carbon metabolism

The year 2024 marks 70 years since the general outline of the carbon pathway in photosynthesis was published. Although several alternative pathways are now known, it is remarkable how many organisms use the reaction sequence described 70 yrs ago, which is now known as the Calvin-Benson cycle or variants such as the Calvin-Benson-Bassham cycle or Benson-Calvin cycle. However, once the carbon has entered the Calvin-Benson cycle and is converted to a 3-carbon sugar, it has many potential fates. This review will examine the last stages of photosynthetic metabolism in leaves. In land plants, this process mostly involves the production of sucrose provided by an endosymbiont (the chloroplast) to its host for use and transport to the rest of the plant. Photosynthetic metabolism also usually involves the synthesis of starch, which helps maintain respiration in the dark and enables the symbiont to supply sugars during both the day and night. Other end products made in the chloroplast are closely tied to photosynthetic CO2 assimilation. These include serine from photorespiration and various amino acids, fatty acids, isoprenoids, and shikimate pathway products. I also describe 2 pathways that can short circuit parts of the Calvin-Benson cycle. These final processes of photosynthetic metabolism play many important roles in plants.

UMAMIT44 is a key player in glutamate export from Arabidopsis chloroplasts

Selective partitioning of amino acids among organelles, cells, tissues, and organs is essential for cellular metabolism and plant growth. Nitrogen assimilation into glutamine and glutamate and de novo biosynthesis of most protein amino acids occur in chloroplasts; therefore, various transport mechanisms must exist to accommodate their directional efflux from the stroma to the cytosol and feed the amino acids into the extraplastidial metabolic and long-distance transport pathways. Yet, Arabidopsis (Arabidopsis thaliana) transporters functioning in plastidial export of amino acids remained undiscovered. Here, USUALLY MULTIPLE ACIDS MOVE IN AND OUT TRANSPORTER 44 (UMAMIT44) was identified and shown to function in glutamate export from Arabidopsis chloroplasts. UMAMIT44 controls glutamate homeostasis within and outside of chloroplasts and influences nitrogen partitioning from leaves to sinks. Glutamate imbalances in chloroplasts and leaves of umamit44 mutants impact cellular redox state, nitrogen and carbon metabolism, and amino acid (AA) and sucrose supply of growing sinks, leading to negative effects on plant growth. Nonetheless, the mutant lines adjust to some extent by upregulating alternative pathways for glutamate synthesis outside the plastids and by mitigating oxidative stress through the production of other amino acids and antioxidants. Overall, this study establishes that the role of UMAMIT44 in glutamate export from chloroplasts is vital for controlling nitrogen availability within source leaf cells and for sink nutrition, with an impact on growth and seed yield.

Mechanisms controlling metabolite concentrations of the Calvin Benson Cycle

Maintaining proper metabolite levels in a complex metabolic network is crucial for maintaining a high flux through the network. In this paper, we discuss major regulatory mechanisms over the Calvin Benson Cycle (CBC) with regard to their roles in conferring homeostasis of metabolite levels in CBC. These include: 1) Redox regulation of enzymes in the CBC on one hand ensures that metabolite levels stay above certain lower bounds under low light while on the other hand increases the flux through the CBC under high light. 2) Metabolite regulations, especially allosteric regulations of major regulatory enzymes, ensure the rapid up-regulation of fluxes to ensure sufficient amount of triose phosphate is available for end product synthesis and concurrently avoid phosphate limitation. 3) A balanced activities of enzymes in the CBC help maintain balanced flux through CBC; some innate product feedback mechanisms, in particular the ADP feedback regulation of GAPDH and F6P feedback regulation of FBPase, exist in CBC to achieve such a balanced enzyme activities and hence flux distribution in the CBC for greater photosynthetic efficiency. Transcriptional regulation and natural variations of enzymes controlling CBC metabolite homeostasis should be further explored to maximize the potential of engineering CBC for greater efficiency.

Natural variation in metabolism of the Calvin-Benson cycle

The Calvin-Benson cycle (CBC) evolved over 2 billion years ago but has been subject to massive selection due to falling atmospheric carbon dioxide, rising atmospheric oxygen and changing nutrient and water availability. In addition, large groups of organisms have evolved carbon-concentrating mechanisms (CCMs) that operate upstream of the CBC. Most previous studies of CBC diversity focused on Rubisco kinetics and regulation. Quantitative metabolite profiling provides a top-down strategy to uncover inter-species diversity in CBC operation. CBC profiles were recently published for twenty species including terrestrial C3 species, terrestrial C4 species that operate a biochemical CCM, and cyanobacteria and green algae that operate different types of biophysical CCM. Distinctive profiles were found for species with different modes of photosynthesis, revealing that evolution of the various CCMs was accompanied by co-evolution of the CBC. Diversity was also found between species that share the same mode of photosynthesis, reflecting lineage-dependent diversity of the CBC. Connectivity analysis uncovers constraints due to pathway and thermodynamic topology, and reveals that cross-species diversity in the CBC is driven by changes in the balance between regulated enzymes and in the balance between the CBC and the light reactions or end-product synthesis.

SIMPEL: using stable isotopes to elucidate dynamics of context specific metabolism

The capacity to leverage high resolution mass spectrometry (HRMS) with transient isotope labeling experiments is an untapped opportunity to derive insights on context-specific metabolism, that is difficult to assess quantitatively. Tools are needed to comprehensively mine isotopologue information in an automated, high-throughput way without errors. We describe a tool, Stable Isotope-assisted Metabolomics for Pathway Elucidation (SIMPEL), to simplify analysis and interpretation of isotope-enriched HRMS datasets. The efficacy of SIMPEL is demonstrated through examples of central carbon and lipid metabolism. In the first description, a dual-isotope labeling experiment is paired with SIMPEL and isotopically nonstationary metabolic flux analysis (INST-MFA) to resolve fluxes in central metabolism that would be otherwise challenging to quantify. In the second example, SIMPEL was paired with HRMS-based lipidomics data to describe lipid metabolism based on a single labeling experiment. Available as an R package, SIMPEL extends metabolomics analyses to include isotopologue signatures necessary to quantify metabolic flux.

Model validation and selection in metabolic flux analysis and flux balance analysis

13C-Metabolic Flux Analysis (13C-MFA) and Flux Balance Analysis (FBA) are widely used to investigate the operation of biochemical networks in both biological and biotechnological research. Both methods use metabolic reaction network models of metabolism operating at steady state so that reaction rates (fluxes) and the levels of metabolic intermediates are constrained to be invariant. They provide estimated (MFA) or predicted (FBA) values of the fluxes through the network in vivo, which cannot be measured directly. These fluxes can shed light on basic biology and have been successfully used to inform metabolic engineering strategies. Several approaches have been taken to test the reliability of estimates and predictions from constraint-based methods and to compare alternative model architectures. Despite advances in other areas of the statistical evaluation of metabolic models, such as the quantification of flux estimate uncertainty, validation and model selection methods have been underappreciated and underexplored. We review the history and state-of-the-art in constraint-based metabolic model validation and model selection. Applications and limitations of the χ2 -test of goodness-of-fit, the most widely used quantitative validation and selection approach in 13C-MFA, are discussed, and complementary and alternative forms of validation and selection are proposed. A combined model validation and selection framework for 13C-MFA incorporating metabolite pool size information that leverages new developments in the field is presented and advocated for. Finally, we discuss how adopting robust validation and selection procedures can enhance confidence in constraint-based modeling as a whole and ultimately facilitate more widespread use of FBA in biotechnology.

Tools for Measuring Photosynthesis at Different Scales

Measurements of in vivo photosynthesis are powerful tools that probe the largest fluxes of carbon and energy in an illuminated leaf, but often the specific techniques used are so varied and specialized that it is difficult for researchers outside the field to select and perform the most useful assays for their research questions. The goal of this chapter is to provide a broad overview of the current tools available for the study of photosynthesis, both in vivo and in vitro, so as to provide a foundation for selecting appropriate techniques, many of which are presented in detail in subsequent chapters. This chapter will also organize current methods into a comparative framework and provide examples of how they have been applied to research questions of broad agronomical, ecological, or biological importance. This chapter closes with an argument that the future of in vivo measurements of photosynthesis lies in the ability to use multiple methods simultaneously and discusses the benefits of this approach to currently open physiological questions. This chapter, combined with the relevant methods chapters, could serve as a laboratory course in methods in photosynthesis research or as part of a more comprehensive laboratory course in general plant physiology methods.

Combining Gas Exchange and Rapid Quenching of Leaf Tissue for Mass Spectrometry Analysis Directly in Gas Exchange Cuvette

Isotopically nonstationary metabolic flux analysis (INST-MFA) is a powerful technique for studying plant central metabolism, which involves introducing a 13CO2 tracer to plant leaves and sampling the labeled metabolic intermediates during the transient period before reaching an isotopic steady state. The metabolic intermediates involved in the C3 cycle have exceptionally fast turnover rates, with some intermediates turning over many times a second. As a result, it is necessary to rapidly introduce the label and then rapidly quench the plant tissue to determine concentrations in the light or capture the labeling kinetics of these intermediates at early labeling time points. Here, we describe a rapid quenching (0.1-0.5 s) system for 13CO2 labeling experiments in plant leaves to minimize metabolic changes during labeling and quenching experiments. This system is integrated into a commercially available gas exchange analyzer to measure initial rates of gas exchange, precisely control ambient conditions, and monitor the conversion from 12CO2 to 13CO2.

Daylength variation affects growth, photosynthesis, leaf metabolism, partitioning, and metabolic fluxes

Daylength, a seasonal and latitudinal variable, exerts a substantial impact on plant growth. However, the relationship between daylength and growth is nonproportional, suggesting the existence of adaptive mechanisms. Thus, our study aimed to comprehensively investigate the adaptive strategies employed by plants in response to daylength variation. We grew false flax (Camelina sativa) plants, a model oilseed crop, under long-day (LD) and short-day (SD) conditions and used growth measurements, gas exchange measurements, and isotopic labeling techniques, including 13C, 14C, and 2H2O, to determine responses to different daylengths. Our findings revealed that daylength influences various growth parameters, photosynthetic physiology, carbon partitioning, metabolic fluxes, and metabolite levels. SD plants employed diverse mechanisms to compensate for reduced CO2 fixation in the shorter photoperiod. These mechanisms included enhanced photosynthetic rates and reduced respiration in the light (RL), leading to increased shoot investment. Additionally, SD plants exhibited reduced rates of the glucose 6-phosphate (G6P) shunt and greater partitioning of sugars into starch, thereby sustaining carbon availability during the longer night. Isotopic labeling results further demonstrated substantial alterations in the partitioning of amino acids and TCA cycle intermediates between rapidly and slowly turning over pools. Overall, the results point to multiple developmental, physiological, and metabolic ways in which plants adapt to different daylengths to maintain growth.

Quantitative Dynamic Analysis of de novo Sphingolipid Biosynthesis in Arabidopsis thaliana

Sphingolipids are pivotal for plant development and stress responses. Growing interest has been directed towards fully comprehending the regulatory mechanisms of the sphingolipid pathway. We explore its de novo biosynthesis and homeostasis in Arabidopsis thaliana cell cultures, shedding light on fundamental metabolic mechanisms. Employing 15N isotope labeling and quantitative dynamic modeling approach, we developed a regularized and constraint-based Dynamic Metabolic Flux Analysis (r-DMFA) framework to predict metabolic shifts due to enzymatic changes. Our analysis revealed key enzymes such as sphingoid-base hydroxylase (SBH) and long-chain-base kinase (LCBK) to be critical for maintaining sphingolipid homeostasis. Disruptions in these enzymes were found to affect cellular viability and increase the potential for programmed cell death (PCD). Thus, this work enhances our understanding of sphingolipid metabolism and demonstrates the utility of dynamic modeling in analyzing complex metabolic pathways.

Metabolic engineering of Synechococcus elongatus 7942 for enhanced sucrose biosynthesis

The capability of cyanobacteria to produce sucrose from CO2 and light has a remarkable societal and biotechnological impact since sucrose can serve as a carbon and energy source for a variety of heterotrophic organisms and can be converted into value-added products. However, most metabolic engineering efforts have focused on understanding local pathway alterations that drive sucrose biosynthesis and secretion in cyanobacteria rather than analyzing the global flux re-routing that occurs following induction of sucrose production by salt stress. Here, we investigated global metabolic flux alterations in a sucrose-secreting (cscB-overexpressing) strain relative to its wild-type Synechococcus elongatus 7942 parental strain. We used targeted metabolomics, 13C metabolic flux analysis (MFA), and genome-scale modeling (GSM) as complementary approaches to elucidate differences in cellular resource allocation by quantifying metabolic profiles of three cyanobacterial cultures - wild-type S. elongatus 7942 without salt stress (WT), wild-type with salt stress (WT/NaCl), and the cscB-overexpressing strain with salt stress (cscB/NaCl) - all under photoautotrophic conditions. We quantified the substantial rewiring of metabolic fluxes in WT/NaCl and cscB/NaCl cultures relative to WT and identified a metabolic bottleneck limiting carbon fixation and sucrose biosynthesis. This bottleneck was subsequently mitigated through heterologous overexpression of glyceraldehyde-3-phosphate dehydrogenase in an engineered sucrose-secreting strain. Our study also demonstrates that combining 13C-MFA and GSM is a useful strategy to both extend the coverage of MFA beyond central metabolism and to improve the accuracy of flux predictions provided by GSM.

Quantifying the impact of dynamic plant-environment interactions on metabolic regulation

A plant's genome encodes enzymes, transporters and many other proteins which constitute metabolism. Interactions of plants with their environment shape their growth, development and resilience towards adverse conditions. Although genome sequencing technologies and applications have experienced triumphantly rapid development during the last decades, enabling nowadays a fast and cheap sequencing of full genomes, prediction of metabolic phenotypes from genotype × environment interactions remains, at best, very incomplete. The main reasons are a lack of understanding of how different levels of molecular organisation depend on each other, and how they are constituted and expressed within a setup of growth conditions. Phenotypic plasticity, e.g., of the genetic model plant Arabidopsis thaliana, has provided important insights into plant-environment interactions and the resulting genotype x phenotype relationships. Here, we summarize previous and current findings about plant development in a changing environment and how this might be shaped and reflected in metabolism and its regulation. We identify current challenges in the study of plant development and metabolic regulation and provide an outlook of how methodological workflows might support the application of findings made in model systems to crops and their cultivation.

Mechanisms of metabolic adaptation in the duckweed Lemna gibba: an integrated metabolic, transcriptomic and flux analysis

BACKGROUND

Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions.

RESULTS

To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts.

CONCLUSIONS

Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.

Predicting plasticity of rosette growth and metabolic fluxes in Arabidopsis thaliana

Plants can rapidly mitigate the effects of suboptimal growth environments by phenotypic plasticity of fitness-traits. While genetic variation for phenotypic plasticity offers the means for breeding climate-resilient crop lines, accurate genomic prediction models for plasticity of fitness-related traits are still lacking. Here, we employed condition- and accession-specific metabolic models for 67 Arabidopsis thaliana accessions to dissect and predict plasticity of rosette growth to changes in nitrogen availability. We showed that specific reactions in photorespiration, linking carbon and nitrogen metabolism, as well as key pathways of central carbon metabolism exhibited substantial genetic variation for flux plasticity. We also demonstrated that, in comparison with a genomic prediction model for fresh weight (FW), genomic prediction of growth plasticity improves the predictability of FW under low nitrogen by 58.9% and by additional 15.4% when further integrating data on plasticity of metabolic fluxes. Therefore, the combination of metabolic and statistical modeling provides a stepping stone in understanding the molecular mechanisms and improving the predictability of plasticity for fitness-related traits.

13CO2-labelling and Sampling in Algae for Flux Analysis of Photosynthetic and Central Carbon Metabolism

The flux in photosynthesis can be studied by performing 13CO2 pulse labelling and analysing the temporal labelling kinetics of metabolic intermediates using gas or liquid chromatography linked to mass spectrometry. Metabolic flux analysis (MFA) is the primary approach for analysing metabolic network function and quantifying intracellular metabolic fluxes. Different MFA approaches differ based on the metabolic state (steady vs. non-steady state) and the use of stable isotope tracers. The main methodology used to investigate metabolic systems is metabolite steady state associated with stable isotope labelling experiments. Specifically, in biological systems like photoautotrophic organisms, isotopic non-stationary 113C metabolic flux analysis at metabolic steady state with transient isotopic labelling (13C-INST-MFA) is required. The common requirement for metabolic steady state, alongside its very short half-timed reactions, complicates robust MFA of photosynthetic metabolism. While custom gas chambers design has addressed these challenges in various model plants, no similar tools were developed for liquid photosynthetic cultures (e.g., algae, cyanobacteria), where diffusion and equilibration of inorganic carbon species in the medium entails a new dimension of complexity. Recently, a novel tailor-made microfluidics labelling system has been introduced, supplying short 13CO2 pulses at steady state, and resolving fluxes across most photosynthetic metabolic pathways in algae. The system involves injecting algal cultures and medium containing pre-equilibrated inorganic 13C into a microfluidic mixer, followed by rapid metabolic quenching, enabling precise seconds-level label pulses. This was complemented by a 13CO2-bubbling-based open labelling system (photobioreactor), allowing long pulses (minutes-hours) required for investigating fluxes into central C metabolism and major products. This combined labelling procedure provides a comprehensive fluxome cover for most algal photosynthetic and central C metabolism pathways, thus allowing comparative flux analyses across algae and plants.

Different photorespiratory mechanisms in conifer leaves, where peroxisomes have intrinsically low catalase activity

Photorespiration is an essential metabolic mechanism associated with photosynthesis; however, little is known about the photorespiratory pathway of conifer gymnosperms. Metabolite analyses of the leaves of 27 tree species showed that the mean glycerate content in conifer leaves was lower than that in angiosperm leaves. We performed experiments where [13 C]-serine was fed to detached shoots of a conifer (Cryptomeria japonica), via the transpiration stream, and compared the labeling patterns of photorespiratory metabolites with those of an angiosperm tree (Populus nigra), because glycerate is produced from serine via hydroxypyruvate in peroxisomes. In P. nigra, hydroxypyruvate, glycerate and glycine were labeled with 13 C, whereas in C. japonica, glycolate and a non-canonical photorespiratory metabolite, formate, were also labeled, suggesting that an H2 O2 -mediated non-enzymatic decarboxylation (NED) reaction occurs in C. japonica. We analyzed changes in the metabolite contents of leaves kept in the dark and leaves exposed to illuminated photorespiration-promoting conditions: a positive relationship between formate and serine levels in C. japonica implied that the active C1 -metabolism pathway synthesizes serine from formate. Leaf gas exchange analyses revealed that CO2 produced through NED was recaptured by chloroplasts. Database analysis of the peroxisomal targeting signal motifs of an H2 O2 -scavenging enzyme, catalase, derived from various species, including nine coniferous species, as well as analyses of peroxisomal fractions isolated from C. japonica and P. nigra leaves indicated that conifer peroxisomes had less catalase activity. These results suggest that NED and the subsequent C1 metabolism are involved in the photorespiratory pathway of conifer leaves, where peroxisomes have intrinsically low catalase activity.

Control of histone demethylation by nuclear-localized α-ketoglutarate dehydrogenase

Methylations on nucleosomal histones play fundamental roles in regulating eukaryotic transcription. Jumonji C domain-containing histone demethylases (JMJs) dynamically control the level of histone methylations. However, how JMJ activity is generally regulated is unknown. We found that the tricarboxylic acid cycle-associated enzyme α-ketoglutarate (α-KG) dehydrogenase (KGDH) entered the nucleus, where it interacted with various JMJs to regulate α-KG-dependent histone demethylations by JMJs, and thus controlled genome-wide gene expression in plants. We show that nuclear targeting is regulated by environmental signals and that KGDH is enriched at thousands of loci in Arabidopsis thaliana. Chromatin-bound KGDH catalyzes α-KG decarboxylation and thus may limit its local availability to KGDH-coupled JMJs, inhibiting histone demethylation. Thus, our results uncover a regulatory mechanism for histone demethylations by JMJs.

Cell-type-specific metabolism in plants

Every plant organ contains tens of different cell types, each with a specialized function. These functions are intrinsically associated with specific metabolic flux distributions that permit the synthesis of the ATP, reducing equivalents and biosynthetic precursors demanded by the cell. Investigating such cell-type-specific metabolism is complicated by the mosaic of different cells within each tissue combined with the relative scarcity of certain types. However, techniques for the isolation of specific cells, their analysis in situ by microscopy, or modeling of their function in silico have permitted insight into cell-type-specific metabolism. In this review we present some of the methods used in the analysis of cell-type-specific metabolism before describing what we know about metabolism in several cell types that have been studied in depth; (i) leaf source and sink cells; (ii) glandular trichomes that are capable of rapid synthesis of specialized metabolites; (iii) guard cells that must accumulate large quantities of the osmolytes needed for stomatal opening; (iv) cells of seeds involved in storage of reserves; and (v) the mesophyll and bundle sheath cells of C4 plants that participate in a CO₂ concentrating cycle. Metabolism is discussed in terms of its principal features, connection to cell function and what factors affect the flux distribution. Demand for precursors and energy, availability of substrates and suppression of deleterious processes are identified as key factors in shaping cell-type-specific metabolism.

Model Validation and Selection in Metabolic Flux Analysis and Flux Balance Analysis

13C-Metabolic Flux Analysis (13C-MFA) and Flux Balance Analysis (FBA) are widely used to investigate the operation of biochemical networks in both biological and biotechnological research. Both of these methods use metabolic reaction network models of metabolism operating at steady state, so that reaction rates (fluxes) and the levels of metabolic intermediates are constrained to be invariant. They provide estimated (MFA) or predicted (FBA) values of the fluxes through the network in vivo, which cannot be measured directly. A number of approaches have been taken to test the reliability of estimates and predictions from constraint-based methods and to decide on and/or discriminate between alternative model architectures. Despite advances in other areas of the statistical evaluation of metabolic models, validation and model selection methods have been underappreciated and underexplored. We review the history and state-of-the-art in constraint-based metabolic model validation and model selection. Applications and limitations of the χ2-test of goodness-of-fit, the most widely used quantitative validation and selection approach in 13C-MFA, are discussed, and complementary and alternative forms of validation and selection are proposed. A combined model validation and selection framework for 13C-MFA incorporating metabolite pool size information that leverages new developments in the field is presented and advocated for. Finally, we discuss how the adoption of robust validation and selection procedures can enhance confidence in constraint-based modeling as a whole and ultimately facilitate more widespread use of FBA in biotechnology in particular.

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.

Expression of malic enzyme reveals subcellular carbon partitioning for storage reserve production in soybeans

Central metabolism produces amino and fatty acids for protein and lipids that establish seed value. Biosynthesis of storage reserves occurs in multiple organelles that exchange central intermediates including two essential metabolites, malate, and pyruvate that are linked by malic enzyme. Malic enzyme can be active in multiple subcellular compartments, partitioning carbon and reducing equivalents for anabolic and catabolic requirements. Prior studies based on isotopic labeling and steady-state metabolic flux analyses indicated malic enzyme provides carbon for fatty acid biosynthesis in plants, though genetic evidence confirming this role is lacking. We hypothesized that increasing malic enzyme flux would alter carbon partitioning and result in increased lipid levels in soybeans. Homozygous transgenic soybean plants expressing Arabidopsis malic enzyme alleles, targeting the translational products to plastid or outside the plastid during seed development, were verified by transcript and enzyme activity analyses, organelle proteomics, and transient expression assays. Protein, oil, central metabolites, cofactors, and acyl-acyl carrier protein (ACPs) levels were quantified overdevelopment. Amino and fatty acid levels were altered resulting in an increase in lipids by 0.5-2% of seed biomass (i.e. 2-9% change in oil). Subcellular targeting of a single gene product in central metabolism impacts carbon and reducing equivalent partitioning for seed storage reserves in soybeans.

Energetic considerations for engineering novel biochemistries in photosynthetic organisms

Humans have been harnessing biology to make valuable compounds for generations. From beer and biofuels to pharmaceuticals, biology provides an efficient alternative to industrial processes. With the continuing advancement of molecular tools to genetically modify organisms, biotechnology is poised to solve urgent global problems related to environment, increasing population, and public health. However, the light dependent reactions of photosynthesis are constrained to produce a fixed stoichiometry of ATP and reducing equivalents that may not match the newly introduced synthetic metabolism, leading to inefficiency or damage. While photosynthetic organisms have evolved several ways to modify the ATP/NADPH output from their thylakoid electron transport chain, it is unknown if the native energy balancing mechanisms grant enough flexibility to match the demands of the synthetic metabolism. In this review we discuss the role of photosynthesis in the biotech industry, and the energetic considerations of using photosynthesis to power synthetic biology.

Dynamic response of photorespiration in fluctuating light environments

Photorespiration is a dynamic process that is intimately linked to photosynthetic carbon assimilation. There is a growing interest in understanding carbon assimilation during dynamic conditions, but the role of photorespiration under such conditions is unclear. In this review, we discuss recent work relevant to the function of photorespiration under dynamic conditions, with a special focus on light transients. This work reveals that photorespiration is a fundamental component of the light induction of assimilation where variable diffusive processes limit CO2 exchange with the atmosphere. Additionally, metabolic interactions between photorespiration and the C3 cycle may help balance fluxes under dynamic light conditions. We further discuss how the energy demands of photorespiration present special challenges to energy balancing during dynamic conditions. We finish the review with an overview of why regulation of photorespiration may be important under dynamic conditions to maintain appropriate fluxes through metabolic pathways related to photorespiration such as nitrogen and one-carbon metabolism.

Isotopically non-stationary metabolic flux analysis of heterotrophic Arabidopsis thaliana cell cultures

Fluxes are the ultimate phenotype of metabolism and their accurate quantification is fundamental to any understanding of metabolic networks. Steady state metabolic flux analysis has been the method of choice for quantifying fluxes in heterotrophic cells, but it is unable to measure fluxes during short-lived metabolic states, such as a transient oxidative load. Isotopically non-stationary metabolic flux analysis (INST-MFA) can be performed over shorter timescales (minutes - hours) and might overcome this limitation. INST-MFA has recently been applied to photosynthesising leaves, but agriculturally important tissues such as roots and storage organs, or plants during the night are heterotrophic. Here we outline the application of INST-MFA to heterotrophic plant cells. Using INST-MFA we were able to identify changes in the fluxes supported by phosphoenolpyruvate carboxylase and malic enzyme under oxidative load, highlighting the potential of INST-MFA to measure fluxes during short-lived metabolic states. We discuss the challenges in applying INST-MFA, and highlight further development required before it can be routinely used to quantify fluxes in heterotrophic plant cells.

Integrated flux and pool size analysis in plant central metabolism reveals unique roles of glycine and serine during photorespiration

Photorespiration is an essential process juxtaposed between plant carbon and nitrogen metabolism that responds to dynamic environments. Photorespiration recycles inhibitory intermediates arising from oxygenation reactions catalysed by Rubisco back into the C₃ cycle, but it is unclear what proportions of its nitrogen-containing intermediates (glycine and serine) are exported into other metabolisms in vivo and how these pool sizes affect net CO₂ gas exchange during photorespiratory transients. Here, to address this uncertainty, we measured rates of amino acid export from photorespiration using isotopically non-stationary metabolic flux analysis. This analysis revealed that ~23-41% of the photorespiratory carbon was exported from the pathway as serine under various photorespiratory conditions. Furthermore, we determined that the build-up and relaxation of glycine pools constrained a large portion of photosynthetic acclimation during photorespiratory transients. These results reveal the unique and important roles of glycine and serine in successfully maintaining various photorespiratory fluxes that occur under environmental fluctuations in nature and providing carbon and nitrogen for metabolism.

Millet-inspired systems metabolic engineering of NUE in crops

The use of nitrogen (N) fertilizers in agriculture has a great ability to increase crop productivity. However, their excessive use has detrimental effects on the environment. Therefore, it is necessary to develop crop varieties with improved nitrogen use efficiency (NUE) that require less N but have substantial yields. Orphan crops such as millets are cultivated in limited regions and are well adapted to lower input conditions. Therefore, they serve as a rich source of beneficial traits that can be transferred into major crops to improve their NUE. This review highlights the tremendous potential of systems biology to unravel the enzymes and pathways involved in the N metabolism of millets, which can open new possibilities to generate transgenic crops with improved NUE.

Metabolic synergy in Camelina reproductive tissues for seed development

Photosynthesis in fruits is well documented, but its contribution to seed development and yield remains largely unquantified. In oilseeds, the pods are green and elevated with direct access to sunlight. With ¹³C labeling in planta and through an intact pod labeling system, a unique multi-tissue comprehensive flux model mechanistically described how pods assimilate up to one-half (33 to 45%) of seed carbon by proximal photosynthesis in Camelina sativa. By capturing integrated tissue metabolism, the studies reveal the contribution of plant architecture beyond leaves, to enable seed filling and maximize the number of viable seeds. The latent capacity of the pod wall in the absence of leaves contributes approximately 79% of seed biomass, supporting greater seed sink capacity and higher theoretical yields that suggest an opportunity for crop productivity gains.

Anthocyanin Accumulation Provides Protection against High Light Stress While Reducing Photosynthesis in Apple Leaves

The photoprotective role of anthocyanin remains controversial. In this study, we explored the effects of anthocyanin on photosynthesis and photoprotection using transgenic ‘Galaxy Gala’ apple plants overexpressing MdMYB10 under high light stress. The overexpression of MdMYB10 dramatically enhanced leaf anthocyanin accumulation, allowing more visible light to be absorbed, particularly in the green region. However, through post-transcriptional regulation, anthocyanin accumulation lowered leaf photosynthesis in both photochemical reaction and CO₂ fixation capacities. Anthocyanin accumulation also led to a decreased de-epoxidation state of the xanthophyll cycle and antioxidant capacities, but this is most likely a response to the light-shielding effect of anthocyanin, as indicated by a higher chlorophyll concentration and lower chlorophyll a/b ratio. Under laboratory conditions when detached leaves lost carbon fixation capacity due to the limitation of CO₂ supply, the photoinhibition of detached transgenic red leaves was less severe under strong white, green, or blue light, but it became more severe in response to strong red light compared with that of the wild type. In field conditions when photosynthesis was performed normally in both green and transgenic red leaves, the degree of photoinhibition was comparable between transgenic red leaves and wild type leaves, but it was less severe in transgenic young shoot bark compared with the wild type. Taken together, these data show that anthocyanin protects plants from high light stress by absorbing excessive visible light despite reducing photosynthesis.

H2 S improves salt-stress recovery via organic acid turn-over in apple seedlings

Signalling roles of hydrogen sulphide (H₂ S) in stress biology are widely reported but not sufficiently established to urge its use in agronomic practice. Our lack of quantitative understanding of the metabolic rewiring in H₂ S signalling makes it difficult to elucidate its functions in stress tolerance on the biochemical level. Here, Malus hupehensis Rehd. var. pingyiensis seedlings were first treated with salt stress for 2 weeks and then treated with four different concentrations of NaHS. Through vigorous investigations, including phenotypic analysis, ¹³ C transient labelling and targeted metabolic and transcriptomic analysis, for the first time in the seedlings of a woody fruit crop, we found out that H₂ S recycles fixed carbons through glycolysis and tricarboxylic acid cycle to inhibit the futile accumulation of carbohydrates, to maintain an efficient CO₂ assimilation, to keep a balanced starch metabolism, to produce sufficient H₂ O₂ , to maintain malate/γ-aminobutyric acid homeostasis via an H₂ O₂ -induced anion channel (aluminium-activated malate transporter) and eventually to improve salt-stress recovery. Our results systematically demonstrate the vital roles of central carbon metabolism in H₂ S signalling and clarify the mode of action of H₂ S in apple seedlings. We conclude that H₂ S signalling interacts with central carbon metabolism in a bottom-up manner to recover plant growth after salt stress.

Combination of long-term 13CO2 labeling and isotopolog profiling allows turnover analysis of photosynthetic pigments in Arabidopsis leaves

BACKGROUND

Living cells maintain and adjust structural and functional integrity by continual synthesis and degradation of metabolites and macromolecules. The maintenance and adjustment of thylakoid membrane involve turnover of photosynthetic pigments along with subunits of protein complexes. Quantifying their turnover is essential to understand the mechanisms of homeostasis and long-term acclimation of photosynthetic apparatus. Here we report methods combining whole-plant long-term ¹³CO₂ labeling and liquid chromatography - mass spectrometry (LC-MS) analysis to determine the size of non-labeled population (NLP) of carotenoids and chlorophylls (Chl) in leaf pigment extracts of partially ¹³C-labeled plants.

RESULTS

The labeling chamber enabled parallel ¹³CO₂ labeling of up to 15 plants of Arabidopsis thaliana with real-time environmental monitoring ([CO₂], light intensity, temperature, relative air humidity and pressure) and recording. No significant difference in growth or photosynthetic pigment composition was found in leaves after 7-d exposure to normal CO₂ (~ 400 ppm) or ¹³CO₂ in the labeling chamber, or in ambient air outside the labeling chamber (control). Following chromatographic separation of the pigments and mass peak assignment by high-resolution Fourier-transform ion cyclotron resonance MS, mass spectra of photosynthetic pigments were analyzed by triple quadrupole MS to calculate NLP. The size of NLP remaining after the 7-d ¹³CO₂ labeling was ~ 10.3% and ~ 11.5% for all-trans- and 9-cis-β-carotene, ~ 21.9% for lutein, ~ 18.8% for Chl a and 33.6% for Chl b, highlighting non-uniform turnover of these pigments in thylakoids. Comparable results were obtained in all replicate plants of the ¹³CO₂ labeling experiment except for three that were showing anthocyanin accumulation and growth impairment due to insufficient water supply (leading to stomatal closure and less ¹³C incorporation).

CONCLUSIONS

Our methods allow ¹³CO₂ labeling and estimation of NLP for photosynthetic pigments with high reproducibility despite potential variations in [¹³CO₂] between the experiments. The results indicate distinct turnover rates of carotenoids and Chls in thylakoid membrane, which can be investigated in the future by time course experiments. Since ¹³C enrichment can be measured in a range of compounds, long-term ¹³CO₂ labeling chamber, in combination with appropriate MS methods, facilitates turnover analysis of various metabolites and macromolecules in plants on a time scale of hours to days.

13C metabolic flux analysis: Classification and characterization from the perspective of mathematical modeling and application in physiological research of neural cell

¹³C metabolic flux analysis (¹³C-MFA) has emerged as a forceful tool for quantifying in vivo metabolic pathway activity of different biological systems. This technology plays an important role in understanding intracellular metabolism and revealing patho-physiology mechanism. Recently, it has evolved into a method family with great diversity in experiments, analytics, and mathematics. In this review, we classify and characterize the various branch of ¹³C-MFA from a unified perspective of mathematical modeling. By linking different parts in the model to each step of its workflow, the specific technologies of ¹³C-MFA are put into discussion, including the isotope labeling model (ILM), isotope pattern measuring technique, optimization algorithm and statistical method. Its application in physiological research in neural cell has also been reviewed.

An automated workflow that generates atom mappings for large-scale metabolic models and its application to Arabidopsis thaliana

Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determines cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA), which relies on the patterns of isotope labeling of metabolites in the network. The application of MFA also requires a stoichiometric model with atom mappings that are currently not available for the majority of large-scale metabolic network models, particularly of plants. While automated approaches such as the Reaction Decoder Toolkit (RDT) can produce atom mappings for individual reactions, tracing the flow of individual atoms of the entire reactions across a metabolic model remains challenging. Here we establish an automated workflow to obtain reliable atom mappings for large-scale metabolic models by refining the outcome of RDT, and apply the workflow to metabolic models of Arabidopsis thaliana. We demonstrate the accuracy of RDT through a comparative analysis with atom mappings from a large database of biochemical reactions, MetaCyc. We further show the utility of our automated workflow by simulating ¹⁵ N isotope enrichment and identifying nitrogen (N)-containing metabolites which show enrichment patterns that are informative for flux estimation in future ¹⁵ N-MFA studies of A. thaliana. The automated workflow established in this study can be readily expanded to other species for which metabolic models have been established and the resulting atom mappings will facilitate MFA and graph-theoretic structural analyses with large-scale metabolic networks.

On the rate of phytoplankton respiration in the light

The rate of algal and cyanobacterial respiration in the light is an important ecophysiological term that remains to be completely characterized and quantified. To address this issue, we exploited process-specific decarboxylation rates from flux balance analysis and isotopically nonstationary metabolic flux analysis. Our study, based on published data, suggested that decarboxylation is about 22% of net CO2 assimilation when the tricarboxylic acid cycle is completely open (characterized by the commitment of alpha ketoglutarate to amino acid synthesis and very low rates of succinate formation). This estimate was supported by calculating the decarboxylation rates required to synthesize the major components of biomass (proteins, lipids, and carbohydrates) at their typical abundance. Of the 22 CO2 molecules produced by decarboxylation (normalized to net assimilation = 100), approximately 13 were from pyruvate and 3 were from isocitrate. The remaining six units of decarboxylation were in the amino acid synthesis pathways outside the tricarboxylic acid cycle. A small additional flux came from photorespiration, decarboxylations of six phosphogluconate in the oxidative pentose phosphate pathway, and decarboxylations in the syntheses of lower-abundance compounds, including pigments and ribonucleic acids. This general approach accounted for the high decarboxylation rates in algae and cyanobacteria compared to terrestrial plants. It prompts a simple speculation for the origin of the Kok effect and helps constrain the photoautotrophic respiration rate, in the light, in the euphotic zone of the ocean and lakes.

13CO2 labeling kinetics in maize reveal impaired efficiency of C4 photosynthesis under low irradiance

C4 photosynthesis allows faster photosynthetic rates and higher water and nitrogen use efficiency than C3 photosynthesis, but at the cost of lower quantum yield due to the energy requirement of its biochemical carbon concentration mechanism. It has also been suspected that its operation may be impaired in low irradiance. To investigate fluxes under moderate and low irradiance, maize (Zea mays) was grown at 550 µmol photons m-2 s-l and 13CO2 pulse-labeling was performed at growth irradiance or several hours after transfer to 160 µmol photons m-2 s-1. Analysis by liquid chromatography/tandem mass spectrometry or gas chromatography/mass spectrometry provided information about pool size and labeling kinetics for 32 metabolites and allowed estimation of flux at many steps in C4 photosynthesis. The results highlighted several sources of inefficiency in low light. These included excess flux at phosphoenolpyruvate carboxylase, restriction of decarboxylation by NADP-malic enzyme, and a shift to increased CO2 incorporation into aspartate, less effective use of metabolite pools to drive intercellular shuttles, and higher relative and absolute rates of photorespiration. The latter provides evidence for a lower bundle sheath CO2 concentration in low irradiance, implying that operation of the CO2 concentration mechanism is impaired in this condition. The analyses also revealed rapid exchange of carbon between the Calvin-Benson cycle and the CO2-concentration shuttle, which allows rapid adjustment of the balance between CO2 concentration and assimilation, and accumulation of large amounts of photorespiratory intermediates in low light that provides a major carbon reservoir to build up C4 metabolite pools when irradiance increases.

An Arabidopsis GCMS chemical ionization technique to quantify adaptive responses in central metabolism

We present a methodology to survey central metabolism in 13CO2-labeled Arabidopsis (Arabidopsis thaliana) rosettes by ammonia positive chemical ionization-gas chromatography-mass spectrometry. This technique preserves the molecular ion cluster of methyloxime/trimethylsilyl-derivatized analytes up to 1 kDa, providing unambiguous nominal mass assignment of >200 central metabolites and 13C incorporation rates into a subset of 111 from the tricarboxylic acid (TCA) cycle, photorespiratory pathway, amino acid metabolism, shikimate pathway, and lipid and sugar metabolism. In short-term labeling assays, we observed plateau labeling of ∼35% for intermediates of the photorespiratory cycle except for glyoxylate, which reached only ∼4% labeling and was also present at molar concentrations several fold lower than other photorespiratory intermediates. This suggests photorespiratory flux may involve alternate intermediate pools besides the generally accepted route through glyoxylate. Untargeted scans showed that in illuminated leaves, noncyclic TCA cycle flux and citrate export to the cytosol revert to a cyclic flux mode following methyl jasmonate (MJ) treatment. MJ also caused a block in the photorespiratory transamination of glyoxylate to glycine. Salicylic acid treatment induced the opposite effects in both cases, indicating the antagonistic relationship of these defense signaling hormones is preserved at the metabolome level. We provide complete chemical ionization spectra for 203 Arabidopsis metabolites from central metabolism, which uniformly feature the unfragmented pseudomolecular ion as the base peak. This unbiased, soft ionization technique is a powerful screening tool to identify adaptive metabolic trends in photosynthetic tissue and represents an important advance in methodology to measure plant metabolic flux.

An integrated isotopic labeling and freeze sampling apparatus (ILSA) to support sampling leaf metabolomics at a centi-second scale

BACKGROUND

Photosynthesis close interacts with respiration and nitrogen assimilation, which determine the photosynthetic efficiency of a leaf. Accurately quantifying the metabolic fluxes in photosynthesis, respiration and nitrogen assimilation benefit the design of photosynthetic efficiency improvement. To accurately estimate metabolic fluxes, time-series data including leaf metabolism and isotopic abundance changes should be collected under precisely controlled environments. But for isotopic labelled leaves under defined environments the, time cost of manually sampling usually longer than the turnover time of several intermediates in photosynthetic metabolism. In this case, the metabolic or physiological status of leaf sample would change during the sampling, and the accuracy of metabolomics data could be compromised.

RESULTS

Here we developed an integrated isotopic labeling and freeze sampling apparatus (ILSA), which could finish freeze sampling automatically in 0.05 s. ILSA can not only be used for sampling of photosynthetic metabolism measurement, but also suit for leaf isotopic labeling experiments under controlled environments ([CO₂] and light). Combined with HPLC-MS/MS as the metabolic measurement method, we demonstrated: (1) how pool-size of photosynthetic metabolites change in dark-accumulated rice leaf, and (2) variation in photosynthetic metabolic flux between rice and Arabidopsis thaliana.

CONCLUSIONS

The development of ILSA supports the photosynthetic research on metabolism and metabolic flux analysis and provides a new tool for the study of leaf physiology.

Recent advances in proteomics and metabolomics in plants

Over the past decade, systems biology and plant-omics have increasingly become the main stream in plant biology research. New developments in mass spectrometry and bioinformatics tools, and methodological schema to integrate multi-omics data have leveraged recent advances in proteomics and metabolomics. These progresses are driving a rapid evolution in the field of plant research, greatly facilitating our understanding of the mechanistic aspects of plant metabolisms and the interactions of plants with their external environment. Here, we review the recent progresses in MS-based proteomics and metabolomics tools and workflows with a special focus on their applications to plant biology research using several case studies related to mechanistic understanding of stress response, gene/protein function characterization, metabolic and signaling pathways exploration, and natural product discovery. We also present a projection concerning future perspectives in MS-based proteomics and metabolomics development including their applications to and challenges for system biology. This review is intended to provide readers with an overview of how advanced MS technology, and integrated application of proteomics and metabolomics can be used to advance plant system biology research.

CeCaFLUX: the first web server for standardized and visual instationary 13C metabolic flux analysis

SUMMARY

The number of instationary 13C-metabolic flux (INST-MFA) studies grows every year, making it more important than ever to ensure the clarity, standardization and reproducibility of each study. We proposed CeCaFLUX, the first user-friendly web server that derives metabolic flux distribution from instationary 13C-labeled data. Flux optimization and statistical analysis are achieved through an evolutionary optimization in a parallel manner. It can visualize the flux optimizing process in real time and the ultimate flux outcome. It will also function as a database to enhance the consistency and to facilitate sharing of flux studies.

AVAILABILITY AND IMPLEMENTATION

CeCaFLUX is freely available at https://www.cecaflux.net, the source code can be downloaded at https://github.com/zhzhd82/CeCaFLUX.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Reprogramming of Plant Central Metabolism in Response to Abiotic Stresses: A Metabolomics View

Abiotic stresses rewire plant central metabolism to maintain metabolic and energy homeostasis. Metabolites involved in the plant central metabolic network serve as a hub for regulating carbon and energy metabolism under various stress conditions. In this review, we introduce recent metabolomics techniques used to investigate the dynamics of metabolic responses to abiotic stresses and analyze the trend of publications in this field. We provide an updated overview of the changing patterns in central metabolic pathways related to the metabolic responses to common stresses, including flooding, drought, cold, heat, and salinity. We extensively review the common and unique metabolic changes in central metabolism in response to major abiotic stresses. Finally, we discuss the challenges and some emerging insights in the future application of metabolomics to study plant responses to abiotic stresses.

An efficient LC-MS method for isomer separation and detection of sugars, phosphorylated sugars, and organic acids

Assessing central carbon metabolism in plants can be challenging due to the dynamic range in pool sizes, with low levels of important phosphorylated sugars relative to more abundant sugars and organic acids. Here, we report a sensitive liquid chromatography-mass spectrometry method for analysing central metabolites on a hybrid column, where both anion-exchange and hydrophilic interaction chromatography (HILIC) ligands are embedded in the stationary phase. The liquid chromatography method was developed for enhanced selectivity of 27 central metabolites in a single run with sensitivity at femtomole levels observed for most phosphorylated sugars. The method resolved phosphorylated hexose, pentose, and triose isomers that are otherwise challenging. Compared with a standard HILIC approach, these metabolites had improved peak areas using our approach due to ion enhancement or low ion suppression in the biological sample matrix. The approach was applied to investigate metabolism in high lipid-producing tobacco leaves that exhibited increased levels of acetyl-CoA, a precursor for oil biosynthesis. The application of the method to isotopologue detection and quantification was considered through evaluating 13C-labeled seeds from Camelina sativa. The method provides a means to analyse intermediates more comprehensively in central metabolism of plant tissues.

Compartment-specific energy requirements of photosynthetic carbon metabolism in Camelina sativa leaves

The oxidative pentose phosphate pathway provides cytosolic NADPH yet reduces carbon and energy use efficiency. Repressing this pathway and introducing cytosolic NADPH-dependent malate dehydrogenase may increase crop yields by ≈5%. Detailed knowledge about plant energy metabolism may aid crop improvements. Using published estimates of flux through central carbon metabolism, we phenotype energy metabolism in illuminated Camelina sativa leaves (grown at 22 °C, 500 µmol photons m^-2 s^-1) and report several findings. First, the oxidative pentose phosphate pathway (OPPP) transfers 3.3% of the NADPH consumed in the Calvin-Benson cycle to the cytosol. NADPH supply proceeds at about 10% of the rate of net carbon assimilation. However, concomitantly respired CO₂ accounts for 4.8% of total rubisco activity. Hence, 4.8% of the flux through the Calvin-Benson cycle and photorespiration is spent on supplying cytosolic NADPH, a significant investment. Associated energy requirements exceed the energy output of the OPPP. Thus, autotrophic carbon metabolism is not simply optimised for flux into carbon sinks but sacrifices carbon and energy use efficiency to support cytosolic energy metabolism. To reduce these costs, we suggest bioengineering plants with a repressed cytosolic OPPP, and an inserted cytosolic NADPH-dependent malate dehydrogenase tuned to compensate for the loss in OPPP activity (if required). Second, sucrose cycling is a minor investment in overall leaf energy metabolism but a significant investment in cytosolic energy metabolism. Third, leaf energy balancing strictly requires oxidative phosphorylation, cofactor export from chloroplasts, and peroxisomal NADH import. Fourth, mitochondria are energetically self-sufficient. Fifth, carbon metabolism has an ATP/NADPH demand ratio of 1.52 which is met if ≤ 21.7% of whole electron flux is cyclic. Sixth, electron transport has a photon use efficiency of ≥ 62%. Last, we discuss interactions between the OPPP and the cytosolic oxidation-reduction cycle in supplying leaf cytosolic NADPH.

Design and fabrication of an improved dynamic flow cuvette for 13CO2 labeling in Arabidopsis plants

BACKGROUND

Stable isotope labeling is a non-invasive, sensitive means of monitoring metabolic flux in plants. The most physiologically meaningful information is obtained from experiments that take advantage of the natural photosynthetic carbon assimilation pathway to introduce a traceable marker with minimal effects on the physiology of the organism. The fundamental substrate in isotopic labeling experiments is ¹³CO₂, which can reveal the earliest events in carbon assimilation and realistically portray downstream metabolism when administered under conditions suitable for making kinetic inferences. Efforts to improve the accuracy and resolution of whole plant labeling techniques have focused on improvements in environmental control, air flow characteristics, and harvesting methods.

RESULTS

Here we present a dynamic flow cuvette designed for single Arabidopsis thaliana labeling experiments. We have also verified its suitability for labeling Nicotiana benthamiana and essential oils in Pelargonium graveolens. Complete plans for fabrication of this device are included. The design includes three important innovations. First, uniform, circular air flow over the rosette surface is accomplished by a fan and deflector that creates a mini-cyclone effect within the chamber interior. Second, a network of circulating canals connected to a water bath provides temperature control to within ± 0.1 ºC under variable irradiance, humidity, and air flow conditions. When photosynthetically active radiation (PAR) was varied over a range of 1000 μEinsteins m^-2 s^-1 with no adjustment to the external temperature control system, the abaxial leaf temperature changed by < 3 ºC/1000 PAR. Third, the device is fully compatible with liquid nitrogen quenching of metabolic activity without perturbation of the light environment. For short labeling experiments (< 10 s), the most critical variable is the half-life (t_1/2) of the atmosphere within the chamber, which determines the maximum resolution of the labeling system. Using an infrared gas analyzer, we monitored the atmospheric half-life during the transition from ¹²CO₂ to ¹³CO₂ air at different flow rates and determined that 3.5 L min^-1 is the optimal flow rate to initiate labeling (t_1/2 ~ 5 s). Under these conditions, we observed linear incorporation of ¹³C into triose phosphate with labeling times as short as 5 s.

CONCLUSIONS

Advances in our ability to conduct short term labeling experiments are critical to understanding of the rates and control of the earliest steps in plant metabolism. Precise kinetic measurements in whole plants using ¹³CO₂ inform metabolic models and reveal control points that can be exploited in agricultural or biotechnological contexts. The dynamic labeling cuvette presented here is suitable for studying early events in carbon assimilation and provides high resolution kinetic data for studies of metabolism in intact plants under physiologically realistic scenarios.

Reimport of carbon from cytosolic and vacuolar sugar pools into the Calvin-Benson cycle explains photosynthesis labeling anomalies

SignificancePhotosynthesis metabolites are quickly labeled when ¹³CO₂ is fed to leaves, but the time course of labeling reveals additional contributing processes involved in the metabolic dynamics of photosynthesis. The existence of three such processes is demonstrated, and a metabolic flux model is developed to explore and characterize them. The model is consistent with a slow return of carbon from cytosolic and vacuolar sugars into the Calvin-Benson cycle through the oxidative pentose phosphate pathway. Our results provide insight into how carbon assimilation is integrated into the metabolic network of photosynthetic cells with implications for global carbon fluxes.

Metabolic profiles in C3, C3-C4 intermediate, C4-like, and C4 species in the genus Flaveria

C4 photosynthesis concentrates CO2 around Rubisco in the bundle sheath, favouring carboxylation over oxygenation and decreasing photorespiration. This complex trait evolved independently in >60 angiosperm lineages. Its evolution can be investigated in genera such as Flaveria (Asteraceae) that contain species representing intermediate stages between C3 and C4 photosynthesis. Previous studies have indicated that the first major change in metabolism probably involved relocation of glycine decarboxylase and photorespiratory CO2 release to the bundle sheath and establishment of intercellular shuttles to maintain nitrogen stoichiometry. This was followed by selection for a CO2-concentrating cycle between phosphoenolpyruvate carboxylase in the mesophyll and decarboxylases in the bundle sheath, and relocation of Rubisco to the latter. We have profiled 52 metabolites in nine Flaveria species and analysed 13CO2 labelling patterns for four species. Our results point to operation of multiple shuttles, including movement of aspartate in C3-C4 intermediates and a switch towards a malate/pyruvate shuttle in C4-like species. The malate/pyruvate shuttle increases from C4-like to complete C4 species, accompanied by a rise in ancillary organic acid pools. Our findings support current models and uncover further modifications of metabolism along the evolutionary path to C4 photosynthesis in the genus Flaveria.