A genome-scale metabolic model accurately predicts fluxes in central carbon metabolism under stress conditions

Integration of metabolomics data into metabolic networks

Metabolite levels together with their corresponding metabolic fluxes are integrative outcomes of biochemical transformations and regulatory processes and they can be used to characterize the response of biological systems to genetic and/or environmental changes. However, while changes in transcript or to some extent protein levels can usually be traced back to one or several responsible genes, changes in fluxes and particularly changes in metabolite levels do not follow such rationale and are often the outcome of complex interactions of several components. The increasing quality and coverage of metabolomics technologies have fostered the development of computational approaches for integrating metabolic read-outs with large-scale models to predict the physiological state of a system. Constraint-based approaches, relying on the stoichiometry of the considered reactions, provide a modeling framework amenable to analyses of large-scale systems and to the integration of high-throughput data. Here we review the existing approaches that integrate metabolomics data in variants of constrained-based approaches to refine model reconstructions, to constrain flux predictions in metabolic models, and to relate network structural properties to metabolite levels. Finally, we discuss the challenges and perspectives in the developments of constraint-based modeling approaches driven by metabolomics data.

Pool size measurements facilitate the determination of fluxes at branching points in non-stationary metabolic flux analysis: the case of Arabidopsis thaliana

Pool size measurements are important for the estimation of absolute intracellular fluxes in particular scenarios based on data from heavy carbon isotope experiments. Recently, steady-state fluxes estimates were obtained for central carbon metabolism in an intact illuminated rosette of Arabidopsis thaliana grown photoautotrophically (Szecowka et al., 2013; Heise et al., 2014). Fluxes were estimated therein by integrating mass-spectrometric data of the dynamics of the unlabeled metabolic fraction, data on metabolic pool sizes, partitioning of metabolic pools between cellular compartments and estimates of photosynthetically inactive pools, with a simplified model of plant central carbon metabolism. However, the fluxes were determined by treating the pool sizes as fixed parameters. Here we investigated whether and, if so, to what extent the treatment of pool sizes as parameters to be optimized in three scenarios may affect the flux estimates. The results are discussed in terms of benchmark values for canonical pathways and reactions, including starch and sucrose synthesis as well as the ribulose-1,5-bisphosphate carboxylation and oxygenation reactions. In addition, we discuss pathways emerging from a divergent branch point for which pool sizes are required for flux estimation, irrespective of the computational approach used for the simulation of the observable labeling pattern. Therefore, our findings indicate the necessity for development of techniques for accurate pool size measurements to improve the quality of flux estimates from non-stationary flux estimates in intact plant cells in the absence of alternative flux measurements.

Predictive sulfur metabolism - a field in flux

The key role of sulfur metabolites in response to biotic and abiotic stress in plants, as well as their importance in diet and health has led to a significant interest and effort in trying to understand and manipulate the production of relevant compounds. Metabolic engineering utilizes a set of theoretical tools to help rationally design modifications that enhance the production of a desired metabolite. Such approaches have proven their value in bacterial systems, however, the paucity of success stories to date in plants, suggests that challenges remain. Here, we review the most commonly used methods for understanding metabolic flux, focusing on the sulfur assimilatory pathway. We highlight known issues with both experimental and theoretical approaches, as well as presenting recent methods for integrating different modeling strategies, and progress toward an understanding of flux at the whole plant level.

Plant metabolic modeling: achieving new insight into metabolism and metabolic engineering

Models are used to represent aspects of the real world for specific purposes, and mathematical models have opened up new approaches in studying the behavior and complexity of biological systems. However, modeling is often time-consuming and requires significant computational resources for data development, data analysis, and simulation. Computational modeling has been successfully applied as an aid for metabolic engineering in microorganisms. But such model-based approaches have only recently been extended to plant metabolic engineering, mainly due to greater pathway complexity in plants and their highly compartmentalized cellular structure. Recent progress in plant systems biology and bioinformatics has begun to disentangle this complexity and facilitate the creation of efficient plant metabolic models. This review highlights several aspects of plant metabolic modeling in the context of understanding, predicting and modifying complex plant metabolism. We discuss opportunities for engineering photosynthetic carbon metabolism, sucrose synthesis, and the tricarboxylic acid cycle in leaves and oil synthesis in seeds and the application of metabolic modeling to the study of plant acclimation to the environment. The aim of the review is to offer a current perspective for plant biologists without requiring specialized knowledge of bioinformatics or systems biology.

Extraction for metabolomics: access to the metabolome

INTRODUCTION

The value of information obtained from a metabolomic study depends on how much of the metabolome is present in analysed samples. Thus, only a comprehensive and reproducible extraction method will provide reliable data because the metabolites that will be measured are those that were extracted and all conclusions will be built around this information.

OBJECTIVE

To discuss the efficiency and reliability of available sample pre-treatment methods and their application in different fields of metabolomics.

METHODS

The review has three sections: the first deals with pre-extraction techniques, the second discusses the choice of extraction solvents and their main features and the third includes a brief description of the most used extraction techniques: microwave-assisted extraction, solid-phase extraction, supercritical fluid extraction, Soxhlet and a new method developed in our laboratory--the comprehensive extraction method.

RESULTS

Examination of over 200 studies showed that sample collection, homogenisation, grinding and storage could affect the yield and reproducibility of results. They also revealed that apart from the solvent used for extraction, the extraction techniques have a decisive role on the metabolites available for analysis.

CONCLUSION

It is essential to evaluate efficacy and reproducibility of sample pre-treatment as a first step to ensure the reliability of a metabolomic study. Among the reviewed methods, the comprehensive extraction method appears to provide a promising approach for extracting diverse types of metabolites.

Bottom-up Metabolic Reconstruction of Arabidopsis and Its Application to Determining the Metabolic Costs of Enzyme Production

Large-scale modeling of plant metabolism provides the possibility to compare and contrast different cellular and environmental scenarios with the ultimate aim of identifying the components underlying the respective plant behavior. The existing models of Arabidopsis (Arabidopsis thaliana) are top-down assembled, whereby the starting point is the annotated genome, in particular, the metabolic genes. Hence, dead-end metabolites and blocked reactions can arise that are subsequently addressed by using gap-filling algorithms in combination with species-unspecific genes. Here, we present a bottom-up-assembled, large-scale model that relies solely on Arabidopsis-specific annotations and results in the inclusion of only manually curated reactions. While the existing models are largely condition unspecific by employing a single biomass reaction, we provide three biomass compositions that pertain to realistic and frequently examined scenarios: carbon-limiting, nitrogen-limiting, and optimal growth conditions. The comparative analysis indicates that the proposed Arabidopsis core model exhibits comparable efficiency in carbon utilization and flexibility to the existing network alternatives. Moreover, the model is utilized to quantify the energy demand of amino acid and enzyme de novo synthesis in photoautotrophic growth conditions. Illustrated by the case of the most abundant protein in the world, Rubisco, we determine its synthesis cost in terms of ATP requirements. This, in turn, allows us to explore the tradeoff between protein synthesis and growth in Arabidopsis. Altogether, the model provides a solid basis for completely species-specific integration of high-throughput data, such as gene expression levels, and for condition-specific investigations of in silico metabolic engineering strategies.

The metabolic flux phenotype of heterotrophic Arabidopsis cells reveals a flexible balance between the cytosolic and plastidic contributions to carbohydrate oxidation in response to phosphate limitation

Understanding the mechanisms that allow plants to respond to variable and reduced availability of inorganic phosphate is of increasing agricultural importance because of the continuing depletion of the rock phosphate reserves that are used to combat inadequate phosphate levels in the soil. Changes in gene expression, protein levels, enzyme activities and metabolite levels all point to a reconfiguration of the central metabolic network in response to reduced availability of inorganic phosphate, but the metabolic significance of these changes can only be assessed in terms of the fluxes supported by the network. Steady-state metabolic flux analysis was used to define the metabolic phenotype of a heterotrophic Arabidopsis thaliana cell culture grown on a Murashige and Skoog medium containing 0, 1.25 or 5 mm inorganic phosphate. Fluxes through the central metabolic network were deduced from the redistribution of (13) C into metabolic intermediates and end products when cells were labelled with [1-(13) C], [2-(13) C], or [(13) C6 ]glucose, in combination with (14) C measurements of the rates of biomass accumulation. Analysis of the flux maps showed that reduced levels of phosphate in the growth medium stimulated flux through phosphoenolpyruvate carboxylase and malic enzyme, altered the balance between cytosolic and plastidic carbohydrate oxidation in favour of the plastid, and increased cell maintenance costs. We argue that plant cells respond to phosphate deprivation by reconfiguring the flux distribution through the pathways of carbohydrate oxidation to take advantage of better phosphate homeostasis in the plastid.

Systems analysis of metabolic phenotypes: what have we learnt?

Flux is one of the most informative measures of metabolic behavior. Its estimation requires integration of experimental and modeling approaches and, thus, is at the heart of metabolic systems biology. In this review, we argue that flux analysis and modeling of a range of plant systems points to the importance of the supply of metabolic inputs and demand for metabolic end-products as key drivers of metabolic behavior. This has implications for metabolic engineering, and the use of in silico models will be important to help design more effective engineering strategies. We also consider the importance of cell type-specific metabolism and the challenges of characterizing metabolism at this resolution. A combination of new measurement technologies and modeling approaches is bringing us closer to integrating metabolic behavior with whole-plant physiology and growth.

Metabolic flux phenotype of tobacco hairy roots engineered for increased geraniol production

The goal of this study was to characterise the metabolic flux phenotype of transgenic tobacco (Nicotiana tabacum) hairy roots engineered for increased biosynthesis of geraniol, an intermediate of the terpenoid indole alkaloid pathway. Steady state, stable isotope labelling was used to determine flux maps of central carbon metabolism for transgenic lines over-expressing (i) plastid-targeted geraniol synthase (pGES) from Valeriana officinalis, and (ii) pGES in combination with plastid-targeted geranyl pyrophosphate synthase from Arabidopsis thaliana (pGES+pGPPS), as well as for wild type and control-vector-transformed roots. Fluxes were constrained by the redistribution of label from [1-¹³C]-, [2-¹³C]- or [¹³C6]glucose into amino acids, sugars and organic acids at isotopic steady state, and by biomass output fluxes determined from the fractionation of [U-¹⁴C]glucose into insoluble polymers. No significant differences in growth and biomass composition were observed between the lines. The pGES line accumulated significant amounts of geraniol/geraniol glycosides (151±24 ng/mg dry weight) and the de novo synthesis of geraniol in pGES was confirmed by ¹³C labelling analysis. The pGES+pGPPS also accumulated geraniol and geraniol glycosides, but to lower levels than the pGES line. Although there was a distinct impact of the transgenes at the level of geraniol synthesis, other network fluxes were unaffected, reflecting the capacity of central metabolism to meet the relatively modest demand for increased precursors in the transgenic lines. It is concluded that re-engineering of the terpenoid indole alkaloid pathway will only require simultaneous manipulation of the steps producing the pathway precursors that originate in central metabolism in tissues engineered to produce at least an order of magnitude more geraniol than has been achieved so far.

Complete proteomic-based enzyme reaction and inhibition kinetics reveal how monolignol biosynthetic enzyme families affect metabolic flux and lignin in Populus trichocarpa

We established a predictive kinetic metabolic-flux model for the 21 enzymes and 24 metabolites of the monolignol biosynthetic pathway using Populus trichocarpa secondary differentiating xylem. To establish this model, a comprehensive study was performed to obtain the reaction and inhibition kinetic parameters of all 21 enzymes based on functional recombinant proteins. A total of 104 Michaelis-Menten kinetic parameters and 85 inhibition kinetic parameters were derived from these enzymes. Through mass spectrometry, we obtained the absolute quantities of all 21 pathway enzymes in the secondary differentiating xylem. This extensive experimental data set, generated from a single tissue specialized in wood formation, was used to construct the predictive kinetic metabolic-flux model to provide a comprehensive mathematical description of the monolignol biosynthetic pathway. The model was validated using experimental data from transgenic P. trichocarpa plants. The model predicts how pathway enzymes affect lignin content and composition, explains a long-standing paradox regarding the regulation of monolignol subunit ratios in lignin, and reveals novel mechanisms involved in the regulation of lignin biosynthesis. This model provides an explanation of the effects of genetic and transgenic perturbations of the monolignol biosynthetic pathway in flowering plants.

Genome-scale models of plant metabolism

A genome-scale model (GSM) is an in silico metabolic model comprising hundreds or thousands of chemical reactions that constitute the metabolic inventory of a cell, tissue, or organism. A complete, accurate GSM, in conjunction with a simulation technique such as flux balance analysis (FBA), can be used to comprehensively predict cellular metabolic flux distributions for a given genotype and given environmental conditions. Apart from enabling a user to quantitatively visualize carbon flow through metabolic pathways, these flux predictions also facilitate the hypothesis of new network properties. By simulating the impacts of environmental stresses or genetic interventions on metabolism, GSMs can aid the formulation of nontrivial metabolic engineering strategies. GSMs for plants and other eukaryotes are significantly more complicated than those for prokaryotes due to their extensive compartmentalization and size. The reconstruction of a GSM involves creating an initial model, curating the model, and then rendering the model ready for FBA. Model reconstruction involves obtaining organism-specific reactions from the annotated genome sequence or organism-specific databases. Model curation involves determining metabolite protonation status or charge, ensuring that reactions are stoichiometrically balanced, assigning reactions to appropriate subcellular compartments, deleting generic reactions or creating specific versions of them, linking dead-end metabolites, and filling of pathway gaps to complete the model. Subsequently, the model requires the addition of transport, exchange, and biomass synthesis reactions to make it FBA-ready. This cycle of editing, refining, and curation has to be performed iteratively to obtain an accurate model. This chapter outlines the reconstruction and curation of GSMs with a focus on models of plant metabolism.

Transcript abundance on its own cannot be used to infer fluxes in central metabolism

An attempt has been made to define the extent to which metabolic flux in central plant metabolism is reflected by changes in the transcriptome and metabolome, based on an analysis of in vitro cultured immature embryos of two oilseed rape (Brassica napus) accessions which contrast for seed lipid accumulation. Metabolic flux analysis (MFA) was used to constrain a flux balance metabolic model which included 671 biochemical and transport reactions within the central metabolism. This highly confident flux information was eventually used for comparative analysis of flux vs. transcript (metabolite). Metabolite profiling succeeded in identifying 79 intermediates within the central metabolism, some of which differed quantitatively between the two accessions and displayed a significant shift corresponding to flux. An RNA-Seq based transcriptome analysis revealed a large number of genes which were differentially transcribed in the two accessions, including some enzymes/proteins active in major metabolic pathways. With a few exceptions, differential activity in the major pathways (glycolysis, TCA cycle, amino acid, and fatty acid synthesis) was not reflected in contrasting abundances of the relevant transcripts. The conclusion was that transcript abundance on its own cannot be used to infer metabolic activity/fluxes in central plant metabolism. This limitation needs to be borne in mind in evaluating transcriptome data and designing metabolic engineering experiments.

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.

Optimization of steady-state ¹³C-labeling experiments for metabolic flux analysis

While steady-state (13)C metabolic flux analysis is a powerful method for deducing multiple fluxes in the central metabolic network of heterotrophic and mixotrophic plant tissues, it is also time-consuming and technically challenging. Key steps in the design and interpretation of steady-state (13)C labeling experiments are illustrated with a generic protocol based on applications to plant cell suspension cultures.

Source/sink interactions underpin crop yield: the case for trehalose 6-phosphate/SnRK1 in improvement of wheat

Considerable interest has been evoked by the analysis of the regulatory pathway in carbohydrate metabolism and cell growth involving the non-reducing disaccharide trehalose (TRE). TRE is at small concentrations in mesophytes such as Arabidopsis thaliana and Triticum aestivum, excluding a role in osmoregulation once suggested for it. Studies of TRE metabolism, and genetic modification of it, have shown a very wide and more important role of the pathway in regulation of many processes in development, growth, and photosynthesis. It has now been established that rather than TRE, it is trehalose 6-phosphate (T6P) which has such profound effects. T6P is the intermediary in TRE synthesis formed from glucose-6-phosphate and UDP-glucose, derived from sucrose, by the action of trehalose phosphate synthase. The concentration of T6P is determined both by the rate of synthesis, which depends on the sucrose concentration, and also by the rate of breakdown by trehalose-6-phosphate phosphatase which produces TRE. Changing T6P concentrations by genetically modifying the enzymes of synthesis and breakdown has altered photosynthesis, sugar metabolism, growth, and development which affect responses to, and recovery from, environmental factors. Many of the effects of T6P on metabolism and growth occur via the interaction of T6P with the SnRK1 protein kinase system. T6P inhibits the activity of SnRK1, which de-represses genes encoding proteins involved in anabolism. Consequently, a large concentration of sucrose increases T6P and thereby inhibits SnRK1, so stimulating growth of cells and their metabolic activity. The T6P/SnRK1 mechanism offers an important new view of how the distribution of assimilates to organs, such as developing grains in cereal plants, is achieved. This review briefly summarizes the factors determining, and limiting, yield of wheat (particularly mass/grain which is highly conserved) and considers how T6P/SnRK1 might function to determine grain yield and might be altered to increase them. Increasing the potential rate of filling and mass/grain are ways in which total crop yield could be increased with good husbandry which maintains crop assimilation Cereal yields globally are not increasing, despite the greater production required to meet human demand. Careful targeting of T6P is showing much promise for optimization of source/sink for yield improvement and offers yet further possibilities for increasing sink demand and grain size in wheat.

Metabolic trade-offs between biomass synthesis and photosynthate export at different light intensities in a genome-scale metabolic model of rice

Previously we have used a genome scale model of rice metabolism to describe how metabolism reconfigures at different light intensities in an expanding leaf of rice. Although this established that the metabolism of the leaf was adequately represented, in the model, the scenario was not that of the typical function of the leaf-to provide material for the rest of the plant. Here we extend our analysis to explore the transition to a source leaf as export of photosynthate increases at the expense of making leaf biomass precursors, again as a function of light intensity. In particular we investigate whether, when the leaf is making a smaller range of compounds for export to the phloem, the same changes occur in the interactions between mitochondrial and chloroplast metabolism as seen in biomass synthesis for growth when light intensity increases. Our results show that the same changes occur qualitatively, though there are slight quantitative differences reflecting differences in the energy and redox requirements for the different metabolic outputs.

Network modeling to understand plant immunity

Deciphering the networks that underpin complex biological processes using experimental data remains a significant, but promising, challenge, a task made all the harder by the added complexity of host-pathogen interactions. The aim of this article is to review the progress in understanding plant immunity made so far by applying network modeling algorithms and to show how this computational/mathematical strategy is facilitating a systems view of plant defense. We review the different types of network modeling that have been used, the data required, and the type of insight that such modeling can provide. We discuss the current challenges in modeling the regulatory networks that underlie plant defense and the future developments that may help address these challenges.

Carbon-flux distribution within Streptomyces coelicolor metabolism: a comparison between the actinorhodin-producing strain M145 and its non-producing derivative M1146

Metabolic Flux Analysis is now viewed as essential to elucidate the metabolic pattern of cells and to design appropriate genetic engineering strategies to improve strain performance and production processes. Here, we investigated carbon flux distribution in two Streptomyces coelicolor A3 (2) strains: the wild type M145 and its derivative mutant M1146, in which gene clusters encoding the four main antibiotic biosynthetic pathways were deleted. Metabolic Flux Analysis and (13)C-labeling allowed us to reconstruct a flux map under steady-state conditions for both strains. The mutant strain M1146 showed a higher growth rate, a higher flux through the pentose phosphate pathway and a higher flux through the anaplerotic phosphoenolpyruvate carboxylase. In that strain, glucose uptake and the flux through the Krebs cycle were lower than in M145. The enhanced flux through the pentose phosphate pathway in M1146 is thought to generate NADPH enough to face higher needs for biomass biosynthesis and other processes. In both strains, the production of NADPH was higher than NADPH needs, suggesting a key role for nicotinamide nucleotide transhydrogenase for redox homeostasis. ATP production is also likely to exceed metabolic ATP needs, indicating that ATP consumption for maintenance is substantial.Our results further suggest a possible competition between actinorhodin and triacylglycerol biosynthetic pathways for their common precursor, acetyl-CoA. These findings may be instrumental in developing new strategies exploiting S. coelicolor as a platform for the production of bio-based products of industrial interest.

Flux balance analysis as an alternative method to estimate fluxes without labeling

The analysis of plant metabolic networks essentially contributes to the understanding of the efficiency of plant systems in terms of their biotechnological usage. Metabolic fluxes are determined by biochemical parameters such as metabolite concentrations as well as enzyme properties and activities, which in turn are the result of various regulatory events at various levels between control of transcription and posttranslational regulation of enzyme protein activity. Thus, knowledge about metabolic fluxes on a large scale provides an integrated view on the functional state of a metabolically active cell, organ, or system. In this chapter, we introduce flux balance analysis as a constraint-based method for the prediction of optimal metabolic fluxes in a given metabolic network. Furthermore, we provide a step-by-step protocol for metabolic network reconstruction and constraint-based analysis using the COBRA Toolbox.

Multiscale metabolic modeling: dynamic flux balance analysis on a whole-plant scale

Plant metabolism is characterized by a unique complexity on the cellular, tissue, and organ levels. On a whole-plant scale, changing source and sink relations accompanying plant development add another level of complexity to metabolism. With the aim of achieving a spatiotemporal resolution of source-sink interactions in crop plant metabolism, a multiscale metabolic modeling (MMM) approach was applied that integrates static organ-specific models with a whole-plant dynamic model. Allowing for a dynamic flux balance analysis on a whole-plant scale, the MMM approach was used to decipher the metabolic behavior of source and sink organs during the generative phase of the barley (Hordeum vulgare) plant. It reveals a sink-to-source shift of the barley stem caused by the senescence-related decrease in leaf source capacity, which is not sufficient to meet the nutrient requirements of sink organs such as the growing seed. The MMM platform represents a novel approach for the in silico analysis of metabolism on a whole-plant level, allowing for a systemic, spatiotemporally resolved understanding of metabolic processes involved in carbon partitioning, thus providing a novel tool for studying yield stability and crop improvement.

Modelling metabolic CO₂ evolution--a fresh perspective on respiration

Respiration is a major contributor to net exchange of CO₂ between plants and the atmosphere and thus an important aspect of the vegetation component of global climate change models. However, a mechanistic model of respiration is lacking, and so here we explore the potential for flux balance analysis (FBA) to predict cellular CO₂ evolution rates. Metabolic flux analysis reveals that respiration is not always the dominant source of CO₂, and that metabolic processes such as the oxidative pentose phosphate pathway (OPPP) and lipid synthesis can be quantitatively important. Moreover, there is considerable variation in the metabolic origin of evolved CO₂ between tissues, species and conditions. Comparison of FBA-predicted CO₂ evolution profiles with those determined from flux measurements reveals that FBA is able to predict the metabolic origin of evolved CO₂ in different tissues/species and under different conditions. However, FBA is poor at predicting flux through certain metabolic processes such as the OPPP and we identify the way in which maintenance costs are accounted for as a major area of improvement for future FBA studies. We conclude that FBA, in its standard form, can be used to predict CO₂ evolution in a range of plant tissues and in response to environment.

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.

Diurnal changes of polysome loading track sucrose content in the rosette of wild-type arabidopsis and the starchless pgm mutant

Growth is driven by newly fixed carbon in the light, but at night it depends on reserves, like starch, that are laid down in the light. Unless plants coordinate their growth with diurnal changes in the carbon supply, they will experience acute carbon starvation during the night. Protein synthesis represents a major component of cellular growth. Polysome loading was investigated during the diurnal cycle, an extended night, and low CO2 in Arabidopsis (Arabidopsis thaliana) Columbia (Col-0) and in the starchless phosphoglucomutase (pgm) mutant. In Col-0, polysome loading was 60% to 70% in the light, 40% to 45% for much of the night, and less than 20% in an extended night, while in pgm, it fell to less than 25% early in the night. Quantification of ribosomal RNA species using quantitative reverse transcription-polymerase chain reaction revealed that polysome loading remained high for much of the night in the cytosol, was strongly light dependent in the plastid, and was always high in mitochondria. The rosette sucrose content correlated with overall and with cytosolic polysome loading. Ribosome abundance did not show significant diurnal changes. However, compared with Col-0, pgm had decreased and increased abundance of plastidic and mitochondrial ribosomes, respectively. Incorporation of label from (13)CO2 into protein confirmed that protein synthesis continues at a diminished rate in the dark. Modeling revealed that a decrease in polysome loading at night is required to balance protein synthesis with the availability of carbon from starch breakdown. Costs are also reduced by using amino acids that accumulated in the previous light period. These results uncover a tight coordination of protein synthesis with the momentary supply of carbon.

Responses to light intensity in a genome-scale model of rice metabolism

We describe the construction and analysis of a genome-scale metabolic model representing a developing leaf cell of rice (Oryza sativa) primarily derived from the annotations in the RiceCyc database. We used flux balance analysis to determine that the model represents a network capable of producing biomass precursors (amino acids, nucleotides, lipid, starch, cellulose, and lignin) in experimentally reported proportions, using carbon dioxide as the sole carbon source. We then repeated the analysis over a range of photon flux values to examine responses in the solutions. The resulting flux distributions show that (1) redox shuttles between the chloroplast, cytosol, and mitochondrion may play a significant role at low light levels, (2) photorespiration can act to dissipate excess energy at high light levels, and (3) the role of mitochondrial metabolism is likely to vary considerably according to the balance between energy demand and availability. It is notable that these organelle interactions, consistent with many experimental observations, arise solely as a result of the need for mass and energy balancing without any explicit assumptions concerning kinetic or other regulatory mechanisms.

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.

The metabolic flux phenotype of heterotrophic Arabidopsis cells reveals a complex response to changes in nitrogen supply

The extent to which individual plants utilise nitrate and ammonium, the two principal nitrogen sources in the rhizosphere, is variable and many species require a balance between the two forms for optimal growth. The effects of nitrate and ammonium on gene expression, enzyme activity and metabolite composition have been documented extensively with the aim of understanding the way in which plant cells respond to the different forms of nitrogen, but ultimately the impact of these changes on the organisation and operation of the central metabolic network can only be addressed by analysing the fluxes supported by the network. Accordingly steady-state metabolic flux analysis was used to define the metabolic phenotype of a heterotrophic Arabidopsis thaliana cell culture grown in Murashige and Skoog and ammonium-free media, treatments that influenced growth and biomass composition. Fluxes through the central metabolic network were deduced from the redistribution of label into metabolic intermediates and end products observed when cells were labelled with [1-(13) C]-, [2-(13) C]- or [(13) C6 ]glucose, in tandem with (14) C-measurements of the net accumulation of biomass. Analysis of the flux maps showed that: (i) flux through the oxidative pentose phosphate pathway varied independently of the reductant demand for biosynthesis, (ii) non-plastidic processes made a significant and variable contribution to the provision of reducing power for the plastid, and (iii) the inclusion of ammonium in the growth medium increased cell maintenance costs, in agreement with the futile cycling model of ammonium toxicity. These conclusions highlight the complexity of the metabolic response to a change in nitrogen nutrition.

A metabolite-centric view on flux distributions in genome-scale metabolic models

Background

Genome-scale metabolic models are important tools in systems biology. They permit the in-silico prediction of cellular phenotypes via mathematical optimisation procedures, most importantly flux balance analysis. Current studies on metabolic models mostly consider reaction fluxes in isolation. Based on a recently proposed metabolite-centric approach, we here describe a set of methods that enable the analysis and interpretation of flux distributions in an integrated metabolite-centric view. We demonstrate how this framework can be used for the refinement of genome-scale metabolic models.

Results

We applied the metabolite-centric view developed here to the most recent metabolic reconstruction of Escherichia coli. By compiling the balance sheets of a small number of currency metabolites, we were able to fully characterise the energy metabolism as predicted by the model and to identify a possibility for model refinement in NADPH metabolism. Selected branch points were examined in detail in order to demonstrate how a metabolite-centric view allows identifying functional roles of metabolites. Fructose 6-phosphate aldolase and the sedoheptulose bisphosphate bypass were identified as enzymatic reactions that can carry high fluxes in the model but are unlikely to exhibit significant activity in vivo. Performing a metabolite essentiality analysis, unconstrained import and export of iron ions could be identified as potentially problematic for the quality of model predictions.

Conclusions

The system-wide analysis of split ratios and branch points allows a much deeper insight into the metabolic network than reaction-centric analyses. Extending an earlier metabolite-centric approach, the methods introduced here establish an integrated metabolite-centric framework for the interpretation of flux distributions in genome-scale metabolic networks that can complement the classical reaction-centric framework. Analysing fluxes and their metabolic context simultaneously opens the door to systems biological interpretations that are not apparent from isolated reaction fluxes. Particularly powerful demonstrations of this are the analyses of the complete metabolic contexts of energy metabolism and the folate-dependent one-carbon pool presented in this work. Finally, a metabolite-centric view on flux distributions can guide the refinement of metabolic reconstructions for specific growth scenarios.

Metabolic fluxes in an illuminated Arabidopsis rosette

Photosynthesis is the basis for life, and its optimization is a key biotechnological aim given the problems of population explosion and environmental deterioration. We describe a method to resolve intracellular fluxes in intact Arabidopsis thaliana rosettes based on time-dependent labeling patterns in the metabolome. Plants photosynthesizing under limiting irradiance and ambient CO2 in a custom-built chamber were transferred into a (13)CO2-enriched environment. The isotope labeling patterns of 40 metabolites were obtained using liquid or gas chromatography coupled to mass spectrometry. Labeling kinetics revealed striking differences between metabolites. At a qualitative level, they matched expectations in terms of pathway topology and stoichiometry, but some unexpected features point to the complexity of subcellular and cellular compartmentation. To achieve quantitative insights, the data set was used for estimating fluxes in the framework of kinetic flux profiling. We benchmarked flux estimates to four classically determined flux signatures of photosynthesis and assessed the robustness of the estimates with respect to different features of the underlying metabolic model and the time-resolved data set.

Metabolic reconstruction and flux analysis of industrial Pichia yeasts

Pichia yeasts have been recognized as important microbial cell factories in the biotechnological industry. Notably, the Pichia pastoris and Pichia stipitis species have attracted much research interest due to their unique cellular physiology and metabolic capability: P. pastoris has the ability to utilize methanol for cell growth and recombinant protein production, while P. stipitis is capable of assimilating xylose to produce ethanol under oxygen-limited conditions. To harness these characteristics for biotechnological applications, it is highly required to characterize their metabolic behavior. Recently, following the genome sequencing of these two Pichia species, genome-scale metabolic networks have been reconstructed to model the yeasts' metabolism from a systems perspective. To date, there are three genome-scale models available for each of P. pastoris and P. stipitis. In this mini-review, we provide an overview of the models, discuss certain limitations of previous studies, and propose potential future works that can be conducted to better understand and engineer Pichia yeasts for industrial applications.

Systematic comparison of C3 and C4 plants based on metabolic network analysis

BACKGROUND

The C4 photosynthetic cycle supercharges photosynthesis by concentrating CO2 around ribulose-1,5-bisphosphate carboxylase and significantly reduces the oxygenation reaction. Therefore engineering C4 feature into C3 plants has been suggested as a feasible way to increase photosynthesis and yield of C3 plants, such as rice, wheat, and potato. To identify the possible transition from C3 to C4 plants, the systematic comparison of C3 and C4 metabolism is necessary.

RESULTS

We compared C3 and C4 metabolic networks using the improved constraint-based models for Arabidopsis and maize. By graph theory, we found the C3 network exhibit more dense topology structure than C4. The simulation of enzyme knockouts demonstrated that both C3 and C4 networks are very robust, especially when optimizing CO2 fixation. Moreover, C4 plant has better robustness no matter the objective function is biomass synthesis or CO2 fixation. In addition, all the essential reactions in C3 network are also essential for C4, while there are some other reactions specifically essential for C4, which validated that the basic metabolism of C4 plant is similar to C3, but C4 is more complex. We also identified more correlated reaction sets in C4, and demonstrated C4 plants have better modularity with complex mechanism coordinates the reactions and pathways than that of C3 plants. We also found the increase of both biomass production and CO2 fixation with light intensity and CO2 concentration in C4 is faster than that in C3, which reflected more efficient use of light and CO2 in C4 plant. Finally, we explored the contribution of different C4 subtypes to biomass production by setting specific constraints.

CONCLUSIONS

All results are consistent with the actual situation, which indicate that Flux Balance Analysis is a powerful method to study plant metabolism at systems level. We demonstrated that in contrast to C3, C4 plants have less dense topology, higher robustness, better modularity, and higher CO2 and radiation use efficiency. In addition, preliminary analysis indicated that the rate of CO2 fixation and biomass production in PCK subtype are superior to NADP-ME and NAD-ME subtypes under enough supply of water and nitrogen.

Network analysis: tackling complex data to study plant metabolism

Incomplete knowledge of biochemical pathways makes the holistic description of plant metabolism a non-trivial undertaking. Sensitive analytical platforms, which are capable of accurately quantifying the levels of the various molecular entities of the cell, can assist in tackling this task. However, the ever-increasing amount of high-throughput data, often from multiple technologies, requires significant computational efforts for integrative analysis. Here we introduce the application of network analysis to study plant metabolism and describe the construction and analysis of correlation-based networks from (time-resolved) metabolomics data. By investigating the interactions between metabolites, network analysis can help to interpret complex datasets through the identification of key network components. The relationship between structural and biological roles of network components can be evaluated and employed to aid metabolic engineering.

Towards multiscale plant models: integrating cellular networks

One of the ambitions of 'crop systems biology' is to combine information from molecular biology with a broader view of plant development and growth. In the context of modeling, this calls for a multiscale perspective that focuses on the interplay between cellular and macroscopic studies. With this in mind, in this review we aim to draw attention to a panel of approaches that were developed in the context of systems biology and are used for analyzing and describing the behavior of cellular networks. Ultimately, insights obtained from these methods can be exploited to refine the description of plant processes, leading to integrated plant-cellular models.

Towards multiscale plant models: integrating cellular networks

One of the ambitions of 'crop systems biology' is to combine information from molecular biology with a broader view of plant development and growth. In the context of modeling, this calls for a multiscale perspective that focuses on the interplay between cellular and macroscopic studies. With this in mind, in this review we aim to draw attention to a panel of approaches that were developed in the context of systems biology and are used for analyzing and describing the behavior of cellular networks. Ultimately, insights obtained from these methods can be exploited to refine the description of plant processes, leading to integrated plant-cellular models.

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.

Plant genome-scale metabolic reconstruction and modelling

Genome-scale metabolic reconstructions are used extensively in the study of microbial metabolism and have proven powerful tools to guide rational pathway design of industrial strains. Generation and curation of plant genome-scale metabolic models has proven far more challenging, not the least of which is our incomplete knowledge of compartmentation and organelle transporters in plants. Conversely, the potential value of modelling is far greater when exploring a complex, multi-organelle and multi-tissue metabolism. The first generation of plant genome-scale metabolic reconstructions have proven surprisingly functional and robust as well as capable of predicting many observed complex phenotypes. With further refinement, the application of these models promises to make important contributions to plant biology and metabolic engineering.

Exploring the limit of metazoan thermal tolerance via comparative proteomics: thermally induced changes in protein abundance by two hydrothermal vent polychaetes

Temperatures around hydrothermal vents are highly variable, ranging from near freezing up to 300°C. Nevertheless, animals thrive around vents, some of which live near the known limits of animal thermotolerance. Paralvinella sulfincola, an extremely thermotolerant vent polychaete, and Paralvinella palmiformis, a cooler-adapted congener, are found along the Juan de Fuca Ridge in the northwestern Pacific. We conducted shipboard high-pressure thermotolerance experiments on both species to characterize the physiological adaptations underlying P. sulfincola's pronounced thermotolerance. Quantitative proteomics, expressed sequence tag (EST) libraries and glutathione assays revealed that P. sulfincola (i) exhibited an upregulation in the synthesis and recycling of glutathione with increasing temperature, (ii) downregulated nicotinamide adenine dinucleotide (NADH) and succinate dehydrogenases (key enzymes in oxidative phosphorylation) with increasing temperature, and (iii) maintained elevated levels of heat shock proteins (HSPs) across all treatments. In contrast, P. palmiformis exhibited more typical responses to increasing temperatures (e.g. increasing HSPs at higher temperatures). These data reveal differences in how a mesotolerant and extremely thermotolerant eukaryote respond to thermal stress, and suggest that P. sulfincola's capacity to mitigate oxidative stress via increased synthesis of antioxidants and decreased flux through the mitochondrial electron transport chain enable pronounced thermotolerance. Ultimately, oxidative stress may be the key factor in limiting all metazoan thermotolerance.

Kinetic modelling of plant metabolic pathways

This paper provides a review of kinetic modelling of plant metabolic pathways as a tool for analysing their control and regulation. An overview of different modelling strategies is presented, starting with those approaches that only require a knowledge of the network stoichiometry; these are referred to as structural. Flux-balance analysis, metabolic flux analysis using isotope labelling, and elementary mode analysis are briefly mentioned as three representative examples. The main focus of this paper, however, is a discussion of kinetic modelling, which requires, in addition to the stoichiometry, a knowledge of the kinetic properties of the constituent pathway enzymes. The different types of kinetic modelling analysis, namely time-course simulation, steady-state analysis, and metabolic control analysis, are explained in some detail. An overview is presented of strategies for obtaining model parameters, as well as software tools available for simulation of such models. The kinetic modelling approach is exemplified with discussion of three models from the general plant physiology literature. With the aid of kinetic modelling it is possible to perform a control analysis of a plant metabolic system, to identify potential targets for biotechnological manipulation, as well as to ascertain the regulatory importance of different enzymes (including isoforms of the same enzyme) in a pathway. Finally, a framework is presented for extending metabolic models to the whole-plant scale by linking biochemical reactions with diffusion and advective flow through the phloem. Future challenges include explicit modelling of subcellular compartments, as well as the integration of kinetic models on the different levels of the cellular and organizational hierarchy.

Pathways and fluxes: exploring the plant metabolic network

The transition from a pathway-centred view of plant metabolism to a network-wide perspective is still incomplete. Further progress in this direction requires tools to facilitate the structural description of the network on the basis of fully annotated genomes, techniques for modelling the properties of the network, and experimental methods for constraining the models and verifying their outputs. It also requires a focus on metabolic flux as the key to understanding the regulation of metabolic activity and the relationship between the inputs and outputs of the network. Progress is being made on several fronts and this Special Issue on 'Pathways and fluxes: exploring the plant metabolic network' describes current developments in the genomic reconstruction of metabolic networks, the application of flux-balance analysis to such networks, kinetic modelling, and both steady-state-and non-steady state isotope-based measurements of multiple fluxes in the network of central carbon metabolism. The papers also highlight insights that can be obtained from pathway analysis, particularly in relation to the thermodynamic and kinetic efficiency of the predicted and observed flux distributions.

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.

Strategies for investigating the plant metabolic network with steady-state metabolic flux analysis: lessons from an Arabidopsis cell culture and other systems

Steady-state (13)C metabolic flux analysis (MFA) is currently the experimental method of choice for generating flux maps of the compartmented network of primary metabolism in heterotrophic and mixotrophic plant tissues. While statistically robust protocols for the application of steady-state MFA to plant tissues have been developed by several research groups, the implementation of the method is still far from routine. The effort required to produce a flux map is more than justified by the information that it contains about the metabolic phenotype of the system, but it remains the case that steady-state MFA is both analytically and computationally demanding. This article provides an overview of principles that underpin the implementation of steady-state MFA, focusing on the definition of the metabolic network responsible for redistribution of the label, experimental considerations relating to data collection, the modelling process that allows a set of metabolic fluxes to be deduced from the labelling data, and the interpretation of flux maps. The article draws on published studies of Arabidopsis cell cultures and other systems, including developing oilseeds, with the aim of providing practical guidance and strategies for handling the issues that arise when applying steady-state MFA to the complex metabolic networks encountered in plants.

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.

Model-based confirmation of alternative substrates of mitochondrial electron transport chain

Discrimination of metabolic models based on high throughput metabolomics data, reflecting various internal and external perturbations, is essential for identifying the components that contribute to the emerging behavior of metabolic processes. Here, we investigate 12 different models of the mitochondrial electron transport chain (ETC) in Arabidopsis thaliana during dark-induced senescence in order to elucidate the alternative substrates to this metabolic pathway. Our findings demonstrate that the coupling of the proposed computational approach, based on dynamic flux balance analysis, with time-resolved metabolomics data results in model-based confirmations of the hypotheses that, during dark-induced senescence in Arabidopsis, (i) under conditions where the main substrate for the ETC are not fully available, isovaleryl-CoA dehydrogenase and 2-hydroxyglutarate dehydrogenase are able to donate electrons to the ETC, (ii) phytanoyl-CoA does not act even as an indirect substrate of the electron transfer flavoprotein/electron-transfer flavoprotein:ubiquinone oxidoreductase complex, and (iii) the mitochondrial γ-aminobutyric acid transporter has functional significance in maintaining mitochondrial metabolism. Our study provides a basic framework for future in silico studies of alternative pathways in mitochondrial metabolism under extended darkness whereby the role of its components can be computationally discriminated based on available molecular profile data.

The upregulation of thiamine (vitamin B1) biosynthesis in Arabidopsis thaliana seedlings under salt and osmotic stress conditions is mediated by abscisic acid at the early stages of this stress response

BACKGROUND

Recent reports suggest that vitamin B1 (thiamine) participates in the processes underlying plant adaptations to certain types of abiotic and biotic stress, mainly oxidative stress. Most of the genes coding for enzymes involved in thiamine biosynthesis in Arabidopsis thaliana have been identified. In our present study, we examined the expression of thiamine biosynthetic genes, of genes encoding thiamine diphosphate-dependent enzymes and the levels of thiamine compounds during the early (sensing) and late (adaptation) responses of Arabidopsis seedlings to oxidative, salinity and osmotic stress. The possible roles of plant hormones in the regulation of the thiamine contribution to stress responses were also explored.

RESULTS

The expression of Arabidopsis genes involved in the thiamine diphosphate biosynthesis pathway, including that of THI1, THIC, TH1 and TPK, was analyzed for 48 h in seedlings subjected to NaCl or sorbitol treatment. These genes were found to be predominantly up-regulated in the early phase (2-6 h) of the stress response. The changes in these gene transcript levels were further found to correlate with increases in thiamine and its diphosphate ester content in seedlings, as well as with the enhancement of gene expression for enzymes which require thiamine diphosphate as a cofactor, mainly α-ketoglutarate dehydrogenase, pyruvate dehydrogenase and transketolase. In the case of the phytohormones including the salicylic, jasmonic and abscisic acids which are known to be involved in plant stress responses, only abscisic acid was found to significantly influence the expression of thiamine biosynthetic genes, the thiamine diphosphate levels, as well as the expression of genes coding for main thiamine diphosphate-dependent enzymes. Using Arabidopsis mutant plants defective in abscisic acid production, we demonstrate that this phytohormone is important in the regulation of THI1 and THIC gene expression during salt stress but that the regulatory mechanisms underlying the osmotic stress response are more complex.

CONCLUSIONS

On the basis of the obtained results and earlier reported data, a general model is proposed for the involvement of the biosynthesis of thiamine compounds and thiamine diphosphate-dependent enzymes in abiotic stress sensing and adaptation processes in plants. A possible regulatory role of abscisic acid in the stress sensing phase is also suggested by these data.

Functional genomics tools applied to plant metabolism: a survey on plant respiration, its connections and the annotation of complex gene functions

The application of post-genomic techniques in plant respiration studies has greatly improved our ability to assign functions to gene products. In addition it has also revealed previously unappreciated interactions between distal elements of metabolism. Such results have reinforced the need to consider plant respiratory metabolism as part of a complex network and making sense of such interactions will ultimately require the construction of predictive and mechanistic models. Transcriptomics, proteomics, metabolomics, and the quantification of metabolic flux will be of great value in creating such models both by facilitating the annotation of complex gene function, determining their structure and by furnishing the quantitative data required to test them. In this review, we highlight how these experimental approaches have contributed to our current understanding of plant respiratory metabolism and its interplay with associated process (e.g., photosynthesis, photorespiration, and nitrogen metabolism). We also discuss how data from these techniques may be integrated, with the ultimate aim of identifying mechanisms that control and regulate plant respiration and discovering novel gene functions with potential biotechnological implications.