Mechanical basis and topological routes to cell elimination

Force transmission is a master regulator of mechanical cell competition

Cell competition is a tissue surveillance mechanism for eliminating unwanted cells, being indispensable in development, infection and tumourigenesis. Although studies have established the role of biochemical mechanisms in this process, due to challenges in measuring forces in these systems, how mechanical forces determine the competition outcome remains unclear. Here we report a form of cell competition that is regulated by differences in force transmission capabilities, selecting for cell types with stronger intercellular adhesion. Direct force measurements in ex vivo tissues and different cell lines reveal that there is an increased mechanical activity at the interface between two competing cell types, which can lead to large stress fluctuations resulting in upward forces and cell elimination. We show how a winning cell type endowed with a stronger intercellular adhesion exhibits higher resistance to elimination and benefiting from efficient force transmission to the neighbouring cells. This cell elimination mechanism could have broad implications for keeping the strong force transmission ability for maintaining tissue boundaries and cell invasion pathology.

Cell-Level Modelling of Homeostasis in Confined Epithelial Monolayers

Tissue homeostasis, the biological process of maintaining a steady state in tissue via control of cell proliferation and death, is essential for the development, growth, maintenance, and proper function of living organisms. Disruptions to this process can lead to serious diseases and even death. In this study, we use the vertex model for the cell-level description of tissue mechanics to investigate the impact of the tissue environment and local mechanical properties of cells on homeostasis in confined epithelial tissues. We find a dynamic steady state, where the balance between cell divisions and removals sustains homeostasis, and characterise the homeostatic state in terms of cell count, tissue area, homeostatic pressure, and the cells' neighbour count distribution. This work, therefore, sheds light on the mechanisms underlying tissue homeostasis and highlights the importance of mechanics in its control.

Assessing mechanical agency during apical apoptotic cell extrusion

Homeostasis is necessary for epithelia to maintain barrier function and prevent the accumulation of defective cells. Unfit, excess, and dying cells in the larval zebrafish tail fin epidermis are removed via controlled cell death and extrusion. Extrusion coincides with oscillations of cell area, both in the extruding cell and its neighbors. Here, we develop a biophysical model of this process to explore the role of autonomous and non-autonomous mechanics. We vary biophysical properties and oscillatory behaviors of extruding cells and their neighbors along with tissue-wide cell density and viscosity. We find that cell autonomous processes are major contributors to the dynamics of extrusion, with the mechanical microenvironment providing a less pronounced contribution. We also find that some cells initially resist extrusion, influencing the duration of the expulsion process. Our model provides insights into the cellular dynamics and mechanics that promote elimination of unwanted cells from epithelia during homeostatic tissue maintenance.

The cellular Potts model on disordered lattices

The cellular Potts model, also known as the Glazier-Graner-Hogeweg model, is a lattice-based approach by which biological tissues at the level of individual cells can be numerically studied. Traditionally, a square or hexagonal underlying lattice structure is assumed for two-dimensional systems, and this is known to introduce artifacts in the structure and dynamics of the model tissues. That is, on regular lattices, cells can assume shapes that are dictated by the symmetries of the underlying lattice. Here, we developed a variant of this method that can be applied to a broad class of (ir)regular lattices. We show that on an irregular lattice deriving from a fluid-like configuration, two types of artifacts can be removed. We further report on the transition between a fluid-like disordered and a solid-like hexagonally ordered phase present for monodisperse confluent cells as a function of their surface tension. This transition shows the hallmarks of a first-order phase transition and is different from the glass/jamming transitions commonly reported for the vertex and active Voronoi models. We emphasize this by analyzing the distribution of shape parameters found in our state space. Our analysis provides a useful reference for the future study of epithelia using the (ir)regular cellular Potts model.

Intercellular friction and motility drive orientational order in cell monolayers

Spatiotemporal patterns in multicellular systems are important to understanding tissue dynamics, for instance, during embryonic development and disease. Here, we use a multiphase field model to study numerically the behavior of a near-confluent monolayer of deformable cells with intercellular friction. Varying friction and cell motility drives a solid-liquid transition, and near the transition boundary, we find the emergence of local nematic order of cell deformation driven by shear-aligning cellular flows. Intercellular friction contributes to the monolayer's viscosity, which significantly increases the spatial correlation in the flow and, concomitantly, the extent of nematic order. We also show that local hexatic and nematic order are tightly coupled and propose a mechanical-geometric model for the colocalization of [Formula: see text] nematic defects and 5-7 disclination pairs, which are the structural defects in the hexatic phase. Such topological defects coincide with regions of high cell-cell overlap, suggesting that they may mediate cellular extrusion from the monolayer, as found experimentally. Our results delineate a mechanical basis for the recent observation of nematic and hexatic order in multicellular collectives in experiments and simulations and pinpoint a generic pathway to couple topological and physical effects in these systems.

Multiphase field model of cells on a substrate: From three dimensional to two dimensional

Multiphase field models have emerged as an important computational tool for understanding biological tissue while resolving single-cell properties. While they have successfully reproduced many experimentally observed behaviors of living tissue, the theoretical underpinnings have not been fully explored. We show that a two-dimensional version of the model, which is commonly employed to study tissue monolayers, can be derived from a three-dimensional version in the presence of a substrate. We also show how viscous forces, which arise from friction between different cells, can be included in the model. Finally, we numerically simulate a tissue monolayer and find that intercellular friction tends to solidify the tissue.

Coexisting mechanisms of luminogenesis in pancreatic cancer-derived organoids

Lumens are crucial features of the tissue architecture in both the healthy exocrine pancreas, where ducts shuttle enzymes from the acini to the intestine, and in the precancerous lesions of the highly lethal pancreatic ductal adenocarcinoma (PDAC), similarly displaying lumens that can further develop into cyst-like structures. Branched pancreatic-cancer derived organoids capture key architectural features of both the healthy and diseased pancreas, including lumens. However, their transition from a solid mass of cells to a hollow tissue remains insufficiently explored. Here, we show that organoids display two orthogonal but complementary lumen formation mechanisms: one relying on fluid intake for multiple microlumen nucleation, swelling and fusion, and the other involving the death of a central cell population, thereby hollowing out cavities. These results shed further light on the processes of luminogenesis, deepening our understanding of the early formation of PDAC precancerous lesions, including cystic neoplasia.

Short-range correlation of stress chains near solid-to-liquid transition in active monolayers

Using a three-dimensional model of cell monolayers, we study the spatial organization of active stress chains as the monolayer transitions from a solid to a liquid state. The critical exponents that characterize this transition map the isotropic stress percolation onto the two-dimensional random percolation universality class, suggesting short-range stress correlations near this transition. This mapping is achieved via two distinct, independent pathways: (i) cell-cell adhesion and (ii) active traction forces. We unify our findings by linking the nature of this transition to high-stress fluctuations, distinctly linked to each pathway. The results elevate the importance of the transmission of mechanical information in dense active matter and provide a new context for understanding the non-equilibrium statistical physics of phase transition in active systems.

Cell sorting by active forces in a phase-field model of cell monolayers

Cell sorting, the segregation of cells with different properties into distinct domains, is a key phenomenon in biological processes such as embryogenesis. We use a phase-field model of a confluent cell layer to study the role of activity in cell sorting. We find that a mixture of cells with extensile or contractile dipolar activity, and which are identical apart from their activity, quickly sort into small, elongated patches which then grow slowly in time. We interpret the sorting as driven by the different diffusivity of the extensile and contractile cells, mirroring the ordering of Brownian particles connected to different hot and cold thermostats. We check that the free energy is not changed by either partial or complete sorting, thus confirming that activity can be responsible for the ordering even in the absence of thermodynamic mechanisms.

Elasticity tunes mechanical stress localization around active topological defects

Mechanical stresses are increasingly found to be associated with various biological functionalities. At the same time, topological defects are being identified across a diverse range of biological systems and are points of localized mechanical stress. It is therefore important to ask how mechanical stress localization around topological defects is controlled. Here, we use continuum simulations of nonequilibrium, fluctuating and active nematics to explore the patterns of stress localization, as well as their extent and intensity around topological defects. We find that by increasing the orientational elasticity of the material, the isotropic stress pattern around topological defects is changed substantially, from a stress dipole characterized by symmetric compression-tension regions around the core of the defect, to a localized stress monopole at the defect position. Moreover, we show that elastic anisotropy alters the extent and intensity of the stresses, and can result in the dominance of tension or compression around defects. Finally, including both nonequilibrium fluctuations and active stress generation, we find that the elastic constant tunes the relative effect of each, leading to the flipping of tension and compression regions around topological defects. This flipping of the tension-compression regions only by changing the elastic constant presents an interesting, simple, way of switching the dynamic behavior in active matter by changing a passive material property. We expect these findings to motivate further exploration tuning stresses in active biological materials by varying material properties of the constituent units.