article: Cells in the Drosophila retina have well-defined morphologies that are attained during tissue morphogenesis. We present a computer simulation of the epithelial tissue in which the global interfacial energy between cells is minimized. Experimental data for both normal cells and mutant cells either lacking or misexpressing the adhesion protein N-cadherin can be explained by a simple model incorporating salient features of morphogenesis that include the timing of N-cadherin expression in cells and its temporal relationship to the remodeling of cell-cell contacts. The simulations reproduce the geometries of wild-type and mutant cells, distinguish features of cadherin dynamics, and emphasize the importance of adhesion protein biogenesis and its timing with respect to cell remodeling. The simulations also indicate that N-cadherin protein is recycled from inactive interfaces to active interfaces, thereby modulating adhesion strengths between cells.
  Tissues in multicellular organisms consist of a variety of cells with specialized functions. The differences between cell types are often manifest in the individual cells' shapes, their relative positions, and their neighbor relations [1]. It has become recognized that mechanical forces within cells control cell shapes and their larger-scale organization within tissues [2], [3], [4]. These forces are generated by specific molecules present within and upon the surfaces of cells, including the actin-myosin cytoskeleton and cell-cell adhesion molecules. A fundamental issue is whether tissue morphology can be described as the sum of cell morphologies, which are individually determined by autonomous force generators. An alternative description of tissue morphology assumes that morphology is passively determined by equilibrium mechanics once a small number of force parameters is established by force-generating molecules within cells. This alternative is rationalized by experimental observations finding that cell packing in the Drosophila retina is a consequence of mechanical equilibrium [5]. Moreover, retinal tissue morphology can be quantitatively modeled by assuming global minimization of interfacial energies that are established by cellular force-generating molecules [6]. Models of other epithelial tissues using similar methods have also successfully reproduced morphological properties [7], [8], [9], giving credence to the approach. The Drosophila retina is a pseudostratified epithelium containing over 800 repeating units called ommatidia (Fig. 1). Each ommatidium is on average D≈9·m across its widest axis, and consists of twenty cells, including eight photoreceptor neurons and twelve accessory cells [10]. Four of these accessory cells (called cone cells) adhere together to form a transparent plate that acts as both the floor of the simple lens and the roof of the underlying pool of photoreceptors. Two primary pigment cells optically insulate each cone-cell group; together with the cone cells, they form the “core” structure of the ommatidium. Secondary and tertiary pigment cells form the ommatidium “frame”. Functionally, the cone cells form an “aperture” for focused light to be transmitted from the lens to the photoreceptors, while the opaque primary pigment cells delineate the aperture stop. Cells in the retina express two kinds of cell-cell adhesion molecules: E-cadherin and N-cadherin [5], [11]. All cells contain E-cadherin; however, only the four cone cells contain N-cadherin. The cadherin protein molecules become localized in a thin band of lateral cell membrane (b≈50 nm wide) corresponding to the adherens junction (AJ), which is the major site of adhesion between cells (Fig. 1B). Cadherins in one membrane bind to cadherins located in the membrane of a neighboring cell across the intercellular gap, affecting the adherence of one cell with its neighbor [5]. The binding interactions appear to be homotypic; E- and N-cadherin molecules do not bind to each other [5], even though heterotypic binding has been found for Xenopus cadherins [12]. Intermolecular binding collectively generates adhesion between facing membranes, thereby decreasing the interfacial energy per AJ area and making expansion of the interfacial AJ domain energetically favorable. Membrane that has established such bonds is considered “active”. However, the expansion of one interfacial domain affects other cell-cell interfaces due to constraints on the overall sizes of the AJ membrane domain and cell volume, . Shape changes in one cell induce shape changes in others, and alteration of the elastic energy of the membranes around all cells. Ultimately, the mechanical energy of the entire ommatidium needs to be minimized globally in order to find an equilibrium configuration for the ommatidium [13], [14]. This configuration describes the entire retinal tissue since ommatidia are arranged in the epithelium as identical tiles with six-fold axes of symmetry. It is well known that cell-cell adhesion can play an important role in cell sorting and morphogenesis of tissues [15], [16], as theoretically described and experimentally demonstrated [17], [18]. More elaborate formalisms have been successfully used to simulate a variety of tissues [19], [20]. Our work indicates that in tissues like the Drosophila retina, the role of cadherins goes beyond cell sorting and in fact determines the details of their geometric shapes [6]. In order to understand how cadherin molecules control cell geometry, , one must consider the distribution and dynamic properties of cadherins. In the present work, we apply these considerations to our mechanical model. We then test the fit of various distribution and dynamic models with experimental data for mutants in which certain cells produce altered levels of N-cadherin. We demonstrate (i) that the model describes characteristic shape changes in such mutants, (ii) that the simulations distinguish between different mechanisms of how cadherin levels are attained and controlled, and (iii) that the model, although conceived as an equilibrium tool, incorporates important dynamical features in morphogenesis, such as the temporal sequence of cadherin expression and cell-cell contact remodeling.
  A cell synthesizes cadherin protein in the cytoplasm, and it is transported to the AJ domain of a cell' s outer membrane [21], [22]. If a molecule locally pairs with another molecule on a neighboring membrane, then it is stabilized both spatially and temporally [23]. Cell biologists have long known that unpaired cadherin molecules are internalized by endocytosis whereas paired molecules are not [22], [24], [25], [26], [27]. Endocytosis is critical for maintenance of epithelial adhesive integrity [28], [29], [30]. Once internalized, cadherins are either recycled back to the cell surface or trafficked to lysosomes for destruction [31], [32]. Trafficking to lysosomes entails passage through Rab5- and Rab7-enriched compartments [33]. In contrast, recycling back to the AJ occurs by two routes: directly from sorting endosomes or after transport to the recycling endosome [34], [35]. Faced with such disparate alternative pathways, what determines which route is taken after cadherins are internalized? If unpaired cadherins follow a route to destruction, then they do not have an opportunity to redistribute to other membrane domains. If unpaired cadherins are recycled, they are free to distribute on other domains of the membrane, and will accumulate along AJ interfaces with neighboring cells that contain the same cadherin. We formulated two models that simulate these different scenarios for N-cadherin. Both models assume a steady-state concentration of N-cadherin proteins in a cell. The Destruction Model assumes high rates of N-cadherin synthesis and turnover; if unpaired, molecules are rapidly endocytosed and degraded (Fig. 3A). Reaction equilibria are established locally, so that the unpaired concentration determines the concentration of N-cadherin pairs if the neighboring membrane contains the same kind of cadherin. Such local reaction equilibria at interfaces have been described in general terms [36], [37], [38]. The resulting dimensionless binding energy per membrane length is independent of the concentrations on other domains of the AJ in the same cell. The Recycling Model, by contrast, assumes a low rate of synthesis and turnover of N-cadherin. The cadherin molecules are long-lived and if unpaired, they recycle via endocytosis to and from the AJ. They then redistribute along the AJ in different ways depending on which interfaces are active, i. e. which neighbors express N-cadherin as well (Fig. 3B). In the dimensionless energy functional (1), an N-cadherin binding term along an AJ of length L is written as −γNL or alternatively −Nact Bdim/ES, where Bdim is the dimensional binding energy of a single N-cadherin pair bond, and Nact is the number of such active bonds along the interface. Thus, if a cell has Nact bonds along a total AJ length Lact, and these bonds are distributed evenly, the effective dimensionless binding strength becomes γ = β Nact/Lact, where we define the constant β = Bdim/ES. Note that Nact is not necessarily the total number of N-cadherin molecules synthesized in a cell, since some molecules may not find a binding partner. Cells contain a defined number N0 of cadherin molecules. We assume this number to be identical for C1 and C2 cells, both for simplicity and because genetic regulation of cadherin expression is not known to differ between cone cells. It is clear, however, that C2 cells have a longer active membrane length than C1 cells (see Figs. 2 and 3). In fact, we can write (2) Only a C2 cell has all of its N-cadherin molecules bound and active (NactC2 = N0), whereas some unbound molecules remain within the C1 cells (Fig. 3B). Quantitatively, the number of those molecules is N0 [1−2LC1C2/ (2LC1C2+Lcen) ]. The binding strength established in the Recycling Model between any two cone cells is thus, γN = βN0/ (2LC1C2+Lcen). Instead of explicitly defining γN as a parameter, the Recycling Model implicitly defines it via the parameter βN0 and the self-consistently determined interfacial AJ lengths between cone cells. The C2 cells can be said to be limiting for the process, because all of their N-cadherin is paired, whereas the C1 cells retain some inactive cadherin. In practice, the simulation protocol iteratively adjusts the cadherin strengths and the interface lengths until convergence is reached. In all cases presented here, accurate convergence takes very few iteration steps. The simulations of normal ommatidia are indistinguishable for the two models, as both of them establish a uniform γN along all edges between cone cells. Hence, in order to determine if either model is correct, we turned to situations in which certain cells synthesize less N-cadherin than their neighbors. If asymmetries arise in an ommatidium so that cone cells have different active lengths of AJ with respect to N-cadherin, the two models predict significantly different equilibrium morphologies. Experiments can be performed in the Drosophila eye to create ommatidia where some cells have a normal gene while other cells are missing the gene [5]. Although the cellular composition of such mosaics is generated in a random manner, it is possible to screen through hundreds of ommatidia and find examples where specific cone cells are missing a gene. This was done for the N-cadherin gene, and a number of mosaic ommatidia were found. The advantage of working with these mutants is that no additional parameter has to be introduced: normal cells retain normal levels of N-cadherin synthesis, while for mutant cells there can be no binding strength at all. We first examined an ommatidium with only one mutant C1 cone cell (Fig. 4A). Three of the five cone/cone interfaces are active for homophilic N-cadherin binding, while the two interfaces juxtaposed to the mutant cell are seen to be shorter than normal. An overall symmetry breaking of the ommatidium results. Simulations using the Destruction Model resulted in a pattern that reproduces the overall deformation of the ommatidium (Fig. 4B). In the Recycling Model, a more elaborate distribution of binding strength results. The normal C1 cell has a longer active interface to populate than the C2 cells, and consequently becomes limiting for cadherin distribution. However, the surplus N-cadherin molecules in the C2 cells can form active bonds with each other across the center edge. We thus obtain γN = 0 on the two inactive interfaces, γN≈0. 038 on the two active C1C2 edges and γN≈0. 053 on the center edge, the latter number a much higher binding strength than in the wildtype. The Surface Evolver simulation according to the Recycling Model is shown in Fig. 4C, and is also successful in describing the pattern seen in the experimental image. The differences between Destruction and Recycling simulations are subtle; we quantified the simulation errors using length of the characteristic C1C2 edges on the left side of the ommatidium, together with the angles around them (we use the two tension ratios belonging to the angles P/C1/C2 and P/C2/C1). The definitions of the quantitative individual errors fe and the total error function Fe are found in the Model section. Table 1 shows that the Recycling Model does better in the total error function Fe (it is smaller by about one third – see Table 1), but there is not enough data to make this finding statistically significant. Therefore, we also analyzed an ommatidium that was missing N-cadherin in the bottom C2 and right C1 cells (Fig. 5A). In the Destruction Model, these two cone cells do not alter the concentration of paired N-cadherin along the single active interface, leading to a normal binding strength γN≈0. 032, while all other interfaces between cone cells have γN = 0. Model simulation resulted in a characteristic asymmetry reflecting what is observed in experiment (Fig. 5B). However, the details of the configuration were poorly reproduced. In particular, the active interface is not long enough and the angle under which it meets with the cone/primary pigment cell edges does not fit the data. In the Recycling Model, the entire N-cadherin molecule population becomes distributed along the active interface between the two normal cone cells. Not only does this redistribute the cadherin molecules from other interfaces, but even the number of paired N-cadherin molecules is higher because there are no unpaired molecules left over in the C1 cell. As a result, the N-cadherin binding strength along the active interface is enhanced to βN0/LC1C2. While LC1C2 adjusts itself during a Surface Evolver simulation (and becomes longer because of the energetic advantage of long interfaces carrying large cadherin strength), it is clear that this binding strength is much larger than in the Destruction Model. Indeed, the final value from the Surface Evolver simulation is γN≈0. 081, almost three times the wild-type binding strength (Fig. 5C). The Recycling Model not only simulates the observed asymmetry, but the extreme length and angle associated with the active interface. The error measures in the simulations with the Recycling Model are all significantly smaller than those with the Destruction Model (Table 1); the total error Fe is more than five times smaller. Thus, a prediction of our mechanical energy model is that cone cells use recycling to redistribute N-cadherin along distinct interfaces. Further evidence supports a recycling/redistribution mechanism in cone cells. In the mosaic experiment shown in Fig. 5, the AJs were visualized using an antibody specific for the β-catenin protein (Fig. 5D). This protein stoichiometrically associates with cadherin protein on the cell membrane, where it helps anchor cadherin to the AJ. Importantly, since β-catenin is produced in vast excess and only bound molecules are stable, its abundance is directly proportional to the abundance of cadherin [22]. The fluorescence intensity of β-catenin staining along C1C2 edges of normal ommatidia is significantly lower than the fluorescence along C1C2 edges of mosaic ommatidia with two mutant cone cells, indicating the presence of much more cadherin on the remaining active edges – a finding that is compatible with the Recycling Model (Fig. 5D), but contradictory to the Destruction Model. The mechanical energy model simulations presented so far implicitly assumed that N-cadherin protein expression is established simultaneously in all four cone cells. Moreover, they assumed equal protein synthesis in all cells. While descriptive fluorescence microscopy supports these assumptions [5], it lacks sufficient resolution and quantitativeness to rigorously demonstrate their veracity. Moreover, the four cone cells do not behave identically from a developmental standpoint. C1 cells begin differentiation several hours before C2 cells, and C1 cells specifically induce cells to a P cell fate whereas C2 cells do not [39], [40]. The mechanical energy model is, however, capable of testing these assumptions in certain mutant ommatidia. The Destruction Model for N-cadherin would be insensitive to variability of N-cadherin expression since the eventual local equilibria would be established no matter the sequence or level of cadherin synthesis. However, the Recycling Model would be very sensitive to N-cadherin expression. For example, if the C2 cells were to synthesize N-cadherin earlier than the C1 cells and at a stage where the C2-C2 center edge contact is established, all of the N-cadherin from C2 cells would go onto the center edge, for a strength γNcen = βN0/Lcen, much greater than observed. There would be no N-cadherin left to bind along the C1C2 edges, rendering the other binding energies γNC1C2 = 0. For similar reasons, the Recycling Model would also be sensitive to changes in the level, rather than the timing, of the expression – changes that the Destruction Model, again, would fail to be affected by. We therefore examined situations in which certain cells synthesize N-cadherin at different times and levels from normal, in order to (i) further test the Recycling versus Destruction Models, and (ii) establish in what ways morphology is affected by these changes. We examined mosaic ommatidia in which some cells misexpressed N-cadherin protein. These ommatidia were made by a random genetic switch that activated expression of a N-cadherin transgene in certain cells. Since the switch can be triggered in cells irrespective of their identity, pigment cells can potentially express N-cadherin. Moreover, cone cells affected by the switch express N-cadherin from two sources - the endogenous N-cadherin gene and the transgene i. e. , the total number of N-cadherin molecules is greater than the normal number N0. Production of N-cadherin from the transgene is defined as N+. Switching on of the N-cadherin transgene does not occur synchronously with the onset of endogenous N-cadherin expression but rather occurs 10 to 30 hours earlier in time [5]. Thus, the transgene not only perturbs the level but also the timing of N-cadherin expression (Fig. 6A). We modeled an ommatidium that had N-cadherin trangene expression in two C2 cells and one P cell. With the Recycling Model, the question of which cell is limiting for N-cadherin depends on the ratio N+/N0. However, except for very extreme ratios, the large active interface of the P cell makes it limiting. The binding strength along the P cell is then γNP = βN+/ (2LC2P+LC1P). The next limiting cell is the right-hand C1 cell, which retains N0−LC1PγNP molecules, distributed along the C1C2 edges, where the binding strength is γNC1C2 = β (N0−LC1PγNP) / (2LC1C2). Thus, the overexpressing C2 cells retain a number of N-cadherin molecules given by This number is to be distributed along the left-hand C1C2 edges and the center edge. Depending on whether this number is greater or smaller than N0/2, the left-hand C1 cell or the C2 cells will be limiting. Because the transgene is expressed before the endogenous gene, the P cell initially distributes its N-cadherin along the interfaces contacting the two C2 cells also expressing the transgene, so that γNP = βN+/ (2LC2P). The remaining N+/2 molecules in each C2 cell distribute to the center edge with a binding strength of γNcen, 1 = βN+/ (2Lcen). Thereafter, the endogenous N-cadherin gene is expressed. We assumed that the transgenic N-cadherin would not influence the binding behavior of the endogenous N-cadherin molecules, so that they would distribute similarly (though not identically) to a normal ommatidium. The result of a simulation for N+/N0 = 0. 9 is shown in Fig. 6B. The cadherin binding strengths along different interfaces vary depending on the ratio of N+/N0, with very strong binding along the center edge (Fig. 7A). Qualitatively, the simulation matched an observed experimental ommatidium with the same configuration of misexpression (Fig. 6D). The right-hand C1 cell is bulged out with an elongated C1P edge, while the misexpressing P cell has become more compact and rounded, curving the PP edges at the top and bottom of the ommatidium. We quantified the error in the simulation compared to the observed shape to allow for an unbiased comparison. All of the errors were evaluated as a function of the ratio N+/N0. The overall penalty function showed a well-defined minimum at N+/N0≈1. 25 (Fig. 7B, C). All in all, we evaluated 15 different error contributions (accounting for the left/right asymmetry, there are 6 different tension ratios, 8 different edge lengths, and the parameter Dx – see the Model section for exact definitions of the error quantities). However, we found that only a handful of these measures contribute significantly to the total error (Fig. 7D). We also tested the shapes and penalty function values assuming a Destruction Model for this example, but total error values were very high in this case (data not shown). Thus, the Recycling Model is strongly supported by simulations of N-cadherin misexpression. We determined the sensitivity of the Recycling Model to timing of N-cadherin expression, using the example already described. We reran the simulation with the assumption that endogenous and transgenic N-cadherin were expressed simultaneously. This resulted in a different pattern from observed (Fig. 6C). The penalty function for the simulation was greater than the penalty for the simulation with real timing (Fig. 7B). Thus, the Recycling Model predicts that tissue morphogenesis is highly dependent on timing of cadherin expression even if levels are unchanged. Our Recycling Model illustrates how gene expression established at different times in different cells can lead to dramatically different morphologies. There is another implication to this; cells frequently change neighbors in a prescribed sequence over time during morphogenesis. Thus, even two cells that initiate N-cadherin expression at the same time might not be able to pair their molecules until such time that they become neighbors. This could have profound consequences on the ultimate distribution of N-cadherin and cell morphologies in a recycling scenario. To test this hypothesis, we exploited the known sequence of neighbor changes that P cells undergo (Fig. 8A). P cells first exclusively contact C1 cells; then P cells form additional contact with C2 cells, and only afterwards do they form contact with each other [10]. In normal ommatidia, this has no consequences for N-cadherin binding, as P cells do not contain N-cadherin. We modeled a mosaic ommatidium in which the N-cadherin transgene is switched on in two neighboring P cells. Although this switching occurs at the same time in the two cells, the N-cadherin is unable to pair because the two P cells are not yet neighbors (cf. Fig. 8A). By the time they do become neighbors, the endogenous N-cadherin gene has been expressed in the cone cells for about 10 hours. Therefore, transgenic N-cadherin first pairs on the C1P and C2P edges (Fig. 8A), increasing adhesion and elongating these edges. The binding strength is γNP = βN+/ (2LC2P+LC1P), assuming that the P cells are limiting, which leaves N0−N+ LC1P/ (2LC2P+LC1P) and N0−N+ LC2P/ (2LC2P+LC1P) molecules in the C1 and C2 cells, respectively. The active interface lengths for the two different cone cell species are different (the C2 cells have to cover the center edge as well), but so are the numbers of cadherins to be distributed. The most important feature of the simulation is that, up to N+/N0≈1. 2, the PP edges do not have paired N-cadherin at all, and therefore have a relatively high tension, leading to a characteristic acute P/C2/P angle seen in the simulation (Fig. 8B). Strikingly, an experimental mosaic ommatidium with misexpression in two P neighbors shows a similar morphology (Fig. 8D). A penalty function of the simulation compared to experimental showed a well-defined minimum at N+/N0≈1. 2 (Fig. 9A–C). While the best-fit simulation yielded the cruciform shape with the acute P/C2/P angles seen in experiment, any significant deviation from this optimal N+/N0 led to large shape changes. Note that this optimal ratio is nearly the same as that obtained from the simulation for the misexpression mutant in Fig. 7. The N-cadherin coverages of the various interfaces are plotted in Fig. 9D; note that any abrupt change in slope of the curves indicates a change in the sequence of limiting cells described above. To determine if the simulation worked because of the timing of gene expression with respect to cell contact, we performed a simulation that disregarded the ordered cell contacts and assumed the P cells are in contact when N-cadherin expression begins. Since transgenic N-cadherin is expressed first, it would only be distributed along the PP edges, i. e. , the binding strength there would be a large value γNPP = βN+/ (2LPP), and the effective tension on these edges quite small. The effect of this change produces a simulation that poorly fits the experiment (Fig. 8C). This failure was seen if N-cadherin was expressed either simultaneously or sequentially from the transgene and endogenous gene (Fig. 9A).
  We have applied and extended the mechanical energy functional modeling on the morphology of Drosophila eye tissue in several crucial ways. This model relies entirely on passive global minimization of a simple interfacial energy functional. However, we have shown how to incorporate new features that, implicitly, describe observed biological features of the tissue. These features include: recycling and redistribution of unpaired N-cadherin molecules from the surface of cells, the temporal dynamics of N-cadherin gene expression, and the temporal dynamics of cell-cell contacts. Note that we established these model features without any additional free parameters. Because the dynamics of morphogenesis and gene expression are experimentally described phenomena, and because the Surface Evolver simulations are capable of varying these features, we incorporated these phenomena into a model that describes the experimental data with a high degree of accuracy. When we incorporate features into the model that are not observed in the tissue, the model does not simulate the experiments to the same degree. Two points can be made. First, it is essential for any morphogenesis model to be flexible enough that critical biological features can be modified. Second, our morphogenesis model not only simulates normal morphology but also morphologies observed in complex genetic mutants. In fact, it is the asymmetries in ommatidial shape and the cadherin expression dynamics introduced through mutants that allowed us to successfully test these features within the framework of the model. The success of these simulations suggests that models employing energy minimization can go beyond the quantitative computation of local equilibrium morphologies, and become useful tools for simulating morphogenetic changes that are quasistatic with respect to the time scale of mechanical equilibration. In the case of changes in gene expression or cell neighbor relations, this condition is fulfilled, as these changes can take hours [41], [42], [43]. In contrast, mechanical equilibria are established in minutes at most, depending on dissipative effects [44], [45]. In the experiments we have described, there is also sufficient separation between the onset of cadherin expression and the changes in cell-cell contacts. We simulated misexpression of a N-cadherin transgene to test many of the biological features of the model. In all of these simulations, the best-fit solutions predicted a level of transgenic N-cadherin expression level to be 1. 2-fold greater than that of the endogenous N-cadherin gene. One would expect that independent simulations predict the same level of transgene expression because its expression is independent of which cell or combination of cells express it. Also encouraging is that the level of transgene expression is predicted by simulations to be comparable to the level of endogenous N-cadherin gene expression. The transgene uses the transcriptional promoter from the actin5c gene to indirectly drive transcription of the N-cadherin coding region [5]. Actin5c is one of six actin isoforms and represents only 5–10% of the actin in pupal head cells [46]. Thus, we estimate that 50,000 to 100,000 actin5c protein molecules would be present in these cells, and presumably a similar number of heterologous proteins produced from an actin5c transgene. About 50,000 endogenous cadherin protein molecules are present in a typical animal cell [47]. Therefore, it is reasonable to find that transgenic and endogenous production of N-cadherin are not different by orders of magnitude. Most computer models of tissue morphogenesis simulate visual changes in cell shape and neighbor relations that occur over time [48]. However, use of qualitative visual information to compare normal and perturbed morphogenesis limits their abilities to find optimum mechanical parameters. Our model introduces a novel method for exploring morphogenesis: the effect of perturbations during morphogenesis can be simulated by the final cell shapes and neighbor relations that are achieved at the end of morphogenesis. The simplicity of the final cell geometry allows for unambiguous quantitative comparisons between different morphogenetic conditions, particularly for highly ordered tissues like the Drosophila retinal epithelium.
  In order to have a more general measure for shape differences, we generalized the penalty function of [6], taking into account errors in (i) the interface lengths in the 2D structure, (ii) the angles between interfaces, and (iii) in cases where the entire unit becomes asymmetric, a measure of that global asymmetry. The latter quantity, Dx, is the dimensionless distance of the PP edges to the nearest vertex of the hexagonal ommatidial frame (see Fig. 2; note that in a symmetric ommatidium with centered PP edges Dx = 1/4). The angles between edges (e. g. θj, θk, θl in Fig. 2) are a measure of the differences in the mechanical tensions τj, τk, τl along these edges (if all tensions are equal, all angles are 120·). In fact, the ratio of any two tensions is given by (3) and corresponding cyclic permutation of the indices (jkl). Even subtle differences in angles can translate into strong differences in tension ratios. Therefore, we adopt the quantities ρl instead of the angles to obtain a more sensitive error measure. Any quantity X contributes to the total error via a squared relative error, (4) Here, the superscript SE denotes the quantity obtained from Surface Evolver simulations, while X without superscript is the experimental value. Our total error function is the sum of all fe. Explicitly, (5) where the indices i, j denote cells and the indices l denote angles. The smaller Fe, the better is the agreement between the measured quantities and those obtained by simulations. In many cases, shape differences are dominated by only one or a few fe, while in others the complete set has to be included in Fe.
summary: Tissues are intricate, heterogeneous systems, consisting of individual cells whose shapes and relative positions are of great importance to the tissue' s function, as well as to its formation during morphogenesis. To make progress in our understanding of the formation of organs, their malfunction, and their therapeutic replacement in regenerative medicine, it is crucial to elucidate the connection between shape and function. We have developed a quantitative mechanical model of an epithelial tissue, the retina of Drosophila, and compare the modeling results with experimental data. The model successfully predicts shape changes induced by different expression levels of cell-cell adhesion molecules. Furthermore, the model gives new insight into the changes a tissue undergoes during morphogenesis. Comparing simulations and experiments, we are able to accept or reject different hypotheses about morphogenetic dynamics. In this way, we can identify the time course of adhesion molecule synthesis and of cell-cell contact, as well as gain new insight into the regulation of adhesion strength. Given the prominent role of adhesion in wound healing, cancer research, and many other fields, our fundamental work introduces a novel modeling tool of universal applicability and importance.
keywords: biotechnology
  tissue mechanics
  cadherins
  geometry
  biomechanics
  animal models
  drosophila melanogaster
  model organisms
  mathematics
  cell mechanics
  biophysics simulations
  bioengineering
  cell adhesion
  biophysics theory
  biology
  biophysics
  surface theory
  cell biology
  molecular cell biology
  engineering
title: Cadherin-Dependent Cell Morphology in an Epithelium: Constructing a Quantitative Dynamical Model