Skip to main content
Advertisement

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

  • Ali Nematbakhsh ,

    Contributed equally to this work with: Ali Nematbakhsh, Wenzhao Sun, Pavel A. Brodskiy, Aboutaleb Amiri

    Affiliations Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, Indiana, United States of America, Department of Mathematics, University of California, Riverside, California, United States of America

  • Wenzhao Sun ,

    Contributed equally to this work with: Ali Nematbakhsh, Wenzhao Sun, Pavel A. Brodskiy, Aboutaleb Amiri

    Affiliation Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, Indiana, United States of America

  • Pavel A. Brodskiy ,

    Contributed equally to this work with: Ali Nematbakhsh, Wenzhao Sun, Pavel A. Brodskiy, Aboutaleb Amiri

    Affiliation Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana, United States of America

  • Aboutaleb Amiri ,

    Contributed equally to this work with: Ali Nematbakhsh, Wenzhao Sun, Pavel A. Brodskiy, Aboutaleb Amiri

    Affiliation Department of Physics, University of Notre Dame, Notre Dame, Indiana, United States of America

  • Cody Narciso,

    Affiliation Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana, United States of America

  • Zhiliang Xu,

    Affiliation Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, Indiana, United States of America

  • Jeremiah J. Zartman ,

    malber@nd.edu (MA); jzartman@nd.edu (JJZ)

    Affiliation Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana, United States of America

  • Mark Alber

    malber@nd.edu (MA); jzartman@nd.edu (JJZ)

    Affiliations Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, Indiana, United States of America, Department of Mathematics, University of California, Riverside, California, United States of America

Abstract

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.

Author summary

Mitotic rounding (MR) during cell division which is critical for the robust segregation of chromosomes into daughter cells, plays important roles in tissue growth and morphogenesis, and is frequently perturbed in cancerous cells. Mechanisms of MR have been investigated in individual cultured cells, but mechanisms regulating MR in tissues are still poorly understood. We developed and calibrated an advanced subcellular element-based computational model called Epi-Scale that enables quantitative testing of hypothesized mechanisms governing epithelial cell behavior within the developing tissue microenvironment. Regression analysis of predictive model simulation results reveals the relative contributions of osmotic pressure, cell-cell adhesion and cortical stiffness to mitotic rounding and establishes a novel mechanism for ensuring robustness in mitotic rounding within densely packed epithelia.

Introduction

Epithelia are tissues composed of tightly adherent cells that provide barriers between internal cells of organs and the environment and are one of the four basic tissue types in the human body [13] (Fig 1). Epithelial expansion driven by cell proliferation is a key feature throughout development, and occurs in hyperplasia, a precursor to cancer. Cell divisions during development must occur robustly, as mis-segregation of chromosomes leads to severe genetic abnormalities such as aneuploidy [4]. Over 90% of human tumors are derived from epithelia [5], and the accumulation of genetic errors during cell division can lead to all of the hallmarks of cancer [6]. Division in epithelia is further complicated by the need for a dividing cell to stay connected to its neighbors [7]. A deeper understanding of the biophysical mechanisms governing the behavior of mitotic cells in epithelia will result in a better understanding of many diseases including cancer.

thumbnail
Fig 1. Epithelial mechanics and workflow outline.

(A) Apical surface of epithelial cells within the Drosophila wing imaginal disc that are marked by E-cadherin tagged with fluorescent GFP (DE-cadherin::GFP). Multiple cells within the displayed region are undergoing mitotic rounding with a noticeable decrease in fluorescent intensities of E-Cadherin. (B) Experimental image of cross-section of wing disc marking levels of actomyosin (Myosin II::GFP). (C) Cartoon abstraction of epithelial cells, which are polarized with apical and basal sides. Actomyosin and mechanical forces during mitotic rounding are primarily localized near the apical surface. (D) At the molecular scale, the boundary between cells consists of a lipid bilayer membrane for each cell, E-cadherin molecules that bind to each other through homophilic interactions, and adaptor proteins that connect the adhesion complexes to an underlying actomyosin cortex that provides tensile forces along the rim of apical areas of cells. (E) The graphical workflow of the computational modeling setup, calibration, verification and predictions. Arrows indicate mitotic cells. Scale bars are 10 micrometers.

https://doi.org/10.1371/journal.pcbi.1005533.g001

Epithelial cells entering mitosis rapidly undergo structural changes that result in the apical surface of the cell becoming larger and rounder, in a process known as mitotic rounding (MR) [8,9]. MR occurs in detached cells, cells adherent to a substrate as well as in epithelial cells within tissues [1012]. MR in epithelia coincides with an increased polymerization of actomyosin at the cell cortex, which results in an increase in cortical stiffness [4,11]. Simultaneously, the intracellular pressure increases [11], and cells partially reduce adhesion to their neighbors and the substrate [4].

However, the roles of cell-cell adhesion, cell stiffness, and intracellular pressure during mitotic rounding are not fully resolved in cultured cells, and even less is known in the tissue context [13]. For example, Stewart et al. [11] indicates that both pressure and the actin-myosin cortex are important for mitotic swelling while Zlotek-Zlotkiewics et al. [14] observe that the actin-myosin cortex is not involved in mitotic swelling. Further, it is technically challenging to modulate the mechanical properties of individual mitotic cells in tissues with small perturbations that do not “break” the system. Thus, this gap-in-knowledge is currently extremely hard to address experimentally.

Recently, computational modeling coupled with experimentation has become a powerful tool for identifying the biophysical mechanisms governing organogenesis [1520]. MR is investigated in this paper by using a novel multi-scale sub-cellular element model (SEM) called Epi-Scale that simulates epithelial cells in growing tissues. New biologically relevant features of the model include: i) separate representations of the apical membrane and cytoplasm, as well as cell-cell interactions at the tissue scale; ii) a systematic calibration of the model parameters to provide accurate biological simulations of cell division and tissue growth; and iii) a detailed description of cell properties during mitotic rounding.

We used multi-scale model simulations and response surface methodology (multiple linear regression) [21,22] to investigate the extent to which a mitotic cell-cell adhesion, cortical stiffness, and internal pressure impact the roundness and apical cross-sectional area of mitotic cells. The quantitative analysis of model simulations demonstrated that increased cytoplasmic pressure is the main driver of the increase of mitotic cell’s apical area, which was balanced by both cortical stiffness and cell adhesivity. Increased cortical stiffness and decreased adhesion is shown to promote cell roundness. Surprisingly, within the range of experimentally observed MR values, the relative roundness of cells was not sensitive to small perturbations in cytoplasmic pressure. Understanding how perturbed mechanical properties such as cytoplasmic pressure, cell-cell adhesion and cortical stiffness impact mitotic rounding have important implications on, amongst others, how tumors progressively become more genetically unstable (chromosomal aneuploidy) and more aggressive.

The paper is organized as follows. The Methods section describes modeling background and new model description. The Results section provides details of calibration of single cell model parameters using quantitative experimental data. Calibrated model simulations are shown to predict emergent properties of epithelial topology without requiring further calibration using tissue-level properties. The model is then used to quantify the relative impacts of cell-cell adhesion, membrane stiffness and intracellular pressure on MR using two separate criteria: apical area and apical roundness. The paper ends with the Discussion section, which puts predictions of the model in more general biological context. It also describes future extensions of the computational model environment for simulating epithelial tissue mechanics in greater biological detail.

Methods

Modeling background

Multiple computational approaches have been utilized to model various aspects of epithelial tissue dynamics, each with its particular focus and applications (see, amongst others, reviews [15,2327]). For example, the Cellular Potts Modeling (CPM) approach has been used successfully to take into account cell adhesivity for studying cell aggregation as well as cell morphogenesis (see, amongst others, [2831]). Finite element models (FEMs) and models based on solving Naiver-Stokes equations have also been implemented to investigate cell growth and division [3235]. Vertex based models (VBM) provided an efficient and fast approach to study regulation of cell topology, tissue-size regulation, tissue morphogenesis, and the role of cell contractility in determining tissue curvature [23,3641]. In VBMs, cellular shapes are defined by the shared vertices of neighboring cells and edges between them.

The Subcellular Elements Model (SEM), developed initially by Newman’s group [42] for simulating multi-cellular systems to encompass multiple length scales, has been now adopted by many groups as a general computational modeling approach. A particular advantage of the SEM approach is that it can provide local representations of mechanical properties of individual cells which can be directly related to the experimental data [43]. Each cell in a SEM consists of a set of nodes representing a coarse-grained representation of subcellular components of biological cells. Node-node interactions are represented by energy potentials. SEMs have been extended to predict how mechanical forces generated by cells are redistributed in a tissue and for studying tissue rheology, blood clot deformation, and cell-cell signaling [4446]. For example, a SEM model with GPU implementation was used to compare multiple mechanisms governing the formation of stratified layers of the epidermis [19] as well as mechanisms governing intestinal crypt homeostasis [47]. Jamali et al. [48] also developed an SEM model to represent the membrane and nucleus of the cell by nodes connected by overdamped springs. Gardiner et al. [49] described a SEM with locally-defined mechanical properties. Christely et al. [45] have developed an efficient computational implementation of the SEM simulating role of Notch signaling in cell growth and division, on GPU clusters to decrease computational time. A SEM model was also used to study aspects of epithelial cell mechanics without making assumptions about cell shapes [50].

Multi-component computational model of epithelia

We describe in this section novel multi-scale SEM computational platform called Epi-Scale which simulates the growth of flat epithelial monolayers. Model simulations focus on representing two-dimensional (2D) planar cell shapes near the apical surfaces of cells of the Drosophila wing imaginal disc, which is popular model to study the biophysics and genetics of epithelial tissue growth. The 2D planar model is a common simplifying approximation that was used in many previous models of wing disc growth [18,38,39,51,52]. It is reasonable to use a 2D model for studying many epithelial processes in the Drosophila wing disc pouch because it consists of a single layer of cells and the essential structural components of those cells, including E-cadherins and actomyosin, are concentrated on the apical surface of the epithelia (Fig 1A–1D). E-cadherin is responsible for adhesion between two neighboring cells, and actomyosin, which is concentrated near the apical surface, drives cell contractility. The nucleus and most of the cytoplasm are pushed up to the apical surface during cell division. Using a 2D approximation also allows us to model a large number of cells with high resolution and with special attention to mechanical cell properties. The future development of the Epi-Scale simulation platform implemented on GPU clusters, will also enable 3D simulations with reasonable computational costs.

In what follows, we first describe different types of the sub-cellular nodes that are used to simulate each cell, and the interactions between them. Then, the equations of motion of each subcellular element are provided. Finally, approaches for modeling cell’s growth, transition to mitotic phase, and division are described. The workflow of the model is shown in Fig 1E.

Sub-cellular elements

Epi-Scale represents individual cells as collections of two types of interacting subcellular elements: internal nodes and membrane nodes (Fig 2). The internal nodes account for the cytoplasm of the cell, and the membrane nodes represent both plasma membrane and associated contractile actomyosin cortex. The internal and membrane nodes are placed on a 2D plane, representing the apical surface of epithelia.

thumbnail
Fig 2. Diagram of the underlying physical basis of model simulations.

(A) Intracellular and intercellular interactions between different elements of the model. Symbols and notations are indicated in the legend. (B) Implementation of the simulation of cell cycle in the model.

https://doi.org/10.1371/journal.pcbi.1005533.g002

Interactions between internal and membrane nodes are modeled using potential energy functions as shown in Fig 2A [45,53]. Combined interactions between pairs of internal nodes (EII) represent the cytoplasmic pressure of a cell. Combined interactions between internal nodes and membrane nodes of the same cell (EMI) represent the pressure from cytoplasm to the membrane. Interactions between membrane nodes of the same cell (EMMS) are used to model the cortical stiffness. Cell-cell adhesion (EAdh) is modeled by combining pairwise interactions between nodes of the membranes of two neighboring cells. EMMD is a repulsive Morse potential function between membrane nodes of neighboring cells that prevents membranes of adjacent cells from overlapping. Epi-Scale utilizes spring and Morse energy potential functions to simulate the interactions between subcellular elements. Linear and torsional springs are represented by energy functions EMMS and EAdh [34,54], while Morse potential functions are used in energy function EMI, EII, and EMMD [46] (see Table 1 and Fig 2).

The Morse potential consists of two terms, generating short-range repulsive and long-range attractive forces [42]. For example, the following expression is a Morse potential function used in EMI to represent an interaction between internal node i and membrane node j: (1) where UMI, WMI, ξMI, and γMI are Morse parameters. The same form of the potential with different sets of parameters is used for EII and EMMD (Table 2). These potential functions govern the motion of internal and membrane nodes inside the cells resulting in the deformation and rearrangement of cells within the tissue. S2 Appendix provides the formulation of potential energy functions used in the Epi-Scale model.

Equations of motion of individual nodes

Displacement of each internal or membrane node is calculated at each moment in time based on the potential energy functions. The model assumes that nodes are in an overdamped regime [20,38,53] so that inertia forces acting on the nodes can be neglected. This leads to the following equations of motion describing movements of internal and membrane nodes, respectively: (2) (3) where η is the damping coefficient, and are positions of internal node and membrane nodes indicated by indices i and j. m is the index for any internal node interacting with the internal node i. k is the index for any membrane node of the same cell interacting with the membrane node j. Finally, l is the index for any membrane node of different cell interacting with the membrane node j. Note that adhesion between membranes of two neighboring cells is represented as pair-wise interaction between membrane nodes. Consequently, no summation with respect to different nodes is needed in Eq 3. Eqs 2 & 3 are solved at the same time for all NI internal nodes and NM membrane nodes.

Eqs 2 and 3 are discretized in time using forward Euler method and positions of nodes and are incremented at discrete times as follows (4) where Δt is the time step size. The same discretization technique is used for the equations of motion of the membrane nodes.

Epi-Scale platform is computationally implemented on a cluster of Graphical Processing Units (GPUs). This enables us to run simulations with subcellular resolution at the micro-scale with a reasonable computational cost and to study the impact of changes in individual cell mechanical properties on the tissue development at the macro-scale. S1 Appendix provides details about the simulation algorithm, GPU implementation and computational cost.

Cell cycle

Model parameters were set based on experimental values determined from studies of Drosophila wing disc development, an established genetically accessible model of organ development [55]. The growth of the wing disc is spatially uniform and decreases over time [56]. The growth rate for cell i is modeled by an exponentially decaying function fit to the experimental data for Drosophila wing disc [56], with a random term representing stochastic variation among cells: (5) where is the average growth rate of cells in the beginning of the simulation and Rnd[−g0, g0] is a random number chosen using a uniform distribution in the range of [−g0, g0]. kg is the decay constant of the growth rate.

Cells cycle through interphase and mitosis phases in the simulation. The variable Cell Progress (CP ϵ [0,1]) describes progress of a cell through the cell cycle from the beginning of the interphase (CP = 0) to the end of the cell division (CP = 1). CP is updated based on cell growth rate as follows: (6)

The number of internal nodes of a cell increases as the cell grows. The number of initial and final nodes can be varied based on the desired resolution of a single cell. Simulations in this work start with 20 internal nodes at CP = 0 and end with 40 internal nodes at CP = 1. So, an internal node is added for every 1/20 increase in CP (Fig 2) to reach the desired 40 internal nodes at the end of cell cycle. The new internal node is randomly placed within a radius 0.2Rc from the center of the cell, where Rc is the radius of the cell. Epithelial cells undergoing mitosis increase their intracellular pressure by adjusting their osmolarity relative to their surroundings [57]. Additionally, the actomyosin cortex is enriched, and cellular adhesion to the substrate and to neighboring cells are downregulated [11,5862]. Since these changes in mitotic cells occur concurrently, the relative impact on a mitotic cell cannot be easily decomposed in experiments into separable effects.

To simulate MR, parameters regulating cell-cell adhesion, cortical stiffness, and internal pressure of cells in the mitotic phase (M phase) are varied linearly from interphase parameter values to mitotic parameter values to represent the changes in cell mechanical properties during mitosis (see Table 2) [10,11,61]. For example, UMI, Morse parameter that determines cytoplasmic pressure on the membrane of the cell (see Section 4 in S3 Appendix), was varied from the interphase value () to the mitotic value (), by using the following function: (7)

Similar linear expressions are used for representing enrichment of the actomyosin cortex and reduction in cell adhesion with neighboring cells in mitotic phase (Table 1).

Cells in the mitotic (M) phase—which lasts approximately 30 minutes—divide into two daughter cells (Fig 2). Cytokinesis occurs when CP approaches 1 and is modeled by separating internal and membrane nodes of the mother cell into two sets representing daughter cells. The axis of division is implemented perpendicular to the cell’s longest axis, following Hertwig’s rule [63], prior to the initiation of mitotic rounding [64]. New membrane nodes are created along the cleavage plane for each daughter cell. After division, parameters for nodes of each daughter cell are set back to calibrated interphase values and CP is set to zero for both daughter cells.

Membrane nodes in the beginning of a simulation are arranged in a circle for each cell, and internal nodes are randomly placed within each cell (Fig 3A). After initialization, internal nodes rapidly rearrange in every cell and cells self-organize into a polygonal network, similar to the experimentally observed cell packing geometry of epithelia (Fig 3B). Cells in a simulation constantly grow, divide and interact with each other resulting in a detailed dynamic representation of the developing epithelial tissue (see Fig 3C and 3D and S1 Video).

thumbnail
Fig 3. Initial conditions and sample simulation output.

(A) Initial condition of a simulation with seven initially non-adherent circular cells. Each cell starts with 100 membrane elements and 20 internal elements. (B) Initial formation of an epithelial sheet after cells adhere to each other. An equilibrium distribution of internal nodes is reached for each cell. (C) Epithelial sheet after 55 hours of proliferation. (D) Enlarged view of the selected region showing different cell shapes and sizes due to interactions between cells. The large cell is undergoing mitotic rounding (MR).

https://doi.org/10.1371/journal.pcbi.1005533.g003

Results

Model calibration

Model parameters were calibrated using experimental data for the third instar Drosophila wing disc, which is a powerful model for studying organ formation [24,65] (Figs 4 and 5). Experimental values for similar cell lines were used to calibrate the model parameters when experimental data for Drosophila wing disc were not available.

thumbnail
Fig 4. Calibration of model parameters through simulations.

(A-A″) Calibration test to determine parameters for cell elasticity, analogous to experimental single cell stretching tests [66], (A) Initial condition t = 0, (A′) 6 minutes after simulation with no force applied, (A″) after 72 minutes cell is completely on tension (B-B″) Cell adhesivity test, analogous to experimental tests [69] for calibrating the level of cell-cell adhesion between adjacent cells. (B) Initial condition t = 0, (B′) 6 minutes after simulation begins with no force applied, (B″) after 72 minutes, 15 nN force is applied. (C) Stress versus strain for single cell calibration (red line) and stress versus strain for calibrating the level of adhesivity between the two cells (blue line) [69,70]. Initial negative strain in adhesivity test is due to strong adhesion between two cells. (D) Force and strain as a function of time for adhesivity test. (E) Tissue growth rate calibration by comparing with the experimental data by Wartlick et al. [56]. The 95% confidence interval for the growth rate results is shown in grey color.

https://doi.org/10.1371/journal.pcbi.1005533.g004

thumbnail
Fig 5. Dynamics of mitotic rounding.

(A-A‴) Time-lapse confocal images of cell undergoing mitosis in the wing disc with E-cadherin:GFP-labeled cell boundaries. Scale bar is 5 μm. Arrows indicate daughter cells. (B-D‴) Time series from Epi-Scale simulation of a cell undergoing mitosis and division with illustration of: (B-B‴) adhesive spring stiffness, (C-C‴) cortical spring stiffness, and (D-D‴) internal pressure, respected to their interphase values. (E-F) Comparison of size and roundness of mitotic cells with experimental data for the Drosophila wing disc. Arrow represents mitotic cell in B-D. A t-test comparing the means of computational simulations and experiments result in p = 0.72 for cell area ratio and p = 0.76 for normalized roundness of mitotic cells.

https://doi.org/10.1371/journal.pcbi.1005533.g005

The mechanical stiffness of the actomyosin cortex () was calibrated using the modulus of elasticity (E) of a single cell [66]. E was experimentally obtained by applying forces to opposite sides of a cell and measuring cell deformation [67,68]. This experiment was reproduced in the Epi-Scale model simulation by applying a linearly increasing force to membrane nodes on both sides of a simulated cell and calculating cell’s deformation (Fig 4A–4A”). The slope of the graph of the stress versus strain (Fig 4C) provides elasticity of the cell. The elasticity of a single cell is calibrated by adjusting linear stiffness of springs representing interactions between membrane nodes of a single cell. We have chosen value of the () so that that E = 19 kPa, which is within the biological range of 10 − 55 kPa measured for epithelial cells [67,68].

The cell-cell adhesive force (Fadh) is experimentally determined by measuring the force needed to detach two adhered cells from each other. This experiment is reproduced in silico by applying forces to membrane nodes on either side of two adhered cells, and measuring the force needed to separate them (Fig 4B–4B”). The strength of the cell-cell adhesion for the Drosophila epithelium has not been measured yet. was calibrated so that 10 nN/μm was required to detach two adhered cells from each other, based on published data for S180 cells transfected to express E-cadherin [69] and data from epithelial MDCK cells, which have adherens junctions similar to those along the apical surface of the Drosophila epithelium [70] (Fig 4D). (More details about cell-cell adhesion calibration are provided in S5 Appendix.)

Cells in the wing disc have spatially-uniform growth-rates that slow down as the tissue approaches its final size [56]. The growth rate in the Epi-Scale model described by Eq (5) was calibrated (Table 3) so that the number of cells in the model simulations matched experimental data for the wing disc pouch [56] (Fig 4E).

During mitosis, apical cell area and roundness increase compared to their interphase values (Fig 5A–5A”‘). This correlates with an observed increase in the cell’s internal pressure, cortical stiffness and decrease of intercellular adhesion marked by noticeable reduction in E-Cadherin. To simulate this, parameters of EII and EMI, which can be changed to increase the cytoplasmic pressure (ΔP) (Section 4 in S3 Appendix), cortical stiffness (kStiff), and adhesivity (kAdh) are varied from their interphase values to their mitotic values (Fig 5b–5d). Values were selected such that the ratio of mitotic cell area to interphase cell area (Amit/Ainter) and cell roundness (Rnorm) were calibrated to data collected on mitotic cells from the wing disc (Fig 5E and 5F). The methods for calculating of area, roundness, and pressure are described in S3 Appendix.

Tissue topology emerges from cell self-organization driven by cellular mechanics

After calibration of the model parameters at the cellular scale, validation simulations were run to determine whether the cellular-scale calibration was sufficient to recapitulate expected topological properties of the tissue (Fig 6A) [71,72]. One metric for tissue topology is the distribution of cell neighbor numbers, or polygon class distribution. The polygon class distribution in Epi-Scale simulations approaches to steady state after 35 hours (Fig 6B). This steady state distribution matches the distributions observed in experiments with the wing disc and other epithelial systems [43] (Fig 6C) as well as obtained using other computational models such as vertex based model [38]. We also confirmed that simulations recapitulate experimental observations [63] that cells entering mitosis on average gain a cell bond, increasing the number of neighbors by one (Fig 6D, inset). Further, we investigated the effects of varying cellular modulus of elasticity on the polygon class distribution in the range reported values of epithelial cells (10 − 55 kPa) [67]. The results show that the polygon class distribution is insensitive to the changes in the elasticity values (Fig 6E). This is a reasonable result since the polygon class distribution is strongly conserved among a wide range of epithelial tissues (Fig 6C). Therefore, cellular modulus of elasticity observed for the range of epithelial cells does not impact the polygon class distribution.

thumbnail
Fig 6. Emergence of tissue-level statistics from model simulations.

(A) Sample simulation output showing cells with different numbers of neighbors as different colors (B) Simulations initiated from seven cells reaches steady-state polygon-class distribution after approximately 35 hours of cell proliferation. (C) Comparison of polygon class distributions obtained by Epi-Scale model with various biological systems (data extracted from [79]) and a vertex based model by Farhadifar et al. [38]. (D) Polygon class distribution of cells at different stages of growth, and comparison of mitotic cells distribution with Drosophila wing disc experimental data [63]. (E) Polygon class distribution of cells at different level of cell’s elasticity. The results do not show sensitivity in the range of reported elasticity of epithelial cells [67]. (F) Average relative area (), and average polygon class of neighboring cells verifying that simulation results satisfy Lewis law and Aboav-Weaire law. A is the apical area of cell and is the average apical area of the population of cells.

https://doi.org/10.1371/journal.pcbi.1005533.g006

We further verified that the polygon class distribution of the simulated tissue satisfies three laws describing topological relationships: Euler’s law, Lewis law, and Aboav-Weaire Law. Euler’s law states that cells forming a packed sheet should be hexagonal on average [72,73]. The Lewis law states that cells with more neighbors should have larger normalized area [73]. The Aboav-Weaire law indicates that the average polygon class of neighbors of each cell decreases as the cell’s polygon class increases [74]. Simulation results obtained using calibrated model, show the average side of cells to be equal to 5.98 for interphase and mitotic cells, 5.80 for interphase cells, and 6.49 for mitotic cells. Model simulations also satisfy two other laws as shown in Fig 6D when interphase cells are counted.

Impacts of adhesion, stiffness, and cytoplasmic pressure on mitotic rounding

The Epi-Scale model is suitable for generating and testing hypotheses regarding mechanical mechanisms of MR because it can represent non-polygonal shapes of cells. Simulations were conducted to predict the relative contributions of different cell properties to the relative area ratio (Amit/Ainter) and normalized roundness (Rnorm) of mitotic cells as calculated in S3 Appendix. Parameter values were selected in a three-level full factorial design (FFD) to investigate the relationships between the mitotic parameters of the model (, , and ΔP) and mitotic rounding (Aratio and Rnorm) (Fig 7A, S9 Appendix). A regression model was fit to the results of the FFD, and showed that for large changes, ΔP was the primary regulator of Aratio, and and were the primary regulators of Rnorm (Fig 7C and 7D).

thumbnail
Fig 7. Response surface method analysis of mechanical properties on regulating mitotic expansion and mitotic rounding.

(A) Schematic of initial full factorial design (FFD) for exploring parameter space, and subsequent central composite design (CCD) for developing the response surface models shown in (C, D). (B) Pareto front indicating computational model parameter values with lowest difference with experimental data for area ratio and normalized roundness. The parameter range defined by the CCD (Run 2) spans parameter variation where the error between experiments and simulations is within the propagated uncertainty of measurements and simulations. Error bars are the standard error of means of the normalized deviation between experiments and simulations. (C-D) Contour plots for FFD experiment where (C) shows the area ratio (Aratio = Amit/Ainter) and (D) shows the normalized roundness (Rnorm). (E-F) Contour plots for CCD experiment where (E) shows the area ratio (Aratio) and (F) shows the normalized roundness (Rnorm).

https://doi.org/10.1371/journal.pcbi.1005533.g007

A region of parameter space was selected where the error in mitotic rounding measurements (Aratio and Rnorm) was minimized as shown in the Pareto front (Fig 7B, 7E and 7F, S9 Appendix). The region of parameter space closest to experimental values of cell area and roundness was explored with a central composite design (CCD) as a second iteration to more precisely determine the relative contribution of each physical parameter on MR within experimentally observed ranges (Fig 7A) [75]. This result quantitatively defines the predicted variation in mitotic cell-cell adhesion, stiffness and pressure that explains the variation in mitotic area ratio and rounding observed in mitotic epithelial cells. To keep mitotic rounding within the range of variation observed, ΔP must be tightly regulated (~19% variation about the calibrated point), whereas the requirements for and are less stringent ~120% and ~67% respectively.

Model reduction (S10 Appendix) revealed that regulation of mitotic rounding is approximated well by linear regression models for parameter evaluation resulting in physiological values of Aratio and Rnorm for Drosophila wing disc cells (Fig 5). This suggests that regulation is in the linear regime, which is a good attribute for tightly controlled processes. Since interaction terms are not significant, cell mechanical properties contribute to cell shape changes independently. Mitotic pressure was found to be the primary regulator of mitotic cell area (Fig 7E), while both cell-cell adhesion and cortical stiffness reduced area expansion slightly. An increase in cell-cell adhesion was shown to reduce roundness whereas increased cortical stiffness promoted roundness for small perturbations (Fig 7F).

To define the relative impacts of mechanical properties on Aratio and Rnorm under physiological or “wild-type” conditions, local sensitivity analysis (Fig 8A and 8B, S11 Appendix) was performed after application of the stepwise model reduction (S10 Appendix). Within the physiologically relevant domain of the parameter space, pressure strongly regulates mitotic area expansion but does not have a strong impact on the shape roundness (Figs 7E and 8A). Stiffness and adhesion are important in tuning the degree of mitotic roundness (Figs 7F and 8B).

thumbnail
Fig 8. Quantitation of relative sensitivity of mitotic area expansion and roundness to adhesion, stiffness and pressure changes within the physiological property space.

Sensitivity estimation of (A) (Amit/Ainter) and (B) Rnorm to small perturbation in the three mitotic parameter set points, , , and ΔP. Sensitivity was estimated from the reduced RSM model described in Fig 7C–7F after stepwise model regression (p-value cutoff of 0.01). (C) Proposed mechanical regulatory network defined for “physiological ranges” within the parameter ranges defined by the CCD (Run 2, Fig 7A) that summarizes the local sensitivity analysis. Cell adhesivity, an increase in , slightly inhibits area expansion and strongly inhibits roundness. Membrane stiffness, inhibits area expansion and promotes roundness. Mitotic area expansion is most sensitive to variation in the mitotic pressure change (ΔP), but pressure has little effect on roundness over the calibrated physiological ranges.

https://doi.org/10.1371/journal.pcbi.1005533.g008

These results are summarized in the form of the mechanical sensitivity model in Fig 8C, analogous to protein interaction networks. This model describes how small variations in each cellular mechanical property impact relative mitotic area expansion and roundness.

Discussion

The roles of pressure, stiffness and adhesion in mitotic cells in single cell culture, in suspension or attached to substrates, are still not resolved in the experimental literature and largely unexplored in the tissue context [11,14]. We described in this paper a novel multi-scale sub-cellular model, called Epi-Scale, for simulating mechanical and adhesive properties of cells in the developing columnar epithelium of the wing disc, which consists of a single layer of cells. The model approximates the tissue as a 2D surface since the majority of the contractile and adhesive forces are localized at the apical surface of the epithelium (Fig 1B).

Parameter ranges for the computational model were obtained by calibrating the model using single cell stretching experiments, experiments on stretching a pair of cells adhered to each other, dynamic experimental measurements of the area and roundness of mitotic cells, and the tissue growth rate of the Drosophila wing disc. Cell-cell adhesion and cell elasticity were calibrated using data from experiments with single cells. The calibrated model was verified by successfully reproducing emergent properties of developing tissue such as the polygon class distributions for both interphase and mitotic cells without additional calibration or parameter tuning.

Epi-Scale enables the systematic generation and testing of new hypotheses about the underlying mechanisms governing mitotic rounding within the developing tissue microenvironment. Regression analysis of predictive simulations provided complete assessment of the quantitative contributions of cytoplasmic pressure, cell-cell adhesion and cortical stiffness to mitotic cell rounding and expansion (Figs 7 and 8). Mitotic cell area expansion was shown to be largely driven by regulation of cytoplasmic pressure. Surprisingly, the variability in mitotic roundness within physiological ranges was shown to be primarily driven by varying cell-cell adhesivity and cortical stiffness, rather than pressure.

It is currently challenging to target only dividing cells in a tissue. One experimental approach that might be used in the future for testing the model predictions would be to regulate the expression of E-cadherin, Myosin-II, and osmotic channel antagonists under a Cyclin B promotor, active during mitosis, resulting in modulation only in dividing cells [76,77]. Alternatively, opto-genetic methods could be employed to selectively regulate individual cell properties [78].

The simulation results have also shown that increases of the mitotic rounding under super-physiological pressure (greater than calibrated values) could result in cell-cell rearrangements (T1 transitions) of the neighboring cells, due to rapid increase of the apical surface of the mitotic cell (S7 Appendix). This indicates that Epi-Scale platform could be used for future detailed studies of epithelial morphogenesis.

We have shown that our model simulations provide new insights into the individual contributions of cell properties to MR. Determining which aspects of mitotic rounding are most sensitive to perturbed cell properties in dense tissues, including solid tumors, can help direct future efforts to identify cellular processes that specifically block mitosis in highly proliferative tumors, but that are not damaging to non-proliferative cells [14]. As a flexible computational modeling platform, Epi-Scale can be extended to simulate a wide range of multi-cellular processes, including epithelial morphogenesis, wounding healing and blood clot formation.

Supporting information

S1 Appendix. Computational implementation and computational cost of the Epi-scale model.

https://doi.org/10.1371/journal.pcbi.1005533.s001

(PDF)

S2 Appendix. Potential energy functions used in the Epi-scale model.

https://doi.org/10.1371/journal.pcbi.1005533.s002

(PDF)

S3 Appendix. Metrics for determining polygon class, roundness of cell shape and internal pressure of a cell.

https://doi.org/10.1371/journal.pcbi.1005533.s003

(PDF)

S5 Appendix. Calibration of the adhesion between neighbouring cells.

https://doi.org/10.1371/journal.pcbi.1005533.s005

(PDF)

S6 Appendix. Transgenic fly lines, microscopy and image processing.

https://doi.org/10.1371/journal.pcbi.1005533.s006

(PDF)

S7 Appendix. Mitotic rounding effects on cell-cell rearrangements.

https://doi.org/10.1371/journal.pcbi.1005533.s007

(PDF)

S2 Video. Mitotic cells at high pressure lead to T1 transitions around the mitotic cells.

https://doi.org/10.1371/journal.pcbi.1005533.s014

(MP4)

Acknowledgments

The authors acknowledge support from the Notre Dame Integrated Imaging Facility. They thank members of the Zartman’s lab for critical feedback on manuscript drafts as well as Luis Lazalde for technical analysis of segmented cells.

Author Contributions

  1. Conceptualization: MA JJZ ZX.
  2. Data curation: AN WS PAB AA CN.
  3. Formal analysis: MA JJZ AN PAB AA CN ZX.
  4. Funding acquisition: MA JJZ ZX.
  5. Investigation: MA JJZ AN WS PAB AA CN ZX.
  6. Methodology: MA JJZ AN WS AA ZX.
  7. Project administration: MA JJZ ZX.
  8. Resources: MA JJZ ZX.
  9. Software: WS AN AA ZX PAB MA.
  10. Supervision: MA JJZ ZX AN.
  11. Validation: MA JJZ AN WS PAB AA CN ZX.
  12. Visualization: MA JJZ AN PAB AA CN ZX.
  13. Writing – original draft: AN JJZ MA WS PAB AA ZX CN.
  14. Writing – review & editing: AN JJZ MA PAB AA ZX CN.

References

  1. 1. Gumbiner B. Structure, biochemistry, and assembly of epithelial tight junctions. Am J Physiol—Cell Physiol. 1987;253: C749–C758.
  2. 2. Fristrom D. The cellular basis of epithelial morphogenesis. A review. Tissue Cell. 1988;20: 645–690. pmid:3068832
  3. 3. Lecuit T, Lenne P-F. Cell surface mechanics and the control of cell shape, tissue patterns and morphogenesis. Nat Rev Mol Cell Biol. 2007;8: 633–644. pmid:17643125
  4. 4. Ramkumar N, Baum B. Coupling changes in cell shape to chromosome segregation. Nat Rev Mol Cell Biol. 2016;advance online publication. pmid:27353479
  5. 5. Miller SJ, Lavker RM, Sun T-T. Interpreting epithelial cancer biology in the context of stem cells: tumor properties and therapeutic implications. Biochim Biophys Acta. 2005;1756: 25–52. pmid:16139432
  6. 6. Hanahan D, Weinberg RA. The Hallmarks of Cancer. Cell. 2000;100: 57–70. pmid:10647931
  7. 7. Cadart C, Zlotek-Zlotkiewicz E, Le Berre M, Piel M, Matthews HK. Exploring the Function of Cell Shape and Size during Mitosis. Dev Cell. 2014;29: 159–169. pmid:24780736
  8. 8. Strangeways TSP. Observations on the Changes Seen in Living Cells during Growth and Division. Proc R Soc Lond Ser B Contain Pap Biol Character. 1922;94: 137–141.
  9. 9. McConnell CH. Mitosis in Hydra. Mitosis in the Ectodermal Epithelio-Muscular Cells of Hydra. Biol Bull. 1933;64: 86–95.
  10. 10. Sorce B, Escobedo C, Toyoda Y, Stewart MP, Cattin CJ, Newton R, et al. Mitotic cells contract actomyosin cortex and generate pressure to round against or escape epithelial confinement. Nat Commun. 2015;6: 8872. pmid:26602832
  11. 11. Stewart MP, Helenius J, Toyoda Y, Ramanathan SP, Muller DJ, Hyman AA. Hydrostatic pressure and the actomyosin cortex drive mitotic cell rounding. Nature. 2011;469: 226–230. pmid:21196934
  12. 12. Stewart MP, Toyoda Y, Hyman AA, Müller DJ. Tracking mechanics and volume of globular cells with atomic force microscopy using a constant-height clamp. Nat Protoc. 2012;7: 143–154. pmid:22222789
  13. 13. Champion L, Linder MI, Kutay U. Cellular reorganization during mitotic entry. Trends Cell Biol. 2017;27: 26–41. pmid:27528558
  14. 14. Zlotek-Zlotkiewicz E, Monnier S, Cappello G, Berre ML, Piel M. Optical volume and mass measurements show that mammalian cells swell during mitosis. J Cell Biol. 2015;211: 765–774. pmid:26598614
  15. 15. Brodland GW. How computational models can help unlock biological systems. Semin Cell Dev Biol. 2015;47–48: 62–73. pmid:26165820
  16. 16. Marin-Riera M, Brun-Usan M, Zimm R, Välikangas T, Salazar-Ciudad I. Computational modeling of development by epithelia, mesenchyme and their interactions: a unified model. Bioinformatics. 2016;32: 219–225. pmid:26342230
  17. 17. Shinbrot T, Chun Y, Caicedo-Carvajal C, Foty R. Cellular morphogenesis in silico. Biophys J. 2009;97: 958–967. pmid:19686642
  18. 18. Mao Y, Tournier AL, Hoppe A, Kester L, Thompson BJ, Tapon N. Differential proliferation rates generate patterns of mechanical tension that orient tissue growth. EMBO J. 2013;32: 2790–2803. pmid:24022370
  19. 19. Gord A, Holmes WR, Dai X, Nie Q. Computational modelling of epidermal stratification highlights the importance of asymmetric cell division for predictable and robust layer formation. J R Soc Interface. 2014;11: 20140631. pmid:25100322
  20. 20. Kursawe J, Brodskiy PA, Zartman JJ, Baker RE, Fletcher AG. Capabilities and Limitations of Tissue Size Control through Passive Mechanical Forces. PLoS Comput Biol. 2015;11: e1004679. pmid:26713738
  21. 21. Anderson MJ. RSM Simplified: Optimizing Processes Using Response Surface Methods for Design of Experiments. New York, New York: Productivity Press; 2005.
  22. 22. Zartman J, Restrepo S, Basler K. A high-throughput template for optimizing Drosophila organ culture with response-surface methods. Development. 2013;140: 667–674. pmid:23293298
  23. 23. Fletcher AG, Osterfield M, Baker RE, Shvartsman SY. Vertex models of epithelial morphogenesis. Biophys J. 2014;106: 2291–2304. pmid:24896108
  24. 24. Buchmann A, Alber M, Zartman JJ. Sizing it up: The mechanical feedback hypothesis of organ growth regulation. Semin Cell Dev Biol. 2014; pmid:25020200
  25. 25. Osborne JM, Fletcher AG, Pitt-Francis JM, Maini PK, Gavaghan DJ. Comparing individual-based approaches to modelling the self-organization of multicellular tissues. bioRxiv. 2016; 074351.
  26. 26. Camley BA, Rappel W-J. Physical models of collective cell motility: from cell to tissue. J Phys Appl Phys. 2017;50: 113002.
  27. 27. Mirams GR, Arthurs CJ, Bernabeu MO, Bordas R, Cooper J, Corrias A, et al. Chaste: An Open Source C++ Library for Computational Physiology and Biology. PLOS Comput Biol. 2013;9: e1002970. pmid:23516352
  28. 28. Chen N, Glazier JA, Izaguirre JA, Alber MS. A parallel implementation of the Cellular Potts Model for simulation of cell-based morphogenesis. Comput Phys Commun. 2007;176: 670–681. pmid:18084624
  29. 29. Chaturvedi R, Huang C, Kazmierczak B, Schneider T, Izaguirre JA, Glimm T, et al. On multiscale approaches to three-dimensional modelling of morphogenesis. J R Soc Interface. 2005;2: 237–253. pmid:16849182
  30. 30. Glazier JA, Graner F. Simulation of the differential adhesion driven rearrangement of biological cells. Phys Rev E. 1993;47: 2128–2154.
  31. 31. Graner F, Glazier JA. Simulation of biological cell sorting using a two-dimensional extended Potts model. Phys Rev Lett. 1992;69: 2013. pmid:10046374
  32. 32. Brodland GW, Viens D, Veldhuis JH. A new cell-based FE model for the mechanics of embryonic epithelia. Comput Methods Biomech Biomed Engin. 2007;10: 121–128. pmid:18651278
  33. 33. Zhao J, Naveed H, Kachalo S, Cao Y, Tian W, Liang J. Dynamic mechanical finite element model of biological cells for studying cellular pattern formation. 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). 2013. pp. 4517–4520. pmid:24110738
  34. 34. Rejniak KA, Dillon RH. A single cell-based model of the ductal tumour microarchitecture. Comput Math Methods Med. 2007;8: 51–69.
  35. 35. Rejniak KA. An immersed boundary framework for modelling the growth of individual cells: An application to the early tumour development. J Theor Biol. 2007;247: 186–204. pmid:17416390
  36. 36. Okuda S, Inoue Y, Adachi T. Three-dimensional vertex model for simulating multicellular morphogenesis. Biophys Physicobiology. 2015;12: 13–20.
  37. 37. Osterfield M, Du X, Schüpbach T, Wieschaus E, Shvartsman SY. Three-dimensional epithelial morphogenesis in the developing Drosophila egg. Dev Cell. 2013;24: 400–410. pmid:23449472
  38. 38. Farhadifar R, Röper J-C, Aigouy B, Eaton S, Jülicher F. The Influence of Cell Mechanics, Cell-Cell Interactions, and Proliferation on Epithelial Packing. Curr Biol. 2007;17: 2095–2104. pmid:18082406
  39. 39. Aegerter-Wilmsen T, Heimlicher MB, Smith AC, de Reuille PB, Smith RS, Aegerter CM, et al. Integrating force-sensing and signaling pathways in a model for the regulation of wing imaginal disc size. Development. 2012;139: 3221–3231. pmid:22833127
  40. 40. Jessica CY, Fernandez-Gonzalez R. Quantitative modelling of epithelial morphogenesis: integrating cell mechanics and molecular dynamics. Seminars in Cell & Developmental Biology. Elsevier; 2016. http://www.sciencedirect.com/science/article/pii/S1084952116302348
  41. 41. Sussman DM. cellGPU: massively parallel simulations of dynamic vertex models. ArXiv170202939 Cond-Mat Physicsphysics. 2017; http://arxiv.org/abs/1702.02939
  42. 42. Newman TJ. Modeling Multicellular Structures Using the Subcellular Element Model. In: Anderson DARA, Chaplain PMAJ, Rejniak DKA, editors. Single-Cell-Based Models in Biology and Medicine. Birkhäuser Basel; 2007. pp. 221–239. http://link.springer.com/chapter/10.1007/978-3-7643-8123-3_10
  43. 43. Sandersius SA, Chuai M, Weijer CJ, Newman TJ. Correlating Cell Behavior with Tissue Topology in Embryonic Epithelia. PLoS ONE. 2011;6: e18081. pmid:21559520
  44. 44. Sweet CR, Chatterjee S, Xu Z, Bisordi K, Rosen ED, Alber M. Modelling platelet—blood flow interaction using the subcellular element Langevin method. J R Soc Interface. 2011;8: 1760–1771. pmid:21593027
  45. 45. Christley S, Lee B, Dai X, Nie Q. Integrative multicellular biological modeling: a case study of 3D epidermal development using GPU algorithms. BMC Syst Biol. 2010;4: 107. pmid:20696053
  46. 46. Sandersius SA, Newman TJ. Modeling cell rheology with the Subcellular Element Model. Phys Biol. 2008;5: 015002. pmid:18403827
  47. 47. Du H, Nie Q, Holmes WR. The interplay between Wnt mediated expansion and negative regulation of growth promotes robust intestinal crypt structure and homeostasis. PLoS Comput Biol. 2015;11: e1004285. pmid:26288152
  48. 48. Jamali Y, Azimi M, Mofrad MRK. A Sub-Cellular Viscoelastic Model for Cell Population Mechanics. PLOS ONE. 2010;5: e12097. pmid:20856895
  49. 49. Gardiner BS, Wong KK, Joldes GR, Rich AJ, Tan CW, Burgess AW, et al. Discrete element framework for modelling extracellular matrix, deformable cells and subcellular components. PLoS Comput Biol. 2015;11: e1004544. pmid:26452000
  50. 50. Sandersius SA, Weijer CJ, Newman TJ. Emergent cell and tissue dynamics from subcellular modeling of active biomechanical processes. Phys Biol. 2011;8: 045007. pmid:21750367
  51. 51. Shraiman BI. Mechanical feedback as a possible regulator of tissue growth. Proc Natl Acad Sci U S A. 2005;102: 3318–3323. pmid:15728365
  52. 52. Buceta J, Herranz H, Canela-Xandri O, Reigada R, Sagués F, Milán M. Robustness and stability of the gene regulatory network involved in DV boundary formation in the Drosophila wing. PLoS ONE. 2007;2: e602. pmid:17622347
  53. 53. Newman TJ. Modeling multi-cellular systems using sub-cellular elements. ArXivq-Bio0504028. 2005; http://arxiv.org/abs/q-bio/0504028
  54. 54. Jamali Y, Azimi M, Mofrad MRK. A Sub-Cellular Viscoelastic Model for Cell Population Mechanics. PLOS ONE. 2010;5: e12097. pmid:20856895
  55. 55. Neto-Silva RM, Wells BS, Johnston LA. Mechanisms of growth and homeostasis in the Drosophila wing. Annu Rev Cell Dev Biol. 2009;25: 197–220. pmid:19575645
  56. 56. Wartlick O, Mumcu P, Kicheva A, Bittig T, Seum C, Jülicher F, et al. Dynamics of Dpp Signaling and Proliferation Control. Science. 2011;331: 1154–1159. pmid:21385708
  57. 57. Fischer-Friedrich E, Hyman AA, Jülicher F, Müller DJ, Helenius J. Quantification of surface tension and internal pressure generated by single mitotic cells. Sci Rep. 2014;4: 6213. pmid:25169063
  58. 58. Clark AG, Paluch E. Mechanics and Regulation of Cell Shape During the Cell Cycle. In: Kubiak JZ, editor. Cell Cycle in Development. Springer Berlin Heidelberg; 2011. pp. 31–73. http://link.springer.com/chapter/10.1007/978-3-642-19065-0_3
  59. 59. Kunda P, Pelling AE, Liu T, Baum B. Moesin Controls Cortical Rigidity, Cell Rounding, and Spindle Morphogenesis during Mitosis. Curr Biol. 2008;18: 91–101. pmid:18207738
  60. 60. Lancaster OM, Le Berre M, Dimitracopoulos A, Bonazzi D, Zlotek-Zlotkiewicz E, Picone R, et al. Mitotic Rounding Alters Cell Geometry to Ensure Efficient Bipolar Spindle Formation. Dev Cell. 2013;25: 270–283. pmid:23623611
  61. 61. Guillot C, Lecuit T. Adhesion Disengagement Uncouples Intrinsic and Extrinsic Forces to Drive Cytokinesis in Epithelial Tissues. Dev Cell. 2013;24: 227–241. pmid:23410938
  62. 62. Matzke R, Jacobson K, Radmacher M. Direct, high-resolution measurement of furrow stiffening during division of adherent cells. Nat Cell Biol. 2001;3: 607–610. pmid:11389447
  63. 63. Gibson WT, Veldhuis JH, Rubinstein B, Cartwright HN, Perrimon N, Brodland GW, et al. Control of the Mitotic Cleavage Plane by Local Epithelial Topology. Cell. 2011;144: 427–438. pmid:21295702
  64. 64. Bosveld F, Markova O, Guirao B, Martin C, Wang Z, Pierre A, et al. Epithelial tricellular junctions act as interphase cell shape sensors to orient mitosis. Nature. 2016;530: 495–498. pmid:26886796
  65. 65. Hariharan IK. Organ Size Control: Lessons from Drosophila. Dev Cell. 2015;34: 255–265. pmid:26267393
  66. 66. Micoulet A, Spatz JP, Ott A. Mechanical Response Analysis and Power Generation by Single-Cell Stretching. ChemPhysChem. 2005;6: 663–670. pmid:15881582
  67. 67. Kuznetsova TG, Starodubtseva MN, Yegorenkov NI, Chizhik SA, Zhdanov RI. Atomic force microscopy probing of cell elasticity. Micron. 2007;38: 824–833. pmid:17709250
  68. 68. Laurent VM, Kasas S, Yersin A, Schäffer TE, Catsicas S, Dietler G, et al. Gradient of Rigidity in the Lamellipodia of Migrating Cells Revealed by Atomic Force Microscopy. Biophys J. 2005;89: 667–675. pmid:15849253
  69. 69. Chu Y-S, Thomas WA, Eder O, Pincet F, Perez E, Thiery JP, et al. Force measurements in E-cadherin—mediated cell doublets reveal rapid adhesion strengthened by actin cytoskeleton remodeling through Rac and Cdc42. J Cell Biol. 2004;167: 1183–1194. pmid:15596540
  70. 70. Sim JY, Moeller J, Hart KC, Ramallo D, Vogel V, Dunn AR, et al. Spatial distribution of cell—cell and cell—ECM adhesions regulates force balance while main­taining E-cadherin molecular tension in cell pairs. Mol Biol Cell. 2015;26: 2456–2465. pmid:25971797
  71. 71. Chiu SN. Aboav-Weaire’s and Lewis’ laws—A review. Mater Charact. 1995;34: 149–165.
  72. 72. Sanchez-Gutierrez D, Tozluoglu M, Barry JD, Pascual A, Mao Y, Escudero LM. Fundamental physical cellular constraints drive self-organization of tissues. EMBO J. 2015; pmid:26598531
  73. 73. Lewis FT. A comparison between the mosaic of polygons in a film of artificial emulsion and the pattern of simple epithelium in surface view (cucumber epidermis and human amnion). Anat Rec. 1931;50: 235–265.
  74. 74. Aboav DA. The arrangement of grains in a polycrystal. Metallography. 1970;3: 383–390.
  75. 75. Pargett M, Rundell AE, Buzzard GT, Umulis DM. Model-Based Analysis for Qualitative Data: An Application in Drosophila Germline Stem Cell Regulation. PLOS Comput Biol. 2014;10: e1003498. pmid:24626201
  76. 76. Gavet O, Pines J. Activation of cyclin B1-Cdk1 synchronizes events in the nucleus and the cytoplasm at mitosis. J Cell Biol. 2010;189: 247–259. pmid:20404109
  77. 77. Charvin G, Cross FR, Siggia ED. Forced periodic expression of G1 cyclins phase-locks the budding yeast cell cycle. Proc Natl Acad Sci. 2009;106: 6632–6637. pmid:19346485
  78. 78. Guglielmi G, Barry JD, Huber W, De Renzis S. An Optogenetic Method to Modulate Cell Contractility during Tissue Morphogenesis. Dev Cell. 2015;35: 646–660. pmid:26777292
  79. 79. Gibson MC, Patel AB, Nagpal R, Perrimon N. The emergence of geometric order in proliferating metazoan epithelia. Nature. 2006;442: 1038–1041. pmid:16900102