Capillary interactions drive the self-organization of bacterial colonies.

Many bacteria inhabit thin layers of water on solid surfaces both naturally in soils or on hosts or textiles and in the lab on agar hydrogels. In these environments, cells experience capillary forces, yet an understanding of how these forces shape bacterial collective behaviors remains elusive. Here, we show that the water menisci formed around bacteria lead to capillary attraction between cells while still allowing them to slide past one another. We develop an experimental apparatus that allows us to control bacterial collective behaviors by varying the strength and range of capillary forces. Combining 3D imaging and cell tracking with agent-based modeling, we demonstrate that capillary attraction organizes rod-shaped bacteria into densely packed, nematic groups, and profoundly influences their collective dynamics and morphologies. Our results suggest that capillary forces may be a ubiquitous physical ingredient in shaping microbial communities in partially hydrated environments.

Local polar order controls mechanical stress and triggers layer formation in developing Myxococcus xanthus colonies.

Colonies of the social bacterium Myxococcus xanthus go through a morphological transition from a thin colony of cells to three-dimensional droplet-like fruiting bodies as a strategy to survive starvation. The biological pathways that control the decision to form a fruiting body have been studied extensively. However, the mechanical events that trigger the creation of multiple cell layers and give rise to droplet formation remain poorly understood. By measuring cell orientation, velocity, polarity, and force with cell-scale resolution, we reveal a stochastic local polar order in addition to the more obvious nematic order. Average cell velocity and active force at topological defects agree with predictions from active nematic theory, but their fluctuations are anomalously large due to polar active forces generated by the self-propelled rod-shaped cells. We find that M. xanthus cells adjust their reversal frequency to tune the magnitude of this local polar order, which in turn controls the mechanical stresses and triggers layer formation in the colonies.

Curvature induces active velocity waves in rotating spherical tissues.

The multicellular organization of diverse systems, including embryos, intestines, and tumors relies on coordinated cell migration in curved environments. In these settings, cells establish supracellular patterns of motion, including collective rotation and invasion. While such collective modes have been studied extensively in flat systems, the consequences of geometrical and topological constraints on collective migration in curved systems are largely unknown. Here, we discover a collective mode of cell migration in rotating spherical tissues manifesting as a propagating single-wavelength velocity wave. This wave is accompanied by an apparently incompressible supracellular flow pattern featuring topological defects as dictated by the spherical topology. Using a minimal active particle model, we reveal that this collective mode arises from the effect of curvature on the active flocking behavior of a cell layer confined to a spherical surface. Our results thus identify curvature-induced velocity waves as a mode of collective cell migration, impacting the dynamical organization of 3D curved tissues.

Self-assembly of tessellated tissue sheets by expansion and collision.

Tissues do not exist in isolation-they interact with other tissues within and across organs. While cell-cell interactions have been intensely investigated, less is known about tissue-tissue interactions. Here, we studied collisions between monolayer tissues with different geometries, cell densities, and cell types. First, we determine rules for tissue shape changes during binary collisions and describe complex cell migration at tri-tissue boundaries. Next, we propose that genetically identical tissues displace each other based on pressure gradients, which are directly linked to gradients in cell density. We present a physical model of tissue interactions that allows us to estimate the bulk modulus of the tissues from collision dynamics. Finally, we introduce TissEllate, a design tool for self-assembling complex tessellations from arrays of many tissues, and we use cell sheet engineering techniques to transfer these composite tissues like cellular films. Overall, our work provides insight into the mechanics of tissue collisions, harnessing them to engineer tissue composites as designable living materials.

Cellular Sensing Governs the Stability of Chemotactic Fronts.

In contexts ranging from embryonic development to bacterial ecology, cell populations migrate chemotactically along self-generated chemical gradients, often forming a propagating front. Here, we theoretically show that the stability of such chemotactic fronts to morphological perturbations is determined by limitations in the ability of individual cells to sense and thereby respond to the chemical gradient. Specifically, cells at bulging parts of a front are exposed to a smaller gradient, which slows them down and promotes stability, but they also respond more strongly to the gradient, which speeds them up and promotes instability. We predict that this competition leads to chemotactic fingering when sensing is limited at too low chemical concentrations. Guided by this finding and by experimental data on E. coli chemotaxis, we suggest that the cells' sensory machinery might have evolved to avoid these limitations and ensure stable front propagation. Finally, as sensing of any stimuli is necessarily limited in living and active matter in general, the principle of sensing-induced stability may operate in other types of directed migration such as durotaxis, electrotaxis, and phototaxis.

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.

Collective durotaxis of cohesive cell clusters on a stiffness gradient.

Many types of motile cells perform durotaxis, namely directed migration following gradients of substrate stiffness. Recent experiments have revealed that cell monolayers can migrate toward stiffer regions even when individual cells do not-a phenomenon known as collective durotaxis. Here, we address the spontaneous motion of finite cohesive cell monolayers on a stiffness gradient. We theoretically analyze a continuum active polar fluid model that has been tested in recent wetting assays of epithelial tissues and includes two types of active forces (cell-substrate traction and cell-cell contractility). The competition between the two active forces determines whether a cell monolayer spreads or contracts. Here, we show that this model generically predicts collective durotaxis, and that it features a variety of dynamical regimes as a result of the interplay between the spreading state and the global propagation, including sequential contraction and spreading of the monolayer as it moves toward higher stiffness. We solve the model exactly in some relevant cases, which provides both physical insights into the mechanisms of tissue durotaxis and spreading as well as a variety of predictions that could guide the design of future experiments.

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.

Size-dependent patterns of cell proliferation and migration in freely-expanding epithelia.

The coordination of cell proliferation and migration in growing tissues is crucial in development and regeneration but remains poorly understood. Here, we find that, while expanding with an edge speed independent of initial conditions, millimeter-scale epithelial monolayers exhibit internal patterns of proliferation and migration that depend not on the current but on the initial tissue size, indicating memory effects. Specifically, the core of large tissues becomes very dense, almost quiescent, and ceases cell-cycle progression. In contrast, initially-smaller tissues develop a local minimum of cell density and a tissue-spanning vortex. To explain vortex formation, we propose an active polar fluid model with a feedback between cell polarization and tissue flow. Taken together, our findings suggest that expanding epithelia decouple their internal and edge regions, which enables robust expansion dynamics despite the presence of size- and history-dependent patterns in the tissue interior.

Cells do not exist in isolation. Instead, they form tissues, where individual cells make contact with their neighbors and form microscopic 'architectures'. Epithelia are a type of tissue where cells are arranged in flat sheets, and are found in organs such as the lining of the kidney or the skin. Tissues need to grow, especially early in life. If tissues are damaged ­ for example, if the skin is cut or grazed ­ cells also need to divide (to create new healthy cells) and move as a group (to close the wound). Such coordinated motions result in cells exhibiting distinct group behaviors, similar to those observed within crowds of people or schools of fish. If coordination breaks down, problems can happen such as uncoordinated tissue growth seen in cancer. However, how cell movements are coordinated is still not fully understand. For example, researchers know that cells' positions within a group can determine how they behave, meaning that even the same type of cell could behave differently at the edge or center of a tissue. This suggests that the initial size and shape of a tissue should influence its subsequent growth and behavior; however, the nature of this influence is still largely unknown. Heinrich et al. therefore wanted to determine the differences in the way larger and smaller tissues grow. Microscope imaging was used to track the growth of circular, artificial tissues made from single-layered sheets of dog kidney cells grown in the laboratory. Comparing how quickly the tissues expanded revealed that the area of tissue circles that started out smaller increased at a much faster rate than that of tissue circles that were larger to begin with. This turned out to be because the edges of the tissues grew at a constant speed, independent of their initial size or shape, but circles with a smaller area have a larger proportion of cells on their edges. The motions of the cells at the center of the tissues had no effect on how the edges of the tissue grew. A final observation was that the way tissues of a given size behaved depended on whether they had grown to be that size, or they started off that big. These results shed light on how groups of cells interact in growing tissues. In the future, this information could be used to predict how different tissues grow over time, potentially helping scientists engineer better artificial tissues or organs for transplantation.

Nonuniform growth and surface friction determine bacterial biofilm morphology on soft substrates.

During development, organisms acquire three-dimensional (3D) shapes with important physiological consequences. While basic mechanisms underlying morphogenesis are known in eukaryotes, it is often difficult to manipulate them in vivo. To circumvent this issue, here we present a study of developing Vibrio cholerae biofilms grown on agar substrates in which the spatiotemporal morphological patterns were altered by varying the agar concentration. Expanding biofilms are initially flat but later undergo a mechanical instability and become wrinkled. To gain mechanistic insights into this dynamic pattern-formation process, we developed a model that considers diffusion of nutrients and their uptake by bacteria, bacterial growth/biofilm matrix production, mechanical deformation of both the biofilm and the substrate, and the friction between them. Our model shows quantitative agreement with experimental measurements of biofilm expansion dynamics, and it accurately predicts two distinct spatiotemporal patterns observed in the experiments-the wrinkles initially appear either in the peripheral region and propagate inward (soft substrate/low friction) or in the central region and propagate outward (stiff substrate/high friction). Our results, which establish that nonuniform growth and friction are fundamental determinants of stress anisotropy and hence biofilm morphology, are broadly applicable to bacterial biofilms with similar morphologies and also provide insight into how other bacterial biofilms form distinct wrinkle patterns. We discuss the implications of forming undulated biofilm morphologies, which may enhance the availability of nutrients and signaling molecules and serve as a "bet hedging" strategy.

Active wetting of epithelial tissues.

Development, regeneration and cancer involve drastic transitions in tissue morphology. In analogy with the behavior of inert fluids, some of these transitions have been interpreted as wetting transitions. The validity and scope of this analogy are unclear, however, because the active cellular forces that drive tissue wetting have been neither measured nor theoretically accounted for. Here we show that the transition between two-dimensional epithelial monolayers and three-dimensional spheroidal aggregates can be understood as an active wetting transition whose physics differs fundamentally from that of passive wetting phenomena. By combining an active polar fluid model with measurements of physical forces as a function of tissue size, contractility, cell-cell and cell-substrate adhesion, and substrate stiffness, we show that the wetting transition results from the competition between traction forces and contractile intercellular stresses. This competition defines a new intrinsic lengthscale that gives rise to a critical size for the wetting transition in tissues, a striking feature that has no counterpart in classical wetting. Finally, we show that active shape fluctuations are dynamically amplified during tissue dewetting. Overall, we conclude that tissue spreading constitutes a prominent example of active wetting - a novel physical scenario that may explain morphological transitions during tissue morphogenesis and tumor progression.

Active Fingering Instability in Tissue Spreading.

During the spreading of epithelial tissues, the advancing tissue front often develops fingerlike protrusions. Their resemblance to traditional viscous fingering patterns in driven fluids suggests that epithelial fingers could arise from an interfacial instability. However, the existence and physical mechanism of such a putative instability remain unclear. Here, based on an active polar fluid model for epithelial spreading, we analytically predict a generic instability of the tissue front. On the one hand, active cellular traction forces impose a velocity gradient that leads to an accelerated front, which is, thus, unstable to long-wavelength perturbations. On the other hand, contractile intercellular stresses typically dominate over surface tension in stabilizing short-wavelength perturbations. Finally, the finite range of hydrodynamic interactions in the tissue selects a wavelength for the fingering pattern, which is, thus, given by the smallest between the tissue size and the hydrodynamic screening length. Overall, we show that spreading epithelia experience an active fingering instability based on a simple kinematic mechanism. Moreover, our results underscore the crucial role of long-range hydrodynamic interactions in the dynamics of tissue morphology.

Role of Substrate Stiffness in Tissue Spreading: Wetting Transition and Tissue Durotaxis.

Living tissues undergo wetting transitions: On a surface, they can either form a dropletlike cell aggregate or spread as a monolayer of migrating cells. Tissue wetting depends not only on the chemical but also on the mechanical properties of the substrate. Here, we study the role of substrate stiffness in tissue spreading, which we describe by means of an active polar fluid model. Taking into account that cells exert larger active traction forces on stiffer substrates, we predict a tissue wetting transition at a critical substrate stiffness that decreases with tissue size. On substrates with a stiffness gradient, we find that the tissue spreads faster on the stiffer side. Furthermore, we show that the tissue can wet the substrate on the stiffer side while dewetting from the softer side. We also show that, by means of viscous forces transmitted across the tissue, the stiffer-side interface can transiently drag the softer-side interface toward increasing stiffness, against its spreading tendency. These two effects result in directed tissue migration up the stiffness gradient. This phenomenon-tissue durotaxis-can thus emerge both from dewetting on the soft side and from hydrodynamic interactions between the tissue interfaces. Overall, our work unveils mechanisms whereby substrate stiffness impacts the collective migration and the active wetting properties of living tissues, which are relevant in development, regeneration, and cancer.

Mixed-order phase transition in a colloidal crystal.

Mixed-order phase transitions display a discontinuity in the order parameter like first-order transitions yet feature critical behavior like second-order transitions. Such transitions have been predicted for a broad range of equilibrium and nonequilibrium systems, but their experimental observation has remained elusive. Here, we analytically predict and experimentally realize a mixed-order equilibrium phase transition. Specifically, a discontinuous solid-solid transition in a 2D crystal of paramagnetic colloidal particles is induced by a magnetic field [Formula: see text] At the transition field [Formula: see text], the energy landscape of the system becomes completely flat, which causes diverging fluctuations and correlation length [Formula: see text] Mean-field critical exponents are predicted, since the upper critical dimension of the transition is [Formula: see text] Our colloidal system provides an experimental test bed to probe the unconventional properties of mixed-order phase transitions.

Fluidization and Active Thinning by Molecular Kinetics in Active Gels.

We derive the constitutive equations of an active polar gel from a model for the dynamics of elastic molecules that link polar elements. Molecular binding kinetics induces the fluidization of the material, giving rise to Maxwell viscoelasticity and, provided that detailed balance is broken, to the generation of active stresses. We give explicit expressions for the transport coefficients of active gels in terms of molecular properties, including nonlinear contributions on the departure from equilibrium. In particular, when activity favors linker unbinding, we predict a decrease of viscosity with activity-active thinning-of kinetic origin, which could explain some experimental results on the cell cortex. By bridging the molecular and hydrodynamic scales, our results could help understand the interplay between molecular perturbations and the mechanics of cells and tissues.

Emergent structures and dynamics of cell colonies by contact inhibition of locomotion.

Cells in tissues can organize into a broad spectrum of structures according to their function. Drastic changes of organization, such as epithelial-mesenchymal transitions or the formation of spheroidal aggregates, are often associated either to tissue morphogenesis or to cancer progression. Here, we study the organization of cell colonies by means of simulations of self-propelled particles with generic cell-like interactions. The interplay between cell softness, cell-cell adhesion, and contact inhibition of locomotion (CIL) yields structures and collective dynamics observed in several existing tissue phenotypes. These include regular distributions of cells, dynamic cell clusters, gel-like networks, collectively migrating monolayers, and 3D aggregates. We give analytical predictions for transitions between noncohesive, cohesive, and 3D cell arrangements. We explicitly show how CIL yields an effective repulsion that promotes cell dispersal, thereby hindering the formation of cohesive tissues. Yet, in continuous monolayers, CIL leads to collective cell motion, ensures tensile intercellular stresses, and opposes cell extrusion. Thus, our work highlights the prominent role of CIL in determining the emergent structures and dynamics of cell colonies.

Formation of metastable phases by spinodal decomposition.

Metastable phases may be spontaneously formed from other metastable phases through nucleation. Here we demonstrate the spontaneous formation of a metastable phase from an unstable equilibrium by spinodal decomposition, which leads to a transient coexistence of stable and metastable phases. This phenomenon is generic within the recently introduced scenario of the landscape-inversion phase transitions, which we experimentally realize as a structural transition in a colloidal crystal. This transition exhibits a rich repertoire of new phase-ordering phenomena, including the coexistence of two equilibrium phases connected by two physically different interfaces. In addition, this scenario enables the control of sizes and lifetimes of metastable domains. Our findings open a new setting that broadens the fundamental understanding of phase-ordering kinetics, and yield new prospects of applications in materials science.

Bleb Nucleation through Membrane Peeling.

We study the nucleation of blebs, i.e., protrusions arising from a local detachment of the membrane from the cortex of a cell. Based on a simple model of elastic linkers with force-dependent kinetics, we show that bleb nucleation is governed by membrane peeling. By this mechanism, the growth or shrinkage of a detached membrane patch is completely determined by the linker kinetics, regardless of the energetic cost of the detachment. We predict the critical nucleation radius for membrane peeling and the corresponding effective energy barrier. These may be typically smaller than those predicted by classical nucleation theory, implying a much faster nucleation. We also perform simulations of a continuum stochastic model of membrane-cortex adhesion to obtain the statistics of bleb nucleation times as a function of the stress on the membrane. The determinant role of membrane peeling changes our understanding of bleb nucleation and opens new directions in the study of blebs.

Model for probing membrane-cortex adhesion by micropipette aspiration and fluctuation spectroscopy.

We propose a model for membrane-cortex adhesion that couples membrane deformations, hydrodynamics, and kinetics of membrane-cortex ligands. In its simplest form, the model gives explicit predictions for the critical pressure for membrane detachment and for the value of adhesion energy. We show that these quantities exhibit a significant dependence on the active acto-myosin stresses. The model provides a simple framework to access quantitative information on cortical activity by means of micropipette experiments. We also extend the model to incorporate fluctuations and show that detailed information on the stability of membrane-cortex coupling can be obtained by a combination of micropipette aspiration and fluctuation spectroscopy measurements.

Landscape-inversion phase transition in dipolar colloids: tuning the structure and dynamics of 2D crystals.

We study the 2D crystalline phases of paramagnetic colloidal particles with dipolar interactions and constrained on a periodic substrate. Combining theory, simulation, and experiments, we demonstrate a new scenario of first-order phase transitions that occurs via a complete inversion of the energy landscape, featuring nonconventional properties that allow for (i) tuning of crystal symmetry, (ii) control of dynamical properties of different crystalline orders via tuning of their relative stability with an external magnetic field, (iii) an equivalent but independent control of the same dynamic properties via temporal modulations of that field, and (iv) nonstandard phase-ordering kinetics involving spontaneous formation of transient metastable domains.