Mechanically-driven phase separation in a growing bacterial colony

Mechanics limits ecological diversity and promotes heterogeneity in confined bacterial communities

Multispecies bacterial populations often inhabit confined and densely packed environments where spatial competition determines the ecological diversity of the community. However, the role of mechanical interactions in shaping the ecology is still poorly understood. Here, we study a model system consisting of two populations of nonmotile Escherichia coli bacteria competing within open, monolayer microchannels. The competitive dynamics is observed to be biphasic: After seeding, either one strain rapidly fixates or both strains orient into spatially stratified, stable communities. We find that mechanical interactions with other cells and local spatial constraints influence the resulting community ecology in unexpected ways, severely limiting the overall diversity of the communities while simultaneously allowing for the establishment of stable, heterogeneous populations of bacteria displaying disparate growth rates. Surprisingly, the populations have a high probability of coexisting even when one strain has a significant growth advantage. A more coccus morphology is shown to provide a selective advantage, but agent-based simulations indicate this is due to hydrodynamic and adhesion effects within the microchannel and not from breaking of the nematic ordering. Our observations are qualitatively reproduced by a simple Pólya urn model, which suggests the generality of our findings for confined population dynamics and highlights the importance of early colonization conditions on the resulting diversity and ecology of bacterial communities. These results provide fundamental insights into the determinants of community diversity in dense confined ecosystems where spatial exclusion is central to competition as in organized biofilms or intestinal crypts.

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.

Agent-based modeling of stress anisotropy driven nematic ordering in growing biofilms

Living active collectives have evolved with remarkable self-patterning capabilities to adapt to the physical and biological constraints crucial for their growth and survival. However, the intricate process by which complex multicellular patterns emerge from a single founder cell remains elusive. In this study, we utilize an agent-based model, validated through single-cell microscopy imaging, to track the three-dimensional (3D) morphodynamics of cells within growing bacterial biofilms encased by agarose gels. The confined growth conditions give rise to a spatiotemporally heterogeneous stress landscape within the biofilm. In the core of the biofilm, where high hydrostatic and low shear stresses prevail, cell packing appears disordered. In contrast, near the gel-cell interface, a state of high shear stress and low hydrostatic stress emerges, driving nematic ordering, albeit with a time delay inherent to shear stress relaxation. Strikingly, we observe a robust spatiotemporal correlation between stress anisotropy and nematic ordering within these confined biofilms. This correlation suggests a mechanism whereby stress anisotropy plays a pivotal role in governing the spatial organization of cells. The reciprocity between stress anisotropy and cell ordering in confined biofilms opens new avenues for innovative 3D mechanically guided patterning techniques for living active collectives, which hold significant promise for a wide array of environmental and biomedical applications.

Substrate geometry affects population dynamics in a bacterial biofilm

Biofilms inhabit a range of environments, such as dental plaques or soil micropores, often characterized by noneven surfaces. However, the impact of surface irregularities on the population dynamics of biofilms remains elusive, as most experiments are conducted on flat surfaces. Here, we show that the shape of the surface on which a biofilm grows influences genetic drift and selection within the biofilm. We culture Escherichia coli biofilms in microwells with a corrugated bottom surface and observe the emergence of clonal sectors whose size corresponds to that of the corrugations, despite no physical barrier separating different areas of the biofilm. The sectors are remarkably stable and do not invade each other; we attribute this stability to the characteristics of the velocity field within the biofilm, which hinders mixing and clonal expansion. A microscopically detailed computer model fully reproduces these findings and highlights the role of mechanical interactions such as adhesion and friction in microbial evolution. The model also predicts clonal expansion to be limited even for clones with a significant growth advantage-a finding which we confirm experimentally using a mixture of antibiotic-sensitive and antibiotic-resistant mutants in the presence of sublethal concentrations of the antibiotic rifampicin. The strong suppression of selection contrasts sharply with the behavior seen in range expansion experiments in bacterial colonies grown on agar. Our results show that biofilm population dynamics can be affected by patterning the surface and demonstrate how a better understanding of the physics of bacterial growth can be used to control microbial evolution.

Morphogenesis of bacterial colonies in polymeric environments

Many bacteria live in polymeric fluids, such as mucus, environmental polysaccharides, and extracellular polymers in biofilms. However, lab studies typically focus on cells in polymer-free fluids. Here, we show that interactions with polymers shape a fundamental feature of bacterial life-how they proliferate in space in multicellular colonies. Using experiments, we find that when polymer is sufficiently concentrated, cells generically and reversibly form large serpentine "cables" as they proliferate. By combining experiments with biophysical theory and simulations, we demonstrate that this distinctive form of colony morphogenesis arises from an interplay between polymer-induced entropic attraction between neighboring cells and their hindered ability to diffusely separate from each other in a viscous polymer solution. Our work thus reveals a pivotal role of polymers in sculpting proliferating bacterial colonies, with implications for how they interact with hosts and with the natural environment, and uncovers quantitative principles governing colony morphogenesis in such complex environments.

The role of filamentous matrix molecules in shaping the architecture and emergent properties of bacterial biofilms

Numerous bacteria naturally occur within spatially organised, multicellular communities called biofilms. Moreover, most bacterial infections proceed with biofilm formation, posing major challenges to human health. Within biofilms, bacterial cells are embedded in a primarily self-produced extracellular matrix, which is a defining feature of all biofilms. The biofilm matrix is a complex, viscous mixture primarily composed of polymeric substances such as polysaccharides, filamentous protein fibres, and extracellular DNA. The structured arrangement of the matrix bestows bacteria with beneficial emergent properties that are not displayed by planktonic cells, conferring protection against physical and chemical stresses, including antibiotic treatment. However, a lack of multi-scale information at the molecular level has prevented a better understanding of this matrix and its properties. Here, we review recent progress on the molecular characterisation of filamentous biofilm matrix components and their three-dimensional spatial organisation within biofilms.

Multi-population dissolution in confined active fluids

Autonomous out-of-equilibrium agents or cells in suspension are ubiquitous in biology and engineering. Turning chemical energy into mechanical stress, they generate activity in their environment, which may trigger spontaneous large-scale dynamics. Often, these systems are composed of multiple populations that may reflect the coexistence of multiple species, differing phenotypes, or chemically varying agents in engineered settings. Here, we present a new method for modeling such multi-population active fluids subject to confinement. We use a continuum multi-scale mean-field approach to represent each phase by its first three orientational moments and couple their evolution with those of the suspending fluid. The resulting coupled system is solved using a parallel adaptive level-set-based solver for high computational efficiency and maximal flexibility in the confinement geometry. Motivated by recent experimental work, we employ our method to study the spatiotemporal dynamics of confined bacterial suspensions and swarms dominated by fluid hydrodynamic effects. Our in silico explorations reproduce observed emergent collective patterns, including features of active dissolution in two-population active-passive swarms, with results clearly suggesting that hydrodynamic effects dominate dissolution dynamics. Our work lays the foundation for a systematic characterization and study of collective phenomena in natural or synthetic multi-population systems such as bacteria colonies, bird flocks, fish schools, colloid swimmers, or programmable active matter.

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.

Dynamic coexistence driven by physiological transitions in microbial communities

Microbial ecosystems are commonly modeled by fixed interactions between species in steady exponential growth states. However, microbes often modify their environments so strongly that they are forced out of the exponential state into stressed or non-growing states. Such dynamics are typical of ecological succession in nature and serial-dilution cycles in the laboratory. Here, we introduce a phenomenological model, the Community State model, to gain insight into the dynamic coexistence of microbes due to changes in their physiological states. Our model bypasses specific interactions (e.g., nutrient starvation, stress, aggregation) that lead to different combinations of physiological states, referred to collectively as "community states", and modeled by specifying the growth preference of each species along a global ecological coordinate, taken here to be the total community biomass density. We identify three key features of such dynamical communities that contrast starkly with steady-state communities: increased tolerance of community diversity to fast growth rates of species dominating different community states, enhanced community stability through staggered dominance of different species in different community states, and increased requirement on growth dominance for the inclusion of late-growing species. These features, derived explicitly for simplified models, are proposed here to be principles aiding the understanding of complex dynamical communities. Our model shifts the focus of ecosystem dynamics from bottom-up studies based on idealized inter-species interaction to top-down studies based on accessible macroscopic observables such as growth rates and total biomass density, enabling quantitative examination of community-wide characteristics.

Dynamic coexistence driven by physiological transitions in microbial communities

Microbial ecosystems are commonly modeled by fixed interactions between species in steady exponential growth states. However, microbes often modify their environments so strongly that they are forced out of the exponential state into stressed or non-growing states. Such dynamics are typical of ecological succession in nature and serial-dilution cycles in the laboratory. Here, we introduce a phenomenological model, the Community State model, to gain insight into the dynamic coexistence of microbes due to changes in their physiological states. Our model bypasses specific interactions (e.g., nutrient starvation, stress, aggregation) that lead to different combinations of physiological states, referred to collectively as "community states", and modeled by specifying the growth preference of each species along a global ecological coordinate, taken here to be the total community biomass density. We identify three key features of such dynamical communities that contrast starkly with steady-state communities: increased tolerance of community diversity to fast growth rates of species dominating different community states, enhanced community stability through staggered dominance of different species in different community states, and increased requirement on growth dominance for the inclusion of late-growing species. These features, derived explicitly for simplified models, are proposed here to be principles aiding the understanding of complex dynamical communities. Our model shifts the focus of ecosystem dynamics from bottom-up studies based on idealized inter-species interaction to top-down studies based on accessible macroscopic observables such as growth rates and total biomass density, enabling quantitative examination of community-wide characteristics.

Interplay of cell motility and self-secreted extracellular polymeric substance induced depletion effects on spatial patterning in a growing microbial colony

Reproducing bacteria self-organize to develop patterned biofilms in various conditions. Various factors contribute to the shaping of a multicellular bacterial organization. Here we investigate how motility force and self-secreted extracellular polymeric substances (EPS) influence bacterial cell aggregation, leading to phase-separated colonies using a particle-based/individual-based model. Our findings highlight the critical role of the interplay between motility force and depletion effects in regulating phase separation within a growing colony under far-from-equilibrium conditions. We observe that increased motility force hinders depletion-induced cell aggregation and phase segregation, necessitating a higher depletion effect for highly motile bacteria to undergo phase separation within a growing biofilm. We present a phase diagram illustrating the systematic variation of motility force and repulsive mechanical force, shedding light on the combined contributions of these two factors: self-propulsive motion and aggregation due to the depletion effect, resulting in the presence of small to large bacterial aggregates. Furthermore, our study reveals the dynamic nature of clustering, marked by changes in cluster size over time. Additionally, our findings suggest that differential dispersion among the components can lead to the localization of EPS at the periphery of a growing colony. Our study enhances the understanding of the collective dynamics of motile bacterial cells within a growing colony, particularly in the presence of a self-secreted polymer-driven depletion effect.

Mechano-induced homotypic patterned domain formation by monocytes

Matrix stiffness and corresponding mechano-signaling play indispensable roles in cellular phenotypes and functions. How tissue stiffness influences the behavior of monocytes, a major circulating leukocyte of the innate system, and how it may promote the emergence of collective cell behavior is less understood. Here, using tunable collagen-coated hydrogels of physiological stiffness, we show that human primary monocytes undergo a dynamic local phase separation to form highly regular, reversible, multicellular, multi-layered domains on soft matrix. Local activation of the β2 integrin initiates inter-cellular adhesion, while global soluble inhibitory factors maintain the steady state domain pattern over days. Patterned domain formation generated by monocytes is unique among other key immune cells, including macrophages, B cells, T cells, and NK cells. While inhibiting their phagocytic capability, domain formation promotes monocytes' survival. We develop a computational model based on the Cahn-Hilliard equation of phase separation, combined with a Turing mechanism of local activation and global inhibition suggested by our experiments, and provides experimentally validated predictions of the role of seeding density and both chemotactic and random cell migration on domain pattern formation. This work reveals that, unlike active matters, cells can generate complex cell phases by exploiting their mechanosensing abilities and combined short-range interactions and long-range signals to enhance their survival.

Distinct types of multicellular aggregates in Pseudomonas aeruginosa liquid cultures

Pseudomonas aeruginosa forms suspended multicellular aggregates when cultured in liquid media. These aggregates may be important in disease, and/or as a pathway to biofilm formation. The polysaccharide Psl and extracellular DNA (eDNA) have both been implicated in aggregation, but previous results depend strongly on the experimental conditions. Here we develop a quantitative microscopy-based method for assessing changes in the size distribution of suspended aggregates over time in growing cultures. For exponentially growing cultures of P. aeruginosa PAO1, we find that aggregation is mediated by cell-associated Psl, rather than by either eDNA or secreted Psl. These aggregates arise de novo within the culture via a growth process that involves both collisions and clonal growth, and Psl non-producing cells do not aggregate with producers. In contrast, we find that stationary phase (overnight) cultures contain a different type of multicellular aggregate, in which both eDNA and Psl mediate cohesion. Our findings suggest that the physical and biological properties of multicellular aggregates may be very different in early-stage vs late-stage bacterial cultures.

A subcellular biochemical model for T6SS dynamics reveals winning competitive strategies

The type VI secretion system (T6SS) is a broadly distributed interbacterial weapon that can be used to eliminate competing bacterial populations. Although unarmed target populations are typically used to study T6SS function in vitro, bacteria most likely encounter other T6SS-armed competitors in nature. However, the connection between subcellular details of the T6SS and the outcomes of such mutually lethal battles is not well understood. Here, we incorporate biological data derived from natural competitors of Vibrio fischeri light organ symbionts to build a biochemical model for T6SS at the single-cell level, which we then integrate into an agent-based model (ABM). Using the ABM, we isolate and experiment with strain-specific physiological differences between competitors in ways not possible with biological samples to identify winning strategies for T6SS-armed populations. Through in vitro experiments, we discover that strain-specific differences exist in T6SS activation speed. ABM simulations corroborate that faster activation is dominant in determining survival during competition. Once competitors are fully activated, the energy required for T6SS creates a tipping point where increased weapon building and firing becomes too costly to be advantageous. Through ABM simulations, we identify the threshold where this transition occurs in the T6SS parameter space. We also find that competitive outcomes depend on the geometry of the battlefield: unarmed target cells survive at the edges of a range expansion where unlimited territory can be claimed. Alternatively, competitions within a confined space, much like the light organ crypts where natural V. fischeri compete, result in the rapid elimination of the unarmed population.

Simulation of various biofilm fractal morphologies by agent-based model

Biofilms are clusters of bacteria wrapped in extracellular matrix and polymers. The study of biofilm morphological transformation has been around for a long time and has attracted widespread attention. In this paper, we present a model for biofilm growth based on the interaction force, in which bacteria are treated as tiny particles and locations of particles are updated by calculating the repulsive forces among particles. We adapt a continuity equation to indicate nutrient concentration variation in the substrate. Based on the above, we study the morphological transformation of biofilms. We find that nutrient concentration and nutrient diffusion rate dominate different biofilm morphological transition processes, in which biofilms would grow into fractal morphology under the conditions of low nutrient concentration and nutrient diffusivity. At the same time, we expand our model by introducing a second particle to mimic extracellular polymeric substances (EPS) in biofilms. We find that the interaction between different particles can lead to phase separation patterns between cells and EPSs, and the adhesion effect of EPS can attenuate this phenomenon. In contrast to single particle system models, branches are inhibited due to EPS filling in dual particle system models, and this invalidation is boosted by the enhancement of the depletion effect.

Collective protection against the type VI secretion system in bacteria

Bacteria commonly face attacks from other strains using the type VI secretion system (T6SS), which acts like a molecular speargun to stab and intoxicate competitors. Here we show how bacteria can work together to collectively defend themselves against these attacks. This project began with an outreach activity: while developing an online computer game of bacterial warfare, we noticed that one strategist ("Slimy") that made extracellular polymeric substances (EPS) was able to resist attacks from another strategist that employed the T6SS ("Stabby"). This observation motivated us to model this scenario more formally, using dedicated agent-based simulations. The model predicts that EPS production can serve as a collective defence mechanism, which protects both producing cells and neighbouring cells that do not make EPS. We then tested our model with a synthetic community that contains a T6SS-wielding attacker (Acinetobacter baylyi), and two T6SS-sensitive target strains (Escherichia coli) that either secrete EPS, or not. As predicted by our modelling, we find that the production of EPS leads to collective protection against T6SS attacks, where EPS producers protect each other and nearby non-producers. We identify two processes that explain this protection: EPS sharing between cells and a second general mechanism whereby groups of resistant cells shield susceptible cells, which we call "flank protection". Our work shows how EPS-producing bacteria can work together to defend themselves from the type VI secretion system.

Active layer dynamics drives a transition to biofilm fingering

The emergence of spatial organisation in biofilm growth is one of the most fundamental topics in biofilm biophysics and microbiology. It has long been known that growing biofilms can adopt smooth or rough interface morphologies, depending on the balance between nutrient supply and microbial growth; this 'fingering' transition has been linked with the average width of the 'active layer' of growing cells at the biofilm interface. Here we use long-time individual-based simulations of growing biofilms to investigate in detail the driving factors behind the biofilm-fingering transition. We show that the transition is associated with dynamical changes in the active layer. Fingering happens when gaps form in the active layer, which can cause local parts of the biofilm interface to pin, or become stationary relative to the moving front. Pinning can be transient or permanent, leading to different biofilm morphologies. By constructing a phase diagram for the transition, we show that the controlling factor is the magnitude of the relative fluctuations in the active layer thickness, rather than the active layer thickness per se. Taken together, our work suggests a central role for active layer dynamics in controlling the pinning of the biofilm interface and hence biofilm morphology.

Structure-Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities─and their impacts on environmental processes─requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca-Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ∼0.6 nm between alginate strands.

A mechanistic understanding of microcolony morphogenesis: coexistence of mobile and sessile aggregates

Most bacteria in the natural environment self-organize into collective phases such as cell clusters, swarms, patterned colonies, or biofilms. Several intrinsic and extrinsic factors, such as growth, motion, and physicochemical interactions, govern the occurrence of different phases and their coexistence. Hence, predicting the conditions under which a collective phase emerges due to individual-level interactions is crucial. Here we develop a particle-based biophysical model of bacterial cells and self-secreted extracellular polymeric substances (EPS) to decipher the interplay of growth, motility-mediated dispersal, and mechanical interactions during microcolony morphogenesis. We show that the microcolony dynamics and architecture significantly vary depending upon the heterogeneous EPS production. In particular, microcolony shows the coexistence of both motile and sessile aggregates rendering a transition towards biofilm formation. We identified that the interplay of differential dispersion and the mechanical interactions among the components of the colony determines the fate of the colony morphology. Our results provide a significant understanding of the mechano-self-regulation during biofilm morphogenesis and open up possibilities of designing experiments to test the predictions.

New understanding of multidrug efflux and permeation in antibiotic resistance, persistence, and heteroresistance

Antibiotics effective against Gram-negative ESKAPE pathogens are a critical area of unmet need. Infections caused by these pathogens are not only difficult to treat but finding new therapies to overcome Gram-negative resistance is also a challenge. There are not enough antibiotics in development that target the most dangerous pathogens and there are not enough novel drugs in the pipeline. The major obstacle in the antibiotic discovery pipeline is the lack of understanding of how to breach antibiotic permeability barriers of Gram-negative pathogens. These barriers are created by active efflux pumps acting across both the inner and the outer membranes. Overproduction of efflux pumps alone or together with either modification of the outer membrane or antibiotic-inactivating enzymes and target mutations contribute to clinical levels of antibiotics resistance. Recent efforts have generated significant advances in the rationalization of compound efflux and permeation across the cell envelopes of Gram-negative pathogens. Combined with earlier studies and novel mathematical models, these efforts have led to a multilevel understanding of how antibiotics permeate these barriers and how multidrug efflux and permeation contribute to the development of antibiotic resistance and heteroresistance. Here, we discuss the new developments in this area.

'Phase transitions' in bacteria - From structural transitions in free living bacteria to phenotypic transitions in bacteria within biofilms

Phase transitions are common in inanimate systems and have been studied extensively in natural sciences. Less explored are the rich transitions that take place at the micro- and nano-scales in biological systems. In conventional phase transitions, large-scale properties of the media change discontinuously in response to continuous changes in external conditions. Such changes play a significant role in the dynamic behaviours of organisms. In this review, we focus on some transitions in both free-living and biofilms of bacteria. Particular attention is paid to the transitions in the flagellar motors and filaments of free-living bacteria, in cellular gene expression during the biofilm growth, in the biofilm morphology transitions during biofilm expansion, and in the cell motion pattern transitions during the biofilm formation. We analyse the dynamic characteristics and biophysical mechanisms of these phase transition phenomena and point out the parallels between these transitions and conventional phase transitions. We also discuss the applications of some theoretical and numerical methods, established for conventional phase transitions in inanimate systems, in bacterial biofilms.

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.

Shared biophysical mechanisms determine early biofilm architecture development across different bacterial species

Bacterial biofilms are among the most abundant multicellular structures on Earth and play essential roles in a wide range of ecological, medical, and industrial processes. However, general principles that govern the emergence of biofilm architecture across different species remain unknown. Here, we combine experiments, simulations, and statistical analysis to identify shared biophysical mechanisms that determine early biofilm architecture development at the single-cell level, for the species Vibrio cholerae, Escherichia coli, Salmonella enterica, and Pseudomonas aeruginosa grown as microcolonies in flow chambers. Our data-driven analysis reveals that despite the many molecular differences between these species, the biofilm architecture differences can be described by only 2 control parameters: cellular aspect ratio and cell density. Further experiments using single-species mutants for which the cell aspect ratio and the cell density are systematically varied, and mechanistic simulations show that tuning these 2 control parameters reproduces biofilm architectures of different species. Altogether, our results show that biofilm microcolony architecture is determined by mechanical cell-cell interactions, which are conserved across different species.

Fractal morphology facilitates Bacillus subtilis biofilm growth

In the late stage of Bacillus subtilis biofilm growth, to adapt the extremely nutrient-lacking environment, the biofilm edge grows into a complex branching structure, which allows the biofilm to expand outward at a faster speed, comparing to the expansion speed of the biofilm edge without branching structure. The fractal analysis shows that the fractal dimension (Fd) decreases along the radius in the biofilm branching structure, as shown in Figs. 1d and 3a. The variation of Fd along the radius is not monotonic, which is because of the texture evolution induced by the bacterial clusters' movement. By using the wide field stereomicroscope and image analysis, we find that the ridges in the mature branching structure are composed of inactive substances, and most of the bacterial clusters move through the valleys. Further analysis shows that bacterial clusters move to the area with the high Succolarity (Suc) value.

Engineered Living Hydrogels

Living biological systems, ranging from single cells to whole organisms, can sense, process information, and actuate in response to changing environmental conditions. Inspired by living biological systems, engineered living cells and nonliving matrices are brought together, which gives rise to the technology of engineered living materials. By designing the functionalities of living cells and the structures of nonliving matrices, engineered living materials can be created to detect variability in the surrounding environment and to adjust their functions accordingly, thereby enabling applications in health monitoring, disease treatment, and environmental remediation. Hydrogels, a class of soft, wet, and biocompatible materials, have been widely used as matrices for engineered living cells, leading to the nascent field of engineered living hydrogels. Here, the interactions between hydrogel matrices and engineered living cells are described, focusing on how hydrogels influence cell behaviors and how cells affect hydrogel properties. The interactions between engineered living hydrogels and their environments, and how these interactions enable versatile applications, are also discussed. Finally, current challenges facing the field of engineered living hydrogels for their applications in clinical and environmental settings are highlighted.

Permeation Increases Biofilm Development in Nanofiltration Membranes Operated with Varying Feed Water Phosphorous Concentrations

Nutrient limitation has been proposed as a biofouling control strategy for membrane systems. However, the impact of permeation on biofilm development under phosphorus-limited and enriched conditions is poorly understood. This study analyzed biofilm development in membrane fouling simulators (MFSs) with and without permeation supplied with water varying dosed phosphorus concentrations (0 and 25 μg P·L−1). The MFSs operated under permeation conditions were run at a constant flux of 15.6 L·m2·h−1 for 4.7 days. Feed channel pressure drop, transmembrane pressure, and flux were used as performance indicators. Optical coherence tomography (OCT) images and biomass quantification were used to analyze the developed biofilms. The total phosphorus concentration that accumulated on the membrane and spacer was quantified by using microwave digestion and inductively coupled plasma atomic emission spectroscopy (ICP-OES). Results show that permeation impacts biofilm development depending on nutrient condition with a stronger impact at low P concentration (pressure drop increase: 282%; flux decline: 11%) compared to a higher P condition (pressure drop increase: 206%; flux decline: 2%). The biofilm that developed at 0 μg P·L−1 under permeation conditions resulted in a higher performance decline due to biofilm localization and spread in the MFS. A thicker biofilm developed on the membrane for biofilms grown at 0 μg P·L−1 under permeation conditions, causing a stronger effect on flux decline (11%) compared to non-permeation conditions (5%). The difference in the biofilm thickness on the membrane was attributed to a higher phosphorus concentration in the membrane biofilm under permeation conditions. Permeation has an impact on biofilm development and, therefore, should not be excluded in biofouling studies.

Chemotactic smoothing of collective migration

Collective migration—the directed, coordinated motion of many self-propelled agents—is a fascinating emergent behavior exhibited by active matter with functional implications for biological systems. However, how migration can persist when a population is confronted with perturbations is poorly understood. Here, we address this gap in knowledge through studies of bacteria that migrate via directed motion, or chemotaxis, in response to a self-generated nutrient gradient. We find that bacterial populations autonomously smooth out large-scale perturbations in their overall morphology, enabling the cells to continue to migrate together. This smoothing process arises from spatial variations in the ability of cells to sense and respond to the local nutrient gradient—revealing a population-scale consequence of the manner in which individual cells transduce external signals. Altogether, our work provides insights to predict, and potentially control, the collective migration and morphology of cellular populations and diverse other forms of active matter.

Flocks of birds, schools of fish and herds of animals are all good examples of collective migration, where individuals co-ordinate their behavior to improve survival. This process also happens on a cellular level; for example, when bacteria consume a nutrient in their surroundings, they will collectively move to an area with a higher concentration of food via a process known as chemotaxis.

Several studies have examined how disturbing collective migration can cause populations to fall apart. However, little is known about how groups withstand these interferences. To investigate, Bhattacharjee, Amchin, Alert et al. studied bacteria called Escherichia coli as they moved through a gel towards nutrients.

The E. coli were injected into the gel using a three-dimensional printer, which deposited the bacteria into a wiggly shape that forces the cells apart, making it harder for them to move as a collective group. However, as the bacteria migrated through the gel, they smoothed out the line and gradually made it straighter so they could continue to travel together over longer distances.

Computer simulations revealed that this smoothing process is achieved by differences in how the cells respond to local nutrient levels based on their position. Bacteria towards the front of the group are exposed to more nutrients, causing them to become oversaturated and respond less effectively to the nutrient gradient. As a result, they move more slowly, allowing the cells behind them to eventually catch-up.

These findings reveal a general mechanism in which limitations in how individuals sense and respond to an external signal (in this case local nutrient concentrations) allows them to continue migrating together. This mechanism may apply to other systems that migrate via chemotaxis, as well as groups whose movement is directed by different external factors, such as temperature and light intensity.

Validated In Silico Model for Biofilm Formation in Escherichia coli

Using Escherichia coli as the representative biofilm former, we report here the development of an in silico model built by simulating events that transform a free-living bacterial entity into self-encased multicellular biofilms. Published literature on ∼300 genes associated with pathways involved in biofilm formation was curated, static maps were created, and suitably interconnected with their respective metabolites using ordinary differential equations. Precise interplay of genetic networks that regulate the transitory switching of bacterial growth pattern in response to environmental changes and the resultant multicomponent synthesis of the extracellular matrix were appropriately represented. Subsequently, the in silico model was analyzed by simulating time-dependent changes in the concentration of components by using the R and python environment. The model was validated by simulating and verifying the impact of key gene knockouts (KOs) and systematic knockdowns on biofilm formation, thus ensuring the outcomes were comparable with the reported literature. Similarly, specific gene KOs in laboratory and pathogenic E. coli were constructed and assessed. MiaA, YdeO, and YgiV were found to be crucial in biofilm development. Furthermore, qRT-PCR confirmed the elevation of expression in biofilm-forming clinical isolates. Findings reported in this study offer opportunities for identifying biofilm inhibitors with applications in multiple industries. The application of this model can be extended to the health care sector specifically to develop novel adjunct therapies that prevent biofilms in medical implants and reduce emergence of biofilm-associated resistant polymicrobial-chronic infections. The in silico framework reported here is open source and accessible for further enhancements.

Mechanical forces drive a reorientation cascade leading to biofilm self-patterning

In growing active matter systems, a large collection of engineered or living autonomous units metabolize free energy and create order at different length scales as they proliferate and migrate collectively. One such example is bacterial biofilms, surface-attached aggregates of bacterial cells embedded in an extracellular matrix that can exhibit community-scale orientational order. However, how bacterial growth coordinates with cell-surface interactions to create distinctive, long-range order during biofilm development remains elusive. Here we report a collective cell reorientation cascade in growing Vibrio cholerae biofilms that leads to a differentially ordered, spatiotemporally coupled core-rim structure reminiscent of a blooming aster. Cell verticalization in the core leads to a pattern of differential growth that drives radial alignment of the cells in the rim, while the growing rim generates compressive stresses that expand the verticalized core. Such self-patterning disappears in nonadherent mutants but can be restored through opto-manipulation of growth. Agent-based simulations and two-phase active nematic modeling jointly reveal the strong interdependence of the driving forces underlying the differential ordering. Our findings offer insight into the developmental processes that shape bacterial communities and provide ways to engineer phenotypes and functions in living active matter.

Cluster and conquer: the morphodynamics of invasion of a compliant substrate by active rods

The colonisation of a soft passive material by motile cells such as bacteria is common in biology. The resulting colonies of the invading cells are often observed to exhibit intricate patterns whose morphology and dynamics can depend on a number of factors, particularly the mechanical properties of the substrate and the motility of the individual cells. We use simulations of a minimal 2D model of self-propelled rods moving through a passive compliant medium consisting of particles that offer elastic resistance before being plastically displaced from their equilibrium positions. It is observed that the clustering of active (self-propelled) particles is crucial for understanding the morphodynamics of colonisation. Clustering enables motile colonies to spread faster than they would have as isolated particles. The colonisation rate depends non-monotonically on substrate stiffness with a distinct maximum at a non-zero value of substrate stiffness. This is observed to be due to a change in the morphology of clusters. Furrow networks created by the active particles have a fractal-like structure whose dimension varies systematically with substrate stiffness but is less sensitive to particle activity. The power-law growth exponent of the furrowed area is smaller than unity, suggesting that, to sustain such extensive furrow networks, colonies must regulate their overall growth rate.

Mechanistic underpinning of cell aspect ratio-dependent emergent collective motions in swarming bacteria

Self-propelled bacteria can exhibit a large variety of non-equilibrium self-organized phenomena. Swarming is one such fascinating dynamical scenario where a number of motile individuals group into dynamical clusters and move in synchronized flows and vortices. While precedent investigations into rod-like particles confirm that an increased aspect-ratio promotes alignment and order, recent experimental studies in bacteria Bacillus subtilis show a non-monotonic dependence of the cell-aspect ratio on their swarming motion. Here, by computer simulations of an agent-based model of self-propelled, mechanically interacting, rod-shaped bacteria under overdamped conditions, we explore the collective dynamics of a bacterial swarm subjected to a variety of cell-aspect ratios. When modeled with an identical self-propulsion speed across a diverse range of cell aspect ratios, simulations demonstrate that both shorter and longer bacteria exhibit slow dynamics whereas the fastest speed is obtained at an intermediate aspect ratio. Our investigation highlights that the origin of this observed non-monotonic trend of bacterial speed and vorticity with the cell-aspect ratio is rooted in the cell-size dependence of motility force. The swarming features remain robust for a wide range of surface density of the cells, whereas asymmetry in friction attributes a distinct effect. Our analysis identifies that at an intermediate aspect ratio, an optimum cell size and motility force promote alignment, which reinforces the mechanical interactions among neighboring cells leading to the overall fastest motion. Mechanistic underpinning of the collective motions reveals that it is a joint venture of the short-range repulsive and the size-dependent motility forces, which determines the characteristics of swarming.

Causes and consequences of pattern diversification in a spatially self-organizing microbial community

Surface-attached microbial communities constitute a vast amount of life on our planet. They contribute to all major biogeochemical cycles, provide essential services to our society and environment, and have important effects on human health and disease. They typically consist of different interacting genotypes that arrange themselves non-randomly across space (referred to hereafter as spatial self-organization). While spatial self-organization is important for the functioning, ecology, and evolution of these communities, the underlying determinants of spatial self-organization remain unclear. Here, we performed a combination of experiments, statistical modeling, and mathematical simulations with a synthetic cross-feeding microbial community consisting of two isogenic strains. We found that two different patterns of spatial self-organization emerged at the same length and time scales, thus demonstrating pattern diversification. This pattern diversification was not caused by initial environmental heterogeneity or by genetic heterogeneity within populations. Instead, it was caused by nongenetic heterogeneity within populations, and we provide evidence that the source of this nongenetic heterogeneity is local differences in the initial spatial positionings of individuals. We further demonstrate that the different patterns exhibit different community-level properties; namely, they have different expansion speeds. Together, our results demonstrate that pattern diversification can emerge in the absence of initial environmental heterogeneity or genetic heterogeneity within populations and can affect community-level properties, thus providing novel insights into the causes and consequences of microbial spatial self-organization.

Polymer-induced microcolony compaction in early biofilms: A computer simulation study

Microscopic organisms, such as bacteria, have the ability of colonizing surfaces and developing biofilms that can determine diseases and infections. Most bacteria secrete a significant amount of extracellular polymer substances that are relevant for biofilm stabilization and growth. In this work, we apply computer simulation and perform experiments to investigate the impact of polymer size and concentration on early biofilm formation and growth. We observe as bacterial cells formed loose, disorganized clusters whenever the effect of diffusion exceeded that of cell growth and division. Addition of model polymeric molecules induced particle self-assembly and aggregation to form compact clusters in a polymer size- and concentration-dependent fashion. We also find that large polymer size or concentration lead to the development of intriguing stripe-like and dendritic colonies. The results obtained by Brownian dynamic simulation closely resemble the morphologies that we experimentally observe in biofilms of a Pseudomonas Putida strain with added polymers. The analysis of the Brownian dynamic simulation results suggests the existence of a threshold polymer concentration that distinguishes between two growth regimes. Below this threshold, the main force driving polymer-induced compaction is the hindrance of bacterial cell diffusion, while collective effects play a minor role. Above this threshold, especially for large polymers, polymer-induced compaction is a collective phenomenon driven by depletion forces. Well above this concentration threshold, severely limited diffusion drives the formation of filaments and dendritic colonies.

Cellular and developmental basis of avian structural coloration

Vivid structural colors in birds are a conspicuous and vital part of their phenotype. They are produced by a rich diversity of integumentary photonic nanostructures in skin and feathers. Unlike pigmentary coloration, whose genetic basis is being elucidated, little is known regarding the pathways underpinning organismal structural coloration. Here, we review available data on the development of avian structural colors. In particular, feather photonic nanostructures are understood to be intracellularly self-assembled by physicochemical forces typically seen in soft colloidal systems. We identify promising avenues for future research that can address current knowledge gaps, which are also highly relevant for the sustainable engineering of advanced bioinspired and biomimetic materials.

How Physical Interactions Shape Bacterial Biofilms

Biofilms are structured communities formed by a single or multiple microbial species. Within biofilms, bacteria are embedded into extracellular matrix, allowing them to build macroscopic objects. Biofilm structure can respond to environmental changes such as the presence of antibiotics or predators. By adjusting expression levels of surface and extracellular matrix components, bacteria tune cell-to-cell interactions. One major challenge in the field is the fact that these components are very diverse among different species. Deciphering how physical interactions within biofilms are affected by changes in gene expression is a promising approach to obtaining a more unified picture of how bacteria modulate biofilms. This review focuses on recent advances in characterizing attractive and repulsive forces between bacteria in correlation with biofilm structure, dynamics, and spreading. How bacteria control physical interactions to maximize their fitness is an emerging theme.

Molecular and Genetic Mechanisms That Mediate Transmission of Yersinia pestis by Fleas

The ability to cause plague in mammals represents only half of the life history of Yersinia pestis. It is also able to colonize and produce a transmissible infection in the digestive tract of the flea, its insect host. Parallel to studies of the molecular mechanisms by which Y. pestis is able to overcome the immune response of its mammalian hosts, disseminate, and produce septicemia, studies of Y. pestis-flea interactions have led to the identification and characterization of important factors that lead to transmission by flea bite. Y. pestis adapts to the unique conditions in the flea gut by altering its metabolic physiology in ways that promote biofilm development, a common strategy by which bacteria cope with a nutrient-limited environment. Biofilm localization to the flea foregut disrupts normal fluid dynamics of blood feeding, resulting in regurgitative transmission. Many of the important genes, regulatory pathways, and molecules required for this process have been identified and are reviewed here.

Haloferax volcanii Immersed Liquid Biofilms Develop Independently of Known Biofilm Machineries and Exhibit Rapid Honeycomb Pattern Formation

The ability to form biofilms is shared by many microorganisms, including archaea. Cells in a biofilm are encased in extracellular polymeric substances that typically include polysaccharides, proteins, and extracellular DNA, conferring protection while providing a structure that allows for optimal nutrient flow. In many bacteria, flagella and evolutionarily conserved type IV pili are required for the formation of biofilms on solid surfaces or floating at the air-liquid interface of liquid media. Similarly, in many archaea it has been demonstrated that type IV pili and, in a subset of these species, archaella are required for biofilm formation on solid surfaces. Additionally, in the model archaeon Haloferax volcanii, chemotaxis and AglB-dependent glycosylation play important roles in this process. H. volcanii also forms immersed biofilms in liquid cultures poured into petri dishes. This study reveals that mutants of this haloarchaeon that interfere with the biosynthesis of type IV pili or archaella, as well as a chemotaxis-targeting transposon and aglB deletion mutants, lack obvious defects in biofilms formed in liquid cultures. Strikingly, we have observed that these liquid-based biofilms are capable of rearrangement into honeycomb-like patterns that rapidly form upon removal of the petri dish lid, a phenomenon that is not dependent on changes in light or oxygen concentration but can be induced by controlled reduction of humidity. Taken together, this study demonstrates that H. volcanii requires novel, unidentified strategies for immersed liquid biofilm formation and also exhibits rapid structural rearrangements.IMPORTANCE This first molecular biological study of archaeal immersed liquid biofilms advances our basic biological understanding of the model archaeon Haloferax volcanii Data gleaned from this study also provide an invaluable foundation for future studies to uncover components required for immersed liquid biofilms in this haloarchaeon and also potentially for liquid biofilm formation in general, which is poorly understood compared to the formation of biofilms on surfaces. Moreover, this first description of rapid honeycomb pattern formation is likely to yield novel insights into the underlying structural architecture of extracellular polymeric substances and cells within immersed liquid biofilms.

Mechanics of biofilms formed of bacteria with fimbriae appendages

Gram-negative bacteria, as well as some Gram-positive bacteria, possess hair-like appendages known as fimbriae, which play an important role in adhesion of the bacteria to surfaces or to other bacteria. Unlike the sex pili or flagellum, the fimbriae are quite numerous, with of order 1000 fimbriae appendages per bacterial cell. In this paper, a recently developed hybrid model for bacterial biofilms is used to examine the role of fimbriae tension force on the mechanics of bacterial biofilms. Each bacterial cell is represented in this model by a spherocylindrical particle, which interact with each other through collision, adhesion, lubrication force, and fimbrial force. The bacterial cells absorb water and nutrients and produce extracellular polymeric substance (EPS). The flow of water and EPS, and nutrient diffusion within these substances, is computed using a continuum model that accounts for important effects such as osmotic pressure gradient, drag force on the bacterial cells, and viscous shear. The fimbrial force is modeled using an outer spherocylinder capsule around each cell, which can transmit tensile forces to neighboring cells with which the fimbriae capsule collides. We find that the biofilm structure during the growth process is dominated by a balance between outward drag force on the cells due to the EPS flow away from the bacterial colony and the inward tensile fimbrial force acting on chains of cells connected by adhesive fimbriae appendages. The fimbrial force also introduces a large rotational motion of the cells and disrupts cell alignment caused by viscous torque imposed by the EPS flow. The current paper characterizes the competing effects of EPS drag and fimbrial force using a series of computations with different values of the ratio of EPS to bacterial cell production rate and different numbers of fimbriae per cell.

Bacterial biomechanics-From individual behaviors to biofilm and the gut flora

Bacteria inhabit a variety of locations and play important roles in the environment and health. Our understanding of bacterial biomechanics has improved markedly in the last decade and has revealed that biomechanics play a significant role in microbial biology. The obtained knowledge has enabled investigation of complex phenomena, such as biofilm formation and the dynamics of the gut flora. A bottom-up strategy, i.e., from the cellular to the macroscale, facilitates understanding of macroscopic bacterial phenomena. In this Review, we first cover the biomechanics of individual bacteria in the bulk liquid and on surfaces as the base of complex phenomena. The collective behaviors of bacteria in simple environments are next introduced. We then introduce recent advances in biofilm biomechanics, in which adhesion force and the flow environment play crucial roles. We also review transport phenomena in the intestine and the dynamics of the gut flora, focusing on that in zebrafish. Finally, we provide an overview of the future prospects for the field.

Interaction-dependent effects of surface structure on microbial spatial self-organization

Surface-attached microbial communities consist of different cell types that, at least to some degree, organize themselves non-randomly across space (referred to as spatial self-organization). While spatial self-organization can have important effects on the functioning, ecology and evolution of communities, the underlying determinants of spatial self-organization remain unclear. Here, we hypothesize that the presence of physical objects across a surface can have important effects on spatial self-organization. Using pairs of isogenic strains of Pseudomonas stutzeri, we performed range expansion experiments in the absence or presence of physical objects and quantified the effects on spatial self-organization. We demonstrate that physical objects create local deformities along the expansion frontier, and these deformities increase in magnitude during range expansion. The deformities affect the densities of interspecific boundaries and diversity along the expansion frontier, and thus affect spatial self-organization, but the effects are interaction-dependent. For competitive interactions that promote sectorized patterns of spatial self-organization, physical objects increase the density of interspecific boundaries and diversity. By contrast, for cross-feeding interactions that promote dendritic patterns, they decrease the density of interspecific boundaries and diversity. These qualitatively different outcomes are probably caused by fundamental differences in the orientations of the interspecific boundaries. Thus, in order to predict the effects of physical objects on spatial self-organization, information is needed regarding the interactions present within a community and the general geometric shapes of spatial self-organization that emerge from those interactions. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Hybrid Model of Bacterial Biofilm Growth

Bacterial biofilms play a critical role in environmental processes, water treatment, human health, and food processing. They exhibit highly complex dynamics due to the interactions between the bacteria and the extracellular polymeric substance (EPS), water, and nutrients and minerals that make up the biofilm. We present a hybrid computational model in which the dynamics of discrete bacterial cells are simulated within a multiphase continuum, consisting of EPS and water as separate interacting phases, through which nutrients and minerals diffuse. Bacterial cells in our model consume water and nutrients in order to grow, divide, and produce EPS. Consequently, EPS flows outward from the bacterial colony, while water flows inward. The model predicts bacterial colony formation as a treelike structure. The distribution of bacterial growth and EPS production is found to be sensitive to the pore spacing between bacteria and the consumption of nutrients within the bacterial colony. Forces that are sometimes neglected in biofilm simulations, such as lubrication force between nearby bacterial cells and osmotic (swelling) pressure force resulting from gradients in EPS concentration, are observed to have an important effect on biofilm growth via their influence on bacteria pore spacing and associated water/nutrient percolation into the bacterial colony.

Spatiotemporal pattern formation in E. coli biofilms explained by a simple physical energy balance

While the biofilm growth mode conveys notable thriving advantages to bacterial populations, the mechanisms of biofilm formation are still strongly debated. Here, we investigate the remarkable spontaneous formation of regular spatial patterns during the growth of an Escherichia coli biofilm. These patterns reported here appear with non-motile bacteria, which excludes both chemotactic origins and other motility-based ones. We demonstrate that a minimal physical model based on phase separation describes them well. To confirm the predictive capacity of our model, we tune the cell-cell and cell-surface interactions using cells expressing different surface appendages. We further explain how F pilus-bearing cells enroll their wild type kindred, poorly piliated, into their typical pattern when mixed together. This work supports the hypothesis that purely physicochemical processes, such as the interplay of cell-cell and cell-surface interactions, can drive the emergence of a highly organized spatial structure that is potentially decisive for community fate and for biological functions.

Mono- to Multilayer Transition in Growing Bacterial Colonies

The transition from monolayers to multilayered structures in bacterial colonies is a fundamental step in biofilm development. Observed across different morphotypes and species, this transition is triggered within freely growing bacterial microcolonies comprising a few hundred cells. Using a combination of numerical simulations and analytical modeling, here we demonstrate that this transition originates from the competition between growth-induced in-plane active stresses and vertical restoring forces, due to the cell-substrate interactions. Using a simple chainlike colony of laterally confined cells, we show that the transition sets when individual cells become unstable to rotations; thus it is localized and mechanically deterministic. Asynchronous cell division renders the process stochastic, so that all the critical parameters that control the onset of the transition are continuously distributed random variables. Here we demonstrate that the occurrence of the first division in the colony can be approximated as a Poisson process in the limit of large cell numbers. This allows us to approximately calculate the probability distribution function of the position and time associated with the first extrusion. The rate of such a Poisson process can be identified as the order parameter of the transition, thus highlighting its mixed deterministic-stochastic nature.

Arrested fungal biofilms as low-modulus structural bio-composites: Water holds the key

Biofilms are self-assembling structures consisting of rigid microbial cells embedded in a soft biopolymeric extracellular matrix (ECM), and have been commonly viewed as being detrimental to health and equipment. In this work, we show that biofilms formed by a non-pathogenic fungus Neurospora discreta, are fungal bio-composites (FBCs) that can be directed to self-organize through active stresses to achieve specific properties. We induced active stresses by systematically varying the agitation rate during the growth of FBCs. By growing FBCs that are strong enough to be conventionally tensile loaded, we find that as agitation rate increases, the elongation strain at which the FBCs break, increases linearly, and their elastic modulus correspondingly decreases. Using results from microstructural imaging and thermogravimetry, we rationalize that agitation increases the production of ECM, which concomitantly increases the water content of agitated FBCs up to 250% more than un-agitated FBCs. Water held in the nanopores of the ECM acts a plasticizer and controls the ductility of FBCs in close analogy with polyelectrolyte complexes. This paradigm shift in viewing biofilms as bio-composites opens up the possibility for their use as sustainable, biodegradable, low-modulus structural materials.

Spatial Organization of Expanding Bacterial Colonies Is Affected by Contact-Dependent Growth Inhibition

Identifying how microbes are able to manipulate, survive, and thrive in complex multispecies communities has expanded our understanding of how microbial ecosystems impact human health and the environment. The ability of bacteria to negatively affect neighbors, through explicit toxin delivery systems, provides them with an opportunity to manipulate the composition of growing microbial communities. Contact-dependent inhibition (CDI) systems (a Type Vb secretion system) are a distinct subset of competition systems whose contribution to shaping the development of spatially structured bacterial communities are yet to be fully understood. Here, we compare the impact of different CDI systems, at both the single-cell and population level, to determine the key drivers of CDI-mediated competition within spatially structured bacterial populations. Through an iterative approach using both an Escherichia coli experimental system and computational modeling, we show that CDI systems have subtle and system-specific effects at the single-cell level, generating single-cell-wide boundaries between CDI-expressing inhibitor cells and their neighboring targets. Despite the subtle effects of CDI at a single-cell level, CDI systems greatly diminished the ability of susceptible targets to expand their range during colony growth. The inoculum density of the population, together with the CDI system-specific variables of the speed of inhibition after contact and biological cost of CDI, strongly affects CDI-mediated competition. In contrast, the magnitude of the toxin-induced growth retardation of target cells only weakly impacts the composition of the population. Our work reveals how distinct CDI systems can differentially affect the composition and spatial arrangement of bacterial populations.

•

CDI causes subtle growth inhibition in a subset of contacted target cells

•

Model describes and predicts observed effects on spatial distribution of strains

•

CDI facilitates success of inhibitor strain increasing population patch size

•

A CDI system’s inhibition rate dominates toxicity in driving competition outcome

In Silico Modeling of Biofilm Formation by Nontypeable Haemophilus influenzae In Vivo

Biofilms formed by nontypeable Haemophilus influenzae (NTHI) bacteria play an important role in multiple respiratory tract diseases. Visual inspection of the morphology of biofilms formed during chronic infections shows distinct differences from biofilms formed in vitro To better understand these differences, we analyzed images of NTHI biofilms formed in the middle ears of Chinchilla lanigera and developed an in silico agent-based model of the formation of NTHI biofilms in vivo We found that, as in vitro, NTHI bacteria are organized in self-similar patterns; however, the sizes of NTHI clusters in vivo are more than 10-fold smaller than their in vitro counterparts. The agent-based model reproduced these patterns and suggested that smaller clusters occur due to elimination of planktonic NTHI cells by the host responses. Estimation of model parameters by fitting simulation results to imaging data showed that the effects of several processes in the model change during the course of the infection.IMPORTANCE Multiple respiratory illnesses are associated with formation of biofilms within the human airway by NTHI. However, a substantial amount of our understanding of the mechanisms that underlie NTHI biofilm formation is obtained from in vitro studies. Our in silico model that describes biofilm formation by NTHI within the middle ears of Chinchilla lanigera will help isolate processes potentially responsible for the differences between the morphologies of biofilms formed in vivo versus those formed in vitro Thus, the in silico model can be used to glean mechanisms that underlie biofilm formation in vivo and connect those mechanisms to those obtained from in vitro experiments. The in silico model developed here can be extended to investigate potential roles of specific host responses (e.g., mucociliary clearance) on NTHI biofilm formation in vivo The developed computational tools can also be used to analyze and describe biofilm formation by other bacterial species in vivo.

Viscoelasticity variation in a biofilm-mediated Bacillus subtilis suspension induced by adding polyethylene glycol

Recent experiments show that synthetic polymers can influence the degree of microbial aggregation and the rheological properties of bacterial suspensions, the study of which can help us control biofilm formation. In this article, we add polyethylene glycol (PEG) with various molecular weights and concentrations into two types Bacillus subtilis cell cultures, Luria Broth (LB) and Minimal Salts glutamate glycerol (MSgg), respectively. We first observe cell clusters in cell suspensions with various concentrations of PEG, and measure cluster size in both static and dynamic fluid environments. We find that cells gather together into big clusters and most of the cells are arranged longitudinally; and the large cell clusters are divided into smaller aggregates under fluid shear. We then use a rheometer to measure the viscoelastic properties of various cell cultures, to represent the degree of aggregation of the bacterial suspensions. We find the storage modulus, the loss modulus and the viscosity of bacterial suspensions not only depend on the cell aggregation but also depend on the directionality of cellular motion.

A multiphase theory for spreading microbial swarms and films

Bacterial swarming and biofilm formation are collective multicellular phenomena through which diverse microbial species colonize and spread over water-permeable tissue. During both modes of surface translocation, fluid uptake and transport play a key role in shaping the overall morphology and spreading dynamics. Here we develop a generalized two-phase thin-film model that couples bacterial growth, extracellular matrix swelling, fluid flow, and nutrient transport to describe the expansion of both highly motile bacterial swarms, and sessile bacterial biofilms. We show that swarm expansion corresponds to steady-state solutions in a nutrient-rich, capillarity dominated regime. In contrast, biofilm colony growth is described by transient solutions associated with a nutrient-limited, extracellular polymer stress driven limit. We apply our unified framework to explain a range of recent experimental observations of steady and unsteady expansion of microbial swarms and biofilms. Our results demonstrate how the physics of flow and transport in slender geometries serve to constrain biological organization in microbial communities.

Bacteria can grow and thrive in many different environments. Although we usually think of bacteria as single-celled organisms, they are not always solitary; they can also form groups containing large numbers of individuals. These aggregates work together as one super-colony, allowing the bacteria to feed and protect themselves more efficiently than they could as isolated cells.

These colonies move and grow in characteristic patterns as they respond to their environment. They can form swarms, like insects, or biofilms, which are thin, flat structures containing both cells and a film-like substance that the cells secrete. Availability of food and water influences the way colonies spread; however, since movement and growth are accompanied by mechanical forces, physical constraints are also important. These include the ability of the bacteria to change the water balance and their local mechanical environment, and the forces they create as they grow and move.

Previous research has used a variety of experimental and theoretical approaches to explain the dynamics of bacterial swarms and biofilms as separate phenomena. However, while they do differ biologically, they also share many physical characteristics.

Srinivasan et al. wanted to exploit these similarities, and use them to predict the growth and shape of biofilms and bacterial swarms under different conditions. To do this, a unified mathematical model for the growth of both swarms and biofilms was created. The model accounted for various factors, such as the transport of nutrients into the colony, the movement of water between the colony and the surface on which it grew, and mechanical changes in the environment (e.g. swelling/softening). The theoretical results were then compared with results from experimental measurements of different bacterial aggregates grown on a soft, hydrated gel. For both swarms and biofilms, the model correctly predicted how fast the colony expanded overall, as well as the shape and location of actively growing regions.

Biofilms and other bacterial aggregates can cause diseases and increase inflammation in tissues, and also hinder industrial processes by damage to submerged surfaces, such as ships and waterpipes. The results described here may open up new approaches to restrict the spreading of bacterial aggregates by focusing on their physical constraints.