Mechanical forces across compartments coordinate cell shape and fate transitions to generate tissue architecture

Morphogenesis and cell state transitions must be coordinated in time and space to produce a functional tissue. An excellent paradigm to understand the coupling of these processes is mammalian hair follicle development, which is initiated by the formation of an epithelial invagination-termed placode-that coincides with the emergence of a designated hair follicle stem cell population. The mechanisms directing the deformation of the epithelium, cell state transitions and physical compartmentalization of the placode are unknown. Here we identify a key role for coordinated mechanical forces stemming from contractile, proliferative and proteolytic activities across the epithelial and mesenchymal compartments in generating the placode structure. A ring of fibroblast cells gradually wraps around the placode cells to generate centripetal contractile forces, which, in collaboration with polarized epithelial myosin activity, promote elongation and local tissue thickening. These mechanical stresses further enhance compartmentalization of Sox9 expression to promote stem cell positioning. Subsequently, proteolytic remodelling locally softens the basement membrane to facilitate a release of pressure on the placode, enabling localized cell divisions, tissue fluidification and epithelial invagination into the underlying mesenchyme. Together, our experiments and modelling identify dynamic cell shape transformations and tissue-scale mechanical cooperation as key factors for orchestrating organ formation.

Differences in boundary behavior in the 3D vertex and Voronoi models

An important open question in the modeling of biological tissues is how to identify the right scale for coarse-graining, or equivalently, the right number of degrees of freedom. For confluent biological tissues, both vertex and Voronoi models, which differ only in their representation of the degrees of freedom, have effectively been used to predict behavior, including fluid-solid transitions and cell tissue compartmentalization, which are important for biological function. However, recent work in 2D has hinted that there may be differences between the two models in systems with heterotypic interfaces between two tissue types, and there is a burgeoning interest in 3D tissue models. Therefore, we compare the geometric structure and dynamic sorting behavior in mixtures of two cell types in both 3D vertex and Voronoi models. We find that while the cell shape indices exhibit similar trends in both models, the registration between cell centers and cell orientation at the boundary are significantly different between the two models. We demonstrate that these macroscopic differences are caused by changes to the cusp-like restoring forces introduced by the different representations of the degrees of freedom at the boundary, and that the Voronoi model is more strongly constrained by forces that are an artifact of the way the degrees of freedom are represented. This suggests that vertex models may be more appropriate for 3D simulations of tissues with heterotypic contacts.

Differences in boundary behavior in the 3D vertex and Voronoi models

An important open question in the modeling of biological tissues is how to identify the right scale for coarse-graining, or equivalently, the right number of degrees of freedom. For confluent biological tissues, both vertex and Voronoi models, which differ only in their representation of the degrees of freedom, have effectively been used to predict behavior, including fluid-solid transitions and cell tissue compartmentalization, which are important for biological function. However, recent work in 2D has hinted that there may be differences between the two models in systems with heterotypic interfaces between two tissue types, and there is a burgeoning interest in 3D tissue models. Therefore, we compare the geometric structure and dynamic sorting behavior in mixtures of two cell types in both 3D vertex and Voronoi models. We find that while the cell shape indices exhibit similar trends in both models, the registration between cell centers and cell orientation at the boundary are significantly different between the two models. We demonstrate that these macroscopic differences are caused by changes to the cusp-like restoring forces introduced by the different representations of the degrees of freedom at the boundary, and that the Voronoi model is more strongly constrained by forces that are an artifact of the way the degrees of freedom are represented. This suggests that vertex models may be more appropriate for 3D simulations of tissues with heterotypic contacts.

Collective chemotaxis in a Voronoi model for confluent clusters

Collective chemotaxis, where single cells cannot climb a biochemical signaling gradient but clusters of cells can, has been observed in different biological contexts, including confluent tissues where there are no gaps or overlaps between cells. Although particle-based models have been developed that predict important features of collective chemotaxis, the mechanisms in those models depend on particle overlaps, and so it remains unclear if they can explain behavior in confluent systems. Here, we develop an open-source code that couples a two-dimensional Voronoi simulation for confluent cell mechanics to a dynamic chemical signal that can diffuse, advect, and/or degrade and use the code to study potential mechanisms for collective chemotaxis in cellular monolayers. We first study the impact of advection on collective chemotaxis and delineate a regime where advective terms are important. Next, we investigate two possible chemotactic mechanisms, contact inhibition of locomotion and heterotypic interfacial tension, and demonstrate that both can drive collective chemotaxis in certain parameter regimes. We further demonstrate that the scaling behavior of cluster motion is well captured by simple analytic theories.

Jamming and arrest of cell motion in biological tissues

Collective cell motility is crucial to many biological processes including morphogenesis, wound healing, and cancer invasion. Recently, the biology and biophysics communities have begun to use the term 'cell jamming' to describe the collective arrest of cell motion in tissues. Although this term is widely used, the underlying mechanisms are varied. In this review, we highlight three independent mechanisms that can potentially drive arrest of cell motion - crowding, tension-driven rigidity, and reduction of fluctuations - and propose a framework that connects all three. Because multiple mechanisms may be operating simultaneously, this emphasizes that experiments should strive to identify which mechanism dominates in a given situation. We also discuss how specific cell-scale and molecular-scale biological processes, such as cell-cell and cell-substrate interactions, control aspects of these underlying physical mechanisms.