Chapter 4 Cell Topology, Geometry, and Morphogenesis in Proliferating Epithelia

Shaping epithelial tissues by stem cell mechanics in development and cancer

Adult stem cells balance self-renewal and differentiation to build, maintain and repair tissues. The role of signalling pathways and transcriptional networks in controlling stem cell function has been extensively studied, but there is increasing appreciation that mechanical forces also have a crucial regulatory role. Mechanical forces, signalling pathways and transcriptional networks must be coordinated across diverse length and timescales to maintain tissue homeostasis and function. Such coordination between stem cells and neighbouring cells dictates when cells divide, migrate and differentiate. Recent advances in measuring and manipulating the mechanical forces that act upon and are produced by stem cells are providing new insights into development and disease. In this Review, we discuss the mechanical forces involved when epithelial stem cells construct their microenvironment and what happens in cancer when stem cell niche mechanics are disrupted or dysregulated. As the skin has evolved to withstand the harsh mechanical pressures from the outside environment, we often use the stem cells of mammalian skin epithelium as a paradigm for adult stem cells shaping their surrounding tissues.

Spatially Resolved Single-Cell Morphometry of Benign Breast Disease Biopsy Images Uncovers Quantitative Cytomorphometric Features Predictive of Subsequent Invasive Breast Cancer Risk

Currently, benign breast disease (BBD) pathologic classification and invasive breast cancer (BC) risk assessment are based on qualitative epithelial changes, with limited utility for BC risk stratification for women with lower-risk category BBD (ie, nonproliferative disease [NPD] and proliferative disease without atypia [PDWA]). Here, machine learning-based single-cell morphometry was used to characterize quantitative changes in epithelial nuclear morphology that reflect functional/structural decline (ie, increasing nuclear size, assessed as epithelial nuclear area and nuclear perimeter), altered DNA chromatin content (ie, increasing nuclear chromasia), and increased cellular crowding/proliferation (ie, increasing nuclear contour irregularity). Cytomorphologic changes reflecting chronic stromal inflammation were assessed using stromal cellular density. Data and pathology materials were obtained from a case-control study (n = 972) nested within a cohort of 15,395 women diagnosed with BBD at Kaiser Permanente Northwest (1971-2012). Odds ratios (ORs) and 95% confidence intervals (CIs) for associations of cytomorphometric features with risk of subsequent BC were assessed using multivariable logistic regression. More than 55 million epithelial and 37 million stromal cells were profiled across 972 BBD images. Cytomorphometric features were individually predictive of subsequent BC risk, independently of BBD histologic classification. However, cytomorphometric features of epithelial functional/structural decline were statistically significantly predictive of low-grade but not high-grade BC following PDWA (OR for low-grade BC per 1-SD increase in nuclear area and nuclear perimeter, 2.10; 95% CI, 1.26-3.49, and 2.22; 95% CI, 1.30-3.78, respectively), whereas stromal inflammation was predictive of high-grade but not low-grade BC following NPD (OR for high-grade BC per 1-SD increase in stromal cellular density, 1.53; 95% CI, 1.13-2.08). Associations of nuclear chromasia and nuclear contour irregularity with subsequent tumor grade were context specific, with both features predicting low-grade BC risk following PDWA (OR per 1-SD, 1.58; 95% CI, 1.06-2.35, and 2.21; 95% CI, 1.25-3.91, for nuclear chromasia and nuclear contour irregularity, respectively) and high-grade BC following NPD (OR per 1-SD, 1.47; 95% CI, 1.11-1.96, and 1.29; 95% CI, 1.00-1.70, for nuclear chromasia and nuclear contour irregularity, respectively). The results indicate that cytomorphometric features on BBD hematoxylin-eosin-stained images might help to refine BC risk estimation and potentially inform BC risk reduction strategies for BBD patients, particularly those currently designated as low risk.

Nanomaterial application for protein delivery in bone regeneration therapy

Bone fractures must undergo a complex healing process involving intricate cellular and molecular mechanisms. They require a suitable biological environment to restore skeletal stability and resolve inflammation. Scaffolds play a vital role in bone regeneration, thus reducing disease burden. Autologous bone graft represents the gold standard of therapy. However, its application is limited due to various reasons. Nanotechnology, in the form of nanomaterials and nano-drug delivery systems, has been proven to increase the potency of active substances in mimicking extracellular matrix (ECM), thereby providing physical support benefits and enhancing therapeutic effectiveness. Various materials, including protein, metal oxide, hydroxyapatite, and silica are modified with nanoparticle technology for the purposes of tissue regeneration therapy. Moreover, the properties of nanomaterials such as size, seta potential, and surface properties will affect their effectiveness in bone regeneration therapy. This review provides insights that deepen the knowledge of the manufacturing and application of nanomaterials as a therapeutic agent for bone regeneration.

YAP/TAZ-associated cell signaling - at the crossroads of cancer and neurodevelopmental disorders

YAP/TAZ (Yes-associated protein/paralog transcriptional co-activator with PDZ-binding domain) are transcriptional cofactors that are the key and major downstream effectors of the Hippo signaling pathway. Both are known to play a crucial role in defining cellular outcomes, including cell differentiation, cell proliferation, and apoptosis. Aside from the canonical Hippo signaling cascade with the key components MST1/2 (mammalian STE20-like kinase 1/2), SAV1 (Salvador homologue 1), MOB1A/B (Mps one binder kinase activator 1A/B) and LATS1/2 (large tumor suppressor kinase 1/2) upstream of YAP/TAZ, YAP/TAZ activation is also influenced by numerous other signaling pathways. Such non-canonical regulation of YAP/TAZ includes well-known growth factor signaling pathways such as the epidermal growth factor receptor (EGFR)/ErbB family, Notch, and Wnt signaling as well as cell-cell adhesion, cell-matrix interactions and mechanical cues from a cell's microenvironment. This puts YAP/TAZ at the center of a complex signaling network capable of regulating developmental processes and tissue regeneration. On the other hand, dysregulation of YAP/TAZ signaling has been implicated in numerous diseases including various cancers and neurodevelopmental disorders. Indeed, in recent years, parallels between cancer development and neurodevelopmental disorders have become apparent with YAP/TAZ signaling being one of these pathways. This review discusses the role of YAP/TAZ in brain development, cancer and neurodevelopmental disorders with a special focus on the interconnection in the role of YAP/TAZ in these different conditions.

The Role of Cytonemes and Diffusive Transport in the Establishment of Morphogen Gradients

Spatial distributions of morphogens provide positional information in developing systems, but how the distributions are established and maintained remains an open problem. Transport by diffusion has been the traditional mechanism, but recent experimental work has shown that cells can also communicate by filopodia-like structures called cytonemes that make direct cell-to-cell contacts. Here we investigate the roles each may play individually in a complex tissue and how they can jointly establish a reliable spatial distribution of a morphogen. To this end, we formulate models that capture fundamental aspects of various cytoneme-based transport mechanisms. In simple cases, exact solutions are attainable, and in more complex cases, we discuss results of numerical simulations.

Piezo regulates epithelial topology and promotes precision in organ size control

Mechanosensitive Piezo channels regulate cell division, cell extrusion, and cell death. However, systems-level functions of Piezo in regulating organogenesis remain poorly understood. Here, we demonstrate that Piezo controls epithelial cell topology to ensure precise organ growth by integrating live-imaging experiments with pharmacological and genetic perturbations and computational modeling. Notably, the knockout or knockdown of Piezo increases bilateral asymmetry in wing size. Piezo's multifaceted functions can be deconstructed as either autonomous or non-autonomous based on a comparison between tissue-compartment-level perturbations or between genetic perturbation populations at the whole-tissue level. A computational model that posits cell proliferation and apoptosis regulation through modulation of the cutoff tension required for Piezo channel activation explains key cell and tissue phenotypes arising from perturbations of Piezo expression levels. Our findings demonstrate that Piezo promotes robustness in regulating epithelial topology and is necessary for precise organ size control.

A mechanical transition from tension to buckling underlies the jigsaw puzzle shape morphogenesis of histoblasts in the Drosophila epidermis

The polygonal shape of cells in proliferating epithelia is a result of the tensile forces of the cytoskeletal cortex and packing geometry set by the cell cycle. In the larval Drosophila epidermis, two cell populations, histoblasts and larval epithelial cells, compete for space as they grow on a limited body surface. They do so in the absence of cell divisions. We report a striking morphological transition of histoblasts during larval development, where they change from a tensed network configuration with straight cell outlines at the level of adherens junctions to a highly folded morphology. The apical surface of histoblasts shrinks while their growing adherens junctions fold, forming deep lobules. Volume increase of growing histoblasts is accommodated basally, compensating for the shrinking apical area. The folded geometry of apical junctions resembles elastic buckling, and we show that the imbalance between the shrinkage of the apical domain of histoblasts and the continuous growth of junctions triggers buckling. Our model is supported by laser dissections and optical tweezer experiments together with computer simulations. Our analysis pinpoints the ability of histoblasts to store mechanical energy to a much greater extent than most other epithelial cell types investigated so far, while retaining the ability to dissipate stress on the hours time scale. Finally, we propose a possible mechanism for size regulation of histoblast apical size through the lateral pressure of the epidermis, driven by the growth of cells on a limited surface. Buckling effectively compacts histoblasts at their apical plane and may serve to avoid physical harm to these adult epidermis precursors during larval life. Our work indicates that in growing nondividing cells, compressive forces, instead of tension, may drive cell morphology.

Evolutionary dynamics of mutants that modify population structure

Natural selection is usually studied between mutants that differ in reproductive rate, but are subject to the same population structure. Here we explore how natural selection acts on mutants that have the same reproductive rate, but different population structures. In our framework, population structure is given by a graph that specifies where offspring can disperse. The invading mutant disperses offspring on a different graph than the resident wild-type. We find that more densely connected dispersal graphs tend to increase the invader's fixation probability, but the exact relationship between structure and fixation probability is subtle. We present three main results. First, we prove that if both invader and resident are on complete dispersal graphs, then removing a single edge in the invader's dispersal graph reduces its fixation probability. Second, we show that for certain island models higher invader's connectivity increases its fixation probability, but the magnitude of the effect depends on the exact layout of the connections. Third, we show that for lattices the effect of different connectivity is comparable to that of different fitness: for large population size, the invader's fixation probability is either constant or exponentially small, depending on whether it is more or less connected than the resident.

Building block aspect ratio controls assembly, architecture, and mechanics of synthetic and natural protein networks

Fibrous networks constructed from high aspect ratio protein building blocks are ubiquitous in nature. Despite this ubiquity, the functional advantage of such building blocks over globular proteins is not understood. To answer this question, we engineered hydrogel network building blocks with varying numbers of protein L domains to control the aspect ratio. The mechanical and structural properties of photochemically crosslinked protein L networks were then characterised using shear rheology and small angle neutron scattering. We show that aspect ratio is a crucial property that defines network architecture and mechanics, by shifting the formation from translationally diffusion dominated to rotationally diffusion dominated. Additionally, we demonstrate that a similar transition is observed in the model living system: fibrin blood clot networks. The functional advantages of this transition are increased mechanical strength and the rapid assembly of homogenous networks above a critical protein concentration, crucial for in vivo biological processes such as blood clotting. In addition, manipulating aspect ratio also provides a parameter in the design of future bio-mimetic and bio-inspired materials.

Deciphering the Role of p60AmotL2 in Epithelial Extrusion and Cell Detachment

Preserving an accurate cell count is crucial for maintaining homeostasis. Apical extrusion, a process in which redundant cells are eliminated by neighboring cells, plays a key role in this regard. Recent studies have revealed that apical extrusion can also be triggered in cells transformed by oncogenes, suggesting it may be a mechanism through which tumor cells escape their microenvironment. In previous work, we demonstrated that p60AmotL2 modulates the E-cadherin function by inhibiting its connection to radial actin filaments. This isoform of AmotL2 is expressed in invasive breast and colon tumors and promotes invasion in vitro and in vivo. Transcriptionally regulated by c-Fos, p60AmotL2 is induced by local stress signals such as severe hypoxia. In this study, we investigated the normal role of p60AmotL2 in epithelial tissues. We found that this isoform is predominantly expressed in the gut, where cells experience rapid turnover. Through time-lapse imaging, we present evidence that cells expressing p60AmotL2 are extruded by their normal neighboring cells. Based on these findings, we hypothesize that tumor cells exploit this pathway to detach from normal epithelia and invade surrounding tissues.

Spatial topology of organelle is a new breast cancer cell classifier

Genomics and proteomics have been central to identify tumor cell populations, but more accurate approaches to classify cell subtypes are still lacking. We propose a new methodology to accurately classify cancer cells based on their organelle spatial topology. Herein, we developed an organelle topology-based cell classification pipeline (OTCCP), which integrates artificial intelligence (AI) and imaging quantification to analyze organelle spatial distribution and inter-organelle topology. OTCCP was used to classify a panel of human breast cancer cells, grown as 2D monolayer or 3D tumor spheroids using early endosomes, mitochondria, and their inter-organelle contacts. Organelle topology allows for a highly precise differentiation between cell lines of different subtypes and aggressiveness. These findings lay the groundwork for using organelle topological profiling as a fast and efficient method for phenotyping breast cancer function as well as a discovery tool to advance our understanding of cancer cell biology at the subcellular level.

Beyond the Genetic Code: A Tissue Code?

The genetic code determines how the precise amino acid sequence of proteins is specified by genomic information in cells. But what specifies the precise histologic organization of cells in plant and animal tissues is unclear. We now hypothesize that another code, the tissue code , exists at an even higher level of complexity which determines how tissue organization is dynamically maintained. Accordingly, we modeled spatial and temporal asymmetries of cell division and established that five simple mathematical laws ("the tissue code") convey a set of biological rules that maintain the specific organization and continuous self-renewal dynamics of cells in tissues. These laws might even help us understand wound healing, and how tissue disorganization leads to birth defects and tissue pathology like cancer.

Entropy as a Geometrical Source of Information in Biological Organizations

Considering both biological and non-biological polygonal shape organizations, in this paper we introduce a quantitative method which is able to determine informational entropy as spatial differences between heterogeneity of internal areas from simulation and experimental samples. According to these data (i.e., heterogeneity), we are able to establish levels of informational entropy using statistical insights of spatial orders using discrete and continuous values. Given a particular state of entropy, we establish levels of information as a novel approach which can unveil general principles of biological organization. Thirty-five geometric aggregates are tested (biological, non-biological, and polygonal simulations) in order to obtain the theoretical and experimental results of their spatial heterogeneity. Geometrical aggregates (meshes) include a spectrum of organizations ranging from cell meshes to ecological patterns. Experimental results for discrete entropy using a bin width of 0.5 show that a particular range of informational entropy (0.08 to 0.27 bits) is intrinsically associated with low rates of heterogeneity, which indicates a high degree of uncertainty in finding non-homogeneous configurations. In contrast, differential entropy (continuous) results reflect negative entropy within a particular range (−0.4 to −0.9) for all bin widths. We conclude that the differential entropy of geometrical organizations is an important source of neglected information in biological systems.

Cell Architecture-Dependent Constraints: Critical Safeguards to Carcinogenesis

Animal cells display great diversity in their shape. These morphological characteristics result from crosstalk between the plasma membrane and the force-generating capacities of the cytoskeleton macromolecules. Changes in cell shape are not merely byproducts of cell fate determinants, they also actively drive cell fate decisions, including proliferation and differentiation. Global and local changes in cell shape alter the transcriptional program by a multitude of mechanisms, including the regulation of physical links between the plasma membrane and the nuclear envelope and the mechanical modulation of cation channels and signalling molecules. It is therefore not surprising that anomalies in cell shape contribute to several diseases, including cancer. In this review, we discuss the possibility that the constraints imposed by cell shape determine the behaviour of normal and pro-tumour cells by organizing the whole interconnected regulatory network. In turn, cell behaviour might stabilize cells into discrete shapes. However, to progress towards a fully transformed phenotype and to acquire plasticity properties, pro-tumour cells might need to escape these cell shape restrictions. Thus, robust controls of the cell shape machinery may represent a critical safeguard against carcinogenesis.

Contact-mediated signaling enables disorder-driven transitions in cellular assemblies

We show that, when cells communicate by contact-mediated interactions, heterogeneity in cell shapes and sizes leads to qualitatively distinct collective behavior in the tissue. For intercellular coupling that implements lateral inhibition, such disorder-driven transitions can substantially alter the asymptotic pattern of differentiated cells by modulating their fate choice through changes in the neighborhood geometry. In addition, when contact-induced signals influence inherent cellular oscillations, disorder leads to the emergence of functionally relevant partially-ordered dynamical states.

Random nature of epithelial cancer cell monolayers

Although the polygonal shape of epithelial cells has been drawing the attention of scientists for several centuries, only a decade and a half ago it was demonstrated that distributions of polygon types (DOPTs) are similar in proliferative epithelia of many different plant and animal species. In this study, we show that hyper-proliferation of cancer cells disrupts this universal paradigm and results in randomly organized epithelial structures. Examining non-synchronized and synchronized HeLa cervix cells, we suppose that the spread of cell sizes is the main parameter controlling the DOPT in the cancer cell monolayers. To test this hypothesis, we develop a theory of morphologically similar random polygonal packings. By analysing differences between tumoural and normal epithelial cell monolayers, we conclude that the latter have more ordered structures because of their lower proliferation rates and, consequently, more effective relaxation of mechanical stress associated with cell division and growth. To explain the structural features of normal proliferative epithelium, we take into account the spread of cell sizes in the monolayer. The proposed theory also rationalizes some highly ordered unconventional post-mitotic epithelia.

A local orientational order parameter for systems of interacting particles

Many physical systems are well modeled as collections of interacting particles. Nevertheless, a general approach to quantifying the absolute degree of order immediately surrounding a particle has yet to be described. Motivated thus, we introduce a quantity E that captures the amount of pairwise informational redundancy among the bonds formed by a particle. Particles with larger E have less diversity in bond angles and thus simpler neighborhoods. We show that E possesses a number of intuitive mathematical properties, such as increasing monotonicity in the coordination number of Platonic polyhedral geometries. We demonstrate analytically that E is, in principle, able to distinguish a wide range of structures and conjecture that it is maximized by the icosahedral geometry under the constraint of equal sphere packing. An algorithm for computing E is described and is applied to the structural characterization of crystals and glasses. The findings of this study are generally consistent with existing knowledge on the structure of such systems. We compare E to the Steinhardt order parameter Q6 and polyhedral template matching (PTM). We observe that E has resolution comparable to Q6 and robustness similar to PTM despite being much simpler than the former and far more informative than the latter.

Grhl3 promotes retention of epidermal cells under endocytic stress to maintain epidermal architecture in zebrafish

Epithelia such as epidermis cover large surfaces and are crucial for survival. Maintenance of tissue homeostasis by balancing cell proliferation, cell size, and cell extrusion ensures epidermal integrity. Although the mechanisms of cell extrusion are better understood, how epithelial cells that round up under developmental or perturbed genetic conditions are reintegrated in the epithelium to maintain homeostasis remains unclear. Here, we performed live imaging in zebrafish embryos to show that epidermal cells that round up due to membrane homeostasis defects in the absence of goosepimples/myosinVb (myoVb) function, are reintegrated into the epithelium. Transcriptome analysis and genetic interaction studies suggest that the transcription factor Grainyhead-like 3 (Grhl3) induces the retention of rounded cells by regulating E-cadherin levels. Moreover, Grhl3 facilitates the survival of MyoVb deficient embryos by regulating cell adhesion, cell retention, and epidermal architecture. Our analyses have unraveled a mechanism of retention of rounded cells and its importance in epithelial homeostasis.

A Random Walk Approach to Transport in Tissues and Complex Media: From Microscale Descriptions to Macroscale Models

The biological processes necessary for the development and continued survival of any organism are often strongly influenced by the transport properties of various biologically active species. The transport phenomena involved vary over multiple temporal and spatial scales, from organism-level behaviors such as the search for food, to systemic processes such as the transport of oxygen from the lungs to distant organs, down to microscopic phenomena such as the stochastic movement of proteins in a cell. Each of these processes is influenced by many interrelated factors. Identifying which factors are the most important, and how they interact to determine the overall result is a problem of great importance and interest. Experimental observations are often fit to relatively simple models, but in reality the observations are the output of complicated functions of the physicochemical, topological, and geometrical properties of a given system. Herein we use multistate continuous-time random walks and generalized master equations to model transport processes involving spatial jumps, immobilization at defined sites, and stochastic internal state changes. The underlying spatial models, which are framed as graphs, may have different classes of nodes, and walkers may have internal states that are governed by a Markov process. A general form of the solutions, using Fourier-Laplace transforms and asymptotic analysis, is developed for several spatially infinite regular lattices in one and two spatial dimensions, and the theory is developed for the analysis of transport and internal state changes on general graphs. The goal in each case is to shed light on how experimentally observable macroscale transport coefficients can be explained in terms of microscale properties of the underlying processes. This work is motivated by problems arising in transport in biological tissues, but the results are applicable to a broad class of problems that arise in other applications.

DeepScratch: Single-cell based topological metrics of scratch wound assays

Changes in tissue architecture and multicellular organisation contribute to many diseases, including cancer and cardiovascular diseases. Scratch wound assay is a commonly used tool that assesses cells' migratory ability based on the area of a wound they cover over a certain time. However, analysis of changes in the organisational patterns formed by migrating cells following genetic or pharmacological perturbations are not well explored in these assays, in part because analysing the resulting imaging data is challenging. Here we present DeepScratch, a neural network that accurately detects the cells in scratch assays based on a heterogeneous set of markers. We demonstrate the utility of DeepScratch by analysing images of more than 232,000 lymphatic endothelial cells. In addition, we propose various topological measures of cell connectivity and local cell density (LCD) to characterise tissue remodelling during wound healing. We show that LCD-based metrics allow classification of CDH5 and CDC42 genetic perturbations that are known to affect cell migration through different biological mechanisms. Such differences cannot be captured when considering only the wound area. Taken together, single-cell detection using DeepScratch allows more detailed investigation of the roles of various genetic components in tissue topology and the biological mechanisms underlying their effects on collective cell migration.

Functionalizable Antifouling Coatings as Tunable Platforms for the Stress-Driven Manipulation of Living Cell Machinery

Cells are continuously sensing their microenvironment and subsequently respond to different physicochemical cues by the activation or inhibition of different signaling pathways. To study a very complex cellular response, it is necessary to diminish background environmental influences and highlight the particular event. However, surface-driven nonspecific interactions of the abundant biomolecules from the environment influence the targeted cell response significantly. Yes-associated protein (YAP) translocation may serve as a marker of human hepatocellular carcinoma (Huh7) cell responses to the extracellular matrix and surface-mediated stresses. Here, we propose a platform of tunable functionable antifouling poly(carboxybetain) (pCB)-based brushes to achieve a molecularly clean background for studying arginine, glycine, and aspartic acid (RGD)-induced YAP-connected mechanotransduction. Using two different sets of RGD-functionalized zwitterionic antifouling coatings with varying compositions of the antifouling layer, a clear correlation of YAP distribution with RGD functionalization concentrations was observed. On the other hand, commonly used surface passivation by the oligo(ethylene glycol)-based self-assembled monolayer (SAM) shows no potential to induce dependency of the YAP distribution on RGD concentrations. The results indicate that the antifouling background is a crucial component of surface-based cellular response studies, and pCB-based zwitterionic antifouling brush architectures may serve as a potential next-generation easily functionable surface platform for the monitoring and quantification of cellular processes.

Isotropic myosin-generated tissue tension is required for the dynamic orientation of the mitotic spindle

The ability of cells to divide along their longest axis has been proposed to play an important role in maintaining epithelial tissue homeostasis in many systems. Because the division plane is largely set by the position of the anaphase spindle, it is important to understand how spindles become oriented. While several molecules have been identified that play key roles in spindle orientation across systems, most notably Mud/NuMA and cortical dynein, the precise mechanism by which spindles detect and align with the long cell axis remain poorly understood. Here, in exploring the dynamics of spindle orientation in mechanically distinct regions of the fly notum, we find that the ability of cells to properly reorient their divisions depends on local tissue tension. Thus, spindles reorient to align with the long cell axis in regions where isotropic tension is elevated, but fail to do so in elongated cells within the crowded midline, where tension is low, or in regions that have been mechanically isolated from the rest of the tissue via laser ablation. Importantly, these differences in spindle behavior outside and inside the midline can be recapitulated by corresponding changes in tension induced by perturbations that alter nonmuscle myosin II activity. These data lead us to propose that isotropic tension within an epithelium provides cells with a mechanically stable substrate upon which localized cortical motor complexes can act on astral microtubules to orient the spindle.

Cuticular reticulation replicates the pattern of epidermal cells in lowermost Cambrian scalidophoran worms

The cuticle of ecdysozoans (Panarthropoda, Scalidophora, Nematoida) is secreted by underlying epidermal cells and renewed via ecdysis. We explore here the relationship between epidermis and external cuticular ornament in stem-group scalidophorans from the early Cambrian of China (Kuanchuanpu Formation; ca 535 Ma) that had two types of microscopic polygonal cuticular networks with either straight or microfolded boundaries. Detailed comparisons with modern scalidophorans (priapulids) indicate that these networks faithfully replicate the cell boundaries of the epidermis. This suggests that the cuticle of early scalidophorans formed through the fusion between patches of extracellular material secreted by epidermal cells, as observed in various groups of present-day ecdysozoans, including arthropods. Key genetic, biochemical and mechanical processes associated with ecdysis and cuticle formation seem to have appeared very early (at least not later than 535 Ma) in the evolution of ecdysozoans. Microfolded reticulation is likely to be a mechanical response to absorbing contraction exerted by underlying muscles. The polygonal reticulation in early and extant ecdysozoans is clearly a by-product of the epidermal cell pavement and interacted with the sedimentary environment.

Cellular Heterogeneity in Pressure and Growth Emerges from Tissue Topology and Geometry

Cell-to-cell heterogeneity prevails in many systems, as exemplified by cell growth, although the origin and function of such heterogeneity are often unclear. In plants, growth is physically controlled by cell wall mechanics and cell hydrostatic pressure, alias turgor pressure. Whereas cell wall heterogeneity has received extensive attention, the spatial variation of turgor pressure is often overlooked. Here, combining atomic force microscopy and a physical model of pressurized cells, we show that turgor pressure is heterogeneous in the Arabidopsis shoot apical meristem, a population of stem cells that generates all plant aerial organs. In contrast with cell wall mechanical properties that appear to vary stochastically between neighboring cells, turgor pressure anticorrelates with cell size and cell neighbor number (local topology), in agreement with the prediction by our model of tissue expansion, which couples cell wall mechanics and tissue hydraulics. Additionally, our model predicts two types of correlations between pressure and cellular growth rate, where high pressure may lead to faster- or slower-than-average growth, depending on cell wall extensibility, yield threshold, osmotic pressure, and hydraulic conductivity. The meristem exhibits one of these two regimes, depending on conditions, suggesting that, in this tissue, water conductivity may contribute to growth control. Our results unravel cell pressure as a source of patterned heterogeneity and illustrate links between local topology, cell mechanical state, and cell growth, with potential roles in tissue homeostasis.

Growth control in the Drosophila wing disk

The regulation of size and shape is a fundamental requirement of biological development and has been a subject of scientific study for centuries, but we still lack an understanding of how organisms know when to stop growing. Imaginal wing disks of the fruit fly Drosophila melanogaster, which are precursors of the adult wings, are an archetypal tissue for studying growth control. The growth of the disks is dependent on many inter- and intra-organ factors such as morphogens, mechanical forces, nutrient levels, and hormones that influence gene expression and cell growth. Extracellular signals are transduced into gene-control signals via complex signal transduction networks, and since cells typically receive many different signals, a mechanism for integrating the signals is needed. Our understanding of the effect of morphogens on tissue-level growth regulation via individual pathways has increased significantly in the last half century, but our understanding of how multiple biochemical and mechanical signals are integrated to determine whether or not a cell decides to divide is still rudimentary. Numerous fundamental questions are involved in understanding the decision-making process, and here we review the major biochemical and mechanical pathways involved in disk development with a view toward providing a basis for beginning to understand how multiple signals can be integrated at the cell level, and how this translates into growth control at the level of the imaginal disk. This article is categorized under: Analytical and Computational Methods > Computational Methods Biological Mechanisms > Cell Signaling Models of Systems Properties and Processes > Cellular Models.

Bridging the gap between single-cell migration and collective dynamics

Motivated by the wealth of experimental data recently available, we present a cellular-automaton-based modeling framework focussing on high-level cell functions and their concerted effect on cellular migration patterns. Specifically, we formulate a coarse-grained description of cell polarity through self-regulated actin organization and its response to mechanical cues. Furthermore, we address the impact of cell adhesion on collective migration in cell cohorts. The model faithfully reproduces typical cell shapes and movements down to the level of single cells, yet allows for the efficient simulation of confluent tissues. In confined circular geometries, we find that specific properties of individual cells (polarizability; contractility) influence the emerging collective motion of small cell cohorts. Finally, we study the properties of expanding cellular monolayers (front morphology; stress and velocity distributions) at the level of extended tissues.

Biomimetic modification of dual porosity poly(2-hydroxyethyl methacrylate) hydrogel scaffolds-porosity and stem cell growth evaluation

The macroporous synthetic poly(2-hydroxyethyl methacrylate) (pHEMA) hydrogels as 3D cellular scaffolds with specific internal morphology, so called dual pore size, were designed and studied. The morphological microstructure of hydrogels was characterized in the gel swollen state and the susceptibility of gels for stem cells was evaluated. The effect of specific chemical groups covalently bound in the hydrogel network by copolymerization on cell adhesion and growth, followed by effect of laminin coating were investigated. The evaluated gels contained either carboxyl groups of the methacrylic acid or quaternary ammonium groups brought by polymerizable ammonium salt or their combinations. The morphology of swollen gel was visualized using the laser scanning confocal microscopy. All hydrogels had very similar porous structures - their matrices contained large pores (up to 10² μm) surrounded with gel walls with small pores (10⁰ μm). The total pore volume in hydrogels swollen in buffer solution ranged between 69 and 86 vol%. Prior to the seeding of the mouse embryonal stem cells, the gels were coated with laminin. The hydrogel with quaternary ammonium groups (with or without laminin) stimulated the cell growth the most. The laminin coating lead to a significant and quaternary ammonium groups. The gel chemical modification influenced also the topology of cell coverage that ranged from individual cell clusters to well dispersed multi cellular structures. Findings in this study point out the laser scanning confocal microscopy as an irreplaceable method for a precise and quick assessment of the hydrogel morphology. In addition, these findings help to optimize the chemical composition of the hydrogel scaffold through the combination of chemical and biological factors leading to intensive cell attachment and proliferation.

Topological traits of a cellular pattern versus growth rate anisotropy in radish roots

The topology of a cellular pattern, which means the spatial arrangement of cells, directly corresponds with cell packing, which is crucial for tissue and organ functioning. The topological features of cells that are typically analyzed are the number of their neighbors and the cell area. To date, the objects of most topological studies have been the growing cells of the surface tissues of plant and animal organs. Some of these researches also provide verification of Lewis's Law concerning the linear correlation between the number of neighboring cells and the cell area. Our aim was to analyze the cellular topology and applicability of Lewis's Law to an anisotropically growing plant organ. The object of our study was the root apex of radish. Based on the tensor description of plant organ growth, we specified the level of anisotropy in specific zones (the root proper, the columella of the cap and the lateral parts of the cap) and in specific types of both external (epidermis) and internal tissues (stele and ground tissue) of the apex. The strongest anisotropy occurred in the root proper, while both zones of the cap showed an intermediate level of anisotropy of growth. Some differences in the topology of the cellular pattern in the zones were also detected; in the root proper, six-sided cells predominated, while in the root cap columella and in the lateral parts of the cap, most cells had five neighbors. The correlation coefficient rL between the number of neighboring cells and the cell area was high in the apex as a whole as well as in all of the zones except the root proper and in all of the tissue types except the ground tissue. In general, Lewis's Law was fulfilled in the anisotropically growing radish root apex. However, the level of the applicability (rL value) of Lewis's Law was negatively correlated with the level of the anisotropy of growth, which may suggest that in plant organs in the regions of anisotropic growth, the number of neighboring cells is less dependent on the cell size.

Multicellular Systems Biology: Quantifying Cellular Patterning and Function in Plant Organs Using Network Science

Organ function is at least partially shaped and constrained by the organization of their constituent cells. Extensive investigation has revealed mechanisms explaining how these patterns are generated, with less being known about their functional relevance. In this paper, a methodology to discretize and quantitatively analyze cellular patterning is described. By performing global organ-scale cellular interaction mapping, the organization of cells can be extracted and analyzed using network science. This provides a means to take the developmental analysis of cellular organization in complex organisms beyond qualitative descriptions and provides data-driven approaches to inferring cellular function. The bridging of a structure-function relationship in hypocotyl epidermal cell patterning through global topological analysis provides support for this approach. The analysis of cellular topologies from patterning mutants further enables the contribution of gene activity toward the organizational properties of tissues to be linked, bridging molecular and tissue scales. This systems-based approach to investigate multicellular complexity paves the way to uncovering the principles of complex organ design and achieving predictive genotype-phenotype mapping.

E-cadherin expression pattern during zebrafish embryonic epidermis development

Background: E-cadherin is the major adhesion receptor in epithelial adherens junctions (AJs). On established epidermis, E-cadherin performs fine-tuned cell-cell contact remodeling to maintain tissue integrity, which is characterized by modulation of cell shape, size and packing density. In zebrafish, the organization and distribution of E-cadherin in AJs during embryonic epidermis development remain scarcely described.

Methods: Combining classical immunofluorescence, deconvolution microscopy and 3D-segmentation of AJs in epithelial cells, a quantitative approach was implemented to assess the spatial and temporal distribution of E-cadherin across zebrafish epidermis between 24 and 72 hpf.

Results: increasing levels of E-cadh protein parallel higher cell density and the appearance of hexagonal cells in the enveloping layer (EVL) as well as the establishments of new cell-cell contacts in the epidermal basal layer (EBL), being significantly between 31 and 48 hpf .

Conclusions: Increasing levels of E-cadherin in AJs correlates with extensive changes in cell morphology towards hexagonal packing during the epidermis morphogenesis.

Emergence of robust patterns from local rules during plant development

The formation of spatial and temporal patterns is an essential component of organismal development. Patterns can be observed on every level from subcellular to organismal and may emerge from local rules that correspond to the interactions between molecules, cells, or tissues. The emergence of robust patterns may seem in contradiction with the prominent heterogeneity at subcellular and cellular scales, however it has become increasingly clear that heterogeneity can be instrumental for pattern formation. Here we review recent examples in plant development, involving genetic regulation, cell arrangement, growth and signal gradient. We discuss how patterns emerge from local rules, whether heterogeneity is stochastic or can be patterned, and whether stochastic noise is amplified or requires filtering for robust patterns to be achieved. We also stress the importance of modelling in investigating such questions.

Global Topological Order Emerges through Local Mechanical Control of Cell Divisions in the Arabidopsis Shoot Apical Meristem

The control of cell position and division act in concert to dictate multicellular organization in tissues and organs. How these processes shape global order and molecular movement across organs is an outstanding problem in biology. Using live 3D imaging and computational analyses, we extracted networks capturing cellular connectivity dynamics across the Arabidopsis shoot apical meristem (SAM) and topologically analyzed the local and global properties of cellular architecture. Locally generated cell division rules lead to the emergence of global tissue-scale organization of the SAM, facilitating robust global communication. Cells that lie upon more shorter paths have an increased propensity to divide, with division plane placement acting to limit the number of shortest paths their daughter cells lie upon. Cell shape heterogeneity and global cellular organization requires KATANIN, providing a multiscale link between cell geometry, mechanical cell-cell interactions, and global tissue order.

Development of epithelial tissues: How are cleavage planes chosen?

The cross-section of a cell in a monolayer epithelial tissue can be modeled mathematically as a k-sided polygon. Empirically studied distributions of the proportions of k-sided cells in epithelia show remarkable similarities in a wide range of evolutionarily distant organisms. A variety of mathematical models have been proposed for explaining this phenomenon. The highly parsimonious simulation model of (Patel et al., PLoS Comput. Biol., 2009) that takes into account only the number of sides of a given cell and cell division already achieves a remarkably good fit with empirical distributions from Drosophila, Hydra, Xenopus, Cucumber, and Anagallis. Within the same modeling framework as in that paper, we introduce additional options for the choice of the endpoints of the cleavage plane that appear to be biologically more realistic. By taking the same data sets as our benchmarks, we found that combinations of some of our new options consistently gave better fits with each of these data sets than previously studied ones. Both our algorithm and simulation data are made available as research tools for future investigations.

Multifaceted regulation and functions of YAP/TAZ in tumors (Review)

The Hippo pathway, initially identified through screenings for mutant tumor suppressors in Drosophila, is an evolutionarily conserved signaling pathway that controls organ size by regulating cell proliferation and apoptosis. Abnormal regulation of the Hippo pathway may lead to cancer in mammals. As the major downstream effectors of the Hippo pathway, unphosphorylated Yes-associated protein (YAP) and its homolog transcriptional co-activator TAZ (also called WWTR1) (hereafter called YAP/TAZ) are translocated into the nucleus. In the nucleus, in order to induce target gene expression, YAP/TAZ bind to the TEA domain (TEAD) proteins, and this binding subsequently promotes cell proliferation and inhibits apoptosis. In contrast, as key regulators of tumorigenesis and development, YAP/TAZ are phosphorylated and regulated by multiple molecules and pathways including Lats1/2 of Hippo, Wnt and G-protein-coupled receptor (GPCR) signaling, with a regulatory role in cell physiology, tumor cell development and pathological abnormalities simultaneously. In particular, the crucial role of YAP/TAZ in tumors ensures their potential as targets in designing anticancer drugs. To date, mounting research has elucidated the suppression of YAP/TAZ via effective inhibitors, which significantly highlights their application in cancer treatment. In the present review, we focus on the functions of YAP/TAZ in cancer, discuss their potential as new therapeutic target for tumor treatment, and provide valuable suggestions for further study in this field.

Mechanical Forces Program the Orientation of Cell Division during Airway Tube Morphogenesis

Oriented cell division plays a key role in controlling organogenesis. The mechanisms for regulating division orientation at the whole-organ level are only starting to become understood. By combining 3D time-lapse imaging, mouse genetics, and mathematical modeling, we find that global orientation of cell division is the result of a combination of two types of spindles with distinct spindle dynamic behaviors in the developing airway epithelium. Fixed spindles follow the classic long-axis rule and establish their division orientation before metaphase. In contrast, rotating spindles do not strictly follow the long-axis rule and determine their division orientation during metaphase. By using both a cell-based mechanical model and stretching-lung-explant experiments, we showed that mechanical force can function as a regulatory signal in maintaining the stable ratio between fixed spindles and rotating spindles. Our findings demonstrate that mechanical forces, cell geometry, and oriented cell division function together in a highly coordinated manner to ensure normal airway tube morphogenesis.

Cytoplasmic Cl- couples membrane remodeling to epithelial morphogenesis

Chloride is the major free anion in the extracellular space (>100 mM) and within the cytoplasm in eukaryotes (10 ∼ 20 mM). Cytoplasmic Cl- level is dynamically regulated by Cl- channels and transporters. It is well established that movement of Cl- across the cell membrane is coupled with cell excitability through changes in membrane potential and with water secretion. However, whether cytoplasmic Cl- plays additional roles in animal development and tissue homeostasis is unknown. Here we use genetics, cell biological and pharmacological tools to demonstrate that TMEM16A, an evolutionarily conserved calcium-activated chloride channel (CaCC), regulates cytoplasmic Cl- homeostasis and promotes plasma membrane remodeling required for mammalian epithelial morphogenesis. We demonstrate that TMEM16A-mediated control of cytoplasmic Cl- regulates the organization of the major phosphoinositide species PtdIns(4,5)P2 into microdomains on the plasma membrane, analogous to processes that cluster soluble and membrane proteins into phase-separated droplets. We further show that an adequate cytoplasmic Cl- level is required for proper endocytic trafficking and membrane supply during early stages of ciliogenesis and adherens junction remodeling. Our study thus uncovers a critical function of CaCC-mediated cytoplasmic Cl- homeostasis in controlling the organization of PtdIns(4,5)P2 microdomains and membrane remodeling. This newly defined role of cytoplasmic Cl- may shed light on the mechanisms of intracellular Cl- signaling events crucial for regulating tissue architecture and organelle biogenesis during animal development.

Modulation of collective cell behaviour by geometrical constraints

Intracellular and extracellular mechanical forces play a crucial role during tissue growth, modulating nuclear shape and function and resulting in complex collective cell behaviour. However, the mechanistic understanding of how the orientation, shape, symmetry and homogeneity of cells are affected by environmental geometry is still lacking. Here we investigate cooperative cell behaviour and patterns under geometric constraints created by topographically patterned substrates. We show how cells cooperatively adopt their geometry, shape, positioning of the nucleus and subsequent proliferation activity. Our findings indicate that geometric constraints induce significant squeezing of cells and nuclei, cytoskeleton reorganization, drastic condensation of chromatin resulting in a change in the cell proliferation rate and the anisotropic growth of cultures. Altogether, this work not only demonstrates complex non-trivial collective cellular responses to geometrical constraints but also provides a tentative explanation of the observed cell culture patterns grown on different topographically patterned substrates. These findings provide important fundamental knowledge, which could serve as a basis for better controlled tissue growth and cell-engineering applications.

Pavement cells and the topology puzzle

D'Arcy Thompson emphasised the importance of surface tension as a potential driving force in establishing cell shape and topology within tissues. Leaf epidermal pavement cells grow into jigsaw-piece shapes, highly deviating from such classical forms. We investigate the topology of developing Arabidopsis leaves composed solely of pavement cells. Image analysis of around 50,000 cells reveals a clear and unique topological signature, deviating from previously studied epidermal tissues. This topological distribution is established early during leaf development, already before the typical pavement cell shapes emerge, with topological homeostasis maintained throughout growth and unaltered between division and maturation zones. Simulating graph models, we identify a heuristic cellular division rule that reproduces the observed topology. Our parsimonious model predicts how and when cells effectively place their division plane with respect to their neighbours. We verify the predicted dynamics through in vivo tracking of 800 mitotic events, and conclude that the distinct topology is not a direct consequence of the jigsaw piece-like shape of the cells, but rather owes itself to a strongly life history-driven process, with limited impact from cell-surface mechanics.

Structural Characterization and Statistical-Mechanical Model of Epidermal Patterns

In proliferating epithelia of mammalian skin, cells of irregular polygon-like shapes pack into complex, nearly flat two-dimensional structures that are pliable to deformations. In this work, we employ various sensitive correlation functions to quantitatively characterize structural features of evolving packings of epithelial cells across length scales in mouse skin. We find that the pair statistics in direct space (correlation function) and Fourier space (structure factor) of the cell centroids in the early stages of embryonic development show structural directional dependence (statistical anisotropy), which is a reflection of the fact that cells are stretched, which promotes uniaxial growth along the epithelial plane. In the late stages, the patterns tend toward statistically isotropic states, as cells attain global polarization and epidermal growth shifts to produce the skin's outer stratified layers. We construct a minimalist four-component statistical-mechanical model involving effective isotropic pair interactions consisting of hard-core repulsion and extra short-range soft-core repulsion beyond the hard core, whose length scale is roughly the same as the hard core. The model parameters are optimized to match the sample pair statistics in both direct and Fourier spaces. By doing this, the parameters are biologically constrained. In contrast with many vertex-based models, our statistical-mechanical model does not explicitly incorporate information about the cell shapes and interfacial energy between cells; nonetheless, our model predicts essentially the same polygonal shape distribution and size disparity of cells found in experiments, as measured by Voronoi statistics. Moreover, our simulated equilibrium liquid-like configurations are able to match other nontrivial unconstrained statistics, which is a testament to the power and novelty of the model. The array of structural descriptors that we deploy enable us to distinguish between normal, mechanically deformed, and pathological skin tissues. Our statistical-mechanical model enables one to generate tissue microstructure at will for further analysis. We also discuss ways in which our model might be extended to better understand morphogenesis (in particular the emergence of planar cell polarity), wound healing, and disease-progression processes in skin, and how it could be applied to the design of synthetic tissues.

A Method to Categorize 2-Dimensional Patterns Using Statistics of Spatial Organization

We developed a measurement framework of spatial organization to categorize 2-dimensional patterns from 2 multiscalar biological architectures. We propose that underlying shapes of biological entities can be approached using the statistical concept of degrees of freedom, defining it through expansion of area variability in a pattern. To help scope this suggestion, we developed a mathematical argument recognizing the deep foundations of area variability in a polygonal pattern (spatial heterogeneity). This measure uses a parameter called eutacticity. Our measuring platform of spatial heterogeneity can assign particular ranges of distribution of spatial areas for 2 biological architectures: ecological patterns of Namibia fairy circles and epithelial sheets. The spatial organizations of our 2 analyzed biological architectures are demarcated by being in a particular position among spatial order and disorder. We suggest that this theoretical platform can give us some insights about the nature of shapes in biological systems to understand organizational constraints.

Rules of tissue packing involving different cell types: human muscle organization

Natural packed tissues are assembled as tessellations of polygonal cells. These include skeletal muscles and epithelial sheets. Skeletal muscles appear as a mosaic composed of two different types of cells: the "slow" and "fast" fibres. Their relative distribution is important for the muscle function but little is known about how the fibre arrangement is established and maintained. In this work we capture the organizational pattern in two different healthy muscles: biceps brachii and quadriceps. Here we show that the biceps brachii muscle presents a particular arrangement, based on the different sizes of slow and fast fibres. By contrast, in the quadriceps muscle an unbiased distribution exists. Our results indicate that the relative size of each cellular type imposes an intrinsic organization into natural tessellations. These findings establish a new framework for the analysis of any packed tissue where two or more cell types exist.

Rac1/RhoA antagonism defines cell-to-cell heterogeneity during epidermal morphogenesis in nematodes

The antagonism between the GTPases Rac1 and RhoA controls cell-to-cell heterogeneity in isogenic populations of cells in vitro and epithelial morphogenesis in vivo. Its involvement in the regulation of cell-to-cell heterogeneity during epidermal morphogenesis has, however, never been addressed. We used a quantitative cell imaging approach to characterize epidermal morphogenesis at a single-cell level during early elongation of Caenorhabditis elegans embryos. This study reveals that a Rac1-like pathway, involving the Rac/Cdc42 guanine-exchange factor β-PIX/PIX-1 and effector PAK1/PAK-1, and a RhoA-like pathway, involving ROCK/LET-502, control the remodeling of apical junctions and the formation of basolateral protrusions in distinct subsets of hypodermal cells. In these contexts, protrusions adopt lamellipodia or an amoeboid morphology. We propose that lamella formation may reduce tension building at cell-cell junctions during morphogenesis. Cell-autonomous antagonism between these pathways enables cells to switch between Rac1- and RhoA-like morphogenetic programs. This study identifies the first case of cell-to-cell heterogeneity controlled by Rac1/RhoA antagonism during epidermal morphogenesis.

Fundamental physical cellular constraints drive self-organization of tissues

Morphogenesis is driven by small cell shape changes that modulate tissue organization. Apical surfaces of proliferating epithelial sheets have been particularly well studied. Currently, it is accepted that a stereotyped distribution of cellular polygons is conserved in proliferating tissues among metazoans. In this work, we challenge these previous findings showing that diverse natural packed tissues have very different polygon distributions. We use Voronoi tessellations as a mathematical framework that predicts this diversity. We demonstrate that Voronoi tessellations and the very different tissues analysed share an overriding restriction: the frequency of polygon types correlates with the distribution of cell areas. By altering the balance of tensions and pressures within the packed tissues using disease, genetic or computer model perturbations, we show that as long as packed cells present a balance of forces within tissue, they will be under a physical constraint that limits its organization. Our discoveries establish a new framework to understand tissue architecture in development and disease.

Coupling between apical tension and basal adhesion allow epithelia to collectively sense and respond to substrate topography over long distances

Epithelial sheets fold into complex topographies that contribute to their function in vivo. Cells can sense and respond to substrate topography in their immediate vicinity by modulating their interfacial mechanics, but the extent to which these mechanical properties contribute to their ability to sense substrate topography across length scales larger than a single cell has not been explored in detail. To study the relationship between the interfacial mechanics of single cells and their collective behavior as tissues, we grew cell-sheets on substrates engraved with surface features spanning macroscopic length-scales. We found that many epithelial cell-types sense and respond to substrate topography, even when it is locally nearly planar. Cells clear or detach from regions of local negative curvature, but not from regions with positive or no curvature. We investigated this phenomenon using a finite element model where substrate topography is coupled to epithelial response through a balance of tissue contractility and adhesive forces. The model correctly predicts the focal sites of cell-clearing and epithelial detachment. Furthermore, the model predicts that local tissue response to substrate curvature is a function of the surrounding topography of the substrate across long distances. Analysis of cell-cell and cell-substrate contact angles suggests a relationship between these single-cell interfacial properties, epithelial interfacial properties, and collective epithelial response to substrate topography. Finally, we show that contact angles change upon activation of oncogenes or inhibition of cell-contractility, and that these changes correlate with collective epithelial response. Our results demonstrate that in mechanically integrated epithelial sheets, cell contractility can be transmitted through multiple cells and focused by substrate topography to affect a behavioral response at distant sites.

The non-proliferative nature of ascidian folliculogenesis as a model of highly ordered cellular topology distinct from proliferative epithelia

Previous studies have addressed why and how mono-stratified epithelia adopt a polygonal topology. One major additional, and yet unanswered question is how the frequency of different cell shapes is achieved and whether the same distribution applies between non-proliferative and proliferative epithelia. We compared different proliferative and non-proliferative epithelia from a range of organisms as well as Drosophila melanogaster mutants, deficient for apoptosis or hyperproliferative. We show that the distribution of cell shapes in non-proliferative epithelia (follicular cells of five species of tunicates) is distinctly, and more stringently organized than proliferative ones (cultured epithelial cells and Drosophila melanogaster imaginal discs). The discrepancy is not supported by geometrical constraints (spherical versus flat monolayers), number of cells, or apoptosis events. We have developed a theoretical model of epithelial morphogenesis, based on the physics of divided media, that takes into account biological parameters such as cell-cell contact adhesions and tensions, cell and tissue growth, and which reproduces the effects of proliferation by increasing the topological heterogeneity observed experimentally. We therefore present a model for the morphogenesis of epithelia where, in a proliferative context, an extended distribution of cell shapes (range of 4 to 10 neighbors per cell) contrasts with the narrower range of 5-7 neighbors per cell that characterizes non proliferative epithelia.

Cell size versus body size in geophilomorph centipedes

Variation in animal body size is the result of a complex interplay between variation in cell number and cell size, but the latter has seldom been considered in wide-ranging comparative studies, although distinct patterns of variation have been described in the evolution of different lineages. We investigated the correlation between epidermal cell size and body size in a sample of 29 geophilomorph centipede species, representative of a wide range of body sizes, from 6 mm dwarf species to gigantic species more than 200 mm long, exploiting the marks of epidermal cells on the overlying cuticle in the form of micro-sculptures called scutes. We found conspicuous and significant variation in average scute area, both between suprageneric taxa and between genera, while the within-species range of variation is comparatively small. This supports the view that the average epidermal cell size is to some extent taxon specific. However, regression analyses show that neither body size nor the number of leg-bearing segments explain this variation, which suggests that cell size is not an usual target of change for body size evolution in this group of arthropods, although there is evidence of its correlation with other morphological variables, like cuticle thickness. Scute sizes of miniaturized geophilomorph species are well within the range of the lineage to which the species belong, suggesting recent evolutionary transitions to smaller body size.

Heterogeneity for IGF-II production maintained by public goods dynamics in neuroendocrine pancreatic cancer

The extensive intratumor heterogeneity revealed by sequencing cancer genomes is an essential determinant of tumor progression, diagnosis, and treatment. What maintains heterogeneity remains an open question because competition within a tumor leads to a strong selection for the fittest subclone. Cancer cells also cooperate by sharing molecules with paracrine effects, such as growth factors, and heterogeneity can be maintained if subclones depend on each other for survival. Without strict interdependence between subclones, however, nonproducer cells can free-ride on the growth factors produced by neighboring producer cells, a collective action problem known in game theory as the "tragedy of the commons," which has been observed in microbial cell populations. Here, we report that similar dynamics occur in cancer cell populations. Neuroendocrine pancreatic cancer (insulinoma) cells that do not produce insulin-like growth factor II (IGF-II) grow slowly in pure cultures but have a proliferation advantage in mixed cultures, where they can use the IGF-II provided by producer cells. We show that, as predicted by evolutionary game theory, producer cells do not go extinct because IGF-II acts as a nonlinear public good, creating negative frequency-dependent selection that leads to a stable coexistence of the two cell types. Intratumor cell heterogeneity can therefore be maintained even without strict interdependence between cell subclones. Reducing the amount of growth factors available within a tumor may lead to a reduction in growth followed by a new equilibrium, which may explain relapse in therapies that target growth factors.