AraGEM, a genome-scale reconstruction of the primary metabolic network in Arabidopsis

Leaf day respiration: More than just catabolic CO2 production in the light

Day respiration is a net flux resulting from several CO2‐generating and CO2‐fixing reactions, not only related to catabolism but also to anabolism. We review pieces of evidence that decarboxylating reactions are partly fed by carbon sources disconnected from current photosynthesis and how they reflect various metabolic pathways.

Resource allocation modeling for autonomous prediction of plant cell phenotypes

Predicting the plant cell response in complex environmental conditions is a challenge in plant biology. Here we developed a resource allocation model of cellular and molecular scale for the leaf photosynthetic cell of Arabidopsis thaliana, based on the Resource Balance Analysis (RBA) constraint-based modeling framework. The RBA model contains the metabolic network and the major macromolecular processes involved in the plant cell growth and survival and localized in cellular compartments. We simulated the model for varying environmental conditions of temperature, irradiance, partial pressure of CO2 and O2, and compared RBA predictions to known resource distributions and quantitative phenotypic traits such as the relative growth rate, the C:N ratio, and finally to the empirical characteristics of CO2 fixation given by the well-established Farquhar model. In comparison to other standard constraint-based modeling methods like Flux Balance Analysis, the RBA model makes accurate quantitative predictions without the need for empirical constraints. Altogether, we show that RBA significantly improves the autonomous prediction of plant cell phenotypes in complex environmental conditions, and provides mechanistic links between the genotype and the phenotype of the plant cell.

Sugar transporters spatially organize microbiota colonization along the longitudinal root axis of Arabidopsis

Plant roots are functionally heterogeneous in cellular architecture, transcriptome profile, metabolic state, and microbial immunity. We hypothesized that axial differentiation may also impact spatial colonization by root microbiota along the root axis. We developed two growth systems, ArtSoil and CD-Rhizotron, to grow and then dissect Arabidopsis thaliana roots into three segments. We demonstrate that distinct endospheric and rhizosphere bacterial communities colonize the segments, supporting the hypothesis of microbiota differentiation along the axis. Root metabolite profiling of each segment reveals differential metabolite enrichment and specificity. Bioinformatic analyses and GUS histochemistry indicate microbe-induced accumulation of SWEET2, 4, and 12 sugar uniporters. Profiling of root segments from sweet mutants shows altered spatial metabolic profiles and reorganization of endospheric root microbiota. This work reveals the interdependency between root metabolites and microbial colonization and the contribution of SWEETs to spatial diversity and stability of microbial ecosystem.

A multi-organ maize metabolic model connects temperature stress with energy production and reducing power generation

Climate change has adversely affected maize productivity. Thereby, a holistic understanding of metabolic crosstalk among its organs is important to address this issue. Thus, we reconstructed the first multi-organ maize metabolic model, iZMA6517, and contextualized it with heat and cold stress transcriptomics data using expression distributed reaction flux measurement (EXTREAM) algorithm. Furthermore, implementing metabolic bottleneck analysis on contextualized models revealed differences between these stresses. While both stresses had reducing power bottlenecks, heat stress had additional energy generation bottlenecks. We also performed thermodynamic driving force analysis, revealing thermodynamics-reducing power-energy generation axis dictating the nature of temperature stress responses. Thus, a temperature-tolerant maize ideotype can be engineered by leveraging the proposed thermodynamics-reducing power-energy generation axis. We experimentally inoculated maize root with a beneficial mycorrhizal fungus, Rhizophagus irregularis, and as a proof-of-concept demonstrated its efficacy in alleviating temperature stress. Overall, this study will guide the engineering effort of temperature stress-tolerant maize ideotypes.

Quantifying the impact of dynamic plant-environment interactions on metabolic regulation

A plant's genome encodes enzymes, transporters and many other proteins which constitute metabolism. Interactions of plants with their environment shape their growth, development and resilience towards adverse conditions. Although genome sequencing technologies and applications have experienced triumphantly rapid development during the last decades, enabling nowadays a fast and cheap sequencing of full genomes, prediction of metabolic phenotypes from genotype × environment interactions remains, at best, very incomplete. The main reasons are a lack of understanding of how different levels of molecular organisation depend on each other, and how they are constituted and expressed within a setup of growth conditions. Phenotypic plasticity, e.g., of the genetic model plant Arabidopsis thaliana, has provided important insights into plant-environment interactions and the resulting genotype x phenotype relationships. Here, we summarize previous and current findings about plant development in a changing environment and how this might be shaped and reflected in metabolism and its regulation. We identify current challenges in the study of plant development and metabolic regulation and provide an outlook of how methodological workflows might support the application of findings made in model systems to crops and their cultivation.

The first multi-tissue genome-scale metabolic model of a woody plant highlights suberin biosynthesis pathways in Quercus suber

Over the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behaviour at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871). The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyse the pathways associated with the synthesis of suberin monomers, namely the acyl-lipids, phenylpropanoids, isoprenoids, and flavonoids production. The models developed in this work provide a systematic overview of the metabolism of Q. suber, including its secondary metabolism pathways and cork formation.

Reconciliation and evolution of Penicillium rubens genome-scale metabolic networks-What about specialised metabolism?

In recent years, genome sequencing of filamentous fungi has revealed a high proportion of specialised metabolites with growing pharmaceutical interest. However, detecting such metabolites through in silico genome analysis does not necessarily guarantee their expression under laboratory conditions. However, one plausible strategy for enabling their production lies in modifying the growth conditions. Devising a comprehensive experimental design testing in different culture environments is time-consuming and expensive. Therefore, using in silico modelling as a preliminary step, such as Genome-Scale Metabolic Network (GSMN), represents a promising approach to predicting and understanding the observed specialised metabolite production in a given organism. To address these questions, we reconstructed a new high-quality GSMN for the Penicillium rubens Wisconsin 54-1255 strain, a commonly used model organism. Our reconstruction, iPrub22, adheres to current convention standards and quality criteria, incorporating updated functional annotations, orthology searches with different GSMN templates, data from previous reconstructions, and manual curation steps targeting primary and specialised metabolites. With a MEMOTE score of 74% and a metabolic coverage of 45%, iPrub22 includes 5,192 unique metabolites interconnected by 5,919 reactions, of which 5,033 are supported by at least one genomic sequence. Of the metabolites present in iPrub22, 13% are categorised as belonging to specialised metabolism. While our high-quality GSMN provides a valuable resource for investigating known phenotypes expressed in P. rubens, our analysis identifies bottlenecks related, in particular, to the definition of what is a specialised metabolite, which requires consensus within the scientific community. It also points out the necessity of accessible, standardised and exhaustive databases of specialised metabolites. These questions must be addressed to fully unlock the potential of natural product production in P. rubens and other filamentous fungi. Our work represents a foundational step towards the objective of rationalising the production of natural products through GSMN modelling.

Genome-scale metabolic model led engineering of Nothapodytes nimmoniana plant cells for high camptothecin production

Camptothecin (CPT) is a vital monoterpene indole alkaloid used in anti-cancer therapeutics. It is primarily derived from Camptotheca acuminata and Nothapodytes nimmoniana plants that are indigenous to Southeast Asia. Plants have intricate metabolic networks and use them to produce secondary metabolites such as CPT, which is a prerequisite for rational metabolic engineering design to optimize their production. By reconstructing metabolic models, we can predict plant metabolic behavior, facilitating the selection of suitable approaches and saving time, cost, and energy, over traditional hit and trial experimental approaches. In this study, we reconstructed a genome-scale metabolic model for N. nimmoniana (NothaGEM iSM1809) and curated it using experimentally obtained biochemical data. We also used in silico tools to identify and rank suitable enzyme targets for overexpression and knockout to maximize camptothecin production. The predicted over-expression targets encompass enzymes involved in the camptothecin biosynthesis pathway, including strictosidine synthase and geraniol 10-hydroxylase, as well as targets related to plant metabolism, such as amino acid biosynthesis and the tricarboxylic acid cycle. The top-ranked knockout targets included reactions responsible for the formation of folates and serine, as well as the conversion of acetyl CoA and oxaloacetate to malate and citrate. One of the top-ranked overexpression targets, strictosidine synthase, was chosen to generate metabolically engineered cell lines of N. nimmoniana using Agrobacterium tumefaciens-mediated transformation. The transformed cell line showed a 5-fold increase in camptothecin production, with a yield of up to 5 µg g-1.

A Protocol for the Automatic Construction of Highly Curated Genome-Scale Models of Human Metabolism

Genome-scale metabolic models (GEMs) have emerged as a tool to understand human metabolism from a holistic perspective with high relevance in the study of many diseases and in the metabolic engineering of human cell lines. GEM building relies on either automated processes that lack manual refinement and result in inaccurate models or manual curation, which is a time-consuming process that limits the continuous update of reliable GEMs. Here, we present a novel algorithm-aided protocol that overcomes these limitations and facilitates the continuous updating of highly curated GEMs. The algorithm enables the automatic curation and/or expansion of existing GEMs or generates a highly curated metabolic network based on current information retrieved from multiple databases in real time. This tool was applied to the latest reconstruction of human metabolism (Human1), generating a series of the human GEMs that improve and expand the reference model and generating the most extensive and comprehensive general reconstruction of human metabolism to date. The tool presented here goes beyond the current state of the art and paves the way for the automatic reconstruction of a highly curated, up-to-date GEM with high potential in computational biology as well as in multiple fields of biological science where metabolism is relevant.

Toward mechanistic modeling and rational engineering of plant respiration

Plant respiration not only provides energy to support all cellular processes, including biomass production, but also plays a major role in the global carbon cycle. Therefore, modulation of plant respiration can be used to both increase the plant yield and mitigate the effects of global climate change. Mechanistic modeling of plant respiration at sufficient biochemical detail can provide key insights for rational engineering of this process. Yet, despite its importance, plant respiration has attracted considerably less modeling effort in comparison to photosynthesis. In this update review, we highlight the advances made in modeling of plant respiration, emphasizing the gradual but important change from phenomenological to models based on first principles. We also provide a detailed account of the existing resources that can contribute to resolving the challenges in modeling plant respiration. These resources point at tangible improvements in the representation of cellular processes that contribute to CO2 evolution and consideration of kinetic properties of underlying enzymes to facilitate mechanistic modeling. The update review emphasizes the need to couple biochemical models of respiration with models of acclimation and adaptation of respiration for their effective usage in guiding breeding efforts and improving terrestrial biosphere models tailored to future climate scenarios.

Antioxidant Green Factories: Toward Sustainable Production of Vitamin E in Plant In Vitro Cultures

Vitamin E is a dietary supplement synthesized only by photosynthetic organisms and, hence, is an essential vitamin for human well-being. Because of the ever-increasing demand for natural vitamin E and limitations in existing synthesis modes, attempts to improve its yield using plant in vitro cultures have gained traction in recent years. With inflating industrial production costs, integrative approaches to conventional bioprocess optimization is the need of the hour for multifold vitamin E productivity enhancement. In this review, we briefly discuss the structure, isomers, and important metabolic routes of biosynthesis for vitamin E in plants. We then emphasize its vital role in human health and its industrial applications and highlight the market demand and supply. We illustrate the advantages of in vitro plant cell/tissue culture cultivation as an alternative to current commercial production platforms for natural vitamin E. We touch upon the conventional vitamin E metabolic pathway engineering strategies, such as single/multigene overexpression and chloroplast engineering. We highlight the recent progress in plant systems biology to rationally identify metabolic bottlenecks and knockout targets in the vitamin E biosynthetic pathway. We then discuss bioprocess optimization strategies for sustainable vitamin E production, including media/process optimization, precursor/elicitor addition, and scale-up to bioreactors. We culminate the review with a short discussion on kinetic modeling to predict vitamin E production in plant cell cultures and suggestions on sustainable green extraction methods of vitamin E for reduced environmental impact. This review will be of interest to a wider research fraternity, including those from industry and academia working in the field of plant cell biology, plant biotechnology, and bioprocess engineering for phytochemical enhancement.

Integration of physiologically relevant photosynthetic energy flows into whole genome models of light-driven metabolism

Characterizing photosynthetic productivity is necessary to understand the ecological contributions and biotechnology potential of plants, algae, and cyanobacteria. Light capture efficiency and photophysiology have long been characterized by measurements of chlorophyll fluorescence dynamics. However, these investigations typically do not consider the metabolic network downstream of light harvesting. By contrast, genome-scale metabolic models capture species-specific metabolic capabilities but have yet to incorporate the rapid regulation of the light harvesting apparatus. Here, we combine chlorophyll fluorescence parameters defining photosynthetic and non-photosynthetic yield of absorbed light energy with a metabolic model of the pennate diatom Phaeodactylum tricornutum. This integration increases the model predictive accuracy regarding growth rate, intracellular oxygen production and consumption, and metabolic pathway usage. Through the quantification of excess electron transport, we uncover the sequential activation of non-radiative energy dissipation processes, cross-compartment electron shuttling, and non-photochemical quenching as the rapid photoacclimation strategy in P. tricornutum. Interestingly, the photon absorption thresholds that trigger the transition between these mechanisms were consistent at low and high incident photon fluxes. We use this understanding to explore engineering strategies for rerouting cellular resources and excess light energy towards bioproducts in silico. Overall, we present a methodology for incorporating a common, informative data type into computational models of light-driven metabolism and show its utilization within the design-build-test-learn cycle for engineering of photosynthetic organisms.

Genome-scale modeling of the primary-specialized metabolism interface

Environmental challenges and development require plants to reallocate resources between primary and specialized metabolites to survive. Genome-scale metabolic models, which map carbon flux through metabolic pathways, are a valuable tool in the study of tradeoffs that arise at this interface. Due to annotation gaps, models that characterize all the enzymatic steps in individual specialized pathways and their linkages to each other and to central carbon metabolism are difficult to construct. Recent studies have successfully curated subsystems of specialized metabolism and characterized the interfaces where flux is diverted to the precursors of glucosinolates, terpenes, and anthocyanins. Although advances in metabolite profiling can help to constrain models at this interface, quantitative analysis remains challenging because of the different timescales on which specialized metabolites from constitutive and reactive pathways accumulate.

The topology of genome-scale metabolic reconstructions unravels independent modules and high network flexibility

The topology of metabolic networks is recognisably modular with modules weakly connected apart from sharing a pool of currency metabolites. Here, we defined modules as sets of reversible reactions isolated from the rest of metabolism by irreversible reactions except for the exchange of currency metabolites. Our approach identifies topologically independent modules under specific conditions associated with different metabolic functions. As case studies, the E.coli iJO1366 and Human Recon 2.2 genome-scale metabolic models were split in 103 and 321 modules respectively, displaying significant correlation patterns in expression data. Finally, we addressed a fundamental question about the metabolic flexibility conferred by reversible reactions: “Of all Directed Topologies (DTs) defined by fixing directions to all reversible reactions, how many are capable of carrying flux through all reactions?”. Enumeration of the DTs for iJO1366 model was performed using an efficient depth-first search algorithm, rejecting infeasible DTs based on mass-imbalanced and loopy flux patterns. We found the direction of 79% of reversible reactions must be defined before all directions in the network can be fixed, granting a high degree of flexibility.

Genome-scale metabolic reconstructions represent all biochemical reactions that an organism can accomplish. These reconstructions are complex and often difficult to study in great detail. A way to overcome this limitation is to focus on specific pathways or subsystems. We present a novel method to identify metabolic modules based on the network topology. The method relies on reaction directions and ignores currency metabolites, which artificially connect distant metabolic reactions. In this way, topologically independent modules are built, where inputs and outputs are controlled by irreversible reactions. The method is automatic and unbiased, and, the result is a set of condition specific modules with defined metabolic functions. As a proof-of-concept we generated biologically relevant modules for the E.coli and Human genome-scale metabolic reconstructions supported by transcriptomic data. Finally, we applied the novel approach to study the network flexibility conferred by reversible reactions. In the case of the E. coli model, we found that the direction of 79% of structurally reversible reactions (those not directionally constrained by surrounding irreversible reactions) must be fixed to determine all the reaction directions in the network. Therefore, reversible reactions operate practically independent of each other.

Flavonoid deficiency disrupts redox homeostasis and terpenoid biosynthesis in glandular trichomes of tomato

Glandular trichomes (GTs) are epidermal structures that provide the first line of chemical defense against arthropod herbivores and other biotic threats. The most conspicuous structure on leaves of cultivated tomato (Solanum lycopersicum) is the type-VI GT (tVI-GT), which accumulates both flavonoids and volatile terpenoids. Although these classes of specialized metabolites are derived from distinct metabolic pathways, previous studies with a chalcone isomerase 1 (CHI1)-deficient mutant called anthocyanin free (af) showed that flavonoids are required for terpenoid accumulation in tVI-GTs. Here, we combined global transcriptomic and proteomic analyses of isolated trichomes as a starting point to show that the lack of CHI1 is associated with reduced levels of terpenoid biosynthetic transcripts and enzymes. The flavonoid deficiency in af trichomes also resulted in the upregulation of abiotic stress-responsive genes associated with DNA damage and repair. Several lines of biochemical and genetic evidence indicate that the terpenoid defect in af mutants is specific for the tVI-GT and is associated with the absence of bulk flavonoids rather than loss of CHI1 per se. A newly developed genome-scale model of metabolism in tomato tVI-GTs helped identify metabolic imbalances caused by the loss of flavonoid production. We provide evidence that flavonoid deficiency in this cell type leads to increased production of reactive oxygen species (ROS), which may impair terpenoid biosynthesis. Collectively, our findings support a role for flavonoids as ROS-scavenging antioxidants in GTs.

Critical assessment of genome-scale metabolic models of Arabidopsis thaliana

Genome-scale metabolic models (GEMs) have enabled researchers to perform systems-level studies of living organisms. Flux balance analysis (FBA), as a constraint-based technique, enables computation of reaction fluxes and prediction of the metabolic phenotypes of a cell under a set of specified conditions. The quality of a GEM is important for obtaining accurate predictions. In this study, we evaluated the quality of five available GEMs for Arabidopsis thaliana from various points of views. To do this, we inspected some of their important features, including the number of reactions with well-defined gene-protein-reaction rules, number of blocked reactions, mass-unbalanced reactions, prediction accuracy in the simulation of key metabolic functions and existence of erroneous energy generating cycles (EGCs). All of the models were found to include some mass-unbalanced reactions. Moreover, four out of five models were found to include EGCs. However, Aracell includes the maximum number of blocked reactions, which suggests the presence of several incomplete pathways. These results clearly show that simulation by using these models may result in erroneous predictions and all of the publicly available GEMs for A. thaliana require extensive curations before being applied in practice.

Dissecting the metabolic reprogramming of maize root under nitrogen-deficient stress conditions

The growth and development of maize (Zea mays L.) largely depends on its nutrient uptake through the root. Hence, studying its growth, response, and associated metabolic reprogramming to stress conditions is becoming an important research direction. A genome-scale metabolic model (GSM) for the maize root was developed to study its metabolic reprogramming under nitrogen stress conditions. The model was reconstructed based on the available information from KEGG, UniProt, and MaizeCyc. Transcriptomics data derived from the roots of hydroponically grown maize plants were used to incorporate regulatory constraints in the model and simulate nitrogen-non-limiting (N+) and nitrogen-deficient (N-) condition. Model-predicted flux-sum variability analysis achieved 70% accuracy compared with the experimental change of metabolite levels. In addition to predicting important metabolic reprogramming in central carbon, fatty acid, amino acid, and other secondary metabolism, maize root GSM predicted several metabolites (l-methionine, l-asparagine, l-lysine, cholesterol, and l-pipecolate) playing a regulatory role in the root biomass growth. Furthermore, this study revealed eight phosphatidylcholine and phosphatidylglycerol metabolites which, even though not coupled with biomass production, played a key role in the increased biomass production under N-deficient conditions. Overall, the omics-integrated GSM provides a promising tool to facilitate stress condition analysis for maize root and engineer better stress-tolerant maize genotypes.

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

Data Management and Modeling in Plant Biology

The study of plant-environment interactions is a multidisciplinary research field. With the emergence of quantitative large-scale and high-throughput techniques, amount and dimensionality of experimental data have strongly increased. Appropriate strategies for data storage, management, and evaluation are needed to make efficient use of experimental findings. Computational approaches of data mining are essential for deriving statistical trends and signatures contained in data matrices. Although, current biology is challenged by high data dimensionality in general, this is particularly true for plant biology. Plants as sessile organisms have to cope with environmental fluctuations. This typically results in strong dynamics of metabolite and protein concentrations which are often challenging to quantify. Summarizing experimental output results in complex data arrays, which need computational statistics and numerical methods for building quantitative models. Experimental findings need to be combined by computational models to gain a mechanistic understanding of plant metabolism. For this, bioinformatics and mathematics need to be combined with experimental setups in physiology, biochemistry, and molecular biology. This review presents and discusses concepts at the interface of experiment and computation, which are likely to shape current and future plant biology. Finally, this interface is discussed with regard to its capabilities and limitations to develop a quantitative model of plant-environment interactions.

Modeling photosynthetic resource allocation connects physiology with evolutionary environments

The regulation of resource allocation in biological systems observed today is the cumulative result of natural selection in ancestral and recent environments. To what extent are observed resource allocation patterns in different photosynthetic types optimally adapted to current conditions, and to what extent do they reflect ancestral environments? Here, we explore these questions for C₃, C₄, and C₃–C₄ intermediate plants of the model genus Flaveria. We developed a detailed mathematical model of carbon fixation, which accounts for various environmental parameters and for energy and nitrogen partitioning across photosynthetic components. This allows us to assess environment-dependent plant physiology and performance as a function of resource allocation patterns. Models of C₄ plants optimized for conditions experienced by evolutionary ancestors perform better than models accounting for experimental growth conditions, indicating low phenotypic plasticity. Supporting this interpretation, the model predicts that C₄ species need to re-allocate more nitrogen between photosynthetic components than C₃ species to adapt to new environments. We thus hypothesize that observed resource distribution patterns in C₄ plants still reflect optimality in ancestral environments, allowing the quantitative inference of these environments from today’s plants. Our work allows us to quantify environmental effects on photosynthetic resource allocation and performance in the light of evolutionary history.

Elucidation of trophic interactions in an unusual single-cell nitrogen-fixing symbiosis using metabolic modeling

Marine nitrogen-fixing microorganisms are an important source of fixed nitrogen in oceanic ecosystems. The colonial cyanobacterium Trichodesmium and diatom symbionts were thought to be the primary contributors to oceanic N₂ fixation until the discovery of the unusual uncultivated symbiotic cyanobacterium UCYN-A (Candidatus Atelocyanobacterium thalassa). UCYN-A has atypical metabolic characteristics lacking the oxygen-evolving photosystem II, the tricarboxylic acid cycle, the carbon-fixation enzyme RuBisCo and de novo biosynthetic pathways for a number of amino acids and nucleotides. Therefore, it is obligately symbiotic with its single-celled haptophyte algal host. UCYN-A receives fixed carbon from its host and returns fixed nitrogen, but further insights into this symbiosis are precluded by both UCYN-A and its host being uncultured. In order to investigate how this syntrophy is coordinated, we reconstructed bottom-up genome-scale metabolic models of UCYN-A and its algal partner to explore possible trophic scenarios, focusing on nitrogen fixation and biomass synthesis. Since both partners are uncultivated and only the genome sequence of UCYN-A is available, we used the phylogenetically related Chrysochromulina tobin as a proxy for the host. Through the use of flux balance analysis (FBA), we determined the minimal set of metabolites and biochemical functions that must be shared between the two organisms to ensure viability and growth. We quantitatively investigated the metabolic characteristics that facilitate daytime N₂ fixation in UCYN-A and possible oxygen-scavenging mechanisms needed to create an anaerobic environment to allow nitrogenase to function. This is the first application of an FBA framework to examine the tight metabolic coupling between uncultivated microbes in marine symbiotic communities and provides a roadmap for future efforts focusing on such specialized systems.

Reduction of dinitrogen gas to biologically useful forms via nitrogen fixation is a key component of the biogeochemical cycle. In the marine environment, the cyanobacteria UCYN-A (Candidatus Atelocyanobacterium thalassa) has been found to be a primary contributor to biological nitrogen fixation at a global scale. UCYN-A exhibits a highly streamlined genome which lacks genes coding for essential cyanobacterial processes such as the energy-generating TCA cycle, oxygen-producing photosystem II, the carbon-fixing RuBisCo and de novo production pathways for numerous amino acids and nucleotides. Thus, it exists in a symbiosis with unicellular planktonic algae where it exchanges fixed nitrogen for fixed carbon with its host. However, both UCYN-A and its symbiotic partner remain uncultured under laboratory conditions. This necessitates implementing a computational approach to glean insights into UCYN-A’s unique physiology and metabolic processes governing the symbiotic association. To this end, we develop a constraints-based framework that infers all possible trophic scenarios consistent with the observed data. Possible mechanisms employed by UCYN-A to accommodate diazotrophy with daytime carbon fixation by the host (i.e., two mutually incompatible processes) are also elucidated. We envision that the developed framework using UCYN-A and its algal host will be used as a roadmap and motivate the study of similarly unique microbial systems in the future.

Characterization of effects of genetic variants via genome-scale metabolic modelling

Genome-scale metabolic networks for model plants and crops in combination with approaches from the constraint-based modelling framework have been used to predict metabolic traits and design metabolic engineering strategies for their manipulation. With the advances in technologies to generate large-scale genotyping data from natural diversity panels and other populations, genome-wide association and genomic selection have emerged as statistical approaches to determine genetic variants associated with and predictive of traits. Here, we review recent advances in constraint-based approaches that integrate genetic variants in genome-scale metabolic models to characterize their effects on reaction fluxes. Since some of these approaches have been applied in organisms other than plants, we provide a critical assessment of their applicability particularly in crops. In addition, we further dissect the inferred effects of genetic variants with respect to reaction rate constants, abundances of enzymes, and concentrations of metabolites, as main determinants of reaction fluxes and relate them with their combined effects on complex traits, like growth. Through this systematic review, we also provide a roadmap for future research to increase the predictive power of statistical approaches by coupling them with mechanistic models of metabolism.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Transcriptome integrated metabolic modeling of carbon assimilation underlying storage root development in cassava

The existing genome-scale metabolic model of carbon metabolism in cassava storage roots, rMeCBM, has proven particularly resourceful in exploring the metabolic basis for the phenotypic differences between high and low-yield cassava cultivars. However, experimental validation of predicted metabolic fluxes by carbon labeling is quite challenging. Here, we incorporated gene expression data of developing storage roots into the basic flux-balance model to minimize infeasible metabolic fluxes, denoted as rMeCBMx, thereby improving the plausibility of the simulation and predictive power. Three different conceptual algorithms, GIMME, E-Flux, and HPCOF were evaluated. The rMeCBMx-HPCOF model outperformed others in predicting carbon fluxes in the metabolism of storage roots and, in particular, was highly consistent with transcriptome of high-yield cultivars. The flux prediction was improved through the oxidative pentose phosphate pathway in cytosol, as has been reported in various studies on root metabolism, but hardly captured by simple FBA models. Moreover, the presence of fluxes through cytosolic glycolysis and alanine biosynthesis pathways were predicted with high consistency with gene expression levels. This study sheds light on the importance of prediction power in the modeling of complex plant metabolism. Integration of multi-omics data would further help mitigate the ill-posed problem of constraint-based modeling, allowing more realistic simulation.

Environment-coupled models of leaf metabolism

The plant leaf is the main site of photosynthesis. This process converts light energy and inorganic nutrients into chemical energy and organic building blocks for the biosynthesis and maintenance of cellular components and to support the growth of the rest of the plant. The leaf is also the site of gas–water exchange and due to its large surface, it is particularly vulnerable to pathogen attacks. Therefore, the leaf's performance and metabolic modes are inherently determined by its interaction with the environment. Mathematical models of plant metabolism have been successfully applied to study various aspects of photosynthesis, carbon and nitrogen assimilation and metabolism, aided suggesting metabolic intervention strategies for optimized leaf performance, and gave us insights into evolutionary drivers of plant metabolism in various environments. With the increasing pressure to improve agricultural performance in current and future climates, these models have become important tools to improve our understanding of plant–environment interactions and to propel plant breeders efforts. This overview article reviews applications of large-scale metabolic models of leaf metabolism to study plant–environment interactions by means of flux-balance analysis. The presented studies are organized in two ways — by the way the environment interactions are modelled — via external constraints or data-integration and by the studied environmental interactions — abiotic or biotic.

Biotechnological strategies for enhancing heavy metal tolerance in neglected and underutilized legume crops: A comprehensive review

Contamination of agricultural land and water by heavy metals due to rapid industrialization and urbanization including various natural processes have become one of the major constraints to crop growth and productivity. Several studies have reported that to counteract heavy metal stress, plants should be able to maneuver various physiological, biochemical and molecular processes to improve their growth and development under heavy metal stress. With the advent of modern biotechnological tools and techniques it is now possible to tailor legume and other plants overexpressing stress-induced genes, transcription factors, proteins, and metabolites that are directly involved in heavy metal stress tolerance. This review provides an in-depth overview of various biotechnological approaches and/or strategies that can be used for enhancing detoxification of the heavy metals by stimulating phytoremediation processes. Synthetic biology tools involved in the engineering of legume and other crop plants against heavy metal stress tolerance are also discussed herewith some pioneering examples where synthetic biology tools that have been used to modify plants for specific traits. Also, CRISPR based genetic engineering of plants, including their role in modulating the expression of several genes/ transcription factors in the improvement of abiotic stress tolerance and phytoremediation ability using knockdown and knockout strategies has also been critically discussed.

Dissecting the regulatory roles of ORM proteins in the sphingolipid pathway of plants

Sphingolipids are a vital component of plant cellular endomembranes and carry out multiple functional and regulatory roles. Different sphingolipid species confer rigidity to the membrane structure, facilitate trafficking of secretory proteins, and initiate programmed cell death. Although the regulation of the sphingolipid pathway is yet to be uncovered, increasing evidence has pointed to orosomucoid proteins (ORMs) playing a major regulatory role and potentially interacting with a number of components in the pathway, including both enzymes and sphingolipids. However, experimental exploration of new regulatory interactions is time consuming and often infeasible. In this work, a computational approach was taken to address this challenge. A metabolic network of the sphingolipid pathway in plants was reconstructed. The steady-state rates of reactions in the network were then determined through measurements of growth and cellular composition of the different sphingolipids in Arabidopsis seedlings. The Ensemble modeling framework was modified to accurately account for activation mechanisms and subsequently used to generate sets of kinetic parameters that converge to the measured steady-state fluxes in a thermodynamically consistent manner. In addition, the framework was appended with an additional module to automate screening the parameters and to output models consistent with previously reported network responses to different perturbations. By analyzing the network's response in the presence of different combinations of regulatory mechanisms, the model captured the experimentally observed repressive effect of ORMs on serine palmitoyltransferase (SPT). Furthermore, predictions point to a second regulatory role of ORM proteins, namely as an activator of class II (or LOH1 and LOH3) ceramide synthases. This activating role was found to be modulated by the concentration of free ceramides, where an accumulation of these sphingolipid species dampened the activating effect of ORMs on ceramide synthase. The predictions pave the way for future guided experiments and have implications in engineering crops with higher biotic stress tolerance.

Automated engineering of synthetic metabolic pathways for efficient biomanufacturing

Metabolic engineering involves the engineering and optimization of processes from single-cell to fermentation in order to increase production of valuable chemicals for health, food, energy, materials and others. A systems approach to metabolic engineering has gained traction in recent years thanks to advances in strain engineering, leading to an accelerated scaling from rapid prototyping to industrial production. Metabolic engineering is nowadays on track towards a truly manufacturing technology, with reduced times from conception to production enabled by automated protocols for DNA assembly of metabolic pathways in engineered producer strains. In this review, we discuss how the success of the metabolic engineering pipeline often relies on retrobiosynthetic protocols able to identify promising production routes and dynamic regulation strategies through automated biodesign algorithms, which are subsequently assembled as embedded integrated genetic circuits in the host strain. Those approaches are orchestrated by an experimental design strategy that provides optimal scheduling planning of the DNA assembly, rapid prototyping and, ultimately, brings forward an accelerated Design-Build-Test-Learn cycle and the overall optimization of the biomanufacturing process. Achieving such a vision will address the increasingly compelling demand in our society for delivering valuable biomolecules in an affordable, inclusive and sustainable bioeconomy.

Plant Biosystems Design Research Roadmap 1.0

Human life intimately depends on plants for food, biomaterials, health, energy, and a sustainable environment. Various plants have been genetically improved mostly through breeding, along with limited modification via genetic engineering, yet they are still not able to meet the ever-increasing needs, in terms of both quantity and quality, resulting from the rapid increase in world population and expected standards of living. A step change that may address these challenges would be to expand the potential of plants using biosystems design approaches. This represents a shift in plant science research from relatively simple trial-and-error approaches to innovative strategies based on predictive models of biological systems. Plant biosystems design seeks to accelerate plant genetic improvement using genome editing and genetic circuit engineering or create novel plant systems through de novo synthesis of plant genomes. From this perspective, we present a comprehensive roadmap of plant biosystems design covering theories, principles, and technical methods, along with potential applications in basic and applied plant biology research. We highlight current challenges, future opportunities, and research priorities, along with a framework for international collaboration, towards rapid advancement of this emerging interdisciplinary area of research. Finally, we discuss the importance of social responsibility in utilizing plant biosystems design and suggest strategies for improving public perception, trust, and acceptance.

Characterization of maximal enzyme catalytic rates in central metabolism of Arabidopsis thaliana

Availability of plant-specific enzyme kinetic data is scarce, limiting the predictive power of metabolic models and precluding identification of genetic factors of enzyme properties. Enzyme kinetic data are measured in vitro, often under non-physiological conditions, and conclusions elicited from modeling warrant caution. Here we estimate maximal in vivo catalytic rates for 168 plant enzymes, including photosystems I and II, cytochrome-b6f complex, ATP-citrate synthase, sucrose-phosphate synthase as well as enzymes from amino acid synthesis with previously undocumented enzyme kinetic data in BRENDA. The estimations are obtained by integrating condition-specific quantitative proteomics data, maximal rates of selected enzymes, growth measurements from Arabidopsis thaliana rosette with and fluxes through canonical pathways in a constraint-based model of leaf metabolism. In comparison to findings in Escherichia coli, we demonstrate weaker concordance between the plant-specific in vitro and in vivo enzyme catalytic rates due to a low degree of enzyme saturation. This is supported by the finding that concentrations of nicotinamide adenine dinucleotide (phosphate), adenosine triphosphate and uridine triphosphate, calculated based on our maximal in vivo catalytic rates, and available quantitative metabolomics data are below reported KM values and, therefore, indicate undersaturation of respective enzymes. Our findings show that genome-wide profiling of enzyme kinetic properties is feasible in plants, paving the way for understanding resource allocation.

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.

The importance of thiamine (vitamin B1) in plant health: From crop yield to biofortification

Ensuring that people have access to sufficient and nutritious food is necessary for a healthy life and the core tenet of food security. With the global population set to reach 9.8 billion by 2050, and the compounding effects of climate change, the planet is facing challenges that necessitate significant and rapid changes in agricultural practices. In the effort to provide food in terms of calories, the essential contribution of micronutrients (vitamins and minerals) to nutrition is often overlooked. Here, we focus on the importance of thiamine (vitamin B1) in plant health and discuss its impact on human health. Vitamin B1 is an essential dietary component, and deficiencies in this micronutrient underlie several diseases, notably nervous system disorders. The predominant source of dietary vitamin B1 is plant-based foods. Moreover, vitamin B1 is also vital for plants themselves, and its benefits in plant health have received less attention than in the human health sphere. In general, vitamin B1 is well-characterized for its role as a coenzyme in metabolic pathways, particularly those involved in energy production and central metabolism, including carbon assimilation and respiration. Vitamin B1 is also emerging as an important component of plant stress responses, and several noncoenzyme roles of this vitamin are being characterized. We summarize the importance of vitamin B1 in plants from the perspective of food security, including its roles in plant disease resistance, stress tolerance, and crop yield, and review the potential benefits of biofortification of crops with increased vitamin B1 content to improve human health.

The genome of Ectocarpus subulatus - A highly stress-tolerant brown alga

Brown algae are multicellular photosynthetic stramenopiles that colonize marine rocky shores worldwide. Ectocarpus sp. Ec32 has been established as a genomic model for brown algae. Here we present the genome and metabolic network of the closely related species, Ectocarpus subulatus Kützing, which is characterized by high abiotic stress tolerance. Since their separation, both strains show new traces of viral sequences and the activity of large retrotransposons, which may also be related to the expansion of a family of chlorophyll-binding proteins. Further features suspected to contribute to stress tolerance include an expanded family of heat shock proteins, the reduction of genes involved in the production of halogenated defence compounds, and the presence of fewer cell wall polysaccharide-modifying enzymes. Overall, E. subulatus has mainly lost members of gene families down-regulated in low salinities, and conserved those that were up-regulated in the same condition. However, 96% of genes that differed between the two examined Ectocarpus species, as well as all genes under positive selection, were found to encode proteins of unknown function. This underlines the uniqueness of brown algal stress tolerance mechanisms as well as the significance of establishing E. subulatus as a comparative model for future functional studies.

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.

In silico prediction of enzymatic reactions catalyzed by acid phosphatases

In present work, we describe a methodology for prediction of an enzymatic reaction for which no experimental data are available except for a gene sequence. As a challenging case, we have developed the method for identifying the putative substrates of monoester phosphatases, commonly known as acid phosphatase enzymes, which have no strong substrate specificity. Finding a preferable substrate for each one is an important task to unravel pathways involved in plant phosphate metabolism. Having used an Arabidopsis thaliana haloacid dehalogenase (HAD)-related acid phosphatases, HRP9, with an experimentally known structure and preferred substrate as an instance, we firstly predicted the 3 D-structure of HRP1 for subsequent analysis. Then, molecular docking was used to find the best protein interaction with a ligand existing in a set of possible substrates compiled from genome scale metabolic networks of A. thaliana based on binding energy, binding mode as well as the distance between phosphoric ester and cofactor, Mg²⁺, localized in the active site of HRP1. Molecular dynamics simulation ratified stable protein-ligand complex model. Our analysis predicted HRP1 preferably bind to pyridoxamine-5'-phosphate (PMP). Thus, it is deduced that the conversion of PMP to pyridoxamine must be catalyzed by HRP1. This procedure is expected to make a reliable pipeline to predict the enzymatic reactions catalyzed by acid phosphatases. Taken as a whole, it could be applicable for discovery of the interacting ligands, inhibitors as well as interacting proteins which limits lab works or used for gap filling in biosystems.Communicated by Ramaswamy H. Sarma.

SNPeffect: identifying functional roles of SNPs using metabolic networks

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome-wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single-nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype-to-phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme-coding genes. To this end, we also constructed a genome-scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth-determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.

Using automated reasoning to explore the metabolism of unconventional organisms: a first step to explore host-microbial interactions

Systems modelled in the context of molecular and cellular biology are difficult to represent with a single calibrated numerical model. Flux optimisation hypotheses have shown tremendous promise to accurately predict bacterial metabolism but they require a precise understanding of metabolic reactions occurring in the considered species. Unfortunately, this information may not be available for more complex organisms or non-cultured microorganisms such as those evidenced in microbiomes with metagenomic techniques. In both cases, flux optimisation techniques may not be applicable to elucidate systems functioning. In this context, we describe how automatic reasoning allows relevant features of an unconventional biological system to be identified despite a lack of data. A particular focus is put on the use of Answer Set Programming, a logic programming paradigm with combinatorial optimisation functionalities. We describe its usage to over-approximate metabolic responses of biological systems and solve gap-filling problems. In this review, we compare steady-states and Boolean abstractions of metabolic models and illustrate their complementarity via applications to the metabolic analysis of macro-algae. Ongoing applications of this formalism explore the emerging field of systems ecology, notably elucidating interactions between a consortium of microbes and a host organism. As the first step in this field, we will illustrate how the reduction in microbiotas according to expected metabolic phenotypes can be addressed with gap-filling problems.

Introducing an Optimization- and explicit Runge-Kutta- based Approach to Perform Dynamic Flux Balance Analysis

In this work we introduce the generalized Optimization- and explicit Runge-Kutta-based Approach (ORKA) to perform dynamic Flux Balance Analysis (dFBA), which is numerically more accurate and computationally tractable than existing approaches. ORKA is applied to a four-tissue (leaf, root, seed, and stem) model of Arabidopsis thaliana, p-ath773, uniquely capturing the core-metabolism of several stages of growth from seedling to senescence at hourly intervals. Model p-ath773 has been designed to show broad agreement with published plant-scale properties such as mass, maintenance, and senescence, yet leaving reaction-level behavior unconstrainted. Hence, it serves as a framework to study the reaction-level behavior necessary for observed plant-scale behavior. Two such case studies of reaction-level behavior include the lifecycle progression of sulfur metabolism and the diurnal flow of water throughout the plant. Specifically, p-ath773 shows how transpiration drives water flow through the plant and how water produced by leaf tissue metabolism may contribute significantly to transpired water. Investigation of sulfur metabolism elucidates frequent cross-compartment exchange of a standing pool of amino acids which is used to regulate the proton flow. Overall, p-ath773 and ORKA serve as scaffolds for dFBA-based lifecycle modeling of plants and other systems to further broaden the scope of in silico metabolic investigation.

Drought Stress Responses in Context-Specific Genome-Scale Metabolic Models of Arabidopsis thaliana

Drought perturbs metabolism in plants and limits their growth. Because drought stress on crops affects their yields, understanding the complex adaptation mechanisms evolved by plants against drought will facilitate the development of drought-tolerant crops for agricultural use. In this study, we examined the metabolic pathways of Arabidopsis thaliana which respond to drought stress by omics-based in silico analyses. We proposed an analysis pipeline to understand metabolism under specific conditions based on a genome-scale metabolic model (GEM). Context-specific GEMs under drought and well-watered control conditions were reconstructed using transcriptome data and examined using metabolome data. The metabolic fluxes throughout the metabolic network were estimated by flux balance analysis using the context-specific GEMs. We used in silico methods to identify an important reaction contributing to biomass production and clarified metabolic reaction responses under drought stress by comparative analysis between drought and control conditions. This proposed pipeline can be applied in other studies to understand metabolic changes under specific conditions using Arabidopsis GEM or other available plant GEMs.

BioHackathon 2015: Semantics of data for life sciences and reproducible research

We report on the activities of the 2015 edition of the BioHackathon, an annual event that brings together researchers and developers from around the world to develop tools and technologies that promote the reusability of biological data. We discuss issues surrounding the representation, publication, integration, mining and reuse of biological data and metadata across a wide range of biomedical data types of relevance for the life sciences, including chemistry, genotypes and phenotypes, orthology and phylogeny, proteomics, genomics, glycomics, and metabolomics. We describe our progress to address ongoing challenges to the reusability and reproducibility of research results, and identify outstanding issues that continue to impede the progress of bioinformatics research. We share our perspective on the state of the art, continued challenges, and goals for future research and development for the life sciences Semantic Web.

Inferring Biochemical Reactions and Metabolite Structures to Understand Metabolic Pathway Drift

Inferring genome-scale metabolic networks in emerging model organisms is challenged by incomplete biochemical knowledge and partial conservation of biochemical pathways during evolution. Therefore, specific bioinformatic tools are necessary to infer biochemical reactions and metabolic structures that can be checked experimentally. Using an integrative approach combining genomic and metabolomic data in the red algal model Chondrus crispus, we show that, even metabolic pathways considered as conserved, like sterols or mycosporine-like amino acid synthesis pathways, undergo substantial turnover. This phenomenon, here formally defined as "metabolic pathway drift," is consistent with findings from other areas of evolutionary biology, indicating that a given phenotype can be conserved even if the underlying molecular mechanisms are changing. We present a proof of concept with a methodological approach to formalize the logical reasoning necessary to infer reactions and molecular structures, abstracting molecular transformations based on previous biochemical knowledge.

Multi-tissue to whole plant metabolic modelling

Genome-scale metabolic models have been successfully applied to study the metabolism of multiple plant species in the past decade. While most existing genome-scale modelling studies have focussed on studying the metabolic behaviour of individual plant metabolic systems, there is an increasing focus on combining models of multiple tissues or organs to produce multi-tissue models that allow the investigation of metabolic interactions between tissues and organs. Multi-tissue metabolic models were constructed for multiple plants including Arabidopsis, barley, soybean and Setaria. These models were applied to study various aspects of plant physiology including the division of labour between organs, source and sink tissue relationship, growth of different tissues and organs and charge and proton balancing. In this review, we outline the process of constructing multi-tissue genome-scale metabolic models, discuss the strengths and challenges in using multi-tissue models, review the current status of plant multi-tissue and whole plant metabolic models and explore the approaches for integrating genome-scale metabolic models into multi-scale plant models.

Computational systems biology of cellular processes in Arabidopsis thaliana: an overview

Systems biology strives for gaining an understanding of biological phenomena by studying the interactions of different parts of a system and integrating the knowledge obtained into the current view of the underlying processes. This is achieved by a tight combination of quantitative experimentation and computational modeling. While there is already a large quantity of systems biology studies describing human, animal and especially microbial cell biological systems, plant biology has been lagging behind for many years. However, in the case of the model plant Arabidopsis thaliana, the steadily increasing amount of information on the levels of its genome, proteome and on a variety of its metabolic and signalling pathways is progressively enabling more researchers to construct models for cellular processes for the plant, which in turn encourages more experimental data to be generated, showing also for plant sciences how fruitful systems biology research can be. In this review, we provide an overview over some of these recent studies which use different systems biological approaches to get a better understanding of the cell biology of A. thaliana. The approaches used in these are genome-scale metabolic modeling, as well as kinetic modeling of metabolic and signalling pathways. Furthermore, we selected several cases to exemplify how the modeling approaches have led to significant advances or new perspectives in the field.

Reconstruction and analysis of a carbon-core metabolic network for Dunaliella salina

BACKGROUND

The green microalga Dunaliella salina accumulates a high proportion of β-carotene during abiotic stress conditions. To better understand the intracellular flux distribution leading to carotenoid accumulation, this work aimed at reconstructing a carbon core metabolic network for D. salina CCAP 19/18 based on the recently published nuclear genome and its validation with experimental observations and literature data.

RESULTS

The reconstruction resulted in a network model with 221 reactions and 212 metabolites within three compartments: cytosol, chloroplast and mitochondrion. The network was implemented in the MATLAB toolbox CellNetAnalyzer and checked for feasibility. Furthermore, a flux balance analysis was carried out for different light and nutrient uptake rates. The comparison of the experimental knowledge with the model prediction revealed that the results of the stoichiometric network analysis are plausible and in good agreement with the observed behavior. Accordingly, our model provides an excellent tool for investigating the carbon core metabolism of D. salina.

CONCLUSIONS

The reconstructed metabolic network of D. salina presented in this work is able to predict the biological behavior under light and nutrient stress and will lead to an improved process understanding for the optimized production of high-value products in microalgae.

Comparative Metabolic Network Flux Analysis to Identify Differences in Cellular Metabolism

Metabolic network flux analysis uses genome-scale metabolic reconstructions to integrate transcriptomics, proteomics, and/or metabolomics data to allow for comprehensive interpretation of genotype to metabolic phenotype relationships. The compilation of many Constraint-based model analysis methods into one MATLAB package, the COBRAtoolbox, has opened the possibility of using these methods to the many biologists with some knowledge of the commonly used statistical program, MATLAB. Here we outline the steps required to take a published genome-scale metabolic reconstruction and interrogate its consistency and biological feasibility. Subsequently, we demonstrate how mRNA expression data and metabolomics data, relating to one or more cell types or biological contexts, can be applied to constrain and generate metabolic models descriptive of metabolic flux phenotypes. Finally, we describe the comparison of the resulting models and model outputs with the aim of identifying metabolic biomarkers and changes in cellular metabolism.

Using Genome-Scale Metabolic Networks for Analysis, Visualization, and Integration of Targeted Metabolomics Data

Interpretation of metabolomics data in the context of biological pathways is important to gain knowledge about underlying metabolic processes. In this chapter we present methods to analyze genome-scale models (GSMs) and metabolomics data together. This includes reading and mining of GSMs using the SBTab format to retrieve information on genes, reactions, and metabolites. Furthermore, the chapter showcases the generation of metabolic pathway maps using the Escher tool, which can be used for data visualization. Lastly, approaches to constrain flux balance analysis (FBA) by metabolomics data are presented.

Evolution of C4 photosynthesis predicted by constraint-based modelling

Constraint-based modelling (CBM) is a powerful tool for the analysis of evolutionary trajectories. Evolution, especially evolution in the distant past, is not easily accessible to laboratory experimentation. Modelling can provide a window into evolutionary processes by allowing the examination of selective pressures which lead to particular optimal solutions in the model. To study the evolution of C4 photosynthesis from a ground state of C3 photosynthesis, we initially construct a C3 model. After duplication into two cells to reflect typical C4 leaf architecture, we allow the model to predict the optimal metabolic solution under various conditions. The model thus identifies resource limitation in conjunction with high photorespiratory flux as a selective pressure relevant to the evolution of C4. It also predicts that light availability and distribution play a role in guiding the evolutionary choice of possible decarboxylation enzymes. The data shows evolutionary CBM in eukaryotes predicts molecular evolution with precision.

Genome-Scale Metabolic Networks Shed Light on the Carotenoid Biosynthesis Pathway in the Brown Algae Saccharina japonica and Cladosiphon okamuranus

Understanding growth mechanisms in brown algae is a current scientific and economic challenge that can benefit from the modeling of their metabolic networks. The sequencing of the genomes of Saccharina japonica and Cladosiphon okamuranus has provided the necessary data for the reconstruction of Genome-Scale Metabolic Networks (GSMNs). The same in silico method deployed for the GSMN reconstruction of Ectocarpus siliculosus to investigate the metabolic capabilities of these two algae, was used. Integrating metabolic profiling data from the literature, we provided functional GSMNs composed of an average of 2230 metabolites and 3370 reactions. Based on these GSMNs and previously published work, we propose a model for the biosynthetic pathways of the main carotenoids in these two algae. We highlight, on the one hand, the reactions and enzymes that have been preserved through evolution and, on the other hand, the specificities related to brown algae. Our data further indicate that, if abscisic acid is produced by Saccharina japonica, its biosynthesis pathway seems to be different in its final steps from that described in land plants. Thus, our work illustrates the potential of GSMNs reconstructions for formalizing hypotheses that can be further tested using targeted biochemical approaches.