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

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

The metabolic response of cultured tomato cells to low oxygen stress

The storage of fruits and vegetables under a controlled atmosphere can induce low oxygen stress, which can lead to post‐harvest losses through the induction of disorders such as core breakdown and browning. To gain better understanding of the metabolic response of plant organs to low oxygen, cultured tomato cells (Lycopersicum esculentum) were used as a model system to study the metabolic stress response to low oxygen (0 and 1 kPa O₂). By adding ¹³C labelled glucose, changes in the levels of polar metabolites and their ¹³C label accumulation were quantified. Low oxygen stress altered the metabolite profile of tomato cells, with the accumulation of the intermediates of glycolysis in addition to increases in lactate and sugar alcohols. ¹³C label data showed reduced label accumulation in almost all metabolites except lactate and some sugar alcohols. The results showed that low oxygen stress in tomato cell culture activated fermentative metabolism and sugar alcohol synthesis while inhibiting the activity of the TCA cycle and the biosynthesis of metabolites whose precursors are derived from central metabolism, including fluxes to most organic acids, amino acids and sugars.     

Metabolic profiling of transgenic tomato plants overexpressing hexokinase reveals that the influence of hexose phosphorylation diminishes during fruit development

We have conducted a comprehensive metabolic profiling on tomato (Lycopersicon esculentum) leaf and developing fruit tissue using a recently established gas chromatography-mass spectrometry profiling protocol alongside conventional spectrophotometric and liquid chromatographic methodologies. Applying a combination of these techniques, we were able to identify in excess of 70 small-M(r) metabolites and to catalogue the metabolite composition of developing tomato fruit. In addition to comparing differences in metabolite content between source and sink tissues of the tomato plant and after the change in metabolite pool sizes through fruit development, we have assessed the influence of hexose phosphorylation through fruit development by analyzing transgenic plants constitutively overexpressing Arabidopsis hexokinase AtHXK1. Analysis of the total hexokinase activity in developing fruits revealed that both wild-type and transgenic fruits exhibit decreasing hexokinase activity with development but that the relative activity of the transgenic lines with respect to wild type increases with development. Conversely, both point-by-point and principal component analyses suggest that the metabolic phenotype of these lines becomes less distinct from wild type during development. In summary, the data presented in this paper demonstrate that the influence of hexose phosphorylation diminishes during fruit development and highlights the importance of greater temporal resolution of metabolism.     

Differences in regulation of carbohydrate metabolism during early fruit development between domesticated tomato and two wild relatives

Early development and growth of fruit in the domesticated tomato Solanum lycopersicum cultivar Money Maker and two of its wild relatives, S. peruvianum LA0385 and S. habrochaites LA1777, were studied. Although small differences exist, the processes involved and the sequence of events in fruit development are similar in all three species. The growth of developing fruits is exponential and the relative growth rate accelerates from 5 days after pollination (DAP 5) to DAP 8, followed by a decline during further development. Growth is positively correlated to the standard “Brix plus starch’’ in the period DAP 8–DAP 20. Carbohydrate composition and levels of sugars and organic acids differ in fruits of the wild accessions compared to domesticated tomato. The wild accessions accumulate sucrose instead of glucose and fructose, and ripe fruits contain higher levels of malate and citrate. The enzymes responsible for the accumulation of glucose and fructose in domesticated tomatoes are soluble invertase and sucrose synthase. The regulation of initial carbohydrate metabolism in the domesticated tomato differs from that in the wild species, as could be concluded from measuring activities of enzymes involved in primary carbohydrate metabolism. Furthermore, changes in the activity of several enzymes, e.g., cell wall invertase, soluble invertase, fructokinase and phosphoglucomutase, could be attributed to changes in gene expression level. For other enzymes, additional control mechanisms play a role in the developing tomato fruits. Localization by in-situ activity staining of enzymes showed comparable results for fruits of domesticated tomato and the wild accessions. However, in the pericarp of S. peruvianum, less activity staining of phosphogluco-isomerase, phosphoglucomutase and UDP-glucosepyrophosphorylase was observed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Changes in Structural Features of Free N-Glycan and Endoglycosidase Activity during Tomato Fruit Ripening

In developing plants, free N-glycans occur ubiquitously at micromolar concentrations. Such oligosaccharides have been proposed to be signaling molecules in plant development. As a part of a study to elucidate the physiological roles of de-N-glycosylation machinery involved in fruit ripening, we analyzed changes in the amounts and structural features of free N-glycans in tomato fruits at four ripening stages. The amount of high-mannose type free N-glycans increased significantly in accordance with fruit ripening, and the relative amounts of high-molecular size N-glycans, such as Man_8-9GlcNAc₁, became predominant. These observations suggest that the de-N-glycosylation machinery, including endo-β-N-acetylglucosaminidase (ENGase) activity, is stimulated in the later stages of fruit ripening. But contrary to expectation, we found that total ENGase activities in the tomato fruits did not vary significantly with the ripening process, suggesting that ENGase activity must be maintained at a certain level, and that the expression of α-mannosidase involved in the clearance of free N-glycans decreases during tomato fruit ripening.     

Effect of amino acid supplementation on titer and glycosylation distribution in hybridoma cell cultures—Systems biology‐based interpretation using genome‐scale metabolic flux balance model and multivariate data analysis

Genome‐scale flux balance analysis (FBA) is a powerful systems biology tool to characterize intracellular reaction fluxes during cell cultures. FBA estimates intracellular reaction rates by optimizing an objective function, subject to the constraints of a metabolic model and media uptake/excretion rates. A dynamic extension to FBA, dynamic flux balance analysis (DFBA), can calculate intracellular reaction fluxes as they change during cell cultures. In a previous study by Read et al. (2013), a series of informed amino acid supplementation experiments were performed on twelve parallel murine hybridoma cell cultures, and this data was leveraged for further analysis (Read et al., Biotechnol Prog. 2013;29:745–753). In order to understand the effects of media changes on the model murine hybridoma cell line, a systems biology approach is applied in the current study. Dynamic flux balance analysis was performed using a genome‐scale mouse metabolic model, and multivariate data analysis was used for interpretation. The calculated reaction fluxes were examined using partial least squares and partial least squares discriminant analysis. The results indicate media supplementation increases product yield because it raises nutrient levels extending the growth phase, and the increased cell density allows for greater culture performance. At the same time, the directed supplementation does not change the overall metabolism of the cells. This supports the conclusion that product quality, as measured by glycoform assays, remains unchanged because the metabolism remains in a similar state. Additionally, the DFBA shows that metabolic state varies more at the beginning of the culture but less by the middle of the growth phase, possibly due to stress on the cells during inoculation. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1163–1173, 2016     

Model-assisted analysis of tomato fruit growth in relation to carbon and water fluxes

This work proposed a model of tomato growth adapted from the Fishman and Génard model developed to predict carbon and water accumulation in peach fruit. The main adaptations relied on the literature on tomato and mainly concerned: (i) the decrease in cell wall extensibility coefficient during fruit development; (ii) the increase in the membrane reflection coefficient to solute from 0 to 1, which accounted for the switch from symplasmic to apoplasmic phloem unloading, and (iii) the negative influence of the initial fruit weight on the maximum rate of active carbon uptake based on the assumption of higher competition for carbon among cells in large fruits containing more cells. A sensitivity analysis was performed and the model was calibrated and evaluated with satisfaction on 17 experimental datasets obtained under contrasting environmental (temperature, air vapour pressure deficit) and plant (plant fruit load and fruit position) conditions. Then the model was used to analyse the variations in the main fluxes involved in tomato fruit growth and accumulation of carbon in response to virtual carbon and water stresses. The conclusions are that this model, integrating simple biophysical laws, was able to simulate the complex fruit behaviour in response to external or internal factors and thus it may be a powerful tool for managing fruit growth and quality.     

Insights into plant metabolic networks from steady-state metabolic flux analysis

Steady-state metabolic flux analysis (MFA) is an experimental approach that allows the measurement of multiple fluxes in the core network of primary carbon metabolism. It is based on isotopic labelling experiments, and although well established in the analysis of micro-organisms, and some mammalian systems, the extension of the method to plant cells has been challenging because of the extensive subcellular compartmentation of the metabolic network. Despite this difficulty there has been substantial progress in developing robust protocols for the analysis of heterotrophic plant metabolism by steady-state MFA, and flux maps have now been published that reflect the metabolic phenotypes of excised root tips, developing embryos and cotyledons, hairy root cultures, and cell suspensions under a variety of physiological conditions. There has been a steady improvement in the quality, extent and statistical reliability of these analyses, and new information is emerging on the performance of the plant metabolic network and the contributions of specific pathways.     

Effect of Gibberellin and Auxin on Parthenocarpic Fruit Growth Induction in the cv Micro-Tom of Tomato

The effect of applied gibberellin (GA) and auxin on fruit-set and growth has been investigated in tomato (Solanum lycopersicum L.) cv Micro-Tom. It was found that to prevent competition between developing fruits only one fruit per truss should be left on the plant. Unpollinated ovaries responded to GA₃ and to different auxins [indol-3-acetic acid, naphthaleneacetic acid, and 2,4-dichlorophenoxyacetic acid (2,4-D)], 2,4-D being the most efficient. GA₃- and 2,4-D-induced fruits had different internal morphology, with poor locular tissue development in the case of GA, and pseudoembryos development in the case of 2,4-D. Also, GA₃ produced larger cells in the internal region of the mesocarp (IM) associated with higher mean C values, whereas 2,4-D produced more cell layers in the pericarp than pollinated fruits. The smaller size of GA₃- compared with 2,4-D-induced fruits was due to them having fewer cells, only partially compensated by the larger size of IM cells. Simultaneous application of GA₃ and 2,4-D produced parthenocarpic fruits similar to pollinated fruits, but for the absence of seeds, suggesting that both kinds of hormones are involved in the induction of fruit development upon pollination. It is concluded that Micro-Tom constitutes a convenient model system, compared to tall cultivars, to investigate the hormonal regulation of fruit development in tomato.     

A hybrid kinetic and constraint-based model of leaf metabolism allows predictions of metabolic fluxes in different environments

While flux balance analysis (FBA) provides a framework for predicting steady-state leaf metabolic network fluxes, it does not readily capture the response to environmental variables without being coupled to other modelling formulations. To address this, we coupled an FBA model of 903 reactions of soybean (Glycine max) leaf metabolism with e-photosynthesis, a dynamic model that captures the kinetics of 126 reactions of photosynthesis and associated chloroplast carbon metabolism. Successful coupling was achieved in an iterative formulation in which fluxes from e-photosynthesis were used to constrain the FBA model and then, in turn, fluxes computed from the FBA model used to update parameters in e-photosynthesis. This process was repeated until common fluxes in the two models converged. Coupling did not hamper the ability of the kinetic module to accurately predict the carbon assimilation rate, photosystem II electron flux, and starch accumulation of field-grown soybean at two CO₂ concentrations. The coupled model also allowed accurate predictions of additional parameters such as nocturnal respiration, as well as analysis of the effect of light intensity and elevated CO₂ on leaf metabolism. Predictions included an unexpected decrease in the rate of export of sucrose from the leaf at high light, due to altered starch–sucrose partitioning, and altered daytime flux modes in the tricarboxylic acid cycle at elevated CO₂. Mitochondrial fluxes were notably different between growing and mature leaves, with greater anaplerotic, tricarboxylic acid cycle and mitochondrial ATP synthase fluxes predicted in the former, primarily to provide carbon skeletons and energy for protein synthesis.     

Quantitative approaches for analysing fluxes through plant metabolic networks using NMR and stable isotope labelling

The quantitative analysis of metabolic networks is a prerequisite for understanding the integration and regulation of plant metabolism and for devising rational approaches for manipulating resource allocation in plants. The analysis of steady state stable isotope labelling experiments using nuclear magnetic resonance (NMR) spectroscopy has developed into a powerful method for determining these fluxes in micro-organisms and its application to heterotrophic plant metabolism is increasing. After an introductory discussion of the well known role of stable isotopes in pathway delineation, the review considers their application to metabolic flux analysis in plants. These applications are divided into two groups – small scale analyses of fluxes through particular pathways and large scale analyses of multiple fluxes through metabolic networks – and the problems caused by the complexity of intermediary metabolism in plants are discussed. It is concluded that metabolic flux analysis provides a powerful method for defining the metabolic phenotype of wild type, mutant and transgenic plants and that its development should be pursued.     

Trichome regulator SlMIXTA‐like directly manipulates primary metabolism in tomato fruit

Trichomes are storage compartments for specialized metabolites in many plant species. In trichome, plant primary metabolism is significantly changed, providing substrates for downstream secondary metabolism. However, little is known of how plants coordinate trichome formation and primary metabolism regulation. In this report, tomato (Solanum lycopersicum) trichome regulator SlMIXTA‐like is indicated as a metabolic regulation gene by mGWAS analysis. Overexpression of SlMIXTA‐like in tomato fruit enhances trichome formation. In addition, SlMIXTA‐like can directly bind to the promoter region of gene encoding 3‐deoxy‐7‐phosphoheptulonate synthase (SlDAHPS) to activate its expression. Induction of SlDAHPS expression enhances shikimate pathway activities and provides substrates for downstream secondary metabolism. Our data provide direct evidence that trichome regulator can directly manipulate primary metabolism, in which way plants can coordinate metabolic regulation and the formation of storage compartments for specialized metabolites. The newly identified SlMIXTA‐like can be used for future metabolic engineering.     

番茄果实早期发育的分子生理机制研究进展

番茄(Solanum lycopersicum)是目前世界上种植面积最广且最受欢迎的蔬菜作物之一, 也是肉果及茄科的重要模式植物。番茄果实发育主要分为早期果实发育和果实成熟2个时期, 但果实形态结构和大小主要决定于早期果实发育时期。该文围绕番茄早期果实发育时期植物激素、细胞周期、表观遗传和源库代谢等多方面调控的分子机制进行了综述, 旨在认识植物生长与发育的基本生物学问题及促进基础理论研究成果在生产中应用。     

A method for accounting for maintenance costs in flux balance analysis improves the prediction of plant cell metabolic phenotypes under stress conditions

Flux balance models of metabolism generally utilize synthesis of biomass as the main determinant of intracellular fluxes. However, the biomass constraint alone is not sufficient to predict realistic fluxes in central heterotrophic metabolism of plant cells because of the major demand on the energy budget due to transport costs and cell maintenance. This major limitation can be addressed by incorporating transport steps into the metabolic model and by implementing a procedure that uses Pareto optimality analysis to explore the trade‐off between ATP and NADPH production for maintenance. This leads to a method for predicting cell maintenance costs on the basis of the measured flux ratio between the oxidative steps of the oxidative pentose phosphate pathway and glycolysis. We show that accounting for transport and maintenance costs substantially improves the accuracy of fluxes predicted from a flux balance model of heterotrophic Arabidopsis cells in culture, irrespective of the objective function used in the analysis. Moreover, when the new method was applied to cells under control, elevated temperature and hyper‐osmotic conditions, only elevated temperature led to a substantial increase in cell maintenance costs. It is concluded that the hyper‐osmotic conditions tested did not impose a metabolic stress, in as much as the metabolic network is not forced to devote more resources to cell maintenance.