Branching instability in expanding bacterial colonies

Quantifying Geobacter sulfurreducens growth: A mathematical model based on acetate concentration as an oxidizing substrate

We developed a mathematical model to simulate dynamics associated with the proliferation of Geobacter and ultimately optimize cellular operation by analyzing the interaction of its components. The model comprises two segments: an initial part comprising a logistic form and a subsequent segment that incorporates acetate oxidation as a saturation term for the microbial nutrient medium. Given that four parameters can be obtained by minimizing the square root of the mean square error between experimental Geobacter growth and the mathematical model, the model underscores the importance of incorporating nonlinear terms. The determined parameter values closely align with experimental data, providing insights into the mechanisms that govern Geobacter proliferation. Furthermore, the model has been transformed into a scaleless equation with only two parameters to simplify the exploration of qualitative properties. This allowed us to conduct stability analysis of the fixed point and construct a co-dimension two bifurcation diagram.

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.

Parallel evolution of alternate morphotypes of Chryseobacterium gleum during experimental evolution with Caenorhabditis elegans

Microbial evolution within polymicrobial communities is a complex process. Here, we report within-species diversification within multispecies microbial communities during experimental evolution with the nematode Caenorhabditis elegans. We describe morphological diversity in the target species Chryseobacterium gleum, which developed a novel colony morphotype in a small number of replicate communities. Alternate morphotypes coexisted with original morphotypes in communities, as well as in single-species experiments using evolved isolates. We found that the original and alternate morphotypes differed in motility and in spatial expansion in the presence of C. elegans. This study provides insight into the emergence and maintenance of intraspecies diversity in the context of microbial communities.

Analysis of biofilm expansion rate of Bacillus subtilis (MTC871) on agar substrates with different stiffness

The surface morphology of mature biofilms is heterogeneous and can be divided into concentric rings wrinkles (I), labyrinthine networks wrinkles (II), radial ridges wrinkles (III), and branches wrinkles (IV), according to surface wrinkle structure and distribution characteristics. Due to the wrinkle structures, channels are formed between the biofilm and substrate and transport nutrients, water, metabolic products, etc. We find that expansion rate variations of biofilms growing on substrates with high and low agar concentrations (1.5, 2.0, 2.5 wt.%) are not in the same phase. In the first 3 days' growth, the interaction stress between biofilm and each agar substrate increases, which makes the biofilm expansion rate decreases before wrinkle pattern IV (branches) comes up. After 3 days, in the later growth stage after wrinkle pattern IV appears, the biofilm has larger expansion rate growing on 2.0 wt.% agar concentration, which has the larger wrinkle distance in wrinkle pattern IV reducing energy consumption. Our study shows that the stiff substrate does not always inhibit the biofilm expansion, although it does in the earlier stage; after that, mature biofilms acquire larger expansion rate by adjusting the growth mode through the wrinkle evolution even in nutrient extremely depletion.

Swarming of P. aeruginosa: Through the lens of biophysics

Swarming is a collective flagella-dependent movement of bacteria across a surface that is observed across many species of bacteria. Due to the prevalence and diversity of this motility modality, multiple models of swarming have been proposed, but a consensus on a general mechanism for swarming is still lacking. Here, we focus on swarming by Pseudomonas aeruginosa due to the abundance of experimental data and multiple models for this species, including interpretations that are rooted in biology and biophysics. In this review, we address three outstanding questions about P. aeruginosa swarming: what drives the outward expansion of a swarm, what causes the formation of dendritic patterns (tendrils), and what are the roles of flagella? We review models that propose biologically active mechanisms including surfactant sensing as well as fluid mechanics-based models that consider swarms as thin liquid films. Finally, we reconcile recent observations of P. aeruginosa swarms with early definitions of swarming. This analysis suggests that mechanisms associated with sliding motility have a critical role in P. aeruginosa swarm formation.

BMX: Biological modelling and interface exchange

High performance computing has a great potential to provide a range of significant benefits for investigating biological systems. These systems often present large modelling problems with many coupled subsystems, such as when studying colonies of bacteria cells. The aim to understand cell colonies has generated substantial interest as they can have strong economic and societal impacts through their roles in in industrial bioreactors and complex community structures, called biofilms, found in clinical settings. Investigating these communities through realistic models can rapidly exceed the capabilities of current serial software. Here, we introduce BMX, a software system developed for the high performance modelling of large cell communities by utilising GPU acceleration. BMX builds upon the AMRex adaptive mesh refinement package to efficiently model cell colony formation under realistic laboratory conditions. Using simple test scenarios with varying nutrient availability, we show that BMX is capable of correctly reproducing observed behavior of bacterial colonies on realistic time scales demonstrating a potential application of high performance computing to colony modelling. The open source software is available from the zenodo repository https://doi.org/10.5281/zenodo.8084270 under the BSD-2-Clause licence.

Geometrical control of interface patterning underlies active matter invasion

Interaction between active materials and the boundaries of geometrical confinement is key to many emergent phenomena in active systems. For living active matter consisting of animal cells or motile bacteria, the confinement boundary is often a deformable interface, and it has been unclear how activity-induced interface dynamics might lead to morphogenesis and pattern formation. Here, we studied the evolution of bacterial active matter confined by a deformable boundary. We found that an ordered morphological pattern emerged at the interface characterized by periodically spaced interfacial protrusions; behind the interfacial protrusions, bacterial swimmers self-organized into multicellular clusters displaying +1/2 nematic defects. Subsequently, a hierarchical sequence of transitions from interfacial protrusions to creeping branches allowed the bacterial active drop to rapidly invade surrounding space with a striking self-similar branch pattern. We found that this interface patterning is geometrically controlled by the local curvature of the interface, a phenomenon we denote as collective curvature sensing. Using a continuum active model, we revealed that the collective curvature sensing arises from enhanced active stresses near high-curvature regions, with the active length scale setting the characteristic distance between the interfacial protrusions. Our findings reveal a protrusion-to-branch transition as a unique mode of active matter invasion and suggest a strategy to engineer pattern formation of active materials.

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.

Long-range alteration of the physical environment mediates cooperation between Pseudomonas aeruginosa swarming colonies

Pseudomonas aeruginosa makes and secretes massive amounts of rhamnolipid surfactants that enable swarming motility over biogel surfaces. But how these rhamnolipids interact with biogels to assist swarming remains unclear. Here, I use a combination of optical techniques across scales and genetically engineered strains to demonstrate that rhamnolipids can induce agar gel swelling over distances >10,000× the body size of an individual cell. The swelling front is on the micrometric scale and is easily visible using shadowgraphy. Rhamnolipid transport is not restricted to the surface of the gel but occurs through the whole thickness of the plate and, consequently, the spreading dynamics depend on the local thickness. Surprisingly, rhamnolipids can cross the whole gel and induce swelling on the opposite side of a two-face Petri dish. The swelling front delimits an area where the mechanical properties of the surface properties are modified: water wets the surface more easily, which increases the motility of individual bacteria and enables collective motility. A genetically engineered mutant unable to secrete rhamnolipids (ΔrhlA), and therefore unable to swarm, is rescued from afar with rhamnolipids produced by a remote colony. These results exemplify the remarkable capacity of bacteria to change the physical environment around them and its ecological consequences.

'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.

Individual-Based Modeling of Spatial Dynamics of Chemotactic Microbial Populations

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

Morphological instability and roughening of growing 3D bacterial colonies

The morphogenesis of two-dimensional bacterial colonies has been well studied. However, little is known about the colony morphologies of bacteria growing in three dimensions, despite the prevalence of three-dimensional environments (e.g., soil, inside hosts) as natural bacterial habitats. Using experiments on bacteria in granular hydrogel matrices, we find that dense multicellular colonies growing in three dimensions undergo a common morphological instability and roughen, adopting a characteristic broccoli-like morphology when they exceed a critical size. Analysis of a continuum “active fluid” model of the expanding colony reveals that this behavior originates from an interplay of competition for nutrients with growth-driven colony expansion, both of which vary spatially. These results shed light on the fundamental biophysical principles underlying growth in three dimensions.

How do growing bacterial colonies get their shapes? While colony morphogenesis is well studied in two dimensions, many bacteria grow as large colonies in three-dimensional (3D) environments, such as gels and tissues in the body or subsurface soils and sediments. Here, we describe the morphodynamics of large colonies of bacteria growing in three dimensions. Using experiments in transparent 3D granular hydrogel matrices, we show that dense colonies of four different species of bacteria generically become morphologically unstable and roughen as they consume nutrients and grow beyond a critical size—eventually adopting a characteristic branched, broccoli-like morphology independent of variations in the cell type and environmental conditions. This behavior reflects a key difference between two-dimensional (2D) and 3D colonies; while a 2D colony may access the nutrients needed for growth from the third dimension, a 3D colony inevitably becomes nutrient limited in its interior, driving a transition to unstable growth at its surface. We elucidate the onset of the instability using linear stability analysis and numerical simulations of a continuum model that treats the colony as an “active fluid” whose dynamics are driven by nutrient-dependent cellular growth. We find that when all dimensions of the colony substantially exceed the nutrient penetration length, nutrient-limited growth drives a 3D morphological instability that recapitulates essential features of the experimental observations. Our work thus provides a framework to predict and control the organization of growing colonies—as well as other forms of growing active matter, such as tumors and engineered living materials—in 3D environments.

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.

Dependency of active pressure and equation of state on stiffness of wall

Autonomous motion and motility are hallmarks of active matter. Active agents, such as biological cells and synthetic colloidal particles, consume internal energy or extract energy from the environment to generate self-propulsion and locomotion. These systems are persistently out of equilibrium due to continuous energy consumption. It is known that pressure is not always a state function for generic active matter. Torque interaction between active constituents and confinement renders the pressure of the system a boundary-dependent property. The mechanical pressure of anisotropic active particles depends on their microscopic interactions with a solid wall. Using self-propelled dumbbells confined by solid walls as a model system, we perform numerical simulations to explore how variations in the wall stiffness influence the mechanical pressure of dry active matter. In contrast to previous findings, we find that mechanical pressure can be independent of the interaction of anisotropic active particles with walls, even in the presence of intrinsic torque interaction. Particularly, the dependency of pressure on the wall stiffness vanishes when the stiffness is above a critical level. In such a limit, the dynamics of dumbbells near the walls are randomized due to the large torque experienced by the dumbbells, leading to the recovery of pressure as a state variable of density.

Programming cell growth into different cluster shapes using diffusible signals

Advances in genetic engineering technologies have allowed the construction of artificial genetic circuits, which have been used to generate spatial patterns of differential gene expression. However, the question of how cells can be programmed, and how complex the rules need to be, to achieve a desired tissue morphology has received less attention. Here, we address these questions by developing a mathematical model to study how cells can collectively grow into clusters with different structural morphologies by secreting diffusible signals that can influence cellular growth rates. We formulate how growth regulators can be used to control the formation of cellular protrusions and how the range of achievable structures scales with the number of distinct signals. We show that a single growth inhibitor is insufficient for the formation of multiple protrusions but may be achieved with multiple growth inhibitors, and that other types of signals can regulate the shape of protrusion tips. These examples illustrate how our approach could potentially be used to guide the design of regulatory circuits for achieving a desired target structure.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

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.

The spatial organization of microbial communities during range expansion

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

Collective colony growth is optimized by branching pattern formation in Pseudomonas aeruginosa

Branching pattern formation is common in many microbes. Extensive studies have focused on addressing how such patterns emerge from local cell-cell and cell-environment interactions. However, little is known about whether and to what extent these patterns play a physiological role. Here, we consider the colonization of bacteria as an optimization problem to find the colony patterns that maximize colony growth efficiency under different environmental conditions. We demonstrate that Pseudomonas aeruginosa colonies develop branching patterns with characteristics comparable to the prediction of modeling; for example, colonies form thin branches in a nutrient-poor environment. Hence, the formation of branching patterns represents an optimal strategy for the growth of Pseudomonas aeruginosa colonies. The quantitative relationship between colony patterns and growth conditions enables us to develop a coarse-grained model to predict diverse colony patterns under more complex conditions, which we validated experimentally. Our results offer new insights into branching pattern formation as a problem-solving social behavior in microbes and enable fast and accurate predictions of complex spatial patterns in branching colonies.

An expanding bacterial colony forms a depletion zone with growing droplets

Many species of bacteria have developed effective means to spread on solid surfaces. This study focuses on the expansion of Pseudomonas aeruginosa on an agar gel surface under conditions of minimal evaporation. We report the occurrence and spread of a depletion zone within an expanded colony, where the bacteria laden film becomes thinner. The depletion zone is colocalized with a higher concentration of rhamnolipids, the biosurfactants that are produced by the bacteria and accumulate in the older region of the colony. With continued growth in population, dense bacterial droplets occur and coalesce in the depletion zone, displaying remarkable fluid dynamic behavior. Whereas expansion of a central depletion zone requires activities of live bacteria, new zones can be seeded elsewhere by adding rhamnolipids. These depletion zones due to the added surfactants expand quickly, even on plates covered by bacteria that have been killed by ultraviolet light. We explain the observed properties based on considerations of bacterial growth and secretion, osmotic swelling, fluid volume expansion, interfacial fluid dynamics involving Marangoni and capillary flows, and cell-cell cohesion.

Geodesic Loewner paths with varying boundary conditions

Equations of the Loewner class subject to non-constant boundary conditions along the real axis are formulated and solved giving the geodesic paths of slits growing in the upper half complex plane. The problem is motivated by Laplacian growth in which the slits represent thin fingers growing in a diffusion field. A single finger follows a curved path determined by the forcing function appearing in Loewner's equation. This function is found by solving an ordinary differential equation whose terms depend on curvature properties of the streamlines of the diffusive field in the conformally mapped 'mathematical' plane. The effect of boundary conditions specifying either piecewise constant values of the field variable along the real axis, or a dipole placed on the real axis, reveal a range of behaviours for the growing slit. These include regions along the real axis from which no slit growth is possible, regions where paths grow to infinity, or regions where paths curve back toward the real axis terminating in finite time. Symmetric pairs of paths subject to the piecewise constant boundary condition along the real axis are also computed, demonstrating that paths which grow to infinity evolve asymptotically toward an angle of bifurcation of π/5.

THE EVOLVING WRINKLE PATTERN OF THE BACILLUS SUBTILIS BIOFILM PROVIDING MORE LIVING SPACE FOR CELLS

The biofilm wrinkle evolution is the growth mechanism by which bacteria regulate their physiological state in response to the environmental change. We use the parameter of surface complexity to describe different wrinkle patterns. The surface complexity is defined that the biofilm surface area contact with the air is divided by the projected area of the biofilm. We find that the biofilm surface complexity variation is positively proportional to the number of spores. Although each wrinkle pattern has various wrinkle thickness and width, surface complexities of some patterns are almost same, which guarantees cells have enough living space. Through the observation of the growth of the damaged biofilm, we further find that the biofilm expansion along the circumferential direction is faster than that along radial direction, which means that the internal stress along the circumferential direction contributes the wrinkle formation. Our work provides a new perspective to study biofilm morphologies, and relates the morphology evolution with phenotypes in the Bacillus subtilis biofilm.

Bacterial surface motility is modulated by colony-scale flow and granular jamming

Microbes routinely face the challenge of acquiring territory and resources on wet surfaces. Cells move in large groups inside thin, surface-bound water layers, often achieving speeds of 30 µm s^-1 within this environment, where viscous forces dominate over inertial forces (low Reynolds number). The canonical Gram-positive bacterium Bacillus subtilis is a model organism for the study of collective migration over surfaces with groups exhibiting motility on length-scales three orders of magnitude larger than themselves within a few doubling times. Genetic and chemical studies clearly show that the secretion of endogenous surfactants and availability of free surface water are required for this fast group motility. Here, we show that: (i) water availability is a sensitive control parameter modulating an abiotic jamming-like transition that determines whether the group remains fluidized and therefore collectively motile, (ii) groups self-organize into discrete layers as they travel, (iii) group motility does not require proliferation, rather groups are pulled from the front, and (iv) flow within expanding groups is capable of moving material from the parent colony into the expanding tip of a cellular dendrite with implications for expansion into regions of varying nutrient content. Together, these findings illuminate the physical structure of surface-motile groups and demonstrate that physical properties, like cellular packing fraction and flow, regulate motion from the scale of individual cells up to length scales of centimetres.

Application of the Schwarz-Christoffel map to the Laplacian growth of needles and fingers

A numerical procedure based on the Schwarz-Christoffel map suitable for the study of the Laplacian growth of thin two-dimensional protrusions is presented. The protrusions take the form of either straight needles or curved fingers satisfying Loewner's equation, and are represented by slits in the complex plane. Particular use is made of Driscoll's numerical procedure, the SC Toolbox, for computing the Schwarz-Christoffel map from a half plane to a slit half plane. Since the Schwarz-Christoffel map applies only to polygonal regions, the growth of curved fingers is approximated by an increasing number of short straight line segments. The growth rate is given by a fixed power η of the harmonic measure at the finger or needle tips and so includes the possibility of "screening" as the needles of fingers interact with themselves and with boundaries. The method is illustrated with examples of multiple needle and finger growth in half-plane and channel geometries. The effect of η on the trajectories of asymmetric bifurcating fingers is also studied.

A thin-film extensional flow model for biofilm expansion by sliding motility

In the presence of glycoproteins, bacterial and yeast biofilms are hypothesized to expand by sliding motility. This involves a sheet of cells spreading as a unit, facilitated by cell proliferation and weak adhesion to the substratum. In this paper, we derive an extensional flow model for biofilm expansion by sliding motility to test this hypothesis. We model the biofilm as a two-phase (living cells and an extracellular matrix) viscous fluid mixture, and model nutrient depletion and uptake from the substratum. Applying the thin-film approximation simplifies the model, and reduces it to one-dimensional axisymmetric form. Comparison with Saccharomyces cerevisiae mat formation experiments reveals good agreement between experimental expansion speed and numerical solutions to the model withO ( 1 ) parameters estimated from experiments. This confirms that sliding motility is a possible mechanism for yeast biofilm expansion. Having established the biological relevance of the model, we then demonstrate how the model parameters affect expansion speed, enabling us to predict biofilm expansion for different experimental conditions. Finally, we show that our model can explain the ridge formation observed in some biofilms. This is especially true if surface tension is low, as hypothesized for sliding motility.

Front Instabilities Can Reverse Desertification

Degradation processes in living systems often take place gradually by front propagation. An important context of such processes is loss of biological productivity in drylands or desertification. Using a dryland-vegetation model, we analyze the stability and dynamics of desertification fronts, identify linear and nonlinear front instabilities, and highlight the significance of these instabilities in inducing self-recovery. The results are based on the derivation and analysis of a universal amplitude equation for pattern-forming living systems for which nonuniform instabilities cannot emerge from the nonviable (zero) state. The results may therefore be applicable to other contexts of animate matter where degradation processes occur by front propagation.

Nonuniversality of front fluctuations for compact colonies of nonmotile bacteria

The front of a compact bacterial colony growing on a Petri dish is a paradigmatic instance of non-equilibrium fluctuations in the celebrated Eden, or Kardar-Parisi-Zhang (KPZ), universality class. While in many experiments the scaling exponents crucially differ from the expected KPZ values, the source of this disagreement has remained poorly understood. We have performed growth experiments with B. subtilis 168 and E. coli ATCC 25922 under conditions leading to compact colonies in the classically alleged Eden regime, where individual motility is suppressed. Non-KPZ scaling is indeed observed for all accessible times, KPZ asymptotics being ruled out for our experiments due to the monotonic increase of front branching with time. Simulations of an effective model suggest the occurrence of transient nonuniversal scaling due to diffusive morphological instabilities, agreeing with expectations from detailed models of the relevant biological reaction-diffusion processes.

Quantifying the dominant growth mechanisms of dimorphic yeast using a lattice-based model

A mathematical model is presented for the growth of yeast that incorporates both dimorphic behaviour and nutrient diffusion. The budding patterns observed in the standard and pseudohyphal growth modes are represented by a bias in the direction of cell proliferation. A set of spatial indices is developed to quantify the morphology and compare the relative importance of the directional bias to nutrient concentration and diffusivity on colony shape. It is found that there are three different growth modes: uniform growth, diffusion-limited growth (DLG) and an intermediate region in which the bias determines the morphology. The dimorphic transition due to nutrient limitation is investigated by relating the directional bias to the nutrient concentration, and this is shown to replicate the behaviour observed in vivo Comparisons are made with experimental data, from which it is found that the model captures many of the observed features. Both DLG and pseudohyphal growth are found to be capable of generating observed experimental morphologies.

Diffusion-Limited Growth of Microbial Colonies

The emergence of diffusion-limited growth (DLG) within a microbial colony on a solid substrate is studied using a combination of mathematical modelling and experiments. Using an agent-based model of the interaction between microbial cells and a diffusing nutrient, it is shown that growth directed towards a nutrient source may be used as an indicator that DLG is influencing the colony morphology. A continuous reaction–diffusion model for microbial growth is employed to identify the parameter regime in which DLG is expected to arise. Comparisons between the model and experimental data are used to argue that the bacterium Bacillus subtilis can undergo DLG, while the yeast Saccharomyces cerevisiae cannot, and thus the non-uniform growth exhibited by this yeast must be caused by the pseudohyphal growth mode rather than limited nutrient availability. Experiments testing directly for DLG features in yeast colonies are used to confirm this hypothesis.

Capillary flow and mechanical buckling in a growing annular bacterial colony

A growing bacterial colony is a dense suspension of an increasing number of cells capable of individual as well as collective motion. After inoculating Pseudomonas aeruginosa over an annular area on an agar plate, we observe the growth and spread of the bacterial population, and model the process by considering the physical effects that account for the features observed. Over a course of 10-12 hours, the majority of bacteria migrate to and accumulate at the edges. We model the capillary flow induced by imbalanced evaporation flux as the cause for the accumulation, much like the well-known coffee stain phenomenon. Simultaneously, periodic buckles or protrusions occur at the inner edge. These buckles indicate that the crowding bacteria produce a jam, transforming the densely packed population at the inner edge to a solid state. The continued bacterial growth produces buckles. Subsequently, a ring of packed bacteria behind the inner edge detach from it and break into pieces, forming bacterial droplets. These droplets slowly coalesce while they continually grow and collectively surf on the agar surface in the region where the colony had previously spread over. Our study shows a clear example of how fluid dynamics and elasto-mechanics together govern the bacterial colony pattern evolution.

Spreading of nonmotile bacteria on a hard agar plate: Comparison between agent-based and stochastic simulations

We study spreading of a nonmotile bacteria colony on a hard agar plate by using agent-based and continuum models. We show that the spreading dynamics depends on the initial nutrient concentration, the motility, and the inherent demographic noise. Population fluctuations are inherent in an agent-based model, whereas for the continuum model we model them by using a stochastic Langevin equation. We show that the intrinsic population fluctuations coupled with nonlinear diffusivity lead to a transition from a diffusion limited aggregation type of morphology to an Eden-like morphology on decreasing the initial nutrient concentration.

Morphodynamics of a growing microbial colony driven by cell death

Bacterial cells can often self-organize into multicellular structures with complex spatiotemporal morphology. In this work, we study the spatiotemporal dynamics of a growing microbial colony in the presence of cell death. We present an individual-based model of nonmotile bacterial cells which grow and proliferate by consuming diffusing nutrients on a semisolid two-dimensional surface. The colony spreads by growth forces and sliding motility of cells and undergoes cell death followed by subsequent disintegration of the dead cells in the medium. We model cell death by considering two possible situations: In one of the cases, cell death occurs in response to the limitation of local nutrients, while the other case corresponds to an active death process, known as apoptotic or programmed cell death. We demonstrate how the colony morphology is influenced by the presence of cell death. Our results show that cell death facilitates transitions from roughly circular to highly branched structures at the periphery of an expanding colony. Interestingly, our results also reveal that for the colonies which are growing in higher initial nutrient concentrations, cell death occurs much earlier compared to the colonies which are growing in lower initial nutrient concentrations. This work provides new insights into the branched patterning of growing bacterial colonies as a consequence of complex interplay among the biochemical and mechanical effects.

Experimental and theoretical studies on the morphogenesis of bacterial biofilms

Biofilm morphogenesis not only reflects the physiological state of bacteria but also serves as a strategy to sustain bacterial survival. In this paper, we take the Bacillus subtilis colony as a model system to explore the morphomechanics of growing biofilms confined in a defined geometry. We find that the growth-induced stresses may drive the occurrence of both surface wrinkling and interface delamination in the biofilm, leading to the formation of a labyrinthine network on its surface. The wrinkles are perpendicular to the boundary of the constraint region. The variation in the surface undulations is attributed to the spatial stress field, which is isotropic in the inner regime but anisotropic in the vicinity of the boundary. Our experiments show that the directional surface wrinkles can confer biofilms with anisotropic wetting properties. This study not only highlights the role of mechanics in sculpturing organisms within the morphogenetic context but also suggests a promising route toward desired surfaces at the interface between synthetic biology and materials sciences.

Influence of Physical Effects on the Swarming Motility of Pseudomonas aeruginosa

Many species of bacteria can spread over a moist surface via a particular form of collective motion known as "surface swarming". This form of motility is typically studied by inoculating bacteria on a gel formed by 0.4-1.5% agar, which contains essential nutrients for their growth and proliferation. Using Pseudomonas aeruginosa and its pili-less mutant, ΔPilA, we investigate physical factors that either facilitate or restrict the swarming motility, measured by the rate of increase in area covered by a spreading bacterial colony, i.e., a swarm. The wild-type colony spreads over the agar surface in highly branched structures. The pili-less mutant fills up the area more fully as it spreads, but it also produces numerous and fragmented branches, or tendrils, at the swarm front. Whereas additional surfactants enhance swarming, increasing the agar percentage, adding extra salt or sugar or incorporating viscous agents in the agar matrix all decrease swarming, supporting the conclusion that swarming motility is restricted by the surface tension at the swarm front and swarm growth is limited by the rate of water supply from within the agar gel. The physical basis elaborated through this study provides a useful framework for understanding the swarming behavior of numerous species of bacteria.

On the morphological stability of multicellular tumour spheroids growing in porous media

Multicellular tumour spheroids (MCTSs) are extensively used as in vitro system models for investigating the avascular growth phase of solid tumours. In this work, we propose a continuous growth model of heterogeneous MCTSs within a porous material, taking into account a diffusing nutrient from the surrounding material directing both the proliferation rate and the mobility of tumour cells. At the time scale of interest, the MCTS behaves as an incompressible viscous fluid expanding inside a porous medium. The cell motion and proliferation rate are modelled using a non-convective chemotactic mass flux, driving the cell expansion in the direction of the external nutrients' source. At the early stages, the growth dynamics is derived by solving the quasi-stationary problem, obtaining an initial exponential growth followed by an almost linear regime, in accordance with experimental observations. We also perform a linear-stability analysis of the quasi-static solution in order to investigate the morphological stability of the radially symmetric growth pattern. We show that mechano-biological cues, as well as geometric effects related to the size of the MCTS subdomains with respect to the diffusion length of the nutrient, can drive a morphological transition to fingered structures, thus triggering the formation of complex shapes that might promote tumour invasiveness. The results also point out the formation of a retrograde flow in the MCTS close to the regions where protrusions form, that could describe the initial dynamics of metastasis detachment from the in vivo tumour mass. In conclusion, the results of the proposed model demonstrate that the integration of mathematical tools in biological research could be crucial for better understanding the tumour's ability to invade its host environment.

Tumour angiogenesis as a chemo-mechanical surface instability

The hypoxic conditions within avascular solid tumours may trigger the secretion of chemical factors, which diffuse to the nearby vasculature and promote the formation of new vessels eventually joining the tumour. Mathematical models of this process, known as tumour angiogenesis, have mainly investigated the formation of the new capillary networks using reaction-diffusion equations. Since angiogenesis involves the growth dynamics of the endothelial cells sprouting, we propose in this work an alternative mechanistic approach, developing a surface growth model for studying capillary formation and network dynamics. The model takes into account the proliferation of endothelial cells on the pre-existing capillary surface, coupled with the bulk diffusion of the vascular endothelial growth factor (VEGF). The thermo-dynamical consistency is imposed by means of interfacial and bulk balance laws. Finite element simulations show that both the morphology and the dynamics of the sprouting vessels are controlled by the bulk diffusion of VEGF and the chemo-mechanical and geometric properties at the capillary interface. Similarly to dendritic growth processes, we suggest that the emergence of tree-like vessel structures during tumour angiogenesis may result from the free boundary instability driven by competition between chemical and mechanical phenomena occurring at different length-scales.