Bacterial social interactions drive the emergence of differential spatial colony structures

Synthetic microbial ecology: engineering habitats for modular consortia

Microbiomes, the complex networks of micro-organisms and the molecules through which they interact, play a crucial role in health and ecology. Over at least the past two decades, engineering biology has made significant progress, impacting the bio-based industry, health, and environmental sectors; but has only recently begun to explore the engineering of microbial ecosystems. The creation of synthetic microbial communities presents opportunities to help us understand the dynamics of wild ecosystems, learn how to manipulate and interact with existing microbiomes for therapeutic and other purposes, and to create entirely new microbial communities capable of undertaking tasks for industrial biology. Here, we describe how synthetic ecosystems can be constructed and controlled, focusing on how the available methods and interaction mechanisms facilitate the regulation of community composition and output. While experimental decisions are dictated by intended applications, the vast number of tools available suggests great opportunity for researchers to develop a diverse array of novel microbial ecosystems.

Deconstructing synthetic biology across scales: a conceptual approach for training synthetic biologists

Synthetic biology allows us to reuse, repurpose, and reconfigure biological systems to address society's most pressing challenges. Developing biotechnologies in this way requires integrating concepts across disciplines, posing challenges to educating students with diverse expertise. We created a framework for synthetic biology training that deconstructs biotechnologies across scales-molecular, circuit/network, cell/cell-free systems, biological communities, and societal-giving students a holistic toolkit to integrate cross-disciplinary concepts towards responsible innovation of successful biotechnologies. We present this framework, lessons learned, and inclusive teaching materials to allow its adaption to train the next generation of synthetic biologists.

Tuning interdomain conjugation to enable in situ population modification in yeasts

The ability to modify and control natural and engineered microbiomes is essential for biotechnology and biomedicine. Fungi are critical members of most microbiomes, yet technology for modifying the fungal members of a microbiome has lagged far behind that for bacteria. Interdomain conjugation (IDC) is a promising approach, as DNA transfer from bacterial cells to yeast enables in situ modification. While such genetic transfers have been known to naturally occur in a wide range of eukaryotes and are thought to contribute to their evolution, IDC has been understudied as a technique to control fungal or fungal-bacterial consortia. One major obstacle to the widespread use of IDC is its limited efficiency. In this work, we manipulated metabolic and physical interactions between genetically tractable Escherichia coli and Saccharomyces cerevisiae to control the incidence of IDC. We test the landscape of population interactions between the bacterial donors and yeast recipients to find that bacterial commensalism leads to maximized IDC, both in culture and in mixed colonies. We demonstrate the capacity of cell-to-cell binding via mannoproteins to assist both IDC incidence and bacterial commensalism in culture and model how these tunable controls can predictably yield a range of IDC outcomes. Furthermore, we demonstrate that these controls can be utilized to irreversibly alter a recipient yeast population, by both "rescuing" a poor-growing recipient population and collapsing a stable population via a novel IDC-mediated CRISPR/Cas9 system.IMPORTANCEFungi are important but often unaddressed members of most natural and synthetic microbial communities. This work highlights opportunities for modifying yeast microbiome populations through bacterial conjugation. While conjugation has been recognized for its capacity to deliver engineerable DNA to a range of cells, its dependence on cell contact has limited its efficiency. Here, we find "knobs" to control DNA transfer, by engineering the metabolic dependence between bacterial donors and yeast recipients and by changing their ability to physically adhere to each other. Importantly, we functionally validate these "knobs" by irreversibly altering yeast populations. We use these controls to "rescue" a failing yeast population, demonstrate the capacity of conjugated CRISPR/Cas9 to depress or collapse populations, and show that conjugation can be easily interrupted by disrupting cell-to-cell binding. These results offer building blocks toward in situ mycobiome editing, with significant implications for clinical treatments of fungal pathogens and other fungal system engineering.

Resource limitation and population fluctuation drive spatiotemporal order in microbial communities

Microbial communities display complex spatiotemporal behaviors leading to spatially-structured and ordered organization driven by species interactions and environmental factors. Resource availability plays a pivotal role in shaping the dynamics of bacterial colonies. In this study, we delve into the intricate interplay between resource limitation and the emergent properties of a growing colony of two visually distinct bacterial strains having similar growth and mechanical properties. Employing an agent-based modeling and computer simulations, we analyze the resource-driven effect on segregation and sectoring, cell length regulation and nematic ordering within a growing colony. We introduce a dimensionless parameter referred to as the active layer thickness, derived from nutrient diffusion equations, indicating effective population participation due to local resource availability. Our results reveal that lower values of active layer thickness arising from decreased resource abundance lead to rougher colony fronts, fostering heightened population fluctuations within the colony and faster spatial genetic diversity loss. Our temporal analyses unveil the dynamics of mean cell length and fluctuations, showcasing how initial disturbances evolve as colonies are exposed to nutrients and subsequently settle. Furthermore, examining microscopic details, we find that lower resource levels yield diverse cell lengths and enhanced nematic ordering, driven by the increased prevalence of longer rod-shaped cells. Our investigation sheds light on the multifaceted relationship between resource constraints and bacterial colony dynamics, revealing insights into their spatiotemporal organization.

Active pH regulation facilitates Bacillus subtilis biofilm development in a minimally buffered environment

Biofilms provide individual bacteria with many advantages, yet dense cellular proliferation can also create intrinsic metabolic challenges including excessive acidification. Because such pH stress can be masked in buffered laboratory media-such as MSgg commonly used to study Bacillus subtilis biofilms-it is not always clear how such biofilms cope with minimally buffered natural environments. Here, we report how B. subtilis biofilms overcome this intrinsic metabolic challenge through an active pH regulation mechanism. Specifically, we find that these biofilms can modulate their extracellular pH to the preferred neutrophile range, even when starting from acidic and alkaline initial conditions, while planktonic cells cannot. We associate this behavior with dynamic interplay between acetate and acetoin biosynthesis and show that this mechanism is required to buffer against biofilm acidification. Furthermore, we find that buffering-deficient biofilms exhibit dysregulated biofilm development when grown in minimally buffered conditions. Our findings reveal an active pH regulation mechanism in B. subtilis biofilms that could lead to new targets to control unwanted biofilm growth.IMPORTANCEpH is known to influence microbial growth and community dynamics in multiple bacterial species and environmental contexts. Furthermore, in many bacterial species, rapid cellular proliferation demands the use of overflow metabolism, which can often result in excessive acidification. However, in the case of bacterial communities known as biofilms, these acidification challenges can be masked when buffered laboratory media are employed to stabilize the pH environment for optimal growth. Our study reveals that B. subtilis biofilms use an active pH regulation mechanism to mitigate both growth-associated acidification and external pH challenges. This discovery provides new opportunities for understanding microbial communities and could lead to new methods for controlling biofilm growth outside of buffered laboratory conditions.

Flux sampling in genome-scale metabolic modeling of microbial communities

BACKGROUND

Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates.

RESULTS

In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux.

CONCLUSIONS

Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

Spatial constraints and stochastic seeding subvert microbial arms race

Surface attached communities of microbes grow in a wide variety of environments. Often, the size of these microbial community is constrained by their physical surroundings. However, little is known about how size constraints of a colony impact the outcome of microbial competitions. Here, we use individual-based models to simulate contact killing between two bacterial strains with different killing rates in a wide range of community sizes. We found that community size has a substantial impact on outcomes; in fact, in some competitions the identity of the most fit strain differs in large and small environments. Specifically, when at a numerical disadvantage, the strain with the slow killing rate is more successful in smaller environments than in large environments. The improved performance in small spaces comes from finite size effects; stochastic fluctuations in the initial relative abundance of each strain in small environments lead to dramatically different outcomes. However, when the slow killing strain has a numerical advantage, it performs better in large spaces than in small spaces, where stochastic fluctuations now aid the fast killing strain in small communities. Finally, we experimentally validate these results by confining contact killing strains of Vibrio cholerae in transmission electron microscopy grids. The outcomes of these experiments are consistent with our simulations. When rare, the slow killing strain does better in small environments; when common, the slow killing strain does better in large environments. Together, this work demonstrates that finite size effects can substantially modify antagonistic competitions, suggesting that colony size may, at least in part, subvert the microbial arms race.

Genetic mixing and demixing on expanding spherical frontiers

Genetic fluctuation during range expansion is a key process driving evolution. When a bacterial population is expanding on a 2D surface, random fluctuations in the growth of the pioneers at the front line cause a strong demixing of genotypes. Even when there is no selective advantage, sectors of low genetic diversity are formed. Experimental studies of range expansions in surface-attached colonies of fluorescently labelled micro-organisms have contributed significantly to our understanding of fundamental evolutionary dynamics. However, experimental studies on genetic fluctuations in 3D range expansions have been sparse, despite their importance for tumour or biofilm development. We encapsulated populations of two fluorescent Escherichia coli strains in inoculation droplets (volumes [Formula: see text] nl). The confined ensemble of cells grew when embedded in a hydrogel-with nutrients-and developed 3D colonies with well-defined, sector-like regions. Using confocal laser scanning microscopy, we imaged the development of 3D colonies and the emergence of sectors. We characterized how cell concentration in the inoculation droplet controls sectors, growth rate, and the transition from branched colonies to quasi-spherical colonies. We further analysed how sectors on the surface change over time. We complement these experimental results with a modified 3D Eden growth model. The model in 3D spherical growth predicts a phase, where sectors are merging, followed by a steady increase (constant rate), and the experimentally analysed sectors were consistent with this prediction. Therefore, our results demonstrate qualitative differences between radial (2D) and spherical (3D) range expansions and their importance in gene fixation processes.

Tuning Interdomain Conjugation Toward in situ Population Modification in Yeast

The ability to modify and control natural and engineered microbiomes is essential for biotechnology and biomedicine. Fungi are critical members of most microbiomes, yet technology for modifying the fungal members of a microbiome has lagged far behind that for bacteria. Interdomain conjugation (IDC) is a promising approach, as DNA transfer from bacterial cells to yeast enables in situ modification. While such genetic transfers have been known to naturally occur in a wide range of eukaryotes, and are thought to contribute to their evolution, IDC has been understudied as a technique to control fungal or fungal-bacterial consortia. One major obstacle to widespread use of IDC is its limited efficiency. In this work, we utilize interactions between genetically tractable Escherichia coli and Saccharomyces cerevisiae to control the incidence of IDC. We test the landscape of population interactions between the bacterial donors and yeast recipients to find that bacterial commensalism leads to maximized IDC, both in culture and in mixed colonies. We demonstrate the capacity of cell-to-cell binding via mannoproteins to assist both IDC incidence and bacterial commensalism in culture, and model how these tunable controls can predictably yield a range of IDC outcomes. Further, we demonstrate that these lessons can be utilized to lastingly alter a recipient yeast population, by both "rescuing" a poor-growing recipient population and collapsing a stable population via a novel IDC-mediated CRISPR/Cas9 system.

Spatial biology of Ising-like synthetic genetic networks

BACKGROUND

Understanding how spatial patterns of gene expression emerge from the interaction of individual gene networks is a fundamental challenge in biology. Developing a synthetic experimental system with a common theoretical framework that captures the emergence of short- and long-range spatial correlations (and anti-correlations) from interacting gene networks could serve to uncover generic scaling properties of these ubiquitous phenomena.

RESULTS

Here, we combine synthetic biology, statistical mechanics models, and computational simulations to study the spatial behavior of synthetic gene networks (SGNs) in Escherichia coli quasi-2D colonies growing on hard agar surfaces. Guided by the combined mechanisms of the contact process lattice simulation and two-dimensional Ising model (CPIM), we describe the spatial behavior of bi-stable and chemically coupled SGNs that self-organize into patterns of long-range correlations with power-law scaling or short-range anti-correlations. These patterns, resembling ferromagnetic and anti-ferromagnetic configurations of the Ising model near critical points, maintain their scaling properties upon changes in growth rate and cell shape.

CONCLUSIONS

Our findings shed light on the spatial biology of coupled and bistable gene networks in growing cell populations. This emergent spatial behavior could provide insights into the study and engineering of self-organizing gene patterns in eukaryotic tissues and bacterial consortia.

Flux Sampling in Genome-scale Metabolic Modeling of Microbial Communities

Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model. However, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling may capture additional heterogeneity across cells, especially when cells exhibit sub-maximal growth rates. In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. We find significant differences in the predicted metabolism with sampling, including increased cooperative interactions and pathway-specific changes in predicted flux. Our results suggest the importance of sampling-based and objective function-independent approaches to evaluate metabolic interactions and emphasize their utility in quantitatively studying interactions between cells and organisms.

Impact of direct physical association and motility on fitness of a synthetic interkingdom microbial community

Mutualistic exchange of metabolites can play an important role in microbial communities. Under natural environmental conditions, such exchange may be compromised by the dispersal of metabolites and by the presence of non-cooperating microorganisms. Spatial proximity between members during sessile growth on solid surfaces has been shown to promote stabilization of cross-feeding communities against these challenges. Nonetheless, many natural cross-feeding communities are not sessile but rather pelagic and exist in turbulent aquatic environments, where partner proximity is often achieved via direct cell-cell adhesion, and cooperation occurs between physically associated cells. Partner association in aquatic environments could be further enhanced by motility of individual planktonic microorganisms. In this work, we establish a model bipartite cross-feeding community between bacteria and yeast auxotrophs to investigate the impact of direct adhesion between prokaryotic and eukaryotic partners and of bacterial motility in a stirred mutualistic co-culture. We demonstrate that adhesion can provide fitness benefit to the bacterial partner, likely by enabling local metabolite exchange within co-aggregates, and that it counteracts invasion of the community by a non-cooperating cheater strain. In a turbulent environment and at low cell densities, fitness of the bacterial partner and its competitiveness against a non-cooperating strain are further increased by motility that likely facilitates partner encounters and adhesion. These results suggest that, despite their potential fitness costs, direct adhesion between partners and its enhancement by motility may play key roles as stabilization factors for metabolic communities in turbulent aquatic environments.

Assembly strategies for polyethylene-degrading microbial consortia based on the combination of omics tools and the "Plastisphere"

Numerous microorganisms and other invertebrates that are able to degrade polyethylene (PE) have been reported. However, studies on PE biodegradation are still limited due to its extreme stability and the lack of explicit insights into the mechanisms and efficient enzymes involved in its metabolism by microorganisms. In this review, current studies of PE biodegradation, including the fundamental stages, important microorganisms and enzymes, and functional microbial consortia, were examined. Considering the bottlenecks in the construction of PE-degrading consortia, a combination of top-down and bottom-up approaches is proposed to identify the mechanisms and metabolites of PE degradation, related enzymes, and efficient synthetic microbial consortia. In addition, the exploration of the plastisphere based on omics tools is proposed as a future principal research direction for the construction of synthetic microbial consortia for PE degradation. Combining chemical and biological upcycling processes for PE waste could be widely applied in various fields to promote a sustainable environment.

Mechanisms driving spatial distribution of residents in colony biofilms: an interdisciplinary perspective

Biofilms are consortia of microorganisms that form collectives through the excretion of extracellular matrix compounds. The importance of biofilms in biological, industrial and medical settings has long been recognized due to their emergent properties and impact on surrounding environments. In laboratory situations, one commonly used approach to study biofilm formation mechanisms is the colony biofilm assay, in which cell communities grow on solid-gas interfaces on agar plates after the deposition of a population of founder cells. The residents of a colony biofilm can self-organize to form intricate spatial distributions. The assay is ideally suited to coupling with mathematical modelling due to the ability to extract a wide range of metrics. In this review, we highlight how interdisciplinary approaches have provided deep insights into mechanisms causing the emergence of these spatial distributions from well-mixed inocula.

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.

Individual-Based Modeling of Spatial Dynamics of Chemotactic Microbial Populations

One important direction of synthetic biology is to establish desired spatial structures from microbial populations. Underlying this structural development process are different driving factors, among which bacterial motility and chemotaxis serve as a major force. Here, we present an individual-based, biophysical computational framework for mechanistic and multiscale simulation of the spatiotemporal dynamics of motile and chemotactic microbial populations. The framework integrates cellular movement with spatial population growth, mechanical and chemical cellular interactions, and intracellular molecular kinetics. It is validated by a statistical comparison of single-cell chemotaxis simulations with reported experiments. The framework successfully captures colony range expansion of growing isogenic populations and also reveals chemotaxis-modulated, spatial patterns of a two-species amensal community. Partial differential equation-based models subsequently validate these simulation findings. This study provides a versatile computational tool to uncover the fundamentals of microbial spatial ecology as well as to facilitate the design of synthetic consortia for desired spatial patterns.

Emergence of sector and spiral patterns from a two-species mutualistic cross-feeding model

The ubiquitous existence of microbial communities marks the importance of understanding how species interact within the community to coexist and their spatial organization. We study a two-species mutualistic cross-feeding model through a stochastic cellular automaton on a square lattice using kinetic Monte Carlo simulation. Our model encapsulates the essential dynamic processes such as cell growth, and nutrient excretion, diffusion and uptake. Focusing on the interplay among nutrient diffusion and individual cell division, we discover three general classes of colony morphology: co-existing sectors, co-existing spirals, and engulfment. When the cross-feeding nutrient is widely available, either through high excretion or fast diffusion, a stable circular colony with alternating species sector emerges. When the consumer cells rely on being spatially close to the producers, we observe a stable spiral. We also see one species being engulfed by the other when species interfaces merge due to stochastic fluctuation. By tuning the diffusion rate and the growth rate, we are able to gain quantitative insights into the structures of the sectors and the spirals.

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.

Effect of Bacillus subtilis fortified inoculation on the microbial communities in different niches of Daqu

Daqu is the fermentation starter of Baijiu, which is one of the six most renowned distilled spirits. Studies have found that Bacillus is one of the dominant microbial genera in Daqu, and Bacillus subtilis is known to produce amylase, an important enzyme that influences the quality of Daqu. This study aims to explore the influence of B. subtilis inoculation on the microbial community structure in different niches of Daqu. We studied the microbial community structure of the natural inoculated Daqu (i.e., control Daqu) and the fortified Daqu (i.e., B. subtilis-inoculated Daqu) by amplicon sequencing. Our results showed that compared with the control Daqu rind microbial community, the relative abundance of Bacillus, Aspergillus, Thermomyces, and Rasamsonia in the fortified Daqu rind microbial community increased, and the relative abundance of Weissella, Lactobacillus, Leuconostoc, and Pichia decreased. Compared with the control Daqu core microbial community, the relative abundance of Bacillus in the fortified Daqu core microbial community also increased, but the relative abundance of Pseudomonas and Paecilomyces decreased. The effect of B. subtilis inoculation on the rind microbial community of Daqu was more significant. In addition, the bacterial community of Daqu was more susceptible to the effect of the B. subtilis inoculation than was the fungal community of Daqu. The correlation between the bacterial community of Daqu and the fungal community of Daqu increased significantly after the B. subtilis inoculation. These results provide an important theoretical basis for the production of fortified Daqu.

Engineering microbial consortia with rationally designed cellular interactions

Synthetic microbial consortia represent a frontier of synthetic biology that promises versatile engineering of cellular functions. They are primarily developed through the design and construction of cellular interactions that coordinate individual dynamics and generate collective behaviors. Here we review recent advances in the engineering of synthetic communities through cellular-interaction programming. We first examine fundamental building blocks for intercellular communication and unidirectional positive and negative interactions. We then recap the assembly of the building blocks for creating bidirectional interactions in two-species ecosystems, which is followed by the discussion of engineering toward complex communities with increasing species numbers, under spatial contexts, and via model-guided design. We conclude by summarizing major challenges and future opportunities of engineered microbial ecosystems.

Research Progress and Hopeful Strategies of Application of Quorum Sensing in Food, Agriculture and Nanomedicine

Quorum sensing (QS) regulates the expression of several genes including motility, biofilm development, virulence expression, population density detection and plasmid conjugation. It is based on “autoinducers”, small molecules that microorganisms produce and release in the extracellular milieu. The biochemistry of quorum sensing is widely discussed and numerous papers are available to scientists. The main purpose of this research is to understand how knowledge about this mechanism can be exploited for the benefit of humans and the environment. Here, we report the most promising studies on QS and their resulting applications in different fields of global interest: food, agriculture and nanomedicine.

Spatial patterns in ecological systems: from microbial colonies to landscapes

Self-organized spatial patterns are ubiquitous in ecological systems and allow populations to adopt non-trivial spatial distributions starting from disordered configurations. These patterns form due to diverse nonlinear interactions among organisms and between organisms and their environment, and lead to the emergence of new (eco)system-level properties unique to self-organized systems. Such pattern consequences include higher resilience and resistance to environmental changes, abrupt ecosystem collapse, hysteresis loops, and reversal of competitive exclusion. Here, we review ecological systems exhibiting self-organized patterns. We establish two broad pattern categories depending on whether the self-organizing process is primarily driven by nonlinear density-dependent demographic rates or by nonlinear density-dependent movement. Using this organization, we examine a wide range of observational scales, from microbial colonies to whole ecosystems, and discuss the mechanisms hypothesized to underlie observed patterns and their system-level consequences. For each example, we review both the empirical evidence and the existing theoretical frameworks developed to identify the causes and consequences of patterning. Finally, we trace qualitative similarities across systems and propose possible ways of developing a more quantitative understanding of how self-organization operates across systems and observational scales in ecology.

Well-mixed Lotka-Volterra model with random strongly competitive interactions

The random Lotka-Volterra model is widely used to describe the dynamical and thermodynamic features of ecological communities. In this work, we consider random symmetric interactions between species and analyze the strongly competitive interaction case. We investigate different scalings for the distribution of the interactions with the number of species and try to bridge the gap with previous works. Our results show two different behaviors for the mean abundance at zero and finite temperature, respectively, with a continuous crossover between the two. We confirm and extend previous results obtained for weak interactions: at zero temperature, even in the strong competitive interaction limit, the system is in a multiple-equilibria phase, whereas at finite temperature only a unique stable equilibrium can exist. Finally, we establish the qualitative phase diagrams and compare the species abundance distributions in the two cases.

Engineered microbial consortia: strategies and applications

Many applications of microbial synthetic biology, such as metabolic engineering and biocomputing, are increasing in design complexity. Implementing complex tasks in single populations can be a challenge because large genetic circuits can be burdensome and difficult to optimize. To overcome these limitations, microbial consortia can be engineered to distribute complex tasks among multiple populations. Recent studies have made substantial progress in programming microbial consortia for both basic understanding and potential applications. Microbial consortia have been designed through diverse strategies, including programming mutualistic interactions, using programmed population control to prevent overgrowth of individual populations, and spatial segregation to reduce competition. Here, we highlight the role of microbial consortia in the advances of metabolic engineering, biofilm production for engineered living materials, biocomputing, and biosensing. Additionally, we discuss the challenges for future research in microbial consortia.

The spatial organization of microbial communities during range expansion

Microbes in nature often live in dense and diverse communities exhibiting a variety of spatial structures. Microbial range expansion is a universal ecological process that enables populations to form spatial patterns. It can be driven by both passive and active processes, for example, mechanical forces from cell growth and bacterial motility. In this review, we provide a taste of recent creative and sophisticated efforts being made to address basic questions in spatial ecology and pattern formation during range expansion. We especially highlight the role of motility to shape community structures, and discuss the research challenges and future directions.

Repulsive expansion dynamics in colony growth and gene expression

Spatial expansion of a population of cells can arise from growth of microorganisms, plant cells, and mammalian cells. It underlies normal or dysfunctional tissue development, and it can be exploited as the foundation for programming spatial patterns. This expansion is often driven by continuous growth and division of cells within a colony, which in turn pushes the peripheral cells outward. This process generates a repulsion velocity field at each location within the colony. Here we show that this process can be approximated as coarse-grained repulsive-expansion kinetics. This framework enables accurate and efficient simulation of growth and gene expression dynamics in radially symmetric colonies with homogenous z-directional distribution. It is robust even if cells are not spherical and vary in size. The simplicity of the resulting mathematical framework also greatly facilitates generation of mechanistic insights.

Spatiotemporal dynamics are ubiquitous in biology. To understand these phenomena in nature or to program them using synthetic gene circuits, it is critical to resort to mathematical modeling to deduce mechanistic insights or to explore plausible outcomes. Historically, modeling of spatiotemporal dynamics depends on the use of agent-based models or their continuum counterparts consisting of partial differential equations. Here, we show that a class of colony expansion can be treated as being driven by the steric force generated by growing and diving cells. This approximation leads to a drastically simplified framework consisting of only ordinary differential equations. This framework greatly improves the computational efficiency and facilitates development of mechanistic insights into the dynamics of colony growth and pattern formation.

Heterogeneous bacterial swarms with mixed lengths

Heterogeneous systems of active matter exhibit a range of complex emergent dynamical patterns. In particular, it is difficult to predict the properties of the mixed system based on its constituents. These considerations are particularly significant for understanding realistic bacterial swarms, which typically develop heterogeneities even when grown from a single cell. Here, mixed swarms of cells with different aspect ratios are studied both experimentally and in simulations. In contrast with previous theory, there is no macroscopic phase segregation. However, locally, long cells act as nucleation cites, around which aggregates of short, rapidly moving cells can form, resulting in enhanced swarming speeds. On the other hand, high fractions of long cells form a bottleneck for efficient swarming. Our results suggest a physical advantage for the spontaneous heterogeneity of bacterial swarm populations.

Influence of interspecies interactions on the spatial organization of dual species bacterial communities

Interspecies interactions in bacterial biofilms have important impacts on the composition and function of communities in natural and engineered systems. To investigate these interactions, synthetic communities provide experimentally tractable systems. Biofilms grown on agar-surfaces have been used for investigating the eco-evolutionary and biophysical forces that determine community composition and spatial distribution of bacteria. Prior studies have used genetically identical bacterial strains and strains with specific mutations, that express different fluorescent proteins, to investigate intraspecies interactions. Here, we investigated interspecies interactions and, specifically, determined the community composition and spatial distribution in synthetic communities of Pseudomonas aeruginosa, Pseudomonas protegens and Klebsiella pneumoniae. Using quantitative microscopic imaging, we found that interspecies interactions in multispecies colonies were influenced by type IV pilus mediated motility, extracellular matrix secretion, environmental parameters, and these effects were also influenced by the specific partner in the dual species combinations. These results indicate that the patterns observable in mixed species colonies can be used to understand the mechanisms that drive interspecies interactions, which are dependent on the interplay between specific species’ physiology and environmental conditions.

•

Spatial patterns in bacterial colonies are species and interaction dependent.

•

Surface motility and extracellar matrix production affect interspecies interactions.

•

Agar surface colonies show how bacteria interact in biofilms.

Theoretical modelling of competitive microbial range expansion with heterogeneous mechanical interactions

Microbial range expansion experiments provide insight into the complex link between dynamic structure, pattern formation and evolutionary dynamics of growing populations. In this work, we develop a theoretical model in order to investigate the interplay of growth statistics and mechanical interactions which are implemented as division driven pushing and swapping of cells. For the case of the competitive growth of a strongly and a weakly interacting strain we investigate the influence of different mean division times, as well as different mechanical interactions on the development of the colony. Our results show that the susceptibility to cell division induced pushing has a much stronger influence on the structure of the colony than cell sorting towards the colony's perimeter. Motivated by microbial range expansion experiments of Neisseria gonorrhoeae bacteria, we also consider the influence of mutating cells on the structure of the colony. We show that the outgrowth of the three different strains is strongly influenced by the relative strengths of their mechanical interaction. The experimentally observed patterns are reproduced for mechanical interactions of the mutants, which range between those of the strongly and weakly interacting strain.

Intra-colony channels in E. coli function as a nutrient uptake system

The ability of microorganisms to grow as aggregated assemblages has been known for many years, however their structure has remained largely unexplored across multiple spatial scales. The development of the Mesolens, an optical system which uniquely allows simultaneous imaging of individual bacteria over a 36 mm2 field of view, has enabled the study of mature Escherichia coli macro-colony biofilm architecture like never before. The Mesolens enabled the discovery of intra-colony channels on the order of 10 μm in diameter, that are integral to E. coli macro-colony biofilms and form as an emergent property of biofilm growth. These channels have a characteristic structure and re-form after total mechanical disaggregation of the colony. We demonstrate that the channels are able to transport particles and play a role in the acquisition of and distribution of nutrients through the biofilm. These channels potentially offer a new route for the delivery of dispersal agents for antimicrobial drugs to biofilms, ultimately lowering their impact on public health and industry.

The unexplored bacterial lifestyle on leaf surface

Social interactions impact microbial communities and these relationships are mediated by small molecules. The chemical ecology of bacteria on the phylloplane environment is still little explored. The harsh environmental conditions found on leaf surface require high metabolic performances of the bacteria in order to survive. That is interesting both for scientific fields of prospecting natural molecules and for the ecological studies. Important queries about the bacterial lifestyle on leaf surface remain not fully comprehended. Does the hostility of the environment increase the populations' cellular altruism by the production of molecules, which can benefit the whole community? Or does the reverse occur and the production of molecules related to competition between species is increased? Does the phylogenetic distance between the bacterial populations influence the chemical profile during social interactions? Do phylogenetically related bacteria tend to cooperate more than the distant ones? The phylloplane contains high levels of yet uncultivated microorganisms, and understanding the molecular basis of the social networks on this habitat is crucial to gain new insights on the ecology of the mysterious community members due to interspecies molecular dependence. Here, we review and discuss what is known about bacterial social interactions and their chemical lifestyle on leaf surface.

Integrative Circuit-Host Modeling of a Genetic Switch in Varying Environments

Synthetic biology is advancing into a new phase where real-world applications are emphasized. There is hence an urgent need for mathematical modeling that can quantitatively describe the behaviors of genetic devices in natural, fluctuating environments. We utilize an integrative circuit-host modeling framework to examine the dynamics of a genetic switch and its host cell in varying environments. For both steady-state and transient cases, we find increasing nutrient reduces the bistability region of the phase space and eventually drives the switch from bistability to monostability. In response, cellular growth and proteome partitioning experience the same transition. Antibiotic perturbations cause the similar circuit and host responses as nutrient variations. However, one difference is the trend of growth rate, which augments with nutrient but declines with antibiotic levels. The framework provides a mechanistic scheme to account for both the dynamic and static characteristics of the circuit-host system upon environmental perturbations, underscoring the intimacy of gene circuits and their hosts and elucidating the complexity of circuit behaviors arising from environmental variations.

Blue-Light-Switchable Bacterial Cell-Cell Adhesions Enable the Control of Multicellular Bacterial Communities

Although the fundamental importance and biotechnological potential of multibacterial communities, also called biofilms, are well-known, our ability to control them is limited. We present a new way of dynamically controlling bacteria-bacteria adhesions by using blue light and how these photoswitchable adhesions can be used to regulate multicellularity and associated bacterial behavior. To achieve this, the photoswitchable proteins nMagHigh and pMagHigh were expressed on bacterial surfaces as adhesins to allow multicellular clusters to assemble under blue light and reversibly disassemble in the dark. Regulation of the bacterial cell-cell adhesions with visible light provides unique advantages including high spatiotemporal control, tunability, and noninvasive remote regulation. Moreover, these photoswitchable adhesions make it possible to regulate collective bacterial functions including aggregation, quorum sensing, biofilm formation, and metabolic cross-feeding between auxotrophic bacteria with light. Overall, the photoregulation of bacteria-bacteria adhesions provides a new way of studying bacterial cell biology and will enable the design of biofilms for biotechnological applications.

Three-Way Interactions in an Artificial Community of Bacterial Strains Directly Isolated From the Environment and Their Effect on the System Population Dynamics

This work is motivated by previous studies that have analyzed the population ecology of a collection of culturable thermoresistant bacteria, isolated from the Churince lagoon in Cuatro Cienegas, Mexico. In particular, it is aimed at testing a hypothesis from a modeling study, which states that antagonistic and sensitive bacteria co-exist thanks to resistant bacteria that protect sensitive ones by forming physical barriers. We selected three different bacterial strains from the referred collection: one antagonistic, one sensitive, and one resistant, and studied the population dynamics of mixed colonies. Our results show that, although the proposed protective mechanism does not work in this case, the resistant strain confers some kind of protection to sensitive bacteria. Further modeling and experimental results suggest that the presence of resistant bacteria indirectly improves the probability that patches of sensitive bacteria grow in a mixed colony. More precisely, our results suggest that by making antagonistic bacteria produce and secrete an antagonistic substance (with the concomitant metabolic cost and growth rate reduction), resistant bacteria increase the likelihood that sensitive bacteria locally outcompete antagonistic ones.

Bioturbation effect of fortified Daqu on microbial community and flavor metabolite in Chinese strong-flavor liquor brewing microecosystem

Traditional spontaneous fermentation in microecosystem with microbial successions and environmental variables results in inconsistent quality of fermented foods. We therefore propose the directional bioturbation for microbiota regulation and metabolites production in food fermentation. Here, we revealed the bioturbation effect of fortified Daqu on microbial community based on taxonomic composition, co-occurrence network, and metabolic potential, using Chinese strong-flavor liquor fermentation as a microecosystem. According to principal coordinate analysis, microbial communities were obviously influenced by the bioturbation of fortified Daqu. More specifically, bioturbation increased the abundances of Caproiciproducens, Clostridium, Aspergillus, Candida, Methanobacterium, and Methanosarcina, while decreased that of Lactobacillus. Meanwhile, higher abundances of most genes that encoding enzymes involved in interspecies hydrogen transfer between hexanoic acid bacteria and methanogens were observed in the bioturbated ecosystem by PICRUSt approach. Additionally, co-occurrence analysis showed that bioturbation increased the diversity and complexity of interspecies interactions in microecosystem, which contributed to higher production of flavor metabolites such as hexanoic acid, ethyl hexanoate, and hexyl hexanoate. These results indicated that the bioturbation of fortified Daqu is feasible for flavor metabolism by interspecies interactions of functional microbiota in liquor fermentation. Taken together, it is of great importance for regulating Chinese liquor and even other foods fermentation by bioturbation.

Spatial coordination in a mutually beneficial bacterial community enhances its antibiotic resistance

Microbial communities can survive in complex and variable environments by using different cooperative strategies. However, the behaviors of these mutuality formed communities remain poorly understood, particularly with regard to the characteristics of spatial cooperation. Here, we selected two Escherichia coli strains, designated as the nutrition provider and the antibiotic protector, respectively, for construction of a mutually beneficial bacterial community that could be used to study these behaviors. We found that in addition to the functional mutualism, the two strains also cooperated through their spatial distribution. Under antibiotic pressure, the bacterial distribution changed to yield different spatial distributions, which resulted in community growth advantages beyond functional cooperation. The mutualistic behavior of these two strains suggested that similar communities could also use variations in spatial distribution to improve their survival rates in a natural environment or under the action of antibiotics.

Single-species fragmentation: The role of density-dependent feedback

Internal feedback is commonly present in biological populations and can play a crucial role in the emergence of collective behavior. To describe the temporal evolution of the distribution of a single-species population, we consider a generalization of the Fisher-KPP equation. This equation includes the elementary processes of random motion, reproduction, and, importantly, nonlocal interspecific competition, which introduces a spatial scale of interaction. In addition, we take into account feedback mechanisms in diffusion and growth processes, mimicked by power-law density dependencies. This feedback includes, for instance, anomalous diffusion, reaction to overcrowding or to the rarefaction of the population, as well as Allee-like effects. We show that, depending on the kind of feedback that takes place, the population can self-organize splitting into disconnected subpopulations, in the absence of external constraints. Through extensive numerical simulations, we investigate the temporal evolution and the characteristics of the stationary population distribution in the one-dimensional case. We discuss the crucial role that density-dependence has on pattern formation, particularly on fragmentation, which can bring important consequences to processes such as epidemic spread and speciation.

A unifying framework for interpreting and predicting mutualistic systems

Coarse-grained rules are widely used in chemistry, physics and engineering. In biology, however, such rules are less common and under-appreciated. This gap can be attributed to the difficulty in establishing general rules to encompass the immense diversity and complexity of biological systems. Furthermore, even when a rule is established, it is often challenging to map it to mechanistic details and to quantify these details. Here we report a framework that addresses these challenges for mutualistic systems. We first deduce a general rule that predicts the various outcomes of mutualistic systems, including coexistence and productivity. We further develop a standardized machine-learning-based calibration procedure to use the rule without the need to fully elucidate or characterize their mechanistic underpinnings. Our approach consistently provides explanatory and predictive power with various simulated and experimental mutualistic systems. Our strategy can pave the way for establishing and implementing other simple rules for biological systems.

Understanding and Engineering Distributed Biochemical Pathways in Microbial Communities

Microbiomes impact nearly every environment on Earth by modulating the molecular composition of the environment. Temporally changing environmental stimuli and spatial organization are major variables shaping the structure and function of microbiomes. The web of interactions among members of these communities and between the organisms and the environment dictates microbiome functions. Microbial interactions are major drivers of microbiomes and are modulated by spatiotemporal parameters. A mechanistic and quantitative understanding of ecological, molecular, and environmental forces shaping microbiomes could inform strategies to control microbiome dynamics and functions. Major challenges for harnessing the potential of microbiomes for diverse applications include the development of predictive modeling frameworks and tools for precise manipulation of microbiome behaviors.

Spatial interference scale as a determinant of microbial range expansion

In microbial communities, social interactions such as competition occur ubiquitously across multiple spatial scales from local proximity to remote distance. However, it remains unclear how such a spatial variation of interaction contributes to the structural development of microbial populations. Here, we developed synthetic consortia, biophysical theory, and simulations to elucidate the role of spatial interference scale in governing ecosystem organization during range expansion. For consortia with unidirectional interference, we discovered that, at growing fronts, the extinction time of toxin-sensitive species is reciprocal to the spatial interference scale. In contrast, for communities with bidirectional interference, their structures diverge into distinct monoculture colonies under different initial conditions, with the corresponding separatrix set by the spatial scale of interference. Near the separatrix, ecosystem development becomes noise-driven and yields opposite structures. Our results establish spatial interaction scale as a key determinant for microbial range expansion, providing insights into microbial spatial organization and synthetic ecosystem engineering.

Intercellular adhesion promotes clonal mixing in growing bacterial populations

Dense bacterial communities, known as biofilms, can have functional spatial organization driven by self-organizing chemical and physical interactions between cells, and their environment. In this work, we investigated intercellular adhesion, a pervasive property of bacteria in biofilms, to identify effects on the internal structure of bacterial colonies. We expressed the self-recognizing ag43 adhesin protein in Escherichia coli to generate adhesion between cells, which caused aggregation in liquid culture and altered microcolony morphology on solid media. We combined the adhesive phenotype with an artificial colony patterning system based on plasmid segregation, which marked clonal lineage domains in colonies grown from single cells. Engineered E. coli were grown to colonies containing domains with varying adhesive properties, and investigated with microscopy, image processing and computational modelling techniques. We found that intercellular adhesion elongated the fractal-like boundary between cell lineages only when both domains within the colony were adhesive, by increasing the rotational motion during colony growth. Our work demonstrates that adhesive intercellular interactions can have significant effects on the spatial organization of bacterial populations, which can be exploited for biofilm engineering. Furthermore, our approach provides a robust platform to study the influence of intercellular interactions on spatial structure in bacterial populations.

Visualizing Bacterial Colony Morphologies Using Time-Lapse Imaging Chamber MOCHA

Capturing microbial growth on a macroscopic scale is of great importance to further our understanding of microbial life. However, methods for imaging microbial life on a scale of millimeters to centimeters are often limited by designs that have poor environmental control, resulting in dehydration of the agar plate within just a few days. Here, we created MOCHA (microbial chamber), a simple but effective chamber that allows users to study microbial growth for extended periods (weeks) in a stable environment. Agar hydration is maintained with a double-decker design, in which two glass petri dishes are connected by a wick, allowing the lower plate to keep the upper plate hydrated. This flexible chamber allows the observation of a variety of microbiological phenomena, such as the growth and development of single bacterial and fungal colonies, interspecies interactions, swarming motility, and pellicle formation.IMPORTANCE Detailed study of microbial life on the colony scale of millimeters to centimeters has been lagging considerably behind microscopic inspection of microbes. One major reason for this is the lack of inexpensive instrumentation that can reproducibly capture images in a controlled environment. In this study, we present the design and use of a unique chamber that was used to produce several time-lapse movies that aimed to capture the diversity of microbial colony phenotypes over long periods.

The spatiotemporal system dynamics of acquired resistance in an engineered microecology

Great strides have been made in the understanding of complex networks; however, our understanding of natural microecologies is limited. Modelling of complex natural ecological systems has allowed for new findings, but these models typically ignore the constant evolution of species. Due to the complexity of natural systems, unanticipated interactions may lead to erroneous conclusions concerning the role of specific molecular components. To address this, we use a synthetic system to understand the spatiotemporal dynamics of growth and to study acquired resistance in vivo. Our system differs from earlier synthetic systems in that it focuses on the evolution of a microecology from a killer-prey relationship to coexistence using two different non-motile Escherichia coli strains. Using empirical data, we developed the first ecological model emphasising the concept of the constant evolution of species, where the survival of the prey species is dependent on location (distance from the killer) or the evolution of resistance. Our simple model, when expanded to complex microecological association studies under varied spatial and nutrient backgrounds may help to understand the complex relationships between multiple species in intricate natural ecological networks. This type of microecological study has become increasingly important, especially with the emergence of antibiotic-resistant pathogens.

A New Improved and Extended Version of the Multicell Bacterial Simulator gro

gro is a cell programming language developed in Klavins Lab for simulating colony growth and cell-cell communication. It is used as a synthetic biology prototyping tool for simulating multicellular biocircuits and microbial consortia. In this work, we present several extensions made to gro that improve the performance of the simulator, make it easier to use, and provide new functionalities. The new version of gro is between 1 and 2 orders of magnitude faster than the original version. It is able to grow microbial colonies with up to 10⁵ cells in less than 10 min. A new library, CellEngine, accelerates the resolution of spatial physical interactions between growing and dividing cells by implementing a new shoving algorithm. A genetic library, CellPro, based on Probabilistic Timed Automata, simulates gene expression dynamics using simplified and easy to compute digital proteins. We also propose a more convenient language specification layer, ProSpec, based on the idea that proteins drive cell behavior. CellNutrient, another library, implements Monod-based growth and nutrient uptake functionalities. The intercellular signaling management was improved and extended in a library called CellSignals. Finally, bacterial conjugation, another local cell-cell communication process, was added to the simulator. To show the versatility and potential outreach of this version of gro, we provide studies and novel examples ranging from synthetic biology to evolutionary microbiology. We believe that the upgrades implemented for gro have made it into a powerful and fast prototyping tool capable of simulating a large variety of systems and synthetic biology designs.

A mechanistic Individual-based Model of microbial communities

Accurate predictive modelling of the growth of microbial communities requires the credible representation of the interactions of biological, chemical and mechanical processes. However, although biological and chemical processes are represented in a number of Individual-based Models (IbMs) the interaction of growth and mechanics is limited. Conversely, there are mechanically sophisticated IbMs with only elementary biology and chemistry. This study focuses on addressing these limitations by developing a flexible IbM that can robustly combine the biological, chemical and physical processes that dictate the emergent properties of a wide range of bacterial communities. This IbM is developed by creating a microbiological adaptation of the open source Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS). This innovation should provide the basis for “bottom up” prediction of the emergent behaviour of entire microbial systems. In the model presented here, bacterial growth, division, decay, mechanical contact among bacterial cells, and adhesion between the bacteria and extracellular polymeric substances are incorporated. In addition, fluid-bacteria interaction is implemented to simulate biofilm deformation and erosion. The model predicts that the surface morphology of biofilms becomes smoother with increased nutrient concentration, which agrees well with previous literature. In addition, the results show that increased shear rate results in smoother and more compact biofilms. The model can also predict shear rate dependent biofilm deformation, erosion, streamer formation and breakup.

Self-organization of bacterial communities against environmental pH variation: Controlled chemotactic motility arranges cell population structures in biofilms

As with many living organisms, bacteria often live on the surface of solids, such as foods, organisms, buildings and soil. Compared with dispersive behavior in liquid, bacteria on surface environment exhibit significantly restricted mobility. They have access to only limited resources and cannot be liberated from the changing environment. Accordingly, appropriate collective strategies are necessarily required for long-term growth and survival. However, in spite of our deepening knowledge of the structure and characteristics of individual cells, strategic self-organizing dynamics of their community is poorly understood and therefore not yet predictable. Here, we report a morphological change in Bacillus subtilis biofilms due to environmental pH variations, and present a mathematical model for the macroscopic spatio-temporal dynamics. We show that an environmental pH shift transforms colony morphology on hard agar media from notched 'volcano-like' to round and front-elevated 'crater-like'. We discover that a pH-dependent dose-response relationship between nutritional resource level and quantitative bacterial motility at the population level plays a central role in the mechanism of the spatio-temporal cell population structure design in biofilms.

Artificial Symmetry-Breaking for Morphogenetic Engineering Bacterial Colonies

Morphogenetic engineering is an emerging field that explores the design and implementation of self-organized patterns, morphologies, and architectures in systems composed of multiple agents such as cells and swarm robots. Synthetic biology, on the other hand, aims to develop tools and formalisms that increase reproducibility, tractability, and efficiency in the engineering of biological systems. We seek to apply synthetic biology approaches to the engineering of morphologies in multicellular systems. Here, we describe the engineering of two mechanisms, symmetry-breaking and domain-specific cell regulation, as elementary functions for the prototyping of morphogenetic instructions in bacterial colonies. The former represents an artificial patterning mechanism based on plasmid segregation while the latter plays the role of artificial cell differentiation by spatial colocalization of ubiquitous and segregated components. This separation of patterning from actuation facilitates the design-build-test-improve engineering cycle. We created computational modules for CellModeller representing these basic functions and used it to guide the design process and explore the design space in silico. We applied these tools to encode spatially structured functions such as metabolic complementation, RNAPT7 gene expression, and CRISPRi/Cas9 regulation. Finally, as a proof of concept, we used CRISPRi/Cas technology to regulate cell growth by controlling methionine synthesis. These mechanisms start from single cells enabling the study of morphogenetic principles and the engineering of novel population scale structures from the bottom up.