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 highlight 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.

SUMMARYComputational 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.

SUMMARYPlant oils are an important renewable resource, and seed oil content is a key agronomical trait that is in part controlled by the metabolic processes within developing seeds. A large-scale model of cellular metabolism in developing embryos of Brassica napus (bna572) was used to predict biomass formation and to analyze metabolic steady states by flux variability analysis under different physiological conditions. Predicted flux patterns are highly correlated with results from prior 13 C metabolic flux analysis of B. napus developing embryos. Minor differences from the experimental results arose because bna572 always selected only one sugar and one nitrogen source from the available alternatives, and failed to predict the use of the oxidative pentose phosphate pathway. Flux variability, indicative of alternative optimal solutions, revealed alternative pathways that can provide pyruvate and NADPH to plastidic fatty acid synthesis. The nutritional values of different medium substrates were compared based on the overall carbon conversion efficiency (CCE) for the biosynthesis of biomass. Although bna572 has a functional nitrogen assimilation pathway via glutamate synthase, the simulations predict an unexpected role of glycine decarboxylase operating in the direction of NH 4 + assimilation. Analysis of the light-dependent improvement of carbon economy predicted two metabolic phases. At very low light levels small reductions in CO 2 efflux can be attributed to enzymes of the tricarboxylic acid cycle (oxoglutarate dehydrogenase, isocitrate dehydrogenase) and glycine decarboxylase. At higher light levels relevant to the 13 C flux studies, ribulose-1,5-bisphosphate carboxylase activity is predicted to account fully for the light-dependent changes in carbon balance.

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.

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 acetylcoenzyme A for cytosol-localized fatty acid elongation. The phase close to equal oil and protein fractions included an unexpected pathway bypassing a-ketoglutarate-oxidizing steps in the tricarboxylic acid cycle. The in vivo relevance of the findings is discussed based on literature on seed storage metabolism.