## Abstract

How can dense biological tissue maintain sharp boundaries between coexisting cell populations? We explore this question within a simple vertex model for cells, focusing on the role of topology and tissue surface tension. We show that the ability of cells to independently regulate adhesivity and tension, together with neighbor-based interaction rules, lets them support strikingly unusual interfaces. In particular, we show that mechanical- and fluctuation-based measurements of the effective surface tension of a cellular aggregate yield different results, leading to mechanically soft interfaces that are nevertheless extremely sharp.

Original language | English (US) |
---|---|

Article number | 058001 |

Journal | Physical Review Letters |

Volume | 120 |

Issue number | 5 |

DOIs | |

State | Published - Jan 29 2018 |

## ASJC Scopus subject areas

- Physics and Astronomy(all)

## Access to Document

## Other files and links

## Fingerprint Dive into the research topics of 'Soft yet Sharp Interfaces in a Vertex Model of Confluent Tissue'. Together they form a unique fingerprint.

## Cite this

- APA
- Standard
- Harvard
- Vancouver
- Author
- BIBTEX
- RIS

*Physical Review Letters*,

*120*(5), [058001]. https://doi.org/10.1103/PhysRevLett.120.058001

**Soft yet Sharp Interfaces in a Vertex Model of Confluent Tissue.** / Sussman, Daniel M.; Schwarz, J. M.; Marchetti, M. Cristina; Manning, M. Lisa.

Research output: Contribution to journal › Article › peer-review

*Physical Review Letters*, vol. 120, no. 5, 058001. https://doi.org/10.1103/PhysRevLett.120.058001

**Soft yet Sharp Interfaces in a Vertex Model of Confluent Tissue**. In: Physical Review Letters. 2018 ; Vol. 120, No. 5.

}

TY - JOUR

T1 - Soft yet Sharp Interfaces in a Vertex Model of Confluent Tissue

AU - Sussman, Daniel M.

AU - Schwarz, J. M.

AU - Marchetti, M. Cristina

AU - Manning, M. Lisa

N1 - Funding Information: Sussman Daniel M. * Schwarz J. M. Marchetti M. Cristina Manning M. Lisa Department of Physics and Soft Matter Program, Syracuse University , Syracuse, New York 13244, USA * dmsussma@syr.edu 29 January 2018 2 February 2018 120 5 058001 2 October 2017 © 2018 American Physical Society 2018 American Physical Society How can dense biological tissue maintain sharp boundaries between coexisting cell populations? We explore this question within a simple vertex model for cells, focusing on the role of topology and tissue surface tension. We show that the ability of cells to independently regulate adhesivity and tension, together with neighbor-based interaction rules, lets them support strikingly unusual interfaces. In particular, we show that mechanical- and fluctuation-based measurements of the effective surface tension of a cellular aggregate yield different results, leading to mechanically soft interfaces that are nevertheless extremely sharp. National Science Foundation 10.13039/100000001 NSF-POLS-1607416 Simons Foundation 10.13039/100000893 342354 446222 National Science Foundation 10.13039/100000001 DMR-1609208 DMR-1507938 DMR-1352184 Syracuse University 10.13039/100007126 NSF NSF ACI-1541396 Research Corporation for Science Advancement 10.13039/100001309 marker L_SUGG marker SYNOPSIS The process of compartmentalizing different cell populations and maintaining those boundaries is of vital importance in processes ranging from early embryonic development to tumor metastasis [1,2] . A common paradigm, the differential adhesion hypothesis, treats each cell population as an immiscible fluid and suggests that cell sorting and compartmentalization are driven by an effective surface tension [3] , which is in turn governed by a competition between the repulsive and adhesive interactions between cells. The precise cellular mechanisms that govern effective surface tension are still under debate; some investigations suggest it is dominated by adhesive interactions [4] , while others implicate actomyosin contractility [5] or a coregulation of these two effects [1,6,7] . It is not even clear that different methods for measuring the effective surface tension should yield consistent results, which could explain discrepancies between observations and lead to nontrivial and unexpected dynamics for cell sorting and compartmentalization. One hint that something interesting may be happening is a set of experiments demonstrating that many tissues can support extremely sharp boundaries between compartments or coexisting cell populations [1,8–12] . Here we present a possible explanation for these observations based only on the assumption that cells interact mechanically with touching neighbors, and that they might regulate these interactions differently with “unlike” cells. For simplicity, our work focuses on models for single layers of confluent cells, with no cellular gaps or overlaps. 2D vertex models represent confluent monolayers as a polygonal tiling of space where each polygon corresponds to a cell [13–15] ; Voronoi models, which we study here, specialize this idea by taking the cell shapes to be given by a Voronoi tessellation of the cell positions. Vertex and Voronoi models explicitly treat mechanical interactions between neighboring cells, and have successfully been used to model many biophysical processes [12,16–18] , ranging from embryonic development to wound healing to tumor metastasis [19–23] . We include an additional interfacial tension between different cell types to mimic the mechanical changes that are known to occur at so-called “heterotypic” contacts. This extra term naturally leads to a mechanism for robust cell compartmentalization. Surprisingly, we find that this model has a large discrepancy in different macroscopic measurements of the effective surface tension between coexisting cell populations. Specifically, we demonstrate that mechanical measurements and measurements based on the spectrum of interfacial fluctuations, which are equivalent in equilibrium particulate matter, are not equivalent in this system. This difference allows tissues to support mechanically soft interfaces that nevertheless are sharper than a fraction of a cell diameter. This result is a direct consequence of the topological nature of cellular interactions. We begin by writing down a dimensionless form of the vertex model, a commonly used energy functional describing cells in terms of their preferred geometry. Although cells themselves are 3D objects, the 2D vertex model projects these shapes onto the plane and maps the combination of cell volumetric incompressibility and the monolayer’s resistance to height fluctuations into a term quadratic in the cross-sectional area of each cell [19] . e = ∑ i = 1 N [ k A ( a i - a 0 ) 2 + ( p i - p 0 ) 2 ] + ∑ ⟨ i j ⟩ δ [ i ] , [ j ] γ i j l i j . (1) This energy depends on the area a i and perimeter p i of each of the N cells, indexed by i . The unit of length is defined such that the average cell area ⟨ a i ⟩ = 1 . The parameter k A , which we set to unity, controls the cell area stiffness relative to the perimeter stiffness. The preferred values for cell area and perimeter are a 0 and p 0 , respectively. The second sum in Eq. (1) introduces an explicit additional interfacial tension between unlike cells [24] . The sum is over all edges, l i j , between cells i and j ; the delta function is equal to unity if the “type” of cells i and j are different, and zero if they are the same. Biologically, γ i j could arise from spatial reorganization of cadherins and cortical tensions, as well as the generation of acto-myosin cables at cell aggregate boundaries [1] . For simplicity, we assume that the value of the additional tension is the same, γ i j = γ 0 , for all such interfaces. We set a 0 = 1 and choose p 0 to focus on a parameter regime where the bulk has diffusive, fluidlike dynamics ( p 0 > 3.81 ) [25] . To isolate the effect of these heterotypic interfaces, we further assume that “unlike” cells are otherwise identical. There is considerable debate about the dynamical rules that best represent the motion of cells, which are fundamentally out of equilibrium. We emphasize that the unusual interfacial properties we report below do not depend on the precise details of the equations of motion governing the system. To make this point, we carry out simulations [26,27] both in and out of equilibrium. We present equilibrium data using overdamped Brownian dynamics at temperature T . To model out-of-equilibrium dynamics, we adopt those in self-propelled particle models [28] , where each cell tries to move along a polarization vector with self-propulsion speed v 0 . The polarization vector rotates with a diffusion constant D r [24,25,29] . The data we present here are restricted to a modestly out-of-equilibrium regime with 0.025 ≤ v 0 ≤ 0.1 and 0.1 ≤ D r ≤ 10 , and we will define T eff = v 0 2 / 2 D r simply to facilitate comparisons between the dynamical schemes and display them on the same plots. We first probe the effective surface tension of a cell droplet by numerically implementing a parallel-plate compression experiment, a popular technique in biological systems [30,31] . This method measures surface tensions using only the geometry of a deformed droplet and the forces exerted, i.e., without needing to measure fluid viscosities. A schematic diagram is shown in Figs. 1(a) and 1(b) . Although often interpreted in terms of pressure differences across the droplet interface via the Young-Laplace equation [31–33] , the force needed to maintain the plates at fixed spacing can be understood by first assuming an energy given by an effective interfacial tension γ times the perimeter of the droplet, and then taking the derivative of that expression. We fit the droplet shape to an ellipse with major and minor axes given by 2 R 1 and H in Fig. 1(b) , and take an analytic derivative of the perimeter of the ellipse with respect to H . 1 10.1103/PhysRevLett.120.058001.f1 FIG. 1. The effective interfacial tension measured by parallel plate compression matches the imposed microscopic tension. (a) Snapshot of a droplet of cells at the beginning of a parallel plate compression simulation. (b) Geometric quantities used to compute the effective surface tension, as described in the main text or via the Young-Laplace law. (c) Measured surface tension normalized by γ 0 as a function of p 0 . Circles, squares, diamonds, and triangles correspond, respectively, to γ 0 = 0.05 , 0.1, 0.2, and 0.4, with v 0 = 0.1 and D r = 10 . These data are representative of our results for all simulations performed. Inset: Example of the mean force in the y direction on lower plate, shown on a logarithmic scale, as a function of time. A sample initial configuration, showing a droplet of 1250 “type A” cells immersed in fluid of 3400 “type B” cells, is displayed in Fig. 1(a) . The plates, composed of externally forced cells with plate-cell interactions identical to the cell-cell interactions, are moved together at a constant velocity over a brief time window and then held at a fixed distance apart. The transient behavior of the force on the plates during the deformation contains rich physical information, and it is well known that cell aggregates and cellular foams have complicated rheological behavior. Without a detailed viscoelastic model combining viscous effects and plastic cellular rearrangements, one can easily overestimate the value of the surface tension when examining the decay of the transient towards its long-time plateau [34,35] . Thus, we focus on the long-time limit of the external force needed to maintain the droplet in its final ellipsoidal configuration. As seen in Fig. 1(c) , we find an effective surface tension consistent with the value imposed microscopically in the model. We next probe the effective surface tension by measuring the spectrum of fluctuations at an interface. One common method of extracting the surface tension of a fluid phase boundary is to look at the structure and dynamics of capillary waves [36,37] . In equilibrium simulations this is straightforward, as one can directly connect the interfacial roughness to the effective surface tension using standard techniques [36,38] . Assuming that the Fourier spectrum of the interface has independent, Gaussian, equipartitioned modes, each height mode h q has magnitude ⟨ | h q | 2 ⟩ = [ k T / ( γ L y q 2 ) ] , where L y is the linear system size. The size of the periodic box sets the range of accessible wave vectors q = 2 π n / L y for n = 1 to ∞ . The mean interfacial width is w 2 = ∑ q ⟨ | h q | 2 ⟩ , which gives w 2 = w 0 2 + k T L y 12 γ , (2) where w 0 2 captures the contribution of the q = 0 mode. We simulate square domains where L y changes by an order of magnitude, and estimate w for particular configurations by fitting the density profile of the strip to a hyperbolic tangent [36,38] . We then extract the effective surface tension by evaluating the growth of the interfacial width using Eq. (2) . In the cellular models this expected behavior breaks down entirely, both in the scaling with system size and the scaling with wave vector. The majority of our data use the Voronoi model described above, but the open squares in Fig. 2 are simulations of a 2D vertex model; see the Supplemental Material for details [39] . The inset to Fig. 2 shows a sample strip geometry, where visual inspection suggests that the segregation between the cell types is very sharp. This impression is confirmed by Fig. 2 , where we plot the extracted interfacial width as a function of system size normalized by the imposed microscopic tension. For comparison, we also plot a lower bound (with w 0 = 0 ) on the w 2 value expected for an equilibrium system based on our mechanical measurement. Focusing first on the equilibrium system given by the green squares, we see the width is more than an order of magnitude smaller than expected. 2 10.1103/PhysRevLett.120.058001.f2 FIG. 2. Interfaces are much sharper than expected based on the imposed microscopic tension. Normalized growth of the interfacial width w as a function of the size of the periodic simulation domain L y . Solid green squares refer to an equilibrium system at the same T eff as the orange diamonds, which are the self-propelled cell model with v 0 = 0.1 , D r = 1 . Blue circles refer to an active system with v 0 = 0.1 , D r = 0.1 . Open green squares refer to 2D vertex model simulations with the same parameters as the Voronoi simulations represented by solid green squares. The dashed red line is the expected scaling based on the standard capillary wave argument using γ = γ 0 = 0.01 as the interfacial tension and w 0 = 0 . Inset: Sample image of two cell types in a strip geometry. Naively fitting the data would lead to an implied value of γ nearly 2 orders of magnitude larger than γ 0 , but we emphasize that our simulation data do not cleanly fit to the form of Eq. (2) . Simulations with nonequilibrium dynamics are labeled with diamonds and circles. Although these systems are still an order of magnitude sharper than expected, the nonequilibrium dynamics do broaden the interface. This is consistent with observations in self-propelled particle systems [38,43] , and is an interesting avenue for future research. This reduction in the scale of interfacial fluctuations can be understood as a direct consequence of the topological nature of the interaction rules in these cell models. In particular, even though the energy is a continuous function as cell rearrange, there are discontinuous changes in the force whenever cells exchange neighbors. At an interface these discontinuous forces suppress fluctuations, pinning cells to the boundary and sharply compartmentalizing the cells. These forces can be analytically calculated for simplified geometries, as we show in the Supplemental Material for a square lattice of cells [39] . In contrast to homogeneous vertex models, where in the fluid phase fourfold vertices are generically unstable [44] , we predict that in the presence of inhomogeneous line tensions (such as those we consider in this work) some geometric configurations around fourfold vertices are stable. Although not commented upon, fourfold vertices can been seen in, e.g., the vertex-model-simulation images of Ref. [9] . This argument becomes slightly more complicated in the Voronoi model, as at higher values of p 0 the Voronoi constraints can stabilize fourfold vertices even in the bulk [29,45] . Nevertheless, we expect to see an enhancement of fourfold stabilized vertices at the interface in Voronoi models as well. At finite effective temperatures these fourfold vertices will transiently resolve into threefold vertices separated by a short edge. Therefore, we test our predictions by comparing the distribution of edge lengths l i in the bulk and along the interface between two cell types. Figure 3 shows the bulk distribution (black dashed line) and the interfacial distribution for a wide range of imposed tensions (solid lines). For any value of applied tension we see a clear enhancement of very short interfacial edges with respect to the bulk distribution, indicating that the Voronoi tessellation is approximating an excess population of fourfold vertices. As shown in the inset to Fig. 3 , for very short edges we see that the interfacial edge length distribution is approximately exponential, with a decay length l p , consistent with simple energetic considerations based on γ 0 and T eff . 3 10.1103/PhysRevLett.120.058001.f3 FIG. 3. Short interfacial edges become more probable with increasing tension. Probability distribution of interfacial edge lengths, P ( l i ) , for p 0 = 3.95 , v 0 = 0.1 , D r = 1 and γ 0 = 0.005 – 0.64 (light red to dark blue curves). The distribution of cell edge lengths in the bulk is given by the dashed black curve. The l i ≪ 1 distributions are approximately exponential, with characteristic decay length l p . Inset: The length l p , characterizing the short edge distribution for D r = 1 and v 0 = 0.1 , 0.05, 0.025. The dashed red line is a guide to the eye with slope -1, corresponding to a Boltzmann expectation P ( l i ) ∼ exp [ - l i / ( T eff l p ) ] . We next show that the near-fourfold vertices in the disordered geometries in our simulations give rise to discontinuous restoring forces that suppress fluctuations of the interface. We quantify this by systematically displacing each cell at the interface and measuring the restoring force, as described in more detail in the Supplemental Material [39] . This measures the forces in the underlying energy landscape, which we show in the inset to Fig. 4 . At a particulate interface governed by adhesive interactions one would expect a springlike restoring force proportional to the magnitude of the displacement. Instead, we find that the mean force is nearly independent of the displacement over orders of magnitude, and that it is proportional to the applied γ 0 , precisely as predicted by our analysis of the simplified geometry calculated in the Supplemental Material [39] . The main panel of Fig. 4 shows the full distribution of restoring forces for ε x = 10 - 4 , compared to our analytic prediction (dashed lines), indicating that our simple model captures the origin of anomalous behavior in the simulations. 4 10.1103/PhysRevLett.120.058001.f4 FIG. 4. The restoring force changes discontinuously when interfacial cells are displaced. The distribution of restoring forces, P ( f r ) , on cells displaced by ε x = 10 - 4 in the inherent state of a system prepared in a strip geometry with p 0 = 3.95 , D r = 1 , v 0 = 0.1 . From left to right γ 0 = 0.04 , 0.08, 0.16, 0.32, 0.64. The dashed vertical lines show the magnitude of the restoring force predicted from the analytical calculation in the Supplemental Material [39] based on breaking fourfold vertices in a square-lattice geometry for each value of γ 0 (see inset). For each value of γ 0 the mean of the restoring force is nearly independent of the displacement for small displacements, demonstrating the non-Hookean behavior of the interface. The dashed line is the Hookean fit to the γ 0 = 0.04 data at large displacements. How is it possible that a mechanical measurement gives a different answer than the fluctuation-based one? Our works suggests that the nature of the mechanical measurement is important: the parallel plate experiment accesses much larger forces because it is strain controlled, and so it can overcome cusplike pinning forces. The signature of this is likely contained in the transient behavior of the forces shown in the inset of Fig. 1(b) . A small-scale microrheology experiment, for example, might yield a result more similar to the fluctuation-based measurement. We have demonstrated that adding a simple interfacial tension term in a model of topologically interacting cells can lead to highly nontrivial material behavior. We find a strong discrepancy between the effective surface tension defined by mechanical measurements versus those based on fluctuations, even in equilibrium. The roughness of interfaces is almost completely suppressed, leading to strikingly sharp boundaries between fluid domains. The extreme interfacial sharpening is due to the topological nature of the intercellular interactions. Since cells interact with their neighbors, and not according to the distance between cells, a cusplike energy landscape underlies the dynamics. We have confirmed that this behavior is robust to changes to single cell mechanics (e.g., p 0 ), changes to the propulsion mechanism (e.g., thermal vs self-propelled), and changes in the underlying degrees of freedom of the model. Note that in the Voronoi model the topology of the interaction network is set by the geometry, but in models with more degrees of freedom (such as vertex models or cellular Potts models) this is not the case. We believe that exploring the connection between the rules for how cells exchange neighbors and the nature of interfaces is an important avenue of future research. This interfacial sharpening mechanism has obvious implications for cell compartmentalization, but it also may influence the dynamics of cell sorting. It is commonly assumed that both compartmentalizing and sorting proceed as if cells were immiscible fluids; we have seen that many-fold vertices fundamentally alter compartmentalization, and we speculate that they may likewise have profound consequences for the process by which cells sort. Recent work has suggested that surprising consequences can arise from systems interacting via topological interactions rather than purely metric-based ones [46–49] . While much of the work in this direction has focused on explicitly nonequilibrium systems—animal flocking or self-propelled particles interacting with combined Viscek and Voronoi dynamics—here we have shown that surprising interfacial behavior may arise as a generic consequence of the cusplike landscape generated by the topological rules. We speculate that, in the context of real cellular aggregates, epithelial cells may interact topologically whereas mesenchymal or nonconfluent cells may interact metrically (through the surrounding medium or otherwise). In addition to its relevance for confluent cellular aggregates, our findings may point towards a interesting new class of bioinspired materials, where combining mesoscopic interaction units with independent regulation of tension and adhesion may support a diverse set of unusual material properties. This work was primarily supported by NSF-POLS-1607416. Additional support was provided by the Simons Foundation Targeted Grant No. 342354 (M. C. M.) and Investigator Grant No. 446222 (M. L. M.) in the Mathematical Modeling of Living Systems, the National Science Foundation Grants No. DMR-1609208 (M. C. M.), No. DMR-1507938 (J. M. S.), and No. DMR-1352184 (M. L. M.), and a Cottrell Scholar award from the Research Corporation for Science Advancement (M. L. M.). All authors acknowledge support of the Syracuse University Soft Matter Program and computing support through NSF ACI-1541396. The Tesla K40 used for this research was donated by the NVIDIA Corporation. [1] 1 J. D. Amack and M. L. Manning , Science 338 , 212 ( 2012 ). SCIEAS 0036-8075 10.1126/science.1223953 [2] 2 S. Pawlizak , A. W. Fritsch , S. Grosser , D. Ahrens , T. Thalheim , S. Riedel , T. R. Kießling , L. Oswald , M. Zink , M. L. Manning , New J. Phys. 17 , 083049 ( 2015 ). NJOPFM 1367-2630 10.1088/1367-2630/17/8/083049 [3] 3 M. S. Steinberg , Proc. Natl. Acad. Sci. U.S.A. 48 , 1577 ( 1962 ). PNASA6 0027-8424 10.1073/pnas.48.9.1577 [4] 4 M. S. Steinberg , Curr. Opin. Genet. Dev. 17 , 281 ( 2007 ). COGDET 0959-437X 10.1016/j.gde.2007.05.002 [5] 5 A. F. Mertz , S. Banerjee , Y. Che , G. K. German , Y. Xu , C. Hyland , M. C. Marchetti , V. Horsley , and E. R. Dufresne , Phys. Rev. Lett. 108 , 198101 ( 2012 ). PRLTAO 0031-9007 10.1103/PhysRevLett.108.198101 [6] 6 J.-L. Maître , H. Berthoumieux , S. F. G. Krens , G. Salbreux , F. Jülicher , E. Paluch , and C.-P. Heisenberg , Science 338 , 253 ( 2012 ). SCIEAS 0036-8075 10.1126/science.1225399 [7] 7 W. Engl , B. Arasi , L. Yap , J. Thiery , and V. Viasnoff , Nat. Cell Biol. 16 , 584 ( 2014 ). NCBIFN 1465-7392 10.1038/ncb2973 [8] 8 C. Dahmann , A. C. Oates , and M. Brand , Nat. Rev. Genet. 12 , 43 ( 2011 ). 1471-0056 10.1038/nrg2902 [9] 9 K. P. Landsberg , R. Farhadifar , J. Ranft , D. Umetsu , T. J. Widmann , T. Bittig , A. Said , F. Jülicher , and C. Dahmann , Curr. Biol. 19 , 1950 ( 2009 ). CUBLE2 0960-9822 10.1016/j.cub.2009.10.021 [10] 10 B. Monier , A. Pélissier-Monier , and B. Sanson , Cell Mol. Life Sci. 68 , 1897 ( 2011 ). CMLSFI 1420-9071 10.1007/s00018-011-0668-8 [11] 11 K. D. Nnetu , M. Knorr , J. Käs , and M. Zink , New J. Phys. 14 , 115012 ( 2012 ). NJOPFM 1367-2630 10.1088/1367-2630/14/11/115012 [12] 12 M. L. Manning , R. A. Foty , M. S. Steinberg , and E.-M. Schoetz , Proc. Natl. Acad. Sci. U.S.A. 107 , 12517 ( 2010 ). PNASA6 0027-8424 10.1073/pnas.1003743107 [13] 13 H. Honda , J. Theor. Biol. 72 , 523 ( 1978 ). JTBIAP 0022-5193 10.1016/0022-5193(78)90315-6 [14] 14 A. Hočevar and P. Ziherl , Phys. Rev. E 80 , 011904 ( 2009 ). PRESCM 1539-3755 10.1103/PhysRevE.80.011904 [15] 15 T. Vazquez-Faci , R. van Drongelen , M. van der Zee , and T. Idema , arXiv:1705.06205 . [16] 16 G. W. Brodland , Appl. Mech. Rev. 57 , 47 ( 2004 ). AMREAD 0003-6900 10.1115/1.1583758 [17] 17 T. Nagai and H. Honda , Philos. Mag. Part B 81 , 699 ( 2001 ). PMABDJ 1364-2812 10.1080/13642810108205772 [18] 18 D. Bi , J. H. Lopez , J. M. Schwarz , and M. L. Manning , Soft Matter 10 , 1885 ( 2014 ). SMOABF 1744-683X 10.1039/c3sm52893f [19] 19 L. Hufnagel , A. A. Teleman , H. Rouault , S. M. Cohen , and B. I. Shraiman , Proc. Natl. Acad. Sci. U.S.A. 104 , 3835 ( 2007 ). PNASA6 0027-8424 10.1073/pnas.0607134104 [20] 20 R. Farhadifar , J.-C. Röper , B. Aigouy , S. Eaton , and F. Jülicher , Curr. Biol. 17 , 2095 ( 2007 ). CUBLE2 0960-9822 10.1016/j.cub.2007.11.049 [21] 21 P. Friedl and D. Gilmour , Nat. Rev. Mol. Cell Biol. 10 , 445 ( 2009 ). NRMCBP 1471-0072 10.1038/nrm2720 [22] 22 A. Brugués , E. Anon , V. Conte , J. H. Veldhuis , M. Gupta , J. Colombelli , J. J. Muñoz , G. W. Brodland , B. Ladoux , and X. Trepat , Nat. Phys. 10 , 683 ( 2014 ). NPAHAX 1745-2473 10.1038/nphys3040 [23] 23 R. Etournay , M. Popović , M. Merkel , A. Nandi , C. Blasse , B. Aigouy , H. Brandl , G. Myers , G. Salbreux , F. Jülicher , and S. Eaton , eLife 4 , e07090 ( 2015 ). ELIFA8 2050-084X 10.7554/eLife.07090 [24] 24 D. L. Barton , S. Henkes , C. J. Weijer , and R. Sknepnek , PLoS Comput. Biol. 13 , e1005569 ( 2017 ). PCBLBG 1553-7358 10.1371/journal.pcbi.1005569 [25] 25 D. Bi , X. Yang , M. C. Marchetti , and M. L. Manning , Phys. Rev. X 6 , 021011 ( 2016 ). PRXHAE 2160-3308 10.1103/PhysRevX.6.021011 [26] 26 D. M. Sussman , Comput. Phys. Commun. 219 , 400 ( 2017 ). CPHCBZ 0010-4655 10.1016/j.cpc.2017.06.001 [27] 27 https://gitlab.com/dmsussman/cellGPU , 2017 . [28] 28 S. Henkes , Y. Fily , and M. C. Marchetti , Phys. Rev. E 84 , 040301 ( 2011 ). PRESCM 1539-3755 10.1103/PhysRevE.84.040301 [29] 29 D. M. Sussman and M. Merkel , arXiv:1708.03396 . [30] 30 R. A. Foty , G. Forgacs , C. M. Pfleger , and M. S. Steinberg , Phys. Rev. Lett. 72 , 2298 ( 1994 ). PRLTAO 0031-9007 10.1103/PhysRevLett.72.2298 [31] 31 A. Mgharbel , H. Delanoë-Ayari , and J.-P. Rieu , HFSP J. 3 , 213 ( 2009 ). HJFOA5 1955-2068 10.2976/1.3116822 [32] 32 T. Young , Phil. Trans. R. Soc. London 95 , 65 ( 1805 ). PTRSAV 0370-2316 10.1098/rstl.1805.0005 [33] 33 P.-S. Laplace , Courcier, Paris ( 1806 ). [34] 34 P. Marmottant , A. Mgharbel , J. Käfer , B. Audren , J.-P. Rieu , J.-C. Vial , B. Van Der Sanden , A. F. Marée , F. Graner , and H. Delanoë-Ayari , Proc. Natl. Acad. Sci. U.S.A. 106 , 17271 ( 2009 ). PNASA6 0027-8424 10.1073/pnas.0902085106 [35] 35 P. Marmottant and F. Graner , Soft Matter 9 , 9602 ( 2013 ). SMOABF 1744-683X 10.1039/c3sm51220g [36] 36 J. S. Rowlinson and B. Widom , Molecular Theory of Capillarity ( Courier Corporation , New York, 2013 ). [37] 37 M. P. Gelfand and M. E. Fisher , Physica (Amsterdam) 166A , 1 ( 1990 ). PHYADX 0378-4371 10.1016/0378-4371(90)90099-E [38] 38 J. Bialké , J. T. Siebert , H. Löwen , and T. Speck , Phys. Rev. Lett. 115 , 098301 ( 2015 ). PRLTAO 0031-9007 10.1103/PhysRevLett.115.098301 [39] 39 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.120.058001 for additional details, which includes Refs. [40–42]. [40] 40 A. G. Fletcher , J. M. Osborne , P. K. Maini , and D. J. Gavaghan , Prog. Biophys. Molec. Biol. 113 , 299 ( 2013 ). PBIMAC 0079-6107 10.1016/j.pbiomolbio.2013.09.003 [41] 41 A. G. Fletcher , M. Osterfield , R. E. Baker , and S. Y. Shvartsman , Biophys. J. 106 , 2291 ( 2014 ). BIOJAU 0006-3495 10.1016/j.bpj.2013.11.4498 [42] 42 E. Bitzek , P. Koskinen , F. Gähler , M. Moseler , and P. Gumbsch , Phys. Rev. Lett. 97 , 170201 ( 2006 ). PRLTAO 0031-9007 10.1103/PhysRevLett.97.170201 [43] 43 T. Speck , Europhys. Lett. 114 , 30006 ( 2016 ). EULEEJ 0295-5075 10.1209/0295-5075/114/30006 [44] 44 M. A. Spencer , Z. Jabeen , and D. K. Lubensky , Eur. Phys. J. E 40 , 2 ( 2017 ). EPJSFH 1292-8941 10.1140/epje/i2017-11489-4 [45] 45 T. Su and G. Lan , arXiv:1610.04254 . [46] 46 M. Ballerini , N. Cabibbo , R. Candelier , A. Cavagna , E. Cisbani , I. Giardina , V. Lecomte , A. Orlandi , G. Parisi , A. Procaccini , Proc. Natl. Acad. Sci. U.S.A. 105 , 1232 ( 2008 ). PNASA6 0027-8424 10.1073/pnas.0711437105 [47] 47 F. Ginelli and H. Chaté , Phys. Rev. Lett. 105 , 168103 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.105.168103 [48] 48 J. E. Herbert-Read , A. Perna , R. P. Mann , T. M. Schaerf , D. J. Sumpter , and A. J. Ward , Proc. Natl. Acad. Sci. U.S.A. 108 , 18726 ( 2011 ). PNASA6 0027-8424 10.1073/pnas.1109355108 [49] 49 A. Peshkov , S. Ngo , E. Bertin , H. Chaté , and F. Ginelli , Phys. Rev. Lett. 109 , 098101 ( 2012 ). PRLTAO 0031-9007 10.1103/PhysRevLett.109.098101

PY - 2018/1/29

Y1 - 2018/1/29

N2 - How can dense biological tissue maintain sharp boundaries between coexisting cell populations? We explore this question within a simple vertex model for cells, focusing on the role of topology and tissue surface tension. We show that the ability of cells to independently regulate adhesivity and tension, together with neighbor-based interaction rules, lets them support strikingly unusual interfaces. In particular, we show that mechanical- and fluctuation-based measurements of the effective surface tension of a cellular aggregate yield different results, leading to mechanically soft interfaces that are nevertheless extremely sharp.

AB - How can dense biological tissue maintain sharp boundaries between coexisting cell populations? We explore this question within a simple vertex model for cells, focusing on the role of topology and tissue surface tension. We show that the ability of cells to independently regulate adhesivity and tension, together with neighbor-based interaction rules, lets them support strikingly unusual interfaces. In particular, we show that mechanical- and fluctuation-based measurements of the effective surface tension of a cellular aggregate yield different results, leading to mechanically soft interfaces that are nevertheless extremely sharp.

UR - http://www.scopus.com/inward/record.url?scp=85041294281&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85041294281&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.120.058001

DO - 10.1103/PhysRevLett.120.058001

M3 - Article

C2 - 29481188

AN - SCOPUS:85041294281

VL - 120

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 5

M1 - 058001

ER -