Search

Summary: Plants can rapidly mitigate the effects of suboptimal growth environments by phenotypic plasticity of fitness‐traits. While genetic variation for phenotypic plasticity offers the means for breeding climate‐resilient crop lines, accurate genomic prediction models for plasticity of fitness‐related traits are still lacking.Here, we employed condition‐ and accession‐specific metabolic models for 67 Arabidopsis thaliana accessions to dissect and predict plasticity of rosette growth to changes in nitrogen availability.We showed that specific reactions in photorespiration, linking carbon and nitrogen metabolism, as well as key pathways of central carbon metabolism exhibited substantial genetic variation for flux plasticity. We also demonstrated that, in comparison with a genomic prediction model for fresh weight (FW), genomic prediction of growth plasticity improves the predictability of FW under low nitrogen by 58.9% and by additional 15.4% when further integrating data on plasticity of metabolic fluxes.Therefore, the combination of metabolic and statistical modeling provides a stepping stone in understanding the molecular mechanisms and improving the predictability of plasticity for fitness‐related traits. [ABSTRACT FROM AUTHOR]

Lipids play fundamental roles in regulating agronomically important traits. Advances in plant lipid metabolism have until recently largely been based on reductionist approaches, although modulation of its components can have system-wide effects. However, existing models of plant lipid metabolism provide lumped representations, hindering detailed study of component modulation. Here, we present the Plant Lipid Module (PLM) which provides a mechanistic description of lipid metabolism in the Arabidopsis thaliana rosette. We demonstrate that the PLM can be readily integrated in models of A. thaliana Col-0 metabolism, yielding accurate predictions (83%) of single lethal knock-outs and 75% concordance between measured transcript and predicted flux changes under extended darkness. Genome-wide associations with fluxes obtained by integrating the PLM in diel condition- and accession-specific models identify up to 65 candidate genes modulating A. thaliana lipid metabolism. Using mutant lines, we validate up to 40% of the candidates, paving the way for identification of metabolic gene function based on models capturing natural variability in metabolism. The use of automated tools to reconstruct lipid metabolic pathways is not warranted in plants. Here, the authors construct Plant Lipid Module for Arabidopsis rosette using constraint-based modeling, demonstrate its integration in other plant metabolic models, and use it to dissect the genetic architecture of lipid metabolism. [ABSTRACT FROM AUTHOR]

Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determine cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA) that relies on the patterns of isotope labeling of metabolites in the network. For metabolic systems, typical for plants, where all potentially labeled atoms effectively have only one source atom pool, only isotopically nonstationary MFA can provide information about intracellular fluxes. There are several global approaches that implement MFA for an entire metabolic network and estimate, at once, a steady-state flux distribution for all reactions with identifiable fluxes in the network. In contrast, local approaches deal with estimation of fluxes for a subset of reactions, with smaller data demand for flux estimation. Here we present a systematic comparative review and benchmarking of the existing local approaches for isotopically nonstationary MFA. The comparison is conducted with respect to the required data and underlying computational problems solved on a synthetic network example. Furthermore, we benchmark the performance of these approaches in estimating fluxes for a subset of reactions using data obtained from the simulation of nitrogen fluxes in the Arabidopsis thaliana core metabolism. The findings pinpoint practical aspects that need to be considered when applying local approaches for flux estimation in large-scale plant metabolic networks. [ABSTRACT FROM AUTHOR]

Phenotypic plasticity is a heritable trait that provides sessile organisms a strategy to rapidly mitigate negative effects of environmental change. Yet, we have little understanding of the mode of inheritance and genetic architecture of plasticity in different focal traits relevant to agricultural applications. This study builds on our recent discovery of genes controlling temperature-mediated flower size plasticity in Arabidopsis thaliana and focuses on dissecting the mode of inheritance and combining ability of plasticity in the context of plant breeding. We created a full diallel cross using 12 A. thaliana accessions displaying different temperaturemediated flower size plasticities, scored as the fold change between two temperatures. Griffing's analysis of variance in flower size plasticity indicated that non-additive genetic action shapes this trait and pointed at challenges and opportunities when breeding for reduced plasticity. Our findings provide an outlook of flower size plasticity that is important for developing resilient crops for future climates. [ABSTRACT FROM AUTHOR]

Cell division and the resulting changes to the cell organization affect the shape and functionality of all tissues. Thus, understanding the determinants of the tissue-wide changes imposed by cell division is a key question in developmental biology. Here, we use a network representation of live cell imaging data from shoot apical meristems (SAMs) in Arabidopsis thaliana to predict cell division events and their consequences at the tissue level. We show that a support vector machine classifier based on the SAM network properties is predictive of cell division events, with test accuracy of 76%, which matches that based on cell size alone. Furthermore, we demonstrate that the combination of topological and biological properties, including cell size, perimeter, distance and shared cell wall between cells, can further boost the prediction accuracy of resulting changes in topology triggered by cell division. Using our classifiers, we demonstrate the importance of microtubule-mediated cell-to-cell growth coordination in influencing tissue-level topology. Together, the results from our network-based analysis demonstrate a feedback mechanism between tissue topology and cell division in A. thaliana SAMs. [ABSTRACT FROM AUTHOR]

Composition and functions of microbial communities affect important traits in diverse hosts, from crops to humans. Yet, mechanistic understanding of how metabolism of individual microbes is affected by the community composition and metabolite leakage is lacking. Here, we first show that the consensus of automatically generated metabolic reconstructions improves the quality of the draft reconstructions, measured by comparison to reference models. We then devise an approach for gap filling, termed COMMIT, that considers metabolites for secretion based on their permeability and the composition of the community. By applying COMMIT with two soil communities from the Arabidopsis thaliana culture collection, we could significantly reduce the gap-filling solution in comparison to filling gaps in individual reconstructions without affecting the genomic support. Inspection of the metabolic interactions in the soil communities allows us to identify microbes with community roles of helpers and beneficiaries. Therefore, COMMIT offers a versatile fully automated solution for large-scale modelling of microbial communities for diverse biotechnological applications. Author summary: Microbial communities are important in ecology, human health, and crop productivity. However, detailed information on the interactions within natural microbial communities is hampered by the community size, lack of detailed information on the biochemistry of single organisms, and the complexity of interactions between community members. Metabolic models are comprised of biochemical reaction networks based on the genome annotation, and can provide mechanistic insights into community functions. Previous analyses of microbial community models have been performed with high-quality reference models or models generated using a single reconstruction pipeline. However, these models do not contain information on the composition of the community that determines the metabolites exchanged between the community members. In addition, the quality of metabolic models is affected by the reconstruction approach used, with direct consequences on the inferred interactions between community members. Here, we use fully automated consensus reconstructions from four approaches to arrive at functional models with improved genomic support while considering the community composition. We applied our pipeline to two soil communities from the Arabidopsis thaliana culture collection, providing only genome sequences. Finally, we show that the obtained models have 90% genomic support and demonstrate that the derived interactions are corroborated by independent computational predictions. [ABSTRACT FROM AUTHOR]

As autotrophic organisms, plants capture light energy to convert carbon dioxide into ATP, nicotinamide adenine dinucleotide phosphate (NADPH), and sugars, which are essential for the biosynthesis of building blocks, storage, and growth. At night, metabolism and growth can be sustained by mobilizing carbon (C) reserves. In response to changing environmental conditions, such as light-dark cycles, the small-molecule regulation of enzymatic activities is critical for reprogramming cellular metabolism. We have recently demonstrated that proteogenic dipeptides, protein degradation products, act as metabolic switches at the interface of proteostasis and central metabolism in both plants and yeast. Dipeptides accumulate in response to the environmental changes and act via direct binding and regulation of critical enzymatic activities, enabling C flux distribution. Here, we provide evidence pointing to the involvement of dipeptides in the metabolic rewiring characteristics for the day-night cycle in plants. Specifically, we measured the abundance of 13 amino acids and 179 dipeptides over short- (SD) and long-day (LD) diel cycles, each with different light intensities. Of the measured dipeptides, 38 and eight were characterized by day-night oscillation in SD and LD, respectively, reaching maximum accumulation at the end of the day and then gradually falling in the night. Not only the number of dipeptides, but also the amplitude of the oscillation was higher in SD compared with LD conditions. Notably, rhythmic dipeptides were enriched in the glucogenic amino acids that can be converted into glucose. Considering the known role of Target of Rapamycin (TOR) signaling in regulating both autophagy and metabolism, we subsequently investigated whether diurnal fluctuations of dipeptides levels are dependent on the TOR Complex (TORC). The Raptor1b mutant (raptor1b), known for the substantial reduction of TOR kinase activity, was characterized by the augmented accumulation of dipeptides, which is especially pronounced under LD conditions. We were particularly intrigued by the group of 16 dipeptides, which, based on their oscillation under SD conditions and accumulation in raptor1b , can be associated with limited C availability or photoperiod. By mining existing protein-metabolite interaction data, we delineated putative protein interactors for a representative dipeptide Pro-Gln. The obtained list included enzymes of C and amino acid metabolism, which are also linked to the TORC-mediated metabolic network. Based on the obtained results, we speculate that the diurnal accumulation of dipeptides contributes to its metabolic adaptation in response to changes in C availability. We hypothesize that dipeptides would act as alternative respiratory substrates and by directly modulating the activity of the focal enzymes. [ABSTRACT FROM AUTHOR]

Nitrogen (N) is fundamental to plant growth, development and yield. Genes underlying N utilization and assimilation are well‐characterized, but mechanisms underpinning plasticity of different phenotypes in response to N remain elusive. Here, using Arabidopsis thaliana accessions, we dissected the genetic architecture of plasticity in early and late rosette diameter, flowering time and yield, in response to three levels of N in the soil. Furthermore, we found that the plasticity in levels of primary metabolites were related with the plasticities of the studied traits. Genome‐wide association analysis identified three significant associations for phenotypic plasticity, one for early rosette diameter and two for flowering time. We confirmed that the gene At1g19880, hereafter named as PLASTICITY OF ROSETTE TO NITROGEN 1 (PROTON1), encoding for a regulator of chromatin condensation 1 (RCC1) family protein, conferred plasticity of rosette diameter in response to N. Treatment of PROTON1 T‐DNA line with salt implied that the reduced plasticity of early rosette diameter was not a general growth response to stress. We further showed that plasticities of growth and flowering‐related traits differed between environmental cues, indicating decoupled genetic programs regulating these traits. Our findings provide a prospective to identify genes that stabilize performance under fluctuating environments. Genome‐wide and mutant analysis revealed that a RCC1‐family protein controls plasticity of rosette diameter in response to N in soil. Furthermore, our results showed that plasticity of traits in response to different environments are decoupled. [ABSTRACT FROM AUTHOR]

Motivation Genomic selection (GS) is currently deemed the most effective approach to speed up breeding of agricultural varieties. It has been recognized that consideration of multiple traits in GS can improve accuracy of prediction for traits of low heritability. However, since GS forgoes statistical testing with the idea of improving predictions, it does not facilitate mechanistic understanding of the contribution of particular single nucleotide polymorphisms (SNP). Results Here, we propose a L 2 , 1 -norm regularized multivariate regression model and devise a fast and efficient iterative optimization algorithm, called L 2 , 1 -joint, applicable in multi-trait GS. The usage of the L 2 , 1 -norm facilitates variable selection in a penalized multivariate regression that considers the relation between individuals, when the number of SNPs is much larger than the number of individuals. The capacity for variable selection allows us to define master regulators that can be used in a multi-trait GS setting to dissect the genetic architecture of the analyzed traits. Our comparative analyses demonstrate that the proposed model is a favorable candidate compared to existing state-of-the-art approaches. Prediction and variable selection with datasets from Brassica napus , wheat and Arabidopsis thaliana diversity panels are conducted to further showcase the performance of the proposed model. Availability and implementation : The model is implemented using R programming language and the code is freely available from https://github.com/alainmbebi/L21-norm-GS. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Ribosome biogenesis is tightly associated to plant metabolism due to the usage of ribosomes in the synthesis of proteins necessary to drive metabolic pathways. Given the central role of ribosome biogenesis in cell physiology, it is important to characterize the impact of different components involved in this process on plant metabolism. Double mutants of the Arabidopsis thaliana cytosolic 60S maturation factors REIL1 and REIL2 do not resume growth after shift to moderate 10 ∘ C chilling conditions. To gain mechanistic insights into the metabolic effects of this ribosome biogenesis defect on metabolism, we developed TC-iReMet2, a constraint-based modelling approach that integrates relative metabolomics and transcriptomics time-course data to predict differential fluxes on a genome-scale level. We employed TC-iReMet2 with metabolomics and transcriptomics data from the Arabidopsis Columbia 0 wild type and the reil1-1 reil2-1 double mutant before and after cold shift. We identified reactions and pathways that are highly altered in a mutant relative to the wild type. These pathways include the Calvin–Benson cycle, photorespiration, gluconeogenesis, and glycolysis. Our findings also indicated differential NAD(P)/NAD(P)H ratios after cold shift. TC-iReMet2 allows for mechanistic hypothesis generation and interpretation of system biology experiments related to metabolic fluxes on a genome-scale level. [ABSTRACT FROM AUTHOR]

SUMMARY: 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. Significance Statement: By using a constraint‐based modeling approach that integrates publically available proteomics data from Arabidopsis thaliana under 10 different conditions, we estimate maximal in vivo enzyme catalytic rates. Our study provides a comprehensive parameterization of a large‐scale model of Arabidopsis metabolism. [ABSTRACT FROM AUTHOR]

Abiotic stresses cause oxidative damage in plants. Here, we demonstrate that foliar application of an extract from the seaweed Ascophyllum nodosum, SuperFifty (SF), largely prevents paraquat (PQ)-induced oxidative stress in Arabidopsis thaliana. While PQ-stressed plants develop necrotic lesions, plants pre-treated with SF (i.e., primed plants) were unaffected by PQ. Transcriptome analysis revealed induction of reactive oxygen species (ROS) marker genes, genes involved in ROS-induced programmed cell death, and autophagy-related genes after PQ treatment. These changes did not occur in PQ-stressed plants primed with SF. In contrast, upregulation of several carbohydrate metabolism genes, growth, and hormone signaling as well as antioxidant-related genes were specific to SF-primed plants. Metabolomic analyses revealed accumulation of the stress-protective metabolite maltose and the tricarboxylic acid cycle intermediates fumarate and malate in SF-primed plants. Lipidome analysis indicated that those lipids associated with oxidative stress-induced cell death and chloroplast degradation, such as triacylglycerols (TAGs), declined upon SF priming. Our study demonstrated that SF confers tolerance to PQ-induced oxidative stress in A. thaliana, an effect achieved by modulating a range of processes at the transcriptomic, metabolic, and lipid levels. [ABSTRACT FROM AUTHOR]

Summary: Plastid ribosomes are very similar in structure and function to the ribosomes of their bacterial ancestors. Since ribosome biogenesis is not thermodynamically favorable under biological conditions it requires the activity of many assembly factors. Here we have characterized a homolog of bacterial RsgA in Arabidopsis thaliana and show that it can complement the bacterial homolog. Functional characterization of a strong mutant in Arabidopsis revealed that the protein is essential for plant viability, while a weak mutant produced dwarf, chlorotic plants that incorporated immature pre‐16S ribosomal RNA into translating ribosomes. Physiological analysis of the mutant plants revealed smaller, but more numerous, chloroplasts in the mesophyll cells, reduction of chlorophyll a and b, depletion of proplastids from the rib meristem and decreased photosynthetic electron transport rate and efficiency. Comparative RNA sequencing and proteomic analysis of the weak mutant and wild‐type plants revealed that various biotic stress‐related, transcriptional regulation and post‐transcriptional modification pathways were repressed in the mutant. Intriguingly, while nuclear‐ and chloroplast‐encoded photosynthesis‐related proteins were less abundant in the mutant, the corresponding transcripts were increased, suggesting an elaborate compensatory mechanism, potentially via differentially active retrograde signaling pathways. To conclude, this study reveals a chloroplast ribosome assembly factor and outlines the transcriptomic and proteomic responses of the compensatory mechanism activated during decreased chloroplast function. Significance Statement: AtRsgA is an assembly factor necessary for maturation of the small subunit of the chloroplast ribosome. Depletion of AtRsgA leads to dwarfed, chlorotic plants, a decrease of mature 16S rRNA and smaller, but more numerous, chloroplasts. Large‐scale transcriptomic and proteomic analysis revealed that chloroplast‐encoded and ‐targeted proteins were less abundant, while the corresponding transcripts were increased in the mutant. We analyze the transcriptional responses of several retrograde signaling pathways to suggest the mechanism underlying this compensatory response. [ABSTRACT FROM AUTHOR]

Summary: Recent advances in gene function prediction rely on ensemble approaches that integrate results from multiple inference methods to produce superior predictions. Yet, these developments remain largely unexplored in plants. We have explored and compared two methods to integrate 10 gene co‐function networks for Arabidopsis thaliana and demonstrate how the integration of these networks produces more accurate gene function predictions for a larger fraction of genes with unknown function. These predictions were used to identify genes involved in mitochondrial complex I formation, and for five of them, we confirmed the predictions experimentally. The ensemble predictions are provided as a user‐friendly online database, EnsembleNet. The methods presented here demonstrate that ensemble gene function prediction is a powerful method to boost prediction performance, whereas the EnsembleNet database provides a cutting‐edge community tool to guide experimentalists. [ABSTRACT FROM AUTHOR]

Hybrid necrosis is a common type of hybrid incompatibility in plants. This phenomenon is caused by deleterious epistatic interactions, resulting in spontaneous activation of plant defenses associated with leaf necrosis, stunted growth and reduced fertility in hybrids. Specific combinations of alleles of ACCELERATED CELL DEATH 6 ( ACD6) have been shown to be a common cause of hybrid necrosis in Arabidopsis thaliana. Increased ACD6 activity confers broad-spectrum resistance against biotrophic pathogens but reduces biomass production., We generated 996 crosses among individuals derived from a single collection area around Tübingen (Germany) and screened them for hybrid necrosis. Necrotic hybrids were further investigated by genetic linkage, amiRNA silencing, genomic complementation and metabolic profiling. Restriction site associated DNA (RAD)-sequencing was used to understand genetic diversity in the collection sites containing necrosis-inducing alleles., Novel combinations of ACD6 alleles found in neighbouring stands were found to activate the A. thaliana immune system. In contrast to what we observed in controlled conditions, necrotic hybrids did not show reduced fitness in the field. Metabolic profiling revealed changes associated with the activation of the immune system in ACD6-dependent hybrid necrosis., This study expands our current understanding of the active role of ACD6 in mediating trade-offs between defense responses and growth in A. thaliana. [ABSTRACT FROM AUTHOR]

Integrative studies of plant growth require spatially and temporally resolved information from high-throughput imaging systems. However, analysis and interpretation of conventional two-dimensional images is complicated by the three-dimensional nature of shoot architecture and by changes in leaf position over time, termed hyponasty. To solve this problem, Phytotyping4D uses a light-field camera that simultaneously provides a focus image and a depth image, which contains distance information about the object surface. Our automated pipeline segments the focus images, integrates depth information to reconstruct the three-dimensional architecture, and analyses time series to provide information about the relative expansion rate, the timing of leaf appearance, hyponastic movement, and shape for individual leaves and the whole rosette. Phytotyping4D was calibrated and validated using discs of known sizes, and plants tilted at various orientations. Information from this analysis was integrated into the pipeline to allow error assessment during routine operation. To illustrate the utility of Phytotyping4D, we compare diurnal changes in Arabidopsis thaliana wild-type Col-0 and the starchless pgm mutant. Compared to Col-0, pgm showed very low relative expansion rate in the second half of the night, a transiently increased relative expansion rate at the onset of light period, and smaller hyponastic movement including delayed movement after dusk, both at the level of the rosette and individual leaves. Our study introduces light-field camera systems as a tool to accurately measure morphological and growth-related features in plants. [ABSTRACT FROM AUTHOR]

Plants as sessile organisms cannot escape their environment and have to adapt to any changes in the availability of sunlight and nutrients. The quantification of synthesis costs of metabolites, in terms of consumed energy, is a prerequisite to understand trade-offs arising from energetic limitations. Here, we examine the energy consumption of amino acid synthesis in Arabidopsis thaliana. To quantify these costs in terms of the energy equivalent ATP, we introduce an improved cost measure based on flux balance analysis and apply it to three state-of-the-art metabolic reconstructions to ensure robust results. We present the first systematic in silico analysis of the effect of nitrogen supply (nitrate/ammonium) on individual amino acid synthesis costs as well as of the effect of photoautotrophic and heterotrophic growth conditions, integrating day/night-specific regulation. Our results identify nitrogen supply as a key determinant of amino acid costs, in agreement with experimental evidence. In addition, the association of the determined costs with experimentally observed growth patterns suggests that metabolite synthesis costs are involved in shaping regulation of plant growth. Finally, we find that simultaneous uptake of both nitrogen sources can lead to efficient utilization of energy source, which may be the result of evolutionary optimization. [ABSTRACT FROM AUTHOR]

Large-scale modeling of plant metabolism provides the possibility to compare and contrast different cellular and environmental scenarios with the ultimate aim of identifying the components underlying the respective plant behavior. The existing models of Arabidopsis (Arabidopsis thaliana) are top-down assembled, whereby the starting point is the annotated genome, in particular, the metabolic genes. Hence, dead-end metabolites and blocked reactions can arise that are subsequently addressed by using gap-filling algorithms in combination with species-unspecific genes. Here, we present a bottom-up-assembled, large-scale model that relies solely on Arabidopsis-specific annotations and results in the inclusion of only manually curated reactions. While the existing models are largely condition unspecific by employing a single biomass reaction, we provide three biomass compositions that pertain to realistic and frequently examined scenarios: carbon-limiting, nitrogen-limiting, and optimal growth conditions The comparative analysis indicates that the proposed Arabidopsis core model exhibits comparable efficiency in carbon utilization and flexibility to the existing network alternatives. Moreover, the model is utilized to quantify the energy demand of amino acid and enzyme de novo synthesis in photoautotrophic growth conditions. Illustrated by the case of the most abundant protein in the world, Rubisco, we determine its synthesis cost in terms of ATP requirements. This, in turn, allows us to explore the tradeoff between protein synthesis and growth in Arabidopsis. Altogether, the model provides a solid basis for completely species-specific integration of high-throughput data, such as gene expression levels, and for condition-specific investigations of in silico metabolic engineering strategies. [ABSTRACT FROM AUTHOR]

Recent advances in high-throughput omics techniques render it possible to decode the function of genes by using the 'guilt-by-association' principle on biologically meaningful clusters of gene expression data. However, the existing frameworks for biological evaluation of gene clusters are hindered by two bottleneck issues: (1) the choice for the number of clusters, and (2) the external measures which do not take in consideration the structure of the analyzed data and the ontology of the existing biological knowledge. Here, we address the identified bottlenecks by developing a novel framework that allows not only for biological evaluation of gene expression clusters based on existing structured knowledge, but also for prediction of putative gene functions. The proposed framework facilitates propagation of statistical significance at each of the following steps: (1) estimating the number of clusters, (2) evaluating the clusters in terms of novel external structural measures, (3) selecting an optimal clustering algorithm, and (4) predicting gene functions. The framework also includes a method for evaluation of gene clusters based on the structure of the employed ontology. Moreover, our method for obtaining a probabilistic range for the number of clusters is demonstrated valid on synthetic data and available gene expression profiles from Saccharomyces cerevisiae. Finally, we propose a network-based approach for gene function prediction which relies on the clustering of optimal score and the employed ontology. Our approach effectively predicts gene function on the Saccharomyces cerevisiae data set and is also employed to obtain putative gene functions for an Arabidopsis thaliana data set. [ABSTRACT FROM AUTHOR]

Plant behaviors across levels of cellular organization, from biochemical components to tissues and organs, relate and reflect growth habitats. Quantification of the relationship between behaviors captured in various phenotypic characteristics and growth habitats can help reveal molecular mechanisms of plant adaptation. The aim of this article is to introduce the power of using statistics originally developed in the field of geographic variability analysis together with prominent network models in elucidating principles of biological organization. We provide a critical systematic review of the existing statistical and network-based approaches that can be employed to determine patterns of covariation from both uni- and multivariate phenotypic characteristics in plants. We demonstrate that parameter-independent network-based approaches result in robust insights about phenotypic covariation. These insights can be quantified and tested by applying well-established statistics combining the network structure with the phenotypic characteristics. We show that the reviewed network-based approaches are applicable from the level of genes to the study of individuals in a population of Arabidopsis thaliana. Finally, we demonstrate that the patterns of covariation can be generalized to quantifiable biological principles of organization. Therefore, these network-based approaches facilitate not only interpretation of large-scale data sets, but also prediction of biochemical and biological behaviors based on measurable characteristics. [ABSTRACT FROM AUTHOR]

Natural genetic diversity provides a powerful tool to study the complex interrelationship between metabolism and growth. Profiling of metabolic traits combined with network-based and statistical analyses allow the comparison of conditions and identification of sets of traits that predict biomass. However, it often remains unclear why a particular set of metabolites is linked with biomass and to what extent the predictive model is applicable beyond a particular growth condition. A panel of 97 genetically diverse Arabidopsis (Arabidopsis thaliana) accessions was grown in near-optimal carbon and nitrogen supply, restricted carbon supply, and restricted nitrogen supply and analyzed for biomass and 54 metabolic traits. Correlation-based metabolic networks were generated from the genotype-dependent variation in each condition to reveal sets of metabolites that show coordinated changes across accessions. The networks were largely specific for a single growth condition. Partial least squares regression from metabolic traits allowed prediction of biomass within and, slightly more weakly, across conditions (cross-validated Pearson correlations in the range of 0.27-0.58 and 0.21-0.51 and P values in the range of <0.001-<0.13 and <0.001-<0.023, respectively). Metabolic traits that correlate with growth or have a high weighting in the partial least squares regression were mainly condition specific and often related to the resource that restricts growth under that condition. Linear mixed-model analysis using the combined metabolic traits from all growth conditions as an input indicated that inclusion of random effects for the conditions improves predictions of biomass. Thus, robust prediction of biomass across a range of conditions requires condition-specific measurement of metabolic traits to take account of environment-dependent changes of the underlying networks. [ABSTRACT FROM AUTHOR]

Nitrogen (N) is central for plant growth, and metabolic plasticity can provide a strategy to respond to changing N availability. We showed that two local A. thaliana populations exhibited differential plasticity in the compounds of photorespiratory and starch degradation pathways in response to three N conditions. Association of metabolite levels with growth‐related and fitness traits indicated that controlled plasticity in these pathways could contribute to local adaptation and play a role in plant evolution. [ABSTRACT FROM AUTHOR]

Motivation: Comprehensive understanding of cellular processes requires development of approaches which consider the energetic balances in the cell. The existing approaches that address this problem are based on defining energy-equivalent costs which do not include the effects of a changing environment. By incorporating these effects, one could provide a framework for integrating 'omics' data from various levels of the system in order to provide interpretations with respect to the energy state and to elicit conclusions about putative global energy-related response mechanisms in the cell. Results: Here we define a cost measure for amino acid synthesis based on flux balance analysis of a genome-scale metabolic network, and develop methods for its integration with proteomics and metabolomics data. This is a first measure which accounts for the effect of different environmental conditions. We applied this approach to a genome-scale network of Arabidopsis thaliana and calculated the costs for all amino acids and proteins present in the network under light and dark conditions. Integration of function and process ontology terms in the analysis of protein abundances and their costs indicates that, during the night, the cell favors cheaper proteins compared with the light environment. However, this does not imply that there is squandering of resources during the day. The results from the association analysis between the costs, levels and well-defined expenses of amino acid synthesis, indicate that our approach not only captures the adjustment made at the switch of conditions, but also could explain the anticipation of resource usage via a global energy-related regulatory mechanism of amino acid and protein synthesis. [ABSTRACT FROM AUTHOR]

Discrimination of metabolic models based on high throughput metabolomics data, reflecting various internal and external perturbations, is essential for identifying the components that contribute to the emerging behavior of metabolic processes. Here, we investigate 12 different models of the mitochondrial electron transport chain (ETC) in Arabidopsis thaliana during dark-induced senescence in order to elucidate the alternative substrates to this metabolic pathway. Our findings demonstrate that the coupling of the proposed computational approach, based on dynamic flux balance analysis, with time-resolved metabolomics data results in model-based confirmations of the hypotheses that, during dark-induced senescence in Arabidopsis, (i) under conditions where the main substrate for the ETC are not fully available, isovaleryl-CoA dehydrogenase and 2-hydroxyglutarate dehydrogenase are able to donate electrons to the ETC, (ii) phytanoyl-CoA does not act even as an indirect substrate of the electron transfer flavoprotein/electron-transfer flavoprotein:ubiquinone oxidoreductase complex, and (iii) the mitochondrial γ-aminobutyric acid transporter has functional significance in maintaining mitochondrial metabolism. Our study provides a basic framework for future in silico studies of alternative pathways in mitochondrial metabolism under extended darkness whereby the role of its components can be computationally discriminated based on available molecular profile data. [ABSTRACT FROM AUTHOR]