Starch biosynthesis and intermediary metabolism in maize kernels. Quantitative analysis of metabolite flux by nuclear magnetic resonance

Genetic basis of kernel starch content decoded in a maize multi-parent population

Starch is the most abundant storage carbohydrate in maize kernels and provides calories for humans and other animals as well as raw materials for various industrial applications. Decoding the genetic basis of natural variation in kernel starch content is needed to manipulate starch quantity and quality via molecular breeding to meet future needs. Here, we identified 50 unique single quantitative trait loci (QTLs) for starch content with 18 novel QTLs via single linkage mapping, joint linkage mapping and a genome-wide association study in a multi-parent population containing six recombinant inbred line populations. Only five QTLs explained over 10% of phenotypic variation in single populations. In addition to a few large-effect and many small-effect additive QTLs, limited pairs of epistatic QTLs also contributed to the genetic basis of the variation in kernel starch content. A regional association study identified five non-starch-pathway genes that were the causal candidate genes underlying the identified QTLs for starch content. The pathway-driven analysis identified ZmTPS9, which encodes a trehalose-6-phosphate synthase in the trehalose pathway, as the causal gene for the QTL qSTA4-2, which was detected by all three statistical analyses. Knockout of ZmTPS9 increased kernel starch content and, in turn, kernel weight in maize, suggesting potential applications for ZmTPS9 in maize starch and yield improvement. These findings extend our knowledge about the genetic basis of starch content in maize kernels and provide valuable information for maize genetic improvement of starch quantity and quality.

Comparative phosphoproteomic analysis of developing maize seeds suggests a pivotal role for enolase in promoting starch synthesis

Maize (Zea mays) seeds are the major source of starch all over the world and the excellent model for researching starch synthesis. Seed starch content is a typical quantitative phenotype and many reports revealed that the glycolytic enzymes are involved in regulating starch synthesis, however the regulatory mechanism is still unclear. Here, we present a comparative phosphoproteomic study of three maize inbred lines with different seed starch content. It reveals that abundances of 62 proteins and 63 phosphoproteins were regulated during maize seed development. Dynamics of 17 enzymes related to glycolysis and starch synthesis were used to construct a phosphorylation regulatory network of starch synthesis. It shows that starch synthesis and glycolysis in maize seeds utilize the same hexose phosphates pool coming from sorbitol and sucrose as carbon source, and phosphorylation of ZmENO1 are suggested to contribute to increase starch content, because it is positively related to seed starch content in different developmental stages and different lines, and the phosphor-mimic mutant (ZmENO1^S43D) damaged its enzyme activity which is vital in glycolysis. Our results provide a new sight into regulatory process of seed starch synthesis and can be used in maize breeding for high starch content.

Genome-wide analysis of bacterial determinants of plant growth promotion and induced systemic resistance by Pseudomonas fluorescens

Pseudomonas fluorescens strain SS101 (Pf.SS101) promotes growth of Arabidopsis thaliana, enhances greening and lateral root formation, and induces systemic resistance (ISR) against the bacterial pathogen Pseudomonas syringae pv. tomato (Pst). Here, targeted and untargeted approaches were adopted to identify bacterial determinants and underlying mechanisms involved in plant growth promotion and ISR by Pf.SS101. Based on targeted analyses, no evidence was found for volatiles, lipopeptides and siderophores in plant growth promotion by Pf.SS101. Untargeted, genome-wide analyses of 7488 random transposon mutants of Pf.SS101 led to the identification of 21 mutants defective in both plant growth promotion and ISR. Many of these mutants, however, were auxotrophic and impaired in root colonization. Genetic analysis of three mutants followed by site-directed mutagenesis, genetic complementation and plant bioassays revealed the involvement of the phosphogluconate dehydratase gene edd, the response regulator gene colR and the adenylsulfate reductase gene cysH in both plant growth promotion and ISR. Subsequent comparative plant transcriptomics analyses strongly suggest that modulation of sulfur assimilation, auxin biosynthesis and transport, steroid biosynthesis and carbohydrate metabolism in Arabidopsis are key mechanisms linked to growth promotion and ISR by Pf.SS101.

Mapping and comparative proteomic analysis of the starch biosynthetic pathway in rice by 2D PAGE/MS

Our results not only provide a comprehensive overview of the starch biosynthetic pathway in the developing endosperm but also reveal some important protein markers that regulate the synthesis of starch. In human diets, rice (Oryza sativa L.) is an important source of starch, a substantial amount of which is accumulated in developing endosperm. A better understanding of the complicated pathways involved in starch biosynthesis is needed to improve the yield and quality of rice and other cereal crops through breeding. One pure line rice mutant, SA0419, was induced from a wild-type rice, TNG67, by sodium azide mutagenesis; therefore, TNG67 and SA0419 share the same genetic background. SA0419 is, however, a unique glutinous rice with a lower amylose content (8%) than that of TNG67 (20%), and the grains of SA0419 develop earlier and faster than those of TNG67. In this study, we used a comparative proteomic analysis to identify the differentially expressed proteins that may explain the differences in starch biosynthesis and the characteristics of TNG67 and SA0419. A gel-based proteomic approach was applied to profile the expressed proteome in the developing endosperm of these two rice varieties by nano-LC/MS/MS. Several over-expressed proteins were found in SA0419, such as plastidial ADP-glucose pyrophosphorylase (AGPase), phosphoglucomutase (PGM), pyrophosphate-fructose 6-phosphate 1-phosphotransferase (PFP), 6-phosphofructokinase (PFK), pyruvate phosphate dikinase (PPDK), starch branching enzymes (SBE) and starch debranching enzyme (SDBE), with those proteins mainly being involved in the pathways of starch metabolism and PPDK-mediated gluconeogenesis. Those over-expressed enzymes may contribute to the relatively early development, similar starch accumulation and rapid grain filling of SA0419 as compared with TNG67. This study provides a detailed biochemical description of starch biosynthesis and related information regarding a unique starch mutant that may assist future research efforts to improve the yield and quality of grain and starch in rice through breeding.

Nuclear Magnetic Resonance (NMR) Spectroscopy For Metabolic Profiling of Medicinal Plants and Their Products

NMR spectroscopy has multidisciplinary applications, including excellent impact in metabolomics. The analytical capacity of NMR spectroscopy provides information for easy qualitative and quantitative assessment of both endogenous and exogenous metabolites present in biological samples. The complexity of a particular metabolite and its contribution in a biological system are critically important for understanding the functional state that governs the organism's phenotypes. This review covers historical aspects of developments in the NMR field, its applications in chemical profiling, metabolomics, and quality control of plants and their derived medicines, foods, and other products. The bottlenecks of NMR in metabolic profiling are also discussed, keeping in view the future scope and further technological interventions.

The complex isotopologue space of glucose as a framework for the study of human intermediary metabolism

The positional distributions of stable isotopes in metabolites provide specific fingerprints of the pathways and fluxes that have occurred in the organisms under study. In particular, modern nuclear magnetic resonance (NMR) spectroscopy enables the detailed assignment of isotope patterns in natural products, for example, in metabolites obtained from labelling experiments using (13)C-enriched precursors, such as glucose, acetate or CO2. In this study, the transient (13)C-isotopologue composition of blood glucose from an adult human volunteer after intravenous supply of [U-(13)C6]glucose was determined by high-resolution (13)C NMR spectroscopy. The non-linear progression curves displaying the relative amounts of eight (13)C-glucose isotopologues reflected the contributions of glucose metabolism by glycolytic cycling, the pentose phosphate pathway and anaplerotic reactions involving the citric acid cycle. The pilot study suggests that the experimental setting can be useful in analysing under non-invasive conditions the impact of physiological and pharmacological constraints on glucose turnover in humans.

Metabolic cartography: experimental quantification of metabolic fluxes from isotopic labelling studies

For the past decade, flux maps have provided researchers with an in-depth perspective on plant metabolism. As a rapidly developing field, significant headway has been made recently in computation, experimentation, and overall understanding of metabolic flux analysis. These advances are particularly applicable to the study of plant metabolism. New dynamic computational methods such as non-stationary metabolic flux analysis are finding their place in the toolbox of metabolic engineering, allowing more organisms to be studied and decreasing the time necessary for experimentation, thereby opening new avenues by which to explore the vast diversity of plant metabolism. Also, improved methods of metabolite detection and measurement have been developed, enabling increasingly greater resolution of flux measurements and the analysis of a greater number of the multitude of plant metabolic pathways. Methods to deconvolute organelle-specific metabolism are employed with increasing effectiveness, elucidating the compartmental specificity inherent in plant metabolism. Advances in metabolite measurements have also enabled new types of experiments, such as the calculation of metabolic fluxes based on (13)CO(2) dynamic labelling data, and will continue to direct plant metabolic engineering. Newly calculated metabolic flux maps reveal surprising and useful information about plant metabolism, guiding future genetic engineering of crops to higher yields. Due to the significant level of complexity in plants, these methods in combination with other systems biology measurements are necessary to guide plant metabolic engineering in the future.

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

¹⁴C labeling experiments performed with kernel cultures showed that developing maize endosperm is more efficient than other non-photosynthetic tissues such as sunflower and maize embryos at converting maternally supplied substrates into biomass. To characterize the metabolic fluxes in endosperm, maize kernels were labeled to isotopic steady state using ¹³C-labeled glucose. The resultant labeling in free metabolites and biomass was analyzed by NMR and GC-MS. After taking into account the labeling of substrates supplied by the metabolically active cob, the fluxes through central metabolism were quantified by computer-aided modeling. The flux map indicates that 51-69% of the ATP produced is used for biomass synthesis and up to 47% is expended in substrate cycling. These findings point to potential engineering targets for improving yield and increasing oil contents by, respectively, reducing substrate cycling and increasing the commitment of plastidic carbon into fatty acid synthesis at the level of pyruvate kinase.

Understanding fatty acid synthesis in developing maize embryos using metabolic flux analysis

The efficiency with which developing maize embryos convert substrates into seed storage reserves was determined to be 57-71%, by incubating developing maize embryos with uniformly labeled 14C substrates and measuring their conversion to CO2 and biomass products. To map the pattern of metabolic fluxes underlying this efficiency, maize embryos were labeled to isotopic steady state using a combination of labeled 13C-substrates. Intermediary metabolic fluxes were estimated by computer-aided modeling of the central metabolic network using the labeling data collected by NMR and GC-MS and the biomass composition. The resultant flux map reveals that even though 36% of the entering carbon goes through the oxidative pentose-phosphate pathway, this does not fully meet the NADPH demands for fatty acid synthesis. Metabolic flux analysis and enzyme activities highlight the importance of plastidic NADP-dependent malic enzyme, which provides one-third of the carbon and NADPH required for fatty acid synthesis in developing maize embryos.

Metabolic flux analysis in plants: coping with complexity

Theory and experience in metabolic engineering both show that metabolism operates at the network level. In plants, this complexity is compounded by a high degree of compartmentation and the synthesis of a very wide array of secondary metabolic products. A further challenge to understanding and predicting plant metabolic function is posed by our ignorance about the structure of metabolic networks even in well-studied systems. Metabolic flux analysis (MFA) provides tools to measure and model the functioning of metabolism, and is making significant contributions to coping with their complexity. This review gives an overview of different MFA approaches, the measurements required to implement them and the information they yield. The application of MFA methods to plant systems is then illustrated by several examples from the recent literature. Next, the challenges that plant metabolism poses for MFA are discussed together with ways that these can be addressed. Lastly, new developments in MFA are described that can be expected to improve the range and reliability of plant MFA in the coming years.

Impact of the salt stress on the photosynthetic carbon flux and 13C-label distribution within floridoside and digeneaside in Solieria chordalis

The flux of photosynthetic carbon used in the synthesis of low-molecular weight carbohydrates (digeneaside and floridoside) was investigated by (13)C and (1)H NMR spectroscopy in samples of the red seaweed, Solieria chordalis, incubated at different salinities (22, 34 and 50psu). Carbohydrates were labelled, by pulse-chase, with the stable isotope (13)C from NaH(13)CO(3). In vivo NMR analyses carried out with a cryogenic probe optimised for (13)C detection were performed directly on the living algal tissues to evidence the labelling of the carbohydrates with neither preliminary extraction nor purification step(s). The isotopic enrichment of each compound was determined by high-resolution (1)H and (13)C NMR spectroscopy. These analyses evidenced different orientations of the flux of the photosynthetic carbon in the algae according to the salt stress. At normal and low salinities, the photosynthetic carbon flux was responsible of 70% and 67% of the floridoside synthetized during the pulse period, respectively, whereas it was only of 30% in the thalli exposed to the high salinity, meaning a biosynthesis of high floridoside amount from endogen source leading to the osmotic regulation. Under normal and hyper-osmotic conditions, the stock of floridoside was used for cellular needs during the chase period, whereas it was not under hypo-osmotic conditions. The characterization of isotopomers composition of floridoside and digeneaside and the analysis of adjacent (13)C-labelling gives much more details on the effects of salinity on the metabolic pathways leading to the synthesis or the degradation of these molecules. High turnover of floridoside was evidenced at normal salinity during the chase period and products issued from the degradation of floridoside would not be used for the novo biosynthesis. That suggests that synthesis and degradation of floridoside may be realized in two different cellular compartments. The presence of more numerous (13)C-(13)C blocks in the carbon skeleton of the molecules from the up salt stressed thalli than in those from no salt stressed algae, concomitant with a slower degree of isotopic enrichment of the molecule, provided evidence that the two metabolic pathways (endogen and photosynthetic) may not share the precursor molecules involved in the floridoside synthesis and that these two routes may be totally separate until the constitution of floridoside molecule.

The complexities of starch biosynthesis in cereal endosperms

Starch serves not only as an energy source for plants, animals, and humans but also as an environmentally friendly alternative for fossil fuels. Here, we describe recent findings concerning the synthesis of this important molecule in the cereal endosperm. Results from six separate transgenic reports point to the importance of adenosine diphosphate glucose pyrophosphorylase in controlling the amount of starch synthesized. The unexpected cause underlying the contrast in sequence divergence of its two subunits is also described. A major unresolved question concerning the synthesis of starch is the origin of nonrandom or clustered alpha-1,6 branch-points within the major component of starch, amylopectin. Developing evidence that several of the starch biosynthetic enzymes involved in amylopectin synthesis occur in complexes is reviewed. These complexes may provide the specificity for the formation of nonrandom branch-points.

Metabolic flux maps comparing the effect of temperature on protein and oil biosynthesis in developing soybean cotyledons

Metabolic flux maps developed from 13C metabolic flux analysis (13C MFA) are effective tools for assessing the response of biological systems to genetic or environmental perturbations, and for identifying possible metabolic engineering targets. Experimental treatments were designed to distinguish between temperature effects prior to, and during incubation in vitro, on primary metabolism in developing soybeans. Biomass accumulation increased with temperature as did carbon partitioning into lipids. The flux through the plastidic oxidative pentose phosphate pathway (pgl(P)) relative to sucrose intake remained fairly constant [ approximately 56% (+/-24%)] when cotyledons were transferred from an optimum growth temperature to varying temperatures in in vitro culture, signifying a rigid node under these conditions. However, pgl(P) flux ranged from 57 to 77% of sucrose intake when growth temperature in planta varied and were cultured in vitro at the same temperature (as the plant), indicating a flexible node for this case. The carbon flux through the anaplerotic reactions catalysed by plastidic malic enzyme (me(P)), cytosolic phosphoenolpyruvate (PEP) carboxylase and the malate (Mal) transporter from the cytosol to mitochondrion varied dramatically with temperature and had a direct influence on the carbon partitioning into protein and oil from the plastidic pyruvate (Pyr) pool. These results of the in vitro culture indicate that temperature during early stages of development has a dominant effect on establishing capacity for flux through certain components of central carbon metabolism.

13CO2 as a universal metabolic tracer in isotopologue perturbation experiments

A tobacco plant was illuminated for 5h in an atmosphere containing (13)CO(2) and then maintained for 10 days under standard greenhouse conditions. Nicotine, glucose, and amino acids from proteins were isolated chromatographically. Isotopologue abundances of isolated metabolites were determined quantitatively by NMR spectroscopy and mass spectrometry. The observed non-stochastic isotopologue patterns indicate (i) formation of multiply labeled photosynthetic carbohydrates during the (13)CO(2) pulse phase followed by (ii) partial catabolism of the primary photosynthetic products, and (iii) recombination of the (13)C-labeled fragments with unlabeled intermediary metabolites during the chase period. The detected and simulated isotopologue profiles of glucose and amino acids reflect carbon partitioning that is dominated by the Calvin cycle and glycolysis/glucogenesis. Retrobiosynthetic analysis of the nicotine pattern is in line with its known formation from nicotinic acid and putrescine via aspartate, glyceraldehyde phosphate and alpha-ketoglutarate as basic building blocks. The study demonstrates that pulse/chase labeling with (13)CO(2) as precursor is a powerful tool for the analysis of quantitative aspects of plant metabolism in completely unperturbed whole plants.

Experimental and mathematical approaches to modeling plant metabolic networks

To support their sessile and autotrophic lifestyle higher plants have evolved elaborate networks of metabolic pathways. Dynamic changes in these metabolic networks are among the developmental forces underlying the functional differentiation of organs, tissues and specialized cell types. They are also important in the various interactions of a plant with its environment. Further complexity is added by the extensive compartmentation of the various interconnected metabolic pathways in plants. Thus, although being used widely for assessing the control of metabolic flux in microbes, mathematical modeling approaches that require steady-state approximations are of limited utility for understanding complex plant metabolic networks. However, considerable progress has been made when manageable metabolic subsystems were studied. In this article, we will explain in general terms and using simple examples the concepts underlying stoichiometric modeling (metabolic flux analysis and metabolic pathway analysis) and kinetic approaches to modeling (including metabolic control analysis as a special case). Selected studies demonstrating the prospects of these approaches, or combinations of them, for understanding the control of flux through particular plant pathways are discussed. We argue that iterative cycles of (dry) mathematical modeling and (wet) laboratory testing will become increasingly important for simulating the distribution of flux in plant metabolic networks and deriving rational experimental designs for metabolic engineering efforts.

A transient isotopic labeling methodology for 13C metabolic flux analysis of photoautotrophic microorganisms

Metabolic flux analysis is increasingly recognized as an integral component of systems biology. However, techniques for experimental measurement of system-wide metabolic fluxes in purely photoautotrophic systems (growing on CO(2) as the sole carbon source) have not yet been developed due to the unique problems posed by such systems. In this paper, we demonstrate that an approach that balances positional isotopic distributions transiently is the only route to obtaining system-wide metabolic flux maps for purely autotrophic metabolism. The outlined transient (13)C-MFA methodology enables measurement of fluxes at a metabolic steady-state, while following changes in (13)C-labeling patterns of metabolic intermediates as a function of time, in response to a step-change in (13)C-label input. We use mathematical modeling of the transient isotopic labeling patterns of central intermediates to assess various experimental requirements for photoautotrophic MFA. This includes the need for intracellular metabolite concentration measurements and isotopic labeling measurements as a function of time. We also discuss photobioreactor design and operation in order to measure fluxes under precise environmental conditions. The transient MFA technique can be used to measure and compare fluxes under different conditions of light intensity, nitrogen sources or compare strains with various mutations or gene deletions and additions.

Recent applications of NMR spectroscopy in plant metabolomics

Recent research has established NMR as a key method for high-throughput comparative analysis of plant extracts. We discuss recent examples of the use of NMR to provide metabolomic data for various applications in plant science and look forward to the key role that NMR will play in data provision for plant systems biology.

Estimation of metabolic fluxes, expression levels and metabolite dynamics of a secondary metabolic pathway in potato using label pulse-feeding experiments combined with kinetic network modelling and simulation

In this paper we present a method that allows dynamic flux analysis without a priori kinetic knowledge. This method was developed and validated using the pulse-feeding experimental data obtained in our previous study (Matsuda et al., 2005), in which incorporation of exogenously applied l-phenylalanine-d(5) into seven phenylpropanoid metabolites in potato tubers was determined. After identification of the topology of the metabolic network of these biosynthetic pathways, the system was described by dynamic mass balances in combination with power-law kinetics. After the first simulations, some reactions were removed from the network because they were not contributing significantly to network behaviour. As a next step, the exponents of the power-law kinetics were identified and then kept at fixed values during further analysis. The model was tested for statistical reliability using Monte Carlo simulations. Most fluxes could be identified with high accuracy. The two test cases, control and after elicitation, were clearly distinguished, and with elicitation fluxes to N-p-coumaroyloctopamine (pCO) and N-p-coumaroyltyramine (pCT) increased significantly, whereas those for chlorogenic acid (CGA) and p-coumaroylshikimate decreased significantly. According to the model, increases in the first two fluxes were caused by induction/derepression mechanisms. The decreases in the latter two fluxes were caused by decreased concentrations of their substrates, which in turn were caused by increased activity of the pCO- and pCT-producing enzymes. Flux-control analysis showed that, in most cases, flux control was changed after application of elicitor. Thus the results revealed potential targets for improving actions against tissue wounding and pathogen attack.

Metabolic flux analysis: recent advances in carbon metabolism in plants

Isotopic tracers are used to both trace metabolic pathways and quantify fluxes through these pathways. The use of different labeling methods recently led to profound changes in our views of plant metabolism. Examples are taken from primary metabolism, with sugar interconversions, carbon partitioning between glycolysis and the pentose phosphate pathway, or metabolite inputs into the tricarboxylic acid (TCA) cycle, as well as from secondary metabolism with the relative contribution of the plastidial and cytosolic pathways to the biosynthesis of terpenoids. While labeling methods are often distinguished according to the instruments used for label detection, emphasis is put here on labeling duration. Short time labeling is adequate to study limited areas of the metabolic network. Long-term labeling, when designed to obtain metabolic and isotopic steady-state, allows to calculate various fluxes in large areas ofcentral metabolism. After longer labeling periods, large amounts of label accumulate in structural or storage compounds: their detailed study through the retrobiosynthetic method gives access to the biosynthetic pathways of otherwise undetectable precursors. This chapter presents the power and limits of the different methods, and illustrates how they can be associated with each other and with other methods of cell biology, to provide the information needed for a rational approach of metabolic engineering.

Robustness of central carbohydrate metabolism in developing maize kernels

The central carbohydrate metabolism provides the precursors for the syntheses of various storage products in seeds. While the underlying biochemical map is well established, little is known about the organization and flexibility of carbohydrate metabolic fluxes in the face of changing biosynthetic demands or other perturbations. This question was addressed in developing kernels of maize (Zea mays L.), a model system for the study of starch and sugar metabolism. (13)C-labeling experiments were carried out with inbred lines, heterotic hybrids, and starch-deficient mutants that were selected to cover a wide range of performances and kernel phenotypes. In total, 46 labeling experiments were carried out using either [U-(13)C(6)]glucose or [U-(13)C(12)]sucrose and up to three stages of kernel development. Carbohydrate flux distributions were estimated based on glucose isotopologue abundances, which were determined in hydrolysates of starch by using quantitative (13)C-NMR and GC-MS. Similar labeling patterns in all samples indicated robustness of carbohydrate fluxes in maize endosperm, and fluxes were rather stable in response to glucose or sucrose feeding and during development. A lack of ADP-glucose pyrophosphorylase in the bt2 and sh2 mutants triggered significantly increased hexose cycling. In contrast, other mutations with similar kernel phenotypes had no effect. Thus, the distribution of carbohydrate fluxes is stable and not determined by sink strength in maize kernels.

Measuring multiple fluxes through plant metabolic networks

Fluxes through metabolic networks are crucial for cell function, and a knowledge of these fluxes is essential for understanding and manipulating metabolic phenotypes. Labeling provides the key to flux measurement, and in network flux analysis the measurement of multiple fluxes allows a flux map to be superimposed on the metabolic network. The principles and practice of two complementary methods, dynamic and steady-state labeling, are described, emphasizing best practice and illustrating their contribution to network flux analysis with examples taken from the plant and microbial literature. The principal analytical methods for the detection of stable isotopes are also described, as well as the procedures for obtaining flux maps from labeling data. A series of boxes summarizing the key concepts of network flux analysis is provided for convenience.

Studies of the intermediary metabolism in cultured cells of the insect Spodoptera frugiperda using 13C- or 15N-labelled tracers

Background

Insect cells can serve as host systems for the recombinant expression of eukaryotic proteins. Using this platform, the controlled expression of ¹⁵N/¹³C labelled proteins requires the analysis of incorporation paths and rates of isotope-labelled precursors present in the medium into amino acids. For this purpose, Spodoptera frugiperda cells were grown in a complex medium containing [U-¹³C₆]glucose. In a second experiment, cultures of S. frugiperda were grown in the presence of ¹⁵N-phenylalanine.

Results

Quantitative NMR analysis showed incorporation of the proffered [U-¹³C₆]glucose into the ribose moiety of ribonucleosides (40 – 45%) and into the amino acids, alanine (41%), glutamic acid/glutamine (C-4 and C-5, 30%) and aspartate/asparagine (15%). Other amino acids and the purine ring of nucleosides were not formed from exogenous glucose in significant amounts (> 5%). Prior to the incorporation into protein the proffered ¹⁵N-phenylalanine lost about 70% of its label by transamination and the labelled compound was not converted into tyrosine to a significant extent.

Conclusion

Growth of S. frugiperda cells in the presence of [U-¹³C₆]glucose is conducive to the fractional labelling of ribonucleosides, alanine, glutamic acid/glutamine and aspartic acid/asparagine. The isotopolog compositions of the ribonucleosides and of alanine indicate considerable recycling of carbohydrate intermediates in the reductive branch of the pentose phosphate pathway. The incorporation of ¹⁵N-labelled amino acids may be hampered by loss of the ¹⁵N-label by transamination.

Changes in flux pattern of the central carbohydrate metabolism during kernel development in maize

Developing kernels of the inbred maize line W22 were grown in sterile culture and supplied with a mixture of [U-13C6]glucose and unlabeled glucose during three consecutive intervals (11-18, 18-25, or 25-32 days after pollination) within the linear phase of starch formation. At the end of each labeling period, glucose was prepared from starch and analyzed by 13C isotope ratio mass spectrometry and high-resolution (13)C NMR spectroscopy. The abundances of individual glucose isotopologs were calculated by computational deconvolution of the NMR data. [1,2-(13)C2]-, [5,6-(13)C2]-, [2,3-(13)C2]-, [4,5-(13)C2]-, [1,2,3-(13)C3]-, [4,5,6-(13)C3]-, [3,4,5,6-(13)C4]-, and [U-(13)C6]-isotopologs were detected as the major multiple-labeled glucose species, albeit at different normalized abundances in the three intervals. Relative flux contributions by five different pathways in the primary carbohydrate metabolism were determined by computational simulation of the isotopolog space of glucose. The relative fractions of some of these processes in the overall glucose cycling changed significantly during maize kernel development. The simulation showed that cycling via the non-oxidative pentose phosphate pathway was lowest during the middle interval of the experiment. The observed flux pattern could by explained by a low demand for amino acid precursors recruited from the pentose phosphate pathway during the middle interval of kernel development.

A new substrate cycle in plants. Evidence for a high glucose-phosphate-to-glucose turnover from in vivo steady-state and pulse-labeling experiments with [13C]glucose and [14C]glucose

Substrate (futile) cycling involving carbohydrate turnover has been widely reported in plant tissues, although its extent, mechanisms, and functions are not well known. In this study, two complementary approaches, short and steady-state labeling experiments, were used to analyze glucose metabolism in maize (Zea mays) root tips. Unidirectional rates of synthesis for storage compounds (starch, Suc, and cell wall polysaccharides) were determined by short labeling experiments using [U-14C]glucose and compared with net synthesis fluxes to determine the rate of glucose production from these storage compounds. Steady-state labeling with [1-(13)C]glucose and [U-13C]glucose showed that the redistribution of label between carbon C-1 and C-6 in glucose is close to that in cytosolic hexose-P. These results indicate a high resynthesis flux of glucose from hexose-P that is not accounted for by glucose recycling from storage compounds, thus suggesting the occurrence of a direct glucose-P-to-glucose conversion. An enzyme assay confirmed the presence of substantial glucose-6-phosphatase activity in maize root tips. This new glucose-P-to-glucose cycle was shown to consume around 40% of the ATP generated in the cell, whereas Suc cycling consumes at most 3% to 6% of the ATP produced. The rate of glucose-P cycling differs by a factor of 3 between a maize W22 line and the hybrid maize cv Dea, and is significantly decreased by a carbohydrate starvation pretreatment.

Revealing metabolic phenotypes in plants: inputs from NMR analysis

Assessing the performance of the plant metabolic network, with its varied biosynthetic capacity and its characteristic subcellular compartmentation, remains a considerable challenge. The complexity of the network is such that it is not yet possible to build large-scale predictive models of the fluxes it supports, whether on the basis of genomic and gene expression analysis or on the basis of more traditional measurements of metabolites and their interconversions. This limits the agronomic and biotechnological exploitation of plant metabolism, and it undermines the important objective of establishing a rational metabolic engineering strategy. Metabolic analysis is central to removing this obstacle and currently there is particular interest in harnessing high-throughput and/or large-scale analyses to the task of defining metabolic phenotypes. Nuclear magnetic resonance (NMR) spectroscopy contributes to this objective by providing a versatile suite of analytical techniques for the detection of metabolites and the fluxes between them. The principles that underpin the analysis of plant metabolism by NMR are described, including a discussion of the measurement options for the detection of metabolites in vivo and in vitro, and a description of the stable isotope labelling experiments that provide the basis for metabolic flux analysis. Despite a relatively low sensitivity, NMR is suitable for high-throughput system-wide analyses of the metabolome, providing methods for both metabolite fingerprinting and metabolite profiling, and in these areas NMR can contribute to the definition of plant metabolic phenotypes that are based on metabolic composition. NMR can also be used to investigate the operation of plant metabolic networks. Labelling experiments provide information on the operation of specific pathways within the network, and the quantitative analysis of steady-state labelling experiments leads to the definition of large-scale flux maps for heterotrophic carbon metabolism. These maps define multiple unidirectional fluxes between branch-points in the metabolic network, highlighting the existence of substrate cycles and discriminating in favourable cases between fluxes in the cytosol and plastid. Flux maps can be used to define a functionally relevant metabolic phenotype and the extensive application of such maps in microbial systems suggests that they could have important applications in characterising the genotypes produced by plant genetic engineering.

Quantification of compartmented metabolic fluxes in developing soybean embryos by employing biosynthetically directed fractional (13)C labeling, two-dimensional [(13)C, (1)H] nuclear magnetic resonance, and comprehensive isotopomer balancing

Metabolic flux quantification in plants is instrumental in the detailed understanding of metabolism but is difficult to perform on a systemic level. Toward this aim, we report the development and application of a computer-aided metabolic flux analysis tool that enables the concurrent evaluation of fluxes in several primary metabolic pathways. Labeling experiments were performed by feeding a mixture of U-(13)C Suc, naturally abundant Suc, and Gln to developing soybean (Glycine max) embryos. Two-dimensional [(13)C, (1)H] NMR spectra of seed storage protein and starch hydrolysates were acquired and yielded a labeling data set consisting of 155 (13)C isotopomer abundances. We developed a computer program to automatically calculate fluxes from this data. This program accepts a user-defined metabolic network model and incorporates recent mathematical advances toward accurate and efficient flux evaluation. Fluxes were calculated and statistical analysis was performed to obtain sds. A high flux was found through the oxidative pentose phosphate pathway (19.99 +/- 4.39 micromol d(-1) cotyledon(-1), or 104.2 carbon mol +/- 23.0 carbon mol per 100 carbon mol of Suc uptake). Separate transketolase and transaldolase fluxes could be distinguished in the plastid and the cytosol, and those in the plastid were found to be at least 6-fold higher. The backflux from triose to hexose phosphate was also found to be substantial in the plastid (21.72 +/- 5.00 micromol d(-1) cotyledon(-1), or 113.2 carbon mol +/-26.0 carbon mol per 100 carbon mol of Suc uptake). Forward and backward directions of anaplerotic fluxes could be distinguished. The glyoxylate shunt flux was found to be negligible. Such a generic flux analysis tool can serve as a quantitative tool for metabolic studies and phenotype comparisons and can be extended to other plant systems.

Understanding flux in plant metabolic networks

The revolutionary growth in our ability to identify the 'parts list' of cellular infrastructure in plants in detail, and to alter it with precision, challenges us to develop methods to quantify how these parts function. For components of metabolism, this means mapping fluxes at the level of metabolic networks. Advances in experimental, analytical and software tools for metabolic flux analysis now allow maps of the fluxes through central metabolism to be obtained from the results of stable-isotope-labeling experiments. Such maps have led to notable successes in understanding and engineering metabolic function in microorganisms. Recent studies in plants are giving insight into particular fluxes, such as those of the pentose phosphate pathway, and into general phenomena, such as substrate- or futile-cycles and compartmentation. The importance of experimental design and statistical analysis have been illustrated, and analyses of fluxes in heterotrophic plant tissues have been carried out recently.

Isotopolog perturbation techniques for metabolic networks: metabolic recycling of nutritional glucose in Drosophila melanogaster

Drosophila melanogaster strain Oregon-R(*) was grown on standard medium supplemented with [U-(13)C(6)]glucose. One to two days after hatching, flies were extracted with water. Glucose was isolated chromatographically from the extract and was analyzed by (13)C NMR spectroscopy. All (13)C signals of the isolated glucose were multiplets arising by (13)C(13)C coupling. Based on a comprehensive analysis of the coupling constants and heavy isotope shifts in glucose, the integrals of individual (13)C signal patterns afforded the concentrations of certain groups of (13)C isotopologs. These data were deconvoluted by a genetic algorithm affording the abundances of all single-labeled and of 15 multiply labeled isotopologs. Among the latter group, seven isotopologs were found at concentrations >0.1 mol % with [1,2-(13)C(2)]glucose as the most prominent species. The multiply (13)C-labeled glucose isotopologs are caused by metabolic remodeling of the proffered glucose via a complex network of catabolic and anabolic processes involving glycolysis and/or passage through the pentose phosphate, the Cori cycle and/or the citrate cycle. The perturbation method described can be adapted to a wide variety of experimental systems and isotope-labeled precursors.

Biosynthesis of camptothecin. In silico and in vivo tracer study from [1-13C]glucose

Camptothecin derivatives are clinically used antitumor alkaloids that belong to monoterpenoid indole alkaloids. In this study, we investigated the biosynthetic pathway of camptothecin from [1-13C]glucose (Glc) by in silico and in vivo studies. The in silico study measured the incorporation of Glc into alkaloids using the Atomic Reconstruction of Metabolism software and predicted the labeling patterns of successive metabolites from [1-13C]Glc. The in vivo study followed incorporation of [1-13C]Glc into camptothecin with hairy roots of Ophiorrhiza pumila by 13C nuclear magnetic resonance spectroscopy. The 13C-labeling pattern of camptothecin isolated from the hairy roots clearly showed that the monoterpene-secologanin moiety was synthesized via the 2C-methyl-D-erythritol 4-phosphate pathway, not via the mevalonate pathway. This conclusion was supported by differential inhibition of camptothecin accumulation by the pathway-specific inhibitors (fosmidomycin and lovastatin). The quinoline moiety from tryptophan was also labeled as predicted by the Atomic Reconstruction of Metabolism program via the shikimate pathway. These results indicate that camptothecin is formed by the combination of the 2C-methyl-D-erythritol 4-phosphate pathway and the shikimate pathway. This study provides the innovative example for how a computer-aided comprehensive metabolic analysis will refine the experimental design to obtain more precise biological information.

The oxidative pentose phosphate pathway: structure and organisation

The oxidative pentose phosphate pathway is a major source of reducing power and metabolic intermediates for biosynthetic processes. Some, if not all, of the enzymes of the pathway are found in both the cytosol and plastids, although the precise distribution of their activities varies. The apparent absence of sections of the pathway from the cytosol potentially complicates metabolism. These complications are partly offset, however, by exchange of intermediates between the cytosol and the plastids through the activities of a family of plastid phosphate translocators. Molecular analysis is confirming the widespread presence of multiple genes encoding each of the enzymes of the oxidative pentose phosphate pathway. Differential expression of these isozymes may ensure that the kinetic properties of the activity that catalyses a specific reaction match the metabolic requirements of a particular tissue. This hypothesis can be tested thanks to recent developments in the application of 13C-steady-state labelling strategies. These strategies make it possible to quantify flux through metabolic networks and to discriminate between pathways of carbohydrate oxidation in the cytosol and plastids.