Emergence of Active Nematic Behavior in Monolayers of Isotropic Cells

Active flows and deformable surfaces in development

We review progress in active hydrodynamic descriptions of flowing media on curved and deformable manifolds: the state-of-the-art in continuum descriptions of single-layers of epithelial and/or other tissues during development. First, after a brief overview of activity, flows and hydrodynamic descriptions, we highlight the generic challenge of identifying the dependence on dynamical variables of so-called active kinetic coefficients- active counterparts to dissipative Onsager coefficients. We go on to describe some of the subtleties concerning how curvature and active flows interact, and the issues that arise when surfaces are deformable. We finish with a broad discussion around the utility of such theories in developmental biology. This includes limitations to analytical techniques, challenges associated with numerical integration, fitting-to-data and inference, and potential tools for the future, such as discrete differential geometry.

Force-mediated cellular anisotropy and plasticity dictate the elongation dynamics of embryos

We developed a unified dynamic model to explain how cellular anisotropy and plasticity, induced by alignment and severing/rebundling of actin filaments, dictate the elongation dynamics of Caenorhabditis elegans embryos. It was found that the gradual alignment of F-actins must be synchronized with the development of intracellular forces for the embryo to elongate, which is then further sustained by muscle contraction-triggered plastic deformation of cells. In addition, we showed that preestablished anisotropy is essential for the proper onset of the process while defects in the integrity or bundling kinetics of actin bundles result in abnormal embryo elongation, all in good agreement with experimental observations.

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. The patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

Fluctuations can induce local nematic order and extensile stress in monolayers of motile cells

Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix.

Machine learning active-nematic hydrodynamics

Hydrodynamic theories effectively describe many-body systems out of equilibrium in terms of a few macroscopic parameters. However, such parameters are difficult to determine from microscopic information. Seldom is this challenge more apparent than in active matter, where the hydrodynamic parameters are in fact fields that encode the distribution of energy-injecting microscopic components. Here, we use active nematics to demonstrate that neural networks can map out the spatiotemporal variation of multiple hydrodynamic parameters and forecast the chaotic dynamics of these systems. We analyze biofilament/molecular-motor experiments with microtubule/kinesin and actin/myosin complexes as computer vision problems. Our algorithms can determine how activity and elastic moduli change as a function of space and time, as well as adenosine triphosphate (ATP) or motor concentration. The only input needed is the orientation of the biofilaments and not the coupled velocity field which is harder to access in experiments. We can also forecast the evolution of these chaotic many-body systems solely from image sequences of their past using a combination of autoencoders and recurrent neural networks with residual architecture. In realistic experimental setups for which the initial conditions are not perfectly known, our physics-inspired machine-learning algorithms can surpass deterministic simulations. Our study paves the way for artificial-intelligence characterization and control of coupled chaotic fields in diverse physical and biological systems, even in the absence of knowledge of the underlying dynamics.

Epithelial colonies in vitro elongate through collective effects

Epithelial tissues of the developing embryos elongate by different mechanisms, such as neighbor exchange, cell elongation, and oriented cell division. Since autonomous tissue self-organization is influenced by external cues such as morphogen gradients or neighboring tissues, it is difficult to distinguish intrinsic from directed tissue behavior. The mesoscopic processes leading to the different mechanisms remain elusive. Here, we study the spontaneous elongation behavior of spreading circular epithelial colonies in vitro. By quantifying deformation kinematics at multiple scales, we report that global elongation happens primarily due to cell elongations, and its direction correlates with the anisotropy of the average cell elongation. By imposing an external time-periodic stretch, the axis of this global symmetry breaking can be modified and elongation occurs primarily due to orientated neighbor exchange. These different behaviors are confirmed using a vertex model for collective cell behavior, providing a framework for understanding autonomous tissue elongation and its origins.

The vortex-driven dynamics of droplets within droplets

Understanding the fluid-structure interaction is crucial for an optimal design and manufacturing of soft mesoscale materials. Multi-core emulsions are a class of soft fluids assembled from cluster configurations of deformable oil-water double droplets (cores), often employed as building-blocks for the realisation of devices of interest in bio-technology, such as drug-delivery, tissue engineering and regenerative medicine. Here, we study the physics of multi-core emulsions flowing in microfluidic channels and report numerical evidence of a surprisingly rich variety of driven non-equilibrium states (NES), whose formation is caused by a dipolar fluid vortex triggered by the sheared structure of the flow carrier within the microchannel. The observed dynamic regimes range from long-lived NES at low core-area fraction, characterised by a planetary-like motion of the internal drops, to short-lived ones at high core-area fraction, in which a pre-chaotic motion results from multi-body collisions of inner drops, as combined with self-consistent hydrodynamic interactions. The onset of pre-chaotic behavior is marked by transitions of the cores from one vortex to another, a process that we interpret as manifestations of the system to maximize its entropy by filling voids, as they arise dynamically within the capsule.

Integer topological defects of cell monolayers: Mechanics and flows

Monolayers of anisotropic cells exhibit long-ranged orientational order and topological defects. During the development of organisms, orientational order often influences morphogenetic events. However, the linkage between the mechanics of cell monolayers and topological defects remains largely unexplored. This holds specifically at the timescales relevant for tissue morphogenesis. Here, we build on the physics of liquid crystals to determine material parameters of cell monolayers. In particular, we use a hydrodynamical description of an active polar fluid to study the steady-state mechanical patterns at integer topological defects. Our description includes three distinct sources of activity: traction forces accounting for cell-substrate interactions as well as anisotropic and isotropic active nematic stresses accounting for cell-cell interactions. We apply our approach to C2C12 cell monolayers in small circular confinements, which form isolated aster or spiral topological defects. By analyzing the velocity and orientational order fields in spirals as well as the forces and cell number density fields in asters, we determine mechanical parameters of C2C12 cell monolayers. Our work shows how topological defects can be used to fully characterize the mechanical properties of biological active matter.

Reversible elastic phase field approach and application to cell monolayers

Motion and generation of forces by single cells and cell collectives are essential elements of many biological processes, including development, wound healing and cancer cell migration. Quantitative wound healing assays have demonstrated that cell monolayers can be both dynamic and elastic at the same time. However, it is very challenging to model this combination with conventional approaches. Here we introduce an elastic phase field approach that allows us to predict the dynamics of elastic sheets under the action of active stresses and localized forces, e.g. from leader cells. Our method ensures elastic reversibility after release of forces. We demonstrate its potential by studying several paradigmatic situations and geometries relevant for single cells and cell monolayers, including elastic bars, contractile discs and expanding monolayers with leader cells.

Active inter-cellular forces in collective cell motility

The collective behaviour of confluent cell sheets is strongly influenced both by polar forces, arising through cytoskeletal propulsion, and by active inter-cellular forces, which are mediated by interactions across cell-cell junctions. We use a phase-field model to explore the interplay between these two contributions and compare the dynamics of a cell sheet when the polarity of the cells aligns to (i) their main axis of elongation, (ii) their velocity and (iii) when the polarity direction executes a persistent random walk. In all three cases, we observe a sharp transition from a jammed state (where cell rearrangements are strongly suppressed) to a liquid state (where the cells can move freely relative to each other) when either the polar or the inter-cellular forces are increased. In addition, for case (ii) only, we observe an additional dynamical state, flocking (solid or liquid), where the majority of the cells move in the same direction. The flocking state is seen for strong polar forces, but is destroyed as the strength of the inter-cellular activity is increased.

Solid-Liquid Transition of Deformable and Overlapping Active Particles

Experiments and theory have shown that cell monolayers and epithelial tissues exhibit solid-liquid and glass-liquid transitions. These transitions are biologically relevant to our understanding of embryonic development, wound healing, and cancer. Current models of confluent epithelia have focused on the role of cell shape, with less attention paid to cell extrusion, which is key for maintaining homeostasis in biological tissue. Here, we use a multiphase field model to study the solid-liquid transition in a confluent monolayer of deformable cells. Cell overlap is allowed and provides a way for modeling the precursor for extrusion. When cells overlap rather than deform, we find that the melting transition changes from continuous to first order like, and that there is an intermittent regime close to the transition, where solid and liquid states alternate over time. By studying the dynamics of five- and sevenfold disclinations in the hexagonal lattice formed by the cell centers, we observe that these correlate with spatial fluctuations in the cellular overlap, and that cell extrusion tends to initiate near fivefold disclinations.

Active phase field crystal systems with inertial delay and underdamped dynamics

Active matter systems often are well approximated as overdamped, meaning that any inertial momentum is immediately dissipated by the environment. On the other hand, especially for macroscopic systems but also for many mesoscopic ones particle mass can become relevant for the dynamics. For such systems we recently proposed an underdamped continuum model which captures translationally inertial dynamics via two contributions. First, convection and second a damping time scale of inertial motion. In this paper, we ask how both of these features influence the collective behavior compared to overdamped dynamics by studying the example of the active phase field crystal model. We first focus on the case of suppressed convection to study the role of the damping time. We quantify that the relaxation process to the steady collective motion state is considerably prolonged with damping time due to the increasing occurrence of transient groups of circularly moving density peaks. Finally, we illustrate the fully underdamped case with convection. Instead of collective motion of density peaks we then find a coexistence of constant high and low density phases reminiscent of motility-induced phase separation.

Mean field approach of dynamical pattern formation in underdamped active matter with short-ranged alignment and distant anti-alignment interactions

Many active matter systems, especially on the microscopic scale, are well approximated as overdamped, meaning that any inertial momentum is immediately dissipated by the environment. On the other hand, especially for macroscopic active systems but also for many mesoscopic systems the time scale of translational inertial motion can become large enough to be relevant for the dynamics. This raises the question how collective dynamics and the resulting states in active matter are influenced by inertia. Therefore, we propose a coarse-grained continuum model for underdamped active matter based on a mean field description for passive systems. Furthermore, as an example, we apply the model to a system with interactions that support an alignment on short distances and an anti-alignment on longer length scales as known in the context of pattern formation due to orientational interactions. Our numerical calculations of the under- and overdamped dynamics both predict a structured laning state. However, activity induced convective flows that are only present in the underdamped model destabilize this state when the anti-alignment is weakened, leading to a collective motion state which does not occur in the overdamped limit. A turbulent transition regime between the two states can be characterized by strong density fluctuations and the absence of global ordering.

The 2020 motile active matter roadmap

Activity and autonomous motion are fundamental in living and engineering systems. This has stimulated the new field of 'active matter' in recent years, which focuses on the physical aspects of propulsion mechanisms, and on motility-induced emergent collective behavior of a larger number of identical agents. The scale of agents ranges from nanomotors and microswimmers, to cells, fish, birds, and people. Inspired by biological microswimmers, various designs of autonomous synthetic nano- and micromachines have been proposed. Such machines provide the basis for multifunctional, highly responsive, intelligent (artificial) active materials, which exhibit emergent behavior and the ability to perform tasks in response to external stimuli. A major challenge for understanding and designing active matter is their inherent nonequilibrium nature due to persistent energy consumption, which invalidates equilibrium concepts such as free energy, detailed balance, and time-reversal symmetry. Unraveling, predicting, and controlling the behavior of active matter is a truly interdisciplinary endeavor at the interface of biology, chemistry, ecology, engineering, mathematics, and physics. The vast complexity of phenomena and mechanisms involved in the self-organization and dynamics of motile active matter comprises a major challenge. Hence, to advance, and eventually reach a comprehensive understanding, this important research area requires a concerted, synergetic approach of the various disciplines. The 2020 motile active matter roadmap of Journal of Physics: Condensed Matter addresses the current state of the art of the field and provides guidance for both students as well as established scientists in their efforts to advance this fascinating area.

Active matter invasion

Biologically active materials such as bacterial biofilms and eukaryotic cells thrive in confined micro-spaces. Here, we show through numerical simulations that confinement can serve as a mechanical guidance to achieve distinct modes of collective invasion when combined with growth dynamics and the intrinsic activity of biological materials. We assess the dynamics of the growing interface and classify these collective modes of invasion based on the activity of the constituent particles of the growing matter. While at small and moderate activities the active material grows as a coherent unit, we find that blobs of active material collectively detach from the cohort above a well-defined activity threshold. We further characterise the mechanical mechanisms underlying the crossovers between different modes of invasion and quantify their impact on the overall invasion speed.