Topological defects in epithelia govern cell death and extrusion

Voltage-dependent volume regulation controls epithelial cell extrusion and morphology

Epithelial cells work collectively to provide a protective barrier, yet also turn over rapidly by cell death and division. If the number of dying cells does not match those dividing, the barrier would vanish, or tumors can form. Mechanical forces and the stretch-activated ion channel (SAC) Piezo1 link both processes; stretch promotes cell division and crowding triggers cell death by initiating live cell extrusion1,2. However, it was not clear how particular cells within a crowded region are selected for extrusion. Here, we show that individual cells transiently shrink via water loss before they extrude. Artificially inducing cell shrinkage by increasing extracellular osmolarity is sufficient to induce cell extrusion. Pre-extrusion cell shrinkage requires the voltage-gated potassium channels Kv1.1 and Kv1.2 and the chloride channel SWELL1, upstream of Piezo1. Activation of these voltage-gated channels requires the mechano-sensitive Epithelial Sodium Channel, ENaC, acting as the earliest crowd-sensing step. Imaging with a voltage dye indicated that epithelial cells lose membrane potential as they become crowded and smaller, yet those selected for extrusion are markedly more depolarized than their neighbours. Loss of any of these channels in crowded conditions causes epithelial buckling, highlighting an important role for voltage and water regulation in controlling epithelial shape as well as extrusion. Thus, ENaC causes cells with similar membrane potentials to slowly shrink with compression but those with reduced membrane potentials to be eliminated by extrusion, suggesting a chief driver of cell death stems from insufficient energy to maintain cell membrane potential.

Hydrodynamic Enhancement of p-atic Defect Dynamics

We investigate numerically and analytically the effects of hydrodynamics on the dynamics of topological defects in p-atic liquid crystals, i.e., two-dimensional liquid crystals with p-fold rotational symmetry. Importantly, we find that hydrodynamics fuels a generic passive self-propulsion mechanism for defects of winding number s=(p-1)/p and arbitrary p. Strikingly, we discover that hydrodynamics always accelerates the annihilation dynamics of pairs of ±1/p defects and that, contrary to expectations, this effect increases with p. Our Letter paves the way toward understanding cell intercalation and other remodeling events in epithelial layers.

Bioengineering of High Cell Density Tissues with Hierarchical Vascular Networks for Ex Vivo Whole Organs

The development of tissue-like structures such as cell sheets, spheroids, and organoids has contributed to progress in regenerative medicine. Simultaneous achievement of scale up and high cell density of these tissues is challenging because sufficient oxygen cannot be supplied to the inside of large, high cell density tissues. Here, in vitro fabrication of vessels to supply oxygen to the inside of millimeter-sized scaffold-free tissues whose cell density is ≈200 million cells mL^-1 , corresponding to those of native tissues, is shown. Hierarchical vascular networks by anastomosis of capillaries and a large vessel are essential for oxygen supply, whereas a large vessel or capillary networks alone make negligible contributions to the supply. The hierarchical vascular networks are formed by a top-down approach, which offers a new option for ex vivo whole organs without decellularization and 3D-bioprinting.

Curvature in Biological Systems: Its Quantification, Emergence, and Implications across the Scales

Surface curvature both emerges from, and influences the behavior of, living objects at length scales ranging from cell membranes to single cells to tissues and organs. The relevance of surface curvature in biology is supported by numerous experimental and theoretical investigations in recent years. In this review, first, a brief introduction to the key ideas of surface curvature in the context of biological systems is given and the challenges that arise when measuring surface curvature are discussed. Giving an overview of the emergence of curvature in biological systems, its significance at different length scales becomes apparent. On the other hand, summarizing current findings also shows that both single cells and entire cell sheets, tissues or organisms respond to curvature by modulating their shape and their migration behavior. Finally, the interplay between the distribution of morphogens or micro-organisms and the emergence of curvature across length scales is addressed with examples demonstrating these key mechanistic principles of morphogenesis. Overall, this review highlights that curved interfaces are not merely a passive by-product of the chemical, biological, and mechanical processes but that curvature acts also as a signal that co-determines these processes.

Crisscross multilayering of cell sheets

Hydrostatic skeletons such as the Hydra's consist of two stacked layers of muscle cells perpendicularly oriented. In vivo, these bilayers first assemble, and then the muscle fibers of both layers develop and organize with this crisscross orientation. In the present work, we identify an alternative mechanism of crisscross bilayering of myoblasts in vitro, which results from the prior local organization of these active cells in the initial monolayer. The myoblast sheet can be described as a contractile active nematic in which, as expected, most of the +1/2 topological defects associated with this nematic order self-propel. However, as a result of the production of extracellular matrix (ECM) by the cells, a subpopulation of these comet-like defects does not show any self-propulsion. Perpendicular bilayering occurs at these stationary defects. Cells located at the head of these defects converge toward their core where they accumulate until they start migrating on top of the tail of the first layer, while the tail cells migrate in the opposite direction under the head. Since the cells keep their initial orientations, the two stacked layers end up perpendicularly oriented. This concerted process leading to a crisscross bilayering is mediated by the secretion of ECM.

E-cadherin-dependent coordinated epithelial rotation on a two-dimensional discoidal pattern

In vivo, cells collectively migrate in a variety of developmental and pathological contexts. Coordinated epithelial rotation represents a unique type of collective cell migrations, which has been modeled in vitro under spatially confined conditions. Although it is known that the coordinated rotation depends on intercellular interactions, the contribution of E-cadherin, a major cell-cell adhesion molecule, has not been directly addressed on two-dimensional (2D) confined substrates. Here, using well-controlled fibronectin-coated surfaces, we tracked and compared the migratory behaviors of MDCK cells expressing or lacking E-cadherin. We observed that wild-type MDCK II cells exhibited persistent and coordinated rotations on discoidal patterns, while E-cadherin knockout cells migrated in a less coordinated manner without large-scale rotation. Our comparison of the collective dynamics between these two cell types revealed a series of changes in migratory behavior caused by the loss of E-cadherin, including a decreased global migration speed, less regularity in quantified coordination, and increased average density of topological defects. Taken together, these data demonstrate that spontaneous initiation of collective epithelial rotations depends on E-cadherin under 2D discoidal confinements.

Inhomogeneous mechanotransduction defines the spatial pattern of apoptosis-induced compensatory proliferation

The number of cells in tissues is controlled by cell division and cell death, and its misregulation could lead to pathological conditions such as cancer. To maintain the cell numbers, a cell-elimination process called apoptosis also stimulates the proliferation of neighboring cells. This mechanism, apoptosis-induced compensatory proliferation, was originally described more than 40 years ago. Although only a limited number of the neighboring cells need to divide to compensate for the apoptotic cell loss, the mechanisms that select cells to divide have remained elusive. Here, we found that spatial inhomogeneity in Yes-associated protein (YAP)-mediated mechanotransduction in neighboring tissues determines the inhomogeneity of compensatory proliferation in Madin-Darby canine kidney (MDCK) cells. Such inhomogeneity arises from the non-uniform distribution of nuclear size and the non-uniform pattern of mechanical force applied to neighboring cells. Our findings from a mechanical perspective provide additional insight into how tissues precisely maintain homeostasis.

The crucial role of adhesion in the transmigration of active droplets through interstitial orifices

Active fluid droplets are a class of soft materials exhibiting autonomous motion sustained by an energy supply. Such systems have been shown to capture motility regimes typical of biological cells and are ideal candidates as building-block for the fabrication of soft biomimetic materials of interest in pharmacology, tissue engineering and lab on chip devices. While their behavior is well established in unconstrained environments, much less is known about their dynamics under strong confinement. Here, we numerically study the physics of a droplet of active polar fluid migrating within a microchannel hosting a constriction with adhesive properties, and report evidence of a striking variety of dynamic regimes and morphological features, whose properties crucially depend upon droplet speed and elasticity, degree of confinement within the constriction and adhesiveness to the pore. Our results suggest that non-uniform adhesion forces are instrumental in enabling the crossing through narrow orifices, in contrast to larger gaps where a careful balance between speed and elasticity is sufficient to guarantee the transition. These observations may be useful for improving the design of artificial micro-swimmers, of interest in material science and pharmaceutics, and potentially for cell sorting in microfluidic devices.

Small hand-designed convolutional neural networks outperform transfer learning in automated cell shape detection in confluent tissues

Mechanical cues such as stresses and strains are now recognized as essential regulators in many biological processes like cell division, gene expression or morphogenesis. Studying the interplay between these mechanical cues and biological responses requires experimental tools to measure these cues. In the context of large scale tissues, this can be achieved by segmenting individual cells to extract their shapes and deformations which in turn inform on their mechanical environment. Historically, this has been done by segmentation methods which are well known to be time consuming and error prone. In this context however, one doesn't necessarily require a cell-level description and a coarse-grained approach can be more efficient while using tools different from segmentation. The advent of machine learning and deep neural networks has revolutionized the field of image analysis in recent years, including in biomedical research. With the democratization of these techniques, more and more researchers are trying to apply them to their own biological systems. In this paper, we tackle a problem of cell shape measurement thanks to a large annotated dataset. We develop simple Convolutional Neural Networks (CNNs) which we thoroughly optimize in terms of architecture and complexity to question construction rules usually applied. We find that increasing the complexity of the networks rapidly no longer yields improvements in performance and that the number of kernels in each convolutional layer is the most important parameter to achieve good results. In addition, we compare our step-by-step approach with transfer learning and find that our simple, optimized CNNs give better predictions, are faster in training and analysis and don't require more technical knowledge to be implemented. Overall, we offer a roadmap to develop optimized models and argue that we should limit the complexity of such models. We conclude by illustrating this strategy on a similar problem and dataset.

Actin polymerisation and crosslinking drive left-right asymmetry in single cell and cell collectives

Deviations from mirror symmetry in the development of bilateral organisms are common but the mechanisms of initial symmetry breaking are insufficiently understood. The actin cytoskeleton of individual cells self-organises in a chiral manner, but the molecular players involved remain essentially unidentified and the relationship between chirality of an individual cell and cell collectives is unclear. Here, we analysed self-organisation of the chiral actin cytoskeleton in individual cells on circular or elliptical patterns, and collective cell alignment in confined microcultures. Screening based on deep-learning analysis of actin patterns identified actin polymerisation regulators, depletion of which suppresses chirality (mDia1) or reverses chirality direction (profilin1 and CapZβ). The reversed chirality  is mDia1-independent but requires the function of actin-crosslinker α-actinin1. A robust correlation between the effects of a variety of actin assembly regulators on chirality of individual cells and cell collectives is revealed. Thus, actin-driven cell chirality may underlie tissue and organ asymmetry.

Structure and Rheology in Vertex Models under Cell-Shape-Dependent Active Stresses

Biological cells can actively tune their intracellular architecture according to their overall shape. Here we explore the rheological implication of such coupling in a minimal model of a dense cellular material where each cell exerts an active mechanical stress along its axis of elongation. Increasing the active stress amplitude leads to several transitions. An initially hexagonal crystal motif is first destabilized into a solid with anisotropic cells whose shear modulus eventually vanishes at a first critical activity. Increasing activity beyond this first critical value, we find a re-entrant transition to a regime with finite hexatic order and finite shear modulus, in which cells arrange according to a rhombile pattern with periodically arranged rosette structures. The shear modulus vanishes again at a third threshold beyond which spontaneous tissue flows and topological defects of the nematic cell shape field arise. Flow and stress fields around the defects agree with active nematic theory, with either contractile or extensile signs, as also observed in several epithelial tissue experiments.

Mechanobiological approaches to synthetic morphogenesis: learning by building

Tissue morphogenesis occurs in a complex physicochemical microenvironment with limited experimental accessibility. This often prevents a clear identification of the processes that govern the formation of a given functional shape. By applying state-of-the-art methods to minimal tissue systems, synthetic morphogenesis aims to engineer the discrete events that are necessary and sufficient to build specific tissue shapes. Here, we review recent advances in synthetic morphogenesis, highlighting how a combination of microfabrication and mechanobiology is fostering our understanding of how tissues are built.

Activity-driven tissue alignment in proliferating spheroids

We extend the continuum theory of active nematic fluids to study cell flows and tissue dynamics inside multicellular spheroids, spherical, self-assembled aggregates of cells that are widely used as model systems to study tumour dynamics. Cells near the surface of spheroids have better access to nutrients and therefore proliferate more rapidly than those in the resource-depleted core. Using both analytical arguments and three-dimensional simulations, we find that the proliferation gradients result in flows and in gradients of activity both of which can align the orientation axis of cells inside the aggregates. Depending on environmental conditions and the intrinsic tissue properties, we identify three distinct alignment regimes: spheroids in which all the cells align either radially or tangentially to the surface throughout the aggregate and spheroids with angular cell orientation close to the surface and radial alignment in the core. The continuum description of tissue dynamics inside spheroids not only allows us to infer dynamic cell parameters from experimentally measured cell alignment profiles, but more generally motivates novel mechanisms for controlling the alignment of cells within aggregates which has been shown to influence the mechanical properties and invasive capabilities of tumors.

Collective rotational motion of freely expanding T84 epithelial cell colonies

Coordinated rotational motion is an intriguing, yet still elusive mode of collective cell migration, which is relevant in pathological and morphogenetic processes. Most of the studies on this topic have been carried out on epithelial cells plated on micropatterned substrates, where cell motion is confined in regions of well-defined shapes coated with extracellular matrix adhesive proteins. The driver of collective rotation in such conditions has not been clearly elucidated, although it has been speculated that spatial confinement can play an essential role in triggering cell rotation. Here, we study the growth of epithelial cell colonies freely expanding (i.e. with no physical constraints) on the surface of cell culture plates and focus on collective cell rotation in such conditions, a case which has received scarce attention in the literature. One of the main findings of our work is that coordinated cell rotation spontaneously occurs in cell clusters in the free growth regime, thus implying that cell confinement is not necessary to elicit collective rotation as previously suggested. The extent of collective rotation was size and shape dependent: a highly coordinated disc-like rotation was found in small cell clusters with a round shape, while collective rotation was suppressed in large irregular cell clusters generated by merging of different clusters in the course of their growth. The angular motion was persistent in the same direction, although clockwise and anticlockwise rotations were equally likely to occur among different cell clusters. Radial cell velocity was quite low as compared to the angular velocity, in agreement with the free expansion regime where cluster growth is essentially governed by cell proliferation. A clear difference in morphology was observed between cells at the periphery and the ones in the core of the clusters, the former being more elongated and spread out as compared to the latter. Overall, our results, to our knowledge, provide the first quantitative and systematic evidence that coordinated cell rotation does not require a spatial confinement and occurs spontaneously in freely expanding epithelial cell colonies, possibly as a mechanism for the system.

Emergence of macroscopic directional motion of deformable active cells in confined structures

There is now growing evidence of collective turbulentlike motion of cells in dense tissues. However, how to control and harness this collective motion is an open question. We investigate the transport of deformable active cells in a periodically asymmetric channel by using a phase-field model. We demonstrate that collective turbulent-like motion of cells can power and steer the macroscopic directional motion through the ratchet channel. The active intercellular forces proportional to the deformation of cells can break thermodynamical equilibrium and induce the directional motion. This directional motion is caused by the ratchet effect rather than the spontaneous symmetry breaking. The motion direction is determined by the asymmetry of the channel. Remarkably, there exits an optimal nonequilibrium driving (depending on the active strength, the elasticity, and the packing fraction) at which the average velocity reaches the maximum. In addition, the optimized packing fraction and the optimized minimum width of the channel can facilitate the directional motion of cells. Our findings are relevant to understanding how macroscopic directional motion relates to the local force transmission mediated by cell-cell contacts in cellular monolayers.

Alignment rule and geometric confinement lead to stability of a vortex in active flow

Vortices are hallmarks of a wide range of nonequilibrium phenomena in fluids at multiple length scales. In this work, we numerically study the whirling motion of self-propelled soft point particles confined in circular domain, and aim at addressing the stability issue of the coherent vortex structure. By the combination of dynamical and statistical analysis at the individual particle level, we reveal the persistence of the whirling motion resulting from the subtle competition of activity and geometric confinement. In the stable whirling motion, the scenario of the coexistence of the irregular microscopic motions of individual particles and the regular global whirling motion is fundamentally different from the motion of a vortex in passive fluid. Possible orientational order coexisting with the whirling are further explored. This work shows the stability mechanism of vortical dynamics in active media under the alignment rule in confined space and may have implications in creating and harnessing macroscale coherent dynamical states by tuning the confining geometry.

Active Forces in Confluent Cell Monolayers

We use a computational phase-field model together with analytical analysis to study how intercellular active forces can mediate individual cell morphology and collective motion in a confluent cell monolayer. We explore the regime where intercellular forces dominate the tissue dynamics, and polar forces are negligible. Contractile intercellular interactions lead to cell elongation, nematic ordering, and active turbulence characterized by motile topological defects. Extensile interactions result in frustration, and perpendicular cell orientations become more prevalent. Furthermore, we show that contractile behavior can change to extensile behavior if anisotropic fluctuations in cell shape are considered.

Active morphogenesis of patterned epithelial shells

Shape transformations of epithelial tissues in three dimensions, which are crucial for embryonic development or in vitro organoid growth, can result from active forces generated within the cytoskeleton of the epithelial cells. How the interplay of local differential tensions with tissue geometry and with external forces results in tissue-scale morphogenesis remains an open question. Here, we describe epithelial sheets as active viscoelastic surfaces and study their deformation under patterned internal tensions and bending moments. In addition to isotropic effects, we take into account nematic alignment in the plane of the tissue, which gives rise to shape-dependent, anisotropic active tensions and bending moments. We present phase diagrams of the mechanical equilibrium shapes of pre-patterned closed shells and explore their dynamical deformations. Our results show that a combination of nematic alignment and gradients in internal tensions and bending moments is sufficient to reproduce basic building blocks of epithelial morphogenesis, including fold formation, budding, neck formation, flattening, and tubulation.

Mechanical stress driven by rigidity sensing governs epithelial stability

Epithelia act as a barrier against environmental stress and abrasion and in vivo they are continuously exposed to environments of various mechanical properties. The impact of this environment on epithelial integrity remains elusive. By culturing epithelial cells on 2D hydrogels, we observe a loss of epithelial monolayer integrity through spontaneous hole formation when grown on soft substrates. Substrate stiffness triggers an unanticipated mechanical switch of epithelial monolayers from tensile on soft to compressive on stiff substrates. Through active nematic modelling, we find that spontaneous half-integer defect formation underpinning large isotropic stress fluctuations initiate hole opening events. Our data show that monolayer rupture due to high tensile stress is promoted by the weakening of cell-cell junctions that could be induced by cell division events or local cellular stretching. Our results show that substrate stiffness provides feedback on monolayer mechanical state and that topological defects can trigger stochastic mechanical failure, with potential application towards a mechanistic understanding of compromised epithelial integrity during immune response and morphogenesis.

Fluctuations of cell geometry and their nonequilibrium thermodynamics in living epithelial tissue

We measure different contributions to entropy production in a living functional epithelial tissue. We do this by extracting the functional dynamics of development while at the same time quantifying fluctuations. Using the translucent Drosophila melanogaster pupal epithelium as an ideal tissue for high-resolution live imaging, we measure the entropy associated with the stochastic geometry of cells in the epithelium. This is done using a detailed analysis of the dynamics of the shape and orientation of individual cells which enables separation of local and global aspects of the tissue behavior. Intriguingly, we find that we can observe irreversible dynamics in the cell geometries but without a change in the entropy associated with those degrees of freedom, showing that there is a flow of energy into those degrees of freedom. Hence, the living system is controlling how the entropy is being produced and partitioned into its different parts.

Cell polarity and extrusion: How to polarize extrusion and extrude misspolarized cells?

The barrier function of epithelia is one of the cornerstones of the body plan organization of metazoans. It relies on the polarity of epithelial cells which organizes along the apico-basal axis the mechanical properties, signaling as well as transport. This barrier function is however constantly challenged by the fast turnover of epithelia occurring during morphogenesis or adult tissue homeostasis. Yet, the sealing property of the tissue can be maintained thanks to cell extrusion: a series of remodeling steps involving the dying cell and its neighbors leading to seamless cell expulsion. Alternatively, the tissue architecture can also be challenged by local damages or the emergence of mutant cells that may alter its organization. This includes mutants of the polarity complexes which can generate neoplastic overgrowths or be eliminated by cell competition when surrounded by wild type cells. In this review, we will provide an overview of the regulation of cell extrusion in various tissues focusing on the relationship between cell polarity, cell organization and the direction of cell expulsion. We will then describe how local perturbations of polarity can also trigger cell elimination either by apoptosis or by cell exclusion, focusing specifically on how polarity defects can be directly causal to cell elimination. Overall, we propose a general framework connecting the influence of polarity on cell extrusion and its contribution to aberrant cell elimination.

Tilt-induced polar order and topological defects in growing bacterial populations

Rod-shaped bacteria, such as Escherichia coli, commonly live forming mounded colonies. They initially grow two-dimensionally on a surface and finally achieve three-dimensional growth. While it was recently reported that three-dimensional growth is promoted by topological defects of winding number +1/2 in populations of motile bacteria, how cellular alignment plays a role in nonmotile cases is largely unknown. Here, we investigate the relevance of topological defects in colony formation processes of nonmotile E. coli populations, and found that both ±1/2 topological defects contribute to the three-dimensional growth. Analyzing the cell flow in the bottom layer of the colony, we observe that +1/2 defects attract cells and -1/2 defects repel cells, in agreement with previous studies on motile cells, in the initial stage of the colony growth. However, later, cells gradually flow toward -1/2 defects as well, exhibiting a sharp contrast to the existing knowledge. By investigating three-dimensional cell orientations by confocal microscopy, we find that vertical tilting of cells is promoted near the defects. Crucially, this leads to the emergence of a polar order in the otherwise nematic two-dimensional cell orientation. We extend the theory of active nematics by incorporating this polar order and the vertical tilting, which successfully explains the influx toward -1/2 defects in terms of a polarity-induced force. Our work reveals that three-dimensional cell orientations may result in qualitative changes in properties of active nematics, especially those of topological defects, which may be generically relevant in active matter systems driven by cellular growth instead of self-propulsion.

Topological defect-mediated morphodynamics of active-active interfaces

Physical interfaces widely exist in nature and engineering. Although the formation of passive interfaces is well elucidated, the physical principles governing active interfaces remain largely unknown. Here, we combine simulation, theory, and cell-based experiment to investigate the evolution of an active-active interface. We adopt a biphasic framework of active nematic liquid crystals. We find that long-lived topological defects mechanically energized by activity display unanticipated dynamics nearby the interface, where defects perform "U-turns" to keep away from the interface, push the interface to develop local fingers, or penetrate the interface to enter the opposite phase, driving interfacial morphogenesis and cross-interface defect transport. We identify that the emergent interfacial morphodynamics stems from the instability of the interface and is further driven by the activity-dependent defect-interface interactions. Experiments of interacting multicellular monolayers with extensile and contractile differences in cell activity have confirmed our predictions. These findings reveal a crucial role of topological defects in active-active interfaces during, for example, boundary formation and tissue competition that underlie organogenesis and clinically relevant disorders.

Organoids

Organoids have attracted increasing attention because they are simple tissue-engineered cell-based in vitro models that recapitulate many aspects of the complex structure and function of the corresponding in vivo tissue. They can be dissected and interrogated for fundamental mechanistic studies on development, regeneration, and repair in human tissues. Organoids can also be used in diagnostics, disease modeling, drug discovery, and personalized medicine. Organoids are derived from either pluripotent or tissue-resident stem (embryonic or adult) or progenitor or differentiated cells from healthy or diseased tissues, such as tumors. To date, numerous organoid engineering strategies that support organoid culture and growth, proliferation, differentiation and maturation have been reported. This Primer serves to highlight the rationale underlying the selection and development of these materials and methods to control the cellular/tissue niche; and therefore, structure and function of the engineered organoid. We also discuss key considerations for generating robust organoids, such as those related to cell isolation and seeding, matrix and soluble factor selection, physical cues and integration. The general standards for data quality, reproducibility and deposition within the organoid community is also outlined. Lastly, we conclude by elaborating on the limitations of organoids in different applications, and key priorities in organoid engineering for the coming years.

Mechanical regulation of the early stages of angiogenesis

Favouring or thwarting the development of a vascular network is essential in fields as diverse as oncology, cardiovascular disease or tissue engineering. As a result, understanding and controlling angiogenesis has become a major scientific challenge. Mechanical factors play a fundamental role in angiogenesis and can potentially be exploited for optimizing the architecture of the resulting vascular network. Largely focusing on in vitro systems but also supported by some in vivo evidence, the aim of this Highlight Review is dual. First, we describe the current knowledge with particular focus on the effects of fluid and solid mechanical stimuli on the early stages of the angiogenic process, most notably the destabilization of existing vessels and the initiation and elongation of new vessels. Second, we explore inherent difficulties in the field and propose future perspectives on the use of in vitro and physics-based modelling to overcome these difficulties.

Spontaneous multi-scale void formation and closure in densifying epithelial and fibroblast monolayers from the sub-confluent state

Using time-lapse phase contrast microscopy, the formation and closure of spontaneously generated voids in the densifying monolayers of isotropic epithelial cells (ECs) and elongated fibroblast cells (FCs) through proliferation from the sub-confluent state are investigated. It is found that, in both types of monolayers after forming a connected network composed of nematic patches with different orientations, numerous multi-scale voids can be spontaneously formed and gradually close with increasing time. The isotropic fluctuations of deformation and crawling of ECs and the anisotropic axial motion/alignment polarizations of FCs are the two keys leading to the following different generic dynamical behaviors. In EC monolayers, voids exhibit irregular boundary fluctuations and easier cell re-orientation of front layer cells (FLCs) surrounding void boundaries. Void closures are mainly through pinching the gap between the opposite fluctuating void boundaries, and the inward crawling of FLCs to reduce void area associated with topological rearrangement to reduce FLC number. In FC monolayers, large voids have piecewise smooth convex boundaries, and cusp-shaped concave boundaries with cells orienting toward the void at cusp tips. The extension of a thin cell bridge from the cusp tip can bisect a large void into smaller voids. For smaller FC voids dominated by convex boundaries, along which cell alignment prohibits inward crawling, the reduction of FLC number through successive outward squeezing of single FLCs by neighboring FLCs sliding along the void boundary plays an important role for topological rearrangement and void closure. Unlike those surrounding artificial wounds in dense EC monolayers, the absence of ring-like purse-strings surrounding EC and FC voids allows topological rearrangements for reducing void perimeter and void area.

Defect dynamics in active polar fluids vs. active nematics

Topological defects play a key role in two-dimensional active nematics, and a transient role in two-dimensional active polar fluids. Using a variational method, we study both the transient and long-time behavior of defects in two-dimensional active polar fluids in the limit of strong order and overdamped, compressible flow, and compare the defect dynamics with the corresponding active nematics model studied recently. One result is non-central interactions between defect pairs for active polar fluids, and by extending our analysis to allow orientation dynamics of defects, we find that the orientation of +1 defects, unlike that of ±1/2 defects in active nematics, is not locked to defect positions and relaxes to asters. Moreover, using a scaling argument, we explain the transient feature of active polar defects and show that in the steady state, active polar fluids are either devoid of defects or consist of a single aster. We argue that for contractile (extensile) active nematic systems, +1 vortices (asters) should emerge as bound states of a pair of +1/2 defects, which has been recently observed. Moreover, unlike the polar case, we show that for active nematics, a linear chain of equally spaced bound states of pairs of +1/2 defects can screen the activity term. A common feature in both models is the appearance of +1 defects (elementary in polar and composite in nematic) in the steady state.

Active nematic defects in compressible and incompressible flows

We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint. Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, while they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a smecticlike pattern, which induce an associated global polar ordering of +1/2 nematic defects. By contrast, divergence-free flows give rise to a local nematic order of the +1/2 defects, consisting of antialigned pairs of neighboring defects, as established in previous studies.

The Nematic Chiral Liquid Crystal Structure of the Cardiac Myoarchitecture: Disclinations and Topological Singularities

This is our second article devoted to the cardiac myoarchitecture considered as a nematic chiral liquid crystal (NCLC). While the first article focused on the myoarchitecture of the left ventricle (LV), this new article extends to the whole ventricular mass and introduces the concept of disclinations and topological singularities, which characterize the differences and relationships between the left and right ventricles (RV). At the level of the ventricular apices, we constantly observed a vortex shape at the LV apex, corresponding, in the terminology of liquid crystals, to a "+1 disclination"; we never observed this at the RV apex. At the level of the interventricular septum (IVS), we identified "-1/2 disclinations" at the anterior and posterior parts. During the perinatal period, there was a significant difference in their distribution, with more "-1/2 disclinations" in the posterior part of the IVS. After birth, concomitant to major physiological changes, the number of "-1/2 disclinations" significantly decreased, both in the anterior and posterior parts of the IVS. Finally, the description of the disclinations must be considered in any attempt to segment the whole ventricular mass, in biomechanical studies, and, more generally, for the characterization of myocardial remodeling.

Scaling and spontaneous symmetry restoring of topological defect dynamics in liquid crystal

Topological defects-locations of local mismatch of order-are a universal concept playing important roles in diverse systems studied in physics and beyond, including the universe, various condensed matter systems, and recently, even life phenomena. Among these, liquid crystal has been a platform for studying topological defects via visualization, yet it has been a challenge to resolve three-dimensional structures of dynamically evolving singular topological defects. Here, we report a direct confocal observation of nematic liquid crystalline defect lines, called disclinations, relaxing from an electrically driven turbulent state. We focus in particular on reconnections, characteristic of such line defects. We find a scaling law for in-plane reconnection events, by which the distance between reconnecting disclinations decreases by the square root of time to the reconnection. Moreover, we show that apparently asymmetric dynamics of reconnecting disclinations is actually symmetric in a comoving frame, in marked contrast to the two-dimensional counterpart whose asymmetry is established. We argue, with experimental supports, that this is because of energetically favorable symmetric twist configurations that disclinations take spontaneously, thanks to the topology that allows for rotation of the winding axis. Our work illustrates a general mechanism of such spontaneous symmetry restoring that may apply beyond liquid crystal, which can take place if topologically distinct asymmetric defects in lower dimensions become homeomorphic in higher dimensions and if the symmetric intermediate is energetically favorable.

Fluctuation-induced dynamics of nematic topological defects

Topological defects are increasingly being identified in various biological systems, where their characteristic flow fields and stress patterns are associated with continuous active stress generation by biological entities. Here, using numerical simulations of continuum fluctuating nematohydrodynamics, we show that even in the absence of any specific form of active stresses associated with self-propulsion, mesoscopic fluctuations in either orientational alignment or hydrodynamics can independently result in flow patterns around topological defects that resemble the ones observed in active systems. Our simulations further show the possibility of extensile- and contractile-like motion of fluctuation-induced positive half-integer topological defects. Remarkably, isotropic stress fields also reproduce the experimentally measured stress patterns around topological defects in epithelia. Our findings further reveal that extensile- or contractile-like flow and stress patterns around fluctuation-induced defects are governed by passive elastic stresses and flow-aligning behavior of the nematics.

Local Yield and Compliance in Active Cell Monolayers

The rheology of biological tissue plays an important role in many processes, from organ formation to cancer invasion. Here, we use a multiphase field model of motile cells to simulate active microrheology within a tissue monolayer. When unperturbed, the tissue exhibits a transition between a solidlike state and a fluidlike state tuned by cell motility and deformability-the ratio of the energetic costs of steric cell-cell repulsion and cell-edge tension. When perturbed, solid tissues exhibit local yield-stress behavior, with a threshold force for the onset of motion of a probe particle that vanishes upon approaching the solid-to-liquid transition. This onset of motion is qualitatively different in the low and high deformability regimes. At high deformability, the tissue is amorphous when solid, it responds compliantly to deformations, and the probe transition to motion is smooth. At low deformability, the monolayer is more ordered translationally and stiffer, and the onset of motion appears discontinuous. Our results suggest that cellular or nanoparticle transport in different types of tissues can be fundamentally different and point to ways in which it can be controlled.

Dominant geometrical factors of collective cell migration in flexible 3D gelatin tube structures

Collective cell migration is a dynamic and interactive behavior of cell cohorts essential for diverse physiological developments in living organisms. Recent studies have revealed the importance of three-dimensional (3D) topographical confinements to regulate the migration modes of cell cohorts in tubular confinement. However, conventional in vitro assays fail to observe cells' behavior in response to 3D structural changes, which is necessary for examining the geometric regulation factors of collective migration. Here, we introduce a newly developed assay for fabricating flexible 3D structures of capillary microtunnels to examine the behavior of vascular endothelial cells (ECs) as they progress through the successive transition across wide or narrow tube structures. The microtunnels with altered diameters were formed inside gelatin-gel blocks by photo-thermal etching with micrometer-sized spot heating of the focused infrared laser absorption. The ECs migrated and spread two-dimensionally on the inner surface of gelatin capillary microtunnels as a monolayer instead of filling the entire capillary. In the straight cylindrical topographical constraint, leading ECs exhibited no apparent diameter dependence for the maximum peak migration velocity. However, widening the diameter in the narrow-wide structures caused a decrease in migration velocity following in direct proportion to the diameter increase ratio, whereas narrowing the diameter in wide-narrow microtunnels increased the speed without obvious correlation between velocity change and diameter change. The results demonstrated the ability of the newly developed flexible 3D gelatin tube structures for collective cell migration, and the findings provide insights into the dominant geometric factor of the emerging migratory modes for endothelial migration as asymmetric fluid flow-like behavior in the borderless cylindrical cell sheets.

Active Nematic Flows over Curved Surfaces

Cell monolayers are a central model system in the study of tissue biophysics. In vivo, epithelial tissues are curved on the scale of microns, and the curvature's role in the onset of spontaneous tissue flows is still not well understood. Here, we present a hydrodynamic theory for an apical-basal asymmetric active nematic gel on a curved strip. We show that surface curvature qualitatively changes monolayer motion compared with flat space: the resulting flows can be thresholdless, and the transition to motion may change from continuous to discontinuous. Surface curvature, friction, and active tractions are all shown to control the flow pattern selected, from simple shear to vortex chains.

Active Nematic Defects and Epithelial Morphogenesis

Inspired by recent experiments that highlight the role of nematic defects in the morphogenesis of epithelial tissues, we develop a minimal framework to study the dynamics of an active curved surface driven by its nematic texture. Allowing the surface to evolve via relaxational dynamics leads to a theory linking nematic defect dynamics, cellular division rates, and Gaussian curvature. Regions of large positive (negative) curvature and positive (negative) growth are colocalized with the presence of positive (negative) defects. In an ex-vivo setting of cultured murine neural progenitor cells, we show that our framework is consistent with the observed cell accumulation at positive defects and depletion at negative defects. In an in-vivo setting, we show that the defect configuration consisting of a bound +1 defect state, which is stabilized by activity, surrounded by two -1/2 defects can create a stationary ring configuration of tentacles, consistent with observations of a basal marine invertebrate Hydra.

Activity gradients in two- and three-dimensional active nematics

We numerically investigate how spatial variations of extensile or contractile active stress affect bulk active nematic systems in two and three dimensions. In the absence of defects, activity gradients drive flows which re-orient the nematic director field and thus act as an effective anchoring force. At high activity, defects are created and the system transitions into active turbulence, a chaotic flow state characterized by strong vorticity. We find that in two-dimensional (2D) systems active torques robustly align +1/2 defects parallel to activity gradients, with defect heads pointing towards contractile regions. In three-dimensional (3D) active nematics disclination lines preferentially lie in the plane perpendicular to activity gradients due to active torques acting on line segments. The average orientation of the defect structures in the plane perpendicular to the line tangent depends on the defect type, where wedge-like +1/2 defects align parallel to activity gradients, while twist defects are aligned anti-parallel. Understanding the response of active nematic fluids to activity gradients is an important step towards applying physical theories to biology, where spatial variations of active stress impact morphogenetic processes in developing embryos and affect flows and deformations in growing cell aggregates, such as tumours.

Plasticity of body axis polarity in Hydra regeneration under constraints

One of the major events in animal morphogenesis is the emergence of a polar body axis. Here, we combine classic grafting techniques with live imaging to explore the plasticity of polarity determination during whole body regeneration in Hydra. Composite tissues are made by fusing two rings, excised from separate animals, in different configurations that vary in the polarity and original positions of the rings along the body axes of the parent animals. Under frustrating initial configurations, body axis polarity that is otherwise stably inherited from the parent animal, can become labile and even be reversed. Importantly, the site of head regeneration exhibits a strong bias toward the edges of the tissue, even when this involves polarity reversal. In particular, we observe head formation at an originally aboral tissue edge, which is not compatible with models of Hydra regeneration based only on preexisting morphogen gradients or an injury response. The site of the new head invariably contains an aster-like defect in the organization of the supra-cellular ectodermal actin fibers. While a defect is neither required nor sufficient for head formation, we show that the defect at the new head site can arise via different routes, either appearing directly following excision as the tissue seals at its edge or through de novo defect formation at the fusion site. Altogether, our results show that the emergence of a polar body axis depends on the original polarity and position of the excised tissues as well as structural factors, suggesting that axis determination is an integrated process that arises from the dynamic interplay of multiple biochemical and mechanical processes.

Defect absorption and emission for p-atic liquid crystals on cones

We investigate the ground-state configurations of two-dimensional liquid crystals with p-fold rotational symmetry (p-atics) on fixed curved surfaces. We focus on the intrinsic geometry and show that isothermal coordinates are particularly convenient as they explicitly encode a geometric contribution to the elastic potential. In the special case of a cone with half-angle β, the apex develops an effective topological charge of -χ, where 2πχ=2π(1-sinβ) is the deficit angle of the cone, and a topological defect of charge σ behaves as if it had an effective topological charge Q_{eff}=(σ-σ^{2}/2) when interacting with the apex. The effective charge of the apex leads to defect absorption and emission at the cone apex as the deficit angle of the cone is varied. For total topological defect charge 1, e.g., imposed by tangential boundary conditions at the edge, we find that for a disk the ground-state configuration consists of p defects each of charge +1/p lying equally spaced on a concentric ring of radius d=(p-1/3p-1)^{1/2p}R, where R is the radius of the disk. In the case of a cone with tangential boundary conditions at the base, we find three types of ground-state configurations as a function of cone angle: (i) for sharp cones, all of the +1/p defects are absorbed by the apex; (ii) at intermediate cone angles, some of the +1/p defects are absorbed by the apex and the rest lie equally spaced along a concentric ring on the flank; and (iii) for nearly flat cones, all of the +1/p defects lie equally spaced along a concentric ring on the flank. Here the defect positions and the absorption transitions depend intricately on p and the deficit angle, which we analytically compute. We check these results with numerical simulations for a set of commensurate cone angles and find excellent agreement.

Impact of dipole-dipole interactions on motility-induced phase separation

We present a hydrodynamic theory for systems of dipolar active Brownian particles which, in the regime of weak dipolar coupling, predicts the onset of motility-induced phase separation (MIPS), consistent with Brownian dynamics (BD) simulations. The hydrodynamic equations are derived by explicitly coarse-graining the microscopic Langevin dynamics, thus allowing for a mapping of the coarse-grained model and particle-resolved simulations. Performing BD simulations at fixed density, we find that dipolar interactions tend to hinder MIPS, as first reported in [Liao et al., Soft Matter, 2020, 16, 2208]. Here we demonstrate that the theoretical approach indeed captures the suppression of MIPS. Moreover, the analysis of the numerically obtained, angle-dependent correlation functions sheds light into the underlying microscopic mechanisms leading to the destabilization of the homogeneous phase.

Coupling the topological defect phase to the extrinsic curvature in nematic shells

In two dimensional nematics, topological defects are point like singularities with both a charge and a phase. We study the topological defects within curved nematic textures on the surface of a cylinder. This allows us to isolate the effect of extrinsic curvature on the structure of the topological defect. By minimizing the energy associated with distortions in the nematic director around the core of a defect, we show that the phase of the topological defect is coupled to the orientation of the cylinder. This coupling depends on the relative energetic cost associated with splay, bend and twist distortions of the nematic director. We identify a bistability in the phase of the defects when twist deformations dominate. Finally, we show a similar effect for integer charge topological defects.

Nematic order condensation and topological defects in inertial active nematics

Living materials at different length scales manifest active nematic features such as orientational order, nematic topological defects, and active nematic turbulence. Using numerical simulations we investigate the impact of fluid inertia on the collective pattern formation in active nematics. We show that an incremental increase in inertial effects due to reduced viscosity results in gradual melting of nematic order with an increase in topological defect density before a discontinuous transition to a vortex-condensate state. The emergent vortex-condensate state at low enough viscosities coincides with nematic order condensation within the giant vortices and the drop in the density of topological defects. We further show flow field around topological defects is substantially affected by inertial effects. Moreover, we demonstrate the strong dependence of the kinetic energy spectrum on the inertial effects, recover the Kolmogorov scaling within the vortex-condensate phase, but find no evidence of universal scaling at higher viscosities. The findings reveal complexities in active nematic turbulence and emphasize the important cross-talk between active and inertial effects in setting flow and orientational organization of active particles.

Bridging microscopic cell dynamics to nematohydrodynamics of cell monolayers

It is increasingly being realized that liquid-crystalline features can play an important role in the properties and dynamics of cell monolayers. Here, we present a cell-based model of cell layers, based on the phase-field formulation, that connects cell-cell interactions specified at the single cell level to large-scale nematic and hydrodynamic properties of the tissue. In particular, we present a minimal formulation that reproduces the well-known bend and splay hydrodynamic instabilities of the continuum nemato-hydrodynamic formulation of active matter, together with an analytical description of the instability threshold in terms of activity and elasticity of the cells. Furthermore, we provide a quantitative characterisation and comparison of flows and topological defects for extensile and contractile stress generation mechanisms, and demonstrate activity-induced heterogeneity and spontaneous formation of gaps within a confluent monolayer. Together, these results contribute to bridging the gap between cell-scale dynamics and tissue-scale collective cellular organisation.

Mechanical Forces Govern Interactions of Host Cells with Intracellular Bacterial Pathogens

To combat infectious diseases, it is important to understand how host cells interact with bacterial pathogens. Signals conveyed from pathogen to host, and vice versa, may be either chemical or mechanical. While the molecular and biochemical basis of host-pathogen interactions has been extensively explored, relatively less is known about mechanical signals and responses in the context of those interactions. Nevertheless, a wide variety of bacterial pathogens appear to have developed mechanisms to alter the cellular biomechanics of their hosts in order to promote their survival and dissemination, and in turn many host responses to infection rely on mechanical alterations in host cells and tissues to limit the spread of infection. In this review, we present recent findings on how mechanical forces generated by host cells can promote or obstruct the dissemination of intracellular bacterial pathogens. In addition, we discuss how in vivo extracellular mechanical signals influence interactions between host cells and intracellular bacterial pathogens. Examples of such signals include shear stresses caused by fluid flow over the surface of cells and variable stiffness of the extracellular matrix on which cells are anchored. We highlight bioengineering-inspired tools and techniques that can be used to measure host cell mechanics during infection. These allow for the interrogation of how mechanical signals can modulate infection alongside biochemical signals. We hope that this review will inspire the microbiology community to embrace those tools in future studies so that host cell biomechanics can be more readily explored in the context of infection studies.

Thickness of epithelia on wavy substrates: measurements and continuous models

We measured the thickness of MDCK epithelia grown on substrates with a sinusoidal profile. We show that while at long wavelength the profile of the epithelium follows that of the substrate, at short wavelengths cells are thicker in valleys than on ridges. This is reminiscent of the so-called «healing length in the case of a thin liquid film wetting a rough solid substrate. We explore the ability of continuum mechanics models to account for these observations. Modeling the epithelium as a thin liquid film, with surface tension, does not fully account for the measurements. Neither does modeling the epithelium as a thin incompressible elastic film. On the contrary, the addition of an apical active stress gives satisfactory agreement with measurements, with one fitting parameter, the ratio between the active stress and the elastic modulus.

Defect loops in three-dimensional active nematics as active multipoles

We develop a description of defect loops in three-dimensional active nematics based on a multipole expansion of the far-field director and show how this leads to a self-dynamics dependent on the loop's geometric type. The dipole term leads to active stresses that generate a global self-propulsion for splay and bend loops. The quadrupole moment is nonzero only for nonplanar loops and generates a net "active torque," such that defect loops are both self-motile and self-orienting. Our analysis identifies right- and left-handed twist loops as the only force- and torque-free geometries, suggesting a mechanism for generating an excess of twist loops. Finally, we determine the Stokesian flows created by defect loops and describe qualitatively their hydrodynamics.

Quadrupolar active stress induces exotic patterns of defect motion in compressible active nematics

A wide range of living and artificial active matter exists in close contact with substrates and under strong confinement, where in addition to dipolar active stresses, quadrupolar active stresses can become important. Here, we numerically investigate the impact of quadrupolar non-equilibrium stresses on the emergent patterns of self-organisation in non-momentum conserving active nematics. Our results reveal that beyond having stabilising effects, the quadrupolar active forces can induce various modes of topological defect motion in active nematics. In particular, we find the emergence of both polar and nematic ordering of the defects, as well as new patterns of self-organisation that comprise topological defect chains and transient topological defect asters. The results contribute to further understanding of emergent patterns of collective motion and non-equilibrium self-organisation in active matter.

Self-rotation regulates interface evolution in biphasic active matter through taming defect dynamics

Chirality can endow nonequilibrium active matter with unique features and functions. Here, we explore the chiral dynamics in biphasic active nematics composed of self-rotating units that continuously inject energy and angular momentum at the microscale. We show that the self-rotation of units can regularize the boundaries between two phases, rendering sinusoidal-like interfaces, which allow lateral wave propagation and are characterized by chains of ordered antiferromagnetic cross-interface flow vortices. Through the spontaneous coordination of counter-rotating units across the interfaces, topological defects excited by activity are sorted spatiotemporally, where positive defects are locally trapped at the interfaces but, unexpectedly, are transported laterally in a unidirectional rather than wavy mode, whereas inertial negative defects remain spinning in the bulks. Our findings reveal that individual chirality could be harnessed to modulate interfacial morphodynamics in active systems and suggest a potential approach toward controlling topological defects for programmable microfluidics and logic operations.

Integer topological defects organize stresses driving tissue morphogenesis

Tissues acquire function and shape via differentiation and morphogenesis. Both processes are driven by coordinating cellular forces and shapes at the tissue scale, but general principles governing this interplay remain to be discovered. Here we report that self-organization of myoblasts around integer topological defects, namely spirals and asters, suffices to establish complex multicellular architectures. In particular, these arrangements can trigger localized cell differentiation or, alternatively, when differentiation is inhibited, they can drive the growth of swirling protrusions. Both localized differentiation and growth of cellular vortices require specific stress patterns. By analysing the experimental velocity and orientational fields through active gel theory, we show that integer topological defects can generate force gradients that concentrate compressive stresses. We reveal these gradients by assessing spatial changes in nuclear volume and deformations of elastic pillars. We propose integer topological defects as mechanical organizing centres controlling differentiation and morphogenesis.

Bistability of Dielectrically Anisotropic Nematic Crystals and the Adaptation of Endothelial Collectives to Stress Fields

Endothelial monolayers physiologically adapt to flow and flow-induced wall shear stress, attaining ordered configurations in which elongation, orientation, and polarization are coherently organized over many cells. Here, with the flow direction unchanged, a peculiar bi-stable (along the flow direction or perpendicular to it) cell alignment is observed, emerging as a function of the flow intensity alone, while cell polarization is purely instructed by flow directionality. Driven by the experimental findings, the parallelism between endothelia is delineated under a flow field and the transition of dual-frequency nematic liquid crystals under an external oscillatory electric field. The resulting physical model reproduces the two stable configurations and the energy landscape of the corresponding system transitions. In addition, it reveals the existence of a disordered, metastable state emerging upon system perturbation. This intermediate state, experimentally demonstrated in endothelial monolayers, is shown to expose the cellular system to a weakening of cell-to-cell junctions to the detriment of the monolayer integrity. The flow-adaptation of monolayers composed of healthy and senescent endothelia is successfully predicted by the model with adjustable nematic parameters. These results may help to understand the maladaptive response of in vivo endothelial tissues to disturbed hemodynamics and the progressive functional decay of senescent endothelia.