Metabolic network reconstruction and flux variability analysis of storage synthesis in developing oilseed rape (Brassica napus L.) embryos

New insights in metabolism modelling to decipher plant-microbe interactions

Plant disease outbreaks, exacerbated by climate change, threaten food security and environmental sustainability world-wide. Plants interact with a wide range of microorganisms. The quest for resilient agriculture requires a deep insight into the molecular and ecological interplays between plants and their associated microbial communities. Omics methods, by profiling entire molecular sets, have shed light on these complex interactions. Nonetheless, deciphering the relationships among thousands of molecular components remains a formidable challenge, and studies that integrate these components into cohesive biological networks involving plants and associated microbes are still limited. Systems biology has the potential to predict the effects of biotic and abiotic perturbations on these networks. It is therefore a promising framework for addressing the full complexity of plant-microbiome interactions.

Metabolic modeling identifies determinants of thermal growth responses in Arabidopsis thaliana

Temperature is a critical environmental factor affecting nearly all plant processes, including growth, development, and yield. Yet, despite decades of research, we lack the ability to predict plant performance at different temperatures, limiting the development of climate-resilient crops. Further, there is a pressing need to bridge the gap between the prediction of physiological and molecular traits to improve our understanding and manipulation of plant temperature responses. Here, we developed the first enzyme-constrained model of Arabidopsis thaliana's metabolism, facilitating predictions of growth-related phenotypes at different temperatures. We showed that the model can be employed for in silico identification of genes that affect plant growth at suboptimal growth temperature. Using mutant lines, we validated the genes predicted to affect plant growth, demonstrating the potential of metabolic modeling in accurately predicting plant thermal responses. The temperature-dependent enzyme-constrained metabolic model provides a template that can be used for developing sophisticated strategies to engineer climate-resilient crops.

Increasing oil content in Brassica oilseed species

Brassica oilseed species are the third most important in the world, providing approximately 15 % of the total vegetable oils. Three species (Brassica rapa, B. juncea, B. napus) dominate with B. napus being the most common in Canada, China and Europe. Originally, B. napus was a crop producing seed with high erucic acid content, which still persists today, to some extent, and is used for industrial purposes. In contrast, cultivars which produce seed used for food and feed are low erucic acid cultivars which also have reduced glucosinolate content. Because of the limit to agricultural land, recent efforts have been made to increase productivity of oil crops, including Brassica oilseed species. In this article, we have detailed research in this regard. We have covered modern genetic, genomic and metabolic control analysis approaches to identifying potential targets for the manipulation of seed oil content. Details of work on the use of quantitative trait loci, genome-wide association and comparative functional genomics to highlight factors influencing seed oil accumulation are given and functional proteins which can affect this process are discussed. In summary, a wide variety of inputs are proving useful for the improvement of Brassica oilseed species, as major sources of global vegetable oil.

An improved algorithm for flux variability analysis

Flux balance analysis (FBA) is an optimization based approach to find the optimal steady state of a metabolic network, commonly of microorganisms such as yeast strains and Escherichia coli. However, the resulting solution from an FBA is typically not unique, as the optimization problem is, more often than not, degenerate. Flux variability analysis (FVA) is a method to determine the range of possible reaction fluxes that still satisfy, within some optimality factor, the original FBA problem. The resulting range of reaction fluxes can be utilized to determine metabolic reactions of high importance, amongst other analyses. In the literature, this has been done by solving [Formula: see text] linear programs (LPs), with n being the number of reactions in the metabolic network. However, FVA can be solved with less than [Formula: see text] LPs by utilizing the basic feasible solution property of bounded LPs to reduce the number of LPs that are needed to be solved. In this work, a new algorithm is proposed to solve FVA that requires less than [Formula: see text] LPs. The proposed algorithm is benchmarked on a problem set of 112 metabolic network models ranging from single cell organisms (iMM904, ect) to a human metabolic system (Recon3D). Showing a reduction in the number of LPs required to solve the FVA problem and thus the time to solve an FVA problem.

Transcriptomic analysis of rapeseed (Brassica napus. L.) seed development in Xiangride, Qinghai Plateau, reveals how its special eco-environment results in high yield in high-altitude areas

As one of the most important oil crops, rapeseed (Brassica napus) is cultivated worldwide to produce vegetable oil, animal feed, and biodiesel. As the population grows and the need for renewable energy increases, the breeding and cultivation of high-yield rapeseed varieties have become top priorities. The formation of a high rapeseed yield is so complex because it is influenced not only by genetic mechanisms but also by many environmental conditions, such as climatic conditions and different farming practices. Interestingly, many high-yield areas are located in special eco-environments, for example, in the high-altitude Xiangride area of the Qinghai Plateau. However, the molecular mechanisms underlying the formation of high yields in such a special eco-environment area remain largely unknown. Here, we conducted field yield analysis and transcriptome analysis in the Xiangride area. Compared with the yield and environmental factors in the Xinning area (a low-yielding area), we found that the relatively longer daylight length is the key to high rapeseed yield in the Xiangride area, which leads up to a 52.1% increase in rapeseed yield, especially the increase in thousand seed weight and silique number (SN). Combined with transcriptome H-cluster analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses, we can assume that the grain development of rapeseed in the Xiangride area is ahead of schedule and lasts for a long time, leading to the high-yield results in the Xiangride area, confirmed by the expression analysis by quantitative real-time polymerase chain reaction (qRT-PCR) of yield-related genes. Our results provide valuable information for further exploring the molecular mechanism underlying high yield in special ecological environments and provide a helpful reference for studying seed development characteristics in special-producing regions for Brassica napus.

Elucidation of Triacylglycerol Overproduction in the C4 Bioenergy Crop Sorghum bicolor by Constraint-Based Analysis

Upregulation of triacylglycerols (TAGs) in vegetative plant tissues such as leaves has the potential to drastically increase the energy density and biomass yield of bioenergy crops. In this context, constraint-based analysis has the promise to improve metabolic engineering strategies. Here we present a core metabolism model for the C4 biomass crop Sorghum bicolor (iTJC1414) along with a minimal model for photosynthetic CO2 assimilation, sucrose and TAG biosynthesis in C3 plants. Extending iTJC1414 to a four-cell diel model we simulate C4 photosynthesis in mature leaves with the principal photo-assimilatory product being replaced by TAG produced at different levels. Independent of specific pathways and per unit carbon assimilated, energy content and biosynthetic demands in reducing equivalents are about 1.3 to 1.4 times higher for TAG than for sucrose. For plant generic pathways, ATP- and NADPH-demands per CO2 assimilated are higher by 1.3- and 1.5-fold, respectively. If the photosynthetic supply in ATP and NADPH in iTJC1414 is adjusted to be balanced for sucrose as the sole photo-assimilatory product, overproduction of TAG is predicted to cause a substantial surplus in photosynthetic ATP. This means that if TAG synthesis was the sole photo-assimilatory process, there could be an energy imbalance that might impede the process. Adjusting iTJC1414 to a photo-assimilatory rate that approximates field conditions, we predict possible daily rates of TAG accumulation, dependent on varying ratios of carbon partitioning between exported assimilates and accumulated oil droplets (TAG, oleosin) and in dependence of activation of futile cycles of TAG synthesis and degradation. We find that, based on the capacity of leaves for photosynthetic synthesis of exported assimilates, mature leaves should be able to reach a 20% level of TAG per dry weight within one month if only 5% of the photosynthetic net assimilation can be allocated into oil droplets. From this we conclude that high TAG levels should be achievable if TAG synthesis is induced only during a final phase of the plant life cycle.

Plant cell cultures as heterologous bio-factories for secondary metabolite production

Synthetic biology has been developing rapidly in the last decade and is attracting increasing attention from many plant biologists. The production of high-value plant-specific secondary metabolites is, however, limited mostly to microbes. This is potentially problematic because of incorrect post-translational modification of proteins and differences in protein micro-compartmentalization, substrate availability, chaperone availability, product toxicity, and cytochrome p450 reductase enzymes. Unlike other heterologous systems, plant cells may be a promising alternative for the production of high-value metabolites. Several commercial plant suspension cell cultures from different plant species have been used successfully to produce valuable metabolites in a safe, low cost, and environmentally friendly manner. However, few metabolites are currently being biosynthesized using plant platforms, with the exception of the natural pigment anthocyanin. Both Arabidopsis thaliana and Nicotiana tabacum cell cultures can be developed by multiple gene transformations and CRISPR-Cas9 genome editing. Given that the introduction of heterologous biosynthetic pathways into Arabidopsis and N. tabacum is not widely used, the biosynthesis of foreign metabolites is currently limited; however, therein lies great potential. Here, we discuss the exemplary use of plant cell cultures and prospects for using A. thaliana and N. tabacum cell cultures to produce valuable plant-specific metabolites.

Characterization of effects of genetic variants via genome-scale metabolic modelling

Genome-scale metabolic networks for model plants and crops in combination with approaches from the constraint-based modelling framework have been used to predict metabolic traits and design metabolic engineering strategies for their manipulation. With the advances in technologies to generate large-scale genotyping data from natural diversity panels and other populations, genome-wide association and genomic selection have emerged as statistical approaches to determine genetic variants associated with and predictive of traits. Here, we review recent advances in constraint-based approaches that integrate genetic variants in genome-scale metabolic models to characterize their effects on reaction fluxes. Since some of these approaches have been applied in organisms other than plants, we provide a critical assessment of their applicability particularly in crops. In addition, we further dissect the inferred effects of genetic variants with respect to reaction rate constants, abundances of enzymes, and concentrations of metabolites, as main determinants of reaction fluxes and relate them with their combined effects on complex traits, like growth. Through this systematic review, we also provide a roadmap for future research to increase the predictive power of statistical approaches by coupling them with mechanistic models of metabolism.

Transcriptome integrated metabolic modeling of carbon assimilation underlying storage root development in cassava

The existing genome-scale metabolic model of carbon metabolism in cassava storage roots, rMeCBM, has proven particularly resourceful in exploring the metabolic basis for the phenotypic differences between high and low-yield cassava cultivars. However, experimental validation of predicted metabolic fluxes by carbon labeling is quite challenging. Here, we incorporated gene expression data of developing storage roots into the basic flux-balance model to minimize infeasible metabolic fluxes, denoted as rMeCBMx, thereby improving the plausibility of the simulation and predictive power. Three different conceptual algorithms, GIMME, E-Flux, and HPCOF were evaluated. The rMeCBMx-HPCOF model outperformed others in predicting carbon fluxes in the metabolism of storage roots and, in particular, was highly consistent with transcriptome of high-yield cultivars. The flux prediction was improved through the oxidative pentose phosphate pathway in cytosol, as has been reported in various studies on root metabolism, but hardly captured by simple FBA models. Moreover, the presence of fluxes through cytosolic glycolysis and alanine biosynthesis pathways were predicted with high consistency with gene expression levels. This study sheds light on the importance of prediction power in the modeling of complex plant metabolism. Integration of multi-omics data would further help mitigate the ill-posed problem of constraint-based modeling, allowing more realistic simulation.

Model-assisted identification of metabolic engineering strategies for Jatropha curcas lipid pathways

Efficient approaches to increase plant lipid production are necessary to meet current industrial demands for this important resource. While Jatropha curcas cell culture can be used for in vitro lipid production, scaling up the system for industrial applications requires an understanding of how growth conditions affect lipid metabolism and yield. Here we present a bottom-up metabolic reconstruction of J. curcas supported with labeling experiments and biomass characterization under three growth conditions. We show that the metabolic model can accurately predict growth and distribution of fluxes in cell cultures and use these findings to pinpoint energy expenditures that affect lipid biosynthesis and metabolism. In addition, by using constraint-based modeling approaches we identify network reactions whose joint manipulation optimizes lipid production. The proposed model and computational analyses provide a stepping stone for future rational optimization of other agronomically relevant traits in J. curcas.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

Metabolic acclimation-a key to enhancing photosynthesis in changing environments?

Plants adjust their photosynthetic capacity in response to their environment in a way that optimizes their yield and fitness. There is growing evidence that this acclimation is a response to changes in the leaf metabolome, but the extent to which these are linked and how this is optimized remain poorly understood. Using as an example the metabolic perturbations occurring in response to cold, we define the different stages required for acclimation, discuss the evidence for a metabolic temperature sensor, and suggest further work towards designing climate-smart crops. In particular, we discuss how constraint-based and kinetic metabolic modelling approaches can be used to generate targeted hypotheses about relevant pathways, and argue that a stronger integration of experimental and in silico studies will help us to understand the tightly regulated interplay of carbon partitioning and resource allocation required for photosynthetic acclimation to different environmental conditions.

Extended Utilization of Constraint-Based Metabolic Model in a Long-Growing Crop

The constraint-based rMeCBM-KU50 model of cassava storage root growth was analyzed to evaluate its sensitivity, with respect to reaction flux distribution and storage root growth rate, to changes in model inputted data and constraints, including sucrose uptake rate-related data—photosynthetic rate, total leaf area, total photosynthetic rate, storage root dry weight, and biomass function-related data. These mainly varied within ±90% of the model default values, although exceptions were made for the carbohydrate (−90% to 8%) and starch (−90% to 9%) contents. The results indicated that the predicted storage root growth rate was highly affected by specific sucrose uptake rates through the total photosynthetic rate and storage root dry weight variations; whereas the carbon flux distribution, direction and partitioning inclusive, was more sensitive to the variation in biomass content, particularly the carbohydrate content. This study showed that the specific sucrose uptake rate based on the total photosynthetic rate, storage root dry weight, and carbohydrate content were critical to the constraint-based metabolic modeling and deepened our understanding of the input–output relationship—specifically regarding the rMeCBM-KU50 model—providing a valuable platform for the modeling of plant metabolic systems, especially long-growing crops.

Understanding carbon utilization routes between high and low starch-producing cultivars of cassava through Flux Balance Analysis

Analysis of metabolic flux was used for system level assessment of carbon partitioning in Kasetsart 50 (KU50) and Hanatee (HN) cassava cultivars to understand the metabolic routes for their distinct phenotypes. First, the constraint-based metabolic model of cassava storage roots, rMeCBM, was developed based on the carbon assimilation pathway of cassava. Following the subcellular compartmentalization and curation to ensure full network connectivity and reflect the complexity of eukaryotic cells, cultivar specific data on sucrose uptake and biomass synthesis were input, and rMeCBM model was used to simulate storage root growth in KU50 and HN. Results showed that rMeCBM-KU50 and rMeCBM-HN models well imitated the storage root growth. The flux-sum analysis revealed that both cultivars utilized different metabolic precursors to produce energy in plastid. More carbon flux was invested in the syntheses of carbohydrates and amino acids in KU50 than in HN. Also, KU50 utilized less flux for respiration and less energy to synthesize one gram of dry storage root. These results may disclose metabolic potential of KU50 underlying its higher storage root and starch yield over HN. Moreover, sensitivity analysis indicated the robustness of rMeCBM model. The knowledge gained might be useful for identifying engineering targets for cassava yield improvement.

A genome-scale metabolic model of potato late blight suggests a photosynthesis suppression mechanism

BACKGROUND

Phytophthora infestans is a plant pathogen that causes an important plant disease known as late blight in potato plants (Solanum tuberosum) and several other solanaceous hosts. This disease is the main factor affecting potato crop production worldwide. In spite of the importance of the disease, the molecular mechanisms underlying the compatibility between the pathogen and its hosts are still unknown.

RESULTS

To explain the metabolic response of late blight, specifically photosynthesis inhibition in infected plants, we reconstructed a genome-scale metabolic network of the S. tuberosum leaf, PstM1. This metabolic network simulates the effect of this disease in the leaf metabolism. PstM1 accounts for 2751 genes, 1113 metabolic functions, 1773 gene-protein-reaction associations and 1938 metabolites involved in 2072 reactions. The optimization of the model for biomass synthesis maximization in three infection time points suggested a suppression of the photosynthetic capacity related to the decrease of metabolic flux in light reactions and carbon fixation reactions. In addition, a variation pattern in the flux of carboxylation to oxygenation reactions catalyzed by RuBisCO was also identified, likely to be associated to a defense response in the compatible interaction between P. infestans and S. tuberosum.

CONCLUSIONS

In this work, we introduced simultaneously the first metabolic network of S. tuberosum and the first genome-scale metabolic model of the compatible interaction of a plant with P. infestans.

Metabolic engineering to increase crop yield: From concept to execution

Although the return on investment over the last 20 years for mass screening of individual plant genes to improve crop performance has been low, the investment in these activities was essential to establish the infrastructure and tools of modern plant genomics. Complex traits such as crop yield are likely multigenic, and the exhaustive screening of random gene combinations to achieve yield gains is not realistic. Clearly, smart approaches must be developed. In silico analyses of plant metabolism and gene networks can move a trait discovery program beyond trial-and-error approaches and towards rational design strategies. Metabolic models employing flux-balance analysis are useful to determine the contribution of individual genes to a trait, or to compare, optimize, or even design metabolic pathways. Regulatory association networks provide a transcriptome-based view of the plant and can lead to the identification of transcription factors that control expression of multiple genes affecting a trait. In this review, the use of these models from the perspective of an Ag innovation company's trait discovery and development program will be discussed. Important decisions that can have significant impacts on the cost and timeline to develop a commercial trait will also be presented.

Assessing Flux Distribution Associated with Metabolic Specialization of Glandular Trichomes

Many aromatic plants accumulate mixtures of secondary (or specialized) metabolites in anatomical structures called glandular trichomes (GTs). Different GT types may also synthesize different mixtures of secreted metabolites, and this contributes to the enormous chemical diversity reported to occur across species. Over the past two decades, significant progress has been made in characterizing the genes and enzymes that are responsible for the unique metabolic capabilities of GTs in different lineages of flowering plants. Less is known about the processes that regulate flux distribution through precursor pathways toward metabolic end-products. We discuss here the results from a meta-analysis of genome-scale models that were developed to capture the unique metabolic capabilities of different GT types.

Fluxomics links cellular functional analyses to whole-plant phenotyping

Fluxes through metabolic pathways reflect the integration of genetic and metabolic regulations. While it is attractive to measure all the mRNAs (transcriptome), all the proteins (proteome), and a large number of the metabolites (metabolome) in a given cellular system, linking and integrating this information remains difficult. Measurement of metabolome-wide fluxes (termed the fluxome) provides an integrated functional output of the cell machinery and a better tool to link functional analyses to plant phenotyping. This review presents and discusses sets of methodologies that have been developed to measure the fluxome. First, the principles of metabolic flux analysis (MFA), its 'short time interval' version Inst-MFA, and of constraints-based methods, such as flux balance analysis and kinetic analysis, are briefly described. The use of these powerful methods for flux characterization at the cellular scale up to the organ (fruits, seeds) and whole-plant level is illustrated. The added value given by fluxomics methods for unravelling how the abiotic environment affects flux, the process, and key metabolic steps are also described. Challenges associated with the development of fluxomics and its integration with 'omics' for thorough plant and organ functional phenotyping are discussed. Taken together, these will ultimately provide crucial clues for identifying appropriate target plant phenotypes for breeding.

Plant lipidomics at the crossroads: From technology to biology driven science

The identification and quantification of lipids from plant tissues have become commonplace and many researchers now incorporate lipidomics approaches into their experimental studies. Plant lipidomics research continues to involve technological developments such as those in mass spectrometry imaging, but in large part, lipidomics approaches have matured to the point of being accessible to the novice. Here we review some important considerations for those planning to apply plant lipidomics to their biological questions, and offer suggestions for appropriate tools and practices. This article is part of a Special Issue entitled: BBALIP_Lipidomics Opinion Articles edited by Sepp Kohlwein.

Plant Microbe Interactions in Post Genomic Era: Perspectives and Applications

Deciphering plant-microbe interactions is a promising aspect to understand the benefits and the pathogenic effect of microbes and crop improvement. The advancement in sequencing technologies and various 'omics' tool has impressively accelerated the research in biological sciences in this area. The recent and ongoing developments provide a unique approach to describing these intricate interactions and test hypotheses. In the present review, we discuss the role of plant-pathogen interaction in crop improvement. The plant innate immunity has always been an important aspect of research and leads to some interesting information like the adaptation of unique immune mechanisms of plants against pathogens. The development of new techniques in the post - genomic era has greatly enhanced our understanding of the regulation of plant defense mechanisms against pathogens. The present review also provides an overview of beneficial plant-microbe interactions with special reference to Agrobacterium tumefaciens-plant interactions where plant derived signal molecules and plant immune responses are important in pathogenicity and transformation efficiency. The construction of various Genome-scale metabolic models of microorganisms and plants presented a better understanding of all metabolic interactions activated during the interactions. This review also lists the emerging repertoire of phytopathogens and its impact on plant disease resistance. Outline of different aspects of plant-pathogen interactions is presented in this review to bridge the gap between plant microbial ecology and their immune responses.

Flux Balance Analysis of Plant Metabolism: The Effect of Biomass Composition and Model Structure on Model Predictions

The biomass composition represented in constraint-based metabolic models is a key component for predicting cellular metabolism using flux balance analysis (FBA). Despite major advances in analytical technologies, it is often challenging to obtain a detailed composition of all major biomass components experimentally. Studies examining the influence of the biomass composition on the predictions of metabolic models have so far mostly been done on models of microorganisms. Little is known about the impact of varying biomass composition on flux prediction in FBA models of plants, whose metabolism is very versatile and complex because of the presence of multiple subcellular compartments. Also, the published metabolic models of plants differ in size and complexity. In this study, we examined the sensitivity of the predicted fluxes of plant metabolic models to biomass composition and model structure. These questions were addressed by evaluating the sensitivity of predictions of growth rates and central carbon metabolic fluxes to varying biomass compositions in three different genome-/large-scale metabolic models of Arabidopsis thaliana. Our results showed that fluxes through the central carbon metabolism were robust to changes in biomass composition. Nevertheless, comparisons between the predictions from three models using identical modeling constraints and objective function showed that model predictions were sensitive to the structure of the models, highlighting large discrepancies between the published models.

Modeling Rice Metabolism: From Elucidating Environmental Effects on Cellular Phenotype to Guiding Crop Improvement

Crop productivity is severely limited by various biotic and abiotic stresses. Thus, it is highly needed to understand the underlying mechanisms of environmental stress response and tolerance in plants, which could be addressed by systems biology approach. To this end, high-throughput omics profiling and in silico modeling can be considered to explore the environmental effects on phenotypic states and metabolic behaviors of rice crops at the systems level. Especially, the advent of constraint-based metabolic reconstruction and analysis paves a way to characterize the plant cellular physiology under various stresses by combining the mathematical network models with multi-omics data. Rice metabolic networks have been reconstructed since 2013 and currently six such networks are available, where five are at genome-scale. Since their publication, these models have been utilized to systematically elucidate the rice abiotic stress responses and identify agronomic traits for crop improvement. In this review, we summarize the current status of the existing rice metabolic networks and models with their applications. Furthermore, we also highlight future directions of rice modeling studies, particularly stressing how these models can be used to contextualize the affluent multi-omics data that are readily available in the public domain. Overall, we envisage a number of studies in the future, exploiting the available metabolic models to enhance the yield and quality of rice and other food crops.

A genome-scale metabolic network reconstruction of tomato (Solanum lycopersicum L.) and its application to photorespiratory metabolism

Tomato (Solanum lycopersicum L.) has been studied extensively due to its high economic value in the market, and high content in health-promoting antioxidant compounds. Tomato is also considered as an excellent model organism for studying the development and metabolism of fleshy fruits. However, the growth, yield and fruit quality of tomatoes can be affected by drought stress, a common abiotic stress for tomato. To investigate the potential metabolic response of tomato plants to drought, we reconstructed iHY3410, a genome-scale metabolic model of tomato leaf, and used this metabolic network to simulate tomato leaf metabolism. The resulting model includes 3410 genes and 2143 biochemical and transport reactions distributed across five intracellular organelles including cytosol, plastid, mitochondrion, peroxisome and vacuole. The model successfully described the known metabolic behaviour of tomato leaf under heterotrophic and phototrophic conditions. The in silico investigation of the metabolic characteristics for photorespiration and other relevant metabolic processes under drought stress suggested that: (i) the flux distributions through the mevalonate (MVA) pathway under drought were distinct from that under normal conditions; and (ii) the changes in fluxes through core metabolic pathways with varying flux ratio of RubisCO carboxylase to oxygenase may contribute to the adaptive stress response of plants. In addition, we improved on previous studies of reaction essentiality analysis for leaf metabolism by including potential alternative routes for compensating reaction knockouts. Altogether, the genome-scale model provides a sound framework for investigating tomato metabolism and gives valuable insights into the functional consequences of abiotic stresses.

Mathematical models of plant metabolism

Among various modeling approaches in plant metabolic research, applications of Constraint-Based modeling are fast increasing in recent years, apparently driven by current advances in genomics and genome sequencing. Constraint-Based modeling, the functional analysis of metabolic networks at the whole cell or genome scale, is more difficult to apply to plants than to microbes. Here we discuss recent developments in Constraint-Based modeling in plants with focus on issues of model reconstruction and flux prediction. Another topic is the emerging application of integration of Constraint-Based modeling with omics data to increase predictive power. Furthermore, advances in experimental measurements of cellular fluxes by (13)C-Metabolic Flux Analysis are highlighted, including instationary (13)C-MFA used to probe autotrophic metabolism in photosynthetic tissue in the light.

Revisiting the chlorophyll biosynthesis pathway using genome scale metabolic model of Oryza sativa japonica

Chlorophyll is one of the most important pigments present in green plants and rice is one of the major food crops consumed worldwide. We curated the existing genome scale metabolic model (GSM) of rice leaf by incorporating new compartment, reactions and transporters. We used this modified GSM to elucidate how the chlorophyll is synthesized in a leaf through a series of bio-chemical reactions spanned over different organelles using inorganic macronutrients and light energy. We predicted the essential reactions and the associated genes of chlorophyll synthesis and validated against the existing experimental evidences. Further, ammonia is known to be the preferred source of nitrogen in rice paddy fields. The ammonia entering into the plant is assimilated in the root and leaf. The focus of the present work is centered on rice leaf metabolism. We studied the relative importance of ammonia transporters through the chloroplast and the cytosol and their interlink with other intracellular transporters. Ammonia assimilation in the leaves takes place by the enzyme glutamine synthetase (GS) which is present in the cytosol (GS1) and chloroplast (GS2). Our results provided possible explanation why GS2 mutants show normal growth under minimum photorespiration and appear chlorotic when exposed to air.

Fluxes through plant metabolic networks: measurements, predictions, insights and challenges

Although the flows of material through metabolic networks are central to cell function, they are not easy to measure other than at the level of inputs and outputs. This is particularly true in plant cells, where the network spans multiple subcellular compartments and where the network may function either heterotrophically or photoautotrophically. For many years, kinetic modelling of pathways provided the only method for describing the operation of fragments of the network. However, more recently, it has become possible to map the fluxes in central carbon metabolism using the stable isotope labelling techniques of metabolic flux analysis (MFA), and to predict intracellular fluxes using constraints-based modelling procedures such as flux balance analysis (FBA). These approaches were originally developed for the analysis of microbial metabolism, but over the last decade, they have been adapted for the more demanding analysis of plant metabolic networks. Here, the principal features of MFA and FBA as applied to plants are outlined, followed by a discussion of the insights that have been gained into plant metabolic networks through the application of these time-consuming and non-trivial methods. The discussion focuses on how a system-wide view of plant metabolism has increased our understanding of network structure, metabolic perturbations and the provision of reducing power and energy for cell function. Current methodological challenges that limit the scope of plant MFA are discussed and particular emphasis is placed on the importance of developing methods for cell-specific MFA.

Improved evidence-based genome-scale metabolic models for maize leaf, embryo, and endosperm

There is a growing demand for genome-scale metabolic reconstructions for plants, fueled by the need to understand the metabolic basis of crop yield and by progress in genome and transcriptome sequencing. Methods are also required to enable the interpretation of plant transcriptome data to study how cellular metabolic activity varies under different growth conditions or even within different organs, tissues, and developmental stages. Such methods depend extensively on the accuracy with which genes have been mapped to the biochemical reactions in the plant metabolic pathways. Errors in these mappings lead to metabolic reconstructions with an inflated number of reactions and possible generation of unreliable metabolic phenotype predictions. Here we introduce a new evidence-based genome-scale metabolic reconstruction of maize, with significant improvements in the quality of the gene-reaction associations included within our model. We also present a new approach for applying our model to predict active metabolic genes based on transcriptome data. This method includes a minimal set of reactions associated with low expression genes to enable activity of a maximum number of reactions associated with high expression genes. We apply this method to construct an organ-specific model for the maize leaf, and tissue specific models for maize embryo and endosperm cells. We validate our models using fluxomics data for the endosperm and embryo, demonstrating an improved capacity of our models to fit the available fluxomics data. All models are publicly available via the DOE Systems Biology Knowledgebase and PlantSEED, and our new method is generally applicable for analysis transcript profiles from any plant, paving the way for further in silico studies with a wide variety of plant genomes.

Metabolic engineering of higher plants and algae for isoprenoid production

Isoprenoids are a class of compounds derived from the five carbon precursors, dimethylallyl diphosphate, and isopentenyl diphosphate. These molecules present incredible natural chemical diversity, which can be valuable for humans in many aspects such as cosmetics, agriculture, and medicine. However, many terpenoids are only produced in small quantities by their natural hosts and can be difficult to generate synthetically. Therefore, much interest and effort has been directed toward capturing the genetic blueprint for their biochemistry and engineering it into alternative hosts such as plants and algae. These autotrophic organisms are attractive when compared to traditional microbial platforms because of their ability to utilize atmospheric CO2 as a carbon substrate instead of supplied carbon sources like glucose. This chapter will summarize important techniques and strategies for engineering the accumulation of isoprenoid metabolites into higher plants and algae by choosing the correct host, avoiding endogenous regulatory mechanisms, and optimizing potential flux into the target compound. Future endeavors will build on these efforts by fine-tuning product accumulation levels via the vast amount of available "-omic" data and devising metabolic engineering schemes that integrate this into a whole-organism approach. With the development of high-throughput transformation protocols and synthetic biology molecular tools, we have only begun to harness the power and utility of plant and algae metabolic engineering.

Modelling central metabolic fluxes by constraint-based optimization reveals metabolic reprogramming of developing Solanum lycopersicum (tomato) fruit

Modelling of metabolic networks is a powerful tool to analyse the behaviour of developing plant organs, including fruits. Guided by our current understanding of heterotrophic metabolism of plant cells, a medium-scale stoichiometric model, including the balance of co-factors and energy, was constructed in order to describe metabolic shifts that occur through the nine sequential stages of Solanum lycopersicum (tomato) fruit development. The measured concentrations of the main biomass components and the accumulated metabolites in the pericarp, determined at each stage, were fitted in order to calculate, by derivation, the corresponding external fluxes. They were used as constraints to solve the model by minimizing the internal fluxes. The distribution of the calculated fluxes of central metabolism were then analysed and compared with known metabolic behaviours. For instance, the partition of the main metabolic pathways (glycolysis, pentose phosphate pathway, etc.) was relevant throughout fruit development. We also predicted a valid import of carbon and nitrogen by the fruit, as well as a consistent CO2 release. Interestingly, the energetic balance indicates that excess ATP is dissipated just before the onset of ripening, supporting the concept of the climacteric crisis. Finally, the apparent contradiction between calculated fluxes with low values compared with measured enzyme capacities suggest a complex reprogramming of the metabolic machinery during fruit development. With a powerful set of experimental data and an accurate definition of the metabolic system, this work provides important insight into the metabolic and physiological requirements of the developing tomato fruits.

Flux variability analysis: application to developing oilseed rape embryos using toolboxes for constraint-based modeling

Flux variability analysis enables comprehensive exploration of alternate optimal routes in a metabolic network. This method is especially useful with models such as bna572 for the developing oilseed rape embryo which is highly compartmentalized. Here, we describe a protocol for carrying out flux variability analysis on reactions and network projections of bna572 using well-established software (CellNetAnalyzer and COBRA) for constraint-based analysis of stoichiometric network reconstructions.

Integration of a constraint-based metabolic model of Brassica napus developing seeds with (13)C-metabolic flux analysis

The use of large-scale or genome-scale metabolic reconstructions for modeling and simulation of plant metabolism and integration of those models with large-scale omics and experimental flux data is becoming increasingly important in plant metabolic research. Here we report an updated version of bna572, a bottom-up reconstruction of oilseed rape (Brassica napus L.; Brassicaceae) developing seeds with emphasis on representation of biomass-component biosynthesis. New features include additional seed-relevant pathways for isoprenoid, sterol, phenylpropanoid, flavonoid, and choline biosynthesis. Being now based on standardized data formats and procedures for model reconstruction, bna572+ is available as a COBRA-compliant Systems Biology Markup Language (SBML) model and conforms to the Minimum Information Requested in the Annotation of Biochemical Models (MIRIAM) standards for annotation of external data resources. Bna572+ contains 966 genes, 671 reactions, and 666 metabolites distributed among 11 subcellular compartments. It is referenced to the Arabidopsis thaliana genome, with gene-protein-reaction (GPR) associations resolving subcellular localization. Detailed mass and charge balancing and confidence scoring were applied to all reactions. Using B. napus seed specific transcriptome data, expression was verified for 78% of bna572+ genes and 97% of reactions. Alongside bna572+ we also present a revised carbon centric model for (13)C-Metabolic Flux Analysis ((13)C-MFA) with all its reactions being referenced to bna572+ based on linear projections. By integration of flux ratio constraints obtained from (13)C-MFA and by elimination of infinite flux bounds around thermodynamically infeasible loops based on COBRA loopless methods, we demonstrate improvements in predictive power of Flux Variability Analysis (FVA). Using this combined approach we characterize the difference in metabolic flux of developing seeds of two B. napus genotypes contrasting in starch and oil content.

Isotopomer measurement techniques in metabolic flux analysis I: nuclear magnetic resonance

Two-dimensional [(1)H, (13)C] heteronuclear single quantum correlation (HSQC) spectroscopy nuclear magnetic resonance (NMR) is a comprehensive tool in metabolic flux analysis using (13)C-labeling experiments. NMR is particularly relevant when extensive isotopomer measurements are required, such as for plant cells and tissues, which contain multiple cellular compartments. Several isotope isomers (isotopomers) can be detected and their distribution extracted quantitatively from a single 2-D HSQC NMR spectrum. For example, 2-D HSQC detects the labeling patterns of adjacent carbon atoms and provides the enrichment of individual carbon atoms of the amino acids and glucosyl and mannosyl units present in hydrolysates of glycosylated protein. The HSQC analysis can quantitatively distinguish differences between the glucosyl units in the starch hydrolysate and a protein hydrolysate of plant biomass: this specifies crucial information about compartmentalization in the plant system. The peak structures obtained from the HSQC experiment show multiplet patterns that are directly related to the isotopomer abundances. These abundances have a nonlinear relationship to the fluxes via isotopomer balancing. Fluxes are obtained from the numerical solution of these balances and a stoichiometric model that includes biomass composition data as well as consumption rates of carbohydrate and nitrogen sources. Herein, we describe the methods for the experimental measurements for flux analysis, i.e., determination of the biomass composition (lipid, protein, soluble sugar, and starch) as well as detailed procedures of acid hydrolysis of protein and starch samples and NMR sample preparation, using soybean embryo culture as the model plant system. Techniques to obtain the relative intensity of 16 amino acids and glucosyl units for protein hydrolysate and the glucosyl units of starch hydrolysate of soybean embryos in 2-D HSQC NMR spectra also are provided.

Analysis of metabolic flux using dynamic labelling and metabolic modelling

Metabolic fluxes and the capacity to modulate them are a crucial component of the ability of the plant cell to react to environmental perturbations. Our ability to quantify them and to attain information concerning the regulatory mechanisms that control them is therefore essential to understand and influence metabolic networks. For all but the simplest of flux measurements labelling methods have proven to be the most informative. Both steady-state and dynamic labelling approaches have been adopted in the study of plant metabolism. Here the conceptual basis of these complementary approaches, as well as their historical application in microbial, mammalian and plant sciences, is reviewed, and an update on technical developments in label distribution analyses is provided. This is supported by illustrative cases studies involving the kinetic modelling of secondary metabolism. One issue that is particularly complex in the analysis of plant fluxes is the extensive compartmentation of the plant cell. This problem is discussed from both theoretical and experimental perspectives, and the current approaches used to address it are assessed. Finally, current limitations and future perspectives of kinetic modelling of plant metabolism are discussed.

A method for accounting for maintenance costs in flux balance analysis improves the prediction of plant cell metabolic phenotypes under stress conditions

Flux balance models of metabolism generally utilize synthesis of biomass as the main determinant of intracellular fluxes. However, the biomass constraint alone is not sufficient to predict realistic fluxes in central heterotrophic metabolism of plant cells because of the major demand on the energy budget due to transport costs and cell maintenance. This major limitation can be addressed by incorporating transport steps into the metabolic model and by implementing a procedure that uses Pareto optimality analysis to explore the trade-off between ATP and NADPH production for maintenance. This leads to a method for predicting cell maintenance costs on the basis of the measured flux ratio between the oxidative steps of the oxidative pentose phosphate pathway and glycolysis. We show that accounting for transport and maintenance costs substantially improves the accuracy of fluxes predicted from a flux balance model of heterotrophic Arabidopsis cells in culture, irrespective of the objective function used in the analysis. Moreover, when the new method was applied to cells under control, elevated temperature and hyper-osmotic conditions, only elevated temperature led to a substantial increase in cell maintenance costs. It is concluded that the hyper-osmotic conditions tested did not impose a metabolic stress, in as much as the metabolic network is not forced to devote more resources to cell maintenance.

Systems-integration of plant metabolism: means, motive and opportunity

System integration of metabolism is considered in analogy to the investigation of corporate misdemeanour. Motive, or goal-oriented explanation, provides hypotheses that can guide the investigation of network structure. Opportunity can be established by correlative analysis using large-scale omics resources. However, correlative approaches on their own remain inconclusive and seldom identify all the links in a network. Establishment of means, or the ability to act on other network components and contribute to a phenotype, is therefore crucial. This requires functional information. Integration of quantitative data in the context of pathway models provides a powerful approach to establish 'means'. This is illustrated by discussing: first, how protein abundance is regulated by a network including transcript abundance, translation and protein degradation and second, how a combination of experimentation and modelling provides information about pathway flux, an emergent network property that integrates changes in proteins and metabolites and determines composition and biomass.

Seed architecture shapes embryo metabolism in oilseed rape

Constrained to develop within the seed, the plant embryo must adapt its shape and size to fit the space available. Here, we demonstrate how this adjustment shapes metabolism of photosynthetic embryo. Noninvasive NMR-based imaging of the developing oilseed rape (Brassica napus) seed illustrates that, following embryo bending, gradients in lipid concentration became established. These were correlated with the local photosynthetic electron transport rate and the accumulation of storage products. Experimentally induced changes in embryo morphology and/or light supply altered these gradients and were accompanied by alterations in both proteome and metabolome. Tissue-specific metabolic models predicted that the outer cotyledon and hypocotyl/radicle generate the bulk of plastidic reductant/ATP via photosynthesis, while the inner cotyledon, being enclosed by the outer cotyledon, is forced to grow essentially heterotrophically. Under field-relevant high-light conditions, major contribution of the ribulose-1,5-bisphosphate carboxylase/oxygenase-bypass to seed storage metabolism is predicted for the outer cotyledon and the hypocotyl/radicle only. Differences between in vitro- versus in planta-grown embryos suggest that metabolic heterogeneity of embryo is not observable by in vitro approaches. We conclude that in vivo metabolic fluxes are locally regulated and connected to seed architecture, driving the embryo toward an efficient use of available light and space.

Seed-development programs: a systems biology-based comparison between dicots and monocots

Seeds develop differently in dicots and monocots, especially with respect to the major storage organs. High-resolution transcriptome data have provided the first insights into the molecular networks and pathway interactions that function during the development of individual seed compartments. Here, we review mainly recent data obtained by systems biology-based approaches, which have allowed researchers to construct and model complex metabolic networks and fluxes and identify key limiting steps in seed development. Comparative coexpression network analyses define evolutionarily conservative (FUS3/ABI3/LEC1) and divergent (LEC2) networks in dicots and monocots. Finally, we discuss the determination of seed size--an important yield-related characteristic--as mediated by a number of processes (maternal and epigenetic factors, fine-tuned regulation of cell death in distinct seed compartments, and endosperm growth) and underlying genes defined through mutant analyses. Altogether, systems approaches can make important contributions toward a more complete and holistic knowledge of seed biology and thus support strategies for knowledge-based molecular breeding.

Predictive modeling of biomass component tradeoffs in Brassica napus developing oilseeds based on in silico manipulation of storage metabolism

Seed oil content is a key agronomical trait, while the control of carbon allocation into different seed storage compounds is still poorly understood and hard to manipulate. Using bna572, a large-scale model of cellular metabolism in developing embryos of rapeseed (Brassica napus) oilseeds, we present an in silico approach for the analysis of carbon allocation into seed storage products. Optimal metabolic flux states were obtained by flux variability analysis based on minimization of the uptakes of substrates in the natural environment of the embryo. For a typical embryo biomass composition, flux sensitivities to changes in different storage components were derived. Upper and lower flux bounds of each reaction were categorized as oil or protein responsive. Among the most oil-responsive reactions were glycolytic reactions, while reactions related to mitochondrial ATP production were most protein responsive. To assess different biomass compositions, a tradeoff between the fractions of oil and protein was simulated. Based on flux-bound discontinuities and shadow prices along the tradeoff, three main metabolic phases with distinct pathway usage were identified. Transitions between the phases can be related to changing modes of the tricarboxylic acid cycle, reorganizing the usage of organic carbon and nitrogen sources for protein synthesis and acetyl-coenzyme A for cytosol-localized fatty acid elongation. The phase close to equal oil and protein fractions included an unexpected pathway bypassing α-ketoglutarate-oxidizing steps in the tricarboxylic acid cycle. The in vivo relevance of the findings is discussed based on literature on seed storage metabolism.

Metabolic control analysis of developing oilseed rape (Brassica napus cv Westar) embryos shows that lipid assembly exerts significant control over oil accumulation

Metabolic control analysis allows the study of metabolic regulation. We applied both single- and double-manipulation top-down control analysis to examine the control of lipid accumulation in developing oilseed rape (Brassica napus) embryos. The biosynthetic pathway was conceptually divided into two blocks of reactions (fatty acid biosynthesis (Block A), lipid assembly (Block B)) connected by a single system intermediate, the acyl-coenzyme A (acyl-CoA) pool. Single manipulation used exogenous oleate. Triclosan was used to inhibit specifically Block A, whereas diazepam selectively manipulated flux through Block B. Exogenous oleate inhibited the radiolabelling of fatty acids from [1-(14)C]acetate, but stimulated that from [U-14C]glycerol into acyl lipids. The calculation of group flux control coefficients showed that c. 70% of the metabolic control was in the lipid assembly block of reactions. Monte Carlo simulations gave an estimation of the error of the resulting group flux control coefficients as 0.27±0.06 for Block A and 0.73±0.06 for Block B. The two methods of control analysis gave very similar results and showed that Block B reactions were more important under our conditions. This contrasts notably with data from oil palm or olive fruit cultures and is important for efforts to increase oilseed rape lipid yields.

Are we ready for genome-scale modeling in plants?

As it is becoming easier and faster to generate various types of high-throughput data, one would expect that by now we should have a comprehensive systems-level understanding of biology, biochemistry, and physiology at least in major prokaryotic and eukaryotic model systems. Despite the wealth of available data, we only get a glimpse of what is going on at the molecular level from the global perspective. The major reason is the high level of cellular complexity and our limited ability to identify all (or at least important) components and their interactions in virtually infinite number of internal and external conditions. Metabolism can be modeled mathematically by the use of genome-scale models (GEMs). GEMs are in silico metabolic flux models derived from available genome annotation. These models predict the combination of flux values of a defined metabolic network given the influence of internal and external signals. GEMs have been successfully implemented to model bacterial metabolism for over a decade. However, it was not until 2009 when the first GEM for Arabidopsis thaliana cell-suspension cultures was generated. Genome-scale modeling ("GEMing") in plants brings new challenges primarily due to the missing components and complexity of plant cells represented by the existence of: (i) photosynthesis; (ii) compartmentation; (iii) variety of cell and tissue types; and (iv) diverse metabolic responses to environmental and developmental cues as well as pathogens, insects, and competing weeds. This review presents a critical discussion of the advantages of existing plant GEMs, while identifies key targets for future improvements. Plant GEMs tend to be accurate in predicting qualitative changes in selected aspects of central carbon metabolism, while secondary metabolism is largely neglected mainly due to the missing (unknown) genes and metabolites. As such, these models are suitable for exploring metabolism in plants grown in favorable conditions, but not in field-grown plants that have to cope with environmental changes in complex ecosystems. AraGEM is the first GEM describing a photosynthetic and photorespiring plant cell (Arabidopsis thaliana). We demonstrate the use of AraGEM given the current (limited) knowledge of plant metabolism and reveal the unexpected robustness of AraGEM by a series of in silico simulations. The major focus of these simulations is on the assessment of the: (i) network connectivity; (ii) influence of CO₂ and photon uptake rates on cellular growth rates and production of individual biomass components; and (iii) stability of plant central carbon metabolism with internal pH changes.

Insights into metabolic efficiency from flux analysis

The efficiency of carbon and energy flows throughout metabolism defines the potential for growth and reproductive success of plants. Understanding the basis for metabolic efficiency requires relevant definitions of efficiency as well as measurements of biochemical functions through metabolism. Here insights into the basis of efficiency provided by (13)C-based metabolic flux analysis (MFA) as well as the uses and limitations of efficiency in predictive flux balance analysis (FBA) are highlighted. (13)C-MFA studies have revealed unusual features of central metabolism in developing green seeds for the efficient use of light to conserve carbon and identified metabolic inefficiencies in plant metabolism due to dissipation of ATP by substrate cycling. Constraints-based FBA has used efficiency to guide the prediction of the growth and actual internal flux distribution of plant systems. Comparisons in a few cases have been made between flux maps measured by (13)C-based MFA and those predicted by FBA assuming one or more maximal efficiency parameters. These studies suggest that developing plant seeds and photoautotrophic microorganisms may indeed have patterns of metabolic flux that maximize efficiency. MFA and FBA are synergistic toolsets for uncovering and explaining the metabolic basis of efficiencies and inefficiencies in plant systems.

Metabolic cartography: experimental quantification of metabolic fluxes from isotopic labelling studies

For the past decade, flux maps have provided researchers with an in-depth perspective on plant metabolism. As a rapidly developing field, significant headway has been made recently in computation, experimentation, and overall understanding of metabolic flux analysis. These advances are particularly applicable to the study of plant metabolism. New dynamic computational methods such as non-stationary metabolic flux analysis are finding their place in the toolbox of metabolic engineering, allowing more organisms to be studied and decreasing the time necessary for experimentation, thereby opening new avenues by which to explore the vast diversity of plant metabolism. Also, improved methods of metabolite detection and measurement have been developed, enabling increasingly greater resolution of flux measurements and the analysis of a greater number of the multitude of plant metabolic pathways. Methods to deconvolute organelle-specific metabolism are employed with increasing effectiveness, elucidating the compartmental specificity inherent in plant metabolism. Advances in metabolite measurements have also enabled new types of experiments, such as the calculation of metabolic fluxes based on (13)CO(2) dynamic labelling data, and will continue to direct plant metabolic engineering. Newly calculated metabolic flux maps reveal surprising and useful information about plant metabolism, guiding future genetic engineering of crops to higher yields. Due to the significant level of complexity in plants, these methods in combination with other systems biology measurements are necessary to guide plant metabolic engineering in the future.

Experimental flux measurements on a network scale

Metabolic flux is a fundamental property of living organisms. In recent years, methods for measuring metabolic flux in plants on a network scale have evolved further. One major challenge in studying flux in plants is the complexity of the plant's metabolism. In particular, in the presence of parallel pathways in multiple cellular compartments, the core of plant central metabolism constitutes a complex network. Hence, a common problem with the reliability of the contemporary results of (13)C-Metabolic Flux Analysis in plants is the substantial reduction in complexity that must be included in the simulated networks; this omission partly is due to limitations in computational simulations. Here, I discuss recent emerging strategies that will better address these shortcomings.