Confinement-induced self-organization in growing bacterial colonies

Diverse morphology and motility induced emergent order in bacterial collectives

Microbial communities exhibit complex behaviors driven by species interactions and individual characteristics. In this study, we delve into the dynamics of a mixed bacterial population comprising two distinct species with different morphology and motility aspects. Employing agent-based modeling and computer simulations, we analyze the impacts of size ratios and packing fractions on dispersal patterns, aggregate formation, clustering, and spatial ordering. Notably, we find that motility and anisotropy of elongated bacteria significantly influence the distribution and spatial organization of nonmotile spherical species. Passive spherical cells display a superdiffusive behavior, particularly at larger size ratios in the ballistic regime. As the size ratio increases, clustering of passive cells is observed, accompanied by enhanced alignment and closer packing of active cells in the presence of higher passive cell area fractions. In addition, we identify the pivotal role of passive cell area fraction in influencing the response of active cells toward nematicity, with its dependence on size ratio. These findings shed light on the significance of morphology and motility in shaping the collective behavior of microbial communities, providing valuable insights into complex microbial behaviors with implications for ecology, biotechnology, and bioengineering.

Optogenetic Control of Bacterial Cell-Cell Adhesion Dynamics: Unraveling the Influence on Biofilm Architecture and Functionality

The transition of bacteria from an individualistic to a biofilm lifestyle profoundly alters their biology. During biofilm development, the bacterial cell-cell adhesions are a major determinant of initial microcolonies, which serve as kernels for the subsequent microscopic and mesoscopic structure of the biofilm, and determine the resulting functionality. In this study, the significance of bacterial cell-cell adhesion dynamics on bacterial aggregation and biofilm maturation is elucidated. Using photoswitchable adhesins between bacteria, modifying the dynamics of bacterial cell-cell adhesions with periodic dark-light cycles is systematic. Dynamic cell-cell adhesions with liquid-like behavior improve bacterial aggregation and produce more compact microcolonies than static adhesions with solid-like behavior in both experiments and individual-based simulations. Consequently, dynamic cell-cell adhesions give rise to earlier quorum sensing activation, better intermixing of different bacterial populations, improved biofilm maturation, changes in the growth of cocultures, and higher yields in fermentation. The here presented approach of tuning bacterial cell-cell adhesion dynamics opens the door for regulating the structure and function of biofilms and cocultures with potential biotechnological applications.

Topological defects in multi-layered swarming bacteria

Topological defects, which are singular points in a director field, play a major role in shaping active systems. Here, we experimentally study topological defects and the flow patterns around them, that are formed during the highly rapid dynamics of swarming bacteria. The results are compared to the predictions of two-dimensional active nematics. We show that, even though some of the assumptions underlying the theory do not hold, the swarm dynamics is in agreement with two-dimensional nematic theory. In particular, we look into the multi-layered structure of the swarm, which is an important feature of real, natural colonies, and find a strong coupling between layers. Our results suggest that the defect-charge density is hyperuniform, i.e., that long range density-fluctuations are suppressed.

Stress and alignment response to curved obstacles in growing bacterial monolayers

Monolayers of growing bacteria, confined within channel geometries, exhibit self-organization into a highly aligned laminar state along the axis of the channel. Although this phenomenon has been observed in experiments and simulations under various boundary conditions, the underlying physical mechanism driving this alignment remains unclear. In this study, we conduct simulations of growing bacteria in two-dimensional channel geometries perturbed by fixed obstacles, either circular or arc shaped, placed at the channel's center. Our findings reveal that even sizable obstacles cause only short-ranged disruptions to the baseline laminar state. These disruptions arise from a competition between local planar anchoring and bulk laminar alignment. At smaller obstacle sizes, bulk alignment fully dominates, while at larger sizes planar anchoring induces increasing local disruptions. Furthermore, our analysis indicates that the resulting configurations of the bacterial system display a striking resemblance to the arrangement of hard-rod smectic liquid crystals around circular obstacles. This suggests that modeling hard-rod bacterial monolayers as smectic, rather than nematic, liquid crystals may yield successful outcomes. The insights gained from our study contribute to the expanding body of research on bacterial growth in channels. Our work provides perspectives on the stability of the laminar state and extends our understanding to encompass more intricate confinement schemes.

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.

Interplay between environmental yielding and dynamic forcing modulates bacterial growth

Many bacterial habitats-ranging from gels and tissues in the body to cell-secreted exopolysaccharides in biofilms-are rheologically complex, undergo dynamic external forcing, and have unevenly distributed nutrients. How do these features jointly influence how the resident cells grow and proliferate? Here, we address this question by studying the growth of Escherichia coli dispersed in granular hydrogel matrices with defined and highly tunable structural and rheological properties, under different amounts of external forcing imposed by mechanical shaking, and in both aerobic and anaerobic conditions. Our experiments establish a general principle: that the balance between the yield stress of the environment that the cells inhabit, σy, and the external stress imposed on the environment, σ, modulates bacterial growth by altering transport of essential nutrients to the cells. In particular, when σy<σ, the environment is easily fluidized and mixed over large scales, providing nutrients to the cells and sustaining complete cellular growth. By contrast, when σy>σ, the elasticity of the environment suppresses large-scale fluid mixing, limiting nutrient availability and arresting cellular growth. Our work thus reveals a new mechanism, beyond effects that change cellular behavior via local forcing, by which the rheology of the environment may modulate microbial physiology in diverse natural and industrial settings.

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.

Biofilms as self-shaping growing nematics

Active nematics are the nonequilibrium analogue of passive liquid crystals. They consist of anisotropic units that consume free energy to drive emergent behaviour. Like liquid crystal molecules in displays, ordering and dynamics in active nematics are sensitive to boundary conditions. However, unlike passive liquid crystals, active nematics have the potential to regulate their boundaries through self-generated stresses. Here, we show how a three-dimensional, living nematic can actively shape itself and its boundary to regulate its internal architecture through growth-induced stresses, using bacterial biofilms confined by a hydrogel as a model system. We show that biofilms exhibit a sharp transition in shape from domes to lenses upon changing environmental stiffness or cell-substrate friction, which is explained by a theoretical model that considers the competition between confinement and interfacial forces. The growth mode defines the progression of the boundary, which in turn determines the trajectories and spatial distribution of cell lineages. We further demonstrate that the evolving boundary and corresponding stress anisotropy define the orientational ordering of cells and the emergence of topological defects in the biofilm interior. Our findings may provide strategies for the development of programmed microbial consortia with emergent material properties.

Instabilities in a growing system of active particles: scalar and vectorial systems

The physics of micron-scale biological colonies usually benefits from different out-of-equilibrium sources. In bacterial colonies and cellular tissues, the growth process is among the important active sources that determine the dynamics. In this article, we study the generic dynamical instabilities associated with the growth phenomena that may arise in both scalar and vectorial systems. In vectorial systems, where the rotational degrees of particles play a role, a phenomenological growth-mediated torque can affect the rotational dynamics of individual particles. We show that such a growth-mediated torque can result in active traveling waves in the bulk of a growing system. In addition to the bulk properties, we analyze the instabilities in the shape of growing interfaces in both scalar and vectorial systems.

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.

Robust Estimation of Position-Dependent Anisotropic Diffusivity Tensors from Stochastic Trajectories

Materials under confinement can possess properties that deviate considerably from their bulk counterparts. Indeed, confinement makes all physical properties position-dependent and possibly anisotropic, and characterizing such spatial variations and directionality has been an intense area of focus in experimental and computational studies of confined matter. While this task is fairly straightforward for simple mechanical observables, it is far more daunting for transport properties such as diffusivity that can only be estimated from autocorrelations of mechanical observables. For instance, there are well established methods for estimating diffusivity from experimentally observed or computationally generated trajectories in bulk systems. No rigorous generalizations of such methods, however, exist for confined systems. In this work, we present two filtered covariance estimators for computing anisotropic and position-dependent diffusivity tensors and validate them by applying them to stochastic trajectories generated according to known diffusivity profiles. These estimators can accurately capture spatial variations that span over several orders of magnitude and that assume different functional forms. Our kernel-based approach is also very robust to implementation details such as the localization function and time discretization and performs significantly better than estimators that are solely based on local covariance. Moreover, the kernel function does not have to be localized and can instead belong to a dictionary of orthogonal functions. Therefore, the proposed estimator can be readily used to obtain functional estimates of diffusivity rather than a tabulated collection of pointwise estimates. Nonetheless, the susceptibility of the proposed estimators to time discretization is higher at the immediate vicinity of hard boundaries. We demonstrate this heightened susceptibility to be common among all covariance-based estimators.

Microdomains and stress distributions in bacterial monolayers on curved interfaces

Monolayers of growing non-motile rod-shaped bacteria act as active nematic materials composed of hard particles rather than the flexible components of other commonly studied active nematics. The organization of these granular monolayers has been studied on flat surfaces but not on curved surfaces, which are known to change the behavior of other active nematics. We use molecular dynamics simulations to track alignment and stress in growing monolayers fixed to curved surfaces, and investigate how these vary with changing surface curvature and cell aspect ratio. We find that the length scale of alignment (measured by average microdomain size) increases with cell aspect ratio and decreases with curvature. Additionally, we find that alignment controls the distribution of extensile stresses in the monolayer by concentrating stress in negative-order regions. These results connect active nematic physics to bacterial monolayers and can be applied to model bacteria growing on droplets, such as oil-degrading marine bacteria.

Hierarchical encapsulation of bacteria in functional hydrogel beads for inter- and intra- species communication

To sequester prokaryotic cells in a biofilm-like niche, the creation of a pertinent and reliable microenvironment that reflects the heterogeneous nature of biological systems is vital for sustenance. Design of a microenvironment that is conducive for growth and survival of organisms, should account for factors such as mass transport, porosity, stability, elasticity, size, functionality, and biochemical characteristics of the organisms in the confined architecture. In this work we present an artificial long-term confinement model fabricated by natural alginate hydrogels that are structurally stable and can host organisms for over 10 days in physiologically relevant conditions. A unique feature of the confinement platform is the development of stratified habitats wherein bacterial cells can be entrapped in the core as well as in the shell layers, wherein the thickness and the number of shell layers are tunable at fabrication. We show that the hydrogel microenvironment in the beads can host complex subpopulations of organisms similar to that in a biofilm. Dynamic interaction of bacterial colonies encapsulated in different beads or within the core and stratified layers of single beads was demonstrated to show intra- species communication. Inter- species communication between probiotic bacteria and human colorectal carcinoma cells was also demonstrated to highlight a possible bidirectional communication between the organisms in the beads and the environment. STATEMENT OF SIGNIFICANCE: Bacteria confinement in a natural soft hydrogel structure has always been a challenge due to the collapse of hydrogel architectures. Alternative methods have been attempted to encapsulate microorganisms by employing various processes to avoid/minimize rupturing of hydrogel structures. However, most of the past approaches have been unfavorable in balancing cell proliferation and functionality upon confinement. Our study addresses the fundamental gap in knowledge necessary to create favorable and complex 3D biofilm mimics utilizing natural hydrogel for microbial colonization for long-term studies. Our approach represents a cornerstone in the development of 3D functional architectures not only to advance studies in microbial communication, host-microbe interaction but also to address basic and fundamental questions in biology.

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.

Stress anisotropy in confined populations of growing rods

A central feature of living matter is its ability to grow and multiply. The mechanical activity associated with growth produces both macroscopic flows shaped by confinement, and striking self-organization phenomena, such as orientational order and alignment, which are particularly prominent in populations of rod-shaped bacteria due to their nematic properties. However, how active stresses, passive mechanical interactions and flow-induced effects interact to give rise to the observed global alignment patterns remains elusive. Here, we study in silico colonies of growing rod-shaped particles of different aspect ratios confined in channel-like geometries. A spatially resolved analysis of the stress tensor reveals a strong relationship between near-perfect alignment and an inversion of stress anisotropy for particles with large length-to-width ratios. We show that, in quantitative agreement with an asymptotic theory, strong alignment can lead to a decoupling of active and passive stresses parallel and perpendicular to the direction of growth, respectively. We demonstrate the robustness of these effects in a geometry that provides less restrictive confinement and introduces natural perturbations in alignment. Our results illustrate the complexity arising from the inherent coupling between nematic order and active stresses in growing active matter, which is modulated by geometric and configurational constraints due to confinement.

Large-scale orientational order in bacterial colonies during inward growth

During colony growth, complex interactions regulate the bacterial orientation, leading to the formation of large-scale ordered structures, including topological defects, microdomains, and branches. These structures may benefit bacterial strains, providing invasive advantages during colonization. Active matter dynamics of growing colonies drives the emergence of these ordered structures. However, additional biomechanical factors also play a significant role during this process. Here, we show that the velocity profile of growing colonies creates strong radial orientation during inward growth when crowded populations invade a closed area. During this process, growth geometry sets virtual confinement and dictates the velocity profile. Herein, flow-induced alignment and torque balance on the rod-shaped bacteria result in a new stable orientational equilibrium in the radial direction. Our analysis revealed that the dynamics of these radially oriented structures, also known as aster defects, depend on bacterial length and can promote the survival of the longest bacteria around localized nutritional hotspots. The present results indicate a new mechanism underlying structural order and provide mechanistic insights into the dynamics of bacterial growth on complex surfaces.

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.

Topology of Orientational Defects in Confined Smectic Liquid Crystals

We propose a general formalism to characterize orientational frustration of smectic liquid crystals in confinement by interpreting the emerging networks of grain boundaries as objects with a topological charge. In a formal idealization, this charge is distributed in pointlike units of quarter-integer magnitude, which we identify with tetratic disclinations located at the end points and nodes. This coexisting nematic and tetratic order is analyzed with the help of extensive Monte Carlo simulations for a broad range of two-dimensional confining geometries as well as colloidal experiments, showing how the observed defect networks can be universally reconstructed from simple building blocks. We further find that the curvature of the confining wall determines the anchoring behavior of grain boundaries, such that the number of nodes in the emerging networks and the location of their end points can be tuned by changing the number and smoothness of corners, respectively.

Emergent spatiotemporal population dynamics with cell-length control of synthetic microbial consortia

The increased complexity of synthetic microbial biocircuits highlights the need for distributed cell functionality due to concomitant increases in metabolic and regulatory burdens imposed on single-strain topologies. Distributed systems, however, introduce additional challenges since consortium composition and spatiotemporal dynamics of constituent strains must be robustly controlled to achieve desired circuit behaviors. Here, we address these challenges with a modeling-based investigation of emergent spatiotemporal population dynamics using cell-length control in monolayer, two-strain bacterial consortia. We demonstrate that with dynamic control of a strain's division length, nematic cell alignment in close-packed monolayers can be destabilized. We find that this destabilization confers an emergent, competitive advantage to smaller-length strains-but by mechanisms that differ depending on the spatial patterns of the population. We used complementary modeling approaches to elucidate underlying mechanisms: an agent-based model to simulate detailed mechanical and signaling interactions between the competing strains, and a reductive, stochastic lattice model to represent cell-cell interactions with a single rotational parameter. Our modeling suggests that spatial strain-fraction oscillations can be generated when cell-length control is coupled to quorum-sensing signaling in negative feedback topologies. Our research employs novel methods of population control and points the way to programming strain fraction dynamics in consortial synthetic biology.

Morphogenesis and cell ordering in confined bacterial biofilms

Biofilms are aggregates of bacterial cells surrounded by an extracellular matrix. Much progress has been made in studying biofilm growth on solid substrates; however, little is known about the biophysical mechanisms underlying biofilm development in three-dimensional confined environments in which the biofilm-dwelling cells must push against and even damage the surrounding environment to proliferate. Here, combining single-cell imaging, mutagenesis, and rheological measurement, we reveal the key morphogenesis steps of Vibrio cholerae biofilms embedded in hydrogels as they grow by four orders of magnitude from their initial size. We show that the morphodynamics and cell ordering in embedded biofilms are fundamentally different from those of biofilms on flat surfaces. Treating embedded biofilms as inclusions growing in an elastic medium, we quantitatively show that the stiffness contrast between the biofilm and its environment determines biofilm morphology and internal architecture, selecting between spherical biofilms with no cell ordering and oblate ellipsoidal biofilms with high cell ordering. When embedded in stiff gels, cells self-organize into a bipolar structure that resembles the molecular ordering in nematic liquid crystal droplets. In vitro biomechanical analysis shows that cell ordering arises from stress transmission across the biofilm-environment interface, mediated by specific matrix components. Our imaging technique and theoretical approach are generalizable to other biofilm-forming species and potentially to biofilms embedded in mucus or host tissues as during infection. Our results open an avenue to understand how confined cell communities grow by means of a compromise between their inherent developmental program and the mechanical constraints imposed by the environment.