The RAVEN toolbox and its use for generating a genome-scale metabolic model for Penicillium chrysogenum

An industrial perspective on metabolic responses of Penicillium chrysogenum to periodic dissolved oxygen feast-famine cycles in a scale-down system

While traveling through different zones in large-scale bioreactors, microbes are most likely subjected to fluctuating dissolved oxygen (DO) conditions at the timescales of global circulation time. In this study, to mimic industrial-scale spatial DO gradients, we present a scale-down setup based on dynamic feast/famine regime (150 s) that leads to repetitive cycles with rapid changes in DO availability in glucose-limited chemostat cultures of Penicillium chrysogenum. Such DO feast/famine regime induced a stable and repetitive pattern with a reproducible metabolic response in time, and the dynamic response of intracellular metabolites featured specific differences in terms of both coverage and magnitude in comparison to other dynamic conditions, for example, substrate feast/famine cycles. Remarkably, intracellular sugar polyols were considerably increased as the hallmark metabolites along with a dynamic and higher redox state (NADH/NAD+) of the cytosol. Despite the increased availability of NADPH for penicillin production under the oscillatory DO conditions, this positive effect may be counteracted by the decreased ATP supply. Moreover, it is interesting to note that not only the penicillin productivity was reduced under such oscillating DO conditions, but also that of the unrecyclable byproduct ortho-hydroxyphenyl acetic acid and degeneration of penicillin productivity. Furthermore, dynamic flux profiles showed the most pronounced variations in central carbon metabolism, amino acid (AA) metabolism, energy metabolism and fatty acid metabolism upon the DO oscillation. Taken together, the metabolic responses of P. chrysogenum to DO gradients reported here are important for elucidating metabolic regulation mechanisms, improving bioreactor design and scale-up procedures as well as for constructing robust cell strains to cope with heterogenous industrial culture conditions.

Genome-scale community modelling elucidates the metabolic interaction in Indian type-2 diabetic gut microbiota

Type-2 diabetes (T2D) is a rapidly growing multifactorial metabolic disorder that induces the onset of various diseases in the human body. The compositional and metabolic shift of the gut microbiota is a crucial factor behind T2D. Hence, gaining insight into the metabolic profile of the gut microbiota is essential for revealing their role in regulating the metabolism of T2D patients. Here, we have focused on the genome-scale community metabolic model reconstruction of crucial T2D-associated gut microbes. The model-based analysis of biochemical flux in T2D and healthy gut conditions showed distinct biochemical signatures and diverse metabolic interactions in the microbial community. The metabolic interactions encompass cross-feeding of short-chain fatty acids, amino acids, and vitamins among individual microbes within the community. In T2D conditions, a reduction in the metabolic flux of acetate, butyrate, vitamin B5, and bicarbonate was observed in the microbial community model, which can impact carbohydrate metabolism. The decline in butyrate levels is correlated with both insulin resistance and diminished glucose metabolism in T2D patients. Compared to the healthy gut, an overall reduction in glucose consumption and SCFA production flux was estimated in the T2D gut environment. Moreover, the decreased consumption profiles of branch chain amino acids (BCAAs) and aromatic amino acids (AAAs) in the T2D gut microbiota can be a distinct biomarker for T2D. Hence, the flux-level analysis of the microbial community model can provide insights into the metabolic reprogramming in diabetic gut microbiomes, which may be helpful in personalized therapeutics and diet design against T2D.

Oomycete Metabolism Is Highly Dynamic and Reflects Lifestyle Adaptations

The selective pressure of pathogen-host symbiosis drives adaptations. How these interactions shape the metabolism of pathogens is largely unknown. Here, we use comparative genomics to systematically analyze the metabolic networks of oomycetes, a diverse group of eukaryotes that includes saprotrophs as well as animal and plant pathogens, with the latter causing devastating diseases with significant economic and/or ecological impacts. In our analyses of 44 oomycete species, we uncover considerable variation in metabolism that can be linked to lifestyle differences. Comparisons of metabolic gene content reveal that plant pathogenic oomycetes have a bipartite metabolism consisting of a conserved core and an accessory set. The accessory set can be associated with the degradation of defense compounds produced by plants when challenged by pathogens. Obligate biotrophic oomycetes have smaller metabolic networks, and taxonomically distantly related biotrophic lineages display convergent evolution by repeated gene losses in both the conserved as well as the accessory set of metabolisms. When investigating to what extent the metabolic networks in obligate biotrophs differ from those in hemibiotrophic plant pathogens, we observe that the losses of metabolic enzymes in obligate biotrophs are not random and that gene losses predominantly influence the terminal branches of the metabolic networks. Our analyses represent the first metabolism-focused comparison of oomycetes at this scale and will contribute to a better understanding of the evolution of oomycete metabolism in relation to lifestyle adaptation. Numerous oomycete species are devastating plant pathogens that cause major damage in crops and natural ecosystems. Their interactions with hosts are shaped by strong selection, but how selection affects adaptation of the primary metabolism to a pathogenic lifestyle is not yet well established. By pan-genome and metabolic network analyses of distantly related oomycete pathogens and their nonpathogenic relatives, we reveal considerable lifestyle- and lineage-specific adaptations. This study contributes to a better understanding of metabolic adaptations in pathogenic oomycetes in relation to lifestyle, host, and environment, and the findings will help in pinpointing potential targets for disease control. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Rewiring Saccharomyces cerevisiae metabolism for optimised Taxol® precursors production

Saccharomyces cerevisiae has been conveniently used to produce Taxol® anticancer drug early precursors. However, the harmful impact of oxidative stress by the first cytochrome P450-reductase enzymes (CYP725A4-POR) of Taxol® pathway has hampered sufficient progress in yeast. Here, we evolved an oxidative stress-resistant yeast strain with three-fold higher titre of their substrate, taxadiene. The performance of the evolved and parent strains were then evaluated in galactose-limited chemostats before and under the oxidative stress by an oxidising agent. The interaction of evolution and oxidative stress was comprehensively evaluated through transcriptomics and metabolite profiles integration in yeast enzyme-constrained genome scale model. Overall, the evolved strain showed improved respiration, reduced overflow metabolites production and oxidative stress re-induction tolerance. The cross-protection mechanism also potentially contributed to better heme, flavin and NADPH availability, essential for CYP725A4 and POR optimal activity in yeast. The results imply that the evolved strain is a robust cell factory for future efforts towards Taxol© production.

Development and applications of genome-scale metabolic network models

The genome-scale metabolic network model is an effective tool for characterizing the gene-protein-response relationship in the entire metabolic pathway of an organism. By combining various algorithms, the genome-scale metabolic network model can effectively simulate the influence of a specific environment on the physiological state of cells, optimize the culture conditions of strains, and predict the targets of genetic modification to achieve targeted modification of strains. In this review, we summarize the whole process of model building, sort out the various tools that may be involved in the model building process, and explain the role of various algorithms in model analysis. In addition, we also summarized the application of GSMM in network characteristics, cell phenotypes, metabolic engineering, etc. Finally, we discuss the current challenges facing GSMM.

Immune cell-specific and common molecular signatures in rheumatoid arthritis through molecular network approaches

Rheumatoid arthritis (RA) is an autoimmune disorder and common symptom of RA is chronic synovial inflammation. The pathogenesis of RA is not fully understood. Therefore, we aimed to identify underlying common and distinct molecular signatures and pathways among ten types of tissue and cells obtained from patients with RA. In this study, transcriptomic data including synovial tissues, macrophages, blood, T cells, CD4+T cells, CD8+T cells, natural killer T (NKT), cells natural killer (NK) cells, neutrophils, and monocyte cells were analyzed with an integrative and comparative network biology perspective. Each dataset yielded a list of differentially expressed genes as well as a reconstruction of the tissue-specific protein-protein interaction (PPI) network. Molecular signatures were identified by a statistical test using the hypergeometric probability density function by employing the interactions of transcriptional regulators and PPI. Reporter metabolites of each dataset were determined by using genome-scale metabolic networks. It was defined as the common hub proteins, novel molecular signatures, and metabolites in two or more tissue types while immune cell-specific molecular signatures were identified, too. Importantly, miR-155-5p is found as a common miRNA in all tissues. Moreover, NCOA3, PRKDC and miR-3160 might be novel molecular signatures for RA. Our results establish a novel approach for identifying immune cell-specific molecular signatures of RA and provide insights into the role of common tissue-specific genes, miRNAs, TFs, receptors, and reporter metabolites. Experimental research should be used to validate the corresponding genes, miRNAs, and metabolites.

Advances in genome-scale metabolic models of industrially important fungi

Many fungal species have been used industrially for production of biofuels and bioproducts. Developing strains with better performance in biomanufacturing contexts requires a systematic understanding of cellular metabolism. Genome-scale metabolic models (GEMs) offer a comprehensive view of interconnected pathways and a mathematical framework for downstream analysis. Recently, GEMs have been developed or updated for several industrially important fungi. Some of them incorporate enzyme constraints, enabling improved predictions of cell states and proteome allocation. Here, we provide an overview of these newly developed GEMs and computational methods that facilitate construction of enzyme-constrained GEMs and utilize flux predictions from GEMs. Furthermore, we highlight the pivotal roles of these GEMs in iterative design-build-test-learn cycles, ultimately advancing the field of fungal biomanufacturing.

APOLLO: A genome-scale metabolic reconstruction resource of 247,092 diverse human microbes spanning multiple continents, age groups, and body sites

Computational modelling of microbiome metabolism has proved instrumental to catalyse our understanding of diet-host-microbiome-disease interactions through the interrogation of mechanistic, strain- and molecule-resolved metabolic models. We present APOLLO, a resource of 247,092 human microbial genome-scale metabolic reconstructions spanning 19 phyla and accounting for microbial genomes from 34 countries, all age groups, and five body sites. We explored the metabolic potential of the reconstructed strains and developed a machine learning classifier able to predict with high accuracy the taxonomic strain assignments. We also built 14,451 sample-specific microbial community models, which could be stratified by body site, age, and disease states. Finally, we predicted faecal metabolites enriched or depleted in gut microbiomes of people with Crohn's disease, Parkinson disease, and undernourished children. APOLLO is compatible with the human whole-body models, and thus, provide unprecedented opportunities for systems-level modelling of personalised host-microbiome co-metabolism. APOLLO will be freely available under https://www.vmh.life/.

Genome-scale metabolic model analysis of Pichia pastoris for enhancing the production of S-adenosyl-L-methionine

Komagataella phaffii, formerly Pichia pastoris (P. pastoris), is a promising methylotrophic yeast used in industry to produce recombinant protein and valuable metabolites. In this study, a genome-scale metabolic model (GEMs) was reconstructed and used to assess P. pastoris' metabolic capabilities for the production of S-adenosyl-L-methionine (AdoMet or SAM or SAMe) from individual carbon sources along with the addition of L-methionine. In a model-driven P. pastoris strain, the well-established genome-scale metabolic model iAUKM can be implemented to predict high valuable metabolite production. The model, iAUKM, was created by merging the previously published iMT1026 model and the draught model generated using Raven toolbox from the KEGG database which covered 2309 enzymatic reactions associated with 1033 metabolic genes and 1750 metabolites. The highly curated model was successful in capturing P. pastoris growth on various carbon sources, as well as AdoMet production under various growth conditions. Many overexpression gene targets for increasing AdoMet accumulation in the cell have been predicted for various carbon sources. Inorganic phosphatase (IPP) was one of the predicted overexpression targets as revealed from simulations using iAUKM. When IPP gene was integrated into P. pastoris, we found that AdoMet accumulation increased by 16% and 14% using glucose and glycerol as carbon sources, respectively. Our in silico results shed light on the factors limiting AdoMet production, as well as key pathways for rationalized engineering to increase AdoMet yield.

Metabolic engineering of human gut microbiome: Recent developments and future perspectives

Many studies have demonstrated that the gut microbiota is associated with human health and disease. Manipulation of the gut microbiota, e.g. supplementation of probiotics, has been suggested to be feasible, but subject to limited therapeutic efficacy. To develop efficient microbiota-targeted diagnostic and therapeutic strategies, metabolic engineering has been applied to construct genetically modified probiotics and synthetic microbial consortia. This review mainly discusses commonly adopted strategies for metabolic engineering in the human gut microbiome, including the use of in silico, in vitro, or in vivo approaches for iterative design and construction of engineered probiotics or microbial consortia. Especially, we highlight how genome-scale metabolic models can be applied to advance our understanding of the gut microbiota. Also, we review the recent applications of metabolic engineering in gut microbiome studies as well as discuss important challenges and opportunities.

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

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

The genome-scale metabolic model for the purple non-sulfur bacterium Rhodopseudomonas palustris Bis A53 accurately predicts phenotypes under chemoheterotrophic, chemoautotrophic, photoheterotrophic, and photoautotrophic growth conditions

The purple non-sulfur bacterium Rhodopseudomonas palustris is recognized as a critical microorganism in the nitrogen and carbon cycle and one of the most common members in wastewater treatment communities. This bacterium is metabolically extremely versatile. It is capable of heterotrophic growth under aerobic and anaerobic conditions, but also able to grow photoautotrophically as well as mixotrophically. Therefore R. palustris can adapt to multiple environments and establish commensal relationships with other organisms, expressing various enzymes supporting degradation of amino acids, carbohydrates, nucleotides, and complex polymers. Moreover, R. palustris can degrade a wide range of pollutants under anaerobic conditions, e.g., aromatic compounds such as benzoate and caffeate, enabling it to thrive in chemically contaminated environments. However, many metabolic mechanisms employed by R. palustris to breakdown and assimilate different carbon and nitrogen sources under chemoheterotrophic or photoheterotrophic conditions remain unknown. Systems biology approaches, such as metabolic modeling, have been employed extensively to unravel complex mechanisms of metabolism. Previously, metabolic models have been reconstructed to study selected capabilities of R. palustris under limited experimental conditions. Here, we developed a comprehensive metabolic model (M-model) for R. palustris Bis A53 (iDT1294) consisting of 2,721 reactions, 2,123 metabolites, and comprising 1,294 genes. We validated the model using high-throughput phenotypic, physiological, and kinetic data, testing over 350 growth conditions. iDT1294 achieved a prediction accuracy of 90% for growth with various carbon and nitrogen sources and close to 80% for assimilation of aromatic compounds. Moreover, the M-model accurately predicts dynamic changes of growth and substrate consumption rates over time under nine chemoheterotrophic conditions and demonstrated high precision in predicting metabolic changes between photoheterotrophic and photoautotrophic conditions. This comprehensive M-model will help to elucidate metabolic processes associated with the assimilation of multiple carbon and nitrogen sources, anoxygenic photosynthesis, aromatic compound degradation, as well as production of molecular hydrogen and polyhydroxybutyrate.

TranSyT, an innovative framework for identifying transport systems

MOTIVATION

The importance and rate of development of genome-scale metabolic models have been growing for the last few years, increasing the demand for software solutions that automate several steps of this process. However, since TRIAGE's release, software development for the automatic integration of transport reactions into models has stalled.

RESULTS

Here, we present the Transport Systems Tracker (TranSyT). Unlike other transport systems annotation software, TranSyT does not rely on manual curation to expand its internal database, which is derived from highly curated records retrieved from the Transporters Classification Database and complemented with information from other data sources. TranSyT compiles information regarding transporter families and proteins, and derives reactions into its internal database, making it available for rapid annotation of complete genomes. All transport reactions have GPR associations and can be exported with identifiers from four different metabolite databases. TranSyT is currently available as a plugin for merlin v4.0 and an app for KBase.

AVAILABILITY AND IMPLEMENTATION

TranSyT web service: https://transyt.bio.di.uminho.pt/; GitHub for the tool: https://github.com/BioSystemsUM/transyt; GitHub with examples and instructions to run TranSyT: https://github.com/ecunha1996/transyt_paper.

Integrating Omics Data in Genome-Scale Metabolic Modeling: A Methodological Perspective for Precision Medicine

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease.

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

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

Unlocking the magic in mycelium: Using synthetic biology to optimize filamentous fungi for biomanufacturing and sustainability

Filamentous fungi drive carbon and nutrient cycling across our global ecosystems, through its interactions with growing and decaying flora and their constituent microbiomes. The remarkable metabolic diversity, secretion ability, and fiber-like mycelial structure that have evolved in filamentous fungi have been increasingly exploited in commercial operations. The industrial potential of mycelial fermentation ranges from the discovery and bioproduction of enzymes and bioactive compounds, the decarbonization of food and material production, to environmental remediation and enhanced agricultural production. Despite its fundamental impact in ecology and biotechnology, molds and mushrooms have not, to-date, significantly intersected with synthetic biology in ways comparable to other industrial cell factories (e.g. Escherichia coli,Saccharomyces cerevisiae, and Komagataella phaffii). In this review, we summarize a suite of synthetic biology and computational tools for the mining, engineering and optimization of filamentous fungi as a bioproduction chassis. A combination of methods across genetic engineering, mutagenesis, experimental evolution, and computational modeling can be used to address strain development bottlenecks in established and emerging industries. These include slow mycelium growth rate, low production yields, non-optimal growth in alternative feedstocks, and difficulties in downstream purification. In the scope of biomanufacturing, we then detail previous efforts in improving key bottlenecks by targeting protein processing and secretion pathways, hyphae morphogenesis, and transcriptional control. Bringing synthetic biology practices into the hidden world of molds and mushrooms will serve to expand the limited panel of host organisms that allow for commercially-feasible and environmentally-sustainable bioproduction of enzymes, chemicals, therapeutics, foods, and materials of the future.

Genome-scale community modeling for deciphering the inter-microbial metabolic interactions in fungus-farming termite gut microbiome

Specialized microbial communities in the fungus-farming termite gut and fungal comb microbiome help maintain host nutrition through interactive biochemical activities of complex carbohydrate degradation. Numerous research studies have been focused on identifying the microbial species in the termite gut and fungal comb microbiota, but the community-wide metabolic interaction patterns remain obscure. The inter-microbial metabolic interactions in the community environment are essential for executing biochemical processes like complex carbohydrate degradation and maintaining the host's physicochemical homeostasis. Recent progress in high-throughput sequencing techniques and mathematical modeling provides suitable platforms for constructing multispecies genome-scale community metabolic models that can render sound knowledge about microbial metabolic interaction patterns. Here, we have implemented the genome-scale metabolic modeling strategy to map the relationship between genes, proteins, and reactions of 12 key bacterial species from fungal cultivating termite gut and fungal comb microbiota. The resulting individual genome-scale metabolic models (GEMs) have been analyzed using flux balance analysis (FBA) to optimize the metabolic flux distribution pattern. Further, these individual GEMs have been integrated into genome-scale community metabolic models where a heuristics-based computational procedure has been employed to track the inter-microbial metabolic interactions. Two separate genome-scale community metabolic models were reconstructed for the O. badius gut and fungal comb microbiome. Analysis of the community models showed up to ∼167% increased flux range in lignocellulose degradation, amino acid biosynthesis, and nucleotide metabolism pathways. The inter-microbial metabolic exchange of amino acids, SCFAs, and small sugars was also upregulated in the multispecies community for maximum biomass formation. The flux variability analysis (FVA) has also been performed to calculate the feasible flux range of metabolic reactions. Furthermore, based on the calculated metabolic flux values, newly defined parameters, i.e., pairwise metabolic assistance (PMA) and community metabolic assistance (CMA) showed that the microbial species are getting up to 15% higher metabolic benefits in the multispecies community compared to pairwise growth. Assessment of the inter-microbial metabolic interaction patterns through pairwise growth support index (PGSI) indicated an increased mutualistic interaction in the termite gut environment compared to the fungal comb. Thus, this genome-scale community modeling study provides a systematic methodology to understand the inter-microbial interaction patterns with several newly defined parameters like PMA, CMA, and PGSI. The microbial metabolic assistance and interaction patterns derived from this computational approach will enhance the understanding of combinatorial microbial activities and may help develop effective synergistic microcosms to utilize complex plant polymers.

Combined metabolic activators improve metabolic functions in the animal models of neurodegenerative diseases

BACKGROUND

Neurodegenerative diseases (NDDs), including Alzheimer's disease (AD) and Parkinson's disease (PD), are associated with metabolic abnormalities. Integrative analysis of human clinical data and animal studies have contributed to a better understanding of the molecular and cellular pathways involved in the progression of NDDs. Previously, we have reported that the combined metabolic activators (CMA), which include the precursors of nicotinamide adenine dinucleotide and glutathione can be utilized to alleviate metabolic disorders by activating mitochondrial metabolism.

METHODS

We first analysed the brain transcriptomics data from AD patients and controls using a brain-specific genome-scale metabolic model (GEM). Then, we investigated the effect of CMA administration in animal models of AD and PD. We evaluated pathological and immunohistochemical findings of brain and liver tissues. Moreover, PD rats were tested for locomotor activity and apomorphine-induced rotation.

FINDINGS

Analysis of transcriptomics data with GEM revealed that mitochondrial dysfunction is involved in the underlying molecular pathways of AD. In animal models of AD and PD, we showed significant damage in the high-fat diet groups' brain and liver tissues compared to the chow diet. The histological analyses revealed that hyperemia, degeneration and necrosis in neurons were improved by CMA administration in both AD and PD animal models. These findings were supported by immunohistochemical evidence of decreased immunoreactivity in neurons. In parallel to the improvement in the brain, we also observed dramatic metabolic improvement in the liver tissue. CMA administration also showed a beneficial effect on behavioural functions in PD rats.

INTERPRETATION

Overall, we showed that CMA administration significantly improved behavioural scores in parallel with the neurohistological outcomes in the AD and PD animal models and is a promising treatment for improving the metabolic parameters and brain functions in NDDs.

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

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

Cellular Genome-Scale Metabolic Modeling Identifies New Potential Drug Targets Against Hepatocellular Carcinoma

Genome-scale metabolic modeling (GEM) is one of the key approaches to unpack cancer metabolism and for discovery of new drug targets. In this study, we report the Transcriptional Regulated Flux Balance Analysis-CORE (TRFBA-), an algorithm for GEM using key growth-correlated reactions using hepatocellular carcinoma (HCC), an important global health burden, as a case study. We generated a HepG2 cell-specific GEM by integrating this cell line transcriptomic data with a generic human metabolic model to forecast potential drug targets for HCC. A total of 108 essential genes for growth were predicted by the TRFBA-CORE. These genes were enriched for metabolic pathways involved in cholesterol, sterol, and steroid biosynthesis. Furthermore, we silenced a predicted essential gene, 11-beta dehydrogenase hydroxysteroid type 2 (HSD11B2), in HepG2 cells resulting in a reduction in cell viability. To further identify novel potential drug targets in HCC, we examined the effect of nine drugs targeting the essential genes, and observed that most drugs inhibited the growth of HepG2 cells. Some of these drugs in this model performed better than Sorafenib, the first-line therapeutic against HCC. A HepG2 cell-specific GEM highlights sterol metabolism to be essential for cell growth. HSD11B2 downregulation results in lower cell growth. Most of the compounds, selected by drug repurposing approach, show a significant inhibitory effect on cell growth in a wide range of concentrations. These findings offer new molecular leads for drug discovery for hepatic cancer while also illustrating the importance of GEM and drug repurposing in cancer therapeutics innovation.

ChiMera: an easy to use pipeline for bacterial genome based metabolic network reconstruction, evaluation and visualization

BACKGROUND

Genome-scale metabolic reconstruction tools have been developed in the last decades. They have helped to reconstruct eukaryotic and prokaryotic metabolic models, which have contributed to fields, e.g., genetic engineering, drug discovery, prediction of phenotypes, and other model-driven discoveries. However, the use of these programs requires a high level of bioinformatic skills. Moreover, the functionalities required to build models are scattered throughout multiple tools, requiring knowledge and experience for utilizing several tools.

RESULTS

Here we present ChiMera, which combines tools used for model reconstruction, prediction, and visualization. ChiMera uses CarveMe in the reconstruction module, generating a gap-filled draft reconstruction able to produce growth predictions using flux balance analysis for gram-positive and gram-negative bacteria. ChiMera also contains two modules for metabolic network visualization. The first module generates maps for the most important pathways, e.g., glycolysis, nucleotides and amino acids biosynthesis, fatty acid oxidation and biosynthesis and core-metabolism. The second module produces a genome-wide metabolic map, which can be used to retrieve KEGG pathway information for each compound in the model. A module to investigate gene essentiality and knockout is also present.

CONCLUSIONS

Overall, ChiMera uses automation algorithms to combine a variety of tools to automatically perform model creation, gap-filling, flux balance analysis (FBA), and metabolic network visualization. ChiMera models readily provide metabolic insights that can aid genetic engineering projects, prediction of phenotypes, and model-driven discoveries.

Genome-Scale Metabolic Modeling Reveals Sequential Dysregulation of Glutathione Metabolism in Livers from Patients with Alcoholic Hepatitis

Alcoholic hepatitis (AH) is the most severe form of alcoholic liver disease for which there is no efficacious treatment aiding most patients. AH manifests differently in individuals, with some patients showing debilitating symptoms more so than others. Previous studies showed significant metabolic dysregulation associated with AH. Therefore, we sought to analyze how the activity of metabolic pathways differed in the liver of patients with varying degrees of AH severity. We utilized a genome-scale metabolic modeling approach that allowed for integration of a generic human cellular metabolic model with specific RNA-seq data corresponding to healthy and multiple liver disease states to predict the metabolic fluxes within each disease state. Additionally, we performed a systems-level analysis of the transcriptomic data and predicted metabolic flux data to identify the regulatory and functional differences in liver metabolism with increasing severity of AH. Our results provide unique insights into the sequential dysregulation of the solute transport mechanisms underlying the glutathione metabolic pathway with increasing AH disease severity. We propose targeting of the solute transporters in the glutathione pathway to mimic the flux activity of the healthy liver state as a potential therapeutic intervention for AH.

Recent Development of Advanced Biotechnology in the Oleaginous Fungi for Arachidonic Acid Production

Arachidonic acid is an essential ω-6 polyunsaturated fatty acid, which plays a significant role in cardiovascular health and neurological development, leading to its wide use in the food and pharmaceutical industries. Traditionally, ARA is obtained from deep-sea fish oil. However, this source is limited by season and is depleting the already threatened global fish stocks. With the rapid development of synthetic biology in recent years, oleaginous fungi have gradually attracted increasing attention as promising microbial sources for large-scale ARA production. Numerous advanced technologies including metabolic engineering, dynamic regulation of fermentation conditions, and multiomics analysis were successfully adapted to increase ARA synthesis. This review summarizes recent advances in the bioengineering of oleaginous fungi for ARA production. Finally, perspectives for future engineering approaches are proposed to further improve the titer yield and productivity of ARA.

Integrative functional analysis uncovers metabolic differences between Candida species

Candida species are a dominant constituent of the human mycobiome and associated with the development of several diseases. Understanding the Candida species metabolism could provide key insights into their ability to cause pathogenesis. Here, we have developed the BioFung database, providing an efficient annotation of protein-encoding genes. Along, with BioFung, using carbohydrate-active enzyme (CAZymes) analysis, we have uncovered core and accessory features across Candida species demonstrating plasticity, adaption to the environment and acquired features. We show a greater importance of amino acid metabolism, as functional analysis revealed that all Candida species can employ amino acid metabolism. However, metabolomics revealed that only a specific cluster of species (AGAu species—C. albicans, C. glabrata and C. auris) utilised amino acid metabolism including arginine, cysteine, and methionine metabolism potentially improving their competitive fitness in pathogenesis. We further identified critical metabolic pathways in the AGAu cluster with biomarkers and anti-fungal target potential in the CAZyme profile, polyamine, choline and fatty acid biosynthesis pathways. This study, combining genomic analysis, and validation with gene expression and metabolomics, highlights the metabolic diversity with AGAu species that underlies their remarkable ability to dominate they mycobiome and cause disease.

Metabolic differences between Candida species are uncovered using the BioFung database alongside genomic and metabolic analysis.

Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

Metabolic modeling-based drug repurposing in Glioblastoma

The manifestation of intra- and inter-tumor heterogeneity hinders the development of ubiquitous cancer treatments, thus requiring a tailored therapy for each cancer type. Specifically, the reprogramming of cellular metabolism has been identified as a source of potential drug targets. Drug discovery is a long and resource-demanding process aiming at identifying and testing compounds early in the drug development pipeline. While drug repurposing efforts (i.e., inspecting readily available approved drugs) can be supported by a mechanistic rationale, strategies to further reduce and prioritize the list of potential candidates are still needed to facilitate feasible studies. Although a variety of ‘omics’ data are widely gathered, a standard integration method with modeling approaches is lacking. For instance, flux balance analysis is a metabolic modeling technique that mainly relies on the stoichiometry of the metabolic network. However, exploring the network’s topology typically neglects biologically relevant information. Here we introduce Transcriptomics-Informed Stoichiometric Modelling And Network analysis (TISMAN) in a recombinant innovation manner, allowing identification and validation of genes as targets for drug repurposing using glioblastoma as an exemplar.

Virtual Populations for Quantitative Systems Pharmacology Models

Quantitative systems pharmacology (QSP) places an emphasis on dynamic systems modeling, incorporating considerations from systems biology modeling and pharmacodynamics. The goal of QSP is often to quantitatively predict the effects of clinical therapeutics, their combinations, and their doses on clinical biomarkers and endpoints. In order to achieve this goal, strategies for incorporating clinical data into model calibration are critical. Virtual population (VPop) approaches facilitate model calibration while faced with challenges encountered in QSP model application, including modeling a breadth of clinical therapies, biomarkers, endpoints, utilizing data of varying structure and source, capturing observed clinical variability, and simulating with models that may require more substantial computational time and resources than often found in pharmacometrics applications. VPops are frequently developed in a process that may involve parameterization of isolated pathway models, integration into a larger QSP model, incorporation of clinical data, calibration, and quantitative validation that the model with the accompanying, calibrated VPop is suitable to address the intended question or help with the intended decision. Here, we introduce previous strategies for developing VPops in the context of a variety of therapeutic and safety areas: metabolic disorders, drug-induced liver injury, autoimmune diseases, and cancer. We introduce methodological considerations, prior work for sensitivity analysis and VPop algorithm design, and potential areas for future advancement. Finally, we give a more detailed application example of a VPop calibration algorithm that illustrates recent progress and many of the methodological considerations. In conclusion, although methodologies have varied, VPop strategies have been successfully applied to give valid clinical insights and predictions with the assistance of carefully defined and designed calibration and validation strategies. While a uniform VPop approach for all potential QSP applications may be challenging given the heterogeneity in use considerations, we anticipate continued innovation will help to drive VPop application for more challenging cases of greater scale while developing new rigorous methodologies and metrics.

Comparative analyses of parasites with a comprehensive database of genome-scale metabolic models

Protozoan parasites cause diverse diseases with large global impacts. Research on the pathogenesis and biology of these organisms is limited by economic and experimental constraints. Accordingly, studies of one parasite are frequently extrapolated to infer knowledge about another parasite, across and within genera. Model in vitro or in vivo systems are frequently used to enhance experimental manipulability, but these systems generally use species related to, yet distinct from, the clinically relevant causal pathogen. Characterization of functional differences among parasite species is confined to post hoc or single target studies, limiting the utility of this extrapolation approach. To address this challenge and to accelerate parasitology research broadly, we present a functional comparative analysis of 192 genomes, representing every high-quality, publicly-available protozoan parasite genome including Plasmodium, Toxoplasma, Cryptosporidium, Entamoeba, Trypanosoma, Leishmania, Giardia, and other species. We generated an automated metabolic network reconstruction pipeline optimized for eukaryotic organisms. These metabolic network reconstructions serve as biochemical knowledgebases for each parasite, enabling qualitative and quantitative comparisons of metabolic behavior across parasites. We identified putative differences in gene essentiality and pathway utilization to facilitate the comparison of experimental findings and discovered that phylogeny is not the sole predictor of metabolic similarity. This knowledgebase represents the largest collection of genome-scale metabolic models for both pathogens and eukaryotes; with this resource, we can predict species-specific functions, contextualize experimental results, and optimize selection of experimental systems for fastidious species.

Exploring the roles of microbes in facilitating plant adaptation to climate change

Plants benefit from their close association with soil microbes which assist in their response to abiotic and biotic stressors. Yet much of what we know about plant stress responses is based on studies where the microbial partners were uncontrolled and unknown. Under climate change, the soil microbial community will also be sensitive to and respond to abiotic and biotic stressors. Thus, facilitating plant adaptation to climate change will require a systems-based approach that accounts for the multi-dimensional nature of plant-microbe-environment interactions. In this perspective, we highlight some of the key factors influencing plant-microbe interactions under stress as well as new tools to facilitate the controlled study of their molecular complexity, such as fabricated ecosystems and synthetic communities. When paired with genomic and biochemical methods, these tools provide researchers with more precision, reproducibility, and manipulability for exploring plant-microbe-environment interactions under a changing climate.

Flux balance analysis of the ammonia-oxidizing bacterium Nitrosomonas europaea ATCC19718 unravels specific metabolic activities while degrading toxic compounds

The ammonia-oxidizing bacterium Nitrosomonas europaea has been widely recognized as an important player in the nitrogen cycle as well as one of the most abundant members in microbial communities for the treatment of industrial or sewage wastewater. Its natural metabolic versatility and extraordinary ability to degrade environmental pollutants (e.g., aromatic hydrocarbons such as benzene and toluene) enable it to thrive under various harsh environmental conditions. Constraint-based metabolic models constructed from genome sequences enable quantitative insight into the central and specialized metabolism within a target organism. These genome-scale models have been utilized to understand, optimize, and design new strategies for improved bioprocesses. Reduced modeling approaches have been used to elucidate Nitrosomonas europaea metabolism at a pathway level. However, genome-scale knowledge about the simultaneous oxidation of ammonia and pollutant metabolism of N. europaea remains limited. Here, we describe the reconstruction, manual curation, and validation of the genome-scale metabolic model for N. europaea, iGC535. This reconstruction is the most accurate metabolic model for a nitrifying organism to date, reaching an average prediction accuracy of over 90% under several growth conditions. The manually curated model can predict phenotypes under chemolithotrophic and chemolithoorganotrophic conditions while oxidating methane and wastewater pollutants. Calculated flux distributions under different trophic conditions show that several key pathways are affected by the type of carbon source available, including central carbon metabolism and energy production.

Nitrosomonas europaea catalyzes the first step of the nitrification process (ammonia to nitrite). It has been recognized as one of the most important members of microbial communities of wastewater treatment processes. Genome-scale models are powerful tools in process optimization since they can predict the organism’s behavior under different growth conditions. The final genome-scale model of N. europaea ATCC19718, iGC535, can predict growth and oxygen uptake rates with 90.52% accuracy under chemolithotrophic and chemolitoorganotrophic conditions. Moreover, iGC535 can predict the simultaneous oxidation of ammonia and wastewater pollutants, such as benzene, toluene, phenol and, chlorobenzene. iGC535 represents the most comprehensive knowledge-base for a nitrifying organism available to date. The genome-scale model reconstructed in this work brings us closer to understanding N. europaea’s role in a community with other nitrifying bacteria.

Changes in Oxygen Availability during Glucose-Limited Chemostat Cultivations of Penicillium chrysogenum Lead to Rapid Metabolite, Flux and Productivity Responses

Bioreactor scale-up from the laboratory scale to the industrial scale has always been a pivotal step in bioprocess development. However, the transition of a bioeconomy from innovation to commercialization is often hampered by performance loss in titer, rate and yield. These are often ascribed to temporal variations of substrate and dissolved oxygen (for instance) in the environment, experienced by microorganisms at the industrial scale. Oscillations in dissolved oxygen (DO) concentration are not uncommon. Furthermore, these fluctuations can be exacerbated with poor mixing and mass transfer limitations, especially in fermentations with filamentous fungus as the microbial cell factory. In this work, the response of glucose-limited chemostat cultures of an industrial Penicillium chrysogenum strain to different dissolved oxygen levels was assessed under both DO shift-down (60% → 20%, 10% and 5%) and DO ramp-down (60% → 0% in 24 h) conditions. Collectively, the results revealed that the penicillin productivity decreased as the DO level dropped down below 20%, while the byproducts, e.g., 6-oxopiperidine-2-carboxylic acid (OPC) and 6-aminopenicillanic acid (6APA), accumulated. Following DO ramp-down, penicillin productivity under DO shift-up experiments returned to its maximum value in 60 h when the DO was reset to 60%. The result showed that a higher cytosolic redox status, indicated by NADH/NAD+, was observed in the presence of insufficient oxygen supply. Consistent with this, flux balance analysis indicated that the flux through the glyoxylate shunt was increased by a factor of 50 at a DO value of 5% compared to the reference control, favoring the maintenance of redox status. Interestingly, it was observed that, in comparison with the reference control, the penicillin productivity was reduced by 25% at a DO value of 5% under steady state conditions. Only a 14% reduction in penicillin productivity was observed as the DO level was ramped down to 0. Furthermore, intracellular levels of amino acids were less sensitive to DO levels at DO shift-down relative to DO ramp-down conditions; this difference could be caused by different timescales between turnover rates of amino acid pools (tens of seconds to minutes) and DO switches (hours to days at steady state and minutes to hours at ramp-down). In summary, this study showed that changes in oxygen availability can lead to rapid metabolite, flux and productivity responses, and dynamic DO perturbations could provide insight into understanding of metabolic responses in large-scale bioreactors.

Host-mycobiome metabolic interactions in health and disease

Fungal communities (mycobiome) have an important role in sustaining the resilience of complex microbial communities and maintenance of homeostasis. The mycobiome remains relatively unexplored compared to the bacteriome despite increasing evidence highlighting their contribution to host-microbiome interactions in health and disease. Despite being a small proportion of the total species, fungi constitute a large proportion of the biomass within the human microbiome and thus serve as a potential target for metabolic reprogramming in pathogenesis and disease mechanism. Metabolites produced by fungi shape host niches, induce immune tolerance and changes in their levels prelude changes associated with metabolic diseases and cancer. Given the complexity of microbial interactions, studying the metabolic interplay of the mycobiome with both host and microbiome is a demanding but crucial task. However, genome-scale modelling and synthetic biology can provide an integrative platform that allows elucidation of the multifaceted interactions between mycobiome, microbiome and host. The inferences gained from understanding mycobiome interplay with other organisms can delineate the key role of the mycobiome in pathophysiology and reveal its role in human disease.

Metabolic systems biology and multi-omics of cyanobacteria: Perspectives and future directions

Cyanobacteria are oxygenic photoautotrophs whose metabolism contains key biochemical pathways to fix atmospheric CO₂ and synthesize various metabolites. The development of bioengineering tools has enabled the manipulation of cyanobacterial chassis to produce various valuable bioproducts photosynthetically. However, effective utilization of cyanobacteria as photosynthetic cell factories needs a detailed understanding of their metabolism and its interaction with other cellular processes. Implementing systems and synthetic biology tools has generated a wealth of information on various metabolic pathways. However, to design effective engineering strategies for further improvement in growth, photosynthetic efficiency, and enhanced production of target biochemicals, in-depth knowledge of their carbon/nitrogen metabolism, pathway fluxe distribution, genetic regulation and integrative analyses are necessary. In this review, we discuss the recent advances in the development of genome-scale metabolic models (GSMMs), omics analyses (metabolomics, transcriptomics, proteomics, fluxomics), and integrative modeling approaches to showcase the current understanding of cyanobacterial metabolism.

Modelling microbial communities: Harnessing consortia for biotechnological applications

Microbes propagate and thrive in complex communities, and there are many benefits to studying and engineering microbial communities instead of single strains. Microbial communities are being increasingly leveraged in biotechnological applications, as they present significant advantages such as the division of labour and improved substrate utilisation. Nevertheless, they also present some interesting challenges to surmount for the design of efficient biotechnological processes. In this review, we discuss key principles of microbial interactions, followed by a deep dive into genome-scale metabolic models, focussing on a vast repertoire of constraint-based modelling methods that enable us to characterise and understand the metabolic capabilities of microbial communities. Complementary approaches to model microbial communities, such as those based on graph theory, are also briefly discussed. Taken together, these methods provide rich insights into the interactions between microbes and how they influence microbial community productivity. We finally overview approaches that allow us to generate and test numerous synthetic community compositions, followed by tools and methodologies that can predict effective genetic interventions to further improve the productivity of communities. With impending advancements in high-throughput omics of microbial communities, the stage is set for the rapid expansion of microbial community engineering, with a significant impact on biotechnological processes.

Stratification of patients with clear cell renal cell carcinoma to facilitate drug repositioning

Clear cell renal cell carcinoma (ccRCC) is the most common histological type of kidney cancer and has high heterogeneity. Stratification of ccRCC is important since distinct subtypes differ in prognosis and treatment. Here, we applied a systems biology approach to stratify ccRCC into three molecular subtypes with different mRNA expression patterns and prognosis of patients. Further, we developed a set of biomarkers that could robustly classify the patients into each of the three subtypes and predict the prognosis of patients. Then, we reconstructed subtype-specific metabolic models and performed essential gene analysis to identify the potential drug targets. We identified four drug targets, including SOAT1, CRLS1, and ACACB, essential in all the three subtypes and GPD2, exclusively essential to subtype 1. Finally, we repositioned mitotane, an FDA-approved SOAT1 inhibitor, to treat ccRCC and showed that it decreased tumor cell viability and inhibited tumor cell growth based on in vitro experiments.

•

Three consistent molecular ccRCC subtypes were found to guide patients' prognoses

•

REOs-based biomarker was developed to robustly classify patients at individual level

•

SOAT1 is identified as a common drug target for all ccRCC subtypes

•

Mitotane was repositioned treatment of ccRCC via inhibiting SOAT1

Antibiotic resistance: Time of synthesis in a post-genomic age

Antibiotic resistance has been highlighted by international organizations, including World Health Organization, World Bank and United Nations, as one of the most relevant global health problems. Classical approaches to study this problem have focused in infected humans, mainly at hospitals. Nevertheless, antibiotic resistance can expand through different ecosystems and geographical allocations, hence constituting a One-Health, Global-Health problem, requiring specific integrative analytic tools. Antibiotic resistance evolution and transmission are multilayer, hierarchically organized processes with several elements (from genes to the whole microbiome) involved. However, their study has been traditionally gene-centric, each element independently studied. The development of robust-economically affordable whole genome sequencing approaches, as well as other -omic techniques as transcriptomics and proteomics, is changing this panorama. These technologies allow the description of a system, either a cell or a microbiome as a whole, overcoming the problems associated with gene-centric approaches. We are currently at the time of combining the information derived from -omic studies to have a more holistic view of the evolution and spread of antibiotic resistance. This synthesis process requires the accurate integration of -omic information into computational models that serve to analyse the causes and the consequences of acquiring AR, fed by curated databases capable of identifying the elements involved in the acquisition of resistance. In this review, we analyse the capacities and drawbacks of the tools that are currently in use for the global analysis of AR, aiming to identify the more useful targets for effective corrective interventions.

Reconstruction of a Genome-Scale Metabolic Model of Streptomyces albus J1074: Improved Engineering Strategies in Natural Product Synthesis

Streptomyces albus J1074 is recognized as an effective host for heterologous production of natural products. Its fast growth and efficient genetic toolbox due to a naturally minimized genome have contributed towards its advantage in expressing biosynthetic pathways for a diverse repertoire of products such as antibiotics and flavonoids. In order to develop precise model-driven engineering strategies for de novo production of natural products, a genome-scale metabolic model (GEM) was reconstructed for the microorganism based on protein homology to model species Streptomyces coelicolor while drawing annotated data from databases and literature for further curation. To demonstrate its capabilities, the Salb-GEM was used to predict overexpression targets for desirable compounds using flux scanning with enforced objective function (FSEOF). Salb-GEM was also utilized to investigate the effect of a minimized genome on metabolic gene essentialities in comparison to another Streptomyces species, S. coelicolor.

Genome-wide gene expression changes of Pseudomonas veronii 1YdBTEX2 during bioaugmentation in polluted soils

BACKGROUND

Bioaugmentation aims to use the capacities of specific bacterial strains inoculated into sites to enhance pollutant biodegradation. Bioaugmentation results have been mixed, which has been attributed to poor inoculant growth and survival in the field, and, consequently, moderate catalytic performance. However, our understanding of biodegradation activity mostly comes from experiments conducted under laboratory conditions, and the processes occurring during adaptation and invasion of inoculants into complex environmental microbiomes remain poorly known. The main aim of this work was thus to study the specific and different cellular reactions of an inoculant for bioaugmentation during adaptation, growth and survival in natural clean and contaminated non-sterile soils, in order to better understand factors limiting bioaugmentation.

RESULTS

As inoculant we focused on the monoaromatic compound-degrading bacterium Pseudomonas veronii 1YdBTEX2. The strain proliferated in all but one soil types in presence and in absence of exogenously added toluene. RNAseq and differential genome-wide gene expression analysis illustrated both a range of common soil responses such as increased nutrient scavenging and recycling, expression of defense mechanisms, as well as environment-specific reactions, notably osmoprotection and metal homeostasis. The core metabolism of P. veronii remained remarkably constant during exponential growth irrespective of the environment, with slight changes in cofactor regeneration pathways, possibly needed for balancing defense reactions.

CONCLUSIONS

P. veronii displayed a versatile global program, enabling it to adapt to a variety of soil environments in the presence and even in absence of its target pollutant toluene. Our results thus challenge the widely perceived dogma of poor survival and growth of exogenous inoculants in complex microbial ecosystems such as soil and provide a further basis to developing successful bioaugmentation strategies.

Atlas: automatic modeling of regulation of bacterial gene expression and metabolism using rule-based languages

MOTIVATION

Cells are complex systems composed of hundreds of genes whose products interact to produce elaborated behaviors. To control such behaviors, cells rely on transcription factors to regulate gene expression, and gene regulatory networks (GRNs) are employed to describe and understand such behavior. However, GRNs are static models, and dynamic models are difficult to obtain due to their size, complexity, stochastic dynamics and interactions with other cell processes.

RESULTS

We developed Atlas, a Python software that converts genome graphs and gene regulatory, interaction and metabolic networks into dynamic models. The software employs these biological networks to write rule-based models for the PySB framework. The underlying method is a divide-and-conquer strategy to obtain sub-models and combine them later into an ensemble model. To exemplify the utility of Atlas, we used networks of varying size and complexity of Escherichia coli and evaluated in silico modifications, such as gene knockouts and the insertion of promoters and terminators. Moreover, the methodology could be applied to the dynamic modeling of natural and synthetic networks of any bacteria.

AVAILABILITY AND IMPLEMENTATION

Code, models and tutorials are available online (https://github.com/networkbiolab/atlas).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

PhenoMapping: a protocol to map cellular phenotypes to metabolic bottlenecks, identify conditional essentiality, and curate metabolic models

Targeted identification of cellular processes responsible for a phenotype is of major importance in guiding efforts in bioengineering and medicine. Genome-scale metabolic models (GEMs) are widely used to integrate various types of omics data and study the cellular physiology under different conditions. Here, we present PhenoMapping, a protocol that uses GEMs, omics, and phenotypic data to map cellular processes and observed phenotypes. PhenoMapping also classifies genes as conditionally and unconditionally essential and guides a comprehensive curation of GEMs. For complete details on the use and execution of this protocol, please refer to Stanway et al. (2019) and Krishnan et al. (2020).

Deep learning approaches for natural product discovery from plant endophytic microbiomes

Plant microbiomes are not only diverse, but also appear to host a vast pool of secondary metabolites holding great promise for bioactive natural products and drug discovery. Yet, most microbes within plants appear to be uncultivable, and for those that can be cultivated, their metabolic potential lies largely hidden through regulatory silencing of biosynthetic genes. The recent explosion of powerful interdisciplinary approaches, including multi-omics methods to address multi-trophic interactions and artificial intelligence-based computational approaches to infer distribution of function, together present a paradigm shift in high-throughput approaches to natural product discovery from plant-associated microbes. Arguably, the key to characterizing and harnessing this biochemical capacity depends on a novel, systematic approach to characterize the triggers that turn on secondary metabolite biosynthesis through molecular or genetic signals from the host plant, members of the rich 'in planta' community, or from the environment. This review explores breakthrough approaches for natural product discovery from plant microbiomes, emphasizing the promise of deep learning as a tool for endophyte bioprospecting, endophyte biochemical novelty prediction, and endophyte regulatory control. It concludes with a proposed pipeline to harness global databases (genomic, metabolomic, regulomic, and chemical) to uncover and unsilence desirable natural products.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s40793-021-00375-0.

High-quality genome-scale metabolic model of Aurantiochytrium sp. T66

The long-chain, ω-3 polyunsaturated fatty acids (PUFAs) (e.g., eicosapentaenoic acid [EPA] and docosahexaenoic acid [DHA]), are essential for humans and animals, including marine fish species. Presently, the primary source of these PUFAs is fish oils. As the global production of fish oils appears to be reaching its limits, alternative sources of high-quality ω-3 PUFAs is paramount to support the growing aquaculture industry. Thraustochytrids are a group of heterotrophic protists with the capability to synthesize and accrue large amounts of DHA. Thus, the thraustochytrids are prime candidates to solve the increasing demand for ω-3 PUFAs using microbial cell factories. However, a systems-level understanding of their metabolic shift from cellular growth into lipid accumulation is, to a large extent, unclear. Here, we reconstructed a high-quality genome-scale metabolic model of the thraustochytrid Aurantiochytrium sp. T66 termed iVS1191. Through iterative rounds of model refinement and extensive manual curation, we significantly enhanced the metabolic scope and coverage of the reconstruction from that of previously published models, making considerable improvements with stoichiometric consistency, metabolic connectivity, and model annotations. We show that iVS1191 is highly consistent with experimental growth data, reproducing in vivo growth phenotypes as well as specific growth rates on minimal carbon media. The availability of iVS1191 provides a solid framework for further developing our understanding of T66's metabolic properties, as well as exploring metabolic engineering and process-optimization strategies in silico for increased ω-3 PUFA production.

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

Recent advances in systems and synthetic biology approaches for developing novel cell-factories in non-conventional yeasts

Microbial bioproduction of chemicals, proteins, and primary metabolites from cheap carbon sources is currently an advancing area in industrial research. The model yeast, Saccharomyces cerevisiae, is a well-established biorefinery host that has been used extensively for commercial manufacturing of bioethanol from myriad carbon sources. However, its Crabtree-positive nature often limits the use of this organism for the biosynthesis of commercial molecules that do not belong in the fermentative pathway. To avoid extensive strain engineering of S. cerevisiae for the production of metabolites other than ethanol, non-conventional yeasts can be selected as hosts based on their natural capacity to produce desired commodity chemicals. Non-conventional yeasts like Kluyveromyces marxianus, K. lactis, Yarrowia lipolytica, Pichia pastoris, Scheffersomyces stipitis, Hansenula polymorpha, and Rhodotorula toruloides have been considered as potential industrial eukaryotic hosts owing to their desirable phenotypes such as thermotolerance, assimilation of a wide range of carbon sources, as well as ability to secrete high titers of protein and lipid. However, the advanced metabolic engineering efforts in these organisms are still lacking due to the limited availability of systems and synthetic biology methods like in silico models, well-characterised genetic parts, and optimized genome engineering tools. This review provides an insight into the recent advances and challenges of systems and synthetic biology as well as metabolic engineering endeavours towards the commercial usage of non-conventional yeasts. Particularly, the approaches in emerging non-conventional yeasts for the production of enzymes, therapeutic proteins, lipids, and metabolites for commercial applications are extensively discussed here. Various attempts to address current limitations in designing novel cell factories have been highlighted that include the advances in the fields of genome-scale metabolic model reconstruction, flux balance analysis, 'omics'-data integration into models, genome-editing toolkit development, and rewiring of cellular metabolisms for desired chemical production. Additionally, the understanding of metabolic networks using ¹³C-labelling experiments as well as the utilization of metabolomics in deciphering intracellular fluxes and reactions have also been discussed here. Application of cutting-edge nuclease-based genome editing platforms like CRISPR/Cas9, and its optimization towards efficient strain engineering in non-conventional yeasts have also been described. Additionally, the impact of the advances in promising non-conventional yeasts for efficient commercial molecule synthesis has been meticulously reviewed. In the future, a cohesive approach involving systems and synthetic biology will help in widening the horizon of the use of unexplored non-conventional yeast species towards industrial biotechnology.

Linking Gut Microbiome and Lipid Metabolism: Moving beyond Associations

Various studies aiming to elucidate the role of the gut microbiome-metabolome co-axis in health and disease have primarily focused on water-soluble polar metabolites, whilst non-polar microbial lipids have received less attention. The concept of microbiota-dependent lipid biotransformation is over a century old. However, only recently, several studies have shown how microbial lipids alter intestinal and circulating lipid concentrations in the host, thus impacting human lipid homeostasis. There is emerging evidence that gut microbial communities play a particularly significant role in the regulation of host cholesterol and sphingolipid homeostasis. Here, we review and discuss recent research focusing on microbe-host-lipid co-metabolism. We also discuss the interplay of human gut microbiota and molecular lipids entering host systemic circulation, and its role in health and disease.