Geometric constraints during epithelial jamming

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.

Ductile-to-brittle transition and yielding in soft amorphous materials: perspectives and open questions

Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.

Motility driven glassy dynamics in confluent epithelial monolayers

As wounds heal, embryos develop, cancer spreads, or asthma progresses, the cellular monolayer undergoes a glass transition between solid-like jammed and fluid-like flowing states. During some of these processes, the cells undergo an epithelial-to-mesenchymal transition (EMT): they acquire in-plane polarity and become motile. Thus, how motility drives the glassy dynamics in epithelial systems is critical for the EMT process. However, no analytical framework that is indispensable for deeper insights exists. Here, we develop such a theory inspired by a well-known glass theory. One crucial result of this work is that the confluency affects the effective persistence time-scale of active force, described by its rotational diffusivity, Deffr. Deffr differs from the bare rotational diffusivity, Dr, of the motile force due to cell shape dynamics, which acts to rectify the force dynamics: Deffr is equal to Dr when Dr is small and saturates when Dr is large. We test the theoretical prediction of Deffr and how it affects the relaxation dynamics in our simulations of the active Vertex model. This novel effect of Deffr is crucial to understanding the new and previously published simulation data of active glassy dynamics in epithelial monolayers.

Mechanical state transitions in the regulation of tissue form and function

From embryonic development, postnatal growth and adult homeostasis to reparative and disease states, cells and tissues undergo constant changes in genome activity, cell fate, proliferation, movement, metabolism and growth. Importantly, these biological state transitions are coupled to changes in the mechanical and material properties of cells and tissues, termed mechanical state transitions. These mechanical states share features with physical states of matter, liquids and solids. Tissues can switch between mechanical states by changing behavioural dynamics or connectivity between cells. Conversely, these changes in tissue mechanical properties are known to control cell and tissue function, most importantly the ability of cells to move or tissues to deform. Thus, tissue mechanical state transitions are implicated in transmitting information across biological length and time scales, especially during processes of early development, wound healing and diseases such as cancer. This Review will focus on the biological basis of tissue-scale mechanical state transitions, how they emerge from molecular and cellular interactions, and their roles in organismal development, homeostasis, regeneration and disease.

Multicell-Fold: geometric learning in folding multicellular life

During developmental processes such as embryogenesis, how a group of cells fold into specific structures, is a central question in biology. However, it remains a major challenge to understand and predict the behavior of every cell within the living tissue over time during such intricate processes. Here we present a geometric deep-learning model that can accurately capture the highly convoluted interactions among cells. We demonstrate that multicellular data can be represented with both granular and foam-like physical pictures through a unified graph data structure, considering both cellular interactions and cell junction networks. Using this model, we achieve interpretable 4-D morphological sequence alignment, and predicting cell rearrangements before they occur at single-cell resolution. Furthermore, using neural activation map and ablation studies, we demonstrate cell geometries and cell junction networks together regulate morphogenesis at single-cell precision. This approach offers a pathway toward a unified dynamic atlas for a variety of developmental processes.

Learning Dynamics from Multicellular Graphs with Deep Neural Networks

Multicellular self-assembly into functional structures is a dynamic process that is critical in the development and diseases, including embryo development, organ formation, tumor invasion, and others. Being able to infer collective cell migratory dynamics from their static configuration is valuable for both understanding and predicting these complex processes. However, the identification of structural features that can indicate multicellular motion has been difficult, and existing metrics largely rely on physical instincts. Here we show that using a graph neural network (GNN), the motion of multicellular collectives can be inferred from a static snapshot of cell positions, in both experimental and synthetic datasets.

A shape-driven reentrant jamming transition in confluent monolayers of synthetic cell-mimics

Many critical biological processes, like wound healing, require densely packed cell monolayers/tissues to transition from a jammed solid-like to a fluid-like state. Although numerical studies anticipate changes in the cell shape alone can lead to unjamming, experimental support for this prediction is not definitive because, in living systems, fluidization due to density changes cannot be ruled out. Additionally, a cell's ability to modulate its motility only compounds difficulties since even in assemblies of rigid active particles, changing the nature of self-propulsion has non-trivial effects on the dynamics. Here, we design and assemble a monolayer of synthetic cell-mimics and examine their collective behaviour. By systematically increasing the persistence time of self-propulsion, we discovered a cell shape-driven, density-independent, re-entrant jamming transition. Notably, we observed cell shape and shape variability were mutually constrained in the confluent limit and followed the same universal scaling as that observed in confluent epithelia. Dynamical heterogeneities, however, did not conform to this scaling, with the fast cells showing suppressed shape variability, which our simulations revealed is due to a transient confinement effect of these cells by their slower neighbors. Our experiments unequivocally establish a morphodynamic link, demonstrating that geometric constraints alone can dictate epithelial jamming/unjamming.

Vimentin promotes collective cell migration through collagen networks via increased matrix remodeling and spheroid fluidity

The intermediate filament (IF) protein vimentin is associated with many diseases with phenotypes of enhanced cellular migration and aggressive invasion through the extracellular matrix (ECM) of tissues, but vimentin's role in in-vivo cell migration is still largely unclear. Vimentin is important for proper cellular adhesion and force generation, which are critical to cell migration; yet the vimentin cytoskeleton also hinders the ability of cells to squeeze through small pores in ECM, resisting migration. To identify the role of vimentin in collective cell migration, we generate spheroids of wide-type and vimentin-null mouse embryonic fibroblasts (mEFs) and embed them in a 3D collagen matrix. We find that loss of vimentin significantly impairs the ability of the spheroid to collectively expand through collagen networks and remodel the collagen network. Traction force analysis reveals that vimentin null spheroids exert less contractile force than their wild-type counterparts. In addition, spheroids made of mEFs with only vimentin unit length filaments (ULFs) exhibit similar behavior as vimentin-null spheroids, suggesting filamentous vimentin is required to promote 3D collective cell migration. We find the vimentin-mediated collective cell expansion is dependent on matrix metalloproteinase (MMP) degradation of the collagen matrix. Further, 3D vertex model simulation of spheroid and embedded ECM indicates that wild-type spheroids behave more fluid-like, enabling more active pulling and reconstructing the surrounding collagen network. Altogether, these results signify that VIF plays a critical role in enhancing migratory persistence in 3D matrix environments through MMP transportation and tissue fluidity.

Emergence of collective behaviours from local Voronoi topological perception

This article addresses how diverse collective behaviours arise from simple and realistic decisions made entirely at the level of each agent's personal space in the sense of the Voronoi diagram. We present a discrete-time model in two dimensions in which individual agents are aware of their local Voronoi environment and may seek static target locations. In particular, agents only communicate directly with their Voronoi neighbours and make decisions based on the geometry of their own Voronoi cells. With two effective control parameters, it is shown numerically to capture a wide range of collective behaviours in different scenarios. Further, we show that the Voronoi topology facilitates the computation of several novel observables for quantifying discrete collective behaviours. These observables are applicable to all agent-based models and to empirical data.

Geometry of adipocyte packing in subcutaneous tissue contributes to nonlinear tissue properties captured through a Gaussian process surrogate model

Subcutaneous tissue mechanical response is governed by the geometry and mechanical properties at the microscale and drives physiological and clinical processes such as drug delivery. Even though adipocyte packing is known to change with age, disease, and from one individual to another, the link between the geometry of the packing and the overall mechanical response of adipose tissue remains poorly understood. Here we create 1200 periodic representative volume elements (RVEs) that sample the possible space of Laguerre packings describing adipose tissue. RVE mechanics are modeled under tri-axial loading. Equilibrium configuration of RVEs is solved by minimizing an energetic potential that includes volume change contributions from adipocyte expansion, and area change contributions from collagen foam stretching. The resulting mechanical response across all RVE samples is interpolated with the aid of a Gaussian process (GP), revealing how the microscale geometry dictates the overall RVE mechanics. For example, increase in adipocyte size and increase in sphericity lead to adipose tissue softening. We showcase the use of the homogenized model in finite element simulations of drug injection by implementing a Blatz-Ko model, informed by the GP, as a custom material in the popular open-source package FEBio. These simulations show how microscale geometry can lead to vastly different injection dynamics even if the constituent parameters are held constant, highlighting the importance of characterizing individual's adipose tissue structure in the development of personalized therapies.

Cell Division and Motility Enable Hexatic Order in Biological Tissues

Biological tissues transform between solid- and liquidlike states in many fundamental physiological events. Recent experimental observations further suggest that in two-dimensional epithelial tissues these solid-liquid transformations can happen via intermediate states akin to the intermediate hexatic phases observed in equilibrium two-dimensional melting. The hexatic phase is characterized by quasi-long-range (power-law) orientational order but no translational order, thus endowing some structure to an otherwise structureless fluid. While it has been shown that hexatic order in tissue models can be induced by motility and thermal fluctuations, the role of cell division and apoptosis (birth and death) has remained poorly understood, despite its fundamental biological role. Here we study the effect of cell division and apoptosis on global hexatic order within the framework of the self-propelled Voronoi model of tissue. Although cell division naively destroys order and active motility facilitates deformations, we show that their combined action drives a liquid-hexatic-liquid transformation as the motility increases. The hexatic phase is accessed by the delicate balance of dislocation defect generation from cell division and the active binding of disclination-antidisclination pairs from motility. We formulate a mean-field model to elucidate this competition between cell division and motility and the consequent development of hexatic order.

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.

Regulation of Chromatin Modifications through Coordination of Nucleus Size and Epithelial Cell Morphology Heterogeneity

Cell morphology heterogeneity within epithelial collectives is a pervasive phenomenon intertwined with tissue mechanical properties. Despite its widespread occurrence, the underlying mechanisms driving cell morphology heterogeneity and its consequential biological ramifications remain elusive. Here, we investigate the dynamic evolution of epithelial cell morphology and nucleus morphology during crowding, unveiling a consistent correlation between the two. Our investigation reveals a persistent log-normal probability distribution characterizing both cell and nucleus areas across diverse crowding stages and epithelial model systems. We showed that this morphological diversity arises from asymmetric partitioning during cell division and is perpetuated through actomyosin-mediated regulation of cell-nucleus size coordination. Moreover, we provide insights into the impact of nucleus morphology on chromatin dynamics, demonstrating that constraining nucleus area leads to downregulation of the euchromatic mark H3K9ac and upregulation of the heterochromatic mark H3K27me3 through modulation of histone demethylase UTX expression. These findings under-score the significance of cell morphology heterogeneity as a driver of chromatin state diversity, shaping functional variability within epithelial tissues.

Nuclear size-regulated emergence of topological packing order on growing human lung alveolospheres

Within multicellular living systems, cells coordinate their positions with spatiotemporal accuracy to form various structures, setting the clock to control developmental processes and trigger maturation. These arrangements can be regulated by tissue topology, biochemical cues, as well as mechanical perturbations. However, the fundamental rules of how local cell packing order is regulated in forming three-dimensional (3D) multicellular architectures remain unclear. Furthermore, how cellular coordination evolves during developmental processes, and whether this cell patterning behavior is indicative of more complex biological functions, is largely unknown. Here, using human lung alveolospheres as a model system, by combining experiments and numerical simulations, we find that, surprisingly, cell packing behavior on alveolospheres resembles hard-disk packing but with increased randomness; the stiffer cell nuclei act as the hard disks surrounded by deformable cell bodies. Interestingly, we observe the emergence of topological packing order during alveolosphere growth, as a result of increasing nucleus-to-cell size ratio. Specifically, we find more hexagon-concentrated cellular packing with increasing bond orientational order, indicating a topological gas-to-liquid transition. Additionally, by osmotically changing the compactness of cells on alveolospheres, we observe that the variations in packing order align with the change of nucleus-to-cell size ratio. Together, our findings reveal the underlying rules of cell coordination and topological phases during human lung alveolosphere growth. These static packing characteristics are consistent with cell dynamics, together suggesting that better cellular packing stabilizes local cell neighborhoods and may regulate more complex biological functions such as organ development and cellular maturation.

Matrix confinement modulates 3D spheroid sorting and burst-like collective migration

While it is known that cells with differential adhesion tend to segregate and preferentially sort, the physical forces governing sorting and invasion in heterogeneous tumors remain poorly understood. To investigate this, we tune matrix confinement, mimicking changes in the stiffness and confinement of the tumor microenvironment, to explore how physical confinement influences individual and collective cell migration in 3D spheroids. High levels of confinement lead to cell sorting while reducing matrix confinement triggers the collective fluidization of cell motion. Cell sorting, which depends on cell-cell adhesion, is crucial to this phenomenon. Burst-like migration does not occur for spheroids that have not undergone sorting, regardless of the degree of matrix confinement. Using computational Self-Propelled Voronoi modeling, we show that spheroid sorting and invasion into the matrix depend on the balance between cell-generated forces and matrix resistance. The findings support a model where matrix confinement modulates 3D spheroid sorting and unjamming in an adhesion-dependent manner, providing insights into the mechanisms of cell sorting and migration in the primary tumor and toward distant metastatic sites. STATEMENT OF SIGNIFICANCE: The mechanical properties of the tumor microenvironment significantly influence cancer cell migration within the primary tumor, yet how these properties affect intercellular interactions in heterogeneous tumors is not well understood. By utilizing calcium and calcium chelators, we dynamically alter collagen-alginate hydrogel stiffness and investigate tumor cell behavior within co-culture spheroids in response to varying degrees of matrix confinement. High confinement is found to trigger cell sorting while reducing confinement for sorted spheroids facilitates collective cell invasion. Notably, without prior sorting, spheroids do not exhibit burst-like migration, regardless of confinement levels. This work establishes that matrix confinement and intercellular adhesion regulate 3D spheroid dynamics, offering insights into cellular organization and migration within the primary tumor.

Epithelial cells crowded out in asthma

Bronchoconstriction causes epithelial cell extrusion that promotes airway inflammation.

"Patchiness" in mechanical stiffness across a tumor as an early-stage marker for malignancy

Mechanical phenotyping of tumors, either at an individual cell level or tumor cell population level is gaining traction as a diagnostic tool. However, the extent of diagnostic and prognostic information that can be gained through these measurements is still unclear. In this work, we focus on the heterogeneity in mechanical properties of cells obtained from a single source such as a tissue or tumor as a potential novel biomarker. We believe that this heterogeneity is a conventionally overlooked source of information in mechanical phenotyping data. We use mechanics-based in-silico models of cell-cell interactions and cell population dynamics within 3D environments to probe how heterogeneity in cell mechanics drives tissue and tumor dynamics. Our simulations show that the initial heterogeneity in the mechanical properties of individual cells and the arrangement of these heterogenous sub-populations within the environment can dictate overall cell population dynamics and cause a shift towards the growth of malignant cell phenotypes within healthy tissue environments. The overall heterogeneity in the cellular mechanotype and their spatial distributions is quantified by a "patchiness" index, which is the ratio of the global to local heterogeneity in cell populations. We observe that there exists a threshold value of the patchiness index beyond which an overall healthy population of cells will show a steady shift towards a more malignant phenotype. Based on these results, we propose that the "patchiness" of a tumor or tissue sample, can be an early indicator for malignant transformation and cancer occurrence in benign tumors or healthy tissues. Additionally, we suggest that tissue patchiness, measured either by biochemical or biophysical markers, can become an important metric in predicting tissue health and disease likelihood just as landscape patchiness is an important metric in ecology.

Machine learning topological defects in confluent tissues

Active nematics is an emerging paradigm for characterizing biological systems. One aspect of particularly intense focus is the role active nematic defects play in these systems, as they have been found to mediate a growing number of biological processes. Accurately detecting and classifying these defects in biological systems is, therefore, of vital importance to improving our understanding of such processes. While robust methods for defect detection exist for systems of elongated constituents, other systems, such as epithelial layers, are not well suited to such methods. Here, we address this problem by developing a convolutional neural network to detect and classify nematic defects in confluent cell layers. Crucially, our method is readily implementable on experimental images of cell layers and is specifically designed to be suitable for cells that are not rod shaped, which we demonstrate by detecting defects on experimental data using the trained model. We show that our machine learning model outperforms current defect detection techniques and that this manifests itself in our method as requiring less data to accurately capture defect properties. This could drastically improve the accuracy of experimental data interpretation while also reducing costs, advancing the study of nematic defects in biological systems.

A gel-coated air-liquid-interface culture system with tunable substrate stiffness matching healthy and diseased lung tissues

Since its invention in the late 1980s, the air-liquid-interface (ALI) culture system has been the standard in vitro model for studying human airway biology and pulmonary diseases. However, in a conventional ALI system, cells are cultured on a porous plastic membrane that is much stiffer than human airway tissues. Here, we develop a gel-ALI culture system by simply coating the plastic membrane with a thin layer of hydrogel with tunable stiffness matching that of healthy and fibrotic airway tissues. We determine the optimum gel thickness that does not impair the transport of nutrients and biomolecules essential to cell growth. We show that the gel-ALI system allows human bronchial epithelial cells (HBECs) to proliferate and differentiate into pseudostratified epithelium. Furthermore, we discover that HBECs migrate significantly faster on hydrogel substrates with stiffness matching that of fibrotic lung tissues, highlighting the importance of mechanical cues in human airway remodeling. The developed gel-ALI system provides a facile approach to studying the effects of mechanical cues in human airway biology and in modeling pulmonary diseases.NEW & NOTEWORTHY In a conventional ALI system, cells are cultured on a plastic membrane that is much stiffer than human airway tissues. We develop a gel-ALI system by coating the plastic membrane with a thin layer of hydrogel with tunable stiffness matching that of healthy and fibrotic airway tissues. We discover that human bronchial epithelial cells migrate significantly faster on hydrogel substrates with pathological stiffness, highlighting the importance of mechanical cues in human airway remodeling.

Signatures of structural disorder in developing epithelial tissues

Epithelial cells generate functional tissues in developing embryos through collective movements and shape changes. In some morphogenetic events, a tissue dramatically reorganizes its internal structure - often generating high degrees of structural disorder - to accomplish changes in tissue shape. However, the origins of structural disorder in epithelia and what roles it might play in morphogenesis are poorly understood. We study this question in the Drosophila germband epithelium, which undergoes dramatic changes in internal structure as cell rearrangements drive elongation of the embryo body axis. Using two order parameters that quantify volumetric and shear disorder, we show that structural disorder increases during body axis elongation and is strongly linked with specific developmental processes. Both disorder metrics begin to increase around the onset of axis elongation, but then plateau at values that are maintained throughout the process. Notably, the disorder plateau values for volumetric disorder are similar to those for random cell packings, suggesting this may reflect a limit on tissue behavior. In mutant embryos with disrupted external stresses from the ventral furrow, both disorder metrics reach wild-type maximum disorder values with a delay, correlating with delays in cell rearrangements. In contrast, in mutants with disrupted internal stresses and cell rearrangements, volumetric disorder is reduced compared to wild type, whereas shear disorder depends on specific external stress patterns. Together, these findings demonstrate that internal and external stresses both contribute to epithelial tissue disorder and suggest that the maximum values of disorder in a developing tissue reflect physical or biological limits on morphogenesis.

Pulling the strings on solid-to-liquid phase transitions in cell collectives

Cell collectives must dynamically adapt to different biological contexts. For instance, in homeostatic conditions, epithelia must establish a barrier between body compartments and resist external stresses, while during development, wound healing or cancer invasion, these tissues undergo extensive remodeling. Using analogies from inert, passive materials, changes in cellular density, shape, rearrangements and/or migration were shown to result in collective transitions between solid and fluid states. However, what biological mechanisms govern these transitions remains an open question. In particular, the upstream signaling pathways and molecular effectors controlling the key physical axes determining tissue rheology and dynamics remain poorly understood. In this perspective, we focus on emerging evidence identifying the first biological signals determining the collective state of living tissues, with an emphasis on how these mechanisms are exploited for functionality across biological contexts.

Hydrodynamics and multiscale order in confluent epithelia

We formulate a hydrodynamic theory of confluent epithelia: i.e. monolayers of epithelial cells adhering to each other without gaps. Taking advantage of recent progresses toward establishing a general hydrodynamic theory of p-atic liquid crystals, we demonstrate that collectively migrating epithelia feature both nematic (i.e. p = 2) and hexatic (i.e. p = 6) orders, with the former being dominant at large and the latter at small length scales. Such a remarkable multiscale liquid crystal order leaves a distinct signature in the system's structure factor, which exhibits two different power-law scaling regimes, reflecting both the hexagonal geometry of small cells clusters and the uniaxial structure of the global cellular flow. We support these analytical predictions with two different cell-resolved models of epithelia - i.e. the self-propelled Voronoi model and the multiphase field model - and highlight how momentum dissipation and noise influence the range of fluctuations at small length scales, thereby affecting the degree of cooperativity between cells. Our construction provides a theoretical framework to conceptualize the recent observation of multiscale order in layers of Madin-Darby canine kidney cells and pave the way for further theoretical developments.

Differences and similarities in biophysical and biological characteristics between U87 MG glioblastoma and astrocyte cells

Current cancer studies focus on molecular-targeting diagnostics and interactions with surroundings; however, there are still gaps in characterization based on topological differences and elemental composition. Glioblastoma (GBM cells; GBMCs) is an astrocytic aggressive brain tumor. At the molecular level, GBMCs and astrocytes may differ, and cell elemental/topological analysis is critical for identifying potential new cancer targets. Here, we used U87 MG cells for GBMCS. U87 MG cell lines, which are frequently used in glioblastoma research, are an important tool for studying the various features and underlying mechanisms of this aggressive brain tumor. For the first time, atomic force microscopy (AFM), scanning electron microscopy (SEM) accompanied by energy-dispersive X-ray spectroscopy (EDS), and X-ray photoelectron spectroscopy (XPS) are used to report the topology and chemistry of cancer (U87 MG) and healthy (SVG p12) cells. In addition, F-actin staining and cytoskeleton-based gene expression analyses were performed. The degree of gene expression for genes related to the cytoskeleton was similar; however, the intensity of F-actin, anisotropy values, and invasion-related genes were different. Morphologically, GBMCs were longer and narrower while astrocytes were shorter and more disseminated based on AFM. Furthermore, the roughness values of these cells differed slightly between the two call types. In contrast to the rougher astrocyte surfaces in the lamellipodial area, SEM-EDS analysis showed that elongated GBMCs displayed filopodial protrusions. Our investigation provides considerable further insight into rapid cancer cell characterization in terms of a combinatorial spectroscopic and microscopic approach.

Discontinuous Shear Thickening in Biological Tissue Rheology

During embryonic morphogenesis, tissues undergo dramatic deformations in order to form functional organs. Similarly, in adult animals, living cells and tissues are continually subjected to forces and deformations. Therefore, the success of embryonic development and the proper maintenance of physiological functions rely on the ability of cells to withstand mechanical stresses as well as their ability to flow in a collective manner. During these events, mechanical perturbations can originate from active processes at the single-cell level, competing with external stresses exerted by surrounding tissues and organs. However, the study of tissue mechanics has been somewhat limited to either the response to external forces or to intrinsic ones. In this work, we use an active vertex model of a 2D confluent tissue to study the interplay of external deformations that are applied globally to a tissue with internal active stresses that arise locally at the cellular level due to cell motility. We elucidate, in particular, the way in which this interplay between globally external and locally internal active driving determines the emergent mechanical properties of the tissue as a whole. For a tissue in the vicinity of a solid-fluid jamming or unjamming transition, we uncover a host of fascinating rheological phenomena, including yielding, shear thinning, continuous shear thickening, and discontinuous shear thickening. These model predictions provide a framework for understanding the recently observed nonlinear rheological behaviors in vivo.

Minimal vertex model explains how the amnioserosa avoids fluidization during Drosophila dorsal closure

Dorsal closure is a process that occurs during embryogenesis of Drosophila melanogaster . During dorsal closure, the amnioserosa (AS), a one-cell thick epithelial tissue that fills the dorsal opening, shrinks as the lateral epidermis sheets converge and eventually merge. During this process, the aspect ratio of amnioserosa cells increases markedly. The standard 2-dimensional vertex model, which successfully describes tissue sheet mechanics in multiple contexts, would in this case predict that the tissue should fluidize via cell neighbor changes. Surprisingly, however, the amnioserosa remains an elastic solid with no such events. We here present a minimal extension to the vertex model that explains how the amnioserosa can achieve this unexpected behavior. We show that continuous shrinkage of the preferred cell perimeter and cell perimeter polydispersity lead to the retention of the solid state of the amnioserosa. Our model accurately captures measured cell shape and orientation changes and predicts non-monotonic junction tension that we confirm with laser ablation experiments.

Significance Statement

During embryogenesis, cells in tissues can undergo significant shape changes. Many epithelial tissues fluidize, i.e. cells exchange neighbors, when the average cell aspect ratio increases above a threshold value, consistent with the standard vertex model. During dorsal closure in Drosophila melanogaster , however, the amnioserosa tissue remains solid even as the average cell aspect ratio increases well above threshold. We introduce perimeter polydispersity and allow the preferred cell perimeters, usually held fixed in vertex models, to decrease linearly with time as seen experimentally. With these extensions to the standard vertex model, we capture experimental observations quantitatively. Our results demonstrate that vertex models can describe the behavior of the amnioserosa in dorsal closure by allowing normally fixed parameters to vary with time.

Minimal vertex model explains how the amnioserosa avoids fluidization during Drosophila dorsal closure

Dorsal closure is a process that occurs during embryogenesis of Drosophila melanogaster. During dorsal closure, the amnioserosa (AS), a one-cell thick epithelial tissue that fills the dorsal opening, shrinks as the lateral epidermis sheets converge and eventually merge. During this process, the aspect ratio of amnioserosa cells increases markedly. The standard 2-dimensional vertex model, which successfully describes tissue sheet mechanics in multiple contexts, would in this case predict that the tissue should fluidize via cell neighbor changes. Surprisingly, however, the amnioserosa remains an elastic solid with no such events. We here present a minimal extension to the vertex model that explains how the amnioserosa can achieve this unexpected behavior. We show that continuous shrink-age of the preferred cell perimeter and cell perimeter polydispersity lead to the retention of the solid state of the amnioserosa. Our model accurately captures measured cell shape and orientation changes and predicts non-monotonic junction tension that we confirm with laser ablation experiments.

Bridging the gap between collective motility and epithelial-mesenchymal transitions through the active finite voronoi model

We introduce an active version of the recently proposed finite Voronoi model of epithelial tissue. The resultant Active Finite Voronoi (AFV) model enables the study of both confluent and non-confluent geometries and transitions between them, in the presence of active cells. Our study identifies six distinct phases, characterized by aggregation-segregation, dynamical jamming-unjamming, and epithelial-mesenchymal transitions (EMT), thereby extending the behavior beyond that observed in previously studied vertex-based models. The AFV model with rich phase diagram provides a cohesive framework that unifies the well-observed progression to collective motility via unjamming with the intricate dynamics enabled by EMT. This approach should prove useful for challenges in developmental biology systems as well as the complex context of cancer metastasis. The simulation code is also provided.

Using Histologic Image Analysis to Understand Biophysical Regulations of Epithelial Cell Morphology

Epithelial mechanics and mechanobiology have become 2 important research fields in life sciences and bioengineering. These fields investigate how physical factors induced by cell adhesion and collective behaviors can directly regulate biologic processes, such as organ development and disease progression. Cell mechanics and mechanobiology thus make exciting biophysics education topics to illustrate how fundamental physics principles play a role in regulating cell biology. However, the field currently lacks hands-on activities that engage students in learning science and outreach programs in these topics. One such area is the development of robust hands-on modules that allow students to observe features of cell shape and mechanics and connect them to fundamental physics principles. Here, we demonstrate a workflow that engages students in studying epithelial cell mechanics by using commercial histology slides of frog skin. We show that by using recently developed artificial intelligence-based image-segmentation tools, students can easily quantify different cell morphologic features in a high-throughput manner. Using our workflow, students can reproduce 2 essential findings in cell mechanics: the common gamma distribution of normalized cell aspect ratio in jammed epithelia and the constant ratio between the nuclear and cellular area. Importantly, because the only required instrument for this active learning module is a readily available light microscope and a computer, our module is relatively low cost, as well as portable. These features make the module scalable for students at various education levels and outreach programs. This highly accessible education module provides a fun and engaging way to introduce students to the world of epithelial tissue mechanics.

Bayesian traction force estimation using cell boundary-dependent force priors

Understanding the principles of cell migration necessitates measurements of the forces generated by cells. In traction force microscopy (TFM), fluorescent beads are placed on a substrate's surface and the substrate strain caused by the cell traction force is observed as displacement of the beads. Mathematical analysis can estimate traction force from bead displacement. However, most algorithms estimate substrate stresses independently of cell boundary, which results in poor estimation accuracy in low-density bead environments. To achieve accurate force estimation at low density, we proposed a Bayesian traction force estimation (BTFE) algorithm that incorporates cell-boundary-dependent force as a prior. We evaluated the performance of the proposed algorithm using synthetic data generated with mathematical models of cells and TFM substrates. BTFE outperformed other methods, especially in low-density bead conditions. In addition, the BTFE algorithm provided a reasonable force estimation using TFM images from the experiment.

A life off the beaten track in biomechanics: Imperfect elasticity, cytoskeletal glassiness, and epithelial unjamming

Textbook descriptions of elasticity, viscosity, and viscoelasticity fail to account for certain mechanical behaviors that typify soft living matter. Here, we consider three examples. First, strong empirical evidence suggests that within lung parenchymal tissues, the frictional stresses expressed at the microscale are fundamentally not of viscous origin. Second, the cytoskeleton (CSK) of the airway smooth muscle cell, as well as that of all eukaryotic cells, is more solid-like than fluid-like, yet its elastic modulus is softer than the softest of soft rubbers by a factor of 104-105. Moreover, the eukaryotic CSK expresses power law rheology, innate malleability, and fluidization when sheared. For these reasons, taken together, the CSK of the living eukaryotic cell is reminiscent of the class of materials called soft glasses, thus likening it to inert materials such as clays, pastes slurries, emulsions, and foams. Third, the cellular collective comprising a confluent epithelial layer can become solid-like and jammed, fluid-like and unjammed, or something in between. Esoteric though each may seem, these discoveries are consequential insofar as they impact our understanding of bronchospasm and wound healing as well as cancer cell invasion and embryonic development. Moreover, there are reasons to suspect that certain of these phenomena first arose in the early protist as a result of evolutionary pressures exerted by the primordial microenvironment. We have hypothesized, further, that each then became passed down virtually unchanged to the present day as a conserved core process. These topics are addressed here not only because they are interesting but also because they track the journey of one laboratory along a path less traveled by.

Nature-inspired designs for disordered acoustic bandgap materials

We introduce an amorphous mechanical metamaterial inspired by how cells pack in biological tissues. The spatial heterogeneity in the local stiffness of these materials has been recently shown to impact the mechanics of confluent biological tissues and cancer tumor invasion. Here we use this bio-inspired structure as a design template to construct mechanical metamaterials and show that this heterogeneity can give rise to amorphous cellular solids with large, tunable acoustic bandgaps. Unlike acoustic crystals with periodic structures, the bandgaps here are directionally isotropic and robust to defects due to their complete lack of positional order. Possible ways to manipulate bandgaps are explored with a combination of the tissue-level elastic modulus and local stiffness heterogeneity of cells. To further demonstrate the existence of bandgaps, we dynamically perturb the system with an external sinusoidal wave in the perpendicular and horizontal directions. The transmission coefficients are calculated and show valleys that coincide with the location of bandgaps. Experimentally this design should lead to the engineering of self-assembled rigid acoustic structures with full bandgaps that can be controlled via mechanical tuning and promote applications in a broad area from vibration isolations to mechanical waveguides.

Collective heterogeneity of mitochondrial potential in contact inhibition of proliferation

In the epithelium, cell density and cell proliferation are closely connected to each other through contact inhibition of proliferation (CIP). Depending on cell density, CIP proceeds through three distinct stages: the free-growing stage at low density, the pre-epithelial transition stage at medium density, and the post-epithelial transition stage at high density. Previous studies have elucidated how cell morphology, motion, and mechanics vary in these stages. However, it remains unknown whether cellular metabolism also has a density-dependent behavior. By measuring the mitochondrial membrane potential at different cell densities, here we reveal a heterogeneous landscape of metabolism in the epithelium, which appears qualitatively distinct in three stages of CIP and did not follow the trend of other CIP-associated parameters, which increases or decreases monotonically with increasing cell density. Importantly, epithelial cells established a collective metabolic heterogeneity exclusively in the pre-epithelial transition stage, where the multicellular clusters of high- and low-potential cells emerged. However, in the post-epithelial transition stage, the metabolic potential field became relatively homogeneous. Next, to study the underlying dynamics, we constructed a system biology model, which predicted the role of cell proliferation in metabolic potential toward establishing collective heterogeneity. Further experiments indeed revealed that the metabolic pattern spatially correlated with the proliferation capacity of cells, as measured by the nuclear localization of a pro-proliferation protein, YAP. Finally, experiments perturbing the actomyosin contractility revealed that, while metabolic heterogeneity was maintained in the absence of actomyosin contractility, its ab initio emergence depended on the latter. Taken together, our results revealed a density-dependent collective heterogeneity in the metabolic field of a pre-epithelial transition-stage epithelial monolayer, which may have significant implications for epithelial form and function.

HDAC9/p300/F-actin immunoexpression and migration analysis for malignant melanoma stem cell

Melanoma is an aggressive tumor with a poor prognosis that worsens in the metastatic phase. Distruptions of epigenetic mechanisms is known to effect cancer stem cells (CSCs) activity. Malignant melanoma (MM) progression may be promoted by changes in the genetic structure of CSC. Thus, treatments that target epigenetic modifications could be a promising weapon, especially in melanoma. Here, we compared p300, HDAC9, and F-actin proteins in melanoma CSCs (CD133+), non-CSCs (CD133-) and CHL-1 cell line, as well as cell migration and division rates. At 4 and 6 h, P300 protein levels in CHL-1 and CD133 + were remarkably similar, and the CD133- showed increases in expression levels as the incubation period lengthened. HDAC9 protein intensity decreased in CHL-1, increased in the CD133-, and remained relatively unchanged in the CD133+ as the incubation period lengthened. The mean value of F-actin expression level increased in all cell group with time, when the highest increase observed in CHL-1. In conclusion, our studies contribute to the management of metastatic diseases in the future and offer new insight into the molecular basis of the initiation and progression of MM.

Long-term, real-time and label-free live cell image processing and analysis based on a combined algorithm of CellPose and watershed segmentation

Developing a rapid and quantitative method to accurately evaluate the physiological abilities of living cells is critical for tumor control. Many experiments have been conducted in the field of biology in an attempt to measure the proliferation and movement abilities of cells, but existing methods cannot provide real-time and objective data for label-free cells. The quantitative imaging technique, including an automatic segmentation algorithm for individual label-free cells, has been a breakthrough in this regard. In this study, we develop a combined automatic image processing algorithm of CellPose and watershed segmentation for the long-term and real-time imaging of label-free cells. This method shows strong reliability in cell identification regardless of cell densities, allowing us to obtain accurate information about the number and proliferation ability of the target cells. Additionally, our results also suggest that this method is a reliable way to assess real-time data on drug cytotoxicity, cell morphology, and cell movement ability.

Collective cell behaviors manipulated by synthetic DNA nanostructures

Cellular collective motion in confluent epithelial monolayers is involved in many processes such as embryo development, carcinoma invasion, and wound healing. The development of new chemical strategies to achieve large-scale control of cells' collective motion is essential for biomedical applications. Here a series of DNA nanostructures with different dimensions were synthesized and their influences on cells' collective migration and packing behaviors in epithelial monolayers were investigated. We found that the framed DNA nanoassemblies effectively reduced the cells' speed by increasing the rigidity of cells, while the lipid-DNA micelles had a more pronounced effect on cells' projection area and shape factor. These DNA nanostructures all significantly enhanced the dependence of cells' speed on their shape factor. Our results indicate that cells' mobility in monolayers can be manipulated by chemical intercellular interactions without any genetic intervention. This may provide a new chemical strategy for tissue engineering and tumor therapy.

Tissue flows are tuned by actomyosin-dependent mechanics in developing embryos

Rapid epithelial tissue flows are essential to building and shaping developing embryos. However, the mechanical properties of embryonic epithelial tissues and the factors that control these properties are not well understood. Actomyosin generates contractile tensions and contributes to the mechanical properties of cells and cytoskeletal networks in vitro, but it remains unclear how the levels and patterns of actomyosin activity contribute to embryonic epithelial tissue mechanics in vivo. To dissect the roles of cell-generated tensions in the mechanics of flowing epithelial tissues, we use optogenetic tools to manipulate actomyosin contractility with spatiotemporal precision in the Drosophila germband epithelium, which rapidly flows during body axis elongation. We find that manipulating actomyosin-dependent tensions by either optogenetic activation or deactivation of actomyosin alters the solid-fluid mechanical properties of the germband epithelium, leading to changes in cell rearrangements and tissue-level flows. Optogenetically activating actomyosin leads to increases in the overall level but decreases in the anisotropy of tension in the tissue, whereas optogenetically deactivating actomyosin leads to decreases in both the level and anisotropy of tension compared to in wild-type embryos. We find that optogenetically activating actomyosin results in more solid-like (less fluid-like) tissue properties, which is associated with reduced cell rearrangements and tissue flow compared to in wild-type embryos. Optogenetically deactivating actomyosin also results in more solid-like properties than in wild-type embryos but less solid-like properties compared to optogenetically activating actomyosin. Together, these findings indicate that increasing the overall tension level is associated with more solid-like properties in tissues that are relatively isotropic, whereas high tension anisotropy fluidizes the tissue. Our results reveal that epithelial tissue flows in developing embryos involve the coordinated actomyosin-dependent regulation of the mechanical properties of tissues and the tensions driving them to flow in order to achieve rapid tissue remodeling.

Local response and emerging nonlinear elastic length scale in biopolymer matrices

Nonlinear stiffening is a ubiquitous property of major types of biopolymers that make up the extracellular matrices (ECM) including collagen, fibrin, and basement membrane. Within the ECM, many types of cells such as fibroblasts and cancer cells have a spindle-like shape that acts like two equal and opposite force monopoles, which anisotropically stretch their surroundings and locally stiffen the matrix. Here, we first use optical tweezers to study the nonlinear force-displacement response to localized monopole forces. We then propose an effective-probe scaling argument that a local point force application can induce a stiffened region in the matrix, which can be characterized by a nonlinear length scale R* that increases with the increasing force magnitude; the local nonlinear force-displacement response is a result of the nonlinear growth of this effective probe that linearly deforms an increasing portion of the surrounding matrix. Furthermore, we show that this emerging nonlinear length scale R* can be observed around living cells and can be perturbed by varying matrix concentration or inhibiting cell contractility.

Multicellular aligned bands disrupt global collective cell behavior

Mechanical stress patterns emerging from collective cell behavior have been shown to play critical roles in morphogenesis, tissue repair, and cancer metastasis. In our previous work, we constrained valvular interstitial cell (VIC) monolayers on circular protein islands to study emergent behavior in a controlled manner and demonstrated that the general patterns of cell alignment, size, and apoptosis correlate with predicted mechanical stress fields if radially increasing stiffness or contractility are used in the computational models. However, these radially symmetric models did not predict the existence of local regions of dense aligned cells observed in seemingly random locations of individual aggregates. The goal of this study is to determine how the heterogeneities in cell behavior emerge over time and diverge from the predicted collective cell behavior. Cell-cell interactions in circular multicellular aggregates of VICs were studied with time-lapse imaging ranging from hours to days, and migration, proliferation, and traction stresses were measured. Our results indicate that elongated cells create strong local alignment within preconfluent cell populations on the microcontact printed protein islands. These cells influence the alignment of additional cells to create dense, locally aligned bands of cells which disrupt the predicted global behavior. Cells are highly elongated at the endpoints of the bands yet have decreased spread area in the middle and reduced mobility. Although traction stresses at the endpoints of bands are enhanced, even to the point of detaching aggregates from the culture surface, the cells in dense bands exhibit reduced proliferation, less nuclear YAP, and increased apoptotic rates indicating a low stress environment. These findings suggest that strong local cell-cell interactions between primary fibroblastic cells can disrupt the global collective cellular behavior leading to substantial heterogeneity of cell behaviors in constrained monolayers. This local emergent behavior within aggregated fibroblasts may play an important role in development and disease of connective tissues. STATEMENT OF SIGNIFICANCE: Mechanical stress patterns emerging from collective cell behavior play critical roles in morphogenesis, tissue repair, and cancer metastasis. Much has been learned of these collective behaviors by utilizing microcontact printing to constrain cell monolayers (aggregates) into specific shapes. Here we utilize these tools along with long-term video microscopy tracking of individual aggregates to determine how heterogeneous collective behaviors unique to primary fibroblastic cells emerge over time and diverge from computed stress fields. We find that dense multicellular bands form from local collective behavior and disrupt the global collective behavior resulting in heterogeneous patterns of migration, traction stresses, proliferation, and apoptosis. This local emergent behavior within aggregated fibroblasts may play an important role in development and disease of connective tissues.

Micro and Macroscopic Stress-Strain Relations in Disordered Tessellated Networks

We demonstrate that for a rigid and incompressible network in mechanical equilibrium, the microscopic stress and strain follows a simple relation, σ=pE, where σ is the deviatoric stress, E is a mean-field strain tensor, and p is the hydrostatic pressure. This relationship arises as the natural consequence of energy minimization or equivalently, mechanical equilibration. The result suggests not only that the microscopic stress and strain are aligned in the principal directions, but also microscopic deformations are predominantly affine. The relationship holds true regardless of the different (foam or tissue) energy model considered, and directly leads to a simple prediction for the shear modulus, μ=⟨p⟩/2, where ⟨p⟩ is the mean pressure of the tessellation, for general randomized lattices.

Quantification of basal stem cell elongation and stress fiber accumulation in the pseudostratified airway epithelium during the unjamming transition

Under homeostatic conditions, epithelial cells remain non-migratory. However, during embryonic development and pathological conditions, they become migratory. The mechanism underlying the transition of the epithelial layer between non-migratory and migratory phases is a fundamental question in biology. Using well-differentiated primary human bronchial epithelial cells that form a pseudostratified epithelium, we have previously identified that a confluent epithelial layer can transition from a non-migratory to migratory phase through an unjamming transition (UJT). We previously defined collective cellular migration and apical cell elongation as hallmarks of UJT. However, other cell-type-specific changes have not been previously studied in the pseudostratified airway epithelium, which consists of multiple cell types. Here, we focused on the quantifying morphological changes in basal stem cells during the UJT. Our data demonstrate that during the UJT, airway basal stem cells elongated and enlarged, and their stress fibers elongated and aligned. These morphological changes observed in basal stem cells correlated to the previously defined hallmarks of the UJT. Moreover, basal cell and stress fiber elongation were observed prior to apical cell elongation. Together, these morphological changes indicate that basal stem cells in pseudostratified airway epithelium are actively remodeling, presumably through accumulation of stress fibers during the UJT.

Transitions in density, pressure, and effective temperature drive collective cell migration into confining environments

Epithelial cell collectives migrate through tissue interfaces and crevices to orchestrate processes of development, tumor invasion, and wound healing. Naturally, traversal of cell collective through confining environments involves crowding due to the narrowing space, which seems tenuous given the conventional inverse relationship between cell density and migration. However, physical transitions required to overcome such epithelial densification for migration across confinements remain unclear. Here, in contiguous microchannels, we show that epithelial (MCF10A) monolayers accumulate higher cell density before entering narrower channels; however, overexpression of breast cancer oncogene +ErbB2 reduced this need for density accumulation across confinement. While wildtype MCF10A cells migrated faster in narrow channels, this confinement sensitivity reduced after +ErbB2 mutation or with constitutively-active RhoA. The migrating collective developed pressure differentials upon encountering microchannels, like fluid flow into narrowing spaces, and this pressure dropped with their continued migration. These transitions of pressure and density altered cell shapes and increased effective temperature, estimated by treating cells as granular thermodynamic system. While +RhoA cells and those in confined regions were effectively warmer, cancer-like +ErbB2 cells remained cooler. Epithelial reinforcement by metformin treatment increased density and temperature differentials across confinement, indicating that higher cell cohesion could reduce unjamming. Our results provide experimental evidence for previously proposed theories of inverse relationship between density and motility-related effective temperature. Indeed, we show across cell lines that confinement increases pressure and effective temperature, which enable migration by reducing density. This physical interpretation of collective cell migration as granular matter could advance our understanding of complex living systems.

Curvature-Induced Cell Rearrangements in Biological Tissues

On a curved surface, epithelial cells can adapt to geometric constraints by tilting and by exchanging their neighbors from apical to basal sides, known as an apico-basal topological transition 1 (AB-T1). The relationship between cell tilt, AB-T1s, and tissue curvature still lacks a unified understanding. Here, we propose a general framework for cell packing in curved environments and explain the formation of AB-T1s from the perspective of strain anisotropy. We find that steep curvature gradients can lead to cell tilting and induce AB-T1s. Alternatively, pressure differences across the epithelial tissue can drive AB-T1s in regions of large curvature anisotropy. The two mechanisms compete to determine the impact of tissue geometry and mechanics on optimized cell rearrangements in three dimensions.

Analysis of Collective Migration Patterns Within Tumors

Metastasis is a hallmark of cancer and the leading cause of mortality among cancer patients. Cancer, in its most deadly form, is thus not only a disease of uncontrolled cell growth but also a disease of uncontrolled cell migration. The study of tumor cell migration requires both experimental systems that are representative of the complex tumor environment as well as quantitative tools to analyze migration patterns. In this chapter, we focus on experimental and analytical methods to capture and analyze cell migration in live explants from mouse intestinal tumors. We first describe a protocol to extract and perform ex vivo live imaging on intestinal tumors in mice. We then provide a step-by-step image analysis workflow using freely available software and custom analysis scripts for extracting several parameters related to collective cell migration and cell and tissue organization.

Chloroplasts in plant cells show active glassy behavior under low-light conditions

Plants have developed intricate mechanisms to adapt to changing light conditions. Besides phototropism and heliotropism (differential growth toward light and diurnal motion with respect to sunlight, respectively), chloroplast motion acts as a fast mechanism to change the intracellular structure of leaf cells. While chloroplasts move toward the sides of the plant cell to avoid strong light, they accumulate and spread out into a layer on the bottom of the cell at low light to increase the light absorption efficiency. Although the motion of chloroplasts has been studied for over a century, the collective organelle motion leading to light-adapting self-organized structures remains elusive. Here, we study the active motion of chloroplasts under dim-light conditions, leading to an accumulation in a densely packed quasi-2D layer. We observe burst-like rearrangements and show that these dynamics resemble systems close to the glass transition by tracking individual chloroplasts. Furthermore, we provide a minimal mathematical model to uncover relevant system parameters controlling the stability of the dense configuration of chloroplasts. Our study suggests that the meta-stable caging close to the glass transition in the chloroplast monolayer serves a physiological relevance: Chloroplasts remain in a spread-out configuration to increase the light uptake but can easily fluidize when the activity is increased to efficiently rearrange the structure toward an avoidance state. Our research opens questions about the role that dynamical phase transitions could play in self-organized intracellular responses of plant cells toward environmental cues.

Orchestration of tissue shape changes and gene expression patterns in development

In development, tissue shape changes and gene expression patterns give rise to morphogenesis. Understanding tissue shape changes requires the analysis of mechanical properties of the tissue such as tissue rigidity, cell influx from neighboring tissues, cell shape changes and cell proliferation. Local and global gene expression patterns can be influenced by neighbor exchange and tissue shape changes. Here we review recent studies on the mechanisms for tissue elongation and its influences on dynamic gene expression patterns by focusing on vertebrate somitogenesis. We first introduce mechanical and biochemical properties of the segmenting tissue that drive tissue elongation. Then, we discuss patterning in the presence of cell mixing, scaling of signaling gradients, and dynamic phase waves of rhythmic gene expression under tissue shape changes. We also highlight the importance of theoretical approaches to address the relation between tissue shape changes and patterning.

Polyploid/Multinucleated Giant and Slow-Cycling Cancer Cell Enrichment in Response to X-ray Irradiation of Human Glioblastoma Multiforme Cells Differing in Radioresistance and TP53/PTEN Status

Radioresistance compromises the efficacy of radiotherapy for glioblastoma multiforme (GBM), the most devastating and common brain tumor. The present study investigated the relationship between radiation tolerance and formation of polyploid/multinucleated giant (PGCC/MGCC) and quiescent/senescent slow-cycling cancer cells in human U-87, LN-229, and U-251 cell lines differing in TP53/PTEN status and radioresistance. We found significant enrichment in MGCC populations of U-87 and LN-229 cell lines, and generation of numerous small mononuclear (called Raju cells, or RJ cells) U-87-derived cells that eventually form cell colonies, in a process termed neosis, in response to X-ray irradiation (IR) at single acute therapeutic doses of 2-6 Gy. For the first time, single-cell high-content imaging and analysis of Ki-67- and EdU-coupled fluorescence demonstrated that the IR exposure dose-dependently augments two distinct GBM cell populations. Bifurcation of Ki-67 staining suggests fast-cycling and slow-cycling populations with a normal-sized nuclear area, and with an enlarged nuclear area, including one resembling the size of PGCC/MGCCs, that likely underlie the highest radioresistance and propensity for repopulation of U-87 cells. Proliferative activity and anchorage-independent survival of GBM cell lines seem to be related to neosis, low level of apoptosis, fraction of prematurely stress-induced senescent MGCCs, and the expression of p63 and p73, members of p53 family transcription factors, but not to the mutant p53. Collectively, our data support the importance of the TP53wt/PTENmut genotype for the maintenance of cycling radioresistant U-87 cells to produce a significant amount of senescent MGCCs as an IR stress-induced adaptation response to therapeutic irradiation doses.

Cell jamming regulates epithelial chiral morphogenesis

Internal organs such as the heart demonstrate apparent left-right (LR) asymmetric morphology and positioning. Cellular chirality and associated LR biased mechanical behavior such as cell migration have been attributed to LR symmetry breaking during embryonic development. Mathematical models have shown that chiral directional migration can be driven by cellular intrinsic torque. Tissue jamming state (i.e., solid-like vs fluid-like state) strongly regulates collective migratory behavior, but how it might affect chiral morphogenesis is still unknown. Here, we develop a cell vertex model to study the role of tissue rigidity or jamming state on chiral morphogenesis of the cells on a patterned ring-shaped tissue, simulating a previously reported experimental setup for measuring cell chirality. We simulate chirality as torsional forces acting on cell vertices. As expected, the cells undergo bidirectional migration at the opposing (inner and outer) boundaries of the ring-shaped tissue. We discover that more fluid-like tissues (unjammed) demonstrate a stronger chiral cell alignment and elongation than more solid-like (jammed) tissues and maintain a bigger difference in migration velocity between opposing tissue boundaries. Finally, we find that fluid-like tissues undergo more cell-neighbor exchange events. This study reveals that chiral torque is sufficient to achieve a biased cellular alignment as seen in vitro. It further sheds light on the mechanical regulation of chiral morphogenesis of tissues and reveals a role of cell density-independent tissue rigidity in this process.

On the origin of universal cell shape variability in confluent epithelial monolayers

Cell shape is fundamental in biology. The average cell shape can influence crucial biological functions, such as cell fate and division orientation. But cell-to-cell shape variability is often regarded as noise. In contrast, recent works reveal that shape variability in diverse epithelial monolayers follows a nearly universal distribution. However, the origin and implications of this universality remain unclear. Here, assuming contractility and adhesion are crucial for cell shape, characterized via aspect ratio (r), we develop a mean-field analytical theory for shape variability. We find that all the system-specific details combine into a single parameter α that governs the probability distribution function (PDF) of r; this leads to a universal relation between the standard deviation and the average of r. The PDF for the scaled r is not strictly but nearly universal. In addition, we obtain the scaled area distribution, described by the parameter μ. Information of α and μ together can distinguish the effects of changing physical conditions, such as maturation, on different system properties. We have verified the theory via simulations of two distinct models of epithelial monolayers and with existing experiments on diverse systems. We demonstrate that in a confluent monolayer, average shape determines both the shape variability and dynamics. Our results imply that cell shape distribution is inevitable, where a single parameter describes both statics and dynamics and provides a framework to analyze and compare diverse epithelial systems. In contrast to existing theories, our work shows that the universal properties are consequences of a mathematical property and should be valid in general, even in the fluid regime.

Jamming Transitions in Astrocytes and Glioblastoma Are Induced by Cell Density and Tension

Collective behavior of cells emerges from coordination of cell-cell-interactions and is important to wound healing, embryonic and tumor development. Depending on cell density and cell-cell interactions, a transition from a migratory, fluid-like unjammed state to a more static and solid-like jammed state or vice versa can occur. Here, we analyze collective migration dynamics of astrocytes and glioblastoma cells using live cell imaging. Furthermore, atomic force microscopy, traction force microscopy and spheroid generation assays were used to study cell adhesion, traction and mechanics. Perturbations of traction and adhesion were induced via ROCK or myosin II inhibition. Whereas astrocytes resided within a non-migratory, jammed state, glioblastoma were migratory and unjammed. Furthermore, we demonstrated that a switch from an unjammed to a jammed state was induced upon alteration of the equilibrium between cell-cell-adhesion and tension from adhesion to tension dominated, via inhibition of ROCK or myosin II. Such behavior has implications for understanding the infiltration of the brain by glioblastoma cells and may help to identify new strategies to develop anti-migratory drugs and strategies for glioblastoma-treatment.

Picking winners in cell-cell collisions: Wetting, speed, and contact

Groups of eukaryotic cells can coordinate their crawling motion to follow cues more effectively, stay together, or invade new areas. This collective cell migration depends on cell-cell interactions, which are often studied by colliding pairs of cells together. Can the outcome of these collisions be predicted? Recent experiments on trains of colliding epithelial cells suggest that cells with a smaller contact angle to the surface or larger speeds are more likely to maintain their direction ("win") upon collision. When should we expect shape or speed to correlate with the outcome of a collision? To investigate this question, we build a model for two-cell collisions within the phase field framework, which allows for cell shape changes. We can reproduce the observation that cells with high speed and small contact angles are more likely to win with two different assumptions for how cells interact: (1) velocity aligning, in which we hypothesize that cells sense their own velocity and align to it over a finite timescale, and (2) front-front contact repolarization, where cells polarize away from cell-cell contact, akin to contact inhibition of locomotion. Surprisingly, though we simulate collisions between cells with widely varying properties, in each case, the probability of a cell winning is completely captured by a single summary variable: its relative speed (in the velocity-aligning model) or its relative contact angle (in the contact repolarization model). Both models are currently consistent with reported experimental results, but they can be distinguished by varying cell contact angle and speed through orthogonal perturbations.