BacArena: Individual-based metabolic modeling of heterogeneous microbes in complex communities

MetaBiome: a multiscale model integrating agent-based and metabolic networks to reveal spatial regulation in gut mucosal microbial communities

Mucosal microbial communities (MMCs) are complex ecosystems near the mucosal layers of the gut essential for maintaining health and modulating disease states. Despite advances in high-throughput omics technologies, current methodologies struggle to capture the dynamic metabolic interactions and spatiotemporal variations within MMCs. In this work, we present MetaBiome, a multiscale model integrating agent-based modeling (ABM), finite volume methods, and constraint-based models to explore the metabolic interactions within these communities. Integrating ABM allows for the detailed representation of individual microbial agents each governed by rules that dictate cell growth, division, and interactions with their surroundings. Through a layered approach-encompassing microenvironmental conditions, agent information, and metabolic pathways-we simulated different communities to showcase the potential of the model. Using our in-silico platform, we explored the dynamics and spatiotemporal patterns of MMCs in the proximal small intestine and the cecum, simulating the physiological conditions of the two gut regions. Our findings revealed how specific microbes adapt their metabolic processes based on substrate availability and local environmental conditions, shedding light on spatial metabolite regulation and informing targeted therapies for localized gut diseases. MetaBiome provides a detailed representation of microbial agents and their interactions, surpassing the limitations of traditional grid-based systems. This work marks a significant advancement in microbial ecology, as it offers new insights into predicting and analyzing microbial communities.IMPORTANCEOur study presents a novel multiscale model that combines agent-based modeling, finite volume methods, and genome-scale metabolic models to simulate the complex dynamics of mucosal microbial communities in the gut. This integrated approach allows us to capture spatial and temporal variations in microbial interactions and metabolism that are difficult to study experimentally. Key findings from our model include the following: (i) prediction of metabolic cross-feeding and spatial organization in multi-species communities, (ii) insights into how oxygen gradients and nutrient availability shape community composition in different gut regions, and (iii) identification of spatiallyregulated metabolic pathways and enzymes in E. coli. We believe this work represents a significant advance in computational modeling of microbial communities and provides new insights into the spatial regulation of gut microbiome metabolism. The multiscale modeling approach we have developed could be broadly applicable for studying other complex microbial ecosystems.

Integrating network and in-silico simulation insights into the ecological interactions shaped by carbon sources in partial denitrification and anammox system

The underlying ecological mechanism of microbial communities shaped by carbon source in partial denitrification and anammox (PDA) systems remains poorly understood, despite the potential of multiple carbon sources to support the partial denitrification process. Herein, the integrated network and in-silico simulation methods were used to evaluate the considerable impact of carbon sources on the dynamics of ecological interactions. The fluctuation of carbon source (from acetate to glucose and ethanol) significantly destabilized the performance of PDA system (total nitrogen removal efficiency decreased from 96.8% to 69.1%). Glucose simultaneously altered the composition of denitrifying bacteria, resulting in a significant enrichment of the genus Elstera (from 0% to 12.7%). By contrast, genus Thauera re-dominated for partial denitrification with ethanol as carbon source. Importantly, heterotrophic bacteria (e.g., genus Calditrichia) gradually enriched by utilizing ethanol. The presence of acetate in phase IV further enhanced the competitive advantage of heterotrophic bacteria over denitrifying bacteria, thereby resulting in the deteriorated performance of the PDA system. The in-silico simulation of co-culture further revealed that the overgrowth of auxotrophic species Calditrichia utilized amounts of nutrients and limited other functional bacteria. Additionally, the whole co-occurrence network indicated that positive interactions likely improved the adaptability of anammox bacteria under the unsteady conditions. This study provides profound insights into the ecological interactions shaped by carbon sources in PDA systems and underscores the necessity of comprehensive review of the external carbon source to ensure optimal performance.

Agent-based model of the human colon to investigate mechanisms of pathogen colonization resistance

Recent global burden of disease studies have shown that bacterial infections are responsible for over 13 million deaths worldwide, or 1 in every 8 deaths, each year. Enteric diarrheal infections, in particular, pose a significant challenge and strain on healthcare systems as many are difficult to address pharmaceutically, and thus rely primarily on the patient's own immune system and gut microbiome to fight the infection. Nonetheless, the specific mechanisms behind gut microbiome colonization resistance of enteric pathogens are not well defined and microbiome diversity is difficult to represent and study experimentally. To address this gap, we have constructed an agent-based computational model of the colonic epithelium cross section to investigate the colonic invasion of enteric pathogens. The model focuses on three main regions: epithelial layer, mucosal bilayer, and adjacent lumen, and utilizes four main cell types as agents: anaerobic bacteria, facultative anaerobic bacteria, human goblet cells, and pathogens. Utilizing this model, we are able to describe the healthy microbiome cell localization and dynamics from our mucosal bilayer. In addition, we are also able to investigate the impact of host dietary fiber consumption and simulate pathogen invasion. The model exemplifies the possibility and potential to explore key gut microbiome colonization resistance mechanisms and environmental impacts on the gut microbiome using computational methods.

Moving from genome-scale to community-scale metabolic models for the human gut microbiome

Metabolic models of individual microorganisms or small microbial consortia have become standard research tools in the bioengineering and systems biology fields. However, extending metabolic modelling to diverse microbial communities, such as those in the human gut, remains a practical challenge from both modelling and experimental validation perspectives. In complex communities, metabolic models accounting for community dynamics, or those that consider multiple objectives, may provide optimal predictions over simpler steady-state models, but require a much higher computational cost. Here we describe some of the strengths and limitations of microbial community-scale metabolic models and argue for a robust validation framework for developing personalized, mechanistic and accurate predictions of microbial community metabolic behaviours across environmental contexts. Ultimately, quantitatively accurate microbial community-scale metabolic models could aid in the design and testing of personalized prebiotic, probiotic and dietary interventions that optimize for translationally relevant outcomes.

Control of medical digital twins with artificial neural networks

The objective of precision medicine is to tailor interventions to an individual patient's unique characteristics. A key technology for this purpose involves medical digital twins, computational models of human biology that can be personalized and dynamically updated to incorporate patient-specific data. Certain aspects of human biology, such as the immune system, are not easily captured with physics-based models, such as differential equations. Instead, they are often multi-scale, stochastic and hybrid. This poses a challenge to existing control and optimization approaches that cannot be readily applied to such models. Recent advances in neural-network control methods hold promise in addressing complex control problems. However, the application of these approaches to biomedical systems is still in its early stages. This work employs dynamics-informed neural-network controllers as an alternative approach to control of medical digital twins. As a first use case, we focus on the control of agent-based models (ABMs), a versatile and increasingly common modelling platform in biomedicine. The effectiveness of the proposed neural-network control methods is illustrated and benchmarked against other methods with two widely used ABMs. To account for the inherent stochastic nature of the ABMs we aim to control, we quantify uncertainty in relevant model and control parameters.This article is part of the theme issue 'Uncertainty quantification for healthcare and biological systems (Part 1)'.

Harnessing microbial division of labor for biomanufacturing: a review of laboratory and formal modeling approaches

Bioprocess industries aim to meet the increasing demand for product complexity by designing enhanced cellular and metabolic capabilities for the host. Monocultures, standard biomanufacturing workhorses, are often restricted in their capability to meet these demands, and the solution often involves the genetic modification of the host. Synthetic microbial communities are a promising alternative to monocultures because they exhibit division of labor, enabling efficient resource utilization and pathway modularity. This specialization minimizes metabolic burden and enhances robustness to perturbations, providing a competitive advantage. Despite this potential, their utilization in biotechnological or bioprocessing applications remains limited. The recent emergence of new and innovative community design tools and strategies, particularly those harnessing the division of labor, holds promise to change this outlook. Understanding the microbial interactions governing natural microbial communities can be used to identify complementary partners, informing synthetic community design. Therefore, we particularly consider engineering division of labor in synthetic microbial communities as a viable solution to accelerate progress in the field. This review presents the current understanding of how microbial interactions enable division of labor and how this information can be used to design synthetic microbial communities to perform tasks otherwise unfeasible to individual organisms. We then evaluate laboratory and formal modeling approaches specifically developed to: elucidate microbial community physiology, guide experimental design, and improve our understanding of complex community interactions assisting synthetic community design. By synthesizing these insights, we aim to present a comprehensive framework that advances the use of microbial communities in biomanufacturing applications.

Coupling flux balance analysis with reactive transport modeling through machine learning for rapid and stable simulation of microbial metabolic switching

Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.

Genome-scale models in human metabologenomics

Metabologenomics integrates metabolomics with other omics data types to comprehensively study the genetic and environmental factors that influence metabolism. These multi-omics data can be incorporated into genome-scale metabolic models (GEMs), which are highly curated knowledge bases that explicitly account for genes, transcripts, proteins and metabolites. By including all known biochemical reactions catalysed by enzymes and transporters encoded in the human genome, GEMs analyse and predict the behaviour of complex metabolic networks. Continued advancements to the scale and scope of GEMs - from cells and tissues to microbiomes and the whole body - have helped to design effective treatments and develop better diagnostic tools for metabolic diseases. Furthermore, increasing amounts of multi-omics data are incorporated into GEMs to better identify the underlying mechanisms, biomarkers and potential drug targets of metabolic diseases.

Advances in functional analysis of the microbiome: Integrating metabolic modeling, metabolite prediction, and pathway inference with Next-Generation Sequencing data

This review explores current advancements in microbiome functional analysis enabled by next-generation sequencing technologies, which have transformed our understanding of microbial communities from mere taxonomic composition to their functional potential. We examine approaches that move beyond species identification to characterize microbial activities, interactions, and their roles in host health and disease. Genome-scale metabolic models allow for in-depth simulations of metabolic networks, enabling researchers to predict microbial metabolism, growth, and interspecies interactions in diverse environments. Additionally, computational methods for predicting metabolite profiles offer indirect insights into microbial metabolic outputs, which is crucial for identifying biomarkers and potential therapeutic targets. Functional pathway analysis tools further reveal microbial contributions to metabolic pathways, highlighting alterations in response to environmental changes and disease states. Together, these methods offer a powerful framework for understanding the complex metabolic interactions within microbial communities and their impact on host physiology. While significant progress has been made, challenges remain in the accuracy of predictive models and the completeness of reference databases, which limit the applicability of these methods in under-characterized ecosystems. The integration of these computational tools with multi-omic data holds promise for personalized approaches in precision medicine, allowing for targeted interventions that modulate the microbiome to improve health outcomes. This review highlights recent advances in microbiome functional analysis, providing a roadmap for future research and translational applications in human health and environmental microbiology.

Metabolic Flux Modeling in Marine Ecosystems

Ocean metabolism constitutes a complex, multiscale ensemble of biochemical reaction networks harbored within and between the boundaries of a myriad of organisms. Gaining a quantitative understanding of how these networks operate requires mathematical tools capable of solving in silico the resource allocation problem each cell faces in real life. Toward this goal, stoichiometric modeling of metabolism, such as flux balance analysis, has emerged as a powerful computational tool for unraveling the intricacies of metabolic processes in microbes, microbial communities, and multicellular organisms. Here, we provide an overview of this approach and its applications, future prospects, and practical considerations in the context of marine sciences. We explore how flux balance analysis has been employed to study marine organisms, help elucidate nutrient cycling, and predict metabolic capabilities within diverse marine environments, and highlight future prospects for this field in advancing our knowledge of marine ecosystems and their sustainability.

Genome-scale metabolic modelling reveals interactions and key roles of symbiont clades in a sponge holobiont

Sponges harbour complex microbiomes and as ancient metazoans and important ecosystem players are emerging as powerful models to understand the evolution and ecology of symbiotic interactions. Metagenomic studies have previously described the functional features of sponge symbionts, however, little is known about the metabolic interactions and processes that occur under different environmental conditions. To address this issue, we construct here constraint-based, genome-scale metabolic networks for the microbiome of the sponge Stylissa sp. Our models define the importance of sponge-derived nutrients for microbiome stability and discover how different organic inputs can result in net heterotrophy or autotrophy of the symbiont community. The analysis further reveals the key role that a newly discovered bacterial taxon has in cross-feeding activities and how it dynamically adjusts with nutrient inputs. Our study reveals insights into the functioning of a sponge microbiome and provides a framework to further explore and define metabolic interactions in holobionts.

Effects of burning on vegetation, soil physicochemistry and prokaryotic microbial communities in surface and subsurface peat

Prescribed burning is a common management strategy in peatlands that has the potential to affect soil physicochemistry, alter biogeochemical cycles and trigger changes in vegetation structure. How burning affects prokaryotic community composition across different soil profiles is not well understood. This study explored the effects of prescribed burning on the diversity of prokaryotic communities in peat soils. Soil samples were collected from Moor House Nature Reserve, UK, a long-term monitoring site initiated in 1954 subject to three burning treatments: Burning at short rotations every 10 years, burning at long rotations every 20 years and a non-burn control. Observed species richness for archaea was highest in the topsoil of the non-burn control plots and highest for bacteria in the topsoil of the non-burn control and plots under a long rotation regime. Community composition was significantly different between different burn treatments and soil depth. Archaeal community structure was shaped by NH4+ and pH in the topsoil; by Pb, moisture and Al in the 20-40 cm profile; and by total N, total C, Al, Ca, Fe and pH in the 40-60 cm profile. Bacterial community structure was shaped by NH4+, heather cover, pH and Mg in the topsoil; by Fe, K and Pb in the 20-40 cm profile; and by Al, Ca and Fe in the 40-60 cm profile. A co-occurrence network analysis revealed that the topsoil of the non-burn control plots had a more complex network structure with more positive links than those under a rotational burn, but a higher average connectivity with a higher number of negative links was observed in the long rotation 20-40 cm profile. The results provide a new insight into the response processes of soil prokaryotic communities to burning in peatland soils, providing valuable knowledge that can support the evaluation and management of ecosystem services in peatlands.

Review on computer-assisted biosynthetic capacities elucidation to assess metabolic interactions and communication within microbial communities

Microbial communities thrive through interactions and communication, which are challenging to study as most microorganisms are not cultivable. To address this challenge, researchers focus on the extracellular space where communication events occur. Exometabolomics and interactome analysis provide insights into the molecules involved in communication and the dynamics of their interactions. Advances in sequencing technologies and computational methods enable the reconstruction of taxonomic and functional profiles of microbial communities using high-throughput multi-omics data. Network-based approaches, including community flux balance analysis, aim to model molecular interactions within and between communities. Despite these advances, challenges remain in computer-assisted biosynthetic capacities elucidation, requiring continued innovation and collaboration among diverse scientists. This review provides insights into the current state and future directions of computer-assisted biosynthetic capacities elucidation in studying microbial communities.

Bacterial community dynamics as a result of growth-yield trade-off and multispecies metabolic interactions toward understanding the gut biofilm niche

Bacterial communities are ubiquitous, found in natural ecosystems, such as soil, and within living organisms, like the human microbiome. The dynamics of these communities in diverse environments depend on factors such as spatial features of the microbial niche, biochemical kinetics, and interactions among bacteria. Moreover, in many systems, bacterial communities are influenced by multiple physical mechanisms, such as mass transport and detachment forces. One example is gut mucosal communities, where dense, closely packed communities develop under the concurrent influence of nutrient transport from the lumen and fluid-mediated detachment of bacteria. In this study, we model a mucosal niche through a coupled agent-based and finite-volume modeling approach. This methodology enables us to model bacterial interactions affected by nutrient release from various sources while adjusting individual bacterial kinetics. We explored how the dispersion and abundance of bacteria are influenced by biochemical kinetics in different types of metabolic interactions, with a particular focus on the trade-off between growth rate and yield. Our findings demonstrate that in competitive scenarios, higher growth rates result in a larger share of the niche space. In contrast, growth yield plays a critical role in neutralism, commensalism, and mutualism interactions. When bacteria are introduced sequentially, they cause distinct spatiotemporal effects, such as deeper niche colonization in commensalism and mutualism scenarios driven by species intermixing effects, which are enhanced by high growth yields. Moreover, sub-ecosystem interactions dictate the dynamics of three-species communities, sometimes yielding unexpected outcomes. Competitive, fast-growing bacteria demonstrate robust colonization abilities, yet they face challenges in displacing established mutualistic systems. Bacteria that develop a cooperative relationship with existing species typically obtain niche residence, regardless of their growth rates, although higher growth yields significantly enhance their abundance. Our results underscore the importance of bacterial niche dynamics in shaping community properties and succession, highlighting a new approach to manipulating microbial systems.

The marine environmental microbiome mediates physiological outcomes in host nematodes

BACKGROUND

Nematodes are the most abundant metazoans in marine sediments, many of which are bacterivores; however, how habitat bacteria affect physiological outcomes in marine nematodes remains largely unknown.  RESULTS: Here, we used a Litoditis marina inbred line to assess how native bacteria modulate host nematode physiology. We characterized seasonal dynamic bacterial compositions in L. marina habitats and examined the impacts of 448 habitat bacteria isolates on L. marina development, then focused on HQbiome with 73 native bacteria, of which we generated 72 whole genomes sequences. Unexpectedly, we found that the effects of marine native bacteria on the development of L. marina and its terrestrial relative Caenorhabditis elegans were significantly positively correlated. Next, we reconstructed bacterial metabolic networks and identified several bacterial metabolic pathways positively correlated with L. marina development (e.g., ubiquinol and heme b biosynthesis), while pyridoxal 5'-phosphate biosynthesis pathway was negatively associated. Through single metabolite supplementation, we verified CoQ10, heme b, acetyl-CoA, and acetaldehyde promoted L. marina development, while vitamin B6 attenuated growth. Notably, we found that only four development correlated metabolic pathways were shared between L. marina and C. elegans. Furthermore, we identified two bacterial metabolic pathways correlated with L. marina lifespan, while a distinct one in C. elegans. Strikingly, we found that glycerol supplementation significantly extended L. marina but not C. elegans longevity. Moreover, we comparatively demonstrated the distinct gut microbiota characteristics and their effects on L. marina and C. elegans physiology.

CONCLUSIONS

Given that both bacteria and marine nematodes are dominant taxa in sedimentary ecosystems, the resource presented here will provide novel insights to identify mechanisms underpinning how habitat bacteria affect nematode biology in a more natural context. Our integrative approach will provide a microbe-nematodes framework for microbiome mediated effects on host animal fitness.

Metabolic interactions shape emergent biofilm structures in a conceptual model of gut mucosal bacterial communities

The gut microbiome plays a major role in human health; however, little is known about the structural arrangement of microbes and factors governing their distribution. In this work, we present an in silico agent-based model (ABM) to conceptually simulate the dynamics of gut mucosal bacterial communities. We explored how various types of metabolic interactions, including competition, neutralism, commensalism, and mutualism, affect community structure, through nutrient consumption and metabolite exchange. Results showed that, across scenarios with different initial species abundances, cross-feeding promotes species coexistence. Morphologically, competition and neutralism resulted in segregation, while mutualism and commensalism fostered high intermixing. In addition, cooperative relations resulted in community properties with little sensitivity to the selective uptake of metabolites produced by the host. Moreover, metabolic interactions strongly influenced colonization success following the invasion of newcomer species. These results provide important insights into the utility of ABM in deciphering complex microbiome patterns.

Colon or semicolon: gut sampling microdevices for omics insights

Ingestible microdevices represent a breakthrough in non-invasive sampling of the human gastrointestinal (GI) tract. By capturing the native spatiotemporal microbiome and intricate biochemical gradients, these devices allow a non-invasive multi-omic access to the unperturbed host-microbiota crosstalk, immune/nutritional landscapes and gut-organ connections. We present the current progress of GI sampling microdevices towards personalized metabolism and fostering collaboration among clinicians, engineers, and data scientists.

New Drinking Water Genome Catalog Identifies a Globally Distributed Bacterial Genus Adapted to Disinfected Drinking Water Systems

Genome-resolved insights into the structure and function of the drinking water microbiome can advance the effective management of drinking water quality. To enable this, we constructed and curated thousands of metagenome-assembled and isolate genomes from drinking water distribution systems globally to develop a Drinking Water Genome Catalog (DWGC). The current DWGC disproportionately represents disinfected drinking water systems due to a paucity of metagenomes from nondisinfected systems. Using the DWGC, we identify core genera of the drinking water microbiome including a genus (UBA4765) within the order Rhizobiales that is frequently detected and highly abundant in disinfected drinking water systems. We demonstrate that this genus has been widely detected but incorrectly classified in previous amplicon sequencing-based investigations of the drinking water microbiome. Further, we show that a single genome variant (genomovar) within this genus is detected in 75% of drinking water systems included in this study. We propose a name for this uncultured bacterium as "Raskinella chloraquaticus" and describe the genus as "Raskinella" (endorsed by SeqCode). Metabolic annotation and modeling-based predictions indicate that this bacterium is capable of necrotrophic growth, is able to metabolize halogenated compounds, proliferates in a biofilm-based environment, and shows clear indications of disinfection-mediated selection.

DyMMM-LEAPS: An ML-based framework for modulating evenness and stability in synthetic microbial communities

There have been a growing number of computational strategies to aid in the design of synthetic microbial consortia. A framework to identify regions in parametric space to maximize two essential properties, evenness and stability, is critical. In this study, we introduce DyMMM-LEAPS (dynamic multispecies metabolic modeling-locating evenness and stability in large parametric space), an extension of the DyMMM framework. Our method explores the large parametric space of genetic circuits in synthetic microbial communities to identify regions of evenness and stability. Due to the high computational costs of exhaustive sampling, we utilize adaptive sampling and surrogate modeling to reduce the number of simulations required to map the vast space. Our framework predicts engineering targets and computes their operating ranges to maximize the probability of the engineered community to have high evenness and stability. We demonstrate our approach by simulating five cocultures and one three-strain culture with different social interactions (cooperation, competition, and predation) employing quorum-sensing-based genetic circuits. In addition to guiding circuit tuning, our pipeline gives an opportunity for a detailed analysis of pockets of evenness and stability for the circuit under investigation, which can further help dissect the relationship between the two properties. DyMMM-LEAPS is easily customizable and can be expanded to a larger community with more complex interactions.

The Future Exploring of Gut Microbiome-Immunity Interactions: From In Vivo/Vitro Models to In Silico Innovations

Investigating the complex interactions between microbiota and immunity is crucial for a fruitful understanding progress of human health and disease. This review assesses animal models, next-generation in vitro models, and in silico approaches that are used to decipher the microbiome-immunity axis, evaluating their strengths and limitations. While animal models provide a comprehensive biological context, they also raise ethical and practical concerns. Conversely, modern in vitro models reduce animal involvement but require specific costs and materials. When considering the environmental impact of these models, in silico approaches emerge as promising for resource reduction, but they require robust experimental validation and ongoing refinement. Their potential is significant, paving the way for a more sustainable and ethical future in microbiome-immunity research.

Measuring and understanding information storage and transfer in a simulated human gut microbiome

Considering biological systems as information processing entities and analyzing their organizational structure via information-theoretic measures has become an established approach in life sciences. We transfer this framework to a field of broad general interest, the human gut microbiome. We use BacArena, a software combining agent-based modelling and flux-balance analysis, to simulate a simplified human intestinal microbiome (SIHUMI). In a first step, we derive information theoretic measures from the simulated abundance data, and, in a second step, relate them to the metabolic processes underlying the abundance data. Our study provides further evidence on the role of active information storage as an indicator of unexpected structural change in the observed system. Besides, we show that information transfer reflects coherent behavior in the microbial community, both as a reaction to environmental changes and as a result of direct effective interaction. In this sense, purely abundance-based information theoretic measures can provide meaningful insight on metabolic interactions within bacterial communities. Furthermore, we shed light on the important however little noticed technical aspect of distinguishing immediate and delayed effects in the interpretation of local information theoretical measures.

Synthetic microbial communities: A gateway to understanding resistance, resilience, and functionality in spontaneously fermented food microbiomes

This review delves into the intricate traits of microbial communities encountered in spontaneously fermented foods (SFF), contributing to resistance, resilience, and functionality drivers. Traits of SFF microbiomes comprise of fluctuations in community composition, genetic stability, and condition-specific phenotypes. Synthetic microbial communities (SMCs) serve as a portal for mechanistic insights and strategic re-programming of microbial communities. Current literature underscores the pivotal role of microbiomes in SFF in shaping quality attributes and preserving the cultural heritage of their origin. In contrast to starter driven fermentations that tend to be more controlled but lacking the capacity to maintain or reproduce the complex flavors and intricacies found in SFF. SMCs, therefore, become indispensable tools, providing a nuanced understanding and control over fermented food microbiomes. They empower the prediction and engineering of microbial interactions and metabolic pathways with the aim of optimizing outcomes in food processing. Summarizing the current application of SMCs in fermented foods, there is still space for improvement. Challenges in achieving stability and reproducibility in SMCs are identified, stemming from non-standardized approaches. The future direction should involve embracing standardized protocols, advanced monitoring tools, and synthetic biology applications. A holistic, multi-disciplinary approach is paramount to unleashing the full potential of SMCs and fostering sustainable and innovative applications in fermented food systems.

Biobased short chain fatty acid production - Exploring microbial community dynamics and metabolic networks through kinetic and microbial modeling approaches

In recent years, there has been growing interest in harnessing anaerobic digestion technology for resource recovery from waste streams. This approach has evolved beyond its traditional role in energy generation to encompass the production of valuable carboxylic acids, especially volatile fatty acids (VFAs) like acetic acid, propionic acid, and butyric acid. VFAs hold great potential for various industries and biobased applications due to their versatile properties. Despite increasing global demand, over 90% of VFAs are currently produced synthetically from petrochemicals. Realizing the potential of large-scale biobased VFA production from waste streams offers significant eco-friendly opportunities but comes with several key challenges. These include low VFA production yields, unstable acid compositions, complex and expensive purification methods, and post-processing needs. Among these, production yield and acid composition stand out as the most critical obstacles impacting economic viability and competitiveness. This paper seeks to offer a comprehensive view of combining complementary modeling approaches, including kinetic and microbial modeling, to understand the workings of microbial communities and metabolic pathways in VFA production, enhance production efficiency, and regulate acid profiles through the integration of omics and bioreactor data.

Genome-resolved metagenomics: a game changer for microbiome medicine

Recent substantial evidence implicating commensal bacteria in human diseases has given rise to a new domain in biomedical research: microbiome medicine. This emerging field aims to understand and leverage the human microbiota and derivative molecules for disease prevention and treatment. Despite the complex and hierarchical organization of this ecosystem, most research over the years has relied on 16S amplicon sequencing, a legacy of bacterial phylogeny and taxonomy. Although advanced sequencing technologies have enabled cost-effective analysis of entire microbiota, translating the relatively short nucleotide information into the functional and taxonomic organization of the microbiome has posed challenges until recently. In the last decade, genome-resolved metagenomics, which aims to reconstruct microbial genomes directly from whole-metagenome sequencing data, has made significant strides and continues to unveil the mysteries of various human-associated microbial communities. There has been a rapid increase in the volume of whole metagenome sequencing data and in the compilation of novel metagenome-assembled genomes and protein sequences in public depositories. This review provides an overview of the capabilities and methods of genome-resolved metagenomics for studying the human microbiome, with a focus on investigating the prokaryotic microbiota of the human gut. Just as decoding the human genome and its variations marked the beginning of the genomic medicine era, unraveling the genomes of commensal microbes and their sequence variations is ushering us into the era of microbiome medicine. Genome-resolved metagenomics stands as a pivotal tool in this transition and can accelerate our journey toward achieving these scientific and medical milestones.

Microbiome modeling: a beginner's guide

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

Emerging methods for genome-scale metabolic modeling of microbial communities

Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels.

Surrogate modeling and control of medical digital twins

The vision of personalized medicine is to identify interventions that maintain or restore a person's health based on their individual biology. Medical digital twins, computational models that integrate a wide range of health-related data about a person and can be dynamically updated, are a key technology that can help guide medical decisions. Such medical digital twin models can be high-dimensional, multi-scale, and stochastic. To be practical for healthcare applications, they often need to be simplified into low-dimensional surrogate models that can be used for optimal design of interventions. This paper introduces surrogate modeling algorithms for the purpose of optimal control applications. As a use case, we focus on agent-based models (ABMs), a common model type in biomedicine for which there are no readily available optimal control algorithms. By deriving surrogate models that are based on systems of ordinary differential equations, we show how optimal control methods can be employed to compute effective interventions, which can then be lifted back to a given ABM. The relevance of the methods introduced here extends beyond medical digital twins to other complex dynamical systems.

Gut-associated functions are favored during microbiome assembly across a major part of C. elegans life

The microbiome expresses a variety of functions that influence host biology. The range of functions depends on the microbiome's composition, which can change during the host's lifetime due to neutral assembly processes, host-mediated selection, and environmental conditions. To date, the exact dynamics of microbiome assembly, the underlying determinants, and the effects on host-associated functions remain poorly understood. Here, we used the nematode Caenorhabditis elegans and a defined community of fully sequenced, naturally associated bacteria to study microbiome dynamics and functions across a major part of the worm's lifetime of hosts under controlled experimental conditions. Bacterial community composition initially shows strongly declining levels of stochasticity, which increases during later time points, suggesting selective effects in younger animals as opposed to more random processes in older animals. The adult microbiome is enriched in genera Ochrobactrum and Enterobacter compared to the direct substrate and a host-free control environment. Using pathway analysis, metabolic, and ecological modeling, we further find that the lifetime assembly dynamics increase competitive strategies and gut-associated functions in the host-associated microbiome, indicating that the colonizing bacteria benefit the worm. Overall, our study introduces a framework for studying microbiome assembly dynamics based on stochastic, ecological, and metabolic models, yielding new insights into the processes that determine host-associated microbiome composition and function.

IMPORTANCE

The microbiome plays a crucial role in host biology. Its functions depend on the microbiome composition that can change during a host's lifetime. To date, the dynamics of microbiome assembly and the resulting functions still need to be better understood. This study introduces a new approach to characterize the functional consequences of microbiome assembly by modeling both the relevance of stochastic processes and metabolic characteristics of microbial community changes. The approach was applied to experimental time-series data obtained for the microbiome of the nematode Caenorhabditis elegans across the major part of its lifetime. Stochastic processes played a minor role, whereas beneficial bacteria as well as gut-associated functions enriched in hosts. This indicates that the host might actively shape the composition of its microbiome. Overall, this study provides a framework for studying microbiome assembly dynamics and yields new insights into C. elegans microbiome functions.

BN-BacArena: Bayesian network extension of BacArena for the dynamic simulation of microbial communities

MOTIVATION

Simulating gut microbial dynamics is extremely challenging. Several computational tools, notably the widely used BacArena, enable modeling of dynamic changes in the microbial environment. These methods, however, do not comprehensively account for microbe-microbe stimulant or inhibitory effects or for nutrient-microbe inhibitory effects, typically observed in different compounds present in the daily diet.

RESULTS

Here, we present BN-BacArena, an extension of BacArena consisting on the incorporation within the native computational framework of a Bayesian network model that accounts for microbe-microbe and nutrient-microbe interactions. Using in vitro experiments, 16S rRNA gene sequencing data and nutritional composition of 55 foods, the output Bayesian network showed 23 significant nutrient-bacteria interactions, suggesting the importance of compounds such as polyols, ascorbic acid, polyphenols and other phytochemicals, and 40 bacteria-bacteria significant relationships. With test data, BN-BacArena demonstrates a statistically significant improvement over BacArena to predict the time-dependent relative abundance of bacterial species involved in the gut microbiota upon different nutritional interventions. As a result, BN-BacArena opens new avenues for the dynamic modeling and simulation of the human gut microbiota metabolism.

AVAILABILITY AND IMPLEMENTATION

MATLAB and R code are available in https://github.com/PlanesLab/BN-BacArena.

Rational construction of synthetic consortia: Key considerations and model-based methods for guiding the development of a novel biosynthesis platform

The rapid development of synthetic biology has significantly improved the capabilities of mono-culture systems in converting different substrates into various value-added bio-chemicals through metabolic engineering. However, overexpression of biosynthetic pathways in recombinant strains can impose a heavy metabolic burden on the host, resulting in imbalanced energy distribution and negatively affecting both cell growth and biosynthesis capacity. Synthetic consortia, consisting of two or more microbial species or strains with complementary functions, have emerged as a promising and efficient platform to alleviate the metabolic burden and increase product yield. However, research on synthetic consortia is still in its infancy, with numerous challenges regarding the design and construction of stable synthetic consortia. This review provides a comprehensive comparison of the advantages and disadvantages of mono-culture systems and synthetic consortia. Key considerations for engineering synthetic consortia based on recent advances are summarized, and simulation and computational tools for guiding the advancement of synthetic consortia are discussed. Moreover, further development of more efficient and cost-effective synthetic consortia with emerging technologies such as artificial intelligence and machine learning is highlighted.

Spatial transcriptome-guided multi-scale framework connects P. aeruginosa metabolic states to oxidative stress biofilm microenvironment

With the generation of spatially resolved transcriptomics of microbial biofilms, computational tools can be used to integrate this data to elucidate the multi-scale mechanisms controlling heterogeneous biofilm metabolism. This work presents a Multi-scale model of Metabolism In Cellular Systems (MiMICS) which is a computational framework that couples a genome-scale metabolic network reconstruction (GENRE) with Hybrid Automata Library (HAL), an existing agent-based model and reaction-diffusion model platform. A key feature of MiMICS is the ability to incorporate multiple -omics-guided metabolic models, which can represent unique metabolic states that yield different metabolic parameter values passed to the extracellular models. We used MiMICS to simulate Pseudomonas aeruginosa regulation of denitrification and oxidative stress metabolism in hypoxic and nitric oxide (NO) biofilm microenvironments. Integration of P. aeruginosa PA14 biofilm spatial transcriptomic data into a P. aeruginosa PA14 GENRE generated four PA14 metabolic model states that were input into MiMICS. Characteristic of aerobic, denitrification, and oxidative stress metabolism, the four metabolic model states predicted different oxygen, nitrate, and NO exchange fluxes that were passed as inputs to update the agent's local metabolite concentrations in the extracellular reaction-diffusion model. Individual bacterial agents chose a PA14 metabolic model state based on a combination of stochastic rules, and agents sensing local oxygen and NO. Transcriptome-guided MiMICS predictions suggested microscale denitrification and oxidative stress metabolic heterogeneity emerged due to local variability in the NO biofilm microenvironment. MiMICS accurately predicted the biofilm's spatial relationships between denitrification, oxidative stress, and central carbon metabolism. As simulated cells responded to extracellular NO, MiMICS revealed dynamics of cell populations heterogeneously upregulating reactions in the denitrification pathway, which may function to maintain NO levels within non-toxic ranges. We demonstrated that MiMICS is a valuable computational tool to incorporate multiple -omics-guided metabolic models to mechanistically map heterogeneous microbial metabolic states to the biofilm microenvironment.

Control of Medical Digital Twins with Artificial Neural Networks

The objective of personalized medicine is to tailor interventions to an individual patient's unique characteristics. A key technology for this purpose involves medical digital twins, computational models of human biology that can be personalized and dynamically updated to incorporate patient-specific data collected over time. Certain aspects of human biology, such as the immune system, are not easily captured with physics-based models, such as differential equations. Instead, they are often multi-scale, stochastic, and hybrid. This poses a challenge to existing model-based control and optimization approaches that cannot be readily applied to such models. Recent advances in automatic differentiation and neural-network control methods hold promise in addressing complex control problems. However, the application of these approaches to biomedical systems is still in its early stages. This work introduces dynamics-informed neural-network controllers as an alternative approach to control of medical digital twins. As a first use case for this method, the focus is on agent-based models, a versatile and increasingly common modeling platform in biomedicine. The effectiveness of the proposed neural-network control method is illustrated and benchmarked against other methods with two widely-used agent-based model types. The relevance of the method introduced here extends beyond medical digital twins to other complex dynamical systems.

Progress on network modeling and analysis of gut microecology: a review

The gut microecological network is a complex microbial community within the human body that plays a key role in linking dietary nutrition and host physiology. To understand the complex relationships among microbes and their functions within this community, network analysis has emerged as a powerful tool. By representing the interactions between microbes and their associated omics data as a network, we can gain a comprehensive understanding of the ecological mechanisms that drive the human gut microbiota. In addition, the network-based approach provides a more intuitive analysis of the gut microbiota, simplifying the study of its complex dynamics and interdependencies. This review provides a comprehensive overview of the methods used to construct and analyze networks in the context of gut microecological background. We discuss various types of network modeling approaches, including co-occurrence networks, causal networks, dynamic networks, and multi-omics networks, and describe the analytical techniques used to identify important network properties. We also highlight the challenges and limitations of network modeling in this area, such as data scarcity and heterogeneity, and provide future research directions to overcome these limitations. By exploring these network-based methods, researchers can gain valuable insights into the intricate relationships and functional roles of microbial communities within the gut, ultimately advancing our understanding of the gut microbiota's impact on human health.

Biophysical metabolic modeling of complex bacterial colony morphology

Microbial colony growth is shaped by the physics of biomass propagation and nutrient diffusion, and by the metabolic reactions that organisms activate as a function of the surrounding environment. While microbial colonies have been explored using minimal models of growth and motility, full integration of biomass propagation and metabolism is still lacking. Here, building upon our framework for Computation of Microbial Ecosystems in Time and Space (COMETS), we combine dynamic flux balance modeling of metabolism with collective biomass propagation and demographic fluctuations to provide nuanced simulations of E. coli colonies. Simulations produced realistic colony morphology, consistent with our experiments. They characterize the transition between smooth and furcated colonies and the decay of genetic diversity. Furthermore, we demonstrate that under certain conditions, biomass can accumulate along "metabolic rings" that are reminiscent of coffee-stain rings, but have a completely different origin. Our approach is a key step towards predictive microbial ecosystems modeling.

Life at the borderlands: microbiomes of interfaces critical to One Health

Microbiomes are foundational components of the environment that provide essential services relating to food security, carbon sequestration, human health, and the overall well-being of ecosystems. Microbiota exert their effects primarily through complex interactions at interfaces with their plant, animal, and human hosts, as well as within the soil environment. This review aims to explore the ecological, evolutionary, and molecular processes governing the establishment and function of microbiome-host relationships, specifically at interfaces critical to One Health-a transdisciplinary framework that recognizes that the health outcomes of people, animals, plants, and the environment are tightly interconnected. Within the context of One Health, the core principles underpinning microbiome assembly will be discussed in detail, including biofilm formation, microbial recruitment strategies, mechanisms of microbial attachment, community succession, and the effect these processes have on host function and health. Finally, this review will catalogue recent advances in microbiology and microbial ecology methods that can be used to profile microbial interfaces, with particular attention to multi-omic, advanced imaging, and modelling approaches. These technologies are essential for delineating the general and specific principles governing microbiome assembly and functions, mapping microbial interconnectivity across varying spatial and temporal scales, and for the establishment of predictive frameworks that will guide the development of targeted microbiome-interventions to deliver One Health outcomes.

Metabolic model predictions enable targeted microbiome manipulation through precision prebiotics

While numerous health-beneficial interactions between host and microbiota have been identified, there is still a lack of targeted approaches for modulating these interactions. Thus, we here identify precision prebiotics that specifically modulate the abundance of a microbiome member species of interest. In the first step, we show that defining precision prebiotics by compounds that are only taken up by the target species but no other species in a community is usually not possible due to overlapping metabolic niches. Subsequently, we use metabolic modeling to identify precision prebiotics for a two-member Caenorhabditis elegans microbiome community comprising the immune-protective target species Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71. We experimentally confirm four of the predicted precision prebiotics, L-serine, L-threonine, D-mannitol, and γ-aminobutyric acid, to specifically increase the abundance of MYb11. L-serine was further assessed in vivo, leading to an increase in MYb11 abundance also in the worm host. Overall, our findings demonstrate that metabolic modeling is an effective tool for the design of precision prebiotics as an important cornerstone for future microbiome-targeted therapies.IMPORTANCEWhile various mechanisms through which the microbiome influences disease processes in the host have been identified, there are still only few approaches that allow for targeted manipulation of microbiome composition as a first step toward microbiome-based therapies. Here, we propose the concept of precision prebiotics that allow to boost the abundance of already resident health-beneficial microbial species in a microbiome. We present a constraint-based modeling pipeline to predict precision prebiotics for a minimal microbial community in the worm Caenorhabditis elegans comprising the host-beneficial Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71 with the aim to boost the growth of MYb11. Experimentally testing four of the predicted precision prebiotics, we confirm that they are specifically able to increase the abundance of MYb11 in vitro and in vivo. These results demonstrate that constraint-based modeling could be an important tool for the development of targeted microbiome-based therapies against human diseases.

Changes in Bacterial Gut Composition in Parkinson's Disease and Their Metabolic Contribution to Disease Development: A Gut Community Reconstruction Approach

Parkinson's disease (PD) is a chronic and progressive neurodegenerative disease with the major symptoms comprising loss of movement coordination (motor dysfunction) and non-motor dysfunction, including gastrointestinal symptoms. Alterations in the gut microbiota composition have been reported in PD patients vs. controls. However, it is still unclear how these compositional changes contribute to disease etiology and progression. Furthermore, most of the available studies have focused on European, Asian, and North American cohorts, but the microbiomes of PD patients in Latin America have not been characterized. To address this problem, we obtained fecal samples from Colombian participants (n = 25 controls, n = 25 PD idiopathic cases) to characterize the taxonomical community changes during disease via 16S rRNA gene sequencing. An analysis of differential composition, diversity, and personalized computational modeling was carried out, given the fecal bacterial composition and diet of each participant. We found three metabolites that differed in dietary habits between PD patients and controls: carbohydrates, trans fatty acids, and potassium. We identified six genera that changed significantly in their relative abundance between PD patients and controls, belonging to the families Lachnospiraceae, Lactobacillaceae, Verrucomicrobioaceae, Peptostreptococcaceae, and Streptococcaceae. Furthermore, personalized metabolic modeling of the gut microbiome revealed changes in the predicted production of seven metabolites (Indole, tryptophan, fructose, phenylacetic acid, myristic acid, 3-Methyl-2-oxovaleric acid, and N-Acetylneuraminic acid). These metabolites are associated with the metabolism of aromatic amino acids and their consumption in the diet. Therefore, this research suggests that each individual's diet and intestinal composition could affect host metabolism. Furthermore, these findings open the door to the study of microbiome-host interactions and allow us to contribute to personalized medicine.

Is it selfish to be filamentous in biofilms? Individual-based modeling links microbial growth strategies with morphology using the new and modular iDynoMiCS 2.0

Microbial communities are found in all habitable environments and often occur in assemblages with self-organized spatial structures developing over time. This complexity can only be understood, predicted, and managed by combining experiments with mathematical modeling. Individual-based models are particularly suited if individual heterogeneity, local interactions, and adaptive behavior are of interest. Here we present the completely overhauled software platform, the individual-based Dynamics of Microbial Communities Simulator, iDynoMiCS 2.0, which enables researchers to specify a range of different models without having to program. Key new features and improvements are: (1) Substantially enhanced ease of use (graphical user interface, editor for model specification, unit conversions, data analysis and visualization and more). (2) Increased performance and scalability enabling simulations of up to 10 million agents in 3D biofilms. (3) Kinetics can be specified with any arithmetic function. (4) Agent properties can be assembled from orthogonal modules for pick and mix flexibility. (5) Force-based mechanical interaction framework enabling attractive forces and non-spherical agent morphologies as an alternative to the shoving algorithm. The new iDynoMiCS 2.0 has undergone intensive testing, from unit tests to a suite of increasingly complex numerical tests and the standard Benchmark 3 based on nitrifying biofilms. A second test case was based on the "biofilms promote altruism" study previously implemented in BacSim because competition outcomes are highly sensitive to the developing spatial structures due to positive feedback between cooperative individuals. We extended this case study by adding morphology to find that (i) filamentous bacteria outcompete spherical bacteria regardless of growth strategy and (ii) non-cooperating filaments outcompete cooperating filaments because filaments can escape the stronger competition between themselves. In conclusion, the new substantially improved iDynoMiCS 2.0 joins a growing number of platforms for individual-based modeling of microbial communities with specific advantages and disadvantages that we discuss, giving users a wider choice.

A genome-wide association study reveals the relationship between human genetic variation and the nasal microbiome

The nasal cavity harbors diverse microbiota that contributes to human health and respiratory diseases. However, whether and to what extent the host genome shapes the nasal microbiome remains largely unknown. Here, by dissecting the human genome and nasal metagenome data from 1401 healthy individuals, we demonstrated that the top three host genetic principal components strongly correlated with the nasal microbiota diversity and composition. The genetic association analyses identified 63 genome-wide significant loci affecting the nasal microbial taxa and functions, of which 2 loci reached study-wide significance (p < 1.7 × 10-10): rs73268759 within CAMK2A associated with genus Actinomyces and family Actinomycetaceae; and rs35211877 near POM121L12 with Gemella asaccharolytica. In addition to respiratory-related diseases, the associated loci are mainly implicated in cardiometabolic or neuropsychiatric diseases. Functional analysis showed the associated genes were most significantly expressed in the nasal airway epithelium tissue and enriched in the calcium signaling and hippo signaling pathway. Further observational correlation and Mendelian randomization analyses consistently suggested the causal effects of Serratia grimesii and Yokenella regensburgei on cardiometabolic biomarkers (cystine, glutamic acid, and creatine). This study suggested that the host genome plays an important role in shaping the nasal microbiome.

Predicting microbial interactions with approaches based on flux balance analysis: an evaluation

BACKGROUND

Given a genome-scale metabolic model (GEM) of a microorganism and criteria for optimization, flux balance analysis (FBA) predicts the optimal growth rate and its corresponding flux distribution for a specific medium. FBA has been extended to microbial consortia and thus can be used to predict interactions by comparing in-silico growth rates for co- and monocultures. Although FBA-based methods for microbial interaction prediction are becoming popular, a systematic evaluation of their accuracy has not yet been performed.

RESULTS

Here, we evaluate the accuracy of FBA-based predictions of human and mouse gut bacterial interactions using growth data from the literature. For this, we collected 26 GEMs from the semi-curated AGORA database as well as four previously published curated GEMs. We tested the accuracy of three tools (COMETS, Microbiome Modeling Toolbox and MICOM) by comparing growth rates predicted in mono- and co-culture to growth rates extracted from the literature and also investigated the impact of different tool settings and media. We found that except for curated GEMs, predicted growth rates and their ratios (i.e. interaction strengths) do not correlate with growth rates and interaction strengths obtained from in vitro data.

CONCLUSIONS

Prediction of growth rates with FBA using semi-curated GEMs is currently not sufficiently accurate to predict interaction strengths reliably.

Simple Porifera holobiont reveals complex interactions between the host, an archaeon, a bacterium, and a phage

The basal metazoan phylum Porifera (sponges) is increasingly used as a model to investigate ecological and evolutionary features of microbe-animal symbioses. However, sponges often host complex microbiomes, which has hampered our understanding of their interactions with their microbial symbionts. Here, we describe the discovery and characterization of the simplest sponge holobiont reported to date, consisting of the deep-sea glass sponge Aphrocallistes beatrix and two newly-described microbial symbionts: an autotrophic ammonia-oxidizing archaeon and a bacterial heterotroph. Omics analyses and metabolic modeling revealed the dependency of the ammonia-oxidizing archaea on sponge-derived ammonia to drive primary production, which in turn supports the bacterium's growth by providing the dicarboxylate fumarate. Furthermore, virus-mediated archaeal lysis appears crucial to overcome the bacterium's vitamin B12 auxotrophy. These findings reveal that the exchanges of vitamin B12 and dicarboxylate may be evolutionarily conserved features of symbiosis as they can also be found in interactions between free-living marine bacteria, and between microbes and plants or diatoms.

Microbial interactions from a new perspective: reinforcement learning reveals new insights into microbiome evolution

MOTIVATION

Microbes are essential part of all ecosystems, influencing material flow and shaping their surroundings. Metabolic modeling has been a useful tool and provided tremendous insights into microbial community metabolism. However, current methods based on flux balance analysis (FBA) usually fail to predict metabolic and regulatory strategies that lead to long-term survival and stability especially in heterogenous communities.

RESULTS

Here, we introduce a novel reinforcement learning algorithm, Self-Playing Microbes in Dynamic FBA, which treats microbial metabolism as a decision-making process, allowing individual microbial agents to evolve by learning and adapting metabolic strategies for enhanced long-term fitness. This algorithm predicts what microbial flux regulation policies will stabilize in the dynamic ecosystem of interest in the presence of other microbes with minimal reliance on predefined strategies. Throughout this article, we present several scenarios wherein our algorithm outperforms existing methods in reproducing outcomes, and we explore the biological significance of these predictions.

AVAILABILITY AND IMPLEMENTATION

The source code for this article is available at: https://github.com/chan-csu/SPAM-DFBA.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Acetate and succinate benefit host muscle energetics as exercise-associated post-biotics

Recently, the gut microbiome has emerged as a potent modulator of exercise-induced systemic adaptation and appears to be crucial for mediating some of the benefits of exercise. This study builds upon previous evidence establishing a gut microbiome-skeletal muscle axis, identifying exercise-induced changes in microbiome composition. Metagenomics sequencing of fecal samples from non-exercise-trained controls or exercise-trained mice was conducted. Biodiversity indices indicated exercise training did not change alpha diversity. However, there were notable differences in beta-diversity between trained and untrained microbiomes. Exercise significantly increased the level of the bacterial species Muribaculaceae bacterium DSM 103720. Computation simulation of bacterial growth was used to predict metabolites that accumulate under in silico culture of exercise-responsive bacteria. We identified acetate and succinate as potential gut microbial metabolites that are produced by Muribaculaceae bacterium, which were then administered to mice during a period of mechanical overload-induced muscle hypertrophy. Although no differences were observed for the overall muscle growth response to succinate or acetate administration during the first 5 days of mechanical overload-induced hypertrophy, acetate and succinate increased skeletal muscle mitochondrial respiration. When given as post-biotics, succinate or acetate treatment may improve oxidative metabolism during muscle hypertrophy.

PhysiBoSS 2.0: a sustainable integration of stochastic Boolean and agent-based modelling frameworks

In systems biology, mathematical models and simulations play a crucial role in understanding complex biological systems. Different modelling frameworks are employed depending on the nature and scales of the system under study. For instance, signalling and regulatory networks can be simulated using Boolean modelling, whereas multicellular systems can be studied using agent-based modelling. Herein, we present PhysiBoSS 2.0, a hybrid agent-based modelling framework that allows simulating signalling and regulatory networks within individual cell agents. PhysiBoSS 2.0 is a redesign and reimplementation of PhysiBoSS 1.0 and was conceived as an add-on that expands the PhysiCell functionalities by enabling the simulation of intracellular cell signalling using MaBoSS while keeping a decoupled, maintainable and model-agnostic design. PhysiBoSS 2.0 also expands the set of functionalities offered to the users, including custom models and cell specifications, mechanistic submodels of substrate internalisation and detailed control over simulation parameters. Together with PhysiBoSS 2.0, we introduce PCTK, a Python package developed for handling and processing simulation outputs, and generating summary plots and 3D renders. PhysiBoSS 2.0 allows studying the interplay between the microenvironment, the signalling pathways that control cellular processes and population dynamics, suitable for modelling cancer. We show different approaches for integrating Boolean networks into multi-scale simulations using strategies to study the drug effects and synergies in models of cancer cell lines and validate them using experimental data. PhysiBoSS 2.0 is open-source and publicly available on GitHub with several repositories of accompanying interoperable tools.

Degenerative Cervical Myelopathy induces sex-specific dysbiosis in mice

Degenerative Cervical Myelopathy (DCM) is the most common cause of spinal cord impairment in elderly populations. It describes a spectrum of disorders that cause progressive spinal cord compression, neurological impairment, loss of bladder and bowel functions, and gastrointestinal dysfunction. The gut microbiota has been recognized as an environmental factor that can modulate both the function of the central nervous system and the immune response through the microbiota-gut-brain axis. Changes in gut microbiota composition or microbiota-producing factors have been linked to the progression and development of several pathologies. However, little is known about the potential role of the gut microbiota in the pathobiology of DCM. Here, DCM was induced in C57BL/6 mice by implanting an aromatic polyether material underneath the C5-6 laminae. The extent of DCM-induced changes in microbiota composition was assessed by 16S rRNA sequencing of the fecal samples. The immune cell composition was assessed using flow cytometry. To date, several bacterial members have been identified using BLAST against the largest collection of metagenome-derived genomes from the mouse gut. In both, female and males DCM caused gut dysbiosis compared to the sham group. However, dysbiosis was more pronounced in males than in females, and several bacterial members of the families Lachnospiraceae and Muribaculaceae were significantly altered in the DCM group. These changes were also associated with altered microbe-derived metabolic changes in propionate-, butyrate-, and lactate-producing bacterial members. Our results demonstrate that DCM causes dynamic changes over time in the gut microbiota, reducing the abundance of butyrate-producing bacteria, and lactate-producing bacteria to a lesser extent. Genome-scale metabolic modeling using gapseq successfully identified pyruvate-to-butanoate and pyruvate-to-propionate reactions involving genes such as Buk and ACH1, respectively. These results provide a better understanding of the sex-specific molecular effects of changes in the gut microbiota on DCM pathobiology.

A structured evaluation of genome-scale constraint-based modeling tools for microbial consortia

Harnessing the power of microbial consortia is integral to a diverse range of sectors, from healthcare to biotechnology to environmental remediation. To fully realize this potential, it is critical to understand the mechanisms behind the interactions that structure microbial consortia and determine their functions. Constraint-based reconstruction and analysis (COBRA) approaches, employing genome-scale metabolic models (GEMs), have emerged as the state-of-the-art tool to simulate the behavior of microbial communities from their constituent genomes. In the last decade, many tools have been developed that use COBRA approaches to simulate multi-species consortia, under either steady-state, dynamic, or spatiotemporally varying scenarios. Yet, these tools have not been systematically evaluated regarding their software quality, most suitable application, and predictive power. Hence, it is uncertain which tools users should apply to their system and what are the most urgent directions that developers should take in the future to improve existing capacities. This study conducted a systematic evaluation of COBRA-based tools for microbial communities using datasets from two-member communities as test cases. First, we performed a qualitative assessment in which we evaluated 24 published tools based on a list of FAIR (Findability, Accessibility, Interoperability, and Reusability) features essential for software quality. Next, we quantitatively tested the predictions in a subset of 14 of these tools against experimental data from three different case studies: a) syngas fermentation by C. autoethanogenum and C. kluyveri for the static tools, b) glucose/xylose fermentation with engineered E. coli and S. cerevisiae for the dynamic tools, and c) a Petri dish of E. coli and S. enterica for tools incorporating spatiotemporal variation. Our results show varying performance levels of the best qualitatively assessed tools when examining the different categories of tools. The differences in the mathematical formulation of the approaches and their relation to the results were also discussed. Ultimately, we provide recommendations for refining future GEM microbial modeling tools.

Regulated bacterial interaction networks: A mathematical framework to describe competitive growth under inclusion of metabolite cross-feeding

When bacterial species with the same resource preferences share the same growth environment, it is commonly believed that direct competition will arise. A large variety of competition and more general 'interaction' models have been formulated, but what is currently lacking are models that link monoculture growth kinetics and community growth under inclusion of emerging biological interactions, such as metabolite cross-feeding. In order to understand and mathematically describe the nature of potential cross-feeding interactions, we design experiments where two bacterial species Pseudomonas putida and Pseudomonas veronii grow in liquid medium either in mono- or as co-culture in a resource-limited environment. We measure population growth under single substrate competition or with double species-specific substrates (substrate 'indifference'), and starting from varying cell ratios of either species. Using experimental data as input, we first consider a mean-field model of resource-based competition, which captures well the empirically observed growth rates for monocultures, but fails to correctly predict growth rates in co-culture mixtures, in particular for skewed starting species ratios. Based on this, we extend the model by cross-feeding interactions where the consumption of substrate by one consumer produces metabolites that in turn are resources for the other consumer, thus leading to positive feedback in the species system. Two different cross-feeding options were considered, which either lead to constant metabolite cross-feeding, or to a regulated form, where metabolite utilization is activated with rates according to either a threshold or a Hill function, dependent on metabolite concentration. Both mathematical proof and experimental data indicate regulated cross-feeding to be the preferred model to constant metabolite utilization, with best co-culture growth predictions in case of high Hill coefficients, close to binary (on/off) activation states. This suggests that species use the appearing metabolite concentrations only when they are becoming high enough; possibly as a consequence of their lower energetic content than the primary substrate. Metabolite sharing was particularly relevant at unbalanced starting cell ratios, causing the minority partner to proliferate more than expected from the competitive substrate because of metabolite release from the majority partner. This effect thus likely quells immediate substrate competition and may be important in natural communities with typical very skewed relative taxa abundances and slower-growing taxa. In conclusion, the regulated bacterial interaction network correctly describes species substrate growth reactions in mixtures with few kinetic parameters that can be obtained from monoculture growth experiments.

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.

Metabolic interaction models recapitulate leaf microbiota ecology

Resource allocation affects the structure of microbiomes, including those associated with living hosts. Understanding the degree to which this dependency determines interspecies interactions may advance efforts to control host-microbiome relationships. We combined synthetic community experiments with computational models to predict interaction outcomes between plant-associated bacteria. We mapped the metabolic capabilities of 224 leaf isolates from Arabidopsis thaliana by assessing the growth of each strain on 45 environmentally relevant carbon sources in vitro. We used these data to build curated genome-scale metabolic models for all strains, which we combined to simulate >17,500 interactions. The models recapitulated outcomes observed in planta with >89% accuracy, highlighting the role of carbon utilization and the contributions of niche partitioning and cross-feeding in the assembly of leaf microbiomes.