Three-dimensional vertex model for simulating multicellular morphogenesis

Effective mechanical potential of cell-cell interaction in tissues harboring cavity and in cell sheet toward morphogenesis

Measuring mechanical forces of cell-cell interactions is important for studying morphogenesis in multicellular organisms. We previously reported an image-based statistical method for inferring effective mechanical potentials of pairwise cell-cell interactions by fitting cell tracking data with a theoretical model. However, whether this method is applicable to tissues with non-cellular components such as cavities remains elusive. Here we evaluated the applicability of the method to cavity-harboring tissues. Using synthetic data generated by simulations, we found that the effect of expanding cavities was added to the pregiven potentials used in the simulations, resulting in the inferred effective potentials having an additional repulsive component derived from the expanding cavities. Interestingly, simulations by using the effective potentials reproduced the cavity-harboring structures. Then, we applied our method to the mouse blastocysts, and found that the inferred effective potentials can reproduce the cavity-harboring structures. Pairwise potentials with additional repulsive components were also detected in two-dimensional cell sheets, by which curved sheets including tubes and cups were simulated. We conclude that our inference method is applicable to tissues harboring cavities and cell sheets, and the resultant effective potentials are useful to simulate the morphologies.

Poissonian Cellular Potts Models Reveal Nonequilibrium Kinetics of Cell Sorting

Cellular Potts models are broadly applied across developmental biology and cancer research. We overcome limitations of the traditional approach, which reinterprets a modified Metropolis sampling as ad hoc dynamics, by introducing a physical timescale through Poissonian kinetics and by applying principles of stochastic thermodynamics to separate thermal and relaxation effects from athermal noise and nonconservative forces. Our method accurately describes cell-sorting dynamics in mouse-embryo development and identifies the distinct contributions of nonequilibrium processes, e.g., cell growth and active fluctuations.

Balancing competing effects of tissue growth and cytoskeletal regulation during Drosophila wing disc development

How a developing organ robustly coordinates the cellular mechanics and growth to reach a final size and shape remains poorly understood. Through iterations between experiments and model simulations that include a mechanistic description of interkinetic nuclear migration, we show that the local curvature, height, and nuclear positioning of cells in the Drosophila wing imaginal disc are defined by the concurrent patterning of actomyosin contractility, cell-ECM adhesion, ECM stiffness, and interfacial membrane tension. We show that increasing cell proliferation via different growth-promoting pathways results in two distinct phenotypes. Triggering proliferation through insulin signaling increases basal curvature, but an increase in growth through Dpp signaling and Myc causes tissue flattening. These distinct phenotypic outcomes arise from differences in how each growth pathway regulates the cellular cytoskeleton, including contractility and cell-ECM adhesion. The coupled regulation of proliferation and cytoskeletal regulators is a general strategy to meet the multiple context-dependent criteria defining tissue morphogenesis.

General, open-source vertex modeling in biological applications using Tissue Forge

Vertex models are a widespread approach for describing the biophysics and behaviors of multicellular systems, especially of epithelial tissues. Vertex models describe a wide variety of developmental scenarios and behaviors like cell rearrangement and tissue folding. Often, these models are implemented as single-use or closed-source software, which inhibits reproducibility and decreases accessibility for researchers with limited proficiency in software development and numerical methods. We developed a physics-based vertex model methodology in Tissue Forge, an open-source, particle-based modeling and simulation environment. Our methodology describes the properties and processes of vertex model objects on the basis of vertices, which allows integration of vertex modeling with the particle-based formalism of Tissue Forge, enabling an environment for developing mixed-method models of multicellular systems. Our methodology in Tissue Forge inherits all features provided by Tissue Forge, delivering open-source, extensible vertex modeling with interactive simulation, real-time simulation visualization and model sharing in the C, C++ and Python programming languages and a Jupyter Notebook. Demonstrations show a vertex model of cell sorting and a mixed-method model of cell migration combining vertex- and particle-based models. Our methodology provides accessible vertex modeling for a broad range of scientific disciplines, and we welcome community-developed contributions to our open-source software implementation.

Generative models of morphogenesis in developmental biology

Understanding the mechanism by which cells coordinate their differentiation and migration is critical to our understanding of many fundamental processes such as wound healing, disease progression, and developmental biology. Mathematical models have been an essential tool for testing and developing our understanding, such as models of cells as soft spherical particles, reaction-diffusion systems that couple cell movement to environmental factors, and multi-scale multi-physics simulations that combine bottom-up rule-based models with continuum laws. However, mathematical models can often be loosely related to data or have so many parameters that model behaviour is weakly constrained. Recent methods in machine learning introduce new means by which models can be derived and deployed. In this review, we discuss examples of mathematical models of aspects of developmental biology, such as cell migration, and how these models can be combined with these recent machine learning methods.

General, Open-Source Vertex Modeling in Biological Applications Using Tissue Forge

Vertex models are a widespread approach for describing the biophysics and behaviors of multicellular systems, especially of epithelial tissues. Vertex models describe a wide variety of developmental scenarios and behaviors like cell rearrangement and tissue folding. Often, these models are implemented as single-use or closed-source software, which inhibits reproducibility and decreases accessibility for researchers with limited proficiency in software development and numerical methods. We developed a physics-based vertex model methodology in Tissue Forge, an open-source, particle-based modeling and simulation environment. Our methodology describes the properties and processes of vertex model objects on the basis of vertices, which allows integration of vertex modeling with the particle-based formalism of Tissue Forge, enabling an environment for developing mixed-method models of multicellular systems. Our methodology in Tissue Forge inherits all features provided by Tissue Forge, delivering opensource, extensible vertex modeling with interactive simulation, real-time simulation visualization and model sharing in the C , C + + and Python programming languages and a Jupyter Notebook. Demonstrations show a vertex model of cell sorting and a mixed-method model of cell migration combining vertex- and particle-based models. Our methodology provides accessible vertex modeling for a broad range of scientific disciplines, and we welcome community-developed contributions to our open-source software implementation.

Precise alternating cellular pattern in the inner ear by coordinated hopping intercalations and delaminations

The mammalian hearing organ, the organ of Corti, is one of the most organized tissues in mammals. It contains a precisely positioned array of alternating sensory hair cells (HCs) and nonsensory supporting cells. How such precise alternating patterns emerge during embryonic development is not well understood. Here, we combine live imaging of mouse inner ear explants with hybrid mechano-regulatory models to identify the processes that underlie the formation of a single row of inner hair cells (IHCs). First, we identify a previously unobserved morphological transition, termed "hopping intercalation," that allows cells differentiating toward IHC fate to "hop" under the apical plane into their final position. Second, we show that out-of-row cells with low levels of the HC marker Atoh1 delaminate. Last, we show that differential adhesion between cell types contributes to straightening of the IHC row. Our results support a mechanism for precise patterning based on coordination between signaling and mechanical forces that is likely relevant for many developmental processes.

Modelling of Tissue Invasion in Epithelial Monolayers

Mathematical and computational models are used to describe biomechanical processes in multicellular systems. Here, we develop a model to analyse how two types of epithelial cell layers interact during tissue invasion depending on their cellular properties, i.e., simulating cancer cells expanding into a region of normal cells. We model the tissue invasion process using the cellular Potts model and implement our two-dimensional computational simulations in the software package CompuCell3D. The model predicts that differences in mechanical properties of cells can lead to tissue invasion, even if the division rates and death rates of the two cell types are the same. We also show how the invasion speed varies depending on the cell division and death rates and the mechanical properties of the cells.

Novel Ground-Up 3D Multicellular Simulators for Synthetic Biology CAD Integrating Stochastic Gillespie Simulations Benchmarked with Topologically Variable SBML Models

The elevation of Synthetic Biology from single cells to multicellular simulations would be a significant scale-up. The spatiotemporal behavior of cellular populations has the potential to be prototyped in silico for computer assisted design through ergonomic interfaces. Such a platform would have great practical potential across medicine, industry, research, education and accessible archiving in bioinformatics. Existing Synthetic Biology CAD systems are considered limited regarding population level behavior, and this work explored the in silico challenges posed from biological and computational perspectives. Retaining the connection to Synthetic Biology CAD, an extension of the Infobiotics Workbench Suite was considered, with potential for the integration of genetic regulatory models and/or chemical reaction networks through Next Generation Stochastic Simulator (NGSS) Gillespie algorithms. These were executed using SBML models generated by in-house SBML-Constructor over numerous topologies and benchmarked in association with multicellular simulation layers. Regarding multicellularity, two ground-up multicellular solutions were developed, including the use of Unreal Engine 4 contrasted with CPU multithreading and Blender visualization, resulting in a comparison of real-time versus batch-processed simulations. In conclusion, high-performance computing and client-server architectures could be considered for future works, along with the inclusion of numerous biologically and physically informed features, whilst still pursuing ergonomic solutions.

Organoid Cultures In Silico: Tools or Toys?

The implementation of stem-cell-based organoid culture more than ten years ago started a development that created new avenues for diagnostic analyses and regenerative medicine. In parallel, computational modelling groups realized the potential of this culture system to support their theoretical approaches to study tissues in silico. These groups developed computational organoid models (COMs) that enabled testing consistency between cell biological data and developing theories of tissue self-organization. The models supported a mechanistic understanding of organoid growth and maturation and helped linking cell mechanics and tissue shape in general. What comes next? Can we use COMs as tools to complement the equipment of our biological and medical research? While these models already support experimental design, can they also quantitatively predict tissue behavior? Here, we review the current state of the art of COMs and discuss perspectives for their application.

Computational approaches for simulating luminogenesis

Lumens, liquid-filled cavities surrounded by polarized tissue cells, are elementary units involved in the morphogenesis of organs. Theoretical modeling and computations, which can integrate various factors involved in biophysics of morphogenesis of cell assembly and lumens, may play significant roles to elucidate the mechanisms in formation of such complex tissue with lumens. However, up to present, it has not been documented well what computational approaches or frameworks can be applied for this purpose and how we can choose the appropriate approach for each problem. In this review, we report some typical lumen morphologies and basic mechanisms for the development of lumens, focusing on three keywords - mechanics, hydraulics and geometry - while outlining pros and cons of the current main computational strategies. We also describe brief guidance of readouts, i.e., what we should measure in experiments to make the comparison with the model's assumptions and predictions.

Appropriate tension sensitivity of α-catenin ensures rounding morphogenesis of epithelial spheroids

The adherens junction (AJ) is an actin filament-anchoring junction. It plays a central role in epithelial morphogenesis through cadherin-based recognition and adhesion among cells. The stability and plasticity of AJs are required for the morphogenesis. An actin-binding α-catenin is an essential component of the cadherin-catenin complex and functions as a tension transducer that changes its conformation and induces AJ development in response to tension. Despite much progress in understanding molecular mechanisms of tension sensitivity of α-catenin, its significance on epithelial morphogenesis is still unknown. Here we show that the tension sensitivity of α-catenin is essential for epithelial cells to form round spheroids through proper multicellular rearrangement. Using a novel in vitro suspension culture model, we found that epithelial cells form round spheroids even from rectangular-shaped cell masses with high aspect ratios without using high tension and that increased tension sensitivity of α-catenin affected this morphogenesis. Analyses of AJ formation and cellular tracking during rounding morphogenesis showed cellular rearrangement, probably through AJ remodeling. The rearrangement occurs at the cell mass level, but not single-cell level. Hypersensitive α-catenin mutant-expressing cells did not show cellular rearrangement at the cell mass level, suggesting that the appropriate tension sensitivity of α-catenin is crucial for the coordinated round morphogenesis.Key words: α-catenin, vinculin, adherens junction, morphogenesis, mechanotransduction.

Arrested coalescence of multicellular aggregates

Multicellular aggregates are known to exhibit liquid-like properties. The fusion process of two cell aggregates is commonly studied as the coalescence of two viscous drops. However, tissues are complex materials and can exhibit viscoelastic behaviour. It is known that elastic effects can prevent the complete fusion of two drops, a phenomenon known as arrested coalescence. Here we study this phenomenon in stem cell aggregates and provide a theoretical framework which agrees with the experiments. In addition, agent-based simulations show that active cell fluctuations can control a solid-to-fluid phase transition, revealing that arrested coalescence can be found in the vicinity of an unjamming transition. By analysing the dynamics of the fusion process and combining it with nanoindentation measurements, we obtain the effective viscosity, shear modulus and surface tension of the aggregates. More generally, our work provides a simple, fast and inexpensive method to characterize the mechanical properties of viscoelastic materials.

Seven challenges in the multiscale modeling of multicellular tissues

The growth and dynamics of multicellular tissues involve tightly regulated and coordinated morphogenetic cell behaviors, such as shape changes, movement, and division, which are governed by subcellular machinery and involve coupling through short- and long-range signals. A key challenge in the fields of developmental biology, tissue engineering and regenerative medicine is to understand how relationships between scales produce emergent tissue-scale behaviors. Recent advances in molecular biology, live-imaging and ex vivo techniques have revolutionized our ability to study these processes experimentally. To fully leverage these techniques and obtain a more comprehensive understanding of the causal relationships underlying tissue dynamics, computational modeling approaches are increasingly spanning multiple spatial and temporal scales, and are coupling cell shape, growth, mechanics, and signaling. Yet such models remain challenging: modeling at each scale requires different areas of technical skills, while integration across scales necessitates the solution to novel mathematical and computational problems. This review aims to summarize recent progress in multiscale modeling of multicellular tissues and to highlight ongoing challenges associated with the construction, implementation, interrogation, and validation of such models. This article is categorized under: Reproductive System Diseases > Computational Models Metabolic Diseases > Computational Models Cancer > Computational Models.

Extracellular hyaluronate pressure shaped by cellular tethers drives tissue morphogenesis

How tissues acquire complex shapes is a fundamental question in biology and regenerative medicine. Zebrafish semicircular canals form from invaginations in the otic epithelium (buds) that extend and fuse to form the hubs of each canal. We find that conventional actomyosin-driven behaviors are not required. Instead, local secretion of hyaluronan, made by the enzymes uridine 5'-diphosphate dehydrogenase (ugdh) and hyaluronan synthase 3 (has3), drives canal morphogenesis. Charged hyaluronate polymers osmotically swell with water and generate isotropic extracellular pressure to deform the overlying epithelium into buds. The mechanical anisotropy needed to shape buds into tubes is conferred by a polarized distribution of actomyosin and E-cadherin-rich membrane tethers, which we term cytocinches. Most work on tissue morphogenesis ascribes actomyosin contractility as the driving force, while the extracellular matrix shapes tissues through differential stiffness. Our work inverts this expectation. Hyaluronate pressure shaped by anisotropic tissue stiffness may be a widespread mechanism for powering morphological change in organogenesis and tissue engineering.

A rigid body framework for multicellular modeling

Off-lattice models are a well-established approach in multicellular modeling, where cells are represented as points that are free to move in space. The representation of cells as point objects is useful in a wide range of settings, particularly when large populations are involved; however, a purely point-based representation is not naturally equipped to deal with objects that have length, such as cell boundaries or external membranes. Here we introduce an off-lattice modeling framework that exploits rigid body mechanics to represent objects using a collection of conjoined one-dimensional edges in a viscosity-dominated system. This framework can be used to represent cells as free moving polygons, to allow epithelial layers to smoothly interact with themselves, to model rod-shaped cells such as bacteria and to robustly represent membranes. We demonstrate that this approach offers solutions to the problems that limit the scope of current off-lattice multicellular models.

Dissecting Organismal Morphogenesis by Bridging Genetics and Biophysics

Multicellular organisms develop complex shapes from much simpler, single-celled zygotes through a process commonly called morphogenesis. Morphogenesis involves an interplay between several factors, ranging from the gene regulatory networks determining cell fate and differentiation to the mechanical processes underlying cell and tissue shape changes. Thus, the study of morphogenesis has historically been based on multidisciplinary approaches at the interface of biology with physics and mathematics. Recent technological advances have further improved our ability to study morphogenesis by bridging the gap between the genetic and biophysical factors through the development of new tools for visualizing, analyzing, and perturbing these factors and their biochemical intermediaries. Here, we review how a combination of genetic, microscopic, biophysical, and biochemical approaches has aided our attempts to understand morphogenesis and discuss potential approaches that may be beneficial to such an inquiry in the future. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Cell fate coordinates mechano-osmotic forces in intestinal crypt formation

Intestinal organoids derived from single cells undergo complex crypt-villus patterning and morphogenesis. However, the nature and coordination of the underlying forces remains poorly characterized. Here, using light-sheet microscopy and large-scale imaging quantification, we demonstrate that crypt formation coincides with a stark reduction in lumen volume. We develop a 3D biophysical model to computationally screen different mechanical scenarios of crypt morphogenesis. Combining this with live-imaging data and multiple mechanical perturbations, we show that actomyosin-driven crypt apical contraction and villus basal tension work synergistically with lumen volume reduction to drive crypt morphogenesis, and demonstrate the existence of a critical point in differential tensions above which crypt morphology becomes robust to volume changes. Finally, we identified a sodium/glucose cotransporter that is specific to differentiated enterocytes that modulates lumen volume reduction through cell swelling in the villus region. Together, our study uncovers the cellular basis of how cell fate modulates osmotic and actomyosin forces to coordinate robust morphogenesis.

A hybrid model of intercellular tension and cell-matrix mechanical interactions in a multicellular geometry

Epithelial cells form continuous sheets of cells that exist in tensional homeostasis. Homeostasis is maintained through cell-to-cell junctions that distribute tension and balance forces between cells and their underlying matrix. Disruption of tensional homeostasis can lead to epithelial-mesenchymal transition (EMT), a transdifferentiation process in which epithelial cells adopt a mesenchymal phenotype, losing cell-cell adhesion and enhancing cellular motility. This process is critical during embryogenesis and wound healing, but is also dysregulated in many disease states. To further understand the role of intercellular tension in spatial patterning of epithelial cell monolayers, we developed a multicellular computational model of cell-cell and cell-substrate forces. This work builds on a hybrid cellular Potts model (CPM)-finite element model to evaluate cell-matrix mechanical feedback of an adherent multicellular cluster. Cellular movement is governed by thermodynamic constraints from cell volume, cell-cell and cell-matrix contacts, and durotaxis, which arises from cell-generated traction forces on a finite element substrate. Junction forces at cell-cell contacts balance these traction forces, thereby producing a mechanically stable epithelial monolayer. Simulations were compared to in vitro experiments using fluorescence-based junction force sensors in clusters of cells undergoing EMT. Results indicate that the multicellular CPM model can reproduce many aspects of EMT, including epithelial monolayer formation dynamics, changes in cell geometry, and spatial patterning of cell-cell forces in an epithelial tissue.

Mechanochemical Coupling and Junctional Forces during Collective Cell Migration

Cell migration, a fundamental physiological process in which cells sense and move through their surrounding physical environment, plays a critical role in development and tissue formation, as well as pathological processes, such as cancer metastasis and wound healing. During cell migration, dynamics are governed by the bidirectional interplay between cell-generated mechanical forces and the activity of Rho GTPases, a family of small GTP-binding proteins that regulate actin cytoskeleton assembly and cellular contractility. These interactions are inherently more complex during the collective migration of mechanically coupled cells because of the additional regulation of cell-cell junctional forces. In this study, we adapted a recent minimal modeling framework to simulate the interactions between mechanochemical signaling in individual cells and interactions with cell-cell junctional forces during collective cell migration. We find that migration of individual cells depends on the feedback between mechanical tension and Rho GTPase activity in a biphasic manner. During collective cell migration, waves of Rho GTPase activity mediate mechanical contraction/extension and thus synchronization throughout the tissue. Further, cell-cell junctional forces exhibit distinct spatial patterns during collective cell migration, with larger forces near the leading edge. Larger junctional force magnitudes are associated with faster collective cell migration and larger tissue size. Simulations of heterogeneous tissue migration exhibit a complex dependence on the properties of both leading and trailing cells. Computational predictions demonstrate that collective cell migration depends on both the emergent dynamics and interactions between cellular-level Rho GTPase activity and contractility and multicellular-level junctional forces.

Balance of mechanical forces drives endothelial gap formation and may facilitate cancer and immune-cell extravasation

The formation of gaps in the endothelium is a crucial process underlying both cancer and immune cell extravasation, contributing to the functioning of the immune system during infection, the unfavorable development of chronic inflammation and tumor metastasis. Here, we present a stochastic-mechanical multiscale model of an endothelial cell monolayer and show that the dynamic nature of the endothelium leads to spontaneous gap formation, even without intervention from the transmigrating cells. These gaps preferentially appear at the vertices between three endothelial cells, as opposed to the border between two cells. We quantify the frequency and lifetime of these gaps, and validate our predictions experimentally. Interestingly, we find experimentally that cancer cells also preferentially extravasate at vertices, even when they first arrest on borders. This suggests that extravasating cells, rather than initially signaling to the endothelium, might exploit the autonomously forming gaps in the endothelium to initiate transmigration.

Mechanics of epithelial tissue formation

A key process in the life of any multicellular organism is its development from a single egg into a full grown adult. The first step in this process often consists of forming a tissue layer out of randomly placed cells on the surface of the egg. We present a model for generating such a tissue, based on mechanical interactions between the cells, and find that the resulting cellular pattern corresponds to the Voronoi tessellation of the nuclei of the cells. Experimentally, we obtain the same result in both fruit flies and flour beetles, with a distribution of cell shapes that matches that of the model, without any adjustable parameters. Finally, we show that this pattern is broken when the cells grow at different rates.

Topological graph description of multicellular dynamics based on vertex model

Vertex models are generally powerful tools for exploring biological insights into multicellular dynamics. In these models, a multicellular structure is represented by a network, which is dynamically rearranged using topological operations. Remarkably, the topological dynamics of the network are important in guaranteeing the results from the models and their biological implications. However, it remains unclear whether the entire pattern of multicellular topological dynamics can be accurately expressed by a set of operators in the models. Surprisingly, vertex models have been empirically used for several decades without any mathematical verification. In this study, we propose a rigorous two-/three-dimensional (2D/3D) vertex model to describe multicellular topological dynamics. To do this, we classify several types of vertex models from a graph-theoretic perspective. Based on the classification, mathematical analyses reveal several conditions that enable us to apply the operators accurately without topological errors. Under these conditions, the operators can completely express the entire pattern of multicellular topological dynamics. From these results, we newly propose rigorous 2D/3D vertex models that can be applied to general multicellular dynamics, and we clarify several points to verify the results obtained from previous models.

Turing Instability-Driven Biofabrication of Branching Tissue Structures: A Dynamic Simulation and Analysis Based on the Reaction⁻Diffusion Mechanism †

Four-dimensional (4D) biofabrication techniques aim to dynamically produce and control three-dimensional (3D) biological structures that would transform their shapes or functionalities with time, when a stimulus is imposed or cell post-printing self-assembly occurs. The evolution of 3D branching patterns via self-assembly of cells is critical for the 4D biofabrication of artificial organs or tissues with branched geometry. However, it is still unclear how the formation and evolution of these branching patterns are biologically encoded. Here, we study the biofabrication of lung branching structures utilizing a simulation model based on Turing instability that raises a dynamic reaction⁻diffusion (RD) process of the biomolecules and cells. The simulation model incorporates partial differential equations of four variables, describing the tempo-spatial distribution of the variables in 3D over time. The simulation results present the formation and evolution process of 3D branching patterns over time and also interpret both the behaviors of side-branching and tip-splitting as the stalk grows and the fabrication style under an external concentration gradient of morphogen, through 3D visualization. This provides a theoretical framework for rationally guiding the 4D biofabrication of lung airway grafts via cellular self-organization, which would potentially reduce the complexity of future experimental research and number of trials.

Combining Turing and 3D vertex models reproduces autonomous multicellular morphogenesis with undulation, tubulation, and branching

This study demonstrates computational simulations of multicellular deformation coupled with chemical patterning in the three-dimensional (3D) space. To address these aspects, we proposes a novel mathematical model, where a reaction–diffusion system is discretely expressed at a single cell level and combined with a 3D vertex model. To investigate complex phenomena emerging from the coupling of patterning and deformation, as an example, we employed an activator–inhibitor system and converted the activator concentration of individual cells into their growth rate. Despite the simplicity of the model, by growing a monolayer cell vesicle, the coupling system provided rich morphological dynamics such as undulation, tubulation, and branching. Interestingly, the morphological variety depends on the difference in time scales between patterning and deformation, and can be partially understood by the intrinsic hysteresis in the activator-inhibitor system with domain growth. Importantly, the model can be applied to 3D multicellular dynamics that couple the reaction–diffusion patterning with various cell behaviors, such as deformation, rearrangement, division, apoptosis, differentiation, and proliferation. Thus, the results demonstrate the significant advantage of the proposed model as well as the biophysical importance of exploring spatiotemporal dynamics of the coupling phenomena of patterning and deformation in 3D space.

Linking stem cell function and growth pattern of intestinal organoids

Intestinal stem cells (ISCs) require well-defined signals from their environment in order to carry out their specific functions. Most of these signals are provided by neighboring cells that form a stem cell niche, whose shape and cellular composition self-organize. Major features of this self-organization can be studied in ISC-derived organoid culture. In this system, manipulation of essential pathways of stem cell maintenance and differentiation results in well-described growth phenotypes. We here provide an individual cell-based model of intestinal organoids that enables a mechanistic explanation of the observed growth phenotypes. In simulation studies of the 3D structure of expanding organoids, we investigate interdependences between Wnt- and Notch-signaling which control the shape of the stem cell niche and, thus, the growth pattern of the organoids. Similar to in vitro experiments, changes of pathway activities alter the cellular composition of the organoids and, thereby, affect their shape. Exogenous Wnt enforces transitions from branched into a cyst-like growth pattern; known to occur spontaneously during long term organoid expansion. Based on our simulation results, we predict that the cyst-like pattern is associated with biomechanical changes of the cells which assign them a growth advantage. The results suggest ongoing stem cell adaptation to in vitro conditions during long term expansion by stabilizing Wnt-activity. Our study exemplifies the potential of individual cell-based modeling in unraveling links between molecular stem cell regulation and 3D growth of tissues. This kind of modeling combines experimental results in the fields of stem cell biology and cell biomechanics constituting a prerequisite for a better understanding of tissue regeneration as well as developmental processes.

Epithelial Folding Driven by Apical or Basal-Lateral Modulation: Geometric Features, Mechanical Inference, and Boundary Effects

During embryonic development, epithelial sheets fold into complex structures required for tissue and organ functions. Although substantial efforts have been devoted to identifying molecular mechanisms underlying epithelial folding, far less is understood about how forces deform individual cells to sculpt the overall sheet morphology. Here we describe a simple and general theoretical model for the autonomous folding of monolayered epithelial sheets. We show that active modulation of intracellular mechanics along the basal-lateral as well as the apical surfaces is capable of inducing fold formation in the absence of buckling instability. Apical modulation sculpts epithelia into shallow and V-shaped folds, whereas basal-lateral modulation generates deep and U-shaped folds. These characteristic tissue shapes remain unchanged when subject to mechanical perturbations from the surroundings, illustrating that the autonomous folding is robust against environmental variabilities. At the cellular scale, how cells change shape depends on their initial aspect ratios and the modulation mechanisms. Such cell deformation characteristics are verified via experimental measurements for a canonical folding process driven by apical modulation, indicating that our theory could be used to infer the underlying folding mechanisms based on experimental data. The mechanical principles revealed in our model could potentially guide future studies on epithelial folding in diverse systems.

Contractile actin belt and mesh structures provide the opposite dependence of epithelial stiffness on the spontaneous curvature of constituent cells

Actomyosin generates contractile forces within cells, which have a crucial role in determining the macroscopic mechanical properties of epithelial tissues. Importantly, actin cytoskeleton, which propagates actomyosin contractile forces, forms several characteristic structures in a 3D intracellular space, such as a circumferential actin belt lining adherence junctions and an actin mesh beneath the apical membrane. However, little is known about how epithelial mechanical property depends on the intracellular contractile structures. We performed computational simulations using a 3D vertex model, and demonstrated the longitudinal tensile test of an epithelial tube, whose inside and outside are defined as the apical and basal surfaces, respectively. As a result, these subcellular structures provide the contrary dependence of epithelial stiffness and fracture force on the spontaneous curvature of constituent cells; the epithelial stiffness increases with increasing the spontaneous curvature in the case of belt, meanwhile it decreases in the case of mesh. This qualitative difference emerges from the different anisotropic deformability of apical cell surfaces; while belt preserves isotropic apical cell shapes, mesh does not. Moreover, the difference in the anisotropic deformability determines the frequency of cell rearrangements, which in turn effectively decrease the tube stiffness. These results illustrate the importance of the intracellular contractile structures, which may be regulated to optimize mechanical functions of individual epithelial tissues.

Active Vertex Model for cell-resolution description of epithelial tissue mechanics

We introduce an Active Vertex Model (AVM) for cell-resolution studies of the mechanics of confluent epithelial tissues consisting of tens of thousands of cells, with a level of detail inaccessible to similar methods. The AVM combines the Vertex Model for confluent epithelial tissues with active matter dynamics. This introduces a natural description of the cell motion and accounts for motion patterns observed on multiple scales. Furthermore, cell contacts are generated dynamically from positions of cell centres. This not only enables efficient numerical implementation, but provides a natural description of the T1 transition events responsible for local tissue rearrangements. The AVM also includes cell alignment, cell-specific mechanical properties, cell growth, division and apoptosis. In addition, the AVM introduces a flexible, dynamically changing boundary of the epithelial sheet allowing for studies of phenomena such as the fingering instability or wound healing. We illustrate these capabilities with a number of case studies.

Role of molecular turnover in dynamic deformation of a three-dimensional cellular membrane

In cells, the molecular constituents of membranes are dynamically turned over by transportation from one membrane to another. This molecular turnover causes the membrane to shrink or expand by sensing the stress state within the cell, changing its morphology. At present, little is known as to how this turnover regulates the dynamic deformation of cellular membranes. In this study, we propose a new physical model by which molecular turnover is coupled with three-dimensional membrane deformation to explore mechanosensing roles of turnover in cellular membrane deformations. In particular, as an example of microscopic machinery, based on a coarse-graining description, we suppose that molecular turnover depends on the local membrane strain. Using the proposed model, we demonstrate computational simulations of a single vesicle. The results show that molecular turnover adaptively facilitates vesicle deformation, owing to its stress dependence; while the vesicle drastically expands in the case with low bending rigidity, it shrinks in that with high bending rigidity. Moreover, localized active tension on the membrane causes cellular migration by driving the directional transport of molecules within the cell. These results illustrate the use of the proposed model as well as the role of turnover in the dynamic deformations of cellular membranes.

Multi-scale computational study of the mechanical regulation of cell mitotic rounding in epithelia

Mitotic rounding during cell division is critical for preventing daughter cells from inheriting an abnormal number of chromosomes, a condition that occurs frequently in cancer cells. Cells must significantly expand their apical area and transition from a polygonal to circular apical shape to achieve robust mitotic rounding in epithelial tissues, which is where most cancers initiate. However, how cells mechanically regulate robust mitotic rounding within packed tissues is unknown. Here, we analyze mitotic rounding using a newly developed multi-scale subcellular element computational model that is calibrated using experimental data. Novel biologically relevant features of the model include separate representations of the sub-cellular components including the apical membrane and cytoplasm of the cell at the tissue scale level as well as detailed description of cell properties during mitotic rounding. Regression analysis of predictive model simulation results reveals the relative contributions of osmotic pressure, cell-cell adhesion and cortical stiffness to mitotic rounding. Mitotic area expansion is largely driven by regulation of cytoplasmic pressure. Surprisingly, mitotic shape roundness within physiological ranges is most sensitive to variation in cell-cell adhesivity and stiffness. An understanding of how perturbed mechanical properties impact mitotic rounding has important potential implications on, amongst others, how tumors progressively become more genetically unstable due to increased chromosomal aneuploidy and more aggressive.