Histone retention preserves epigenetic marks during heat stress-induced transcriptional memory in plants.

Plants often experience recurrent stressful events, for example, during heat waves. They can be primed by heat stress (HS) to improve the survival of more severe heat stress conditions. At certain genes, sustained expression is induced for several days beyond the initial heat stress. This transcriptional memory is associated with hyper-methylation of histone H3 lysine 4 (H3K4me3), but it is unclear how this is maintained for extended periods. Here, we determined histone turnover by measuring the chromatin association of HS-induced histone H3.3. Genome-wide histone turnover was not homogenous; in particular, H3.3 was retained longer at heat stress memory genes compared to HS-induced non-memory genes during the memory phase. While low nucleosome turnover retained H3K4 methylation, methylation loss did not affect turnover, suggesting that low nucleosome turnover sustains H3K4 methylation, but not vice versa. Together, our results unveil the modulation of histone turnover as a mechanism to retain environmentally mediated epigenetic modifications.

Tyr-Asp inhibition of glyceraldehyde 3-phosphate dehydrogenase affects plant redox metabolism.

How organisms integrate metabolism with the external environment is a central question in biology. Here, we describe a novel regulatory small molecule, a proteogenic dipeptide Tyr-Asp, which improves plant tolerance to oxidative stress by directly interfering with glucose metabolism. Specifically, Tyr-Asp inhibits the activity of a key glycolytic enzyme, glyceraldehyde 3-phosphate dehydrogenase (GAPC), and redirects glucose toward pentose phosphate pathway (PPP) and NADPH production. In line with the metabolic data, Tyr-Asp supplementation improved the growth performance of both Arabidopsis and tobacco seedlings subjected to oxidative stress conditions. Moreover, inhibition of Arabidopsis phosphoenolpyruvate carboxykinase (PEPCK) activity by a group of branched-chain amino acid-containing dipeptides, but not by Tyr-Asp, points to a multisite regulation of glycolytic/gluconeogenic pathway by dipeptides. In summary, our results open the intriguing possibility that proteogenic dipeptides act as evolutionarily conserved small-molecule regulators at the nexus of stress, protein degradation, and metabolism.

Topological properties accurately predict cell division events and organization of shoot apical meristem in Arabidopsis thaliana.

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.

Fluorescent cytoskeletal markers reveal associations between the actin and microtubule cytoskeleton in rice cells.

Rice (Oryza sativa) is one of our main food crops, feeding ∼3.5 billion people worldwide. An increasing number of studies note the importance of the cytoskeleton, including actin filaments and microtubules, in rice development and environmental responses. Yet, reliable in vivo cytoskeleton markers are lacking in rice, which limits our knowledge of cytoskeletal functions in living cells. Therefore, we generated bright fluorescent marker lines of the actin and microtubule cytoskeletons in rice, suitable for live-cell imaging in a wide variety of rice tissues. Using these lines, we show that actin bundles and microtubules engage and co-function during pollen grain development, how the cytoskeletal components are coordinated during root cell development, and that the actin cytoskeleton is robust and facilitates microtubule responses during salt stress. Hence, we conclude that our cytoskeletal marker lines, highlighted by our findings of cytoskeletal associations and dynamics, will substantially further future investigations in rice biology.

Genome-wide association of the metabolic shifts underpinning dark-induced senescence in Arabidopsis.

Dark-induced senescence provokes profound metabolic shifts to recycle nutrients and to guarantee plant survival. To date, research on these processes has largely focused on characterizing mutants deficient in individual pathways. Here, we adopted a time-resolved genome-wide association-based approach to characterize dark-induced senescence by evaluating the photochemical efficiency and content of primary and lipid metabolites at the beginning, or after 3 or 6 days in darkness. We discovered six patterns of metabolic shifts and identified 215 associations with 81 candidate genes being involved in this process. Among these associations, we validated the roles of four genes associated with glycine, galactinol, threonine, and ornithine levels. We also demonstrated the function of threonine and galactinol catabolism during dark-induced senescence. Intriguingly, we determined that the association between tyrosine contents and TYROSINE AMINOTRANSFERASE 1 influences enzyme activity of the encoded protein and transcriptional activity of the gene under normal and dark conditions, respectively. Moreover, the single-nucleotide polymorphisms affecting the expression of THREONINE ALDOLASE 1 and the amino acid transporter gene AVT1B, respectively, only underlie the variation in threonine and glycine levels in the dark. Taken together, these results allow us to present a very detailed model of the metabolic aspects of dark-induced senescence, as well as the process itself.

The genome of Haberlea rhodopensis provides insights into the mechanisms for tolerance to multiple extreme environments.

Haberlea rhodopensis, a resurrection species, is the only plant known to be able to survive multiple extreme environments, including desiccation, freezing temperatures, and long-term darkness. However, the molecular mechanisms underlying tolerance to these stresses are poorly studied. Here, we present a high-quality genome of Haberlea and found that ~ 23.55% of the 44,306 genes are orphan. Comparative genomics analysis identified 89 significantly expanded gene families, of which 25 were specific to Haberlea. Moreover, we demonstrated that Haberlea preserves its resurrection potential even in prolonged complete darkness. Transcriptome profiling of plants subjected to desiccation, darkness, and low temperatures revealed both common and specific footprints of these stresses, and their combinations. For example, PROTEIN PHOSPHATASE 2C (PP2C) genes were substantially induced in all stress combinations, while PHYTOCHROME INTERACTING FACTOR 1 (PIF1) and GROWTH RESPONSE FACTOR 4 (GRF4) were induced only in darkness. Additionally, 733 genes with unknown functions and three genes encoding transcription factors specific to Haberlea were specifically induced/repressed upon combination of stresses, rendering them attractive targets for future functional studies. The study provides a comprehensive understanding of the genomic architecture and reports details of the mechanisms of multi-stress tolerance of this resurrection species that will aid in developing strategies that allow crops to survive extreme and multiple abiotic stresses.

Prediction and integration of metabolite-protein interactions with genome-scale metabolic models.

Metabolites, as small molecules, can act not only as substrates to enzymes, but also as effectors of activity of proteins with different functions, thereby affecting various cellular processes. While several experimental techniques have started to catalogue the metabolite-protein interactions (MPIs) present in different cellular contexts, characterizing the functional relevance of MPIs remains a challenging problem. Computational approaches from the constrained-based modeling framework allow for predicting MPIs and integrating their effects in the in silico analysis of metabolic and physiological phenotypes, like cell growth. Here, we provide a classification of all existing constraint-based approaches that predict and integrate MPIs using genome-scale metabolic networks as input. In addition, we benchmark the performance of the approaches to predict MPIs in a comparative study using different features extracted from the model structure and predicted metabolic phenotypes with the state-of-the-art metabolic networks of Escherichia coli and Saccharomyces cerevisiae. Lastly, we provide an outlook for future, feasible directions to expand the consideration of MPIs in constraint-based modeling approaches with wide biotechnological applications.

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.

Protein constraints in genome-scale metabolic models: Data integration, parameter estimation, and prediction of metabolic phenotypes.

Genome-scale metabolic models provide a valuable resource to study metabolism and cell physiology. These models are employed with approaches from the constraint-based modeling framework to predict metabolic and physiological phenotypes. The prediction performance of genome-scale metabolic models can be improved by including protein constraints. The resulting protein-constrained models consider data on turnover numbers (kcat ) and facilitate the integration of protein abundances. In this systematic review, we present and discuss the current state-of-the-art regarding the estimation of kinetic parameters used in protein-constrained models. We also highlight how data-driven and constraint-based approaches can aid the estimation of turnover numbers and their usage in improving predictions of cellular phenotypes. Finally, we identify standing challenges in protein-constrained metabolic models and provide a perspective regarding future approaches to improve the predictive performance.

Gene regulatory network inference using mixed-norms regularized multivariate model with covariance selection.

Despite extensive research efforts, reconstruction of gene regulatory networks (GRNs) from transcriptomics data remains a pressing challenge in systems biology. While non-linear approaches for reconstruction of GRNs show improved performance over simpler alternatives, we do not yet have understanding if joint modelling of multiple target genes may improve performance, even under linearity assumptions. To address this problem, we propose two novel approaches that cast the GRN reconstruction problem as a blend between regularized multivariate regression and graphical models that combine the L2,1-norm with classical regularization techniques. We used data and networks from the DREAM5 challenge to show that the proposed models provide consistently good performance in comparison to contenders whose performance varies with data sets from simulation and experiments from model unicellular organisms Escherichia coli and Saccharomyces cerevisiae. Since the models' formulation facilitates the prediction of master regulators, we also used the resulting findings to identify master regulators over all data sets as well as their plasticity across different environments. Our results demonstrate that the identified master regulators are in line with experimental evidence from the model bacterium E. coli. Together, our study demonstrates that simultaneous modelling of several target genes results in improved inference of GRNs and can be used as an alternative in different applications.

Maximizing multi-reaction dependencies provides more accurate and precise predictions of intracellular fluxes than the principle of parsimony.

Intracellular fluxes represent a joint outcome of cellular transcription and translation and reflect the availability and usage of nutrients from the environment. While approaches from the constraint-based metabolic framework can accurately predict cellular phenotypes, such as growth and exchange rates with the environment, accurate prediction of intracellular fluxes remains a pressing problem. Parsimonious flux balance analysis (pFBA) has become an approach of choice to predict intracellular fluxes by employing the principle of efficient usage of protein resources. Nevertheless, comparative analyses of intracellular flux predictions from pFBA against fluxes estimated from labeling experiments remain scarce. Here, we posited that steady-state flux distributions derived from the principle of maximizing multi-reaction dependencies are of improved accuracy and precision than those resulting from pFBA. To this end, we designed a constraint-based approach, termed complex-balanced FBA (cbFBA), to predict steady-state flux distributions that support the given specific growth rate and exchange fluxes. We showed that the steady-state flux distributions resulting from cbFBA in comparison to pFBA show better agreement with experimentally measured fluxes from 17 Escherichia coli strains and are more precise, due to the smaller space of alternative solutions. We also showed that the same principle holds in eukaryotes by comparing the predictions of pFBA and cbFBA against experimentally derived steady-state flux distributions from 26 knock-out mutants of Saccharomyces cerevisiae. Furthermore, our results showed that intracellular fluxes predicted by cbFBA provide better support for the principle of minimizing metabolic adjustment between mutants and wild types. Together, our findings point that other principles that consider the dynamics and coordination of steady states may govern the distribution of intracellular fluxes.

PARROT: Prediction of enzyme abundances using protein-constrained metabolic models.

Protein allocation determines the activity of cellular pathways and affects growth across all organisms. Therefore, different experimental and machine learning approaches have been developed to quantify and predict protein abundance and how they are allocated to different cellular functions, respectively. Yet, despite advances in protein quantification, it remains challenging to predict condition-specific allocation of enzymes in metabolic networks. Here, using protein-constrained metabolic models, we propose a family of constrained-based approaches, termed PARROT, to predict how much of each enzyme is used based on the principle of minimizing the difference between a reference and an alternative growth condition. To this end, PARROT variants model the minimization of enzyme reallocation using four different (combinations of) distance functions. We demonstrate that the PARROT variant that minimizes the Manhattan distance between the enzyme allocation of a reference and an alternative condition outperforms existing approaches based on the parsimonious distribution of fluxes or enzymes for both Escherichia coli and Saccharomyces cerevisiae. Further, we show that the combined minimization of flux and enzyme allocation adjustment leads to inconsistent predictions. Together, our findings indicate that minimization of protein allocation rather than flux redistribution is a governing principle determining steady-state pathway activity for microorganism grown in alternative growth conditions.

An automated workflow that generates atom mappings for large-scale metabolic models and its application to Arabidopsis thaliana.

Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determines cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA), which relies on the patterns of isotope labeling of metabolites in the network. The application of MFA also requires a stoichiometric model with atom mappings that are currently not available for the majority of large-scale metabolic network models, particularly of plants. While automated approaches such as the Reaction Decoder Toolkit (RDT) can produce atom mappings for individual reactions, tracing the flow of individual atoms of the entire reactions across a metabolic model remains challenging. Here we establish an automated workflow to obtain reliable atom mappings for large-scale metabolic models by refining the outcome of RDT, and apply the workflow to metabolic models of Arabidopsis thaliana. We demonstrate the accuracy of RDT through a comparative analysis with atom mappings from a large database of biochemical reactions, MetaCyc. We further show the utility of our automated workflow by simulating 15 N isotope enrichment and identifying nitrogen (N)-containing metabolites which show enrichment patterns that are informative for flux estimation in future 15 N-MFA studies of A. thaliana. The automated workflow established in this study can be readily expanded to other species for which metabolic models have been established and the resulting atom mappings will facilitate MFA and graph-theoretic structural analyses with large-scale metabolic networks.

Accurate prediction of in vivo protein abundances by coupling constraint-based modelling and machine learning.

Quantification of how different environmental cues affect protein allocation can provide important insights for understanding cell physiology. While absolute quantification of proteins can be obtained by resource-intensive mass-spectrometry-based technologies, prediction of protein abundances offers another way to obtain insights into protein allocation. Here we present CAMEL, a framework that couples constraint-based modelling with machine learning to predict protein abundance for any environmental condition. This is achieved by building machine learning models that leverage static features, derived from protein sequences, and condition-dependent features predicted from protein-constrained metabolic models. Our findings demonstrate that CAMEL results in excellent prediction of protein allocation in E. coli (average Pearson correlation of at least 0.9), and moderate performance in S. cerevisiae (average Pearson correlation of at least 0.5). Therefore, CAMEL outperformed contending approaches without using molecular read-outs from unseen conditions and provides a valuable tool for using protein allocation in biotechnological applications.

Optogenetic control of Cdc48 for dynamic metabolic engineering in yeast.

Dynamic metabolic engineering is a strategy to switch key metabolic pathways in microbial cell factories from biomass generation to accumulation of target products. Here, we demonstrate that optogenetic intervention in the cell cycle of budding yeast can be used to increase production of valuable chemicals, such as the terpenoid ß-carotene or the nucleoside analog cordycepin. We achieved optogenetic cell-cycle arrest in the G2/M phase by controlling activity of the ubiquitin-proteasome system hub Cdc48. To analyze the metabolic capacities in the cell cycle arrested yeast strain, we studied their proteomes by timsTOF mass spectrometry. This revealed widespread, but highly distinct abundance changes of metabolic key enzymes. Integration of the proteomics data in protein-constrained metabolic models demonstrated modulation of fluxes directly associated with terpenoid production as well as metabolic subsystems involved in protein biosynthesis, cell wall synthesis, and cofactor biosynthesis. These results demonstrate that optogenetically triggered cell cycle intervention is an option to increase the yields of compounds synthesized in a cellular factory by reallocation of metabolic resources.

Predicting plasticity of rosette growth and metabolic fluxes in Arabidopsis thaliana.

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.

Models and molecular mechanisms for trade-offs in the context of metabolism.

Accumulating evidence for trade-offs involving metabolic traits has demonstrated their importance in the evolution of organisms. Metabolic models with different levels of complexity have already been considered when investigating mechanisms that explain various metabolic trade-offs. Here we provide a systematic review of modelling approaches that have been used to study and explain trade-offs between: (i) the kinetic properties of individual enzymes, (ii) rates of metabolic reactions, (iii) the rate and yield of metabolic pathways and networks, (iv) different metabolic objectives in single organisms and in metabolic communities, and (v) metabolic concentrations. In providing insights into the mechanisms underlying these five types of metabolic trade-offs obtained from constraint-based metabolic modelling, we emphasize the relationship of metabolic trade-offs to the classical black box Y-model that provides a conceptual explanation for resource acquisition-allocation trade-offs. In addition, we identify several pressing concerns and offer a perspective for future research in the identification and manipulation of metabolic trade-offs by relying on the toolbox provided by constraint-based metabolic modelling for single organisms and microbial communities.

COMMIT: Consideration of metabolite leakage and community composition improves microbial community reconstructions.

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.

GeneReg: a constraint-based approach for design of feasible metabolic engineering strategies at the gene level.

MOTIVATION: Large-scale metabolic models are widely used to design metabolic engineering strategies for diverse biotechnological applications. However, the existing computational approaches focus on alteration of reaction fluxes and often neglect the manipulations of gene expression to implement these strategies. RESULTS: Here, we find that the association of genes with multiple reactions leads to infeasibility of engineering strategies at the flux level, since they require contradicting manipulations of gene expression. Moreover, we identify that all of the existing approaches to design gene knockout strategies do not ensure that the resulting design may also require other gene alterations, such as up- or downregulations, to match the desired flux distribution. To address these issues, we propose a constraint-based approach, termed GeneReg, that facilitates the design of feasible metabolic engineering strategies at the gene level and that is readily applicable to large-scale metabolic networks. We show that GeneReg can identify feasible strategies to overproduce ethanol in Escherichia coli and lactate in Saccharomyces cerevisiae, but overproduction of the TCA cycle intermediates is not feasible in five organisms used as cell factories under default growth conditions. Therefore, GeneReg points at the need to couple gene regulation and metabolism to design rational metabolic engineering strategies. AVAILABILITY AND IMPLEMENTATION: https://github.com/MonaRazaghi/GeneReg. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PC2P: parameter-free network-based prediction of protein complexes.

MOTIVATION: Prediction of protein complexes from protein-protein interaction (PPI) networks is an important problem in systems biology, as they control different cellular functions. The existing solutions employ algorithms for network community detection that identify dense subgraphs in PPI networks. However, gold standards in yeast and human indicate that protein complexes can also induce sparse subgraphs, introducing further challenges in protein complex prediction. RESULTS: To address this issue, we formalize protein complexes as biclique spanned subgraphs, which include both sparse and dense subgraphs. We then cast the problem of protein complex prediction as a network partitioning into biclique spanned subgraphs with removal of minimum number of edges, called coherent partition. Since finding a coherent partition is a computationally intractable problem, we devise a parameter-free greedy approximation algorithm, termed Protein Complexes from Coherent Partition (PC2P), based on key properties of biclique spanned subgraphs. Through comparison with nine contenders, we demonstrate that PC2P: (i) successfully identifies modular structure in networks, as a prerequisite for protein complex prediction, (ii) outperforms the existing solutions with respect to a composite score of five performance measures on 75% and 100% of the analyzed PPI networks and gold standards in yeast and human, respectively, and (iii,iv) does not compromise GO semantic similarity and enrichment score of the predicted protein complexes. Therefore, our study demonstrates that clustering of networks in terms of biclique spanned subgraphs is a promising framework for detection of complexes in PPI networks. AVAILABILITY AND IMPLEMENTATION: https://github.com/SaraOmranian/PC2P. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.