Advancing microbial sciences by individual-based modelling

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.

Role of biogeochemical and hydrodynamic characteristics in simulating nitrogen dynamics in river confluence

The confluence area is the link of different river systems, whose specific hydrodynamic characteristics can significantly influence mass transport and distribution, which can further make a difference to microorganism growth and biogeochemical processes. However, the specific influences of hydrodynamic characteristics in confluence on formation processes of microbial communities and the biogeochemical processes remain unclear. To this end, the present study established an indoor self-circulation confluence flume and conducted 28-day culture experiment to thoroughly investigate the characteristics of microbial communities and nitrogen dynamics in sediment of confluence area. Results illustrated that the initial homogenous microbial communities gradually emerged differences among varied hydrodynamic zones with experiment going on. Concentrations of nitrogenous materials also changed at different experiment period, NO3- concentrations peaked at day 14, and then exhibited significant downtrend. The mean NO3- concentrations decreased the most in flow separation zone, with a 62 % decrease from day 14 to day 28. A numerical model was further established following the thermodynamics of enzyme catalysis reactions to simulate nitrogen transformation rates based on abundances of associated functional genes (gene-centric model). The average relative deviation between simulated and measured N2 production rates was 32 %. To further investigate the influence of hydrodynamic characteristics on nitrogen dynamics, DamKöhler numbers were calculated as the ratio of characteristic residence time to reaction time. DamKöhler numbers were better fitted with measured N2 production rates than simulated results of gene-centric model, signifying the importance of hydrodynamic characteristics in simulating nitrogen dynamics in confluence area.

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.

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.

Integrative analysis of yeast colony growth

Yeast colonies are routinely grown on agar plates in everyday experimental settings to understand basic molecular processes, produce novel drugs, improve health, and so on. Standardized conditions ensure these colonies grow in a reproducible fashion, while in nature microbes are under a constantly changing environment. Here we combine the power of computational simulations and laboratory experiments to investigate the impact of non-standard environmental factors on colony growth. We present the developement and parameterization of a quantitative agent-based model for yeast colony growth to reproduce measurements on colony size and cell number in a colony at non-standard environmental conditions. Specifically, we establish experimental conditions that mimic the effects of humidity changes and nutrient gradients. Our results show how colony growth is affected by moisture changes, nutrient availability, and initial colony inoculation conditions. We show that initial colony spread, not initial cell number have higher impact on the final size and cell number of colonies. Parameters of the model were identified by fitting these experiments and the fitted model gives guidance to establish conditions which enable unlimited growth of yeast colonies.

ABCD of IA: A multi-scale agent-based model of T cell activation in inflammatory arthritis

Biomaterial-based agents have been demonstrated to regulate the function of immune cells in models of autoimmunity. However, the complexity of the kinetics of immune cell activation can present a challenge in optimizing the dose and frequency of administration. Here, we report a model of autoreactive T cell activation, which are key drivers in autoimmune inflammatory joint disease. The model is termed a multi-scale Agent-Based, Cell-Driven model of Inflammatory Arthritis (ABCD of IA). Using kinetic rate equations and statistical theory, ABCD of IA simulated the activation and presentation of autoantigens by dendritic cells, interactions with cognate T cells and subsequent T cell proliferation in the lymph node and IA-affected joints. The results, validated with in vivo data from the T cell driven SKG mouse model, showed that T cell proliferation strongly correlated with the T cell receptor (TCR) affinity distribution (TCR-ad), with a clear transition state from homeostasis to an inflammatory state. T cell proliferation was strongly dependent on the amount of antigen in antigenic stimulus event (ASE) at low concentrations. On the other hand, inflammation driven by Th17-inducing cytokine mediated T cell phenotype commitment was influenced by the initial level of Th17-inducing cytokines independent of the amount of arthritogenic antigen. The introduction of inhibitory artificial antigen presenting cells (iaAPCs), which locally suppress T cell activation, reduced T cell proliferation in a dose-dependent manner. The findings in this work set up a framework based on theory and modeling to simulate personalized therapeutic strategies in IA.

Cooperation between autotrophic and heterotrophic denitrifiers under low C/N ratios revealed by individual-based modelling

Denitrifying biofilms, in which autotrophic denitrifiers (AD) and heterotrophic denitrifiers (HD) coexist, play a crucial role in removing nitrate from water or wastewater. However, it is difficult to elucidate the interactions between HD and AD through sequencing-based experimental methods. Here, we developed an individual-based model to describe the interspecies dynamics and priority effects between sulfur-based AD (Thiobacillus denitrificans) and HD (Thauera phenylcarboxya) under different C/N ratios. In test I (coexistence simulation), AD and HD were initially inoculated at a ratio of 1:1. The simulation results showed excellent denitrification performance and a coaggregation pattern of denitrifiers, indicating that cooperation was the predominant interaction at a C/N ratio of 0.25 to 1.5. In test II (invasion simulation), in which only one type of denitrifier was initially inoculated and the other was added at the invasion time, denitrifiers exhibited a stratification pattern in biofilms. When HD invaded AD, the final HD abundance decreased with increasing invasion time, indicating an enhanced priority effect. When AD invaded HD, insufficient organic carbon sources weakened the priority effect by limiting the growth of HD populations. This study reveals the interaction between autotrophic and heterotrophic denitrifiers, providing guidance for optimizing wastewater treatment process.

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.

Residual cells and nutrient availability guide wound healing in bacterial biofilms

Biofilms are multicellular heterogeneous bacterial communities characterized by social-like division of labor, and remarkable robustness with respect to external stresses. Increasingly often an analogy between biofilms and arguably more complex eukaryotic tissues is being drawn. One illustrative example of where this analogy can be practically useful is the process of wound healing. While it has been extensively studied in eukaryotic tissues, the mechanism of wound healing in biofilms is virtually unexplored. Combining experiments in Bacillus subtilis bacteria, a model organism for biofilm formation, and a lattice-based theoretical model of biofilm growth, we studied how biofilms recover after macroscopic damage. We suggest that nutrient gradients and the abundance of proliferating cells are key factors augmenting wound closure. Accordingly, in the model, cell quiescence, nutrient fluxes, and biomass represented by cells and self-secreted extracellular matrix are necessary to qualitatively recapitulate the experimental results for damage repair. One of the surprising experimental findings is that residual cells, persisting in a damaged area after removal of a part of the biofilm, prominently affect the healing process. Taken together, our results outline the important roles of nutrient gradients and residual cells on biomass regrowth on macroscopic scales of the whole biofilm. The proposed combined experiment-simulation framework opens the way to further investigate the possible relation between wound healing, cell signaling and cell phenotype alternation in the local microenvironment of the wound.

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.

Simultaneous spatiotemporal transcriptomics and microscopy of Bacillus subtilis swarm development reveal cooperation across generations

Development of microbial communities is a complex multiscale phenomenon with wide-ranging biomedical and ecological implications. How biological and physical processes determine emergent spatial structures in microbial communities remains poorly understood due to a lack of simultaneous measurements of gene expression and cellular behaviour in space and time. Here we combined live-cell microscopy with a robotic arm for spatiotemporal sampling, which enabled us to simultaneously acquire phenotypic imaging data and spatiotemporal transcriptomes during Bacillus subtilis swarm development. Quantitative characterization of the spatiotemporal gene expression patterns revealed correlations with cellular and collective properties, and phenotypic subpopulations. By integrating these data with spatiotemporal metabolome measurements, we discovered a spatiotemporal cross-feeding mechanism fuelling swarm development: during their migration, earlier generations deposit metabolites which are consumed by later generations that swarm across the same location. These results highlight the importance of spatiotemporal effects during the emergence of phenotypic subpopulations and their interactions in bacterial communities.

Community composition and the environment modulate the population dynamics of type VI secretion in human gut bacteria

Understanding the relationship between the composition of the human gut microbiota and the ecological forces shaping it is of great importance; however, knowledge of the biogeographical and ecological relationships between physically interacting taxa is limited. Interbacterial antagonism may play an important role in gut community dynamics, yet the conditions under which antagonistic behaviour is favoured or disfavoured by selection in the gut are not well understood. Here, using genomics, we show that a species-specific type VI secretion system (T6SS) repeatedly acquires inactivating mutations in Bacteroides fragilis in the human gut. This result implies a fitness cost to the T6SS, but we could not identify laboratory conditions under which such a cost manifests. Strikingly, experiments in mice illustrate that the T6SS can be favoured or disfavoured in the gut depending on the strains and species in the surrounding community and their susceptibility to T6SS antagonism. We use ecological modelling to explore the conditions that could underlie these results and find that community spatial structure modulates interaction patterns among bacteria, thereby modulating the costs and benefits of T6SS activity. Our findings point towards new integrative models for interrogating the evolutionary dynamics of type VI secretion and other modes of antagonistic interaction in microbiomes.

In Situ and Individual-Based Analysis of the Influence of Polystyrene Microplastics on Escherichia coli Conjugative Gene Transfer at the Single-Cell Level

The impact of microplastic particles of micro- and nanometer sizes on microbial horizontal gene transfer (HGT) remains a controversial topic. Existing studies rely on traditional approaches, which analyze population behavior, leading to conflicting conclusions and a limited understanding. The present study addressed these limitations by employing a novel microfluidic chamber system for in situ visualization and precise quantification of the effects of different concentrations of polystyrene (PS) microbeads on microbial HGT at the single-cell level. The statistical analysis indicated no significant difference in the division times of both the donor and recipient bacteria across different PS microbead concentrations. However, as the concentration of PS microbeads increased from 0 to 2000 mg L-1, the average conjugation frequency of Escherichia coli decreased from 0.028 ± 0.015 to 0.004 ± 0.003. Our observations from the microfluidic experiments revealed that 500 nm PS microbeads created a barrier effect on bacterial conjugative transfer. The presence of microbeads resulted in reduced contact and interaction between the donor and recipient strains, thereby causing a decrease in the conjugation transfer frequency. These findings were validated by an individual-based modeling framework parameterized by the data from the individual-level microfluidic experiments. Overall, this study offers a fresh perspective and strategy for investigating the risks associated with the dissemination of antibiotic resistance genes related to microplastics.

Advances in automated and reactive flow cytometry for synthetic biotechnology

Automated flow cytometry (FC) has been initially considered for bioprocess monitoring and optimization. More recently, new physical and software interfaces have been made available, facilitating the access to this technology for labs and industries. It also comes with new capabilities, such as being able to act on the cultivation conditions based on population data. This approach, known as reactive FC, extended the range of applications of automated FC to bioprocess control and the stabilization of cocultures, but also to the broad field of synthetic and systems biology for the characterization of gene circuits. However, several issues must be addressed before automated and reactive FC can be considered standard and modular technologies.

Corrected and Republished from: "Understanding Lactobacillus paracasei and Streptococcus oralis Biofilm Interactions through Agent-Based Modeling"

As common commensals residing on mucosal tissues, Lactobacillus species are known to promote health, while some Streptococcus species act to enhance the pathogenicity of other organisms in those environments. In this study we used a combination of in vitro imaging of live biofilms and computational modeling to explore biofilm interactions between Streptococcus oralis, an accessory pathogen in oral candidiasis, and Lactobacillus paracasei, an organism with known probiotic properties. A computational agent-based model was created where the two species interact only by competing for space, oxygen, and glucose. Quantification of bacterial growth in live biofilms indicated that S. oralis biomass and cell numbers were much lower than predicted by the model. Two subsequent models were then created to examine more complex interactions between these species, one where L. paracasei secretes a surfactant and another where L. paracasei secretes an inhibitor of S. oralis growth. We observed that the growth of S. oralis could be affected by both mechanisms. Further biofilm experiments support the hypothesis that L. paracasei may secrete an inhibitor of S. oralis growth, although they do not exclude that a surfactant could also be involved. This contribution shows how agent-based modeling and experiments can be used in synergy to address multiple-species biofilm interactions, with important roles in mucosal health and disease. IMPORTANCE We previously discovered a role of the oral commensal Streptococcus oralis as an accessory pathogen. S. oralis increases the virulence of Candida albicans infections in murine oral candidiasis and epithelial cell models through mechanisms which promote the formation of tissue-damaging biofilms. Lactobacillus species have known inhibitory effects on biofilm formation of many microbes, including Streptococcus species. Agent-based modeling has great advantages as a means of exploring multifaceted relationships between organisms in complex environments such as biofilms. Here, we used an iterative collaborative process between experimentation and modeling to reveal aspects of the mostly unexplored relationship between S. oralis and L. paracasei in biofilm growth. The inhibitory nature of L. paracasei on S. oralis in biofilms may be exploited as a means of preventing or alleviating mucosal fungal infections.

Community composition and the environment modulate the population dynamics of type VI secretion in human gut bacteria

Understanding the relationship between the composition of the human gut microbiota and the ecological forces shaping it is of high importance as progress towards therapeutic modulation of the microbiota advances. However, given the inaccessibility of the gastrointestinal tract, our knowledge of the biogeographical and ecological relationships between physically interacting taxa has been limited to date. It has been suggested that interbacterial antagonism plays an important role in gut community dynamics, but in practice the conditions under which antagonistic behavior is favored or disfavored by selection in the gut environment are not well known. Here, using phylogenomics of bacterial isolate genomes and analysis of infant and adult fecal metagenomes, we show that the contact-dependent type VI secretion system (T6SS) is repeatedly lost from the genomes of Bacteroides fragilis in adults compare to infants. Although this result implies a significant fitness cost to the T6SS, but we could not identify in vitro conditions under which such a cost manifests. Strikingly, however, experiments in mice illustrated that the B. fragilis T6SS can be favored or disfavored in the gut environment, depending on the strains and species in the surrounding community and their susceptibility to T6SS antagonism. We use a variety of ecological modeling techniques to explore the possible local community structuring conditions that could underlie the results of our larger scale phylogenomic and mouse gut experimental approaches. The models illustrate robustly that the pattern of local community structuring in space can modulate the extent of interactions between T6SS-producing, sensitive, and resistant bacteria, which in turn control the balance of fitness costs and benefits of performing contact-dependent antagonistic behavior. Taken together, our genomic analyses, in vivo studies, and ecological theory point toward new integrative models for interrogating the evolutionary dynamics of type VI secretion and other predominant modes of antagonistic interaction in diverse microbiomes.

Controlling the human microbiome

We coexist with a vast number of microbes that live in and on our bodies. Those microbes and their genes are collectively known as the human microbiome, which plays important roles in human physiology and diseases. We have acquired extensive knowledge of the organismal compositions and metabolic functions of the human microbiome. However, the ultimate proof of our understanding of the human microbiome is reflected in our ability to manipulate it for health benefits. To facilitate the rational design of microbiome-based therapies, there are many fundamental questions to be addressed at the systems level. Indeed, we need a deep understanding of the ecological dynamics associated with such a complex ecosystem before we rationally design control strategies. In light of this, this review discusses progress from various fields, e.g., community ecology, network science, and control theory, that are helping us make progress toward the ultimate goal of controlling the human microbiome.

Model-aided targeted volatile fatty acid production from food waste using a defined co-culture microbial community

The production of volatile fatty acids (VFA) is gaining momentum due to their central role in the emerging carboxylate platform. Particularly, the production of the longest VFA (from butyrate to caproate) is desired due to their increased economic value and easier downstream processing. While the use of undefined microbial cultures is usually preferred with organic waste streams, the use of defined microbial co-culture processes could tackle some of their drawbacks such as poor control over the process outcome, which often leads to low selectivity for the desired products. However, the extensive experimentation needed to design a co-culture system hinders the use of this technology. In this work, a workflow based on the combined use of mathematical models and wet experimentation is proposed to accelerate the design of novel bioprocesses. In particular, a co-culture consisting of Pediococcus pentosaceus and Megaphaera cerevisiae is used to target the production of high-value odd- and even‑carbon VFA. An unstructured kinetic model was developed, calibrated and used to design experiments with the goal of increasing the selectivity for the desired VFA, which were experimentally validated. In the case of even‑carbon VFA, the experimental validation showed an increase of 38 % in caproate yield and, in the case of enhanced odd‑carbon VFA experiments, the yield of butyrate and caproate diminished by 62 % and 94 %, respectively, while propionate became one of the main end products and valerate yield value increased from 0.007 to 0.085 g_valearte per g_{consumed sugar.} The workflow followed in this work proved to be a sound tool for bioprocess design due to its capacity to explore and design new experiments in silico in a fast way and ability to quickly adapt to new scenarios.

Computational Design of Antimicrobial Active Surfaces via Automated Bayesian Optimization

Biofilms pose significant problems for engineers in diverse fields, such as marine science, bioenergy, and biomedicine, where effective biofilm control is a long-term goal. The adhesion and surface mechanics of biofilms play crucial roles in generating and removing biofilm. Designing customized nanosurfaces with different surface topologies can alter the adhesive properties to remove biofilms more easily and greatly improve long-term biofilm control. To rapidly design such topologies, we employ individual-based modeling and Bayesian optimization to automate the design process and generate different active surfaces for effective biofilm removal. Our framework successfully generated optimized functional nanosurfaces for improved biofilm removal through applied shear and vibration. Densely distributed short pillar topography is the optimal geometry to prevent biofilm formation. Under fluidic shearing, the optimal topography is to sparsely distribute tall, slim, pillar-like structures. When subjected to either vertical or lateral vibrations, thick trapezoidal cones are found to be optimal. Optimizing the vibrational loading indicates a small vibration magnitude with relatively low frequencies is more efficient in removing biofilm. Our results provide insights into various engineering fields that require surface-mediated biofilm control. Our framework can also be applied to more general materials design and optimization.

Novel Ground-Up 3D Multicellular Simulators for Synthetic Biology CAD Integrating Stochastic Gillespie Simulations Benchmarked with Topologically Variable SBML Models

The elevation of Synthetic Biology from single cells to multicellular simulations would be a significant scale-up. The spatiotemporal behavior of cellular populations has the potential to be prototyped in silico for computer assisted design through ergonomic interfaces. Such a platform would have great practical potential across medicine, industry, research, education and accessible archiving in bioinformatics. Existing Synthetic Biology CAD systems are considered limited regarding population level behavior, and this work explored the in silico challenges posed from biological and computational perspectives. Retaining the connection to Synthetic Biology CAD, an extension of the Infobiotics Workbench Suite was considered, with potential for the integration of genetic regulatory models and/or chemical reaction networks through Next Generation Stochastic Simulator (NGSS) Gillespie algorithms. These were executed using SBML models generated by in-house SBML-Constructor over numerous topologies and benchmarked in association with multicellular simulation layers. Regarding multicellularity, two ground-up multicellular solutions were developed, including the use of Unreal Engine 4 contrasted with CPU multithreading and Blender visualization, resulting in a comparison of real-time versus batch-processed simulations. In conclusion, high-performance computing and client-server architectures could be considered for future works, along with the inclusion of numerous biologically and physically informed features, whilst still pursuing ergonomic solutions.

Environmental and ecological controls of the spatial distribution of microbial populations in aggregates

In microbial communities, the ecological interactions between species of different populations are responsible for the spatial distributions observed in aggregates (granules, biofilms or flocs). To explore the underlying mechanisms that control these processes, we have developed a mathematical modelling framework able to describe, label and quantify defined spatial structures that arise from microbial and environmental interactions in communities. An artificial system of three populations collaborating or competing in an aggregate is simulated using individual-based modelling under different environmental conditions. In this study, neutralism, competition, commensalism and concurrence of commensalism and competition have been considered. We were able to identify interspecific segregation of communities that appears in competitive environments (columned stratification), and a layered distribution of populations that emerges in commensal (layered stratification). When different ecological interactions were considered in the same aggregate, the resultant spatial distribution was identified as the one controlled by the most limiting substrate. A theoretical modulus was defined, with which we were able to quantify the effect of environmental conditions and ecological interactions to predict the most probable spatial distribution. The specific microbial patterns observed in our results allowed us to identify the optimal spatial organizations for bacteria to thrive when building a microbial community and how this permitted co-existence of populations at different growth rates. Our model reveals that although ecological relationships between different species dictate the distribution of bacteria, the environment controls the final spatial distribution of the community.

Guided by the principles of microbiome engineering: Accomplishments and perspectives for environmental use

Although the accomplishments of microbiome engineering highlight its significance for the targeted manipulation of microbial communities, knowledge and technical gaps still limit the applications of microbiome engineering in biotechnology, especially for environmental use. Addressing the environmental challenges of refractory pollutants and fluctuating environmental conditions requires an adequate understanding of the theoretical achievements and practical applications of microbiome engineering. Here, we review recent cutting-edge studies on microbiome engineering strategies and their classical applications in bioremediation. Moreover, a framework is summarized for combining both top-down and bottom-up approaches in microbiome engineering toward improved applications. A strategy to engineer microbiomes for environmental use, which avoids the build-up of toxic intermediates that pose a risk to human health, is suggested. We anticipate that the highlighted framework and strategy will be beneficial for engineering microbiomes to address difficult environmental challenges such as degrading multiple refractory pollutants and sustain the performance of engineered microbiomes in situ with indigenous microorganisms under fluctuating conditions.

Agent-Based Modeling of Microbial Communities

Microbial communities are complex living systems that populate the planet with diverse functions and are increasingly harnessed for practical human needs. To deepen the fundamental understanding of their organization and functioning as well as to facilitate their engineering for applications, mathematical modeling has played an increasingly important role. Agent-based models represent a class of powerful quantitative frameworks for investigating microbial communities because of their individualistic nature in describing cells, mechanistic characterization of molecular and cellular processes, and intrinsic ability to produce emergent system properties. This review presents a comprehensive overview of recent advances in agent-based modeling of microbial communities. It surveys the state-of-the-art algorithms employed to simulate intracellular biomolecular events, single-cell behaviors, intercellular interactions, and interactions between cells and their environments that collectively serve as the driving forces of community behaviors. It also highlights three lines of applications of agent-based modeling, namely, the elucidation of microbial range expansion and colony ecology, the design of synthetic gene circuits and microbial populations for desired behaviors, and the characterization of biofilm formation and dispersal. The review concludes with a discussion of existing challenges, including the computational cost of the modeling, and potential mitigation strategies.

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.

Do microbes have a memory? History-dependent behavior in the adaptation to variable environments

Microbes are constantly confronted with changes and challenges in their environment. A proper response to these environmental cues is needed for optimal cellular functioning and fitness. Interestingly, past exposure to environmental cues can accelerate or boost the response when this condition returns, even in daughter cells that have not directly encountered the initial cue. Moreover, this behavior is mostly epigenetic and often goes hand in hand with strong heterogeneity in the strength and speed of the response between isogenic cells of the same population, which might function as a bet-hedging strategy. In this review, we discuss examples of history-dependent behavior (HDB) or “memory,” with a specific focus on HDB in fluctuating environments. In most examples discussed, the lag time before the response to an environmental change is used as an experimentally measurable proxy for HDB. We highlight different mechanisms already implicated in HDB, and by using HDB in fluctuating carbon conditions as a case study, we showcase how the metabolic state of a cell can be a key determining factor for HDB. Finally, we consider possible evolutionary causes and consequences of such HDB.

Impacts of atmospheric stilling and climate warming on cyanobacterial blooms: An individual-based modelling approach

Harmful algal blooms of the freshwater cyanobacteria genus Microcystis are a global problem and are expected to intensify with climate change. In studies of climate change impacts on Microcystis blooms, atmospheric stilling has not been considered. Stilling is expected to occur in some regions of the world with climate warming, and it will affect lake stratification regimes. We tested if stilling could affect water column Microcystis distributions using a novel individual-based model (IBM). Using the IBM coupled to a three-dimensional hydrodynamic model, we assessed responses of colonial Microcystis biomass to wind speed decrease and air temperature increase projected under a future climate. The IBM altered Microcystis colony size using relationships with turbulence from the literature, and included light, temperature, and nutrient effects on Microcystis growth using input data from a shallow urban lake. The model results show that dynamic variations in colony size are critical for accurate prediction of cyanobacterial bloom development and decay. Colony size (mean and variability) increased more than six-fold for a 20% decrease in wind speed compared with a 2 °C increase in air temperature. Our results suggest that atmospheric stilling needs to be included in projections of changes in the frequency, distribution and magnitude of blooms of buoyant, colony-forming cyanobacteria under climate change.

Ecological modelling approaches for predicting emergent properties in microbial communities

Recent studies have brought forward the critical role of emergent properties in shaping microbial communities and the ecosystems of which they are a part. Emergent properties-patterns or functions that cannot be deduced linearly from the properties of the constituent parts-underlie important ecological characteristics such as resilience, niche expansion and spatial self-organization. While it is clear that emergent properties are a consequence of interactions within the community, their non-linear nature makes mathematical modelling imperative for establishing the quantitative link between community structure and function. As the need for conservation and rational modulation of microbial ecosystems is increasingly apparent, so is the consideration of the benefits and limitations of the approaches to model emergent properties. Here we review ecosystem modelling approaches from the viewpoint of emergent properties. We consider the scope, advantages and limitations of Lotka-Volterra, consumer-resource, trait-based, individual-based and genome-scale metabolic models. Future efforts in this research area would benefit from capitalizing on the complementarity between these approaches towards enabling rational modulation of complex microbial ecosystems.

A Self-Controlled and Self-Healing Model of Bacterial Cells

A new kind of self-assembly model, morphogenetic (M) systems, assembles spatial units into larger structures through local interactions of simpler components and enables discovery of new principles for cellular membrane assembly, development, and its interface function. The model is based on interactions among three kinds of constitutive objects such as tiles and protein-like elements in discrete time and continuous 3D space. It was motivated by achieving a balance between three conflicting goals: biological, physical-chemical, and computational realism. A recent example is a unified model of morphogenesis of a single biological cell, its membrane and cytoskeleton formation, and finally, its self-reproduction. Here, a family of dynamic M systems (Mbac) is described with similar characteristics, modeling the process of bacterial cell formation and division that exhibits bacterial behaviors of living cells at the macro-level (including cell growth that is self-controlled and sensitive to the presence/absence of nutrients transported through membranes), as well as self-healing properties. Remarkably, it consists of only 20 or so developmental rules. Furthermore, since the model exhibits membrane formation and septic mitosis, it affords more rigorous definitions of concepts such as injury and self-healing that enable quantitative analyses of these kinds of properties. Mbac shows that self-assembly and interactions of living organisms with their environments and membrane interfaces are critical for self-healing, and that these properties can be defined and quantified more rigorously and precisely, despite their complexity.

Substrate availability and toxicity shape the structure of microbial communities engaged in metabolic division of labor

Metabolic division of labor (MDOL) represents a widespread natural phenomenon, whereby a complex metabolic pathway is shared between different strains within a community in a mutually beneficial manner. However, little is known about how the composition of such a microbial community is regulated. We hypothesized that when degradation of an organic compound is carried out via MDOL, the concentration and toxicity of the substrate modulate the benefit allocation between the two microbial populations, thus affecting the structure of this community. We tested this hypothesis by combining modeling with experiments using a synthetic consortium. Our modeling analysis suggests that the proportion of the population executing the first metabolic step can be simply estimated by Monod-like formulas governed by substrate concentration and toxicity. Our model and the proposed formula were able to quantitatively predict the structure of our synthetic consortium. Further analysis demonstrates that our rule is also applicable in estimating community structures in spatially structured environments. Together, our work clearly demonstrates that the structure of MDOL communities can be quantitatively predicted using available information on environmental factors, thus providing novel insights into how to manage artificial microbial systems for the wide application of the bioindustry.

Models predict planned phosphorus load reduction will make Lake Erie more toxic

Harmful cyanobacteria are a global environmental problem, yet we lack actionable understanding of toxigenic versus nontoxigenic strain ecology and toxin production. We performed a large-scale meta-analysis including 103 papers and used it to develop a mechanistic, agent-based model of Microcystis growth and microcystin production. Simulations for Lake Erie suggest that the observed toxigenic-to-nontoxigenic strain succession during the 2014 Toledo drinking water crisis was controlled by different cellular oxidative stress mitigation strategies (protection by microcystin versus degradation by enzymes) and the different susceptibility of those mechanisms to nitrogen limitation. This model, as well as a simpler empirical one, predicts that the planned phosphorus load reduction will lower biomass but make nitrogen and light more available, which will increase toxin production, favor toxigenic cells, and increase toxin concentrations.

Leveraging the model-experiment loop: Examples from cellular slime mold chemotaxis

Interplay between models and experimental data advances discovery and understanding in biology, particularly when models generate predictions that allow well-designed experiments to distinguish between alternative mechanisms. To illustrate how this feedback between models and experiments can lead to key insights into biological mechanisms, we explore three examples from cellular slime mold chemotaxis. These examples include studies that identified chemotaxis as the primary mechanism behind slime mold aggregation, discovered that cells likely measure chemoattractant gradients by sensing concentration differences across cell length, and tested the role of cell-associated chemoattractant degradation in shaping chemotactic fields. Although each study used a different model class appropriate to their hypotheses - qualitative, mathematical, or simulation-based - these examples all highlight the utility of modeling to formalize assumptions and generate testable predictions. A central element of this framework is the iterative use of models and experiments, specifically: matching experimental designs to the models, revising models based on mismatches with experimental data, and validating critical model assumptions and predictions with experiments. We advocate for continued use of this interplay between models and experiments to advance biological discovery.

Roadmap from Microbial Communities to Individuality Modeling for Anaerobic Digestion of Sewage Sludge

Biological models describing anaerobic digestion (AD) of sewage sludge have been widely applied to test various control and operation strategies. Anaerobic digestion model 1 (ADM1) provides a generic platform that includes the main processes of AD, excluding homoacetogenesis and the microbial structure. Homoacetogenic bacteria have been identified as important competitors for hydrogen consumption and acetate production. Although recent advances in meta-omics techniques have improved our characterization of AD microbial communities, conventional models treat functional groups as homogeneous and overlook the physiology and behavior of microbial individuality, limiting insights into mechanisms governing process performance. A novel microbial individuality model (MIM) that integrates kinetics, energetics, and agent-based modeling to describe a microbiome's behavior as heterogenic populations, including homoacetogenesis, was developed. The MIM was validated with two datasets from previous studies through daily biogas production, methane content, compound concentrations, and microbial relative abundance changes. The MIM identified the emergence of Methanosaeta at low concentrations of acetate. Moreover, this simulation supports experimental studies confirming that the overlooked homoacetogenesis is an important hydrogen sink in AD. Validated MIMs are expected to provide insights into syntrophic and competitive interactions among microbiomes in AD systems while testing different operational parameters in a virtual environment. The MIM offers a methodological framework to other biological treatment systems and their microbial community dynamics.

Agent-based modelling of iron cycling bacteria provides a framework for testing alternative environmental conditions and modes of action

Iron-reducing and iron-oxidizing bacteria are of interest in a variety of environmental and industrial applications. Such bacteria often co-occur at oxic-anoxic gradients in aquatic and terrestrial habitats. In this paper, we present the first computational agent-based model of microbial iron cycling, between the anaerobic ferric iron (Fe³⁺)-reducing bacteria Shewanella spp. and the microaerophilic ferrous iron (Fe²⁺)-oxidizing bacteria Sideroxydans spp. By including the key processes of reduction/oxidation, movement, adhesion, Fe²⁺-equilibration and nanoparticle formation, we derive a core model which enables hypothesis testing and prediction for different environmental conditions including temporal cycles of oxic and anoxic conditions. We compared (i) combinations of different Fe³⁺-reducing/Fe²⁺-oxidizing modes of action of the bacteria and (ii) system behaviour for different pH values. We predicted that the beneficial effect of a high number of iron-nanoparticles on the total Fe³⁺ reduction rate of the system is not only due to the faster reduction of these iron-nanoparticles, but also to the nanoparticles' additional capacity to bind Fe²⁺ on their surfaces. Efficient iron-nanoparticle reduction is confined to pH around 6, being twice as high than at pH 7, whereas at pH 5 negligible reduction takes place. Furthermore, in accordance with experimental evidence our model showed that shorter oxic/anoxic periods exhibit a faster increase of total Fe³⁺ reduction rate than longer periods.

McComedy: A user-friendly tool for next-generation individual-based modeling of microbial consumer-resource systems

Individual-based modeling is widely applied to investigate the ecological mechanisms driving microbial community dynamics. In such models, the population or community dynamics emerge from the behavior and interplay of individual entities, which are simulated according to a predefined set of rules. If the rules that govern the behavior of individuals are based on generic and mechanistically sound principles, the models are referred to as next-generation individual-based models. These models perform particularly well in recapitulating actual ecological dynamics. However, implementation of such models is time-consuming and requires proficiency in programming or in using specific software, which likely hinders a broader application of this powerful method. Here we present McComedy, a modeling tool designed to facilitate the development of next-generation individual-based models of microbial consumer-resource systems. This tool allows flexibly combining pre-implemented building blocks that represent physical and biological processes. The ability of McComedy to capture the essential dynamics of microbial consumer-resource systems is demonstrated by reproducing and furthermore adding to the results of two distinct studies from the literature. With this article, we provide a versatile tool for developing next-generation individual-based models that can foster understanding of microbial ecology in both research and education.

Microorganisms such as bacteria and fungi can be found in virtually any natural environment. To better understand the ecology of these microorganisms–which is important for several research fields including medicine, biotechnology, and conservation biology–researchers often use computer models to simulate and predict the behavior of microbial communities. Commonly, a particular technique called individual-based modeling is used to generate structurally realistic models of these communities by explicitly simulating each individual microorganism. Here we developed a tool called McComedy that helps researchers applying individual-based modeling efficiently without having to program low-level processes, thus allowing them to focus on their actual research questions. To test whether McComedy is not only convenient to use but also generates meaningful models, we used it to reproduce previously reported findings of two other research groups. Given that our results could well recapitulate and furthermore extend the original findings, we are confident that McComedy is a powerful and versatile tool that can help to address fundamental questions in microbial ecology.

Kinetics and mathematical model of killer/sensitive interaction under different physicochemical conditions of must/wine: a study from a biological point of view

Fermentation of grape must to wine is carried out by a complex microbial mixture, which also involves spoilage yeasts of wine. The latter yeasts produce organoleptic changes that cause significant economic losses to the wine industry. SO₂ is traditionally used to control this spoilage populations, but because of its harmful effects on human health, biocontrol has emerged as an alternative treatment. Although studies have been carried out to select biocontroller yeasts and examine their underlying mechanisms of action, reports on their application have not been published yet. To better understand the interaction and the successful application of biocontrol, the use of mathematical models, among other methods, is important, as they facilitate the prediction of success or failure of the antagonist. The objective of the present study was to use an existing mathematical model to obtain information about the yeast's interaction assayed and to validate its predictive use under different physicochemical conditions during the wine fermentation, and eventually predict biocontrol kinetics. The mathematical model was applied to the fermentation conditions and provided information on the kinetic parameters of the biocontrol interaction and allowed interpretations about other parameters. The model was applied in the different physicochemical conditions for the biocontrol and did not fit correctly to experimental data, and therefore an improvement was proposed which was successful and presented new hypotheses.

Understanding Lactobacillus paracasei and Streptococcus oralis Biofilm Interactions through Agent-Based Modeling

As common commensals residing on mucosal tissues, Lactobacillus species are known to promote health, while some Streptococcus species act to enhance the pathogenicity of other organisms in those environments. In this study, we used a combination of in vitro imaging of live biofilms and computational modeling to explore biofilm interactions between Streptococcus oralis, an accessory pathogen in oral candidiasis, and Lactobacillus paracasei, an organism with known probiotic properties. A computational agent-based model was created where the two species interact only by competing for space, oxygen and glucose. Quantification of bacterial growth in live biofilms indicated that S. oralis biomass and cell numbers were much lower than predicted by the model. Two subsequent models were then created to examine more complex interactions between these species, one where L. paracasei secretes a surfactant, and another where L. paracasei secretes an inhibitor of S. oralis growth. We observed that the growth of S. oralis could be affected by both mechanisms. Further biofilm experiments support the hypothesis that L. paracasei may secrete an inhibitor of S. oralis growth, although they do not exclude that a surfactant could also be involved. This contribution shows how agent-based modeling and experiments can be used in synergy to address multiple species biofilm interactions, with important roles in mucosal health and disease. IMPORTANCE We previously discovered a role of the oral commensal Streptococcus oralis as an accessory pathogen. S. oralis increases the virulence of Candida albicans infections in murine oral candidiasis and epithelial cell models through mechanisms which promote the formation of tissue-damaging biofilms. Lactobacillus species have known inhibitory effects on biofilm formation of many microbes, including Streptococcus species. Agent-based modeling has great advantages as a means of exploring multifaceted relationships between organisms in complex environments such as biofilms. Here, we used an iterative collaborative process between experimentation and modeling to reveal aspects of the mostly unexplored relationship between S. oralis and L. paracasei in biofilm growth. The inhibitory nature of L. paracasei on S. oralis in biofilms may be exploited as a means of preventing or alleviating mucosal fungal infections.

Derivation of continuum models from discrete models of mechanical forces in cell populations

In certain discrete models of populations of biological cells, the mechanical forces between the cells are center based or vertex based on the microscopic level where each cell is individually represented. The cells are circular or spherical in a center based model and polygonal or polyhedral in a vertex based model. On a higher, macroscopic level, the time evolution of the density of the cells is described by partial differential equations (PDEs). We derive relations between the modelling on the micro and macro levels in one, two, and three dimensions by regarding the micro model as a discretization of a PDE for conservation of mass on the macro level. The forces in the micro model correspond on the macro level to a gradient of the pressure scaled by quantities depending on the cell geometry. The two levels of modelling are compared in numerical experiments in one and two dimensions.

Symbiosis in Digital Evolution: Past, Present, and Future

Symbiosis, the living together of unlike organisms as symbionts, is ubiquitous in the natural world. Symbioses occur within and across all scales of life, from microbial to macro-faunal systems. Further, the interactions between symbionts are multimodal in both strength and type, can span from parasitic to mutualistic within one partnership, and persist over generations. Studying the ecological and evolutionary dynamics of symbiosis in natural or laboratory systems poses a wide range of challenges, including the long time scales at which symbioses evolve de novo, the limited capacity to experimentally control symbiotic interactions, the weak resolution at which we can quantify interactions, and the idiosyncrasies of current model systems. These issues are especially challenging when seeking to understand the ecological effects and evolutionary pressures on and of a symbiosis, such as how a symbiosis may shift between parasitic and mutualistic modes and how that shift impacts the dynamics of the partner population. In digital evolution, populations of computational organisms compete, mutate, and evolve in a virtual environment. Digital evolution features perfect data tracking and allows for experimental manipulations that are impractical or impossible in natural systems. Furthermore, modern computational power allows experimenters to observe thousands of generations of evolution in minutes (as opposed to several months or years), which greatly expands the range of possible studies. As such, digital evolution is poised to become a keystone technique in our methodological repertoire for studying the ecological and evolutionary dynamics of symbioses. Here, we review how digital evolution has been used to study symbiosis, and we propose a series of open questions that digital evolution is well-positioned to answer.

Mathematical Models of Plasmid Population Dynamics

With plasmid-mediated antibiotic resistance thriving and threatening to become a serious public health problem, it is paramount to increase our understanding of the forces that enable the spread and maintenance of drug resistance genes encoded in mobile genetic elements. The relevance of plasmids as vehicles for the dissemination of antibiotic resistance genes, in addition to the extensive use of plasmid-derived vectors for biotechnological and industrial purposes, has promoted the in-depth study of the molecular mechanisms controlling multiple aspects of a plasmids' life cycle. This body of experimental work has been paralleled by the development of a wealth of mathematical models aimed at understanding the interplay between transmission, replication, and segregation, as well as their consequences in the ecological and evolutionary dynamics of plasmid-bearing bacterial populations. In this review, we discuss theoretical models of plasmid dynamics that span from the molecular mechanisms of plasmid partition and copy-number control occurring at a cellular level, to their consequences in the population dynamics of complex microbial communities. We conclude by discussing future directions for this exciting research topic.

Individual-based modelling of cyanobacteria blooms: Physical and physiological processes

Lakes and reservoirs throughout the world are increasingly adversely affected by cyanobacterial harmful algal blooms (CyanoHABs). The development and spatiotemporal distributions of blooms are governed by complex physical mixing and transport processes that interact with physiological processes affecting the growth and loss of bloom-forming species. Individual-based models (IBMs) can provide a valuable tool for exploring and integrating some of these processes. Here we contend that the advantages of IBMs have not been fully exploited. The main reasons for the lack of progress in mainstreaming IBMs in numerical modelling are their complexity and high computational demand. In this review, we identify gaps and challenges in the use of IBMs for modelling CyanoHABs and provide an overview of the processes that should be considered for simulating the spatial and temporal distributions of cyanobacteria. Notably, important processes affecting cyanobacteria distributions, in particular their vertical passive movement, have not been considered in many existing lake ecosystem models. We identify the following research gaps that should be addressed in future studies that use IBMs: 1) effects of vertical movement and physiological processes relevant to cyanobacteria growth and accumulations, 2) effects and feedbacks of CyanoHABs on their environment; 3) inter and intra-specific competition of cyanobacteria species for nutrients and light; 4) use of high resolved temporal-spatial data for calibration and verification targets for IBMs; and 5) climate change impacts on the frequency, intensity and duration of CyanoHABs. IBMs are well adapted to incorporate these processes and should be considered as the next generation of models for simulating CyanoHABs.

Application of Machine Learning Techniques to an Agent-Based Model of Pantoea

Agent-based modeling (ABM) is a powerful simulation technique which describes a complex dynamic system based on its interacting constituent entities. While the flexibility of ABM enables broad application, the complexity of real-world models demands intensive computing resources and computational time; however, a metamodel may be constructed to gain insight at less computational expense. Here, we developed a model in NetLogo to describe the growth of a microbial population consisting of Pantoea. We applied 13 parameters that defined the model and actively changed seven of the parameters to modulate the evolution of the population curve in response to these changes. We efficiently performed more than 3,000 simulations using a Python wrapper, NL4Py. Upon evaluation of the correlation between the active parameters and outputs by random forest regression, we found that the parameters which define the depth of medium and glucose concentration affect the population curves significantly. Subsequently, we constructed a metamodel, a dense neural network, to predict the simulation outputs from the active parameters and found that it achieves high prediction accuracy, reaching an R² coefficient of determination value up to 0.92. Our approach of using a combination of ABM with random forest regression and neural network reduces the number of required ABM simulations. The simplified and refined metamodels may provide insights into the complex dynamic system before their transition to more sophisticated models that run on high-performance computing systems. The ultimate goal is to build a bridge between simulation and experiment, allowing model validation by comparing the simulated data to experimental data in microbiology.

A Review on Occurrence and Spread of Antibiotic Resistance in Wastewaters and in Wastewater Treatment Plants: Mechanisms and Perspectives

This paper reviews current knowledge on sources, spread and removal mechanisms of antibiotic resistance genes (ARGs) in microbial communities of wastewaters, treatment plants and downstream recipients. Antibiotic is the most important tool to cure bacterial infections in humans and animals. The over- and misuse of antibiotics have played a major role in the development, spread, and prevalence of antibiotic resistance (AR) in the microbiomes of humans and animals, and microbial ecosystems worldwide. AR can be transferred and spread amongst bacteria via intra- and interspecies horizontal gene transfer (HGT). Wastewater treatment plants (WWTPs) receive wastewater containing an enormous variety of pollutants, including antibiotics, and chemicals from different sources. They contain large and diverse communities of microorganisms and provide a favorable environment for the spread and reproduction of AR. Existing WWTPs are not designed to remove micropollutants, antibiotic resistant bacteria (ARB) and ARGs, which therefore remain present in the effluent. Studies have shown that raw and treated wastewaters carry a higher amount of ARB in comparison to surface water, and such reports have led to further studies on more advanced treatment processes. This review summarizes what is known about AR removal efficiencies of different wastewater treatment methods, and it shows the variations among different methods. Results vary, but the trend is that conventional activated sludge treatment, with aerobic and/or anaerobic reactors alone or in series, followed by advanced post treatment methods like UV, ozonation, and oxidation removes considerably more ARGs and ARB than activated sludge treatment alone. In addition to AR levels in treated wastewater, it examines AR levels in biosolids, settled by-product from wastewater treatment, and discusses AR removal efficiency of different biosolids treatment procedures. Finally, it puts forward key-points and suggestions for dealing with and preventing further increase of AR in WWTPs and other aquatic environments, together with a discussion on the use of mathematical models to quantify and simulate the spread of ARGs in WWTPs. Mathematical models already play a role in the analysis and development of WWTPs, but they do not consider AR and challenges remain before models can be used to reliably study the dynamics and reduction of AR in such systems.

Assessment of Bacillus subtilis Plasmid pLS20 Conjugation in the Absence of Quorum Sensing Repression

Bacillus subtilis conjugative plasmid pLS20 uses a quorum-sensing mechanism to control expression levels of its conjugation genes, involving the repressor Rco_pLS20, the anti-repressor Rap_pLS20, and the signaling peptide Phr*_pLS20. In previous studies, artificial overexpression of rap_pLS20 in the donor cells was shown to enhance conjugation efficiency. However, we found that the overexpression of rap_pLS20 led to various phenotypic traits, including cell aggregation and death, which might have affected the correct determination of the conjugation efficiency when determined by colony formation assay. In the current study, conjugation efficiencies were determined under different conditions using a two-color fluorescence-activated flow cytometry method and measuring a single-round of pLS20-mediated transfer of a mobilizable plasmid. Under standard conditions, the conjugation efficiency obtained by fluorescence-activated flow cytometry was 23-fold higher than that obtained by colony formation. Furthermore, the efficiency difference increased to 45-fold when rap_pLS20 was overexpressed.

The Inclusion of the Food Microstructural Influence in Predictive Microbiology: State-of-the-Art

Predictive microbiology has steadily evolved into one of the most important tools to assess and control the microbiological safety of food products. Predictive models were traditionally developed based on experiments in liquid laboratory media, meaning that food microstructural effects were not represented in these models. Since food microstructure is known to exert a significant effect on microbial growth and inactivation dynamics, the applicability of predictive models is limited if food microstructure is not taken into account. Over the last 10-20 years, researchers, therefore, developed a variety of models that do include certain food microstructural influences. This review provides an overview of the most notable microstructure-including models which were developed over the years, both for microbial growth and inactivation.

Predicting plasmid persistence in microbial communities by coarse-grained modeling

Plasmids are a major type of mobile genetic elements (MGEs) that mediate horizontal gene transfer. The stable maintenance of plasmids plays a critical role in the functions and survival for microbial populations. However, predicting and controlling plasmid persistence and abundance in complex microbial communities remain challenging. Computationally, this challenge arises from the combinatorial explosion associated with the conventional modeling framework. Recently, a plasmid-centric framework (PCF) has been developed to overcome this computational bottleneck. This framework enables the derivation of a simple metric, the persistence potential, to predict plasmid persistence and abundance. Here, we discuss how PCF can be extended to account for plasmid interactions. We also discuss how such model-guided predictions of plasmid fates can benefit from the development of new experimental tools and data-driven computational methods.

Lattice Boltzmann Method in Modeling Biofilm Formation, Growth and Detachment

Biofilms are a complex and heterogeneous aggregation of multiple populations of microorganisms linked together by their excretion of extracellular polymer substances (EPS). Biofilms can cause many serious problems, such as chronic infections, food contamination and equipment corrosion, although they can be useful for constructive purposes, such as in wastewater treatment, heavy metal removal from hazardous waste sites, biofuel production, power generation through microbial fuel cells and microbially enhanced oil recovery; however, biofilm formation and growth are complex due to interactions among physicochemical and biological processes under operational and environmental conditions. Advanced numerical modeling techniques using the lattice Boltzmann method (LBM) are enabling the prediction of biofilm formation and growth and microbial community structures. This study is the first attempt to perform a general review on major contributions to LBM-based biofilm models, ranging from pioneering efforts to more recent progress. We present our understanding of the modeling of biofilm formation, growth and detachment using LBM-based models and present the fundamental aspects of various LBM-based biofilm models. We describe how the LBM couples with cellular automata (CA) and individual-based model (IbM) approaches and discuss their applications in assessing the spatiotemporal distribution of biofilms and their associated parameters and evaluating bioconversion efficiency. Finally, we discuss the main features and drawbacks of LBM-based biofilm models from ecological and biotechnological perspectives and identify current knowledge gaps and future research priorities.