In tissue development, wound healing, and cancer invasion, coordinated cell motion arises from active forces produced by the cells. The relationship between force and motion remains unclear, however, because the forces result from a sum of contributions from activity and the constitutive response of the cell collective. Here, we develop a method to decouple the forces due to activity from those due to constitutive response. As a model of an epithelial tissue, we use a monolayer of epithelial cells in the fluid state, for which the constitutive behavior is that of a viscous fluid. By careful study of the distribution of the ratio between shear stress and strain rate, we show that the order of magnitude of viscosity within the epithelial tissue is 100 Pa-hr and that increasing (decreasing) the actomyosin cytoskeleton and cell-cell adhesions increase (decrease) the magnitude of tissue viscosity. These results establish tissue viscosity as a meaningful way to describe the mechanical behavior of epithelial tissues, and demonstrate a direct relationship between tissue microstructure and material properties. By providing the first experimental measurement of tissue viscosity, our study is
Cell migration plays a fundamental role in numerous physiological processes, including embryonic development, wound healing, and cancer metastasis. While cell-cell adhesion is known to regulate motion by shaping cell morphology and intercellular force balance, its dynamic, rate-dependent contributions to tissue behavior remain poorly understood. In this study, we examine how the dissipative nature of cell-cell adhesion influences tissue dynamics and collective migration using an extended vertex model with explicit junctional viscosity. Our findings reveal a nontrivial interplay between two distinct components of adhesion: an interfacial adhesion energy (energetic, rate-independent) contribution, which sets the effective junctional tension, and a dissipative (rate-dependent) contribution, which controls resistance to relative motion during cell rearrangements. We show that increasing the energetic component promotes migration by modifying cell shape and lowering the barrier to neighbor exchanges, whereas strengthening the dissipative component induces jamming and suppresses cell motion. Linear rheological analysis further demonstrates that, in the unjammed regime, vertex-model tissu
Local stresses in a tissue, a collective property, regulate cell division and apoptosis. In turn, cell growth and division induce active stresses in the tissue. As a consequence, there is a feedback between cell growth and local stresses. However, how the cell dynamics depend on local stress-dependent cell division and the feedback strength is not fully understood. Here, we probe the consequences of stress-mediated growth and cell division on cell dynamics using agent-based simulations of a two-dimensional growing tissue. We discover a rich dynamical behavior of individual cells, ranging from jamming (mean square displacement, $Δ(t) \sim t^α$ with $α$ less than unity), to hyperdiffusion ($α> 2$) depending on cell division rate and the strength of the mechanical feedback. Strikingly, $Δ(t)$ is determined by the tissue growth law, which quantifies cell proliferation (number of cells $N(t)$ as a function of time). The growth law ($N(t) \sim t^λ$ at long times) is regulated by the critical pressure that controls the strength of the mechanical feedback and the ratio between cell division-apoptosis rates. We show that $λ\sim α$, which implies that higher growth rate leads to a greater
The application of single-cell molecular profiling coupled with spatial technologies has enabled charting cellular heterogeneity in reference tissues and in disease. This new wave of molecular data has highlighted the expected diversity of single-cell dynamics upon shared external queues and spatial organizations. However, little is known about the relationship between single cell heterogeneity and the emergence and maintenance of robust multicellular processes in developed tissues and its role in (patho)physiology. Here, we present emerging computational modeling strategies that use increasingly available large-scale cross-condition single cell and spatial datasets, to study multicellular organization in tissues and complement cell taxonomies. This perspective should enable us to better understand how cells within tissues collectively process information and adapt synchronized responses in disease contexts and to bridge the gap between structural changes and functions in tissues.
Cell detection is a fundamental task in computational pathology that can be used for extracting high-level medical information from whole-slide images. For accurate cell detection, pathologists often zoom out to understand the tissue-level structures and zoom in to classify cells based on their morphology and the surrounding context. However, there is a lack of efforts to reflect such behaviors by pathologists in the cell detection models, mainly due to the lack of datasets containing both cell and tissue annotations with overlapping regions. To overcome this limitation, we propose and publicly release OCELOT, a dataset purposely dedicated to the study of cell-tissue relationships for cell detection in histopathology. OCELOT provides overlapping cell and tissue annotations on images acquired from multiple organs. Within this setting, we also propose multi-task learning approaches that benefit from learning both cell and tissue tasks simultaneously. When compared against a model trained only for the cell detection task, our proposed approaches improve cell detection performance on 3 datasets: proposed OCELOT, public TIGER, and internal CARP datasets. On the OCELOT test set in partic
The spatiotemporal coordination and regulation of cell proliferation is fundamental in many aspects of development and tissue maintenance. Cells have the ability to adapt their division rates in response to mechanical constraints, yet we do not fully understand how cell proliferation regulation impacts cell migration phenomena. Here, we present a minimal continuum model of cell migration with cell cycle dynamics, which includes density-dependent effects and hence can account for cell proliferation regulation. By combining minimal mathematical modelling, Bayesian inference, and recent experimental data, we quantify the impact of tissue crowding across different cell cycle stages in epithelial tissue expansion experiments. Our model suggests that cells sense local density and adapt cell cycle progression in response, during G1 and the combined S/G2/M phases, providing an explicit relationship between each cell cycle stage duration and local tissue density, which is consistent with several experimental observations. Finally, we compare our mathematical model predictions to different experiments studying cell cycle regulation and present a quantitative analysis on the impact of density
Multicellular tissues are the building blocks of many biological systems and organs. These tissues are not static, but dynamically change over time. Even if the overall structure remains the same there is a turnover of cells within the tissue. This dynamic homeostasis is maintained by numerous governing mechanisms which are finely tuned in such a way that the tissue remains in a homeostatic state, even across large timescales. Some of these governing mechanisms include cell motion, and cell fate selection through inter cellular signalling. However, it is not yet clear how to link these two processes, or how they may affect one another across the tissue. In this paper, we present a multicellular, multiscale model, which brings together the two phenomena of cell motility, and inter cellular signalling, to describe cell fate selection on a dynamic tissue. We find that the affinity for cellular signalling to occur greatly influences a cells ability to differentiate. We also find that our results support claims that cell differentiation is a finely tuned process within dynamic tissues at homeostasis, with excessive cell turnover rates leading to unhealthy (undifferentiated and unpattern
Tissue mechanical properties such as rigidity and fluidity, and changes in these properties driven by jamming-unjamming transitions (UJT), have come under recent highlight as mechanical markers of health and disease in various biological processes including cancer. However, most analysis of these mechanical properties and UJT have sidestepped the effect of cellular death and division in these systems. Cellular apoptosis (programmed cell death) and mitosis (cell division) can drive significant changes in tissue properties. The balance between the two is crucial in maintaining tissue function, and an imbalance between the two is seen in situations such as cancer progression, wound healing and necrosis. In this work we investigate the impact of cell death and division on tissue mechanical properties, by incorporating specific mechanosensitive triggers of cell death and division based on the size and geometry of the cell within in silico models of tissue dynamics. Specifically, we look at cell migration, tissue response to external stress, tissue extrusion propensity and self-organization of different cell types within the tissue, as a function of cell death and division and the rules
Coordinated and cooperative motion of cells is essential for embryonic development, tissue morphogenesis, wound healing and cancer invasion. A predictive understanding of the emergent mechanical behaviors in collective cell motion is challenging due to the complex interplay between cell-cell interactions, cell-matrix adhesions and active cell behaviors. To overcome this challenge, we develop a predictive cellular vertex model that can delineate the relative roles of substrate rigidity, tissue mechanics and active cell properties on the movement of cell collectives. We apply the model to the specific case of collective motion in cell aggregates as they spread into a two-dimensional cell monolayer adherent to a soft elastic matrix. Consistent with recent experiments, we find that substrate stiffness regulates the driving forces for the spreading of cellular monolayer, which can be pressure-driven or crawling-based depending on substrate rigidity. On soft substrates, cell monolayer spreading is driven by an active pressure due to the influx of cells coming from the aggregate, whereas on stiff substrates, cell spreading is driven primarily by active crawling forces. Our model predicts
Formulating quantitative and predictive models for tissue development requires consideration of the complex, stochastic gene expression dynamics, its regulation via cell-to-cell interactions, and cell proliferation. Including all of these processes into a practical mathematical framework requires complex expressions that are difficult to interpret and apply. We construct a simple theory that incorporates intracellular stochastic gene expression dynamics, signaling chemicals that influence these dynamics and mediate cell-cell interactions, and cell proliferation and its accompanying differentiation. Cellular states (genetic and epigenetic) are described by a Waddington vector field that allows for non-gradient dynamics (cycles, entropy production, loss of detailed balance) which is precluded in Waddington potential landscape representations of gene expression dynamics. We define an epigenetic fitness landscape that describes the proliferation of different cell types, and elucidate how this fitness landscape is related to Waddington's vector field. We illustrate the applicability of our framework by analyzing two model systems: an interacting two-gene differentiation process and a sp
Cell-based, mathematical modeling of collective cell behavior has become a prominent tool in developmental biology. Cell-based models represent individual cells as single particles or as sets of interconnected particles, and predict the collective cell behavior that follows from a set of interaction rules. In particular, vertex-based models are a popular tool for studying the mechanics of confluent, epithelial cell layers. They represent the junctions between three (or sometimes more) cells in confluent tissues as point particles, connected using structural elements that represent the cell boundaries. A disadvantage of these models is that cell-cell interfaces are represented as straight lines. This is a suitable simplification for epithelial tissues, where the interfaces are typically under tension, but this simplification may not be appropriate for mesenchymal tissues or tissues that are under compression, such that the cell-cell boundaries can buckle. In this paper we introduce a variant of VMs in which this and two other limitations of VMs have been resolved. The new model can also be seen as on off-the-lattice generalization of the Cellular Potts Model. It is an extension of t
Multiplexed immuno-fluorescence tissue imaging, allowing simultaneous detection of molecular properties of cells, is an essential tool for characterizing the complex cellular mechanisms in translational research and clinical practice. New image analysis approaches are needed because tissue section stained with a mixture of protein, DNA and RNA biomarkers are introducing various complexities, including spurious edges due to fluorescent staining artifacts between touching or overlapping cells. We have developed the RRScell method harnessing the stochastic random-reaction-seed (RRS) algorithm and deep neural learning U-net to extract single-cell resolution profiling-map of gene expression over a million cells tissue section accurately and automatically. Furthermore, with the use of manifold learning technique UMAP for cell phenotype cluster analysis, the AI-driven RRScell has equipped with a marker-based image cytometry analysis tool (markerUMAP) in quantifying spatial distribution of cell phenotypes from tissue images with a mixture of biomarkers. The results achieved in this study suggest that RRScell provides a robust enough way for extracting cytometric single cell morphology as w
The growth of several biological tissues is known to be controlled in part by local geometrical features, such as the curvature of the tissue interface. This control leads to changes in tissue shape that in turn can affect the tissue's evolution. Understanding the cellular basis of this control is highly significant for bioscaffold tissue engineering, the evolution of bone microarchitecture, wound healing, and tumour growth. While previous models have proposed geometrical relationships between tissue growth and curvature, the role of cell density and cell vigor remains poorly understood. We propose a cell-based mathematical model of tissue growth to investigate the systematic influence of curvature on the collective crowding or spreading of tissue-synthesising cells induced by changes in local tissue surface area during the motion of the interface. Depending on the strength of diffusive damping, the model exhibits complex growth patterns such as undulating motion, efficient smoothing of irregularities, and the generation of cusps. We compare this model with in-vitro experiments of tissue deposition in bioscaffolds of different geometries. By accounting for the depletion of active c
It is known that the orientation of cellulose microfibrils within plant cell walls has an important impact on the morphogenesis of plant cells and tissues. Viewing the shape of a plant cell as a square prism or cylinder with the axis aligning with the primary direction of expansion and growth, the orientation of the microfibrils within the cell wall on the sides of the cell is known. However, not much is known about their orientation at the ends of the cell. Here we investigate the impact of the orientation of cellulose microfibrils within a plant cell wall at the ends of the cell by solving the equations of linear elasticity numerically. Three different scenarios for the orientation of the microfibrils are considered. The macroscopic elastic properties of the cell wall are obtained using homogenization theory from the microscopic description of the elastic properties of the cell wall microfibrils and wall matrix. It is found that the orientation of the microfibrils in the upper and lower parts of cell walls do not affect the expansion of the cell in the direction of its axis but do affect its expansion in the lateral directions. The arrangement of the microfibrils in the upper and
Cell-cell communication is essential for tissue development, regeneration and function, and its disruption can lead to diseases and developmental abnormalities. The revolution of single-cell genomics technologies offers unprecedented insights into cellular identities, opening new avenues to resolve the intricate cellular interactions present in tissue niches. CellPhoneDB is a bioinformatics toolkit designed to infer cell-cell communication by combining a curated repository of bona fide ligand-receptor interactions with a set of computational and statistical methods to integrate them with single-cell genomics data. Importantly, CellPhoneDB captures the multimeric nature of molecular complexes, thus representing cell-cell communication biology faithfully. Here we present CellPhoneDB v5, an updated version of the tool, which offers several new features. Firstly, the repository has been expanded by one-third with the addition of new interactions. These encompass interactions mediated by non-protein ligands such as endocrine hormones and GPCR ligands. Secondly, it includes a differentially expression-based methodology for more tailored interaction queries. Thirdly, it incorporates novel
During animal development and homeostasis, the structure of tissues, including muscles, blood vessels and connective tissues adapts to mechanical strains in the extracellular matrix (ECM). These strains originate from the differential growth of tissues or forces due to muscle contraction or gravity. Here we show using a computational model that by amplifying local strain cues, active cell contractility can facilitate and accelerate the reorientation of single cells to static strains. At the collective cell level, the model simulations show that active cell contractility can facilitate the formation of strings along the orientation of stretch. The computational model is based on a hybrid cellular Potts and finite-element simulation framework describing a mechanical cell-substrate feedback, where: 1) cells apply forces on the ECM, such that 2) local strains are generated in the ECM, and 3) cells preferentially extend protrusions along the strain orientation. In accordance with experimental observations, simulated cells align and form string-like structures parallel to static uniaxial stretch. Our model simulations predict that the magnitude of the uniaxial stretch and the strength of
Cell neighbor exchanges are integral to tissue rearrangements in biology, including development and repair. Often these processes occur via topological T1 transitions analogous to those observed in foams, grains and colloids. However, in contrast to in non-living materials the T1 transitions in biological tissues are rate-limited and cannot occur instantaneously due to the finite time required to remodel complex structures at cell-cell junctions. Here we study how this rate-limiting process affects the mechanics and collective behavior of cells in a tissue by introducing this important biological constraint in a theoretical vertex-based model as an intrinsic single-cell property. We report in the absence of this time constraint, the tissue undergoes a motility-driven glass transition characterized by a sharp increase in the intermittency of cell-cell rearrangements. Remarkably, this glass transition disappears as T1 transitions are temporally limited. As a unique consequence of limited rearrangements, we also find that the tissue develops spatially correlated streams of fast and slow cells, in which the fast cells organize into stream-like patterns with leader-follower interactions
Regulation of cell proliferation is a crucial aspect of tissue development and homeostasis and plays a major role in morphogenesis, wound healing, and tumor invasion. A phenomenon of such regulation is contact inhibition, which describes the dramatic slowing of proliferation, cell migration and individual cell growth when multiple cells are in contact with each other. While many physiological, molecular and genetic factors are known, the mechanism of contact inhibition is still not fully understood. In particular, the relevance of cellular signaling due to interfacial contact for contact inhibition is still debated. Cellular automata (CA) have been employed in the past as numerically efficient mathematical models to study the dynamics of cell ensembles, but they are not suitable to explore the origins of contact inhibition as such agent-based models assume fixed cell sizes. We develop a minimal, data-driven model to simulate the dynamics of planar cell cultures by extending a probabilistic CA to incorporate size changes of individual cells during growth and cell division. We successfully apply this model to previous in-vitro experiments on contact inhibition in epithelial tissue: A
In vivo and in vitro cells rely on the support of an underlying biocompatible substrate, such as the extracellular matrix or a culture substrate, to spread and proliferate. The mechanical and chemical properties of such structures play a central role in the dynamical and statistical properties of the tissue. At the cell scale, these substrates are highly disordered. Here, we investigate how spatial heterogeneities of the cell-substrate interaction influence the motility of the cells in a model confluent tissue. We use the Self-Propelled Voronoi model and describe the disorder as a spatially dependent preferred geometry of the individual cells. We found that when the characteristic length scale of the preferred geometry is smaller than the cell size, the tissue is less rigid than its homogeneous counterpart, with a consequent increase in cell motility. This result is in sharp contrast to what has been reported for tissues with heterogeneity in the mechanical properties of the individual cells, where the disorder favors rigidity. Using the fraction of rigid cells, we observe a collapse of the motility data for different model parameters and provide evidence that the rigidity transiti
Measurements on embryonic epithelial tissues in a diverse range of organisms have shown that the statistics of cell neighbor numbers are universal in tissues where cell proliferation is the primary cell activity. Highly simplified non-spatial models of proliferation are claimed to accurately reproduce these statistics. Using a systematic critical analysis, we show that non-spatial models are not capable of robustly describing the universal statistics observed in proliferating epithelia, indicating strong spatial correlations between cells. Furthermore we show that spatial simulations using the Subcellular Element Model are able to robustly reproduce the universal histogram. In addition these simulations are able to unify ostensibly divergent experimental data in the literature. We also analyze cell neighbor statistics in early stages of chick embryo development in which cell behaviors other than proliferation are important. We find from experimental observation that cell neighbor statistics in the primitive streak region, where cell motility and ingression are also important, show a much broader distribution. A non-spatial Markov process model provides excellent agreement with this