Search

Secondary metabolites in Cyclocarya paliurus (Batalin) Iljinsk. leaves are beneficial for human health. The synthesis and accumulation of secondary metabolites form a complex process that is influenced by the trade-off between primary and secondary metabolism and by the biosynthetic pathways themselves. In this study, we explored the relationship between secondary metabolite accumulation and the activity of metabolic networks in leaves of C. paliurus. Leaves at three different growth stages were subjected to transcriptomic and non-targeted metabolomic analyses. The results revealed that nitrogen assimilation increased and carbon assimilation decreased as leaves matured, and the patterns of secondary metabolite accumulation and gene expression differed among the leaves at different growth stages. Mature green leaves had higher nitrogen assimilation and lower carbon assimilation, which were correlated with variations in secondary metabolite accumulation. As a major source of carbon and nitrogen, glutamine accumulated in the mature green leaves of C. paliurus. The accumulation of glutamine inhibited phenylalanine biosynthesis by modulating the pentose phosphate pathway but promoted acetyl-CoA biosynthesis through the tricarboxylic acid cycle. These changes led to decreased flavonoid contents and increased triterpenoid contents in mature leaves. These metabolomic and transcriptomic data reveal the differential expression of metabolic regulatory networks during three stages of leaf development and highlight the trade-off between primary and secondary metabolism. Our results provide a comprehensive picture of the metabolic pathways that are active in the leaves of C. paliurus at different growth stages.

The chemistry underpinning microbial interactions provides an integrative framework for linking the activities of individual microbes, microbial communities, plants, and their environments. Currently, we know very little about the functions of genes and metabolites within these communities because genome annotations and functions are derived from the minority of microbes that have been propagated in the laboratory. Yet the diversity, complexity, inaccessibility, and irreproducibility of native microbial consortia limit our ability to interpret chemical signaling and map metabolic networks. In this perspective, we contend that standardized laboratory ecosystems are needed to dissect the chemistry of soil microbiomes. We argue that dissemination and application of standardized laboratory ecosystems will be transformative for the field, much like how model organisms have played critical roles in advancing biochemistry and molecular and cellular biology. Community consensus on fabricated ecosystems ("EcoFABs") along with protocols and data standards will integrate efforts and enable rapid improvements in our understanding of the biochemical ecology of microbial communities.

Abstract The metabolic networks underlying skotomorphogenesis in seedlings remain relatively unknown. On the basis of our previous study on the folate metabolism in seedlings grown in darkness, the plastidial folylpolyglutamate synthetase gene (AtDFB) T‐DNA insertion Arabidopsis thaliana mutant (atdfb‐3) was examined. Under the nitrate‐sufficient condition, the mutant exhibited deficient folate metabolism and hypocotyl elongation, which affected skotomorphogenesis. Further analyses revealed changes to multiple intermediate metabolites related to carbon and nitrogen metabolism in the etiolated atdfb‐3 seedlings. Specifically, the sugar, polyol, and fatty acid contents decreased in the atdfb‐3 mutant under the nitrate‐sufficient condition, whereas the abundance of various organic acids and amino acids increased. In response to nitrate‐limited stress, multiple metabolites, including sugars, polyols, fatty acids, organic acids, and amino acids, accumulated more in the mutant than in the wild‐type control. The differences in the contents of multiple metabolites between the atdfb‐3 and wild‐type seedlings decreased following the addition of exogenous 5‐F‐THF under both nitrogen conditions. Additionally, the mutant accumulated high levels of one‐carbon metabolites, such as Cys, S‐adenosylmethionine, and S‐adenosylhomocysteine, under both nitrogen conditions. Thus, our data demonstrated that the perturbed folate metabolism in the atdfb‐3 seedlings, which was caused by the loss‐of‐function mutation to AtDFB, probably altered carbon and nitrogen metabolism, thereby modulating skotomorphogenesis. Furthermore, the study findings provide new evidence of the links among folate metabolism, metabolic networks, and skotomorphogenesis.

Background: Recent studies have provided insights into the persistence and succession of microbes aboard the International Space Station (ISS), notably the dominance of Klebsiella pneumoniae. However, the interactions between the various microbes aboard the ISS and how they shape the microbiome remain to be clearly understood. In this study, we apply a computational approach to predict possible metabolic interactions in the ISS microbiome and shed further light on its organization. Results: Through a combination of a systems-based graph-theoretical approach, and a constraint-based community metabolic modeling approach, we demonstrated several key interactions in the ISS microbiome. These complementary approaches provided insights into the metabolic interactions and dependencies present amongst various microbes in a community, highlighting key interactions and keystone species. Our results showed that the presence of K. pneumoniae is beneficial to many other microorganisms it coexists with, notably those from the Pantoea genus. Species belonging to the Enterobacteriaceae family were often found to be the most beneficial for the survival of other microorganisms in the ISS microbiome. However, K. pneumoniae was found to exhibit parasitic and amensalistic interactions with Aspergillus and Penicillium species, respectively. To prove this metabolic prediction, K. pneumoniae and Aspergillus fumigatus were co-cultured under normal and simulated microgravity, where K. pneumoniae cells showed parasitic characteristics to the fungus. The electron micrography revealed that the presence of K. pneumoniae compromised the morphology of fungal conidia and degenerated its biofilm-forming structures. Conclusion: Our study underscores the importance of K. pneumoniae in the ISS, and its potential positive and negative interactions with other microbes, including potential pathogens. This integrated modeling approach, combined with experiments, demonstrates the potential for understanding the organization of other such microbiomes, unravelling key organisms and their interdependencies. -HMRmrRFAwPDF1qD73q3qY Video Abstract [ABSTRACT FROM AUTHOR]

Two novel natural products, the polyketide cuniculene and the peptide antibiotic aquimarin, were recently discovered from the marine bacterial genus Aquimarina. However, the diversity of the secondary metabolite biosynthetic gene clusters (SM-BGCs) in Aquimarina genomes indicates a far greater biosynthetic potential. In this study, nine representative Aquimarina strains were tested for antimicrobial activity against diverse human-pathogenic and marine microorganisms and subjected to metabolomic and genomic profiling. We found an inhibitory activity of most Aquimarina strains against Candida glabrata and marine Vibrio and Alphaproteobacteria species. Aquimarina sp. Aq135 and Aquimarina muelleri crude extracts showed particularly promising antimicrobial activities, amongst others against methicillin-resistant Staphylococcus aureus. The metabolomic and functional genomic profiles of Aquimarina spp. followed similar patterns and were shaped by phylogeny. SM-BGC and metabolomics networks suggest the presence of novel polyketides and peptides, including cyclic depsipeptide-related compounds. Moreover, exploration of the 'Sponge Microbiome Project' dataset revealed that Aquimarina spp. possess low-abundance distributions worldwide across multiple marine biotopes. Our study emphasizes the relevance of this member of the microbial rare biosphere as a promising source of novel natural products. We predict that future metabologenomics studies of Aquimarina species will expand the spectrum of known secondary metabolites and bioactivities from marine ecosystems. [ABSTRACT FROM AUTHOR]

Genome-scale metabolism can best be described as a highly interconnected network of biochemical reactions and metabolites. The flow of metabolites, i.e., flux, throughout these networks can be predicted and analyzed using approaches such as flux balance analysis (FBA). By knowing the network topology and employing only a few simple assumptions, FBA can efficiently predict metabolic functions at the genome scale as well as microbial phenotypes. The network topology is represented in the form of genome-scale metabolic models (GEMs), which provide a direct mapping between network structure and function via the enzyme-coding genes and corresponding metabolic capacity. Recently, the role of protein limitations in shaping metabolic phenotypes have been extensively studied following the reconstruction of enzyme-constrained GEMs. This framework has been shown to significantly improve the accuracy of predicting microbial phenotypes, and it has demonstrated that a global limitation in protein availability can prompt the ubiquitous metabolic strategy of overflow metabolism. Being one of the most abundant and differentially expressed proteome sectors, metabolic proteins constitute a major cellular demand on proteinogenic amino acids. However, little is known about the impact and sensitivity of amino acid availability with regards to genome-scale metabolism. Here, we explore these aspects by extending on the enzyme-constrained GEM framework by also accounting for the usage of amino acids in expressing the metabolic proteome. Including amino acids in an enzyme-constrained GEM of Saccharomyces cerevisiae, we demonstrate that the expanded model is capable of accurately reproducing experimental amino acid levels. We further show that the metabolic proteome exerts variable demands on amino acid supplies in a condition-dependent manner, suggesting that S. cerevisiae must have evolved to efficiently fine-tune the synthesis of amino acids for expressing its metabolic proteins in response to changes in the external environment. Finally, our results demonstrate how the metabolic network of S. cerevisiae is robust towards perturbations of individual amino acids, while simultaneously being highly sensitive when the relative amino acid availability is set to mimic a priori distributions of both yeast and non-yeast origins.

Normal aging causes profound changes of structural degeneration and glucose hypometabolism in the human brain, even in the absence of disease. In recent years, with the extensive exploration of the topological characteristics of the human brain, related studies in rats have begun to investigate. However, age-related alterations of topological properties in individual brain metabolic network of rats remain unknown. In this study, a total of 48 healthy female Sprague–Dawley (SD) rats were used, including 24 young rats and 24 aged rats. We used Jensen-Shannon Divergence Similarity Estimation (JSSE) method for constructing individual metabolic networks to explore age-related topological properties and rich-club organization changes. Compared with the young rats, the aged rats showed significantly decreased clustering coefficient (Cp) and local efficiency (E loc ) across the whole-brain metabolic network. In terms of changes in local network measures, degree (D) and nodal efficiency (E nod ) of left posterior dorsal hippocampus, and E nod of left olfactory tubercle were higher in the aged rats than in the young rats. About the rich-club analysis, the existence of rich-club organization in individual brain metabolic networks of rats was demonstrated. In addition, our findings further confirmed that rich-club connections were susceptible to aging. Relative to the young rats, the overall strength of rich-club connections was significantly reduced in the aged rats, while the overall strength of feeder and local connections was significantly increased. These findings demonstrated the age-related reorganization principle of the brain structure and improved our understanding of brain alternations during aging. [ABSTRACT FROM AUTHOR]

Mathematical modeling of metabolic networks is a powerful approach to investigate the underlying principles of metabolism and growth. Such approaches include, among others, differential-equation-based modeling of metabolic systems, constraint-based modeling and metabolic network expansion of metabolic networks. Most of these methods are well established and are implemented in numerous software packages, but these are scattered between different programming languages, packages and syntaxes. This complicates establishing straight forward pipelines integrating model construction and simulation. We present a Python package moped that serves as an integrative hub for reproducible construction, modification, curation and analysis of metabolic models. moped supports draft reconstruction of models directly from genome/proteome sequences and pathway/genome databases utilizing GPR annotations, providing a completely reproducible model construction and curation process within executable Python scripts. Alternatively, existing models published in SBML format can be easily imported. Models are represented as Python objects, for which a wide spectrum of easy-to-use modification and analysis methods exist. The model structure can be manually altered by adding, removing or modifying reactions, and gap-filling reactions can be found and inspected. This greatly supports the development of draft models, as well as the curation and testing of models. Moreover, moped provides several analysis methods, in particular including the calculation of biosynthetic capacities using metabolic network expansion. The integration with other Python-based tools is facilitated through various model export options. For example, a model can be directly converted into a CobraPy object for constraint-based analyses. moped is a fully documented and expandable Python package. We demonstrate the capability to serve as a hub for integrating reproducible model construction and curation, database import, metabolic network expansion and export for constraint-based analyses. [ABSTRACT FROM AUTHOR]

Abstract Background Computational mathematical models of biological and biomedical systems have been successfully applied to advance our understanding of various regulatory processes, metabolic fluxes, effects of drug therapies, and disease evolution and transmission. Unfortunately, despite community efforts leading to the development of SBML and the BioModels database, many published models have not been fully exploited, largely due to a lack of proper documentation or the dependence on proprietary software. To facilitate the reuse and further development of systems biology and systems medicine models, an open-source toolbox that makes the overall process of model construction more consistent, understandable, transparent, and reproducible is desired. Results and discussion We provide an update on the development of modelbase, a free, expandable Python package for constructing and analysing ordinary differential equation-based mathematical models of dynamic systems. It provides intuitive and unified methods to construct and solve these systems. Significantly expanded visualisation methods allow for convenient analysis of the structural and dynamic properties of models. After specifying reaction stoichiometries and rate equations modelbase can automatically assemble the associated system of differential equations. A newly provided library of common kinetic rate laws reduces the repetitiveness of the computer programming code. modelbase is also fully compatible with SBML. Previous versions provided functions for the automatic construction of networks for isotope labelling studies. Now, using user-provided label maps, modelbase v1.2.3 streamlines the expansion of classic models to their isotope-specific versions. Finally, the library of previously published models implemented in modelbase is growing continuously. Ranging from photosynthesis to tumour cell growth to viral infection evolution, all these models are now available in a transparent, reusable and unified format through modelbase. Conclusion With this new Python software package, which is written in currently one of the most popular programming languages, the user can develop new models and actively profit from the work of others. modelbase enables reproducing and replicating models in a consistent, tractable and expandable manner. Moreover, the expansion of models to their isotopic label-specific versions enables simulating label propagation, thus providing quantitative information regarding network topology and metabolic fluxes.

Normal aging causes profound changes of structural degeneration and glucose hypometabolism in the human brain, even in the absence of disease. In recent years, with the extensive exploration of the topological characteristics of the human brain, related studies in rats have begun to investigate. However, age-related alterations of topological properties in individual brain metabolic network of rats remain unknown. In this study, a total of 48 healthy female Sprague–Dawley (SD) rats were used, including 24 young rats and 24 aged rats. We used Jensen-Shannon Divergence Similarity Estimation (JSSE) method for constructing individual metabolic networks to explore age-related topological properties and rich-club organization changes. Compared with the young rats, the aged rats showed significantly decreased clustering coefficient (Cp) and local efficiency (Eloc) across the whole-brain metabolic network. In terms of changes in local network measures, degree (D) and nodal efficiency (Enod) of left posterior dorsal hippocampus, and Enod of left olfactory tubercle were higher in the aged rats than in the young rats. About the rich-club analysis, the existence of rich-club organization in individual brain metabolic networks of rats was demonstrated. In addition, our findings further confirmed that rich-club connections were susceptible to aging. Relative to the young rats, the overall strength of rich-club connections was significantly reduced in the aged rats, while the overall strength of feeder and local connections was significantly increased. These findings demonstrated the age-related reorganization principle of the brain structure and improved our understanding of brain alternations during aging.

Abstract Genome-scale metabolic models of microorganisms are powerful frameworks to predict phenotypes from an organism’s genotype. While manual reconstructions are laborious, automated reconstructions often fail to recapitulate known metabolic processes. Here we present gapseq ( https://github.com/jotech/gapseq ), a new tool to predict metabolic pathways and automatically reconstruct microbial metabolic models using a curated reaction database and a novel gap-filling algorithm. On the basis of scientific literature and experimental data for 14,931 bacterial phenotypes, we demonstrate that gapseq outperforms state-of-the-art tools in predicting enzyme activity, carbon source utilisation, fermentation products, and metabolic interactions within microbial communities.

Abstract Background Significant efforts have been made in building large-scale kinetic models of cellular metabolism in the past two decades. However, most kinetic models published to date, remain focused around central carbon pathways or are built around ad hoc reduced models without clear justification on their derivation and usage. Systematic algorithms exist for reducing genome-scale metabolic reconstructions to build thermodynamically feasible and consistently reduced stoichiometric models. However, it is important to study how network complexity affects conclusions derived from large-scale kinetic models built around consistently reduced models before we can apply them to study biological systems. Results We reduced the iJO1366 Escherichia Coli genome-scale metabolic reconstruction systematically to build three stoichiometric models of different size. Since the reduced models are expansions around the core subsystems for which the reduction was performed, the models are nested. We present a method for scaling up the flux profile and the concentration vector reference steady-states from the smallest model to the larger ones, whilst preserving maximum equivalency. Populations of kinetic models, preserving similarity in kinetic parameters, were built around the reference steady-states and their metabolic sensitivity coefficients (MSCs) were computed. The MSCs were sensitive to the model complexity. We proposed a metric for measuring the sensitivity of MSCs to these structural changes. Conclusions We proposed for the first time a workflow for scaling up the size of kinetic models while preserving equivalency between the kinetic models. Using this workflow, we demonstrate that model complexity in terms of networks size has significant impact on sensitivity characteristics of kinetic models. Therefore, it is essential to account for the effects of network complexity when constructing kinetic models. The presented metric for measuring MSC sensitivity to structural changes can guide modelers and experimentalists in improving model quality and guide synthetic biology and metabolic engineering. Our proposed workflow enables the testing of the suitability of a kinetic model for answering certain study-specific questions. We argue that the model-based metabolic design targets that are common across models of different size are of higher confidence, while those that are different could be the objective of investigations for model improvement.

Human serine hydroxymethyltransferase (SHMT) regulates the serine-glycine one carbon metabolism and plays a role in cancer metabolic reprogramming. Two SHMT isozymes are acting in the cell: SHMT1 encoding the cytoplasmic isozyme, and SHMT2 encoding the mitochondrial one. Here we present a molecular model built on experimental data reporting the interaction between SHMT1 protein and SHMT2 mRNA, recently discovered in lung cancer cells. Using a stochastic dynamic model, we show that RNA moieties dynamically regulate serine and glycine concentration, shaping the system behaviour. For the first time we observe an active functional role of the RNA in the regulation of the serine-glycine metabolism and availability, which unravels a complex layer of regulation that cancer cells exploit to fine tune amino acids availability according to their metabolic needs. The quantitative model, complemented by an experimental validation in the lung adenocarcinoma cell line H1299, exploits RNA molecules as metabolic switches of the SHMT1 activity. Our results pave the way for the development of RNA-based molecules able to unbalance serine metabolism in cancer cells.

Mycobacterium tuberculosis, the causative agent of tuberculosis, is composed of several lineages characterized by a genome identity higher than 99%. Although the majority of the lineages are associated with humans, at least four lineages are adapted to other mammals, including different M. tuberculosis ecotypes. Host specificity is associated with higher virulence in its preferred host in ecotypes such as M. bovis. Deciphering what determines the preference of the host can reveal host-specific virulence patterns. However, it is not clear which genomic determinants might be influencing host specificity. In this study, we apply a combination of unsupervised and supervised classification methods on genomic data of ~27,000 M. tuberculosis clinical isolates to decipher host-specific genomic determinants. Host-specific genomic signatures are scarce beyond known lineage-specific mutations. Therefore, we integrated lineage-specific mutations into the iEK1011 2.0 genome-scale metabolic model to obtain lineage-specific versions of it. Flux distributions sampled from the solution spaces of these models can be accurately separated according to host association. This separation correlated with differences in cell wall processes, lipid, amino acid and carbon metabolic subsystems. These differences were observable when more than 95% of the samples had a specific growth rate significantly lower than the maximum achievable by the models. This suggests that these differences might manifest at low growth rate settings, such as the restrictive conditions M. tuberculosis suffers during macrophage infection. [ABSTRACT FROM AUTHOR]

The ever-increasing importance of structured data in different applications, especially in the biomedical field, has driven the need for reducing its complexity through projections into a more manageable space. The latest methods for learning features on graphs focus mainly on the neighborhood of nodes and edges. Methods capable of providing a representation that looks beyond the single node neighborhood are kernel graphs. However, they produce handcrafted features unaccustomed with a generalized model. To reduce this gap, in this work we propose a neural embedding framework, based on probability distribution representations of graphs, named Netpro2vec. The goal is to look at basic node descriptions other than the degree, such as those induced by the Transition Matrix and Node Distance Distribution. Netpro2vec provides embeddings completely independent from the task and nature of the data. The framework is evaluated on synthetic and various real biomedical network datasets through a comprehensive experimental classification phase and is compared to well-known competitors. [ABSTRACT FROM AUTHOR]

Motivation: The increasing availability of metabolomic data and their analysis are improving the understanding of cellular mechanisms and how biological systems respond to different perturbations. Currently, there is a need for novel computational methods that facilitate the analysis and integration of increasing volume of available data. Results: In this paper, we present Totoro a new constraint-based approach that integrates quantitative non-targeted metabolomic data of two different metabolic states into genome-wide metabolic models and predicts reactions that were most likely active during the transient state. We applied Totoro to real data of three different growth experiments (pulses of glucose, pyruvate, succinate) from Escherichia coli and we were able to predict known active pathways and gather new insights on the different metabolisms related to each substrate. We used both the E. coli core and the iJO1366 models to demonstrate that our approach is applicable to both smaller and larger networks. Availability: Totoro is an open source method (available at https://gitlab.inria.fr/erable/totoro) suitable for any organism with an available metabolic model. It is implemented in C++ and depends on IBM CPLEX which is freely available for academic purposes. [ABSTRACT FROM AUTHOR]

Subject of Research. The joint clustering method in network and correlation context is designed to identify active modules in metabolic graphs based on transcriptomic data represented by a large number of samples. The active modules obtained by this method describe the dynamic metabolic regulation in all samples of the analyzed dataset. The paper presents modifications of the proposed method for application on real data. Method. For results stability study the modified method was repeatedly run on real data with small variations of the initial parameters. For result analysis, several metrics were formulated that display modules similarity and representation under various start-up conditions. Main Results. The analysis results are sufficiently robust. For the most modules, their profiles are detected well in the noisy data, and the most genes are also preserved. Practical Relevance. The results of the presented study have shown that the modified method analyzes successfully real data by producing active modules that are stable and easy in interpretation.

Equilibrium bifurcations arise from sign changes of Jacobian determinants, as parameters are varied. Therefore we address the Jacobian determinant for metabolic networks with general reaction kinetics. Our approach is based on the concept of Child Selections: each (mother) metabolite is mapped, injectively, to one of those (child) reactions that it drives as an input. Our analysis distinguishes reaction network Jacobians with constant sign from the bifurcation case, where that sign depends on specific reaction rates. In particular, we distinguish "good" Child Selections, which do not affect the sign, from more interesting and mischievous "bad" children, which gang up towards sign changes, instability, and bifurcation.

Abstract Background The concept of minimal cut sets (MCS) has become an important mathematical framework for analyzing and (re)designing metabolic networks. However, the calculation of MCS in genome-scale metabolic models is a complex computational problem. The development of duality-based algorithms in the last years allowed the enumeration of thousands of MCS in genome-scale networks by solving mixed-integer linear problems (MILP). A recent advancement in this field was the introduction of the MCS2 approach. In contrast to the Farkas-lemma-based dual system used in earlier studies, the MCS2 approach employs a more condensed representation of the dual system based on the nullspace of the stoichiometric matrix, which, due to its reduced dimension, holds promise to further enhance MCS computations. Results In this work, we introduce several new variants and modifications of duality-based MCS algorithms and benchmark their effects on the overall performance. As one major result, we generalize the original MCS2 approach (which was limited to blocking the operation of certain target reactions) to the most general case of MCS computations with arbitrary target and desired regions. Building upon these developments, we introduce a new MILP variant which allows maximal flexibility in the formulation of MCS problems and fully leverages the reduced size of the nullspace-based dual system. With a comprehensive set of benchmarks, we show that the MILP with the nullspace-based dual system outperforms the MILP with the Farkas-lemma-based dual system speeding up MCS computation with an averaged factor of approximately 2.5. We furthermore present several simplifications in the formulation of constraints, mainly related to binary variables, which further enhance the performance of MCS-related MILP. However, the benchmarks also reveal that some highly condensed formulations of constraints, especially on reversible reactions, may lead to worse behavior when compared to variants with a larger number of (more explicit) constraints and involved variables. Conclusions Our results further enhance the algorithmic toolbox for MCS calculations and are of general importance for theoretical developments as well as for practical applications of the MCS framework.

Motivation: The increasing availability of metabolomic data and their analysis are improving the understanding of cellular mechanisms and how biological systems respond to different perturbations. Currently, there is a need for novel computational methods that facilitate the analysis and integration of increasing volume of available data.Results: In this paper, we present Totoro a new constraint-based approach that integrates quantitative non-targeted metabolomic data of two different metabolic states into genome-wide metabolic models and predicts reactions that were most likely active during the transient state. We applied Totoro to real data of three different growth experiments (pulses of glucose, pyruvate, succinate) from Escherichia coli and we were able to predict known active pathways and gather new insights on the different metabolisms related to each substrate. We used both the E. coli core and the iJO1366 models to demonstrate that our approach is applicable to both smaller and larger networks.Availability:Totoro is an open source method (available at https://gitlab.inria.fr/erable/totoro) suitable for any organism with an available metabolic model. It is implemented in C++ and depends on IBM CPLEX which is freely available for academic purposes.

Studying the correlation between microbiome metabolism and flavor of fermented foods has garnered significant attention recently. Understanding the contribution of metabolic regulation and environmental stress to microecosystems is essential for exploring the mechanisms of action of traditional fermented foods. Here, the interaction between microbial communities was investigated using a Xiaoqu fermentation system, processed as “simulative microecosystems,” in which starters were composed of Rhizopus-specific species/strains, Meyerozyma guilliermondii, and Bacillus licheniformis. The differences between community succession and metabolites were also explored. The results indicated that Rhizopus species/strain specificity affected starch hydrolyzation, resulting in a remarkable difference in the type and content of organic acids. This further suggested that the differences in nutrient abundance and organic acids influenced the colonization of microorganisms in the fermentation system, thereby influencing the succession of their communities. The fungi in the community predominantly originated from starters, whereas the bacteria were derived from both the environment and starter. Environmentally colonized microbes were the major contributors to the co-occurrence network and were strongly correlated with network. Regional characteristics of fermented foods were closely related to environmental microbes. These results contribute to the understanding of microbial assembly and flavor metabolism in fermented foods and provide strategies for quality regulation.

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs. [ABSTRACT FROM AUTHOR]

• Abstract interpretation allows gene knock out predictions based on elementary modes. • Qualitative analysis of metabolic networks realized through abstract interpretation. • Elementary modes provide Boolean analysis of reaction networks. Abstract interpretation was proposed for predicting changes of reaction networks with partial kinetic information in systems biology. This requires to compute the set of difference abstractions of a system of linear equations under nonlinear constraints. We present the first practical algorithm that can compute the difference abstractions of linear equation systems exactly. We also present a new heuristics based on minimal support consequences for overapproximating the set of difference abstractions. Our algorithms rely on elementary modes, first-order definitions, and finite domain constraint programming. We implemented our algorithms and applied them to change prediction in systems biology. It turns out experimentally that the new heuristics is often exact in practice, while outperforming the exact algorithm. This journal article extends on a paper published at the 17th International Conference on Computational Methods in Systems Biology (CMSB'2019) [1]. [ABSTRACT FROM AUTHOR]

Metabolism is the set of biochemical reactions of an organism that enables it to assimilate nutrients from its environment and to generate building blocks for growth and proliferation. It forms a complex network that is intertwined with the many molecular and cellular processes that take place within cells. Systems biology aims to capture the complexity of cells, organisms, or communities by reconstructing models based on information gathered by high-throughput analyses (omics data) and prior knowledge. One type of model is a genome-scale metabolic model (GEM) that allows studying the distributions of metabolic fluxes, i.e., the "mass-flow" through the network of biochemical reactions. GEMs are nowadays widely applied and have been reconstructed for various microbial pathogens, either in a free-living state or in interaction with their hosts, with the aim to gain insight into mechanisms of pathogenicity. In this review, we first introduce the principles of systems biology and GEMs. We then describe how metabolic modeling can contribute to unraveling microbial pathogenesis and host–pathogen interactions, with a specific focus on oomycete plant pathogens and in particular Phytophthora infestans. Subsequently, we review achievements obtained so far and identify and discuss potential pitfalls of current models. Finally, we propose a workflow for reconstructing high-quality GEMs and elaborate on the resources needed to advance a system biology approach aimed at untangling the intimate interactions between plants and pathogens. [ABSTRACT FROM AUTHOR]

The study of microbiomes by sequencing has revealed a plethora of correlations between microbial community composition and various life-history characteristics of the corresponding host species. However, inferring causation from correlation is often hampered by the sheer compositional complexity of microbiomes, even in simple organisms. Synthetic communities offer an effective approach to infer cause-effect relationships in host-microbiome systems. Yet the available communities suffer from several drawbacks, such as artificial (thus non-natural) choice of microbes, microbe-host mismatch (e.g., human microbes in gnotobiotic mice), or hosts lacking genetic tractability. Here we introduce CeMbio, a simplified natural Caenorhabditis elegans microbiota derived from our previous meta-analysis of the natural microbiome of this nematode. The CeMbio resource is amenable to all strengths of the C. elegans model system, strains included are readily culturable, they all colonize the worm gut individually, and comprise a robust community that distinctly affects nematode life-history. Several tools have additionally been developed for the CeMbio strains, including diagnostic PCR primers, completely sequenced genomes, and metabolic network models. With CeMbio, we provide a versatile resource and toolbox for the in-depth dissection of naturally relevant host-microbiome interactions in C. elegans.

Abstract Background Biological networks are representative of the diverse molecular interactions that occur within cells. Some of the commonly studied biological networks are modeled through protein-protein interactions, gene regulatory, and metabolic pathways. Among these, metabolic networks are probably the most studied, as they directly influence all physiological processes. Exploration of biochemical pathways using multigraph representation is important in understanding complex regulatory mechanisms. Feature extraction and clustering of these networks enable grouping of samples obtained from different biological specimens. Clustering techniques separate networks depending on their mutual similarity. Results We present a clustering analysis on tissue-specific metabolic networks for single samples from three primary tumor sites: breast, lung, and kidney cancer. The metabolic networks were obtained by integrating genome scale metabolic models with gene expression data. We performed network simplification to reduce the computational time needed for the computation of network distances. We empirically proved that networks clustering can characterize groups of patients in multiple conditions. Conclusions We provide a computational methodology to explore and characterize the metabolic landscape of tumors, thus providing a general methodology to integrate analytic metabolic models with gene expression data. This method represents a first attempt in clustering large scale metabolic networks. Moreover, this approach gives the possibility to get valuable information on what are the effects of different conditions on the overall metabolism.

Increased technological methods have enabled the investigation of biology at nanoscale levels. Such systems require the use of computational methods to comprehend the complex interactions that occur. The dynamics of metabolic systems have been traditionally described utilizing differential equations without fully capturing the heterogeneity of biological systems. Stochastic modeling approaches have recently emerged with the capacity to incorporate the statistical properties of such systems. However, the processing of stochastic algorithms is a computationally intensive task with intrinsic limitations. Alternatively, the queueing theory approach, historically used in the evaluation of telecommunication networks, can significantly reduce the computational power required to generate simulated results while simultaneously reducing the expansion of errors. We present here the application of queueing theory to simulate stochastic metabolic networks with high efficiency. With the use of glycolysis as a well understood biological model, we demonstrate the power of the proposed modeling methods discussed herein. Furthermore, we describe the simulation and pharmacological inhibition of glycolysis to provide an example of modeling capabilities.

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs.

Metabolism is the set of biochemical reactions of an organism that enables it to assimilate nutrients from its environment and to generate building blocks for growth and proliferation. It forms a complex network that is intertwined with the many molecular and cellular processes that take place within cells. Systems biology aims to capture the complexity of cells, organisms, or communities by reconstructing models based on information gathered by high-throughput analyses (omics data) and prior knowledge. One type of model is a genome-scale metabolic model (GEM) that allows studying the distributions of metabolic fluxes, i.e., the “mass-flow” through the network of biochemical reactions. GEMs are nowadays widely applied and have been reconstructed for various microbial pathogens, either in a free-living state or in interaction with their hosts, with the aim to gain insight into mechanisms of pathogenicity. In this review, we first introduce the principles of systems biology and GEMs. We then describe how metabolic modeling can contribute to unraveling microbial pathogenesis and host–pathogen interactions, with a specific focus on oomycete plant pathogens and in particular Phytophthora infestans. Subsequently, we review achievements obtained so far and identify and discuss potential pitfalls of current models. Finally, we propose a workflow for reconstructing high-quality GEMs and elaborate on the resources needed to advance a system biology approach aimed at untangling the intimate interactions between plants and pathogens.

Metabolic conversions allow organisms to produce essential metabolites from the available nutrients in an environment, frequently requiring metabolic exchanges among co-inhabiting organisms. Here, we applied genomic-based simulations for exploring tri-trophic interactions among the sap-feeding insect whitefly (Bemisia tabaci), its host-plants, and symbiotic bacteria. The simplicity of this ecosystem allows capturing the interacting organisms (based on genomic data) and the environmental content (based on metabolomics data). Simulations explored the metabolic capacities of insect-symbiont combinations under environments representing natural phloem. Predictions were correlated with experimental data on the dynamics of symbionts under different diets. Simulation outcomes depict a puzzle of three-layer origins (plant-insect-symbionts) for the source of essential metabolites across habitats and stratify interactions enabling the whitefly to feed on diverse hosts. In parallel to simulations, natural and artificial feeding experiments provide supporting evidence for an environment-based effect on symbiont dynamics. Based on simulations, a decrease in the relative abundance of a symbiont can be associated with a loss of fitness advantage due to an environmental excess in amino-acids whose production in a deprived environment used to depend on the symbiont. The study demonstrates that genomic-based predictions can bridge environment and community dynamics and guide the design of symbiont manipulation strategies. [ABSTRACT FROM AUTHOR]

Pregnancy generally alters the balance of maternal metabolism, but the molecular profiles in early pregnancy and associated factors of folate supplementation in pregnant women remains incompletely understood. Untargeted metabonomics based on high-performance liquid chromatography-high-resolution mass spectrometry integrated with multivariate metabolic pathway analysis were applied to characterize metabolite profiles and associated factors of folate supplements in early pregnancy. The metabolic baseline of early pregnancy was determined by metabolic analysis of 510 serum samples from 131 non-pregnant and 379 pregnant healthy Chinese women. The pathophysiology of adaptive reactions and metabolic challenges induced by folate supplementation in early pregnancy was further compared between pregnant women with (n = 168) and without (n = 184) folate supplements. Compared with non-pregnant participants, 106 metabolites, majority of which are related to amino acids and lysophosphatidylcholine/phosphatidylcholine, and 13 metabolic pathways were significantly changed in early pregnancy. The supplementation of folate in early pregnancy induced marked changes in N-acyl ethanolamine 22:0, N-acyl taurine 18:2, glycerophosphoserine 44:1 and 8,11,14-eicosatrienoate, proline, and aminoimidazole ribotide levels. During early pregnancy, the metabolism of amino acids significantly changes to meet the physiological requirements of pregnant women. Folate intake may change glucose and lipid metabolism. These findings provide a comprehensive landscape for understanding the basic characteristics and gestational metabolic networks of early pregnancy and folate supplementation. This study provides a basis for further research into the relationship between metabolic markers and pregnancy diseases. This study protocol was registered on www.ClinicalTrials.gov , NCT03651934, on August 29, 2018 (prior to recruitment). [ABSTRACT FROM AUTHOR]

The nutritional diversity of resources can affect the adaptive evolution of consumer metabolism and consumer diversification. The omega‐3 long‐chain polyunsaturated fatty acids eicosapentaenoic acid (EPA; 20:5n‐3) and docosahexaenoic acid (DHA; 22:6n‐3) have a high potential to affect consumer fitness, through their widespread effects on reproduction, growth and survival. However, few studies consider the evolution of fatty acid metabolism within an ecological context. In this review, we first document the extensive diversity in both primary producer and consumer fatty acid distributions amongst major ecosystems, between habitats and amongst species within habitats. We highlight some of the key nutritional contrasts that can shape behavioural and/or metabolic adaptation in consumers, discussing how consumers can evolve in response to the spatial, seasonal and community‐level variation of resource quality. We propose a hierarchical trait‐based approach for studying the evolution of consumers' metabolic networks and review the evolutionary genetic mechanisms underpinning consumer adaptation to EPA and DHA distributions. In doing so, we consider how the metabolic traits of consumers are hierarchically structured, from cell membrane function to maternal investment, and have strongly environment‐dependent expression. Finally, we conclude with an outlook on how studying the metabolic adaptation of consumers within the context of nutritional landscapes can open up new opportunities for understanding evolutionary diversification. [ABSTRACT FROM AUTHOR]

Model development and simulation of biological networks is recognized as a key task in Systems Biology. Integrated with in vitro and in vivo experimental data, network simulation allows for the discovery of the dynamics that regulate biological systems. Stochastic Petri Nets (SPNs) have become a widespread and reference formalism to model metabolic networks thanks to their natural expressiveness to represent metabolites, reactions, molecule interactions, and simulation randomness due to system fluctuations and environmental noise. In the literature, starting from the network model and the complete set of system parameters, there exist frameworks that allow for dynamic system simulation. Nevertheless, they do not allow for automatic model parameterization, which is a crucial task to identify, in silico, the network configurations that lead the model to satisfy specific temporal properties. To cover such a gap, this work first presents a framework to implement SPN models into SystemC code. Then, it shows how the framework allows for automatic parameterization of the networks. The user formally defines the network properties to be observed and the framework automatically extrapolates, through Assertion-based Verification (ABV), the parameter configurations that satisfy such properties. We present the results obtained by applying the proposed framework to model the complex metabolic network of the purine metabolism. We show how the automatic extrapolation of the system parameters allowed us to simulate the model under different conditions, which led to the understanding of behavioral differences in the regulation of the entire purine network. We also show the scalability of the approach through the modeling and simulation of four biological networks, each one with different structural characteristics. [ABSTRACT FROM AUTHOR]

URAT1 (slc22a12) was identified as the transporter responsible for renal reabsorption of the medically important compound, uric acid. However, subsequent studies have indicated that other transporters make contributions to this process, and that URAT1 transports other organic anions besides urate (including several in common with the closely related multi-specific renal organic anion transporters, OAT1 (slc22a6) and OAT3 (slc22a8)). These findings raise the possibility that urate transport is not the sole physiological function of URAT1. We previously characterized mice null for the murine ortholog of URAT1 (mURAT1; previously cloned as RST), finding a relatively modest decrement in urate reabsorptive capacity. Nevertheless, there were shifts in the plasma and urinary concentrations of multiple small molecules, suggesting significant metabolic changes in the knockouts. Although these molecules remain unidentified, here we have computationally delineated the biochemical networks consistent with transcriptomic data from the null mice. These analyses suggest alterations in the handling of not only urate but also other putative URAT1 substrates comprising intermediates in nucleotide, carbohydrate, and steroid metabolism. Moreover, the analyses indicate changes in multiple other pathways, including those relating to the metabolism of glycosaminoglycans, methionine, and coenzyme A, possibly reflecting downstream effects of URAT1 loss. Taken together with the available substrate and metabolomic data for the other OATs, our findings suggest that the transport and biochemical functions of URAT1 overlap those of OAT1 and OAT3, and could contribute to our understanding of the relationship between uric acid and the various metabolic disorders to which it has been linked.

High fructose consumption has been implicated with deleterious effects on human health, including hyperlipidemia elicited through de novo lipogenesis. However, more global effects of fructose on cellular metabolism have not been elucidated. In order to explore the metabolic impact of fructose-containing nutrients, we applied both GC-TOF and HILIC-QTOF mass spectrometry metabolomic strategies using extracts from cultured HepG2 cells exposed to fructose, glucose, or fructose + glucose. Cellular responses were analyzed in a time-dependent manner, incubated in media containing 5.5 mM glucose + 5.0 mM fructose in comparison to controls incubated in media containing either 5.5 mM glucose or 10.5 mM glucose. Mass spectrometry identified 156 unique known metabolites and a large number of unknown compounds, which revealed metabolite changes due to both utilization of fructose and high-carbohydrate loads independent of hexose structure. Fructose was shown to be partially converted to sorbitol, and generated higher levels of fructose-1-phosphate as a precursor for glycolytic intermediates. Differentially regulated ratios of 3-phosphoglycerate to serine pathway intermediates in high fructose media indicated a diversion of carbon backbones away from energy metabolism. Additionally, high fructose conditions changed levels of complex lipids toward phosphatidylethanolamines. Patterns of acylcarnitines in response to high hexose exposure (10.5 mM glucose or glucose/fructose combination) suggested a reduction in mitochondrial beta-oxidation.