Central metabolic fluxes in the endosperm of developing maize seeds and their implications for metabolic engineering

13C-labeling reveals non-conventional pathways providing carbon for hydroxy fatty acid synthesis in Physaria fendleri

Physaria fendleri is a member of the Brassicaceae that produces in its embryos hydroxy fatty acids, constituents of oils that are very valuable and widely used by industry for cosmetics, lubricants, biofuels, etc. Free of toxins and rich in hydroxy fatty acids, Physaria provides a promising alternative to imported castor oil and is on the verge of being commercialized. This study aims to identify important biochemical step(s) for oil synthesis in Physaria, which may serve as target(s) for future crop improvement. To advance towards this goal, the endosperm composition was analysed by LC-MS/MS to develop and validate culture conditions that mimic the development of the embryos in planta. Using developing Physaria embryos in culture and 13C-labeling, our studies revealed that: (i) Physaria embryos metabolize carbon into biomass with an efficiency significantly lower than other photosynthetic embryos; (ii) the plastidic malic enzyme provides 42% of the pyruvate used for de novo fatty acid synthesis, which is the highest measured so far in developing 'green' oilseed embryos; and (iii) Physaria uses non-conventional pathways to channel carbon into oil, namely the Rubisco shunt, which fixes CO2 released in the plastid, and the reversibility of isocitrate dehydrogenase, which provides additional carbon for fatty acid elongation.

Progress in understanding and improving oil content and quality in seeds

The world's population is projected to increase by two billion by 2050, resulting in food and energy insecurity. Oilseed crops have been identified as key to address these challenges: they produce and store lipids in the seeds as triacylglycerols that can serve as a source of food/feed, renewable fuels, and other industrially-relevant chemicals. Therefore, improving seed oil content and composition has generated immense interest. Research efforts aiming to unravel the regulatory pathways involved in fatty acid synthesis and to identify targets for metabolic engineering have made tremendous progress. This review provides a summary of the current knowledge of oil metabolism and discusses how photochemical activity and unconventional pathways can contribute to high carbon conversion efficiency in seeds. It also highlights the importance of ¹³C-metabolic flux analysis as a tool to gain insights on the pathways that regulate oil biosynthesis in seeds. Finally, a list of key genes and regulators that have been recently targeted to enhance seed oil production are reviewed and additional possible targets in the metabolic pathways are proposed to achieve desirable oil content and quality.

Effects of Fruit Shading on Gene and Protein Expression During Starch and Oil Accumulation in Developing Styrax tonkinensis Kernels

Styrax tonkinensis has great potential as a biofuel feedstock source having industrial oilseeds with excellent fatty acids (FAs) composition and good fuel properties. Photosynthesis in the developing pericarp could affect the carbon distribution in kernel. During kernel development, more carbon sources are allocated to starch rather than lipid, when the pericarp photosynthesis is reduced by fruit shading treatment. After shading the fruits at 50 days after flowering (DAF), samples of shaded fruit (FSK) and controls (CK) were collected at 80 DAF and analyzed using the proteomic method. We identified 3,181 proteins, of which 277 were differentially expressed proteins, all downregulated in the FSK group. There were 56 proteins found involved in carbohydrate metabolism and lipid biosynthesis leading to oil accumulation with their iTRAQ ratios of FSK/CK ranging from 0.7123 to 1.1075. According to the qRT-PCR analyses, the key genes related to FA and triacylglycerol (TAG) biosynthesis were significantly downregulated between 60 and 90 DAF especially at 80 DAF, while the key genes involved in starch biosynthesis and FA desaturase had no significant difference between the two groups at 80 DAF. Fruit shading is a negative treatment for lipid accumulation but not starch accumulation by restraining enzymic protein expression involved in FA and TAG biosynthesis during S. tonkinensis kernel development.

Establishment of a GC-MS-based 13 C-positional isotopomer approach suitable for investigating metabolic fluxes in plant primary metabolism

¹³ C-Metabolic flux analysis (¹³ C-MFA) has greatly contributed to our understanding of plant metabolic regulation. However, the generation of detailed in vivo flux maps remains a major challenge. Flux investigations based on nuclear magnetic resonance have resolved small networks with high accuracy. Mass spectrometry (MS) approaches have broader potential, but have hitherto been limited in their power to deduce flux information due to lack of atomic level position information. Herein we established a gas chromatography (GC) coupled to MS-based approach that provides ¹³ C-positional labelling information in glucose, malate and glutamate (Glu). A map of electron impact (EI)-mediated MS fragmentation was created and validated by ¹³ C-positionally labelled references via GC-EI-MS and GC-atmospheric pressure chemical ionization-MS technologies. The power of the approach was revealed by analysing previous ¹³ C-MFA data from leaves and guard cells, and ¹³ C-HCO₃ labelling of guard cells harvested in the dark and after the dark-to-light transition. We demonstrated that the approach is applicable to established GC-EI-MS-based ¹³ C-MFA without the need for experimental adjustment, but will benefit in the future from paired analyses by the two GC-MS platforms. We identified specific glucose carbon atoms that are preferentially labelled by photosynthesis and gluconeogenesis, and provide an approach to investigate the phosphoenolpyruvate carboxylase (PEPc)-derived ¹³ C-incorporation into malate and Glu. Our results suggest that gluconeogenesis and the PEPc-mediated CO₂ assimilation into malate are activated in a light-independent manner in guard cells. We further highlight that the fluxes from glycolysis and PEPc toward Glu are restricted by the mitochondrial thioredoxin system in illuminated leaves.

Comparative transcriptome and metabolome profiling in the maturing seeds of contrasting cluster bean (Cyamopsis tetragonoloba L. Taub) cultivars identified key molecular variations leading to increased gum accumulation

Cluster bean (Guar) is the major source of industrial gum. Knowledge on the molecular events regulating galactomannan gum accumulation in guar will pave way for accelerated development of gummy guar genotypes. RNA Seq analysis in the immature seeds of contrasting cluster bean genotypes HGS 563 (gum type) and Pusa Navbahar (vegetable type) resulted in the generation of 19,855,490 and 21,488,472 quality reads. Data analysis identified 4938 differentially expressed genes between the gummy vs vegetable genotypes. A set of 2241 genes were up-regulated and 2587 genes were down-regulated in gummy guar. Significant up-regulation of genes involved in the biosynthesis of galactomannan and cell wall storage polysaccharides was observed in the gummy HGS 563. Genes involved in carotenoids, flavonoids, non mevalonic acid, terpenoids, and wax metabolism were also up-regulated in HGS 563. Mannose and galactose were the major nucleotide sugars in Pusa Navbahar and HGS 563 immature seeds. Immature seeds of HGS 563 showed high concentration of mannose and galactose accumulation compared to Pusa Navbahar. qRT-PCR analysis of selected genes confirmed the findings of transcriptome data.

Sugar modulation of anaerobic-response networks in maize root tips

Sugar supply is a key component of hypoxia tolerance and acclimation in plants. However, a striking gap remains in our understanding of mechanisms governing sugar impacts on low-oxygen responses. Here, we used a maize (Zea mays) root-tip system for precise control of sugar and oxygen levels. We compared responses to oxygen (21 and 0.2%) in the presence of abundant versus limited glucose supplies (2.0 and 0.2%). Low-oxygen reconfigured the transcriptome with glucose deprivation enhancing the speed and magnitude of gene induction for core anaerobic proteins (ANPs). Sugar supply also altered profiles of hypoxia-responsive genes carrying G4 motifs (sources of regulatory quadruplex structures), revealing a fast, sugar-independent class followed more slowly by feast-or-famine-regulated G4 genes. Metabolite analysis showed that endogenous sugar levels were maintained by exogenous glucose under aerobic conditions and demonstrated a prominent capacity for sucrose re-synthesis that was undetectable under hypoxia. Glucose abundance had distinctive impacts on co-expression networks associated with ANPs, altering network partners and aiding persistence of interacting networks under prolonged hypoxia. Among the ANP networks, two highly interconnected clusters of genes formed around Pyruvate decarboxylase 3 and Glyceraldehyde-3-phosphate dehydrogenase 4. Genes in these clusters shared a small set of cis-regulatory elements, two of which typified glucose induction. Collective results demonstrate specific, previously unrecognized roles of sugars in low-oxygen responses, extending from accelerated onset of initial adaptive phases by starvation stress to maintenance and modulation of co-expression relationships by carbohydrate availability.

Model-assisted identification of metabolic engineering strategies for Jatropha curcas lipid pathways

Efficient approaches to increase plant lipid production are necessary to meet current industrial demands for this important resource. While Jatropha curcas cell culture can be used for in vitro lipid production, scaling up the system for industrial applications requires an understanding of how growth conditions affect lipid metabolism and yield. Here we present a bottom-up metabolic reconstruction of J. curcas supported with labeling experiments and biomass characterization under three growth conditions. We show that the metabolic model can accurately predict growth and distribution of fluxes in cell cultures and use these findings to pinpoint energy expenditures that affect lipid biosynthesis and metabolism. In addition, by using constraint-based modeling approaches we identify network reactions whose joint manipulation optimizes lipid production. The proposed model and computational analyses provide a stepping stone for future rational optimization of other agronomically relevant traits in J. curcas.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Non-conventional pathways enable pennycress (Thlaspi arvense L.) embryos to achieve high efficiency of oil biosynthesis

Pennycress (Thlaspi arvense L.) accumulates oil up to 35% of the total seed biomass, and its overall fatty acid composition is suitable for aviation fuel. However, for this plant to become economically viable, its oil production needs to be improved. In vivo culture conditions that resemble the development of pennycress embryos in planta were developed based on the composition of the liquid endosperm. Then, substrate uptake rates and biomass accumulation were measured from cultured pennycress embryos, revealing a biosynthetic efficiency of 93%, which is one of the highest in comparison with other oilseeds to date. Additionally, the ratio of carbon in oil to CO2 indicated that non-conventional pathways are likely to be responsible for such a high carbon conversion efficiency. To identify the reactions enabling this phenomenon, parallel labeling experiments with 13C-labeled substrates were conducted in pennycress embryos. The main findings of these labeling experiments include: (i) the occurrence of the oxidative reactions of the pentose phosphate pathway in the cytosol; (ii) the reversibility of isocitrate dehydrogenase; (iii) the operation of the plastidic NADP-dependent malic enzyme; and (iv) the refixation of CO2 by Rubisco. These reactions are key providers of carbon and reductant for fatty acid synthesis and elongation.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

Liquid Chromatography Tandem Mass Spectrometry Quantification of 13C-Labeling in Sugars

Subcellular compartmentation has been challenging in plant ¹³C-metabolic flux analysis. Indeed, plant cells are highly compartmented: they contain vacuoles and plastids in addition to the regular organelles found in other eukaryotes. The distinction of reactions between compartments is possible when metabolites are synthesized in a particular compartment or by a unique pathway. Sucrose is an example of such a metabolite: it is specifically produced in the cytosol from glucose 6-phosphate (G6P) and fructose 6-phosphate (F6P). Therefore, determining the ¹³C-labeling in the fructosyl and glucosyl moieties of sucrose directly informs about the labeling of cytosolic F6P and G6P, respectively. To date, the most commonly used method to monitor sucrose labeling is by nuclear magnetic resonance, which requires substantial amounts of biological sample. This study describes a new methodology that accurately measures the labeling in free sugars using liquid chromatography tandem mass spectrometry (LC-MS/MS). For this purpose, maize embryos were pulsed with [U-¹³C]-fructose, intracellular sugars were extracted, and their time-course labeling was analyzed by LC-MS/MS. Additionally, extracts were enzymatically treated with hexokinase to remove the soluble hexoses, and then invertase to cleave sucrose into fructose and glucose. Finally, the labeling in the glucosyl and fructosyl moieties of sucrose was determined by LC-MS/MS.

Metabolite Profiling and Classification of Developing Styrax tonkinensis Kernels

Background: Styrax tonkinensis is an economic tree species with high timber, medicine, oil, and ornamental value. Its seed, containing a particularly high oil content, are widely studied for their biodiesel properties by nutritional components and oil body ultrastructure. However, their comprehensive biochemical compositions have not been studied. Methods: During S. tonkinensis kernel development, we collected samples from four time points for metabolite profiling and classification through gas chromatography-mass spectrometry and liquid chromatography-mass spectrometry. Results: A total of 187 and 1556 metabolites were obtained, respectively. All of the metabolites were grouped into 19 and 21 classes by their chemical properties and into 8 clusters based on their change trends, respectively. Among all the metabolites, carboxylic acids and derivatives, flavonoids, fatty acyls, glycerophospholipids, organooxygen compounds, prenol lipids, and steroids and steroid derivatives were the main components. Alanine, glutamine, tryptophan, tyrosine and valine were the five most abundant amino acids. Palmitic acid, stearic acid, oleic acid and linoleic acid were the four major free fatty acids. Flavans, flavonoid glycosides and o-methylated flavonoids were the three major flavonoids. The differential metabolites distributions between different time points were identified. A pathway enrichment was performed, which was mainly focused on three groups, amino acids metabolism, carbon flow from sucrose to lipid and secondary metabolites biosynthesis. Conclusions: It’s the first time to analyze the metabolite fingerprinting for developing S. tonkinensis kernels and identify varied kinds of flavonoids. We performed metabolite profiling, classification and pathway enrichment to assess the comprehensive biochemical compositions. Our results described the change in major metabolites and main metabolic processes during S. tonkinensis kernel development and provided a variety of bases for seed applications as biofuel or medicine.

Analysis of the Barley Malt Rootlet Proteome

Barley seeds are one of the main ingredients of the malting industry for brewing beer. The barley rootlets that are separated from the kilned seeds at the end of the malting process and used as animal feed are one of the byproducts of this industry. In this study, the proteome of rootlets derived from two stages of the malting process, germination and kilning, from a popular malting barley variety were analyzed. A label-free shotgun proteomics strategy was used to identify more than 800 proteins from the barley rootlets. A high coverage and high confidence Gene Ontology annotations of the barley genome was used to facilitate the functional annotation of the proteins that were identified in the rootlets. An analysis of these proteins using Kellogg Encyclopedia of Genes and Genomes (KEGG) and Plant Reactome databases indicated the enrichment of pathways associated with phytohormones, protein biosynthesis, secondary metabolism, and antioxidants. Increased levels of jasmonic acid and auxin in the rootlets further supported the in silico analysis. As a rich source of proteins and amino acids use of these by-products of the malting industry for animal feed is validated. This study also indicates rootlets as a potential source of naturally occurring phenylpropanoids and antioxidants that can be further exploited in the development of functional foods.

Shoot tip culture: a step towards 13C metabolite flux analysis of sink leaf metabolism

BACKGROUND

Better understanding of the physiological and metabolic status of plants can only be obtained when metabolic fluxes are accurately assessed in a growing plant. Steady state ¹³C-MFA has been established as a routine method for analysis of fluxes in plant primary metabolism. However, the experimental system needs to be improved for continuous carbon enrichment from labelled sugars into metabolites for longer periods until complex secondary metabolism reaches steady state.

RESULTS

We developed an in vitro plant culture strategy by using peppermint as a model plant with minimizing unlabelled carbon fixation where growing shoot tip was strongly dependent on labelled glucose for their carbon necessity. We optimized the light condition and detected the satisfactory phenotypical growth under the lower light intensity. Total volatile terpenes were also highest at the same light. Analysis of label incorporation into pulegone monoterpene after continuous U-¹³C₆ glucose feeding revealed nearly 100% ¹³C, even at 15 days after first leaf visibility (DALV). Label enrichment was gradually scrambled with increasing light intensity and leaf age. This study was validated by showing similar levels of label enrichment in proteinogenic amino acids. The efficiency of this method was also verified in oregano.

CONCLUSIONS

Our shoot tip culture depicted a method in achieving long term, stable and a high percentage of label accumulation in secondary metabolites within a fully functional growing plant system. It recommends the potential application for the investigations of various facets of plant metabolism by steady state ¹³C-MFA. The system also provides a greater potential to study sink leaf metabolism. Overall, we propose a system to accurately describe complex metabolic phenotypes in a growing plant.

Sustained substrate cycles between hexose phosphates and free sugars in phosphate-deficient potato (Solanum tuberosum) cell cultures

Futile cycling between free sugars and hexose phosphates occurring under phosphate deficiency could be involved in the maintenance of a threshold level of free cellular phosphate to preserve respiratory metabolism. We studied the metabolic response of potato cell cultures growing in Pi sufficient (2.5 mM, +Pi) or deficient (125 µM, -Pi) conditions. Under Pi deficiency, cellular growth was severely affected, however -Pi cells were able to maintain a low but steady level of free Pi. We surveyed the activities of 33 primary metabolic enzymes during the course of a 12 days Pi deficiency period. Our results show that many of these enzymes had higher specific activity in -Pi cells. Among these, we found typical markers of Pi deficiency such as phosphoenolpyruvate phosphatase and phosphoenolpyruvate carboxylase as well as enzymes involved in the biosynthesis of organic acids. Intriguingly, several ATP-consuming enzymes such as hexokinase (HK) and phosphofructokinase also displayed increased activity in -Pi condition. For HK, this was associated with an increase in the steady state of a specific HK polypeptide. Quantification of glycolytic intermediates showed a pronounced decrease in phosphate esters under Pi deficiency. Adenylate levels also decreased in -Pi cells, but the Adenylate Energy Charge was not affected by the treatment. To investigate the significance of HK induction under low Pi, [U-¹⁴C]-glucose tracer studies were conducted. We found in vivo evidence of futile cycling between pools of hexose phosphates and free sugars under Pi deficiency. Our study suggests that the futile cycling between hexose phosphates and free sugars which is active under +Pi conditions is sustained under Pi deficiency. The possibility that this process represents a metabolic adaptation to Pi deficiency is discussed with respect to Pi homeostasis in Pi-deficient conditions.

Gluconeogenesis and nitrogen metabolism in maize

Two pathways can be used by gluconeogenesis in plants: one employs phosphoenolpyruvate carboxykinase (PEPCK) and the other pyruvate orthophosphate dikinase (PPDK). The occurrence-location of these enzymes was determined in developing kernels of maize. PPDK was much more abundant than PEPCK in extracts of whole kernels. However, their location within the kernel was different. PPDK was particularly abundant in the peripheral endosperm (in which alanine is abundant), whereas PEPCK was localised in the pedicel and basal endosperm transfer cells (where asparagine is metabolised). The abundance of these enzymes was also determined in maize roots where there was a massive increase in abundance of PEPCK and a small increase in abundance of PPDK when they were fed ammonium; PEPCK was located in the pericycle and various cell types associated with the vasculature. On the other hand, there was a large increase in abundance of PPDK in roots subjected to anoxia (which induces an accumulation of alanine), whereas the abundance of PEPCK was decreased. These results show: firstly, that gluconeogenesis can potentially occur in many different tissues of maize. Secondly, within one organ PPDK can be abundant in some tissues and PEPCK in others. Thirdly, the abundance of PPDK and PEPCK is often associated with the metabolism of certain nitrogenous compounds and can be dramatically altered by factors related to nitrogen metabolism. In maize roots and developing kernels PPDK was associated with alanine metabolism. By contrast, the presence of PEPCK in maize roots and kernels was associated with either ammonium or asparagine metabolism. We propose that gluconeogenesis is often a component of a widespread mechanism that is used in coordinating the import/mobilisation of nitrogenous compounds with their utilisation. Further, potentially component of this mechanism may have provided building blocks that were used in the evolution of processes such as C₄ photosynthesis, Crassulacean acid metabolism, stomatal metabolism and the biochemical pH stat.

High-Resolution Magic-Angle-Spinning NMR and Magnetic Resonance Imaging Spectroscopies Distinguish Metabolome and Structural Properties of Maize Seeds from Plants Treated with Different Fertilizers and Arbuscular mycorrhizal fungi

Both high-resolution magic-angle-spinning (HRMAS) and magnetic resonance imaging (MRI) NMR spectroscopies were applied here to identify the changes of metabolome, morphology, and structural properties induced in seeds (caryopses) of maize plants grown at field level under either mineral or compost fertilization in combination with the inoculation by arbuscular mycorrhizal fungi (AMF). The metabolome of intact caryopses was examined by HRMAS-NMR, while the morphological aspects, endosperm properties and seed water distribution were investigated by MRI. Principal component analysis (PCA) was applied to evaluate ¹H CPMG (Carr-Purcel-Meiboom-Gill) HRMAS spectra as well as several MRI-derived parameters ( T₁, T₂, and self-diffusion coefficients) of intact maize caryopses. PCA score-plots from spectral results indicated that both seeds metabolome and structural properties depended on the specific field treatment undergone by maize plants. Our findings show that a combination of multivariate statistical analyses with advanced and nondestructive NMR techniques, such as HRMAS and MRI, enables the evaluation of the effects induced on maize caryopses by different fertilization and management practices at field level. The spectroscopic approach adopted here may become useful for the objective appraisal of the quality of seeds produced under a sustainable agriculture.

Gas Chromatography-Mass Spectrometry-Based 13C-Labeling Studies in Plant Metabolomics

Stable-isotope labeling analysis has been used to discover new metabolic pathways and their key regulatory points in a wide range of organisms. Given the complexity of the plant metabolic network, this analysis provides information complementary to that obtained from metabolite profiling that can be used to understand how plants cope with adverse conditions, and how metabolism varies between different cells, tissues, and organs. Here we describe the experimental procedures from sample harvesting and extraction to mass spectral analysis and interpretation that allow the researcher to perform ¹³C-labeling experiments. A wide range of plant material, from single cells to whole plants, can be used to investigate the metabolic fate of the ¹³C from a predefined tracer. Thus, a key point of this analysis is to choose the correct biological system, the substrate and the condition to be investigated; all of which implicitly relies on the biological question to be investigated. Rapid sample quenching and a careful data analysis are also critical points in such studies. By contrast to other metabolomic approaches, stable-isotope labeling can provide information concerning the fluxes through metabolic networks, which is essential for understanding and manipulating metabolic phenotypes and therefore of pivotal importance for both systems biology and plant metabolic engineering.

Starch Biosynthesis in the Developing Endosperms of Grasses and Cereals

The starch-rich endosperms of the Poaceae, which includes wild grasses and their domesticated descendents the cereals, have provided humankind and their livestock with the bulk of their daily calories since the dawn of civilization up to the present day. There are currently unprecedented pressures on global food supplies, largely resulting from population growth, loss of agricultural land that is linked to increased urbanization, and climate change. Since cereal yields essentially underpin world food and feed supply, it is critical that we understand the biological factors contributing to crop yields. In particular, it is important to understand the biochemical pathway that is involved in starch biosynthesis, since this pathway is the major yield determinant in the seeds of six out of the top seven crops grown worldwide. This review outlines the critical stages of growth and development of the endosperm tissue in the Poaceae, including discussion of carbon provision to the growing sink tissue. The main body of the review presents a current view of our understanding of storage starch biosynthesis, which occurs inside the amyloplasts of developing endosperms.

Cell-Type Specific Metabolic Flux Analysis: A Challenge for Metabolic Phenotyping and a Potential Solution in Plants

Stable isotope labelling experiments are used routinely in metabolic flux analysis (MFA) to determine the metabolic phenotype of cells and tissues. A complication arises in multicellular systems because single cell measurements of transcriptomes, proteomes and metabolomes in multicellular organisms suggest that the metabolic phenotype will differ between cell types. In silico analysis of simulated metabolite isotopomer datasets shows that cellular heterogeneity confounds conventional MFA because labelling data averaged over multiple cell types does not necessarily yield averaged flux values. A potential solution to this problem—the use of cell-type specific reporter proteins as a source of cell-type specific labelling data—is proposed and the practicality of implementing this strategy in the roots of Arabidopsis thaliana seedlings is explored. A protocol for the immunopurification of ectopically expressed green fluorescent protein (GFP) from Arabidopsis thaliana seedlings using a GFP-binding nanobody is developed, and through GC-MS analysis of protein hydrolysates it is established that constitutively expressed GFP reports accurately on the labelling of total protein in root tissues. It is also demonstrated that the constitutive expression of GFP does not perturb metabolism. The principal obstacle to the implementation of the method in tissues with cell-type specific GFP expression is the sensitivity of the GC-MS system.

Plant metabolism as studied by NMR spectroscopy

The study of plant metabolism impacts a broad range of domains such as plant cultural practices, plant breeding, human or animal nutrition, phytochemistry and green biotechnologies. Plant metabolites are extremely diverse in terms of structure or compound families as well as concentrations. This review attempts to illustrate how NMR spectroscopy, with its broad variety of experimental approaches, has contributed widely to the study of plant primary or specialized metabolism in very diverse ways. The review presents recent developments of one-dimensional and multi-dimensional NMR methods to study various aspects of plant metabolism. Through recent examples, it highlights how NMR has proved to be an invaluable tool for the global characterization of sample composition within metabolomic studies, and shows some examples of use for targeted phytochemistry, with a special focus on compound identification and quantitation. In such cases, NMR approaches are often used to provide snapshots of the plant sample composition. The review also covers dynamic aspects of metabolism, with a description of NMR techniques to measure metabolic fluxes - in most cases after stable isotope labelling. It is mainly intended for NMR specialists who would be interested to learn more about the potential of their favourite technique in plant sciences and about specific details of NMR approaches in this field. Therefore, as a practical guide, a paragraph on the specific precautions that should be taken for sample preparation is also included. In addition, since the quality of NMR metabolic studies is highly dependent on approaches to data processing and data sharing, a specific part is dedicated to these aspects. The review concludes with perspectives on the emerging methods that could change significantly the role of NMR in the field of plant metabolism by boosting its sensitivity. The review is illustrated throughout with examples of studies selected to represent diverse applications of liquid-state or HR-MAS NMR.

Metabolic flux analysis of heterotrophic growth in Chlamydomonas reinhardtii

Despite the wealth of knowledge available for C. reinhardtii, the central metabolic fluxes of growth on acetate have not yet been determined. In this study, ¹³C-metabolic flux analysis (¹³C-MFA) was used to determine and quantify the metabolic pathways of primary metabolism in C. reinhardtii cells grown under heterotrophic conditions with acetate as the sole carbon source. Isotopic labeling patterns of compartment specific biomass derived metabolites were used to calculate the fluxes. It was found that acetate is ligated with coenzyme A in the three subcellular compartments (cytosol, mitochondria and plastid) included in the model. Two citrate synthases were found to potentially be involved in acetyl-coA metabolism; one localized in the mitochondria and the other acting outside the mitochondria. Labeling patterns demonstrate that Acetyl-coA synthesized in the plastid is directly incorporated in synthesis of fatty acids. Despite having a complete TCA cycle in the mitochondria, it was also found that a majority of the malate flux is shuttled to the cytosol and plastid where it is converted to oxaloacetate providing reducing equivalents to these compartments. When compared to predictions by flux balance analysis, fluxes measured with ¹³C-MFA were found to be suboptimal with respect to biomass yield; C. reinhardtii sacrifices biomass yield to produce ATP and reducing equivalents.

Fluxomics links cellular functional analyses to whole-plant phenotyping

Fluxes through metabolic pathways reflect the integration of genetic and metabolic regulations. While it is attractive to measure all the mRNAs (transcriptome), all the proteins (proteome), and a large number of the metabolites (metabolome) in a given cellular system, linking and integrating this information remains difficult. Measurement of metabolome-wide fluxes (termed the fluxome) provides an integrated functional output of the cell machinery and a better tool to link functional analyses to plant phenotyping. This review presents and discusses sets of methodologies that have been developed to measure the fluxome. First, the principles of metabolic flux analysis (MFA), its 'short time interval' version Inst-MFA, and of constraints-based methods, such as flux balance analysis and kinetic analysis, are briefly described. The use of these powerful methods for flux characterization at the cellular scale up to the organ (fruits, seeds) and whole-plant level is illustrated. The added value given by fluxomics methods for unravelling how the abiotic environment affects flux, the process, and key metabolic steps are also described. Challenges associated with the development of fluxomics and its integration with 'omics' for thorough plant and organ functional phenotyping are discussed. Taken together, these will ultimately provide crucial clues for identifying appropriate target plant phenotypes for breeding.

Plant lipidomics at the crossroads: From technology to biology driven science

The identification and quantification of lipids from plant tissues have become commonplace and many researchers now incorporate lipidomics approaches into their experimental studies. Plant lipidomics research continues to involve technological developments such as those in mass spectrometry imaging, but in large part, lipidomics approaches have matured to the point of being accessible to the novice. Here we review some important considerations for those planning to apply plant lipidomics to their biological questions, and offer suggestions for appropriate tools and practices. This article is part of a Special Issue entitled: BBALIP_Lipidomics Opinion Articles edited by Sepp Kohlwein.

High-throughput quantification of the levels and labeling abundance of free amino acids by liquid chromatography tandem mass spectrometry

Accurate assessment of mass isotopomer distributions (MIDs) of intracellular metabolites, such as free amino acids (AAs), is crucial for quantifying in vivo fluxes. To date, the majority of studies that measured AA MIDs have relied on the analysis of proteinogenic rather than free AAs by: i) GC-MS, which involved cumbersome process of derivatization, or ii) NMR, which requires large quantities of biological sample. In this work, the development and validation of a high-throughput LC-MS/MS method allowing the quantification of the levels and labeling of free AAs is described. Sensitivity in the order of the femtomol was achieved using multiple reaction monitoring mode (MRM). The MIDs of all free AAs were assessed without the need of derivatization, and were validated (except for Trp) on a mixture of unlabeled AA standards. Finally, this method was applied to the determination of the ¹³C-labeling abundance in free AAs extracted from maize embryos cultured with ¹³C-glutamine or ¹³C-glucose. Although Cys was below the limit of detection in these biological samples, the MIDs of a total of 18 free AAs were successfully determined. Due to the increased application of tandem mass spectrometry for ¹³C-Metabolic Flux Analysis, this novel method will enable the assessment of more complete and accurate labeling information of intracellular AAs, and therefore a better definition of the fluxes.

T-6b allocates more assimilation product for oil synthesis and less for polysaccharide synthesis during the seed development of Arabidopsis thaliana

BACKGROUND

As an Agrobacterium tumefaciens T-DNA oncogene, T-6b induces the development of tumors and the enation syndrome in vegetative tissues of transgenic plants. Most of these effects are related to increases in soluble sugar contents. To verify the potential roles of T-6b in the distribution of carbon in developing seeds, not in vegetative tissues, we fused an endosperm-specific promoter to the T-6b gene for expression in transgenic Arabidopsis thaliana plants.

RESULTS

The expression of T-6b in reproductive organs did not induce the development of the enation syndrome, and moreover, promoted endosperm expansion, which increased the total seed biomass by more than 10%. Additionally, T-6b also increased oil content in mature seeds by more than 10% accompanied with the decrease of starch and mucilage content at the same time.

CONCLUSIONS

T-6b enhances seed biomass and helps oil biosynthesis but not polysaccharides in reproductive organs without disturbing vegetative growth and development. Our findings suggest T-6b may be very useful for increasing oil production in biodiesel plants.

Understanding the control of acyl flux through the lipid metabolic network of plant oil biosynthesis

Plant oil biosynthesis involves a complex metabolic network with multiple subcellular compartments, parallel pathways, cycles, and pathways that have a dual function to produce essential membrane lipids and triacylglycerol. Modern molecular biology techniques provide tools to alter plant oil compositions through bioengineering, however with few exceptions the final composition of triacylglycerol cannot be predicted. One reason for limited success in oilseed bioengineering is the inadequate understanding of how to control the flux of fatty acids through various fatty acid modification, and triacylglycerol assembly pathways of the lipid metabolic network. This review focuses on the mechanisms of acyl flux through the lipid metabolic network, and highlights where uncertainty resides in our understanding of seed oil biosynthesis. This article is part of a Special Issue entitled: Plant Lipid Biology edited by Kent D. Chapman and Ivo Feussner.

Assessing compartmentalized flux in lipid metabolism with isotopes

Metabolism in plants takes place across multiple cell types and within distinct organelles. The distributions equate to spatial heterogeneity; though the limited means to experimentally assess metabolism frequently involve homogenizing tissues and mixing metabolites from different locations. Most current isotope investigations of metabolism therefore lack the ability to resolve spatially distinct events. Recognition of this limitation has resulted in inspired efforts to advance metabolic flux analysis and isotopic labeling techniques. Though a number of these efforts have been applied to studies in central metabolism; recent advances in instrumentation and techniques present an untapped opportunity to make similar progress in lipid metabolism where the use of stable isotopes has been more limited. These efforts will benefit from sophisticated radiolabeling reports that continue to enrich our knowledge on lipid biosynthetic pathways and provide some direction for stable isotope experimental design and extension of MFA. Evidence for this assertion is presented through the review of several elegant stable isotope studies and by taking stock of what has been learned from radioisotope investigations when spatial aspects of metabolism were considered. The studies emphasize that glycerolipid production occurs across several locations with assembly of lipids in the ER or plastid, fatty acid biosynthesis occurring in the plastid, and the generation of acetyl-CoA and glycerol-3-phosphate taking place at multiple sites. Considering metabolism in this context underscores the cellular and subcellular organization that is important to enhanced production of glycerolipids in plants. An attempt is made to unify salient features from a number of reports into a diagrammatic model of lipid metabolism and propose where stable isotope labeling experiments and further flux analysis may help address questions in the field. This article is part of a Special Issue entitled: Plant Lipid Biology edited by Kent D. Chapman and Ivo Feussner.

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.

Combined transcript, proteome, and metabolite analysis of transgenic maize seeds engineered for enhanced carotenoid synthesis reveals pleotropic effects in core metabolism

The aim of this study was to assess whether endosperm-specific carotenoid biosynthesis influenced core metabolic processes in maize embryo and endosperm and how global seed metabolism adapted to this expanded biosynthetic capacity. Although enhancement of carotenoid biosynthesis was targeted to the endosperm of maize kernels, a concurrent up-regulation of sterol and fatty acid biosynthesis in the embryo was measured. Targeted terpenoid analysis, and non-targeted metabolomic, proteomic, and transcriptomic profiling revealed changes especially in carbohydrate metabolism in the transgenic line. In-depth analysis of the data, including changes of metabolite pools and increased enzyme and transcript concentrations, gave a first insight into the metabolic variation precipitated by the higher up-stream metabolite demand by the extended biosynthesis capacities for terpenoids and fatty acids. An integrative model is put forward to explain the metabolic regulation for the increased provision of terpenoid and fatty acid precursors, particularly glyceraldehyde 3-phosphate and pyruvate or acetyl-CoA from imported fructose and glucose. The model was supported by higher activities of fructokinase, glucose 6-phosphate isomerase, and fructose 1,6-bisphosphate aldolase indicating a higher flux through the glycolytic pathway. Although pyruvate and acetyl-CoA utilization was higher in the engineered line, pyruvate kinase activity was lower. A sufficient provision of both metabolites may be supported by a by-pass in a reaction sequence involving phosphoenolpyruvate carboxylase, malate dehydrogenase, and malic enzyme.

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.

Tracking the metabolic pulse of plant lipid production with isotopic labeling and flux analyses: Past, present and future

Metabolism is comprised of networks of chemical transformations, organized into integrated biochemical pathways that are the basis of cellular operation, and function to sustain life. Metabolism, and thus life, is not static. The rate of metabolites transitioning through biochemical pathways (i.e., flux) determines cellular phenotypes, and is constantly changing in response to genetic or environmental perturbations. Each change evokes a response in metabolic pathway flow, and the quantification of fluxes under varied conditions helps to elucidate major and minor routes, and regulatory aspects of metabolism. To measure fluxes requires experimental methods that assess the movements and transformations of metabolites without creating artifacts. Isotopic labeling fills this role and is a long-standing experimental approach to identify pathways and quantify their metabolic relevance in different tissues or under different conditions. The application of labeling techniques to plant science is however far from reaching it potential. In light of advances in genetics and molecular biology that provide a means to alter metabolism, and given recent improvements in instrumentation, computational tools and available isotopes, the use of isotopic labeling to probe metabolism is becoming more and more powerful. We review the principal analytical methods for isotopic labeling with a focus on seminal studies of pathways and fluxes in lipid metabolism and carbon partitioning through central metabolism. Central carbon metabolic steps are directly linked to lipid production by serving to generate the precursors for fatty acid biosynthesis and lipid assembly. Additionally some of the ideas for labeling techniques that may be most applicable for lipid metabolism in the future were originally developed to investigate other aspects of central metabolism. We conclude by describing recent advances that will play an important future role in quantifying flux and metabolic operation in plant tissues.

Improved evidence-based genome-scale metabolic models for maize leaf, embryo, and endosperm

There is a growing demand for genome-scale metabolic reconstructions for plants, fueled by the need to understand the metabolic basis of crop yield and by progress in genome and transcriptome sequencing. Methods are also required to enable the interpretation of plant transcriptome data to study how cellular metabolic activity varies under different growth conditions or even within different organs, tissues, and developmental stages. Such methods depend extensively on the accuracy with which genes have been mapped to the biochemical reactions in the plant metabolic pathways. Errors in these mappings lead to metabolic reconstructions with an inflated number of reactions and possible generation of unreliable metabolic phenotype predictions. Here we introduce a new evidence-based genome-scale metabolic reconstruction of maize, with significant improvements in the quality of the gene-reaction associations included within our model. We also present a new approach for applying our model to predict active metabolic genes based on transcriptome data. This method includes a minimal set of reactions associated with low expression genes to enable activity of a maximum number of reactions associated with high expression genes. We apply this method to construct an organ-specific model for the maize leaf, and tissue specific models for maize embryo and endosperm cells. We validate our models using fluxomics data for the endosperm and embryo, demonstrating an improved capacity of our models to fit the available fluxomics data. All models are publicly available via the DOE Systems Biology Knowledgebase and PlantSEED, and our new method is generally applicable for analysis transcript profiles from any plant, paving the way for further in silico studies with a wide variety of plant genomes.

Isotopically nonstationary 13C flux analysis of changes in Arabidopsis thaliana leaf metabolism due to high light acclimation

Improving plant productivity is an important aim for metabolic engineering. There are few comprehensive methods that quantitatively describe leaf metabolism, although such information would be valuable for increasing photosynthetic capacity, enhancing biomass production, and rerouting carbon flux toward desirable end products. Isotopically nonstationary metabolic flux analysis (INST-MFA) has been previously applied to map carbon fluxes in photoautotrophic bacteria, which involves model-based regression of transient (13)C-labeling patterns of intracellular metabolites. However, experimental and computational difficulties have hindered its application to terrestrial plant systems. We performed in vivo isotopic labeling of Arabidopsis thaliana rosettes with (13)CO2 and estimated fluxes throughout leaf photosynthetic metabolism by INST-MFA. Plants grown at 200 µmol m(-2)s(-1) light were compared with plants acclimated for 9 d at an irradiance of 500 µmol⋅m(-2)⋅s(-1). Approximately 1,400 independent mass isotopomer measurements obtained from analysis of 37 metabolite fragment ions were regressed to estimate 136 total fluxes (54 free fluxes) under each condition. The results provide a comprehensive description of changes in carbon partitioning and overall photosynthetic flux after long-term developmental acclimation of leaves to high light. Despite a doubling in the carboxylation rate, the photorespiratory flux increased from 17 to 28% of net CO2 assimilation with high-light acclimation (Vc/Vo: 3.5:1 vs. 2.3:1, respectively). This study highlights the potential of (13)C INST-MFA to describe emergent flux phenotypes that respond to environmental conditions or plant physiology and cannot be obtained by other complementary approaches.

Preferred carbon precursors for lipid labelling in the heterotrophic endosperm of developing oat (Avena sativa L.) grains

Oat (Avena sativa L.) is unusual among the cereal grains in storing high amounts of oil in the endosperm; up to 90% of total grain oil. By using oat as a model species for oil metabolism in the cereal endosperm, we can learn how to develop strategies to redirect carbon from starch to achieve high-oil yielding cereal crops. Carbon precursors for lipid synthesis were compared in two genetically close oat cultivars with different endosperm oil content (about 6% and 10% of grain dw, medium-oil; MO, and high-oil; HO cultivar, respectively) by supplying a variety of (14)C-labelled substrates to the grain from both up- and downstream parts of glycolysis, either through detached oat panicles in vitro or by direct injection in planta. When supplied by direct injection, (14)C from acetate was identified to label the lipid fraction of the grain to the highest extent among substrates tested; 46% of net accumulated (14)C, demonstrating its applicability as a marker for lipids in the endosperm. Time course analyses of injected (14)C acetate during grain development suggested a more efficient transfer of fatty acids from polar lipids to triacylglycerol in the HO as compared to the MO cultivar, and turnover of triacylglycerol was suggested to not play a major role for the final oil content of oat grain endosperm despite the low amount of protective oleosins in this tissue. Moreover, availability of light was shown to drastically affect grain net carbon accumulation from (14)C-sucrose when supplied through detached panicles for the HO cultivar.

¹⁴C pulse labeling to estimate external fluxes and turnovers in primary metabolism

Steady state (13)C-MFA is classically used to measure fluxes in complex metabolic networks. However, the modeling of steady state labeling allows the quantification of internal fluxes only and requires the estimation, by other methods, of the external fluxes, corresponding to substrate uptake (carbon input into the network) and to the production rate of compounds that accumulate within plant cells (network output). Additionally, it is not always possible to discriminate between different pathways that lead to the same label distribution. Methods to measure fluxes, based on direct measurements of pool size and on (14)C short-time labeling experiments, are described in this chapter. To illustrate this approach, we focus on the quantification of sucrose and starch turnovers.

Genetic and molecular bases of yield-associated traits: a translational biology approach between rice and wheat

Transferring the knowledge bases between related species may assist in enlarging the yield potential of crop plants. Being cereals, rice and wheat share a high level of gene conservation; however, they differ at metabolic levels as a part of the environmental adaptation resulting in different yield capacities. This review focuses on the current understanding of genetic and molecular regulation of yield-associated traits in both crop species, highlights the similarities and differences and presents the putative knowledge gaps. We focus on the traits associated with phenology, photosynthesis, and assimilate partitioning and lodging resistance; the most important drivers of yield potential. Currently, there are large knowledge gaps in the genetic and molecular control of such major biological processes that can be filled in a translational biology approach in transferring genomics and genetics informations between rice and wheat.

Brassica napus seed endosperm - metabolism and signaling in a dead end tissue

Oilseeds are an important element of human nutrition and of increasing significance for the production of industrial materials. The development of the seeds is based on a coordinated interplay of the embryo and its surrounding tissue, the endosperm. This study aims to give insights into the physiological role of endosperm for seed development in the oilseed crop Brassica napus. Using protein separation by two-dimensional (2D) isoelectric focusing (IEF)/SDS polyacrylamide gel electrophoresis (PAGE) and protein identification by mass spectrometry three proteome projects were carried out: (i) establishment of an endosperm proteome reference map, (ii) proteomic characterization of endosperm development and (iii) comparison of endosperm and embryo proteomes. The endosperm proteome reference map comprises 930 distinct proteins, including enzymes involved in genetic information processing, carbohydrate metabolism, environmental information processing, energy metabolism, cellular processes and amino acid metabolism. To investigate dynamic changes in protein abundance during seed development, total soluble proteins were extracted from embryo and endosperm fractions at defined time points. Proteins involved in sugar converting and recycling processes, ascorbate metabolism, amino acid biosynthesis and redox balancing were found to be of special importance for seed development in B. napus. Implications for the seed filling process and the function of the endosperm for seed development are discussed.

BIOLOGICAL SIGNIFICANCE

The endosperm is of key importance for embryo development during seed formation in plants. We present a broad study for characterizing endosperm proteins in the oilseed plant B. napus. Furthermore, a project on the biochemical interplay between the embryo and the endosperm during seed development is presented. We provide evidence that the endosperm includes a complete set of enzymes necessary for plant primary metabolism. Combination of our results with metabolome data will further improve systems-level understanding of the seed filling process and provide rational strategies for plant bioengineering.

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

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

Elucidation and in planta reconstitution of the parthenolide biosynthetic pathway

Parthenolide, the main bioactive compound of the medicinal plant feverfew (Tanacetum parthenium), is a promising anti-cancer drug. However, the biosynthetic pathway of parthenolide has not been elucidated yet. Here we report on the isolation and characterization of all the genes from feverfew that are required for the biosynthesis of parthenolide, using a combination of 454 sequencing of a feverfew glandular trichome cDNA library, co-expression analysis and metabolomics. When parthenolide biosynthesis was reconstituted by transient co-expression of all pathway genes in Nicotiana benthamiana, up to 1.4μgg(-1) parthenolide was produced, mostly present as cysteine and glutathione conjugates. These relatively polar conjugates were highly active against colon cancer cells, with only slightly lower activity than free parthenolide. In addition to these biosynthetic genes, another gene encoding a costunolide and parthenolide 3β-hydroxylase was identified opening up further options to improve the water solubility of parthenolide and therefore its potential as a drug.

Targeted metabolomics of Physaria fendleri, an industrial crop producing hydroxy fatty acids

Physaria fendleri (syn. Lesquerella) is a Brassicaceae producing lesquerolic acid, a highly valued hydroxy fatty acid that could be used for several industrial applications, such as cosmetics, lubricating greases, paints, plastics and biofuels. Free of toxins, Physaria oil is an attractive alternative to imported castor (Ricinus communis) oil, and is hence on the verge of commercialization. Gas chromatography-mass spectrometry analysis of fatty acid methyl esters revealed that lesquerolic acid was synthesized and accumulated in the embryos, reaching 60% (w/w) of the total fatty acids. The sequential extraction and characterization of biomass compounds revealed that Physaria embryo metabolism switched from protein to fatty acid biosynthesis between 18 and 24 days post-anthesis (DPA). In order to unravel the metabolic pathways involved in fatty acid synthesis, a targeted metabolomics study was conducted on Physaria embryos at different stages of development. For this purpose, two novel high-throughput liquid chromatography-tandem mass spectrometry methods were developed and validated to quantify sugars, sugar alcohols and amino acids. Specificity was achieved using multiple reaction monitoring, and the limits of quantification were in the pmole-fmole range. The comparative metabolomic study underlined that: (i) the majority of the metabolites accumulate in Physaria embryos between 18 and 27 DPA; (ii) the oxidative pentose phosphate pathway, glycolysis, the tricarboxilic acid cycle and the anaplerotic pathway drain a substantial amount of carbon; and (iii) ribulose-1,5-bisphosphate is present, which specifically indicates that the Calvin cycle is occurring. The importance and the relevance of these findings regarding fatty acid synthesis were discussed.

Metabolic flux analysis using ¹³C peptide label measurements

¹³C metabolic flux analysis (MFA) has become the experimental method of choice to investigate the cellular metabolism of microbes, cell cultures and plant seeds. Conventional steady-state MFA utilizes isotopic labeling measurements of amino acids obtained from protein hydrolysates. To retain spatial information in conventional steady-state MFA, tissues or subcellular fractions must be dissected or biochemically purified. In contrast, peptides retain their identity in complex protein extracts, and may therefore be associated with a specific time of expression, tissue type and subcellular compartment. To enable 'single-sample' spatially and temporally resolved steady-state flux analysis, we investigated the suitability of peptide mass distributions (PMDs) as an alternative to amino acid label measurements. PMDs are the discrete convolution of the mass distributions of the constituent amino acids of a peptide. We investigated the requirements for the unique deconvolution of PMDs into amino acid mass distributions (AAMDs), the influence of peptide sequence length on parameter sensitivity, and how AAMD and flux estimates that are determined through deconvolution compare to estimates from a conventional GC-MS measurement-based approach. Deconvolution of PMDs of the storage protein β-conglycinin of soybean (Glycine max) resulted in good AAMD and flux estimates if fluxes were directly fitted to PMDs. Unconstrained deconvolution resulted in inferior AAMD and flux estimates. PMD measurements do not include amino acid backbone fragments, which increase the information content in GC-MS-derived analyses. Nonetheless, the resulting flux maps were of comparable quality due to the precision of Orbitrap quantification and the larger number of peptide measurements.

High-throughput data pipelines for metabolic flux analysis in plants

In this chapter we illustrate the methodology for high-throughput metabolic flux analysis. Central to this is developing an end to end data pipeline, crucial for integrating the wet lab experiments and analytics, combining hardware and software automation, and standardizing data representation providing importers and exporters to support third party tools. The use of existing software at the start, data extraction from the chromatogram, and the end, MFA analysis, allows for the most flexibility in this workflow. Developing iMS2Flux provided a standard, extensible, platform independent tool to act as the "glue" between these end points. Most importantly this tool can be easily adapted to support different data formats, data verification and data correction steps allowing it to be central to managing the data necessary for high-throughput MFA. An additional tool was needed to automate the MFA software and in particular to take advantage of the course grained parallel nature of high-throughput analysis and available high performance computing facilities.In combination these methods show the development of high-throughput pipelines that allow metabolic flux analysis to join as a full member of the omics family.

Liquid chromatography tandem mass spectrometry for measuring ¹³C-labeling in intermediates of the glycolysis and pentose phosphate pathway

This chapter describes a procedure to analyze (13)C-labeled phosphorylated compounds by liquid chromatography tandem mass spectrometry. Phosphorylated compounds, intermediaries of the glycolysis and pentose phosphate pathway, are separated by anion exchange chromatography and their isotopic labeling is determined by mass spectrometry. A sensitivity in the fmole range is achieved using scheduled multiple reaction monitoring mode.

Isotopomer measurement techniques in metabolic flux analysis I: nuclear magnetic resonance

Two-dimensional [(1)H, (13)C] heteronuclear single quantum correlation (HSQC) spectroscopy nuclear magnetic resonance (NMR) is a comprehensive tool in metabolic flux analysis using (13)C-labeling experiments. NMR is particularly relevant when extensive isotopomer measurements are required, such as for plant cells and tissues, which contain multiple cellular compartments. Several isotope isomers (isotopomers) can be detected and their distribution extracted quantitatively from a single 2-D HSQC NMR spectrum. For example, 2-D HSQC detects the labeling patterns of adjacent carbon atoms and provides the enrichment of individual carbon atoms of the amino acids and glucosyl and mannosyl units present in hydrolysates of glycosylated protein. The HSQC analysis can quantitatively distinguish differences between the glucosyl units in the starch hydrolysate and a protein hydrolysate of plant biomass: this specifies crucial information about compartmentalization in the plant system. The peak structures obtained from the HSQC experiment show multiplet patterns that are directly related to the isotopomer abundances. These abundances have a nonlinear relationship to the fluxes via isotopomer balancing. Fluxes are obtained from the numerical solution of these balances and a stoichiometric model that includes biomass composition data as well as consumption rates of carbohydrate and nitrogen sources. Herein, we describe the methods for the experimental measurements for flux analysis, i.e., determination of the biomass composition (lipid, protein, soluble sugar, and starch) as well as detailed procedures of acid hydrolysis of protein and starch samples and NMR sample preparation, using soybean embryo culture as the model plant system. Techniques to obtain the relative intensity of 16 amino acids and glucosyl units for protein hydrolysate and the glucosyl units of starch hydrolysate of soybean embryos in 2-D HSQC NMR spectra also are provided.

Analysis of proteinogenic amino acid and starch labeling by 2D NMR

Comprehensive analysis of isotopic labeling patterns of metabolites in proteinogenic amino acids and starch for plant systems lay in the powerful tool of 2-Dimensional [(1)H, (13)C] Nuclear Magnetic Resonance (2D NMR) spectroscopy. From (13)C-labeling experiments, 2D NMR provides information on the labeling of particular carbon positions, which contributes to the quantification of positional isotope isomers (isotopomer). 2D Heteronuclear Single Quantum Correlation (HSQC) NMR distinguishes particularly between the labeling patterns of adjacent carbon atoms, and leads to a characteristic enrichment of each carbon atom of amino acids and glucosyl and mannosyl units present in hydrolysates of glycosylated protein. Furthermore, this technique can quantitatively classify differences in glucosyl units of starch hydrolysate and of protein hydrolysate of plant biomass. Therefore, the 2D HSQC NMR method uses proteinogenic amino acids and starch to provide an understanding of carbon distribution of compartmentalization in the plant system. NMR has the advantage of minimal sample handle without separate individual compounds prior to analysis, for example multiple isotopomers can be detected, and their distribution extracted quantitatively from a single 2D HSQC NMR spectrum. The peak structure obtained from the HSQC experiment show multiplet patterns, which are directly related to isotopomer balancing. These abundances can be translated to maximum information on the metabolic flux analysis. Detailed methods for the extractions of protein, oil, soluble sugars, and starch, hydrolysis of proteinogenic amino acid and starch, and NMR preparation using soybean embryos cultured in vitro as a model plant systems are reported in this text. In addition, this chapter includes procedures to obtain the relative intensity of 16 amino acids and glucosyl units from protein hydrolysate and the glucosyl units of starch hydrolysate of soybean embryos in 2D HSQC NMR spectra.