Non-Turing stripes and spots: a novel mechanism for biological cell clustering
Introduction
Mesenchymal condensation, in which mesenchymal cells (i.e., cells that form part of the bulk tissues rather than the epithelial surfaces) condense (i.e., coalesce) into compact clusters) is the earliest stage of organogenesis and is crucial to the development of skeletal and other mesenchymal tissues (e.g. cartilage, kidney, lung, etc.). The morphology of the condensations lays down the shape of future organs [1], [2]. In this paper we propose a novel developmental mechanism for mesenchymal condensation, based on cell–cell and cell–extracellular matrix (ECM, the loosely structured noncellular scaffolding which surrounds and supports the cells) adhesion. We use it to explain the clustering and differentiation of a previously dispersed population of mesenchymal cells during precartilage mesenchymal condensation in chick limb chondrogenesis (bone patterning and formation), a specific case of the more general developmental phenomenon. We compare the results to those the Turing mechanism produces, which Newman and others have proposed as explaining mesenchymal condensation [3].
Turing proposed a reaction-diffusion mechanism for biological pattern formation in 1952 [4] and patterns similar to Turing-type also occur in certain chemical reactions [5]. A Turing pattern results spontaneously from the reaction and diffusion of two initially homogeneous chemical substances: activator (A) and inhibitor (I), provided that I diffuses faster than A. A promotes the production of both A and I, while I inhibits the production of A. A zone of elevated concentration of A increases the concentration of I over a larger area, inhibiting A in the surrounding region and leading to an inhomogeneous spatial distribution of A and I. Basic patterns include spots, simple stripes and labyrinths, while fine tuning the parameters can produce a whole zoo of patterns.
Many explanations of biological patterning invoke the Turing mechanism, e.g. for the patterning of animal coats, feather buds and fish skin [6], [7]. In this class of models, the standard Turing mechanism establishes a chemical pre-pattern and cells interpret the pattern by differentiating only where the chemical (morphogen) concentration is above a certain threshold [1], [2]. The Turning model is incomplete in several respects: the huge inventory of molecules involved in morphogenesis includes very few which behave like classical activators or inhibitors. No activators or inhibitors are known for mesenchymal condensation. The cells do not influence the static chemical ‘prepattern’, making development into a slave process, while real cells actively reshape their environment. The model requires an additional process (e.g. selective differentiation or chemotaxis) to lock in the labile prepattern. These problems led us to look for a simple patterning mechanism that could reproduce the basic Turing patterns using only biological mechanisms proved experimentally to act during mesenchymal condensation.
Section snippets
Biology
Recent experiments show that extracellular matrix molecules and membrane-bound cell adhesion molecules, such as fibronectin (a sticky protein in the ECM) and N-CAM (a cell–cell adhesion molecule) play an important role in initiating cell clustering during precartilage mesenchymal condensation [8].
Precartilage mesenchymal cells are swarming cells that secrete a variety of complex molecules which form the ECM, the porous environment which is largely responsible for the mechanical integrity of
Methods
We implement a Cellular Potts Model (CPM) [21], [22] simulation on a two-dimensional lattice. The two-dimensional model reflects the quasi-two-dimensional micromass experiment. Our simulation produces a density dependent morphology which we compare to experimental data. We also discuss possible extensions to our model to study the formation of pattern and shape during limb chondrogenesis.
Glazier and Graner introduced the CPM to describe surface-energy-driven biological processes at the cellular
Model and simulation
We propose the following model for precartilage mesenchymal condensation in limb chondrogenesis: mesenchymal cells produce nondiffusing fibronectin and deposit it onto the substrate. Cells execute a random walk biased by their binding more strongly to fibronectin than to the substrate (haptotaxis). Fibronectin–cell surface binding upregulates the production of cell–cell adhesion molecules such as N-cadherin, increasing cells’ adhesiveness to each other. The pattern continues to evolve until the
Results
Precartilage mesenchymal condensations at medium and low culture density (Fig. 2) show very different morphologies. When the density is relatively high (Fig. 2(A)), the condensations consist of both spots and ‘finger-like’ stripes, while dots form during low density condensation (Fig. 2(B)). In Fig. 2(C) we have a cell concentration of 35%, while in Fig. 2(D) the concentration is about 15%. Both cases show reasonable agreement with experiment. Figs. 2(D and E) show the fibronectin distribution
Discussion
Our simulations show that spot and stripe patterns can result from a density dependent mechanism unrelated to the Turing mechanism. In our model, the instability arises from positive feedback in the fibronectin production. Cells tend to stay longer in regions with fibronectin, and hence become more adhesive and produce more fibronectin in those regions, further increasing the time the cells spend there. The wavelength grows continuously until cells stick and clusters cannot coalesce any more.
Acknowledgements
Grants NSF-IBN-008365, NASA-NAG2-1619 and CAPES-Brazil. The authors are grateful to Stuart Newman for his experimental instruction and discussions about the biological aspects of the work, though his interpretation of the experiments differs from ours. We thank Mark Alber, Cornelis Weijer and Gemunu Gunaratne for critical reading and comments.
References (34)
- et al.
Role of transforming growth factor-beta in chondrogenic pattern formation in the embryonic limbstimulation of mesenchymal condensation and fibronectin gene expression by exogenous TGF-beta and evidence for endogenous TGF-beta-like activity
Dev. Biol.
(1991) - et al.
Ectodermal FGFs induce perinodular inhibition of limb chondrogenesis in vitro and in vivo via FGF receptor 2
Dev. Biol.
(2002) - et al.
Cell cycle length affects gene expression and pattern formation in limbs
Dev. Biol.
(1997) - et al.
Cell migration and chick limb developmentchemotactic action of FGF-4 and the AER
Dev. Biol.
(1999) Principles of Development
(2001)Developmental Biology
(2003)- et al.
Dynamics of skeletal pattern formation in developing chick limb
Science
(1979) The chemical basis of morphogenesis
Phil. Trans. Roy. Soc. B
(1952)- et al.
Experimental evidence of a sustained standing turing-type nonequilibrium chemical pattern
Phys. Rev. Lett.
(1990) Mathematical Biology
(1989)
A viable reaction-diffusion wave on the skin of Pomacanthus, a marine angelfish
Nature
All for one and one for allcondensations and the initiation of skeletal development
BioEssays
Defects in mesoderm, neural tube and vascular development in mouse embryos lacking fibronectin
Development
Fibronectin gene expression during limb cartilage differentiation
Development
TGFbeta2 acts as an “activator” molecule in reaction–diffusion model and is involved in cell sorting phenomenon in mouse limb micromass culture
Dev. Dyn.
Expression and functional involvement of N-cadherin in embryonic limb chondrogenesis
Development
Cited by (31)
The design principles of discrete turing patterning systems
2021, Journal of Theoretical BiologyCitation Excerpt :Cellular automata constitute one further level of abstraction in that they model concentrations by a finite set of discrete states. Variations of CAs have been used extensively in the context of biological pattern formation (Dormann et al., 2001; Deutsch and Dormann, 2005; Alber et al., 2003; Zeng et al., 2004; Zeng et al., 2002; Boon et al., 1996; Weimar and Boon, 1994). These studies have greatly enhanced our understanding of Turing patterns in different biological scenarios.
Mathematical modeling of chondrogenic pattern formation during limb development: Recent advances in continuous models
2020, Mathematical BiosciencesCitation Excerpt :The selection of models is by no means complete, and we mostly concentrate on models proposed within the last 15 years or so. Many other influential models are not treated in depth here, e.g. earlier reaction-diffusion models such as to Newman and Frisch’s model [109] and the work by Hentschel et al. [69], the mechanochemical models of Murray and Oster [124,125] and Zeng et al. [181] or the boundary model by Meinhardt [95,97]. Among these models, the most molecularly detailed model is that by Hentschel et al. [69], which incorporates a system of activator-inhibitor reaction-diffusion equations together with an explicitly modeled cell density.
Coordinated Action of N-CAM, N-cadherin, EphA4, and ephrinB2 Translates Genetic Prepatterns into Structure during Somitogenesis in Chick
2008, Current Topics in Developmental BiologyCitation Excerpt :To simulate our somite-segmentation model, we implement it using the Glazier–Graner–Hogeweg model (GGH) (also known as the Extended or Cellular Potts Model) (Glazier and Graner, 1992; 1993). Extensive comparisons between experiments and GGH simulations have validated GGH methods (Mareé and Hogeweg, 2001; Zeng et al., 2004; Chaturvedi et al., 2005; Poplawski et al., 2007; Merks et al., 2006; Merks and Glazier, 2006) for multicell morphogenesis modeling. To translate our model into a GGH simulation, we must define the GGH parameters and their temporal changes for each cell.
Adhesion between cells, diffusion of growth factors, and elasticity of the AER produce the paddle shape of the chick limb
2007, Physica A: Statistical Mechanics and its ApplicationsA cell-centered approach to developmental biology
2005, Physica A: Statistical Mechanics and its ApplicationsModeling the extracellular matrix in cell migration and morphogenesis: a guide for the curious biologist
2024, Frontiers in Cell and Developmental Biology
- 1
Permanent address: Instituto de Fı́sica, Universidade Federal do Rio Grande do Sul, Porto Alegre, RS 91501-970, Brazil.