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

Computational simulation of large‐scale biochemical networks can be used to analyze and predict the metabolic behavior of an organism, such as a developing seed. Based on the biochemical literature, pathways databases and decision rules defining reaction directionality we reconstructed bna572, a stoichiometric metabolic network model representing Brassica napus seed storage metabolism. In the highly compartmentalized network about 25% of the 572 reactions are transport reactions interconnecting nine subcellular compartments and the environment. According to known physiological capabilities of developing B. napus embryos, four nutritional conditions were defined to simulate heterotrophy or photoheterotrophy, each in combination with the availability of inorganic nitrogen (ammonia, nitrate) or amino acids as nitrogen sources. Based on mathematical linear optimization the optimal solution space was comprehensively explored by flux variability analysis, thereby identifying for each reaction the range of flux values allowable under optimality. The range and variability of flux values was then categorized into flux variability types. Across the four nutritional conditions, approximately 13% of the reactions have variable flux values and 10–11% are substitutable (can be inactive), both indicating metabolic redundancy given, for example, by isoenzymes, subcellular compartmentalization or the presence of alternative pathways. About one‐third of the reactions are never used and are associated with pathways that are suboptimal for storage synthesis. Fifty‐seven reactions change flux variability type among the different nutritional conditions, indicating their function in metabolic adjustments. This predictive modeling framework allows analysis and quantitative exploration of storage metabolism of a developing B. napus oilseed.     

Proteomes of the barley aleurone layer: A model system for plant signalling and protein secretion

The cereal aleurone layer is of major importance due to its nutritional properties as well as its central role in seed germination and industrial malting. Cereal seed germination involves mobilisation of storage reserves in the starchy endosperm to support seedling growth. In response to gibberellic acid produced by the embryo, the aleurone layer synthesises hydrolases that are secreted to the endosperm for the degradation of storage products. The barley aleurone layer can be separated from the other seed tissues and maintained in culture, allowing the study of the effect of added signalling molecules in an isolated system. These properties have led to its use as a model system for the study of plant signalling and germination. More recently, proteome analysis of the aleurone layer has provided new insight into this unique tissue including identification of plasma membrane proteins and targeted analysis of germination-related changes and the thioredoxin system. Here, analysis of intracellular and secreted proteomes reveals features of the aleurone layer system that makes it promising for investigations of plant protein secretion mechanisms.     

Energy state and its control on seed development: starch accumulation is associated with high ATP and steep oxygen gradients within barley grains

The role of oxygen and energy state in development and storage activity of cereal grains is an important issue, but has remained largely uninvestigated due to the lack of appropriate analytical methods. Metabolic profiling, bioluminescence-based in situ imaging of ATP, and oxygen-sensitive microsensors were combined here to investigate barley seed development. For the first time temporal and spatial maps of O2 and ATP distribution in cereal grains were determined and related to the differentiation pattern. Steep O2 gradients were demonstrated and strongly hypoxic regions were detected within the caryopsis (<0.1% of atmospheric saturation). Growing lateral and peripheral regions of endosperm remained well-supplied with O2 due to pericarp photosynthesis. ATP distribution in the developing grain was coupled to endosperm differentiation. High ATP concentrations were associated with the local onset of starch storage within endosperm, while low ATP overlapped with the hypoxic regions. Temporally, the building of steep gradients in ATP coincided with overall elevating metabolite levels, specific changes in the metabolite profiles (glycolysis and citrate cycle), and channelling of metabolic fluxes towards storage (increase of starch accumulation rate). These findings implicate an inhomogenous spatial arrangement of metabolic activity within the caryopsis. It is suggested that the local onset of starch storage is coupled with the accumulation of ATP and elevated metabolic activity. Thus, the ATP level reflects the metabolic state of storage tissue. On the basis of these findings, a hypothetical model for the regulation of starch storage in barley seeds is proposed.     

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.     

Mitochondrial metabolism in developing embryos of Brassica napus

The metabolism of developing plant seeds is directed toward transforming primary assimilatory products (sugars and amino acids) into seed storage compounds. To understand the role of mitochondria in this metabolism, metabolic fluxes were determined in developing embryos of Brassica napus. After labeling with [1,2-(13)C2]glucose + [U-(13)C6]glucose, [U-(13)C3]alanine, [U-(13)C5]glutamine, [(15)N]alanine, (amino)-[(15)N]glutamine, or (amide)-[(15)N]glutamine, the resulting labeling patterns in protein amino acids and in fatty acids were analyzed by gas chromatography-mass spectrometry. Fluxes through mitochondrial metabolism were quantified using a steady state flux model. Labeling information from experiments using different labeled substrates was essential for model validation and reliable flux estimation. The resulting flux map shows that mitochondrial metabolism in these developing seeds is very different from that in either heterotrophic or autotrophic plant tissues or in most other organisms: (i) flux around the tricarboxylic acid cycle is absent and the small fluxes through oxidative reactions in the mitochondrion can generate (via oxidative phosphorylation) at most 22% of the ATP needed for biosynthesis; (ii) isocitrate dehydrogenase is reversible in vivo; (iii) about 40% of mitochondrial pyruvate is produced by malic enzyme rather than being imported from the cytosol; (iv) mitochondrial flux is largely devoted to providing precursors for cytosolic fatty acid elongation; and (v) the uptake of amino acids rather than anaplerosis via PEP carboxylase determines carbon flow into storage proteins.     

Phloem import and storage metabolism are highly coordinated by the low oxygen concentrations within developing wheat seeds

We studied the influence of the internal oxygen concentration in seeds of wheat (Triticum aestivum) on storage metabolism and its relation to phloem import of nutrients. Wheat seeds that were developing at ambient oxygen (21%) were found to be hypoxic (2.1%). Altering the oxygen supply by decreasing or increasing the external oxygen concentration induced parallel changes in the internal oxygen tension. However, the decrease in internal concentration was proportionally less than the reduction in external oxygen. This indicates that decreasing the oxygen supply induces short-term adaptive responses to reduce oxygen consumption of the seeds. When external oxygen was decreased to 8%, internal oxygen decreased to approximately 0.5% leading to a decrease in energy production via respiration. Conversely, increasing the external oxygen concentration above ambient levels increased the oxygen content as well as the energy status of the seeds, indicating that under normal conditions the oxygen supply is strongly limiting for energy metabolism in developing wheat seeds. The intermediate metabolites of seed storage metabolism were not substantially affected when oxygen was either increased or decreased. However, at subambient external oxygen concentrations (8%) the metabolic flux of carbon into starch and protein, measured by injecting (14)C-Suc into the seeds, was reduced by 17% and 32%, respectively, whereas no significant effect was observed at superambient (40%) oxygen. The observed decrease in biosynthetic fluxes to storage compounds is suggested to be part of an adaptive response to reduce energy consumption preventing excessive oxygen consumption when oxygen supply is limited. Phloem transport toward ears exposed to low (8%) oxygen was significantly reduced within 1 h, whereas exposing ears to elevated oxygen (40%) had no significant effect. This contrasts with the situation where the distribution of assimilates has been modified by removing the lower source leaves from the plant, resulting in less assimilates transported to the ear in favor of transport to the lower parts of the plant. Under these conditions, with two strongly competing sinks, elevated oxygen (40%) did lead to a strong increase in phloem transport to the ear. The results show that sink metabolism is affected by the prevailing low oxygen concentrations in developing wheat seeds, determining the import rate of assimilates via the phloem.     

An in silico compartmentalized metabolic model of Brassica napus enables the systemic study of regulatory aspects of plant central metabolism

Biochemical network reconstructions represent valuable tools for the computational metabolic modeling of organisms that present a great biotechnological interest. An in silico multi-compartmental model of the central metabolism of the plant Brassica napus (Rapeseed) was constructed, aiming to investigate the metabolic properties of the Brassicaceae family. This family comprises many plants with major importance for the energy and nutrition sector, including the model plant Arabidopsis thaliana. The model utilized as objective function to be subsequently optimized, the biomass production of rapeseed developing embryos, which are characterized by a very high, oil content, up to 60% of biomass weight. In order to study global network properties of seed metabolism, various methods were employed, like Flux Balance Analysis, Principal Component Analysis of the flux space and reaction deletion studies, which simulate the effect of gene knock-out experiments. The model successfully simulated seed growth during the stage of oil accumulation and provided insight, regarding certain aspects of network plasticity, with the emphasis given in lipid biosynthesis regulation.     

An assessment of the biotechnological use of hemoglobin modulation in cereals

Non‐symbiotic hemoglobin (nsHb) genes are ubiquitous in plants, but their biological functions have mostly been studied in model plant species rather than in crops. nsHb influences cell signaling and metabolism by modulating the levels of nitric oxide (NO). Class 1 nsHb is upregulated under hypoxia and is involved in various biotic and abiotic stress responses. Ectopic overexpression of nsHb in Arabidopsis thaliana accelerates development, whilst targeted overexpression in seeds can increase seed yield. Such observations suggest that manipulating nsHb could be a valid biotechnological target. We studied the effects of overexpression of class 1 nsHb in the monocotyledonous crop plant barley (Hordeum vulgare cv. Golden Promise). nsHb was shown to be involved in NO metabolism in barley, as ectopic overexpression reduced the amount of NO released during hypoxia. Further, as in Arabidopsis, nsHb overexpression compromised basal resistance toward pathogens in barley. However, unlike Arabidopsis, nsHb ectopic overexpression delayed growth and development in barley, and seed specific overexpression reduced seed yield. Thus, nsHb overexpression in barley does not seem to be an efficient strategy for increasing yield in cereal crops. These findings highlight the necessity for using actual crop plants rather than laboratory model plants when assessing the effects of biotechnological approaches to crop improvement.     

Metabolic pathway analysis of Scheffersomyces (Pichia) stipitis: effect of oxygen availability on ethanol synthesis and flux distributions

Elementary mode analysis (EMA) identifies all possible metabolic states of the cell metabolic network. Investigation of these states can provide a detailed insight into the underlying metabolism in the cell. In this study, the flux states of Scheffersomyces (Pichia) stipitis metabolism were examined. It was shown that increasing oxygen levels led to a decrease of ethanol synthesis. This trend was confirmed by experimental evaluation of S. stipitis in glucose-xylose fermentation. The oxygen transfer rate for an optimal ethanol production was 1.8 mmol/l/h, which gave the ethanol yield of 0.40 g/g and the ethanol productivity of 0.25 g/l/h. For a better understanding of the cell's regulatory mechanism in response to oxygenation levels, EMA was used to examine metabolic flux patterns under different oxygen levels. Up- and downregulation of enzymes in the network during the change of culturing condition from oxygen limitation to oxygen sufficiency were identified. The results indicated the flexibility of S. stipitis metabolism to cope with oxygen availability. In addition, relevant genetic targets towards improved ethanol-producing strains under all oxygenation levels were identified. These targeted genes limited the metabolic functionality of the cell to function according to the most efficient ethanol synthesis pathways. The presented approach is promising and can contribute to the development of culture optimization and strain engineers for improved lignocellulosic ethanol production by S. stipitis.     

Transgenic barley expressing a fungal xylanase gene in the endosperm of the developing grains

The feasibility of producing plant cell wall polysaccharide-hydrolysing feed enzymes in the endosperm of barley grain was investigated. The coding region of a modified xylanase gene (xynA) from the rumen fungus, Neocallimastix patriciarum, linked with an endosperm-specific promoter from cereal storage protein genes was introduced into barley by Agrobacterium-mediated transformation. Twenty-four independently transformed barley lines with the xylanase gene were produced and analysed. The fungal xylanase was produced in the developing endosperm under the control of either the rice glutelin B-1 (GluB-1) or barley B1 hordein (Hor2-4) promoter. The rice GluB-1 promoter provided an apparently higher expression level of recombinant proteins in barley grain than the barley Hor2-4 promoter in both transient and stable expression experiments. In particular, the mean value for the fungal xylanase activity driven by the GluB-1 promoter in the mature grains of transgenic barley was more than twice that with the Hor2-4 promoter. Expression of the xylanase transgene under these endosperm-specific promoters was not observed in the leaf, stem and root tissues. Accumulation of the fungal xylanase in the developing grains of transgenic barley followed the pattern of storage protein deposition. The xylanase was stably maintained in the grain during grain maturation and desiccation and post-harvest storage. These results indicate that the cereal grain expression system may provide an economic means for large scale production of feed enzymes in the future.     

Seed Architecture Shapes Embryo Metabolism in Oilseed Rape

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

Physiological,Biochemical and Molecular Responses of Barley Seedlings to Aluminum Stress

Barley (Hordeum vulgare L.) is one of the most Aluminum (Al) sensitive cereal species. In this study, the physiological, biochemical, and molecular response of barley seedlings to Al treatment was examined to gain insight into Al response and tolerance mechanisms. The results showed that superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activity were inhibited to different degrees following Al exposure. The MDA content also significantly increased with increasing Al concentrations. SRAP results indicated significant differences between Al treatments and controls in terms of SRAP profile, and the genomic template stability (GTS) decreased with increasing Al concentration and duration. These integrative results help to elucidate the underlying mechanisms that the barley response to Al toxicity.     

Efficacy of imidacloprid for control of cereal leaf beetle (Coleoptera: Chrysomelidae) in barley

The toxicity of imidacloprid to the cereal leaf beetle, Oulema melanopus (L.), was measured under laboratory and field conditions. Insect mortality and plant damage were determined from artificial and natural infestations of O. melanopus applied to various growth stages of barley. All rates of imidacloprid formulated and applied as a seed treatment caused >90% mortality to cereal leaf beetle larvae when barley was infested with eggs at the 4-leaf stage, but were ineffective when barley was infested with eggs at the early tillering or flag-leaf stages of barley. This window of susceptibility influenced results obtained in field trials where peak larval emergence did not occur until the early tillering stage of barley. The resulting mortality in plants from treated seeds never exceeded 40% in the field. Foliar imidacloprid, however, caused >90% mortality in the field, and may be another option in the management of the cereal leaf beetle.     

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

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

Dynamic flux balance analysis of the metabolism of Saccharomyces cerevisiae during the shift from fully respirative or respirofermentative metabolic states to anaerobiosis

Dynamic flux balance analysis was utilized to simulate the metabolic behaviour of initially fully respirative and respirofermentative steady-state cultures of Saccharomyces?cerevisiae during sudden oxygen depletion. The hybrid model for the dynamic flux balance analysis included a stoichiometric genome-scale metabolic model as a static part and dynamic equations for the uptake of glucose and the cessation of respirative metabolism. The yeast consensus genome-scale metabolic model [Herrg?rd MJ et?al. (2008) Nat Biotechnol26, 1155-1160; Dobson PD et?al. (2010) BMC Syst Biol4, 145] was refined with respect to oxygen-dependent energy metabolism and further modified to reflect S.?cerevisiae anabolism in the absence of oxygen. Dynamic flux balance analysis captured well the essential features of the dynamic metabolic behaviour of S.?cerevisiae during adaptation to anaerobiosis. Modelling and simulation enabled the identification of short time-scale flux distribution dynamics under the transition to anaerobic metabolism, during which the specific growth rate was reduced, as well as longer time-scale process dynamics when the specific growth rate recovered. Expression of the metabolic genes was set into the context of the identified dynamics. Metabolic gene expression responses associated with the specific growth rate and with the cessation of respirative metabolism were distinguished.     

Cloning and characterization of SOC1 homologs in barley (Hordeum vulgare) and their expression during seed development and in response to vernalization

A number of genes are involved in the vernalization pathway, such as VRN1, VRN2 and VRN3/FT1, whose function has been studied in barley and wheat. However, the function of the flowering and vernalization integrator SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) has not been well studied in Triticeae, and particularly in barley. Herein, we cloned and characterized two barley SOC1-like homologs, HvSOC1-like1 and HvSOC1-like2. Primary sequence analysis of the predicted HvSOC1-like1 and HvSOC1-like2 proteins showed that they are members of the type II MADS-box protein family. Phylogenetic analysis placed the predicted proteins with other SOC1 and SOC1-like proteins from different species neighboring those from other cereal plant species. Primary and secondary structures of the predicted proteins are conserved to each other and more distant to the recently identified barley ODDSOC1 proteins. Genomic organization of HvSOC1-like1 is very similar to the Arabidopsis and Brachypodium SOC1 genes and localized in highly syntenic chromosomal regions. Regulatory cis-acting elements detected in the HvSOC1-like1 promoter include the CArG-box, implicated in the regulation of SOC1 expression in Arabidopsis. Both HvSOC1-like1 and HvSOCI-like2 are expressed in vegetative and reproductive tissues and at different stages of seed development. Both are upregulated in a particular seed developmental stage suggesting their possible implication in seed development. Furthermore, HvSOC1-like1 was induced in two winter barley cultivars after vernalization treatment pointing to its probable involvement in the vernalization process. The study of the SOC1 genes reported here opens the way for a better understanding of both the vernalization process and seed development and germination in this important cereal crop.     

Metabolic redesign of vitamin E biosynthesis in plants for tocotrienol production and increased antioxidant content

Tocotrienols are the primary form of vitamin E in seeds of most monocot plants, including cereals such as rice and wheat. As potent antioxidants, tocotrienols contribute to the nutritive value of cereal grains in human and livestock diets. cDNAs encoding homogentisic acid geranylgeranyl transferase (HGGT), which catalyzes the committed step of tocotrienol biosynthesis, were isolated from barley, wheat and rice seeds. Transgenic expression of the barley HGGT in Arabidopsis thaliana leaves resulted in accumulation of tocotrienols, which were absent from leaves of nontransformed plants, and a 10- to 15-fold increase in total vitamin E antioxidants (tocotrienols plus tocopherols). Overexpression of the barley HGGT in corn seeds resulted in an increase in tocotrienol and tocopherol content of as much as six-fold. These results provide insight into the genetic basis for tocotrienol biosynthesis in plants and demonstrate the ability to enhance the antioxidant content of crops by introduction of an enzyme that redirects metabolic flux.     

Involvement of low molecular mass soluble acyl-CoA-binding protein in seed oil biosynthesis

Acyl-CoA-binding protein (ACBP), a low molecular mass (m) (～ 10 kDa) soluble protein ubiquitous in eukaryotes, plays an important housekeeping role in lipid metabolism by maintaining the intracellular acyl-CoA pool. ACBP is involved in lipid biosynthesis and transport, gene expression, and membrane biogenesis. In plants, low m ACBP and high m ACBPs participate in response mechanisms to biotic and abiotic factors, acyl-CoA transport in phloem, and biosynthesis of structural and storage lipids. In light of current research on the modification of seed oil, insight into mechanisms of substrate trafficking within lipid biosynthetic pathways is crucial for developing rational strategies for the production of specialty oils with the desired alterations in fatty acid composition. In this review, we summarize our knowledge of plant ACBPs with emphasis on the role of low m ACBP in seed oil biosynthesis, based on in vitro studies and analyses of transgenic plants. Future prospects and possible applications of low m ACBP in seed oil modification are discussed.