Systems analysis of metabolic phenotypes: what have we learnt?

Comparative constraint-based modelling of fruit development across species highlights nitrogen metabolism in the growth-defence trade-off

Although primary metabolism is well conserved across species, it is useful to explore the specificity of its network to assess the extent to which some pathways may contribute to particular outcomes. Constraint-based metabolic modelling is an established framework for predicting metabolic fluxes and phenotypes and helps to explore how the plant metabolic network delivers specific outcomes from temporal series. After describing the main physiological traits during fruit development, we confirmed the correlations between fruit relative growth rate (RGR), protein content and time to maturity. Then a constraint-based method is applied to a panel of eight fruit species with a knowledge-based metabolic model of heterotrophic cells describing a generic metabolic network of primary metabolism. The metabolic fluxes are estimated by constraining the model using a large set of metabolites and compounds quantified throughout fruit development. Multivariate analyses showed a clear common pattern of flux distribution during fruit development with differences between fast- and slow-growing fruits. Only the latter fruits mobilise the tricarboxylic acid cycle in addition to glycolysis, leading to a higher rate of respiration. More surprisingly, to balance nitrogen, the model suggests, on the one hand, nitrogen uptake by nitrate reductase to support a high RGR at early stages of cucumber and, on the other hand, the accumulation of alkaloids during ripening of pepper and eggplant. Finally, building virtual fruits by combining 12 biomass compounds shows that the growth-defence trade-off is supported mainly by cell wall synthesis for fast-growing fruits and by total polyphenols accumulation for slow-growing fruits.

Metabolomics: a systems biology approach for enhancing heat stress tolerance in plants

Comprehensive metabolomic investigations provide a large set of stress-related metabolites and metabolic pathways, advancing crops under heat stress conditions. Metabolomics-assisted breeding, including mQTL and mGWAS boosted our understanding of improving numerous quantitative traits under heat stress. During the past decade, metabolomics has emerged as a fascinating scientific field that includes documentation, evaluation of metabolites, and chemical methods for cell monitoring programs in numerous plant species. A comprehensive metabolome profiling allowed the investigator to handle the comprehensive data groups of metabolites and the equivalent metabolic pathways in an extraordinary manner. Metabolomics, together with transcriptomics, plays an influential role in discovering connections between stress and genes/metabolite, phenotyping, and biomarkers documentation. Further, it helps to decode several metabolic systems connected with heat stress (HS) tolerance in plants. Heat stress is a critical environmental factor that is globally affecting the growth and productivity of plants. Thus, there is an urgent need to exploit modern breeding and biotechnological tools like metabolomics to develop cultivars with improved HS tolerance. Several studies have reported that amino acids, carbohydrates, nitrogen metabolisms, etc. and metabolites involved in the biosynthesis and catalyzing actions play a game-changing role in HS response and help plants to cope with the HS. The use of metabolomics-assisted breeding (MAB) allows a well-organized transmission of higher yield and HS tolerance at the metabolome level with specific properties. Progressive metabolomics systematic techniques have accelerated metabolic profiling. Nonetheless, continuous developments in bioinformatics, statistical tools, and databases are allowing us to produce ever-progressing, comprehensive insights into the biochemical configuration of plants and by what means this is inclined by genetic and environmental cues. Currently, assimilating metabolomics with post-genomic platforms has allowed a significant division of genetic-phenotypic connotation in several plant species. This review highlights the potential of a state-of-the-art plant metabolomics approach for the improvement of crops under HS. The development of plants with specific properties using integrated omics (metabolomics and transcriptomics) and MAB can provide new directions for future research to enhance HS tolerance in plants to achieve a goal of "zero hunger".

Influence of starch deficiency on photosynthetic and post-photosynthetic carbon isotope fractionations

Carbon isotope (13C) fractionations occurring during and after photosynthetic CO2 fixation shape the carbon isotope composition (δ13C) of plant material and respired CO2. However, responses of 13C fractionations to diel variation in starch metabolism in the leaf are not fully understood. Here we measured δ13C of organic matter (δ13COM), concentrations and δ13C of potential respiratory substrates, δ13C of dark-respired CO2 (δ13CR), and gas exchange in leaves of starch-deficient plastidial phosphoglucomutase (pgm) mutants and wild-type plants of four species (Arabidopsis thaliana, Mesembryanthemum crystallinum, Nicotiana sylvestris, and Pisum sativum). The strongest δ13C response to the pgm-induced starch deficiency was observed in N. sylvestris, with more negative δ13COM, δ13CR, and δ13C values for assimilates (i.e. sugars and starch) and organic acids (i.e. malate and citrate) in pgm mutants than in wild-type plants during a diel cycle. The genotype differences in δ13C values could be largely explained by differences in leaf gas exchange. In contrast, the PGM-knockout effect on post-photosynthetic 13C fractionations via the plastidic fructose-1,6-bisphosphate aldolase reaction or during respiration was small. Taken together, our results show that the δ13C variations in starch-deficient mutants are primarily explained by photosynthetic 13C fractionations and that the combination of knockout mutants and isotope analyses allows additional insights into plant metabolism.

Infrastructures of systems biology that facilitate functional genomic study in rice

Rice (Oryza sativa L.) is both a major staple food for the worldwide population and a model crop plant for studying the mode of action of agronomically valuable traits, providing information that can be applied to other crop plants. Due to the development of high-throughput technologies such as next generation sequencing and mass spectrometry, a huge mass of multi-omics data in rice has been accumulated. Through the integration of those data, systems biology in rice is becoming more advanced.To facilitate such systemic approaches, we have summarized current resources, such as databases and tools, for systems biology in rice. In this review, we categorize the resources using six omics levels: genomics, transcriptomics, proteomics, metabolomics, integrated omics, and functional genomics. We provide the names, websites, references, working states, and number of citations for each individual database or tool and discuss future prospects for the integrated understanding of rice gene functions.

Using a Multi-compartmental Metabolic Model to Predict Carbon Allocation in Arabidopsis thaliana

The molecular mechanism of loading/unloading of sucrose into/from the phloem plays an important role in sucrose translocation among plant tissues. Perturbation of this mechanism results in growth phenotypes of a plant. In order to better understand the coupling of sucrose translocation with metabolic processes a multi-compartmental metabolic network of Arabidopsis thaliana was reconstructed and optimized with respect to biomass growth, both in light and in dark conditions. The model can be used to perform flux balance analysis of metabolic fluxes through the central carbon metabolism and catabolic and anabolic pathways. Balances and turnover of energy (ATP/ADP) and redox metabolites (NAD(P)H/NAD(P)) as well as proton concentrations in different compartments can be estimated. Importantly, the model can be used to quantify the translocation of sucrose from source to sink tissues through phloem in association with an integral balance of protons, which in turn is defined by the operational modes of the energy metabolism (light and dark conditions). This chapter describes how a multi-compartmental model to predict carbon allocation is constructed and used.

Metabolomics in Plant Priming Research: The Way Forward?

A new era of plant biochemistry at the systems level is emerging, providing detailed descriptions of biochemical phenomena at the cellular and organismal level. This new era is marked by the advent of metabolomics—the qualitative and quantitative investigation of the entire metabolome (in a dynamic equilibrium) of a biological system. This field has developed as an indispensable methodological approach to study cellular biochemistry at a global level. For protection and survival in a constantly-changing environment, plants rely on a complex and multi-layered innate immune system. This involves surveillance of ‘self’ and ‘non-self,’ molecule-based systemic signalling and metabolic adaptations involving primary and secondary metabolites as well as epigenetic modulation mechanisms. Establishment of a pre-conditioned or primed state can sensitise or enhance aspects of innate immunity for faster and stronger responses. Comprehensive elucidation of the molecular and biochemical processes associated with the phenotypic defence state is vital for a better understanding of the molecular mechanisms that define the metabolism of plant–pathogen interactions. Such insights are essential for translational research and applications. Thus, this review highlights the prospects of metabolomics and addresses current challenges that hinder the realisation of the full potential of the field. Such limitations include partial coverage of the metabolome and maximising the value of metabolomics data (extraction of information and interpretation). Furthermore, the review points out key features that characterise both the plant innate immune system and enhancement of the latter, thus underlining insights from metabolomic studies in plant priming. Future perspectives in this inspiring area are included, with the aim of stimulating further studies leading to a better understanding of plant immunity at the metabolome level.

Metabolic model of central carbon and energy metabolisms of growing Arabidopsis thaliana in relation to sucrose translocation

BACKGROUND

Sucrose translocation between plant tissues is crucial for growth, development and reproduction of plants. Systemic analysis of these metabolic and underlying regulatory processes allow a detailed understanding of carbon distribution within the plant and the formation of associated phenotypic traits. Sucrose translocation from 'source' tissues (e.g. mesophyll) to 'sink' tissues (e.g. root) is tightly bound to the proton gradient across the membranes. The plant sucrose transporters are grouped into efflux exporters (SWEET family) and proton-symport importers (SUC, STP families). To better understand regulation of sucrose export from source tissues and sucrose import into sink tissues, there is a need for a metabolic model that takes in account the tissue organisation of Arabidopsis thaliana with corresponding metabolic specificities of respective tissues in terms of sucrose and proton production/utilization. An ability of the model to operate under different light modes ('light' and 'dark') and correspondingly in different energy producing modes is particularly important in understanding regulatory modules.

RESULTS

Here, we describe a multi-compartmental model consisting of a mesophyll cell with plastid and mitochondrion, a phloem cell, as well as a root cell with mitochondrion. In this model, the phloem was considered as a non-growing transport compartment, the mesophyll compartment was considered as both autotrophic (growing on CO₂ under light) and heterotrophic (growing on starch in darkness), and the root was always considered as heterotrophic tissue dependent on sucrose supply from the mesophyll compartment. In total, the model includes 413 balanced compounds interconnected by 400 transformers. The structured metabolic model accounts for central carbon metabolism, photosynthesis, photorespiration, carbohydrate metabolism, energy and redox metabolisms, proton metabolism, biomass growth, nutrients uptake, proton gradient generation and sucrose translocation between tissues. Biochemical processes in the model were associated with gene-products (742 ORFs). Flux Balance Analysis (FBA) of the model resulted in balanced carbon, nitrogen, proton, energy and redox states under both light and dark conditions. The main H⁺-fluxes were reconstructed and their directions matched with proton-dependent sucrose translocation from 'source' to 'sink' under any light condition.

CONCLUSIONS

The model quantified the translocation of sucrose between plant tissues in association with an integral balance of protons, which in turn is defined by operational modes of the energy metabolism.

Mass spectrometry-based plant metabolomics: Metabolite responses to abiotic stress

Metabolomics is one omics approach that can be used to acquire comprehensive information on the composition of a metabolite pool to provide a functional screen of the cellular state. Studies of the plant metabolome include analysis of a wide range of chemical species with diverse physical properties, from ionic inorganic compounds to biochemically derived hydrophilic carbohydrates, organic and amino acids, and a range of hydrophobic lipid-related compounds. This complexitiy brings huge challenges to the analytical technologies employed in current plant metabolomics programs, and powerful analytical tools are required for the separation and characterization of this extremely high compound diversity present in biological sample matrices. The use of mass spectrometry (MS)-based analytical platforms to profile stress-responsive metabolites that allow some plants to adapt to adverse environmental conditions is fundamental in current plant biotechnology research programs for the understanding and development of stress-tolerant plants. In this review, we describe recent applications of metabolomics and emphasize its increasing application to study plant responses to environmental (stress-) factors, including drought, salt, low oxygen caused by waterlogging or flooding of the soil, temperature, light and oxidative stress (or a combination of them). Advances in understanding the global changes occurring in plant metabolism under specific abiotic stress conditions are fundamental to enhance plant fitness and increase stress tolerance. © 2015 Wiley Periodicals, Inc. Mass Spec Rev 35:620-649, 2016.

Plant Metabolomics: An Indispensable System Biology Tool for Plant Science

As genomes of many plant species have been sequenced, demand for functional genomics has dramatically accelerated the improvement of other omics including metabolomics. Despite a large amount of metabolites still remaining to be identified, metabolomics has contributed significantly not only to the understanding of plant physiology and biology from the view of small chemical molecules that reflect the end point of biological activities, but also in past decades to the attempts to improve plant behavior under both normal and stressed conditions. Hereby, we summarize the current knowledge on the genetic and biochemical mechanisms underlying plant growth, development, and stress responses, focusing further on the contributions of metabolomics to practical applications in crop quality improvement and food safety assessment, as well as plant metabolic engineering. We also highlight the current challenges and future perspectives in this inspiring area, with the aim to stimulate further studies leading to better crop improvement of yield and quality.

Reduced mitochondrial malate dehydrogenase activity has a strong effect on photorespiratory metabolism as revealed by 13C labelling

Mitochondrial malate dehydrogenase (mMDH) catalyses the interconversion of malate and oxaloacetate (OAA) in the tricarboxylic acid (TCA) cycle. Its activity is important for redox control of the mitochondrial matrix, through which it may participate in regulation of TCA cycle turnover. In Arabidopsis, there are two isoforms of mMDH. Here, we investigated to which extent the lack of the major isoform, mMDH1 accounting for about 60% of the activity, affected leaf metabolism. In air, rosettes of mmdh1 plants were only slightly smaller than wild type plants although the fresh weight was decreased by about 50%. In low CO2 the difference was much bigger, with mutant plants accumulating only 14% of fresh weight as compared to wild type. To investigate the metabolic background to the differences in growth, we developed a (13)CO2 labelling method, using a custom-built chamber that enabled simultaneous treatment of sets of plants under controlled conditions. The metabolic profiles were analysed by gas- and liquid- chromatography coupled to mass spectrometry to investigate the metabolic adjustments between wild type and mmdh1 The genotypes responded similarly to high CO2 treatment both with respect to metabolite pools and (13)C incorporation during a 2-h treatment. However, under low CO2 several metabolites differed between the two genotypes and, interestingly most of these were closely associated with photorespiration. We found that while the glycine/serine ratio increased, a concomitant altered glutamine/glutamate/α-ketoglutarate relation occurred. Taken together, our results indicate that adequate mMDH activity is essential to shuttle reductants out from the mitochondria to support the photorespiratory flux, and strengthen the idea that photorespiration is tightly intertwined with peripheral metabolic reactions.

Perspectives for a better understanding of the metabolic integration of photorespiration within a complex plant primary metabolism network

Photorespiration is an essential high flux metabolic pathway that is found in all oxygen-producing photosynthetic organisms. It is often viewed as a closed metabolic repair pathway that serves to detoxify 2-phosphoglycolic acid and to recycle carbon to fuel the Calvin-Benson cycle. However, this view is too simplistic since the photorespiratory cycle is known to interact with several primary metabolic pathways, including photosynthesis, nitrate assimilation, amino acid metabolism, C1 metabolism and the Krebs (TCA) cycle. Here we will review recent advances in photorespiration research and discuss future priorities to better understand (i) the metabolic integration of the photorespiratory cycle within the complex network of plant primary metabolism and (ii) the importance of photorespiration in response to abiotic and biotic stresses.

Quantifying plant phenotypes with isotopic labeling & metabolic flux analysis

Analyses of metabolic flux using stable isotopes in plants have traditionally been restricted to tissues with presumed homogeneous cell populations and long metabolic steady states such as developing seeds, cell suspensions, or cultured roots and root tips. It is now possible to describe these and other metabolically more dynamic tissues such as leaves in greater detail using novel methods in mass spectrometry, isotope labeling strategies, and transient labeling-based flux analyses. Such studies are necessary for a systems level description of plant function that more closely represents biological reality, and provides insights into the genes that will need to be modified as natural resources become ever more limited and environments change.

Plant phenomics and the need for physiological phenotyping across scales to narrow the genotype-to-phenotype knowledge gap

Plants are affected by complex genome×environment×management interactions which determine phenotypic plasticity as a result of the variability of genetic components. Whereas great advances have been made in the cost-efficient and high-throughput analyses of genetic information and non-invasive phenotyping, the large-scale analyses of the underlying physiological mechanisms lag behind. The external phenotype is determined by the sum of the complex interactions of metabolic pathways and intracellular regulatory networks that is reflected in an internal, physiological, and biochemical phenotype. These various scales of dynamic physiological responses need to be considered, and genotyping and external phenotyping should be linked to the physiology at the cellular and tissue level. A high-dimensional physiological phenotyping across scales is needed that integrates the precise characterization of the internal phenotype into high-throughput phenotyping of whole plants and canopies. By this means, complex traits can be broken down into individual components of physiological traits. Since the higher resolution of physiological phenotyping by 'wet chemistry' is inherently limited in throughput, high-throughput non-invasive phenotyping needs to be validated and verified across scales to be used as proxy for the underlying processes. Armed with this interdisciplinary and multidimensional phenomics approach, plant physiology, non-invasive phenotyping, and functional genomics will complement each other, ultimately enabling the in silico assessment of responses under defined environments with advanced crop models. This will allow generation of robust physiological predictors also for complex traits to bridge the knowledge gap between genotype and phenotype for applications in breeding, precision farming, and basic research.

Metabolomics-Inspired Insight into Developmental, Environmental and Genetic Aspects of Tomato Fruit Chemical Composition and Quality

Tomato was one of the first plant species to be evaluated using metabolomics and remains one of the best characterized, with tomato fruit being both an important source of nutrition in the human diet and a valuable model system for the development of fleshy fruits. Additionally, given the broad habitat range of members of the tomato clade and the extensive use of exotic germplasm in tomato genetic research, it represents an excellent genetic model system for understanding both metabolism per se and the importance of various metabolites in conferring stress tolerance. This review summarizes technical approaches used to characterize the tomato metabolome to date and details insights into metabolic pathway structure and regulation that have been obtained via analysis of tissue samples taken under different developmental or environmental circumstance as well as following genetic perturbation. Particular attention is paid to compounds of importance for nutrition or the shelf-life of tomatoes. We propose furthermore how metabolomics information can be coupled to the burgeoning wealth of genome sequence data from the tomato clade to enhance further our understanding of (i) the shifts in metabolic regulation occurring during development and (ii) specialization of metabolism within the tomato clade as a consequence of either adaptive evolution or domestication.

Intracellular and cell-to-apoplast compartmentation of carbohydrate metabolism

In most plants, carbohydrates represent the major energy store as well as providing the building blocks for essential structural polymers. Although the major pathways for carbohydrate biosynthesis, degradation, and transport are well characterized, several key steps have only recently been discovered. In addition, several novel minor metabolic routes have been uncovered in the past few years. Here we review current studies of plant carbohydrate metabolism detailing the expanding compendium of functionally characterized transport proteins as well as our deeper comprehension of more minor and conditionally activated metabolic pathways. We additionally explore the pertinent questions that will allow us to enhance our understanding of the response of both major and minor carbohydrate fluxes to changing cellular circumstances.

A Process-Based Model of TCA Cycle Functioning to Analyze Citrate Accumulation in Pre- and Post-Harvest Fruits

Citrate is one of the most important organic acids in many fruits and its concentration plays a critical role in organoleptic properties. The regulation of citrate accumulation throughout fruit development, and the origins of the phenotypic variability of the citrate concentration within fruit species remain to be clarified. In the present study, we developed a process-based model of citrate accumulation based on a simplified representation of the TCA cycle to predict citrate concentration in fruit pulp during the pre- and post-harvest stages. Banana fruit was taken as a reference because it has the particularity of having post-harvest ripening, during which citrate concentration undergoes substantial changes. The model was calibrated and validated on the two stages, using data sets from three contrasting cultivars in terms of citrate accumulation, and incorporated different fruit load, potassium supply, and harvest dates. The model predicted the pre and post-harvest dynamics of citrate concentration with fairly good accuracy for the three cultivars. The model suggested major differences in TCA cycle functioning among cultivars during post-harvest ripening of banana, and pointed to a potential role for NAD-malic enzyme and mitochondrial malate carriers in the genotypic variability of citrate concentration. The sensitivity of citrate accumulation to growth parameters and temperature differed among cultivars during post-harvest ripening. Finally, the model can be used as a conceptual basis to study citrate accumulation in fleshy fruits and may be a powerful tool to improve our understanding of fruit acidity.

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.

Cell-type specific metabolic profiling of Arabidopsis thaliana protoplasts as a tool for plant systems biology

Flow cytometry combined with cell sorting of protoplasts has previously been used successfully for transcript profiling of the Arabidopsis thaliana root. We have developed the technique further, and in this paper we present a robust and reliable method for metabolite profiling in specific cell types isolated from Arabidopsis roots. The method uses a combination of fluorescence-activated cell sorting and gas chromatography-time of flight-mass spectrometry analysis. Cortical and endodermal cells from the green fluorescent protein (GFP)-expressing enhancer trap line J0571 were analysed and compared with non-GFP-expressing cells and intact root tissue. Of the metabolites identified, several showed significant differences in concentration between cell types. Multivariate statistical analysis was used to compare metabolite patterns between cell and tissue types, showing that the patterns differed substantially. Isolation of specific cell populations combined with highly sensitive MS-analysis will be a powerful tool for future studies of plant metabolism, and can also be combined with transcript and protein profiling for in-depth analyses of cellular processes.

Moving toward a comprehensive map of central plant metabolism

Decades of intensive study have led to the discovery of the main pathways involved in central metabolism but only some of the pathways and regulatory networks in which they are embedded. In this review, we discuss techniques used to assemble these pathways into a systems biology framework that can enable accurate modeling of the response of central metabolism to changes, including ways to perturb metabolic systems and assemble the resulting data into a meaningful network. Critically, these networks are of such size and complexity that it is possible to derive them only if data from different groups can be comprehensively and meaningfully combined. We conclude that it is essential to establish common standards for the description of experimental conditions and data collection and to store this information in databases to which the whole community can contribute.

Protein import into plant mitochondria: signals, machinery, processing, and regulation

The majority of more than 1000 proteins present in mitochondria are imported from nuclear-encoded, cytosolically synthesized precursor proteins. This impressive feat of transport and sorting is achieved by the combined action of targeting signals on mitochondrial proteins and the mitochondrial protein import apparatus. The mitochondrial protein import apparatus is composed of a number of multi-subunit protein complexes that recognize, translocate, and assemble mitochondrial proteins into functional complexes. While the core subunits involved in mitochondrial protein import are well conserved across wide phylogenetic gaps, the accessory subunits of these complexes differ in identity and/or function when plants are compared with Saccharomyces cerevisiae (yeast), the model system for mitochondrial protein import. These differences include distinct protein import receptors in plants, different mechanistic operation of the intermembrane protein import system, the location and activity of peptidases, the function of inner-membrane translocases in linking the outer and inner membrane, and the association/regulation of mitochondrial protein import complexes with components of the respiratory chain. Additionally, plant mitochondria share proteins with plastids, i.e. dual-targeted proteins. Also, the developmental and cell-specific nature of mitochondrial biogenesis is an aspect not observed in single-celled systems that is readily apparent in studies in plants. This means that plants provide a valuable model system to study the various regulatory processes associated with protein import and mitochondrial biogenesis.

Knowledge-driven approaches for engineering complex metabolic pathways in plants

Plant metabolic pathways are complex and often feature multiple levels of regulation. Until recently, metabolic engineering in plants relied on the laborious testing of ad hoc modifications to achieve desirable changes in the metabolic profile. However, technological advances in data mining, modeling, multigene engineering and genome editing are now taking away much of the guesswork by allowing the impact of modifications to be predicted more accurately. In this review we discuss recent developments in knowledge-based metabolic engineering strategies, that is the gathering and mining of genomic, transcriptomic, proteomic and metabolomic data to generate models of metabolic pathways that help to define and refine optimal intervention strategies.

A network perspective on nitrogen metabolism from model to crop plants using integrated 'omics' approaches

Nitrogen (N), as an essential element in amino acids, nucleotides, and proteins, is a key factor in plant growth and development. Omics approaches such as metabolomics and transcriptomics have become a promising way to inspect complex network interactions in N metabolism and can be used for monitoring the uptake and regulation, translocation, and remobilization of N. In this review, the authors highlight recent progress in omics approaches, including transcript profiling using microarrays and deep sequencing, and show recent technical developments in metabolite profiling for N studies. Further, network analysis studies including network inference methods with correlations, information-theoretic measures, and a network concept to examine gene expression clusters in relation to N regulatory systems in plants are introduced, and integrating network inference methods and integrated networks using multiple omics data are discussed. Finally, this review summarizes recent omics application examples using metabolite and/or transcript profiling analysis to elucidate the regulation of N metabolism and signalling and the coordination of N and carbon metabolism in model plants (Arabidopsis and rice), crops (tomato, maize, and legumes), and trees (Populus).

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.

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.

Integrated network analysis and effective tools in plant systems biology

One of the ultimate goals in plant systems biology is to elucidate the genotype-phenotype relationship in plant cellular systems. Integrated network analysis that combines omics data with mathematical models has received particular attention. Here we focus on the latest cutting-edge computational advances that facilitate their combination. We highlight (1) network visualization tools, (2) pathway analyses, (3) genome-scale metabolic reconstruction, and (4) the integration of high-throughput experimental data and mathematical models. Multi-omics data that contain the genome, transcriptome, proteome, and metabolome and mathematical models are expected to integrate and expand our knowledge of complex plant metabolisms.