Review and perspective on mathematical modeling of microbial ecosystems

Deciphering and designing microbial communities by genome-scale metabolic modelling

Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.

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.

Advances in modeling techniques for the production and purification of biomolecules: A comprehensive review

In response to the growing demand for therapeutic biomolecules, there is a need for continuous and cost-effective bio-separation techniques to enhance extraction yield and efficiency. Aqueous biphasic extractive fermentation has emerged as an integrated downstream processing technique, offering selective partitioning, high productivity, and preservation of biomolecule integrity. However, the dynamic nature of this technique requires a comprehensive understanding of the underlying separation mechanisms. Unfortunately, the analysis of parameters influencing this dynamic behavior can be challenging due to limited resources and time. To address this, mathematical modeling approaches can be employed to minimize the tedious trial-and-error experimentation process. This review article presents mathematical modeling approaches for both upstream and downstream processing techniques, focusing on the production of biomolecules which can be used in pharmaceutical industries in a cost-effective manner. By leveraging mathematical models, researchers can optimize the production and purification processes, leading to improved efficiency and processing cost reduction in biomolecule production.

A Review of Basic Bioinformatic Techniques for Microbial Community Analysis in an Anaerobic Digester

Biogas production involves various types of intricate microbial populations in an anaerobic digester (AD). To understand the anaerobic digestion system better, a broad-based study must be conducted on the microbial population. Deep understanding of the complete metagenomics including microbial structure, functional gene form, similarity/differences, and relationships between metabolic pathways and product formation, could aid in optimization and enhancement of AD processes. With advancements in technologies for metagenomic sequencing, for example, next generation sequencing and high-throughput sequencing, have revolutionized the study of microbial dynamics in anaerobic digestion. This review includes a brief introduction to the basic process of metagenomics research and includes a detailed summary of the various bioinformatics approaches, viz., total investigation of data obtained from microbial communities using bioinformatics methods to expose metagenomics characterization. This includes (1) methods of DNA isolation and sequencing, (2) investigation of anaerobic microbial communities using bioinformatics techniques, (3) application of the analysis of anaerobic microbial community and biogas production, and (4) restriction and prediction of bioinformatics analysis on microbial metagenomics. The review has been concluded, giving a summarized insight into bioinformatic tools and also promoting the future prospects of integrating humungous data with artificial intelligence and neural network software.

Population Response of Rhizosphere Microbiota of Garden Pea Genotypes to Inoculation with Arbuscular Mycorrhizal Fungi

This study of a legume's rhizosphere in tripartite symbiosis focused on the relationships between the symbionts and less on the overall rhizosphere microbiome. We used an experimental model with different garden pea genotypes inoculated with AM fungi (Rhizophagus irregularis and with a mix of AM species) to study their influence on the population levels of main trophic groups of soil microorganisms as well as their structure and functional relationships in the rhizosphere microbial community. The experiments were carried out at two phenological cycles of the plants. Analyzes were performed according to classical methods: microbial population density defined as CUF/g a.d.s. and root colonization rate with AMF (%). We found a proven dominant effect of AMF on the densities of micromycetes and actinomycetes in the direction of reduction, suggesting antagonism, and on ammonifying, phosphate-solubilizing and free-living diazotrophic Azotobacter bacteria in the direction of stimulation, an indicator of mutualistic relationships. We determined that the genotype was decisive for the formation of populations of bacteria immobilizing mineral NH4+-N and bacteria Rhizobium. We reported significant two-way relationships between trophic groups related associated with soil nitrogen and phosphorus ions availability. The preserved proportions between trophic groups in the microbial communities were indicative of structural and functional stability.

Calibrating spatiotemporal models of microbial communities to microscopy data: A review

Spatiotemporal models that account for heterogeneity within microbial communities rely on single-cell data for calibration and validation. Such data, commonly collected via microscopy and flow cytometry, have been made more accessible by recent advances in microfluidics platforms and data processing pipelines. However, validating models against such data poses significant challenges. Validation practices vary widely between modelling studies; systematic and rigorous methods have not been widely adopted. Similar challenges are faced by the (macrobial) ecology community, in which systematic calibration approaches are often employed to improve quantitative predictions from computational models. Here, we review single-cell observation techniques that are being applied to study microbial communities and the calibration strategies that are being employed for accompanying spatiotemporal models. To facilitate future calibration efforts, we have compiled a list of summary statistics relevant for quantifying spatiotemporal patterns in microbial communities. Finally, we highlight some recently developed techniques that hold promise for improved model calibration, including algorithmic guidance of summary statistic selection and machine learning approaches for efficient model simulation.

Species abundance correlations carry limited information about microbial network interactions

Unraveling the network of interactions in ecological communities is a daunting task. Common methods to infer interspecific interactions from cross-sectional data are based on co-occurrence measures. For instance, interactions in the human microbiome are often inferred from correlations between the abundances of bacterial phylogenetic groups across subjects. We tested whether such correlation-based methods are indeed reliable for inferring interaction networks. For this purpose, we simulated bacterial communities by means of the generalized Lotka-Volterra model, with variation in model parameters representing variability among hosts. Our results show that correlations can be indicative for presence of bacterial interactions, but only when measurement noise is low relative to the variation in interaction strengths between hosts. Indication of interaction was affected by type of interaction network, process noise and sampling under non-equilibrium conditions. The sign of a correlation mostly coincided with the nature of the strongest pairwise interaction, but this is not necessarily the case. For instance, under rare conditions of identical interaction strength, we found that competitive and exploitative interactions can result in positive as well as negative correlations. Thus, cross-sectional abundance data carry limited information on specific interaction types. Correlations in abundance may hint at interactions but require independent validation.

The bacteria in and on our body (the human microbiome) largely determine how our body functions, and whether we stay healthy or get sick. These bacteria do not live on their own, but interact among each other and with their human host. Finding out which bacteria interact with each other is cumbersome, but patterns of joint occurrence between species might provide a clue to their ecological dependencies. We investigated whether correlations in species abundance can be used for the purpose of ecological network reconstruction. We simulated different bacterial communities with known interactions according to a theoretical population model. After having collected virtual samples from our simulated data, we performed a correlation analysis and then compared the correlation network with our known interaction network. We found that correlations can be informative for underlying interactions, but ecological conclusions should be drawn carefully. An obvious limitation of correlation analysis is that direction of interaction cannot be recovered from co-occurrence data, making correlations insensitive for detection of asymmetric interactions. In addition, we found that competitive and exploitative interactions can induce positive as well as negative correlations. We recommend careful interpretation and validation when inferring networks from cross-sectional abundance data.

Complex and unexpected outcomes of antibiotic therapy against a polymicrobial infection

Antibiotics are our primary approach to treating complex infections, yet we have a poor understanding of how these drugs affect microbial communities. To better understand antimicrobial effects on host-associated microbial communities we treated cultured sputum microbiomes from people with cystic fibrosis (pwCF, n = 24) with 11 different antibiotics, supported by theoretical and mathematical modeling-based predictions in a mucus-plugged bronchiole microcosm. Treatment outcomes we identified in vitro that were predicted in silico were: 1) community death, 2) community resistance, 3) pathogen killing, and 4) fermenter killing. However, two outcomes that were not predicted when antibiotics were applied were 5) community profile shifts with little change in total bacterial load (TBL), and 6) increases in TBL. The latter outcome was observed in 17.8% of samples with a TBL increase of greater than 20% and 6.8% of samples with an increase greater than 40%, demonstrating significant increases in community carrying capacity in the presence of an antibiotic. An iteration of the mathematical model showed that TBL increase was due to antibiotic-mediated release of pH-dependent inhibition of pathogens by anaerobe fermentation. These dynamics were verified in vitro when killing of fermenters resulted in a higher community carrying capacity compared to a no antibiotic control. Metagenomic sequencing of sputum samples during antibiotic therapy revealed similar dynamics in clinical samples. This study shows that the complex microbial ecology dictates the outcomes of antibiotic therapy against a polymicrobial infection.

MicNet toolbox: Visualizing and unraveling a microbial network

Applications of network theory to microbial ecology are an emerging and promising approach to understanding both global and local patterns in the structure and interplay of these microbial communities. In this paper, we present an open-source python toolbox which consists of two modules: on one hand, we introduce a visualization module that incorporates the use of UMAP, a dimensionality reduction technique that focuses on local patterns, and HDBSCAN, a clustering technique based on density; on the other hand, we have included a module that runs an enhanced version of the SparCC code, sustaining larger datasets than before, and we couple the resulting networks with network theory analyses to describe the resulting co-occurrence networks, including several novel analyses, such as structural balance metrics and a proposal to discover the underlying topology of a co-occurrence network. We validated the proposed toolbox on 1) a simple and well described biological network of kombucha, consisting of 48 ASVs, and 2) we validate the improvements of our new version of SparCC. Finally, we showcase the use of the MicNet toolbox on a large dataset from Archean Domes, consisting of more than 2,000 ASVs. Our toolbox is freely available as a github repository (https://github.com/Labevo/MicNetToolbox), and it is accompanied by a web dashboard (http://micnetapplb-1212130533.us-east-1.elb.amazonaws.com) that can be used in a simple and straightforward manner with relative abundance data. This easy-to-use implementation is aimed to microbial ecologists with little to no experience in programming, while the most experienced bioinformatics will also be able to manipulate the source code's functions with ease.

Path-integral solution of MacArthur's resource-competition model for large ecosystems with random species-resources couplings

We solve MacArthur's resource-competition model with random species-resource couplings in the "thermodynamic" limit of infinitely many species and resources using dynamical path integrals à la De Domincis. We analyze how the steady state picture changes upon modifying several parameters, including the degree of heterogeneity of metabolic strategies (encoding the preferences of species) and of maximal resource levels (carrying capacities), and discuss its stability. Ultimately, the scenario obtained by other approaches is recovered by analyzing an effective one-species-one-resource ecosystem that is fully equivalent to the original multi-species one. The technique used here can be applied for the analysis of other model ecosystems related to the version of MacArthur's model considered here.

Modelling microbial communities: Harnessing consortia for biotechnological applications

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

Designing microbial communities to maximize the thermodynamic driving force for the production of chemicals

Microbial communities have become a major research focus due to their importance for biogeochemical cycles, biomedicine and biotechnological applications. While some biotechnological applications, such as anaerobic digestion, make use of naturally arising microbial communities, the rational design of microbial consortia for bio-based production processes has recently gained much interest. One class of synthetic microbial consortia is based on specifically designed strains of one species. A common design principle for these consortia is based on division of labor, where the entire production pathway is divided between the different strains to reduce the metabolic burden caused by product synthesis. We first show that classical division of labor does not automatically reduce the metabolic burden when metabolic flux per biomass is analyzed. We then present ASTHERISC (Algorithmic Search of THERmodynamic advantages in Single-species Communities), a new computational approach for designing multi-strain communities of a single-species with the aim to divide a production pathway between different strains such that the thermodynamic driving force for product synthesis is maximized. ASTHERISC exploits the fact that compartmentalization of segments of a product pathway in different strains can circumvent thermodynamic bottlenecks arising when operation of one reaction requires a metabolite with high and operation of another reaction the same metabolite with low concentration. We implemented the ASTHERISC algorithm in a dedicated program package and applied it on E. coli core and genome-scale models with different settings, for example, regarding number of strains or demanded product yield. These calculations showed that, for each scenario, many target metabolites (products) exist where a multi-strain community can provide a thermodynamic advantage compared to a single strain solution. In some cases, a production with sufficiently high yield is thermodynamically only feasible with a community. In summary, the developed ASTHERISC approach provides a promising new principle for designing microbial communities for the bio-based production of chemicals.

Communities of microbes are ubiquitous in nature and also of high relevance for industrial applications, e.g. for the production of biogas. The development and use of non-natural communities for biotechnological applications has become an important subject of research. In this work, we present a new computational method to design synthetic communities with improved capabilities for the synthesis of desired target metabolites. Our method takes a constraint-based metabolic model of an organism as input and searches for a suitable partitioning of the product pathway via different strains of the organism such that the thermodynamic driving force for product synthesis is maximized. Essentially, this approach exploits the fact that having multiple strains allows adjustment of different metabolite concentrations in the different strains by which the thermodynamic driving force for product synthesis can often be increased. We tested this approach with a core and with a genome-scale metabolic network model of Escherichia coli. We found that, for dozens of metabolites, there exist communities with specifically designed strains of E. coli where the maximal thermodynamic driving force can be increased compared to a single E. coli strain. In summary, our presented method provides a new approach, together with a new design principle, for the computational design of microbial communities.

Co-culturing microbial consortia: approaches for applications in biomanufacturing and bioprocessing

The application of microbial co-cultures is now recognized in the fields of biotechnology, ecology, and medicine. Understanding the biological interactions that govern the association of microorganisms would shape the way in which artificial/synthetic co-cultures or consortia are developed. The ability to accurately predict and control cell-to-cell interactions fully would be a significant enabler in synthetic biology. Co-culturing method development holds the key to strategically engineer environments in which the co-cultured microorganism can be monitored. Various approaches have been employed which aim to emulate the natural environment and gain access to the untapped natural resources emerging from cross-talk between partners. Amongst these methods are the use of a communal liquid medium for growth, use of a solid-liquid interface, membrane separation, spatial separation, and use of microfluidics systems. Maximizing the information content of interactions monitored is one of the major challenges that needs to be addressed by these designs. This review critically evaluates the significance and drawbacks of the co-culturing approaches used to this day in biotechnological applications, relevant to biomanufacturing. It is recommended that experimental results for a co-cultured species should be validated with different co-culture approaches due to variations in interactions that could exist as a result of the culturing method selected.

Mechanism Across Scales: A Holistic Modeling Framework Integrating Laboratory and Field Studies for Microbial Ecology

Over the last century, leaps in technology for imaging, sampling, detection, high-throughput sequencing, and -omics analyses have revolutionized microbial ecology to enable rapid acquisition of extensive datasets for microbial communities across the ever-increasing temporal and spatial scales. The present challenge is capitalizing on our enhanced abilities of observation and integrating diverse data types from different scales, resolutions, and disciplines to reach a causal and mechanistic understanding of how microbial communities transform and respond to perturbations in the environment. This type of causal and mechanistic understanding will make predictions of microbial community behavior more robust and actionable in addressing microbially mediated global problems. To discern drivers of microbial community assembly and function, we recognize the need for a conceptual, quantitative framework that connects measurements of genomic potential, the environment, and ecological and physical forces to rates of microbial growth at specific locations. We describe the Framework for Integrated, Conceptual, and Systematic Microbial Ecology (FICSME), an experimental design framework for conducting process-focused microbial ecology studies that incorporates biological, chemical, and physical drivers of a microbial system into a conceptual model. Through iterative cycles that advance our understanding of the coupling across scales and processes, we can reliably predict how perturbations to microbial systems impact ecosystem-scale processes or vice versa. We describe an approach and potential applications for using the FICSME to elucidate the mechanisms of globally important ecological and physical processes, toward attaining the goal of predicting the structure and function of microbial communities in chemically complex natural environments.

Future directions for studying resilience of the oral ecosystem

The oral ecosystem is shaped by complex interactions between systemic health disease and the resident oral microbiota. Research in the last two decades has produced datasets describing the genetics and physiology of the host and the oral microbiome in health and disease. There are inter-individual differences in the ability to tolerate oral disease-promoting challenges. Identification of the key factors that drive a healthy and resilient oral ecosystem is urgently needed. So far, progress is being made towards replicating the host-microbiota interplay in vitro. Clinical studies may shed light on the mechanisms of oral health resilience. However, most clinical studies are cross-sectional and are insufficient for understanding resilience or for identifying biomarkers that correlate with the point of transition from oral health to dysbiosis. Mathematical and computational models, including artificial intelligence approaches, offer an opportunity to inform the design of clinical studies by identifying key biomarkers and interaction networks in complex datasets and predicting important parameters. This paper discusses some of the challenges and opportunities for understanding the biological basis of resilience of the oral ecosystem. It discusses the current status and challenges, and proposes a way forward to better understand resilience towards oral diseases.

Metabolic modelling approaches for describing and engineering microbial communities

Microbes do not live in isolation but in microbial communities. The relevance of microbial communities is increasing due to growing awareness of their influence on a huge number of environmental, health and industrial processes. Hence, being able to control and engineer the output of both natural and synthetic communities would be of great interest. However, most of the available methods and biotechnological applications involving microorganisms, both in vivo and in silico, have been developed in the context of isolated microbes. In vivo microbial consortia development is extremely difficult and costly because it implies replicating suitable environments in the wet-lab. Computational approaches are thus a good, cost-effective alternative to study microbial communities, mainly via descriptive modelling, but also via engineering modelling. In this review we provide a detailed compilation of examples of engineered microbial communities and a comprehensive, historical revision of available computational metabolic modelling methods to better understand, and rationally engineer wild and synthetic microbial communities.

Experimental and computational approaches to unravel microbial community assembly

Microbial communities have a preponderant role in the life support processes of our common home planet Earth. These extremely diverse communities drive global biogeochemical cycles, and develop intimate relationships with most multicellular organisms, with a significant impact on their fitness. Our understanding of their composition and function has enjoyed a significant thrust during the last decade thanks to the rise of high-throughput sequencing technologies. Intriguingly, the diversity patterns observed in nature point to the possible existence of fundamental community assembly rules. Unfortunately, these rules are still poorly understood, despite the fact that their knowledge could spur a scientific, technological, and economic revolution, impacting, for instance, agricultural, environmental, and health-related practices. In this minireview, I recapitulate the most important wet lab techniques and computational approaches currently employed in the study of microbial community assembly, and briefly discuss various experimental designs. Most of these approaches and considerations are also relevant to the study of microbial microevolution, as it has been shown that it can occur in ecological relevant timescales. Moreover, I provide a succinct review of various recent studies, chosen based on the diversity of ecological concepts addressed, experimental designs, and choice of wet lab and computational techniques. This piece aims to serve as a primer to those new to the field, as well as a source of new ideas to the more experienced researchers.

From sequence to information

Today massive amounts of sequenced metagenomic and metatranscriptomic data from different ecological niches and environmental locations are available. Scientific progress depends critically on methods that allow extracting useful information from the various types of sequence data. Here, we will first discuss types of information contained in the various flavours of biological sequence data, and how this information can be interpreted to increase our scientific knowledge and understanding. We argue that a mechanistic understanding of biological systems analysed from different perspectives is required to consistently interpret experimental observations, and that this understanding is greatly facilitated by the generation and analysis of dynamic mathematical models. We conclude that, in order to construct mathematical models and to test mechanistic hypotheses, time-series data are of critical importance. We review diverse techniques to analyse time-series data and discuss various approaches by which time-series of biological sequence data have been successfully used to derive and test mechanistic hypotheses. Analysing the bottlenecks of current strategies in the extraction of knowledge and understanding from data, we conclude that combined experimental and theoretical efforts should be implemented as early as possible during the planning phase of individual experiments and scientific research projects.

This article is part of the theme issue ‘Integrative research perspectives on marine conservation’.

Dynamics of microbial competition, commensalism, and cooperation and its implications for coculture and microbiome engineering

Microbial consortium is a complex adaptive system with higher‐order dynamic characteristics that are not present by individual members. To accurately predict the social interactions, we formulate a set of unstructured kinetic models to quantitatively capture the dynamic interactions of multiple microbial species. By introducing an interaction coefficient, we analytically derived the steady‐state solutions for the interacting species and the substrate‐depleting profile in the chemostat. We analyzed the stability of the possible coexisting states defined by competition, parasitism, amensalism, commensalism, and cooperation. Our model predicts that only parasitism, commensalism, and cooperation could lead to stable coexisting states. We also determined the optimal social interaction criteria of microbial coculture when sequential metabolic reactions are compartmentalized into two distinct species. Coupled with Luedeking–Piret and Michaelis–Menten equations, accumulation of metabolic intermediates in one species and formation of end‐product in another species could be derived and assessed. We discovered that parasitism consortia disfavor the bioconversion of intermediate to final product; and commensalism consortia could efficiently convert metabolic intermediates to final product and maintain metabolic homeostasis with a broad range of operational conditions (i.e., dilution rates); whereas cooperative consortia leads to highly nonlinear pattern of precursor accumulation and end‐product formation. The underlying dynamics and emergent properties of microbial consortia may provide critical knowledge for us to understand ecological coexisting states, engineer efficient bioconversion process, deliver effective gut therapeutics as well as elucidate probiotic‐pathogen or tumor‐host interactions in general.

KEYLINK: towards a more integrative soil representation for inclusion in ecosystem scale models. I. review and model concept

The relatively poor simulation of the below-ground processes is a severe drawback for many ecosystem models, especially when predicting responses to climate change and management. For a meaningful estimation of ecosystem production and the cycling of water, energy, nutrients and carbon, the integration of soil processes and the exchanges at the surface is crucial. It is increasingly recognized that soil biota play an important role in soil organic carbon and nutrient cycling, shaping soil structure and hydrological properties through their activity, and in water and nutrient uptake by plants through mycorrhizal processes. In this article, we review the main soil biological actors (microbiota, fauna and roots) and their effects on soil functioning. We review to what extent they have been included in soil models and propose which of them could be included in ecosystem models. We show that the model representation of the soil food web, the impact of soil ecosystem engineers on soil structure and the related effects on hydrology and soil organic matter (SOM) stabilization are key issues in improving ecosystem-scale soil representation in models. Finally, we describe a new core model concept (KEYLINK) that integrates insights from SOM models, structural models and food web models to simulate the living soil at an ecosystem scale.

Stochasticity in microbiology: managing unpredictability to reach the Sustainable Development Goals

Pure (single) cultures of microorganisms and mixed microbial communities (microbiomes) have been important for centuries in providing renewable energy, clean water and food products to human society and will continue to play a crucial role to pursue the Sustainable Development Goals. To use microorganisms effectively, microbial engineered processes require adequate control. Microbial communities are shaped by manageable deterministic processes, but also by stochastic processes, which can promote unforeseeable variations and adaptations. Here, we highlight the impact of stochasticity in single culture and microbiome engineering. First, we discuss the concepts and mechanisms of stochasticity in relation to microbial ecology of single cultures and microbiomes. Second, we discuss the consequences of stochasticity in relation to process performance and human health, which are reflected in key disadvantages and important opportunities. Third, we propose a suitable decision tool to deal with stochasticity in which monitoring of stochasticity and setting the boundaries of stochasticity by regulators are central aspects. Stochasticity may give rise to some risks, such as the presence of pathogens in microbiomes. We argue here that by taking the necessary precautions and through clever monitoring and interpretation, these risks can be mitigated.

μBialSim: Constraint-Based Dynamic Simulation of Complex Microbiomes

Microbial communities are pervasive in the natural environment, associated with many hosts, and of increasing importance in biotechnological applications. The complexity of these microbial systems makes the underlying mechanisms driving their dynamics difficult to identify. While experimental meta-OMICS techniques are routinely applied to record the inventory and activity of microbiomes over time, it remains difficult to obtain quantitative predictions based on such data. Mechanistic, quantitative mathematical modeling approaches hold the promise to both provide predictive power and shed light on cause-effect relationships driving these dynamic systems. We introduce μbialSim (pronounced "microbial sim"), a dynamic Flux-Balance-Analysis-based (dFBA) numerical simulator which is able to predict the time course in terms of composition and activity of microbiomes containing 100s of species in batch or chemostat mode. Activity of individual species is simulated by using separate FBA models which have access to a common pool of compounds, allowing for metabolite exchange. A novel augmented forward Euler method ensures numerical accuracy by temporarily reducing the time step size when compound concentrations decrease rapidly due to high compound affinities and/or the presence of many consuming species. We present three exemplary applications of μbialSim: a batch culture of a hydrogenotrophic archaeon, a syntrophic methanogenic biculture, and a 773-species human gut microbiome which exhibits a complex and dynamic pattern of metabolite exchange. Focusing on metabolite exchange as the main interaction type, μbialSim allows for the mechanistic simulation of microbiomes at their natural complexity. Simulated trajectories can be used to contextualize experimental meta-OMICS data and to derive hypotheses on cause-effect relationships driving community dynamics based on scenario simulations. μbialSim is implemented in Matlab and relies on the COBRA Toolbox or CellNetAnalyzer for FBA calculations. The source code is available under the GNU General Public License v3.0 at https://git.ufz.de/UMBSysBio/microbialsim.

Modeling population heterogeneity from microbial communities to immune response in cells

Heterogeneity is universally observed in all natural systems and across multiple scales. Understanding population heterogeneity is an intriguing and attractive topic of research in different disciplines, including microbiology and immunology. Microbes and mammalian immune cells present obviously rather different system-specific biological features. Nevertheless, as typically occurs in science, similar methods can be used to study both types of cells. This is particularly true for mathematical modeling, in which key features of a system are translated into algorithms to challenge our mechanistic understanding of the underlying biology. In this review, we first present a broad overview of the experimental developments that allowed observing heterogeneity at the single cell level. We then highlight how this "data revolution" requires the parallel advancement of algorithms and computing infrastructure for data processing and analysis, and finally present representative examples of computational models of population heterogeneity, from microbial communities to immune response in cells.

Integrating Computational Methods to Investigate the Macroecology of Microbiomes

Studies in microbiology have long been mostly restricted to small spatial scales. However, recent technological advances, such as new sequencing methodologies, have ushered an era of large-scale sequencing of environmental DNA data from multiple biomes worldwide. These global datasets can now be used to explore long standing questions of microbial ecology. New methodological approaches and concepts are being developed to study such large-scale patterns in microbial communities, resulting in new perspectives that represent a significant advances for both microbiology and macroecology. Here, we identify and review important conceptual, computational, and methodological challenges and opportunities in microbial macroecology. Specifically, we discuss the challenges of handling and analyzing large amounts of microbiome data to understand taxa distribution and co-occurrence patterns. We also discuss approaches for modeling microbial communities based on environmental data, including information on biological interactions to make full use of available Big Data. Finally, we summarize the methods presented in a general approach aimed to aid microbiologists in addressing fundamental questions in microbial macroecology, including classical propositions (such as “everything is everywhere, but the environment selects”) as well as applied ecological problems, such as those posed by human induced global environmental changes.

Modelling microbial communities using biochemical resource allocation analysis

To understand the functioning and dynamics of microbial communities is a fundamental challenge in current biology. To tackle this challenge, the construction of computational models of interacting microbes is an indispensable tool. There is, however, a large chasm between ecologically motivated descriptions of microbial growth used in many current ecosystems simulations, and detailed metabolic pathway and genome-based descriptions developed in the context of systems and synthetic biology. Here, we seek to demonstrate how resource allocation models of microbial growth offer the potential to advance ecosystem simulations and their parametrization. In particular, recent work on quantitative resource allocation allow us to formulate mechanistic models of microbial growth that are physiologically meaningful while remaining computationally tractable. These models go beyond Michaelis-Menten and Monod-type growth models, and are capable of accounting for emergent properties that underlie the remarkable plasticity of microbial growth. We outline the utility and advantages of using biochemical resource allocation models by considering a coarse-grained model of cyanobacterial growth and demonstrate how the model allows us to address specific questions of relevance for the simulation of marine microbial ecosystems, including the physiological acclimation of protein expression to different environments, the description of co-limitation by several nutrients and the differential use of alternative nutrient sources, as well as the description of metabolic diversity based on our increasing knowledge about quantitative cell physiology.

Quantifying the impact of treatment history on plasmid-mediated resistance evolution in human gut microbiota

To understand how antibiotic use affects the risk of a resistant infection, we present a computational model of the population dynamics of gut microbiota including antibiotic resistance-conferring plasmids. We then describe how this model is parameterized based on published microbiota data. Finally, we investigate how treatment history affects the prevalence of resistance among opportunistic enterobacterial pathogens. We simulate treatment histories and identify which properties of prior antibiotic exposure are most influential in determining the prevalence of resistance. We find that resistance prevalence can be predicted by 3 properties, namely the total days of drug exposure, the duration of the drug-free period after last treatment, and the center of mass of the treatment pattern. Overall this work provides a framework for capturing the role of the microbiome in the selection of antibiotic resistance and highlights the role of treatment history for the prevalence of resistance.

Emergent Subpopulation Behavior Uncovered with a Community Dynamic Metabolic Model of Escherichia coli Diauxic Growth

Microbes have adapted to greatly variable environments in order to survive both short-term perturbations and permanent changes. A classical and yet still actively studied example of adaptation to dynamic environments is the diauxic shift of Escherichia coli, in which cells grow on glucose until its exhaustion and then transition to using previously secreted acetate. Here we tested different hypotheses concerning the nature of this transition by using dynamic metabolic modeling. To reach this goal, we developed an open source modeling framework integrating dynamic models (ordinary differential equation systems) with structural models (metabolic networks) which can take into account the behavior of multiple subpopulations and smooth flux transitions between time points. We used this framework to model the diauxic shift, first with a single E. coli model whose metabolic state represents the overall population average and then with a community of two subpopulations, each growing exclusively on one carbon source (glucose or acetate). After introduction of an environment-dependent transition function that determined the balance between subpopulations, our model generated predictions that are in strong agreement with published data. Our results thus support recent experimental evidence that diauxie, rather than a coordinated metabolic shift, would be the emergent pattern of individual cells differentiating for optimal growth on different substrates. This work offers a new perspective on the use of dynamic metabolic modeling to investigate population heterogeneity dynamics. The proposed approach can easily be applied to other biological systems composed of metabolically distinct, interconverting subpopulations and could be extended to include single-cell-level stochasticity. IMPORTANCE Escherichia coli diauxie is a fundamental example of metabolic adaptation, a phenomenon that is not yet completely understood. Further insight into this process can be achieved by integrating experimental and computational modeling methods. We present a dynamic metabolic modeling approach that captures diauxie as an emergent property of subpopulation dynamics in E. coli monocultures. Without fine-tuning the parameters of the E. coli core metabolic model, we achieved good agreement with published data. Our results suggest that single-organism metabolic models can only approximate the average metabolic state of a population, therefore offering a new perspective on the use of such modeling approaches. The open source modeling framework that we provide can be applied to model general subpopulation systems in more-complex environments and can be extended to include single-cell-level stochasticity.