The oral microbiome is associated with HPA axis response to a psychosocial stressor.

Intense psychosocial stress during early life has a detrimental effect on health-disease balance in later life. Simultaneously, despite its sensitivity to stress, the developing microbiome contributes to long-term health. Following stress exposure, HPA-axis activation regulates the "fight or flight" response with the release of glucose and cortisol. Here, we investigated the interaction between the oral microbiome and the stress response. We used a cohort of 115 adults, mean age 24, who either experienced institutionalisation and adoption (n = 40) or were non-adopted controls (n = 75). Glucose and cortisol measurements were taken from participants following an extended socially evaluated cold pressor test (seCPT) at multiple time points. The cohort´s oral microbiome was profiled via 16S-V4 sequencing on microbial DNA from saliva and buccal samples. Using mixed-effect linear regressions, we identified 12 genera that exhibited an interaction with host's cortisol-glucose response to stress, strongly influencing intensity and clearance of cortisol and glucose following stress exposure. Particularly, the identified taxa influenced the glucose and cortisol release profiles and kinetics following seCPT exposure. In conclusion, our study provided evidence for the oral microbiome modifying the effect of stress on the HPA-axis and human metabolism, as shown in glucose-cortisol time series data.

Microbial Abundances Retrieved from Sequencing data-automated NCBI Taxonomy (MARS): a pipeline to create relative microbial abundance data for the Microbiome Modelling Toolbox and utilizing homosynonyms for efficient mapping to resources.

Motivation: Computational approaches to the functional characterization of the microbiome, such as the Microbiome Modelling Toolbox, require precise information on microbial composition and relative abundances. However, challenges arise from homosynonyms-different names referring to the same taxon, which can hinder the mapping process and lead to missed species mapping when using microbial metabolic reconstruction resources, such as AGORA and APOLLO. Results: We introduce the integrated MARS pipeline, a user-friendly Python-based solution that addresses these challenges. MARS automates the extraction of relative abundances from metagenomic reads, maps species and genera onto microbial metabolic reconstructions, and accounts for alternative taxonomic names. It normalizes microbial reads, provides an optional cut-off for low-abundance taxa, and produces relative abundance tables apt for integration with the Microbiome Modelling Toolbox. A sub-component of the pipeline automates the task of identifying homosynonyms, leveraging web scraping to find taxonomic IDs of given species, searching NCBI for alternative names, and cross-reference them with microbial reconstruction resources. Taken together, MARS streamlines the entire process from processed metagenomic reads to relative abundance, thereby significantly reducing time and effort when working with microbiome data. Availability and implementation: MARS is implemented in Python. It can be found as an interactive application here: https://mars-pipeline.streamlit.app/along with a detailed documentation here: https://github.com/ThieleLab/mars-pipeline.

Personalized metabolic whole-body models for newborns and infants predict growth and biomarkers of inherited metabolic diseases.

Comprehensive whole-body models (WBMs) accounting for organ-specific dynamics have been developed to simulate adult metabolism, but such models do not exist for infants. Here, we present a resource of 360 organ-resolved, sex-specific models of newborn and infant metabolism (infant-WBMs) spanning the first 180 days of life. These infant-WBMs were parameterized to represent the distinct metabolic characteristics of newborns and infants, including nutrition, energy requirements, and thermoregulation. We demonstrate that the predicted infant growth was consistent with the recommendation by the World Health Organization. We assessed the infant-WBMs' reliability and capabilities for personalization by simulating 10,000 newborns based on their blood metabolome and birth weight. Furthermore, the infant-WBMs accurately predicted changes in known biomarkers over time and metabolic responses to treatment strategies for inherited metabolic diseases. The infant-WBM resource holds promise for personalized medicine, as the infant-WBMs could be a first step to digital metabolic twins for newborn and infant metabolism.

Whole-body metabolic modelling reveals microbiome and genomic interactions on reduced urine formate levels in Alzheimer's disease.

In this study, we aimed to understand the potential role of the gut microbiome in the development of Alzheimer's disease (AD). We took a multi-faceted approach to investigate this relationship. Urine metabolomics were examined in individuals with AD and controls, revealing decreased formate and fumarate concentrations in AD. Additionally, we utilised whole-genome sequencing (WGS) data obtained from a separate group of individuals with AD and controls. This information allowed us to create and investigate host-microbiome personalised whole-body metabolic models. Notably, AD individuals displayed diminished formate microbial secretion in these models. Additionally, we identified specific reactions responsible for the production of formate in the host, and interestingly, these reactions were linked to genes that have correlations with AD. This study suggests formate as a possible early AD marker and highlights genetic and microbiome contributions to its production. The reduced formate secretion and its genetic associations point to a complex connection between gut microbiota and AD. This holistic understanding might pave the way for novel diagnostic and therapeutic avenues in AD management.

Replication Protein A, the Main Eukaryotic Single-Stranded DNA Binding Protein, a Focal Point in Cellular DNA Metabolism.

Replication protein A (RPA) is a heterotrimeric protein complex and the main single-stranded DNA (ssDNA)-binding protein in eukaryotes. RPA has key functions in most of the DNA-associated metabolic pathways and DNA damage signalling. Its high affinity for ssDNA helps to stabilise ssDNA structures and protect the DNA sequence from nuclease attacks. RPA consists of multiple DNA-binding domains which are oligonucleotide/oligosaccharide-binding (OB)-folds that are responsible for DNA binding and interactions with proteins. These RPA-ssDNA and RPA-protein interactions are crucial for DNA replication, DNA repair, DNA damage signalling, and the conservation of the genetic information of cells. Proteins such as ATR use RPA to locate to regions of DNA damage for DNA damage signalling. The recruitment of nucleases and DNA exchange factors to sites of double-strand breaks are also an important RPA function to ensure effective DNA recombination to correct these DNA lesions. Due to its high affinity to ssDNA, RPA's removal from ssDNA is of central importance to allow these metabolic pathways to proceed, and processes to exchange RPA against downstream factors are established in all eukaryotes. These faceted and multi-layered functions of RPA as well as its role in a variety of human diseases will be discussed.

A toolbox of machine learning software to support microbiome analysis.

The human microbiome has become an area of intense research due to its potential impact on human health. However, the analysis and interpretation of this data have proven to be challenging due to its complexity and high dimensionality. Machine learning (ML) algorithms can process vast amounts of data to uncover informative patterns and relationships within the data, even with limited prior knowledge. Therefore, there has been a rapid growth in the development of software specifically designed for the analysis and interpretation of microbiome data using ML techniques. These software incorporate a wide range of ML algorithms for clustering, classification, regression, or feature selection, to identify microbial patterns and relationships within the data and generate predictive models. This rapid development with a constant need for new developments and integration of new features require efforts into compile, catalog and classify these tools to create infrastructures and services with easy, transparent, and trustable standards. Here we review the state-of-the-art for ML tools applied in human microbiome studies, performed as part of the COST Action ML4Microbiome activities. This scoping review focuses on ML based software and framework resources currently available for the analysis of microbiome data in humans. The aim is to support microbiologists and biomedical scientists to go deeper into specialized resources that integrate ML techniques and facilitate future benchmarking to create standards for the analysis of microbiome data. The software resources are organized based on the type of analysis they were developed for and the ML techniques they implement. A description of each software with examples of usage is provided including comments about pitfalls and lacks in the usage of software based on ML methods in relation to microbiome data that need to be considered by developers and users. This review represents an extensive compilation to date, offering valuable insights and guidance for researchers interested in leveraging ML approaches for microbiome analysis.

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

Causal inference on microbiome-metabolome relations in observational host-microbiome data via in silico in vivo association pattern analyses.

Understanding the effects of the microbiome on the host's metabolism is core to enlightening the role of the microbiome in health and disease. Herein, we develop the paradigm of in silico in vivo association pattern analyses, combining microbiome metabolome association studies with in silico constraint-based community modeling. Via theoretical dissection of confounding and causal paths, we show that in silico in vivo association pattern analyses allow for causal inference on microbiome-metabolome relations in observational data. We justify the corresponding theoretical criterion by structural equation modeling of host-microbiome systems, integrating deterministic microbiome community modeling into population statistics approaches. We show the feasibility of our approach on a published multi-omics dataset (n = 347), demonstrating causal microbiome-metabolite relations for 26 out of 54 fecal metabolites. In summary, we generate a promising approach for causal inference in metabolic host-microbiome interactions by integrating hypothesis-free screening association studies with knowledge-based in silico modeling.

Cardinality optimization in constraint-based modelling: application to human metabolism.

MOTIVATION: Several applications in constraint-based modelling can be mathematically formulated as cardinality optimization problems involving the minimization or maximization of the number of nonzeros in a vector. These problems include testing for stoichiometric consistency, testing for flux consistency, testing for thermodynamic flux consistency, computing sparse solutions to flux balance analysis problems and computing the minimum number of constraints to relax to render an infeasible flux balance analysis problem feasible. Such cardinality optimization problems are computationally complex, with no known polynomial time algorithms capable of returning an exact and globally optimal solution. RESULTS: By approximating the zero-norm with nonconvex continuous functions, we reformulate a set of cardinality optimization problems in constraint-based modelling into a difference of convex functions. We implemented and numerically tested novel algorithms that approximately solve the reformulated problems using a sequence of convex programs. We applied these algorithms to various biochemical networks and demonstrate that our algorithms match or outperform existing related approaches. In particular, we illustrate the efficiency and practical utility of our algorithms for cardinality optimization problems that arise when extracting a model ready for thermodynamic flux balance analysis given a human metabolic reconstruction. AVAILABILITY AND IMPLEMENTATION: Open source scripts to reproduce the results are here https://github.com/opencobra/COBRA.papers/2023_cardOpt with general purpose functions integrated within the COnstraint-Based Reconstruction and Analysis toolbox: https://github.com/opencobra/cobratoolbox.

Whole-body modelling reveals microbiome and genomic interactions on reduced urine formate levels in Alzheimer's disease.

In this study, we aimed to understand the potential role of the gut microbiome in the development of Alzheimer's disease (AD). We took a multi-faceted approach to investigate this relationship. Urine metabolomics were examined in individuals with AD and controls, revealing decreased formate and fumarate concentrations in AD. Additionally, we utilized whole-genome sequencing (WGS) data obtained from a separate group of individuals with AD and controls. This information allowed us to create and investigate host-microbiome personalized models. Notably, AD individuals displayed diminished formate microbial secretion in these models. Additionally, we identified specific reactions responsible for the production of formate in the host, and interestingly, these reactions were linked to genes that have correlations with AD. This study suggests formate as a possible early AD marker and highlights genetic and microbiome contributions to its production. The reduced formate secretion and its genetic associations point to a complex connection between gut microbiota and AD. This holistic understanding might pave the way for novel diagnostic and therapeutic avenues in AD management.

The Effects of Hospitalisation on the Serum Metabolome in COVID-19 Patients.

COVID-19, a systemic multi-organ disease resulting from infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is known to result in a wide array of disease outcomes, ranging from asymptomatic to fatal. Despite persistent progress, there is a continued need for more accurate determinants of disease outcomes, including post-acute symptoms after COVID-19. In this study, we characterised the serum metabolomic changes due to hospitalisation and COVID-19 disease progression by mapping the serum metabolomic trajectories of 71 newly hospitalised moderate and severe patients in their first week after hospitalisation. These 71 patients were spread out over three hospitals in Switzerland, enabling us to meta-analyse the metabolomic trajectories and filter consistently changing metabolites. Additionally, we investigated differential metabolite-metabolite trajectories between fatal, severe, and moderate disease outcomes to find prognostic markers of disease severity. We found drastic changes in serum metabolite concentrations for 448 out of the 901 metabolites. These results included markers of hospitalisation, such as environmental exposures, dietary changes, and altered drug administration, but also possible markers of physiological functioning, including carboxyethyl-GABA and fibrinopeptides, which might be prognostic for worsening lung injury. Possible markers of disease progression included altered urea cycle metabolites and metabolites of the tricarboxylic acid (TCA) cycle, indicating a SARS-CoV-2-induced reprogramming of the host metabolism. Glycerophosphorylcholine was identified as a potential marker of disease severity. Taken together, this study describes the metabolome-wide changes due to hospitalisation and COVID-19 disease progression. Moreover, we propose a wide range of novel potential biomarkers for monitoring COVID-19 disease course, both dependent and independent of the severity.

Longitudinal flux balance analyses of a patient with episodic colonic inflammation reveals microbiome metabolic dynamics.

We report the first use of constraint-based microbial community modeling on a single individual with episodic inflammation of the gastrointestinal tract, who has a well documented set of colonic inflammatory biomarkers, as well as metagenomically-sequenced fecal time series covering seven dates over 16 months. Between the first two time steps the individual was treated with both steroids and antibiotics. Our methodology enabled us to identify numerous time-correlated microbial species and metabolites. We found that the individual's dynamical microbial ecology in the disease state led to time-varying in silico overproduction, compared to healthy controls, of more than 24 biologically important metabolites, including methane, thiamine, formaldehyde, trimethylamine N-oxide, folic acid, serotonin, histamine, and tryptamine. The microbe-metabolite contribution analysis revealed that some Dialister species changed metabolic pathways according to the inflammation phases. At the first time point, characterized by the highest levels of serum (complex reactive protein) and fecal (calprotectin) inflammation biomarkers, they produced L-serine or formate. The production of the compounds, through a cascade effect, was mediated by the interaction with pathogenic Escherichia coli strains and Desulfovibrio piger. We integrated the microbial community metabolic models of each time point with a male whole-body, organ-resolved model of human metabolism to track the metabolic consequences of dysbiosis at different body sites. The presence of D. piger in the gut microbiome influenced the sulfur metabolism with a domino effect affecting the liver. These results revealed large longitudinal variations in an individual's gut microbiome ecology and metabolite production, potentially impacting other organs in the body. Future simulations with more time points from an individual could permit us to assess how external drivers, such as diet change or medical interventions, drive microbial community dynamics.

Circulating metabolites modulated by diet are associated with depression.

Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.

Systems biology elucidates the distinctive metabolic niche filled by the human gut microbe Eggerthella lenta.

Human gut bacteria perform diverse metabolic functions with consequences for host health. The prevalent and disease-linked Actinobacterium Eggerthella lenta performs several unusual chemical transformations, but it does not metabolize sugars and its core growth strategy remains unclear. To obtain a comprehensive view of the metabolic network of E. lenta, we generated several complementary resources: defined culture media, metabolomics profiles of strain isolates, and a curated genome-scale metabolic reconstruction. Stable isotope-resolved metabolomics revealed that E. lenta uses acetate as a key carbon source while catabolizing arginine to generate ATP, traits which could be recapitulated in silico by our updated metabolic model. We compared these in vitro findings with metabolite shifts observed in E. lenta-colonized gnotobiotic mice, identifying shared signatures across environments and highlighting catabolism of the host signaling metabolite agmatine as an alternative energy pathway. Together, our results elucidate a distinctive metabolic niche filled by E. lenta in the gut ecosystem. Our culture media formulations, atlas of metabolomics data, and genome-scale metabolic reconstructions form a freely available collection of resources to support further study of the biology of this prevalent gut bacterium.

A nutrition algorithm to optimize feed and medium composition using genome-scale metabolic models.

The optimization of animal feeds and cell culture media are problems of interest to a wide range of industries and scientific disciplines. Both problems are dictated by the properties of an organism's metabolism. However, due to the tremendous complexity of metabolic systems, it can be difficult to predict how metabolism will respond to changes in nutrient availability. A common tool used to capture the complexity of metabolism in a computational framework is a genome-scale metabolic model (GEM). GEMs are useful for predicting the fluxes of reactions within an organism's metabolism. To optimize feed or media, in silico experiments can be performed with GEMs by systematically varying nutritional constraints and predicting metabolic activity. In this way, the influence of various nutritional changes on metabolic outcomes can be evaluated. However, this methodology does not guarantee an optimal solution. Here, we develop a nutrition algorithm that utilizes linear programming to search the entire flux solution space of possible dietary intervention strategies to identify the most efficient changes to nutrition for a desirable metabolic outcome. We illustrate the utility of the nutrition algorithm on GEMs of Atlantic salmon (Salmo salar) and Chinese hamster ovary (CHO) cell metabolism and find that the nutrition algorithm makes predictions that not only align with experimental findings but reveal new insights into promising feeding strategies. We show that the nutrition algorithm is highly versatile and customizable to meet the user's needs. For instance, we demonstrate that the nutrition algorithm can be used to predict feed/media compositions that maximize profit margins. While the nutrition algorithm can be used to define an optimal feed/medium ab initio, it can also identify minimal changes to be made to an existing feed/medium to drive the largest metabolic shift. Moreover, the nutrition algorithm can target multiple metabolic pathways simultaneously with only a marginal increase in computational expense. While the nutrition algorithm has its limitations, we believe that this tool can be leveraged in a broad range of biotechnological applications to enhance the feed/medium optimization process.

Genome-scale metabolic reconstruction of 7,302 human microorganisms for personalized medicine.

The human microbiome influences the efficacy and safety of a wide variety of commonly prescribed drugs. Designing precision medicine approaches that incorporate microbial metabolism would require strain- and molecule-resolved, scalable computational modeling. Here, we extend our previous resource of genome-scale metabolic reconstructions of human gut microorganisms with a greatly expanded version. AGORA2 (assembly of gut organisms through reconstruction and analysis, version 2) accounts for 7,302 strains, includes strain-resolved drug degradation and biotransformation capabilities for 98 drugs, and was extensively curated based on comparative genomics and literature searches. The microbial reconstructions performed very well against three independently assembled experimental datasets with an accuracy of 0.72 to 0.84, surpassing other reconstruction resources and predicted known microbial drug transformations with an accuracy of 0.81. We demonstrate that AGORA2 enables personalized, strain-resolved modeling by predicting the drug conversion potential of the gut microbiomes from 616 patients with colorectal cancer and controls, which greatly varied between individuals and correlated with age, sex, body mass index and disease stages. AGORA2 serves as a knowledge base for the human microbiome and paves the way to personalized, predictive analysis of host-microbiome metabolic interactions.

Preterm birth is associated with xenobiotics and predicted by the vaginal metabolome.

Spontaneous preterm birth (sPTB) is a leading cause of maternal and neonatal morbidity and mortality, yet its prevention and early risk stratification are limited. Previous investigations have suggested that vaginal microbes and metabolites may be implicated in sPTB. Here we performed untargeted metabolomics on 232 second-trimester vaginal samples, 80 from pregnancies ending preterm. We find multiple associations between vaginal metabolites and subsequent preterm birth, and propose that several of these metabolites, including diethanolamine and ethyl glucoside, are exogenous. We observe associations between the metabolome and microbiome profiles previously obtained using 16S ribosomal RNA amplicon sequencing, including correlations between bacteria considered suboptimal, such as Gardnerella vaginalis, and metabolites enriched in term pregnancies, such as tyramine. We investigate these associations using metabolic models. We use machine learning models to predict sPTB risk from metabolite levels, weeks to months before birth, with good accuracy (area under receiver operating characteristic curve of 0.78). These models, which we validate using two external cohorts, are more accurate than microbiome-based and maternal covariates-based models (area under receiver operating characteristic curve of 0.55-0.59). Our results demonstrate the potential of vaginal metabolites as early biomarkers of sPTB and highlight exogenous exposures as potential risk factors for prematurity.

MetaboAnnotator: an efficient toolbox to annotate metabolites in genome-scale metabolic reconstructions.

MOTIVATION: Genome-scale metabolic reconstructions have been assembled for thousands of organisms using a wide range of tools. However, metabolite annotations, required to compare and link metabolites between reconstructions, remain incomplete. Here, we aim to further extend metabolite annotation coverage using various databases and chemoinformatic approaches. RESULTS: We developed a COBRA toolbox extension, deemed MetaboAnnotator, which facilitates the comprehensive annotation of metabolites with database independent and dependent identifiers, obtains molecular structure files, and calculates metabolite formula and charge at pH 7.2. The resulting metabolite annotations allow for subsequent cross-mapping between reconstructions and mapping of, e.g., metabolomic data. AVAILABILITY AND IMPLEMENTATION: MetaboAnnotator and tutorials are freely available at https://github.com/opencobra. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Whole-body metabolic modelling predicts isoleucine dependency of SARS-CoV-2 replication.

We aimed at investigating host-virus co-metabolism during SARS-CoV-2 infection. Therefore, we extended comprehensive sex-specific, whole-body organ resolved models of human metabolism with the necessary reactions to replicate SARS-CoV-2 in the lung as well as selected peripheral organs. Using this comprehensive host-virus model, we obtained the following key results: 1. The predicted maximal possible virus shedding rate was limited by isoleucine availability. 2. The supported initial viral load depended on the increase in CD4+ T-cells, consistent with the literature. 3. During viral infection, the whole-body metabolism changed including the blood metabolome, which agreed well with metabolomic studies from COVID-19 patients and healthy controls. 4. The virus shedding rate could be reduced by either inhibition of the guanylate kinase 1 or availability of amino acids, e.g., in the diet. 5. The virus variants differed in their maximal possible virus shedding rates, which could be inversely linked to isoleucine occurrences in the sequences. Taken together, this study presents the metabolic crosstalk between host and virus and emphasises the role of amino acid metabolism during SARS-CoV-2 infection, in particular of isoleucine. As such, it provides an example of how computational modelling can complement more canonical approaches to gain insight into host-virus crosstalk and to identify potential therapeutic strategies.