Mechanisms of metabolic adaptation in the duckweed Lemna gibba: an integrated metabolic, transcriptomic and flux analysis.

BACKGROUND: Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions. RESULTS: To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts. CONCLUSIONS: Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.

Cellular Plasticity in Response to Suppression of Storage Proteins in the Brassica napus Embryo.

The tradeoff between protein and oil storage in oilseed crops has been tested here in oilseed rape (Brassica napus) by analyzing the effect of suppressing key genes encoding protein storage products (napin and cruciferin). The phenotypic outcomes were assessed using NMR and mass spectrometry imaging, microscopy, transcriptomics, proteomics, metabolomics, lipidomics, immunological assays, and flux balance analysis. Surprisingly, the profile of storage products was only moderately changed in RNA interference transgenics. However, embryonic cells had undergone remarkable architectural rearrangements. The suppression of storage proteins led to the elaboration of membrane stacks enriched with oleosin (sixfold higher protein abundance) and novel endoplasmic reticulum morphology. Protein rebalancing and amino acid metabolism were focal points of the metabolic adjustments to maintain embryonic carbon/nitrogen homeostasis. Flux balance analysis indicated a rather minor additional demand for cofactors (ATP and NADPH). Thus, cellular plasticity in seeds protects against perturbations to its storage capabilities and, hence, contributes materially to homeostasis. This study provides mechanistic insights into the intriguing link between lipid and protein storage, which have implications for biotechnological strategies directed at improving oilseed crops.

Waterlogging tolerance rendered by oxylipin-mediated metabolic reprogramming in Arabidopsis.

Environmental stresses induce production of oxylipins synthesized by the two main biosynthetic branches, allene oxide synthase (AOS) and hydroperoxide lyase (HPL). Here, we investigate how waterlogging-mediated alteration of AOS- and HPL-derived metabolic profile results in modulation of central metabolism and ultimately enhanced tolerance to this environmental stress in Arabidopsis thaliana. Waterlogging leads to increased levels of AOS- and HPL-derived metabolites, and studies of genotypes lacking either one or both branches further support the key function of these oxylipins in waterlogging tolerance. Targeted quantitative metabolic profiling revealed oxylipin-dependent alterations in selected primary metabolites, and glycolytic and citric acid cycle intermediates, as well as a prominent shift in sucrose cleavage, hexose activation, the methionine salvage pathway, shikimate pathway, antioxidant system, and energy metabolism in genotypes differing in the presence of one or both functional branches of the oxylipin biosynthesis pathway. Interestingly, despite some distinct metabolic alterations caused specifically by individual branches, overexpression of HPL partially or fully alleviates the majority of altered metabolic profiles observed in AOS-depleted lines. Collectively, these data identify the key role of AOS- and HPL-derived oxylipins in altering central metabolism, and further provide a metabolic platform targeted at identification of gene candidates for enhancing plant tolerance to waterlogging.

Vacuolar processing enzyme 4 contributes to maternal control of grain size in barley by executing programmed cell death in the pericarp.

The angiosperm embryo and endosperm are limited in space because they grow inside maternal seed tissues. The elimination of cell layers of the maternal seed coat by programmed cell death (PCD) could provide space and nutrition to the filial organs. Using the barley (Hordeum vulgare L.) seed as a model, we elucidated the role of vacuolar processing enzyme 4 (VPE4) in cereals by using an RNAi approach and targeting the enzymatic properties of the recombinant protein. A comparative characterization of transgenic versus wild-type plants included transcriptional and metabolic profiling, flow cytometry, histology and nuclear magnetic imaging of grains. The recombinant VPE4 protein exhibited legumain and caspase-1 properties in vitro. Pericarp disintegration was delayed in the transgenic grains. Although the VPE4 gene and enzymatic activity was decreased in the early developing pericarp, storage capacity and the size of the endosperm and embryo were reduced in the mature VPE4-repressed grains. The persistence of the pericarp in the VPE4-affected grains constrains endosperm and embryo growth and leads to transcriptional reprogramming, perturbations in signalling and adjustments in metabolism. We conclude that VPE4 expression executes PCD in the pericarp, which is required for later endosperm filling, and argue for a role of PCD in maternal control of seed size in cereals.

Quantitative Multilevel Analysis of Central Metabolism in Developing Oilseeds of Oilseed Rape during in Vitro Culture.

Seeds provide the basis for many food, feed, and fuel products. Continued increases in seed yield, composition, and quality require an improved understanding of how the developing seed converts carbon and nitrogen supplies into storage. Current knowledge of this process is often based on the premise that transcriptional regulation directly translates via enzyme concentration into flux. In an attempt to highlight metabolic control, we explore genotypic differences in carbon partitioning for in vitro cultured developing embryos of oilseed rape (Brassica napus). We determined biomass composition as well as 79 net fluxes, the levels of 77 metabolites, and 26 enzyme activities with specific focus on central metabolism in nine selected germplasm accessions. Overall, we observed a tradeoff between the biomass component fractions of lipid and starch. With increasing lipid content over the spectrum of genotypes, plastidic fatty acid synthesis and glycolytic flux increased concomitantly, while glycolytic intermediates decreased. The lipid/starch tradeoff was not reflected at the proteome level, pointing to the significance of (posttranslational) metabolic control. Enzyme activity/flux and metabolite/flux correlations suggest that plastidic pyruvate kinase exerts flux control and that the lipid/starch tradeoff is most likely mediated by allosteric feedback regulation of phosphofructokinase and ADP-glucose pyrophosphorylase. Quantitative data were also used to calculate in vivo mass action ratios, reaction equilibria, and metabolite turnover times. Compounds like cyclic 3',5'-AMP and sucrose-6-phosphate were identified to potentially be involved in so far unknown mechanisms of metabolic control. This study provides a rich source of quantitative data for those studying central metabolism.

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.

Combined noninvasive imaging and modeling approaches reveal metabolic compartmentation in the barley endosperm.

The starchy endosperm of cereals is a priori taken as a metabolically uniform tissue. By applying a noninvasive assay based on (13)C/(1)H-magnetic resonance imaging (MRI) to barley (Hordeum vulgare) grains, we uncovered metabolic compartmentation in the endosperm. (13)C-Suc feeding during grain filling showed that the primary site of Ala synthesis was the central region of the endosperm, the part of the caryopsis experiencing the highest level of hypoxia. Region-specific metabolism in the endosperm was characterized by flux balance analysis (FBA) and metabolite profiling. FBA predicts that in the central region of the endosperm, the tricarboxylic acid cycle shifts to a noncyclic mode, accompanied by elevated glycolytic flux and the accumulation of Ala. The metabolic compartmentation within the endosperm is advantageous for the grain's carbon and energy economy, with a prominent role being played by Ala aminotransferase. An investigation of caryopses with a genetically perturbed tissue pattern demonstrated that Ala accumulation is a consequence of oxygen status, rather than being either tissue specific or dependent on the supply of Suc. Hence the (13)C-Ala gradient can be used as an in vivo marker for hypoxia. The combination of MRI and metabolic modeling offers opportunities for the noninvasive analysis of metabolic compartmentation in plants.

Environmental stresses of field growth allow cinnamyl alcohol dehydrogenase-deficient Nicotiana attenuata plants to compensate for their structural deficiencies.

The organized lignocellulosic assemblies of cell walls provide the structural integrity required for the large statures of terrestrial plants. Silencing two CINNAMYL ALCOHOL DEHYDROGENASE (CAD) genes in Nicotiana attenuata produced plants (ir-CAD) with thin, red-pigmented stems, low CAD and sinapyl alcohol dehydrogenase activity, low lignin contents, and rubbery, structurally unstable stems when grown in the glasshouse (GH). However, when planted into their native desert habitat, ir-CAD plants produced robust stems that survived wind storms as well as the wild-type plants. Despite efficient silencing of NaCAD transcripts and enzymatic activity, field-grown ir-CAD plants had delayed and restricted spread of red stem pigmentation, a color change reflecting blocked lignification by CAD silencing, and attained wild-type-comparable total lignin contents. The rubbery GH phenotype was largely restored when field-grown ir-CAD plants were protected from wind, herbivore attack, and ultraviolet B exposure and grown in restricted rooting volumes; conversely, it was lost when ir-CAD plants were experimentally exposed to wind, ultraviolet B, and grown in large pots in growth chambers. Transcript and liquid chromatography-electrospray ionization-time-of-flight analysis revealed that these environmental stresses enhanced the accumulation of various phenylpropanoids in stems of field-grown plants; gas chromatography-mass spectrometry and nuclear magnetic resonance analysis revealed that the lignin of field-grown ir-CAD plants had GH-grown comparable levels of sinapaldehyde and syringaldehyde cross-linked into their lignins. Additionally, field-grown ir-CAD plants had short, thick stems with normal xylem element traits, which collectively enabled field-grown ir-CAD plants to compensate for the structural deficiencies associated with CAD silencing. Environmental stresses play an essential role in regulating lignin biosynthesis in lignin-deficient plants.

A dual role of tobacco hexokinase 1 in primary metabolism and sugar sensing.

Hexokinase (HXK) is present in all virtually living organisms and is central to carbohydrate metabolism catalysing the ATP-dependent phosphorylation of hexoses. In plants, HXKs are supposed to act as sugar sensors and/or to interact with other enzymes directly supplying metabolic pathways such as glycolysis, the nucleotide phosphate monosaccharide (NDP-glucose) pathway and the pentose phosphate pathway. We identified nine members of the tobacco HXK gene family and observed that among RNAi lines of these nine NtHXKs, only RNAi lines of NtHXK1 showed an altered phenotype, namely stunted growth and leaf chlorosis. NtHXK1 was also the isoform with highest relative expression levels among all NtHXKs. GFP-tagging and immunolocalization indicated that NtHXK1 is associated with mitochondrial membranes. Overexpression of NtHXK1 resulted in elevated glucose phosphorylation activity in leaf extracts or chloroplasts. Moreover, NtHXK1 was able to complement the glucose-insensitive Arabidopsis mutant gin2-1 suggesting that NtHXK1 can take over glucose sensing functions. RNAi lines of NtHXK1 showed severely damaged leaf and chloroplast structure, coinciding with an excess accumulation of starch. We conclude that NtHXK1 is not only essential for maintaining glycolytic activity during respiration but also for regulating starch turnover, especially during the night.

Dynamic ¹³C/¹ H NMR imaging uncovers sugar allocation in the living seed.

Seed growth and accumulation of storage products relies on the delivery of sucrose from the maternal to the filial tissues. The transport route is hidden inside the seed and has never been visualized in vivo. Our approach, based on high-field nuclear magnetic resonance and a custom made (13)C/(1) H double resonant coil, allows the non-invasive imaging and monitoring of sucrose allocation within the seed. The new technique visualizes the main stream of sucrose and determines its velocity during the grain filling in barley (Hordeum vulgare L.). Quantifiable dynamic images are provided, which allow observing movement of (13)C-sucrose at a sub-millimetre level of resolution. The analysis of genetically modified barley grains (Jekyll transgenic lines, seg8 and Risø13 mutants) demonstrated that sucrose release via the nucellar projection towards the endosperm provides an essential mean for the control of seed growth by maternal organism. The sucrose allocation was further determined by structural and metabolic features of endosperm. Sucrose monitoring was integrated with an in silico flux balance analysis, representing a powerful platform for non-invasive study of seed filling in crops.

R2R3-NaMYB8 regulates the accumulation of phenylpropanoid-polyamine conjugates, which are essential for local and systemic defense against insect herbivores in Nicotiana attenuata.

Although phenylpropanoid-polyamine conjugates (PPCs) occur ubiquitously in plants, their biological roles remain largely unexplored. The two major PPCs of Nicotiana attenuata plants, caffeoylputrescine (CP) and dicaffeoylspermidine, increase dramatically in local and systemic tissues after herbivore attack and simulations thereof. We identified NaMYB8, a homolog of NtMYBJS1, which in BY-2 cells regulates PPC biosynthesis, and silenced its expression by RNA interference in N. attenuata (ir-MYB8), to understand the ecological role(s) of PPCs. The regulatory role of NaMYB8 in PPC biosynthesis was validated by a microarray analysis, which revealed that transcripts of several key biosynthetic genes in shikimate and polyamine metabolism accumulated in a NaMYB8-dependent manner. Wild-type N. attenuata plants typically contain high levels of PPCs in their reproductive tissues; however, NaMYB8-silenced plants that completely lacked CP and dicaffeoylspermidine showed no changes in reproductive parameters of the plants. In contrast, a defensive role for PPCs was clear; both specialist (Manduca sexta) and generalist (Spodoptera littoralis) caterpillars feeding on systemically preinduced young stem leaves performed significantly better on ir-MYB8 plants lacking PPCs compared with wild-type plants expressing high levels of PPCs. Moreover, the growth of M. sexta caterpillars was significantly reduced when neonates were fed ir-MYB8 leaves sprayed with synthetic CP, corroborating the role of PPCs as direct plant defense. The spatiotemporal accumulation and function of PPCs in N. attenuata are consistent with the predictions of the optimal defense theory: plants preferentially protect their most fitness-enhancing and vulnerable parts, young tissues and reproductive organs, to maximize their fitness.

Polymorphism in jasmonate signaling partially accounts for the variety of volatiles produced by Nicotiana attenuata plants in a native population.

Herbivore- and jasmonate-induced volatile organic compounds (VOCs), which mediate indirect defense, must provide reliable information for predators that frequently learn to associate their release with feeding herbivores. Yet little is known about variation of these cues within populations of native plants, on a scale encountered by predators. We examined variation in herbivore-elicited VOC emissions and patterns of herbivore-induced jasmonate signaling from accessions of Nicotiana attenuata co-occurring in a native population. VOC emissions elicited by herbivore oral secretions (OS) and by methyl jasmonate (MJ) were characterized using gas chromatography-mass spectrometry (GC-MS), high-resolution two-dimensional gas chromatography-time-of-flight mass spectrometry (GCxGC-ToF-MS) and micro-hydrolysis and micro-hydrogenation reactions. Accessions varied in emissions of abundant (trans-alpha-bergamotene, alpha-duprezianene, trans-beta-ocimene, and cis-3-hexenol) and total detectable VOCs, as well as the accumulation of jasmonates, the jasmonate antagonist salicylic acid (SA), abscisic acid (ABA) and jasmonate signaling-related transcripts after OS elicitation. Yet MJ treatment exacerbated differences in VOC emission, suggesting that much variation in VOC emission is caused by processes downstream of jasmonate signaling. Co-occurring N. attenuata plants emit different VOCs following simulated herbivore elicitation as a result in part of differences in jasmonate production and responsiveness, which could reduce the effectiveness of induced indirect defense.

The Role of Persulfide Metabolism During Arabidopsis Seed Development Under Light and Dark Conditions.

The sulfur dioxygenase ETHE1 oxidizes persulfides in the mitochondrial matrix and is involved in the degradation of L-cysteine and hydrogen sulfide. ETHE1 has an essential but as yet undefined function in early embryo development of Arabidopsis thaliana. In leaves, ETHE1 is strongly induced by extended darkness and participates in the use of amino acids as alternative respiratory substrates during carbohydrate starvation. Thus, we tested the effect of darkness on seed development in an ETHE1 deficient mutant in comparison to the wild type. Since ETHE1 knock-out is embryo lethal, the knock-down line ethe1-1 with about 1% residual sulfur dioxygenase activity was used for this study. We performed phenotypic analysis, metabolite profiling and comparative proteomics in order to investigate the general effect of extended darkness on seed metabolism and further define the specific function of the mitochondrial sulfur dioxygenase ETHE1 in seeds. Shading of the siliques had no morphological effect on embryogenesis in wild type plants. However, the developmental delay that was already visible in ethe1-1 seeds under control conditions was further enhanced in the darkness. Dark conditions strongly affected seed quality parameters of both wild type and mutant plants. The effect of ETHE1 knock-down on amino acid profiles was clearly different from that found in leaves indicating that in seeds persulfide oxidation interacts with alanine and glycine rather than branched-chain amino acid metabolism. Sulfur dioxygenase deficiency led to defects in endosperm development possibly due to alterations in the cellularization process. In addition, we provide evidence for a potential role of persulfide metabolism in abscisic acid (ABA) signal transduction in seeds. We conclude that the knock-down of ETHE1 causes metabolic re-arrangements in seeds that differ from those in leaves. Putative mechanisms that cause the aberrant endosperm and embryo development are discussed.

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.

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.

A novel procedure for the quantitative analysis of metabolites, storage products and transcripts of laser microdissected seed tissues of Brassica napus.

BACKGROUND: The biology of the seed is complicated by the extensive non-homogeneity (spatial gradients) in gene expression, metabolic conversions and storage product accumulation. The detailed understanding of the mechanisms underlying seed growth and storage therefore requires the development of means to obtain tissue-specific analyses. This approach also represents an important priority in the context of seed biotechnology. RESULTS: We provide a guideline and detailed procedures towards the quantitative analysis of laser micro-dissected (LM) tissues in oilseed rape (Brassica napus). This includes protocols for laser microdissection of the seed, and the subsequent extraction and quantitative analysis of lipids, starch and metabolites (sugars, sugar phosphates, nucleotides, amino acids, intermediates of glycolysis and citric acid cycle). We have also developed a protocol allowing the parallel analysis of the transcriptome using Brassica-specific microarrays. Some data are presented regarding the compartmentation of metabolites within the oilseed rape embryo. CONCLUSION: The described methodology allows for the rapid, combined analysis of metabolic intermediates, major storage products and transcripts in a tissue-specific manner. The protocols are robust for oilseed rape, and should be readily adjustable for other crop species. The suite of methods applied to LM tissues represents an important step in the context of both the systems biology and the biotechnology of oilseeds.