The role of spatially controlled cell proliferation in limb bud morphogenesis

Downstream branches of receptor tyrosine kinase signaling act interdependently to shape the face

Background –

Previously we found that increasing fibroblast growth factor (FGF) signaling in the neural crest cells within the frontonasal process (FNP) of the chicken embryo caused dysmorphology that was correlated with reduced proliferation, disrupted cellular orientation, and lower MAPK activation but no change in PLCy and PI3K activation. This suggests RTK signaling may drive craniofacial morphogenesis through specific downstream effectors that affect cellular activities. In this study we inhibited three downstream branches of RTK signaling to determine their role in regulating cellular activities and how these changes affect morphogenesis of the FNP.

Results –

Small molecule inhibitors of MEK1/2, PI3K, and PLCy were delivered individually and in tandem to the right FNP of chicken embryos. All treatments caused asymmetric proximodistal truncation on the treated side and a mild expansion on the untreated side compared to DMSO control treated FNPs. Inhibiting each pathway caused similar decreased proliferation and disrupted cellular orientation, but did not affect apoptosis.

Conclusions –

Since RTK signaling is a ubiquitous and tightly regulated biochemical system we conclude that the downstream pathways are robust to developmental perturbation through redundant signaling systems.

Proliferation symmetry breaking in growing tissues

Morphogenesis of developing tissues results from anisotropic growth, typically driven by polarized patterns of gene expression. Here we propose an alternative model of anisotropic growth driven by self-organized feedback between cell polarity, mechanical pressure, and cell division rates. Specifically, cell polarity alignment can induce spontaneous symmetry breaking in proliferation, resulting from the anisotropic distribution of mechanical pressure in the tissue. We show that proliferation anisotropy can be controlled by cellular elasticity, motility and contact inhibition, thereby elucidating the design principles for anisotropic morphogenesis.

Adult caudal fin shape is imprinted in the embryonic fin fold

Appendage shape is formed during development (and re-formed during regeneration) according to spatial and temporal cues that orchestrate local cellular morphogenesis. The caudal fin is the primary appendage used for propulsion in most fish species, and exhibits a range of distinct morphologies adapted for different swimming strategies, however the molecular mechanisms responsible for generating these diverse shapes remain mostly unknown. In zebrafish, caudal fins display a forked shape, with longer supportive bony rays at the periphery and shortest rays at the center. Here, we show that a premature, transient pulse of sonic hedgehog a (shha) overexpression during late embryonic development results in excess proliferation and growth of the central rays, causing the adult caudal fin to grow into a triangular, truncate shape. Both global and regional ectopic shha overexpression are sufficient to alter fin shape, and forked shape may be rescued by subsequent treatment with an antagonist of the canonical Shh pathway. The induced truncate fins show a decreased fin ray number and fail to form the hypural diastema that normally separates the dorsal and ventral fin lobes. While forked fins regenerate their original forked morphology, truncate fins regenerate truncate, suggesting that positional memory of the fin rays can be permanently altered by a transient treatment during embryogenesis. Ray finned fish have evolved a wide spectrum of caudal fin morphologies, ranging from truncate to forked, and the current work offers insights into the developmental mechanisms that may underlie this shape diversity.

A fibronectin gradient remodels mixed-phase mesoderm

Physical processes ultimately shape tissue during development. Two emerging proposals are that cells migrate toward stiffer tissue (durotaxis) and that the extent of cell rearrangements reflects tissue phase, but it is unclear whether and how these concepts are related. Here, we identify fibronectin-dependent tissue stiffness as a control variable that underlies and unifies these phenomena in vivo. In murine limb bud mesoderm, cells are either caged, move directionally, or intercalate as a function of their location along a stiffness gradient. A modified Landau phase equation that incorporates tissue stiffness accurately predicts cell diffusivity upon loss or gain of fibronectin. Fibronectin is regulated by WNT5A-YAP feedback that controls cell movements, tissue shape, and skeletal pattern. The results identify a key determinant of phase transition and show how fibronectin-dependent directional cell movement emerges in a mixed-phase environment in vivo.

Mesenchymal Wnts are required for morphogenetic movements of calvarial osteoblasts during apical expansion

Apical expansion of calvarial osteoblast progenitors from the cranial mesenchyme (CM) above the eye is integral to calvarial growth and enclosure of the brain. The cellular behaviors and signals underlying the morphogenetic process of calvarial expansion are unknown. Time-lapse light-sheet imaging of mouse embryos revealed calvarial progenitors intercalate in 3D in the CM above the eye, and exhibit protrusive and crawling activity more apically. CM cells express non-canonical Wnt/planar cell polarity (PCP) core components and calvarial osteoblasts are bidirectionally polarized. We found non-canonical ligand Wnt5a-/- mutants have less dynamic cell rearrangements and protrusive activity. Loss of CM-restricted Wntless (CM-Wls), a gene required for secretion of all Wnt ligands, led to diminished apical expansion of Osx+ calvarial osteoblasts in the frontal bone primordia in a non-cell autonomous manner without perturbing proliferation or survival. Calvarial osteoblast polarization, progressive cell elongation and enrichment for actin along the baso-apical axis were dependent on CM-Wnts. Thus, CM-Wnts regulate cellular behaviors during calvarial morphogenesis for efficient apical expansion of calvarial osteoblasts. These findings also offer potential insights into the etiologies of calvarial dysplasias.

Current research on mechanisms of limb bud development, and challenges for the next decade

The developmental mechanisms of limb buds have been studied in developmental biology as an excellent model of pattern formation. Chick embryos have contributed to the discovery of new principles in developmental biology, as it is easy to observe live embryos and manipulate embryonic tissues. Herein, I outline recent findings and future issues over the next decade regarding three themes, based on my research: limb positioning, proximal-distal limb elongation and digit identity determination. First, how hindlimb position is determined at the molecular level is described, with a focus on the transforming growth factor-β signaling molecule GDF11. Second, I explain how the cell population in the limb bud deforms with developmental progress, shaping the limb bud with elongation along the proximal-distal axis. Finally, I describe the developmental mechanisms that determine digit identity through the interdigits.

Compressive stress gradients direct mechanoregulation of anisotropic growth in the zebrafish jaw joint

Mechanical stimuli arising from fetal movements are critical factors underlying joint growth. Abnormal fetal movements negatively affect joint shape features with important implications for joint health, but the mechanisms by which mechanical forces from fetal movements influence joint growth are still unclear. In this research, we quantify zebrafish jaw joint growth in 3D in free-to-move and immobilised fish larvae between four and five days post fertilisation. We found that the main changes in size and shape in normally moving fish were in the ventrodorsal axis, while growth anisotropy was lost in the immobilised larvae. We next sought to determine the cell level activities underlying mechanoregulated growth anisotropy by tracking individual cells in the presence or absence of jaw movements, finding that the most dramatic changes in growth rates due to jaw immobility were in the ventrodorsal axis. Finally, we implemented mechanobiological simulations of joint growth with which we tested hypotheses relating specific mechanical stimuli to mechanoregulated growth anisotropy. Different types of mechanical stimulation were incorporated into the simulation to provide the mechanoregulated component of growth, in addition to the baseline (non-mechanoregulated) growth which occurs in the immobilised animals. We found that when average tissue stress over the opening and closing cycle of the joint was used as the stimulus for mechanoregulated growth, joint morphogenesis was not accurately predicted. Predictions were improved when using the stress gradients along the rudiment axes (i.e., the variation in magnitude of compression to magnitude of tension between local regions). However, the most accurate predictions were obtained when using the compressive stress gradients (i.e., the variation in compressive stress magnitude) along the rudiment axes. We conclude therefore that the dominant biophysical stimulus contributing to growth anisotropy during early joint development is the gradient of compressive stress experienced along the growth axes under cyclical loading.

An archetype and scaling of developmental tissue dynamics across species

Morphometric studies have revealed the existence of simple geometric relationships among various animal shapes. However, we have little knowledge of the mathematical principles behind the morphogenetic dynamics that form the organ/body shapes of different species. Here, we address this issue by focusing on limb morphogenesis in Gallus gallus domesticus (chicken) and Xenopus laevis (African clawed frog). To compare the deformation dynamics between tissues with different sizes/shapes as well as their developmental rates, we introduce a species-specific rescaled spatial coordinate and a common clock necessary for cross-species synchronization of developmental times. We find that tissue dynamics are well conserved across species under this spacetime coordinate system, at least from the early stages of development through the phase when basic digit patterning is established. For this developmental period, we also reveal that the tissue dynamics of both species are mapped with each other through a time-variant linear transformation in real physical space, from which hypotheses on a species-independent archetype of tissue dynamics and morphogenetic scaling are proposed.

Mesenchymal Wnts are required for morphogenetic movements of calvarial osteoblasts during apical expansion

Apical expansion of calvarial osteoblast progenitors from the cranial mesenchyme (CM) above the eye is integral for calvarial growth and enclosure of the brain. The cellular behaviors and signals underlying the morphogenetic process of calvarial expansion are unknown. During apical expansion, we found that mouse calvarial primordia have consistent cellular proliferation, density, and survival with complex tissue scale deformations, raising the possibility that morphogenetic movements underlie expansion. Time lapse light sheet imaging of mouse embryos revealed that calvarial progenitors intercalate in 3D to converge supraorbital arch mesenchyme mediolaterally and extend it apically. In contrast, progenitors located further apically exhibited protrusive and crawling activity. CM cells express non-canonical Wnt/Planar Cell Polarity (PCP) core components and calvarial osteoblasts are bidirectionally polarized. We found non-canonical ligand, Wnt5a-/- mutants have less dynamic cell rearrangements, protrusive activity, and a flattened head shape. Loss of cranial mesenchyme-restricted Wntless (CM-Wls), a gene required for secretion of all Wnt ligands, led to diminished apical expansion of OSX+ calvarial osteoblasts in the frontal bone primordia in a non-cell autonomous manner without perturbing proliferation or survival. Calvarial osteoblast polarization, progressive cell elongation and enrichment for actin cytoskeleton protein along the baso-apical axis were dependent on CM-Wnts. Thus, CM-Wnts regulate cellular behaviors during calvarial morphogenesis and provide tissue level cues for efficient apical expansion of calvarial osteoblasts. These findings also offer potential insights into the etiologies of calvarial dysplasias.

Interspersed expression of CUP-SHAPED COTYLEDON2 and REDUCED COMPLEXITY shapes Cardamine hirsuta complex leaf form

How genetically regulated growth shapes organ form is a key problem in developmental biology. Here, we investigate this problem using the leaflet-bearing complex leaves of Cardamine hirsuta as a model. Leaflet development requires the action of two growth-repressing transcription factors: REDUCED COMPLEXITY (RCO), a homeodomain protein, and CUP-SHAPED COTYLEDON2 (CUC2), a NAC-domain protein. However, how their respective growth-repressive actions are integrated in space and time to generate complex leaf forms remains unknown. By using live imaging, we show that CUC2 and RCO are expressed in an interspersed fashion along the leaf margin, creating a distinctive striped pattern. We find that this pattern is functionally important because forcing RCO expression in the CUC2 domain disrupts auxin-based marginal patterning and can abolish leaflet formation. By combining genetic perturbations with time-lapse imaging and cellular growth quantifications, we provide evidence that RCO-mediated growth repression occurs after auxin-based leaflet patterning and in association with the repression of cell proliferation. Additionally, through the use of genetic mosaics, we show that RCO is sufficient to repress both cellular growth and proliferation in a cell-autonomous manner. This mechanism of growth repression is different to that of CUC2, which occurs in proliferating cells. Our findings clarify how the two growth repressors RCO and CUC2 coordinate to subdivide developing leaf primordia into distinct leaflets and generate the complex leaf form. They also indicate different relationships between growth repression and cell proliferation in the patterning and post-patterning stages of organogenesis.

Spatial consistency of cell growth direction during organ morphogenesis requires CELLULOSE SYNTHASE INTERACTIVE1

Extracellular matrices contain fibril-like polymers often organized in parallel arrays. Although their role in morphogenesis has been long recognized, it remains unclear how the subcellular control of fibril synthesis translates into organ shape. We address this question using the Arabidopsis sepal as a model organ. In plants, cell growth is restrained by the cell wall (extracellular matrix). Cellulose microfibrils are the main load-bearing wall component, thought to channel growth perpendicularly to their main orientation. Given the key function of CELLULOSE SYNTHASE INTERACTIVE1 (CSI1) in guidance of cellulose synthesis, we investigate the role of CSI1 in sepal morphogenesis. We observe that sepals from csi1 mutants are shorter, although their newest cellulose microfibrils are more aligned compared to wild-type. Surprisingly, cell growth anisotropy is similar in csi1 and wild-type plants. We resolve this apparent paradox by showing that CSI1 is required for spatial consistency of growth direction across the sepal.

Organization of the body wall in lampreys informs the evolution of the vertebrate paired appendages

Vertebrate paired appendages are one of the most important evolutionary novelties in vertebrates. During embryogenesis, the skeletal elements of paired appendages differentiate from the somatic mesoderm, which is a layer of lateral plate mesoderm. However, the presence of the somatic mesoderm in the common ancestor of vertebrates has been controversial. To address this problem, it is necessary but insufficient to understand the developmental process of lateral plate mesoderm formation in lamprey (jawless vertebrates) embryos. Here, I show the presence of the somatic mesoderm in lamprey (Lethenteron camtschaticum) embryos using plastic sectioning and transmission electron microscopy analysis. During the early pharyngeal stages, the somatic mesoderm transforms from the lateral plate mesoderm in the trunk region. Soon after, when the cardiac structures were morphologically distinct, the somatic mesoderm was recognized through the cardiac to more caudal regions. These findings indicated that the somatic mesoderm evolved before the emergence of paired appendages. I also discuss the developmental changes in the body wall organization in the common ancestor of vertebrates, which is likely related to the evolution of the paired appendages.

Planar cell polarity regulators in asymmetric organogenesis during development and disease

The phenomenon of planar cell polarity is critically required for a myriad of morphogenetic processes in metazoan and is accurately controlled by several conserved modules. Six "core" proteins, including Frizzled, Flamingo (Celsr), Van Gogh (Vangl), Dishevelled, Prickle, and Diego (Ankrd6), are major components of the Wnt/planar cell polarity pathway. The Fat/Dchs protocadherins and the Scrib polarity complex also function to instruct cellular polarization. In vertebrates, all these pathways are essential for tissue and organ morphogenesis, such as neural tube closure, left-right symmetry breaking, heart and gut morphogenesis, lung and kidney branching, stereociliary bundle orientation, and proximal-distal limb elongation. Mutations in planar polarity genes are closely linked to various congenital diseases. Striking advances have been made in deciphering their contribution to the establishment of spatially oriented pattern in developing organs and the maintenance of tissue homeostasis. The challenge remains to clarify the complex interplay of different polarity pathways in organogenesis and the link of cell polarity to cell fate specification. Interdisciplinary approaches are also important to understand the roles of mechanical forces in coupling cellular polarization and differentiation. This review outlines current advances on planar polarity regulators in asymmetric organ formation, with the aim to identify questions that deserve further investigation.

Concentration-Dependent Domain Evolution in Reaction-Diffusion Systems

Pattern formation has been extensively studied in the context of evolving (time-dependent) domains in recent years, with domain growth implicated in ameliorating problems of pattern robustness and selection, in addition to more realistic modelling in developmental biology. Most work to date has considered prescribed domains evolving as given functions of time, but not the scenario of concentration-dependent dynamics, which is also highly relevant in a developmental setting. Here, we study such concentration-dependent domain evolution for reaction-diffusion systems to elucidate fundamental aspects of these more complex models. We pose a general form of one-dimensional domain evolution and extend this to N-dimensional manifolds under mild constitutive assumptions in lieu of developing a full tissue-mechanical model. In the 1D case, we are able to extend linear stability analysis around homogeneous equilibria, though this is of limited utility in understanding complex pattern dynamics in fast growth regimes. We numerically demonstrate a variety of dynamical behaviours in 1D and 2D planar geometries, giving rise to several new phenomena, especially near regimes of critical bifurcation boundaries such as peak-splitting instabilities. For sufficiently fast growth and contraction, concentration-dependence can have an enormous impact on the nonlinear dynamics of the system both qualitatively and quantitatively. We highlight crucial differences between 1D evolution and higher-dimensional models, explaining obstructions for linear analysis and underscoring the importance of careful constitutive choices in defining domain evolution in higher dimensions. We raise important questions in the modelling and analysis of biological systems, in addition to numerous mathematical questions that appear tractable in the one-dimensional setting, but are vastly more difficult for higher-dimensional models.

Resolving morphogenesis into quantifiable cell behaviours

Morphogenesis is extremely diverse, but its systematic quantification to determine the physical mechanisms that produce different phenotypes is possible by quantifying the underlying cell behaviours. These are limited and definable: they consist of cell proliferation, orientation of cell division, cell rearrangement, directional matrix production, cell addition/subtraction and cell size/shape change. Although minor variations in these categories are possible, in sum they capture all possible morphogenetic behaviours. This article summarises these processes, discusses their measurement, and highlights some salient examples.

The emerging significance of splicing in vertebrate development

Splicing is a crucial regulatory node of gene expression that has been leveraged to expand the proteome from a limited number of genes. Indeed, the vast increase in intron number that accompanied vertebrate emergence might have aided the evolution of developmental and organismal complexity. Here, we review how animal models for core spliceosome components have provided insights into the role of splicing in vertebrate development, with a specific focus on neuronal, neural crest and skeletal development. To this end, we also discuss relevant spliceosomopathies, which are developmental disorders linked to mutations in spliceosome subunits. Finally, we discuss potential mechanisms that could underlie the tissue-specific phenotypes often observed upon spliceosome inhibition and identify gaps in our knowledge that, we hope, will inspire further research.

4D reconstruction of murine developmental trajectories using spherical harmonics

Normal organogenesis cannot be recapitulated in vitro for mammalian organs, unlike in species including Drosophila and zebrafish. Available 3D data in the form of ex vivo images only provide discrete snapshots of the development of an organ morphology. Here, we propose a computer-based approach to recreate its continuous evolution in time and space from a set of 3D volumetric images. Our method is based on the remapping of shape data into the space of the coefficients of a spherical harmonics expansion where a smooth interpolation over time is simpler. We tested our approach on mouse limb buds and embryonic hearts. A key advantage of this method is that the resulting 4D trajectory can take advantage of all the available data while also being able to interpolate well through time intervals for which there are little or no data. This allows for a quantitative, data-driven 4D description of mouse limb morphogenesis.

•

Computer-based method recreating a 3D plus time evolution of a set of volumetric images

•

Technique based on the interpolation of the coefficients of spherical harmonics

•

Data-driven quantitative 4D description of limb and heart morphogenesis

•

Quantitatively reliable baseline description of organ development

Growth orientations, rather than heterogeneous growth rates, dominate jaw joint morphogenesis in the larval zebrafish

In early limb embryogenesis, synovial joints acquire specific shapes which determine joint motion and function. The process by which the opposing cartilaginous joint surfaces are moulded into reciprocal and interlocking shapes, called joint morphogenesis, is one of the least understood aspects of joint formation and the cell-level dynamics underlying it are yet to be unravelled. In this research, we quantified key cellular dynamics involved in growth and morphogenesis of the zebrafish jaw joint and synthesised them in a predictive computational simulation of joint development. Cells in larval zebrafish jaw joints labelled with cartilage markers were tracked over a 48-h time window using confocal imaging. Changes in distance and angle between adjacent cell centroids resulting from cell rearrangement, volume expansion and extracellular matrix (ECM) deposition were measured and used to calculate the rate and direction of local tissue deformations. We observed spatially and temporally heterogeneous growth patterns with marked anisotropy over the developmental period assessed. There was notably elevated growth at the level of the retroarticular process of the Meckel's cartilage, a feature known to undergo pronounced shape changes during zebrafish development. Analysis of cell dynamics indicated a dominant role for cell volume expansion in growth, with minor influences from ECM volume increases and cell intercalation. Cell proliferation in the joint was minimal over the timeframe of interest. Synthesising the dynamic cell data into a finite element model of jaw joint development resulted in accurate shape predictions. Our biofidelic computational simulation demonstrated that zebrafish jaw joint growth can be reasonably approximated based on cell positional information over time, where cell positional information derives mainly from cell orientation and cell volume expansion. By modifying the input parameters of the simulation, we were able to assess the relative contributions of heterogeneous growth rates and of growth orientation. The use of uniform rather than heterogeneous growth rates only minorly impacted the shape predictions, whereas isotropic growth fields resulted in altered shape predictions. The simulation results suggest that growth anisotropy is the dominant influence on joint growth and morphogenesis. This study addresses the gap of the cellular processes underlying joint morphogenesis, with implications for understanding the aetiology of developmental joint disorders such as developmental dysplasia of the hip and arthrogryposis.

Mechanical Regulation of Limb Bud Formation

Early limb bud development has been of considerable interest for the study of embryological development and especially morphogenesis. The focus has long been on biochemical signalling and less on cell biomechanics and mechanobiology. However, their importance cannot be understated since tissue shape changes are ultimately controlled by active forces and bulk tissue rheological properties that in turn depend on cell-cell interactions as well as extracellular matrix composition. Moreover, the feedback between gene regulation and the biomechanical environment is still poorly understood. In recent years, novel experimental techniques and computational models have reinvigorated research on this biomechanical and mechanobiological side of embryological development. In this review, we consider three stages of early limb development, namely: outgrowth, elongation, and condensation. For each of these stages, we summarize basic biological regulation and examine the role of cellular and tissue mechanics in the morphogenetic process.

THREE-DIMENSIONAL COMPUTATIONAL MODEL OF EARLY UPPER LIMB DEVELOPMENT

Limb development begins during embryogenesis when a series of biochemical interactions are triggered between a particular region of the mesoderm and the ectoderm. These processes affect the morphogenesis and growth of bones, joints, and all the other constituent elements of limbs; nevertheless, how the biochemical regulation affects mesenchymal condensation is not entirely clear. In this study, a three-dimensional computational model is designed to predict the appearance and location of the mesenchymal condensation in the stylopod and zeugopod; the biochemical events were described with reaction–diffusion equations that were solved using the finite elements method. The result of the gene expression in our model was consistent with the one reported in literature; the obtained patterns of Fgf8, Fgf10, and Wnt3a can predict the shape of the mesenchymal condensation of early upper limb development; the simple diffusive patterns of molecules were suitable to explain the areas where sox9 is expressed. Furthermore, our results suggest that the expression of Tgf-[Formula: see text] in the upper limb could be due to the inhibition of retinoic acid. These results suggest the importance of building computational scenarios where pathologies may be comprehensively examined.

Normal skeletal pattern formation in chick limb bud with a mesenchymal hole is mediated by adjustment of cellular properties along the anterior-posterior axis in the limb bud

The chick limb bud has plasticity to reconstruct a normal skeletal pattern after a part of mesenchymal mass is excised to make a hole in its early stage of development. To understand the details of hole closure and re-establishment of normal limb axes to reconstruct a normal limb skeleton, we focused on cellular and molecular changes during hole repair and limb restoration. We excised a cube-shaped mass of mesenchymal cells from the medial region of chick hindlimb bud (stage 23) and observed the following morphogenesis. The hole had closed by 15 ​h after excision, followed by restoration of the limb bud morphology, and the cartilage pattern was largely restored by 48 ​h. Lineage analysis of the mesenchymal cells showed that cells at the anterior and posterior margins of the hole were adjoined at the hole closure site, whereas cells at the proximal and distal margins were not. To investigate cell polarity during hole repair, we analyzed intracellular positioning of the Golgi apparatus relative to the nuclei. We found that the Golgi apparatus tended to be directed toward the hole among cells at the anterior and posterior margins but not among cells at identical positions in normal limb buds or cells at the proximal and distal hole margins. In the manipulated limb buds, the frequency of cell proliferation was maintained compared with the control side. Tbx3 expression, which was usually restricted to anterior and posterior margins of the limb bud, was temporarily expanded medially and then reverted to a normal pattern as limb reconstruction proceeded, with Tbx3 negative cells reappearing in the medial regions of the limb buds. Thus, mesenchymal hole closure and limb reconstruction are mainly mediated by cells at the anterior and posterior hole margins. These results suggest that adjustment of cellular properties along the anteroposterior axis is crucial to restore limb damage and reconstruct normal skeletal patterns.

Programming cell growth into different cluster shapes using diffusible signals

Advances in genetic engineering technologies have allowed the construction of artificial genetic circuits, which have been used to generate spatial patterns of differential gene expression. However, the question of how cells can be programmed, and how complex the rules need to be, to achieve a desired tissue morphology has received less attention. Here, we address these questions by developing a mathematical model to study how cells can collectively grow into clusters with different structural morphologies by secreting diffusible signals that can influence cellular growth rates. We formulate how growth regulators can be used to control the formation of cellular protrusions and how the range of achievable structures scales with the number of distinct signals. We show that a single growth inhibitor is insufficient for the formation of multiple protrusions but may be achieved with multiple growth inhibitors, and that other types of signals can regulate the shape of protrusion tips. These examples illustrate how our approach could potentially be used to guide the design of regulatory circuits for achieving a desired target structure.

Quantification of cell behaviors and computational modeling show that cell directional behaviors drive zebrafish pectoral fin morphogenesis

MOTIVATION

Understanding the mechanisms by which the zebrafish pectoral fin develops is expected to produce insights on how vertebrate limbs grow from a 2D cell layer to a 3D structure. Two mechanisms have been proposed to drive limb morphogenesis in tetrapods: a growth-based morphogenesis with a higher proliferation rate at the distal tip of the limb bud than at the proximal side, and directed cell behaviors that include elongation, division and migration in a non-random manner. Based on quantitative experimental biological data at the level of individual cells in the whole developing organ, we test the conditions for the dynamics of pectoral fin early morphogenesis.

RESULTS

We found that during the development of the zebrafish pectoral fin, cells have a preferential elongation axis that gradually aligns along the proximodistal (PD) axis of the organ. Based on these quantitative observations, we build a center-based cell model enhanced with a polarity term and cell proliferation to simulate fin growth. Our simulations resulted in 3D fins similar in shape to the observed ones, suggesting that the existence of a preferential axis of cell polarization is essential to drive fin morphogenesis in zebrafish, as observed in the development of limbs in the mouse, but distal tip-based expansion is not.

AVAILABILITYAND IMPLEMENTATION

Upon publication, biological data will be available at http://bioemergences.eu/modelingFin, and source code at https://github.com/guijoe/MaSoFin.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dynamic changes in epithelial cell packing during tissue morphogenesis

Cell packing - the spatial arrangement of cells - determines the shapes of organs. Recently, investigations of organ development in a variety of model organisms have uncovered cellular mechanisms that are used by epithelial tissues to change cell packing, and thereby their shapes, to generate functional architectures. Here, we review these cellular mechanisms across a wide variety of developmental processes in vertebrates and invertebrates and identify a set of common motifs in the morphogenesis toolbox that, in combination, appear to allow any change in tissue shape. We focus on tissue elongation, folding and invagination, and branching. We also highlight how these morphogenetic processes are achieved by cell-shape changes, cell rearrangements, and oriented cell division. Finally, we describe approaches that have the potential to engineer three-dimensional tissues for both basic science and translational purposes. This review provides a framework for future analyses of how tissues are shaped by the dynamics of epithelial cell packing.

Secreted inhibitors drive the loss of regeneration competence in Xenopus limbs

Absence of a specialized wound epidermis is hypothesized to block limb regeneration in higher vertebrates. However, the factors preventing its formation in regeneration-incompetent animals are poorly understood. To characterize the endogenous molecular and cellular regulators of specialized wound epidermis formation in Xenopus laevis tadpoles, and the loss of their regeneration competency during development, we used single-cell transcriptomics and ex vivo regenerating limb cultures. Transcriptomic analysis revealed that the specialized wound epidermis is not a novel cell state, but a re-deployment of the apical-ectodermal-ridge (AER) programme underlying limb development. Enrichment of secreted inhibitory factors, including Noggin, a morphogen expressed in developing cartilage/bone progenitor cells, are identified as key inhibitors of AER cell formation in regeneration-incompetent tadpoles. These factors can be overridden by Fgf10, which operates upstream of Noggin and blocks chondrogenesis. These results indicate that manipulation of the extracellular environment and/or chondrogenesis may provide a strategy to restore regeneration potential in higher vertebrates.

Summary: Secreted inhibitors associated with chondrogenic progression inhibit AER cell formation and restrict limb regeneration potential.

Mechanical control of tissue shape and morphogenetic flows during vertebrate body axis elongation

Shaping embryonic tissues into their functional morphologies requires cells to control the physical state of the tissue in space and time. While regional variations in cellular forces or cell proliferation have been typically assumed to be the main physical factors controlling tissue morphogenesis, recent experiments have revealed that spatial variations in the tissue physical (fluid/solid) state play a key role in shaping embryonic tissues. Here we theoretically study how the regional control of fluid and solid tissue states guides morphogenetic flows to shape the extending vertebrate body axis. Our results show that both the existence of a fluid-to-solid tissue transition along the anteroposterior axis and the tissue surface tension determine the shape of the tissue and its ability to elongate unidirectionally, with large tissue tensions preventing unidirectional elongation and promoting blob-like tissue expansions. We predict both the tissue morphogenetic flows and stresses that enable unidirectional axis elongation. Our results show the existence of a sharp transition in the structure of morphogenetic flows, from a flow with no vortices to a flow with two counter-rotating vortices, caused by a transition in the number and location of topological defects in the flow field. Finally, comparing the theoretical predictions to quantitative measurements of both tissue flows and shape during zebrafish body axis elongation, we show that the observed morphogenetic events can be explained by the existence of a fluid-to-solid tissue transition along the anteroposterior axis. These results highlight the role of spatiotemporally-controlled fluid-to-solid transitions in the tissue state as a physical mechanism of embryonic morphogenesis.

The Story of the Hand

This review describes the Story of the Human Hand. It traces the functional needs that led to evolution of the human hand as well as its embryological development. The various in utero stages of formation of the human hand are covered along with a description of the various molecular and genetic factors that control this process.

FACEts of mechanical regulation in the morphogenesis of craniofacial structures

During embryonic development, organs undergo distinct and programmed morphological changes as they develop into their functional forms. While genetics and biochemical signals are well recognized regulators of morphogenesis, mechanical forces and the physical properties of tissues are now emerging as integral parts of this process as well. These physical factors drive coordinated cell movements and reorganizations, shape and size changes, proliferation and differentiation, as well as gene expression changes, and ultimately sculpt any developing structure by guiding correct cellular architectures and compositions. In this review we focus on several craniofacial structures, including the tooth, the mandible, the palate, and the cranium. We discuss the spatiotemporal regulation of different mechanical cues at both the cellular and tissue scales during craniofacial development and examine how tissue mechanics control various aspects of cell biology and signaling to shape a developing craniofacial organ.

On the evolution and development of morphological complexity: A view from gene regulatory networks

How does morphological complexity evolve? This study suggests that the likelihood of mutations increasing phenotypic complexity becomes smaller when the phenotype itself is complex. In addition, the complexity of the genotype-phenotype map (GPM) also increases with the phenotypic complexity. We show that complex GPMs and the above mutational asymmetry are inevitable consequences of how genes need to be wired in order to build complex and robust phenotypes during development. We randomly wired genes and cell behaviors into networks in EmbryoMaker. EmbryoMaker is a mathematical model of development that can simulate any gene network, all animal cell behaviors (division, adhesion, apoptosis, etc.), cell signaling, cell and tissues biophysics, and the regulation of those behaviors by gene products. Through EmbryoMaker we simulated how each random network regulates development and the resulting morphology (i.e. a specific distribution of cells and gene expression in 3D). This way we obtained a zoo of possible 3D morphologies. Real gene networks are not random, but a random search allows a relatively unbiased exploration of what is needed to develop complex robust morphologies. Compared to the networks leading to simple morphologies, the networks leading to complex morphologies have the following in common: 1) They are rarer; 2) They need to be finely tuned; 3) Mutations in them tend to decrease morphological complexity; 4) They are less robust to noise; and 5) They have more complex GPMs. These results imply that, when complexity evolves, it does so at a progressively decreasing rate over generations. This is because as morphological complexity increases, the likelihood of mutations increasing complexity decreases, morphologies become less robust to noise, and the GPM becomes more complex. We find some properties in common, but also some important differences, with non-developmental GPM models (e.g. RNA, protein and gene networks in single cells).

WNT5A-Ca2+-CaN-NFAT signalling plays a permissive role during cartilage differentiation in embryonic chick digit development

During digit development, the correct balance of chondrogenic signals ensures the recruitment of undifferentiated cells into the cartilage lineage or the maintenance of cells at the undifferentiated stage. WNT/β catenin maintains the pool of progenitor cells, whereas TGFβ signalling promotes cartilage differentiation by inducing Sox9 expression. Moreover, WNT5A promotes the degradation of β catenin during mouse limb development. Although these mechanisms are well established, it is still unknown whether the signalling pathway downstream WNT5A is also involved in early chondrogenesis during digit formation. Thus, the aim of this study was to determine the role of WNT5A during the recruitment of progenitor cells during digit development. Our results showed that WNT5A activated calcium (Ca²⁺) release in the undifferentiated region during digit development. Further, the blockade of Ca²⁺ release or calcineurin (CaN) or nuclear factor of activated T-cells (NFAT) functions resulted in an inhibition of cartilage differentiation. Together, our results demonstrate that non canonical WNT5A-Ca²⁺-CaN-NFAT signalling plays a key role during embryonic digit development in vivo promoting the competence for chondrogenic signals and also acts as a permissive factor for chondrogenesis independently of cell death mechanisms.

Bone Morphogenetic Protein Signaling Restricts Proximodistal Extension of the Ventral Fin Fold

Unpaired fins, which are the most ancient form of locomotory appendages in chordates, had emerged at least 500 million years ago. While it has been suggested that unpaired fins and paired fins share structural similarities, cellular and molecular mechanisms that regulate the outgrowth of the former have not been fully elucidated yet. Using the ventral fin fold in zebrafish as a model, here, we investigate how the outgrowth of the unpaired fin is modulated. We show that Bone Morphogenetic Protein (BMP) signaling restricts extension of the ventral fin fold along the proximodistal axis by modulating diverse aspects of cellular behaviors. We find that lack of BMP signaling, either caused by genetic or chemical manipulation, prolongs the proliferative capacity of epithelial cells and substantially increases the number of cells within the ventral fin fold. In addition, inhibition of BMP signaling attenuates the innate propensity of cell division along the anteroposterior axis and shifts the orientation of cell division toward the proximodistal axis. Moreover, abrogating BMP signaling appears to induce excessive distal migration of cells within the ventral fin fold, and therefore precipitates extension along the proximodistal axis. Taken together, our data suggest that BMP signaling restricts the outgrowth of the ventral fin fold during zebrafish development.

Of mitogens and morphogens: modelling Sonic Hedgehog mechanisms in vertebrate development

Sonic Hedgehog (Shh) Is a critical protein in vertebrate development, orchestrating patterning and growth in many developing systems. First described as a classic morphogen that patterns tissues through a spatial concentration gradient, subsequent studies have revealed a more complex mechanism, in which Shh can also regulate proliferation and differentiation. While the mechanism of action of Shh as a morphogen is well understood, it remains less clear how Shh might integrate patterning, proliferation and differentiation in a given tissue, to ultimately direct its morphogenesis. In tandem with experimental studies, mathematical modelling can help gain mechanistic insights into these processes and bridge the gap between Shh-regulated patterning and growth, by integrating these processes into a common theoretical framework. Here, we briefly review the roles of Shh in vertebrate development, focusing on its functions as a morphogen, mitogen and regulator of differentiation. We then discuss the contributions that modelling has made to our understanding of the action of Shh and highlight current challenges in using mathematical models in a quantitative and predictive way. This article is part of a discussion meeting issue 'Contemporary morphogenesis'.

IRX3/5 regulate mitotic chromatid segregation and limb bud shape

Pattern formation is influenced by transcriptional regulation as well as by morphogenetic mechanisms that shape organ primordia, although factors that link these processes remain under-appreciated. Here we show that, apart from their established transcriptional roles in pattern formation, IRX3/5 help to shape the limb bud primordium by promoting the separation and intercalation of dividing mesodermal cells. Surprisingly, IRX3/5 are required for appropriate cell cycle progression and chromatid segregation during mitosis, possibly in a nontranscriptional manner. IRX3/5 associate with, promote the abundance of, and share overlapping functions with co-regulators of cell division such as the cohesin subunits SMC1, SMC3, NIPBL and CUX1. The findings imply that IRX3/5 coordinate early limb bud morphogenesis with skeletal pattern formation.

Establishing the pattern of the vertebrate limb

The vertebrate limb continues to serve as an influential model of growth, morphogenesis and pattern formation. With this Review, we aim to give an up-to-date picture of how a population of undifferentiated cells develops into the complex pattern of the limb. Focussing largely on mouse and chick studies, we concentrate on the positioning of the limbs, the formation of the limb bud, the establishment of the principal limb axes, the specification of pattern, the integration of pattern formation with growth and the determination of digit number. We also discuss the important, but little understood, topic of how gene expression is interpreted into morphology.

Loss of U11 small nuclear RNA in the developing mouse limb results in micromelia

Disruption of the minor spliceosome due to mutations in RNU4ATAC is linked to primordial dwarfism in microcephalic osteodysplastic primordial dwarfism type 1, Roifman syndrome, and Lowry-Wood syndrome. Similarly, primordial dwarfism in domesticated animals is linked to positive selection in minor spliceosome components. Despite being vital for limb development and size regulation, its role remains unexplored. Here, we disrupt minor spliceosome function in the developing mouse limb by ablating one of its essential components, U11 small nuclear RNA, which resulted in micromelia. Notably, earlier loss of U11 corresponded to increased severity. We find that limb size is reduced owing to elevated minor intron retention in minor intron-containing genes that regulate cell cycle. As a result, limb progenitor cells experience delayed prometaphase-to-metaphase transition and prolonged S-phase. Moreover, we observed death of rapidly dividing, distally located progenitors. Despite cell cycle defects and cell death, the spatial expression of key limb patterning genes was maintained. Overall, we show that the minor spliceosome is required for limb development via size control potentially shared in disease and domestication.

A "Numerical Evo-Devo" Synthesis for the Identification of Pattern-Forming Factors

Animals display extensive diversity in motifs adorning their coat, yet these patterns have reproducible orientation and periodicity within species or groups. Morphological variation has been traditionally used to dissect the genetic basis of evolutionary change, while pattern conservation and stability in both mathematical and organismal models has served to identify core developmental events. Two patterning theories, namely instruction and self-organisation, emerged from this work. Combined, they provide an appealing explanation for how natural patterns form and evolve, but in vivo factors underlying these mechanisms remain elusive. By bridging developmental biology and mathematics, novel frameworks recently allowed breakthroughs in our understanding of pattern establishment, unveiling how patterning strategies combine in space and time, or the importance of tissue morphogenesis in generating positional information. Adding results from surveys of natural variation to these empirical-modelling dialogues improves model inference, analysis, and in vivo testing. In this evo-devo-numerical synthesis, mathematical models have to reproduce not only given stable patterns but also the dynamics of their emergence, and the extent of inter-species variation in these dynamics through minimal parameter change. This integrative approach can help in disentangling molecular, cellular and mechanical interaction during pattern establishment.

Shaping Organs: Shared Structural Principles Across Kingdoms

Development encapsulates the morphogenesis of an organism from a single fertilized cell to a functional adult. A critical part of development is the specification of organ forms. Beyond the molecular control of morphogenesis, shape in essence entails structural constraints and thus mechanics. Revisiting recent results in biophysics and development, and comparing animal and plant model systems, we derive key overarching principles behind the formation of organs across kingdoms. In particular, we highlight how growing organs are active rather than passive systems and how such behavior plays a role in shaping the organ. We discuss the importance of considering different scales in understanding how organs form. Such an integrative view of organ development generates new questions while calling for more cross-fertilization between scientific fields and model system communities.

Self-sustained planar intercalations due to mechanosignaling feedbacks lead to robust axis extension during morphogenesis

Tissue elongation is a necessary process in metazoans to implement their body plans that is not fully understood. Here we propose a mechanism based on the interplay between cellular mechanics and primordia patterning that results in self-sustained planar intercalations. Thus, we show that a location-dependent modulation of the mechanical properties of cells leads to robust axis extension. To illustrate the plausibility of this mechanism, we test it against different patterning models by means of computer simulations of tissues where we implemented mechano-signaling feedbacks. Our results suggest that robust elongation relies on a trade-off between cellular and tissue strains that is orchestrated through the cleavage orientation. In the particular context of axis extension in Turing-patterned tissues, we report that different directional cell activities cooperate synergetically to achieve elongation. Altogether, our findings help to understand how the axis extension phenomenon emerges from the dynamics of individual cells.

On the problem of biological form

Embryonic development, which inspired the first theories of biological form, was eventually excluded from the conceptual framework of the Modern Synthesis as irrelevant. A major question during the last decades has centred on understanding whether new advances in developmental biology are compatible with the standard view or whether they compel a new theory. Here, I argue that the answer to this question depends on which concept of morphogenesis is held. Morphogenesis can be conceived as (1) a chemically driven or (2) a mechanically driven process. According to the first option, genetic regulatory networks drive morphogenesis. According to the second, morphogenesis results from an invariant tendency of embryonic tissues to restore changes in mechanical stress. While chemically driven morphogenesis allows an extension of the standard view, mechanically driven morphogenesis would deeply transform it. Which of these hypotheses has wider explanatory power is unknown. At present, the problem of biological form remains unsolved.

Multiscale modeling of vertebrate limb development

We review the current state of mathematical modeling of cartilage pattern formation in vertebrate limbs. We place emphasis on several reaction-diffusion type models that have been proposed in the last few years. These models are grounded in more detailed knowledge of the relevant regulatory processes than previous ones but generally refer to different molecular aspects of these processes. Considering these models in light of comparative phylogenomics permits framing of hypotheses on the evolutionary order of appearance of the respective mechanisms and their roles in the fin-to-limb transition. This article is categorized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Mechanistic Models Developmental Biology > Developmental Processes in Health and Disease Analytical and Computational Methods > Analytical Methods.

Mathematical modeling of chondrogenic pattern formation during limb development: Recent advances in continuous models

The phenomenon of chondrogenic pattern formation in the vertebrate limb is one of the best studied examples of organogenesis. Many different models, mathematical as well as conceptual, have been proposed for it in the last fifty years or so. In this review, we give a brief overview of the fundamental biological background, then describe in detail several models which aim to describe qualitatively and quantitatively the corresponding biological phenomena. We concentrate on several new models that have been proposed in recent years, taking into account recent experimental progress. The major mathematical tools in these approaches are ordinary and partial differential equations. Moreover, we discuss models with non-local flux terms used to account for cell-cell adhesion forces and a structured population model with diffusion. We also include a detailed list of gene products and potential morphogens which have been identified to play a role in the process of limb formation and its growth.

Finite element analysis reveals an important role for cell morphology in response to mechanical compression

Mechanical loading naturally controls cell phenotype, development, motility and various other biological functions; however, prolonged or substantial loading can cause cell damage and eventual death. Loading-induced mechanobiological and mechanostructural responses of different cell types affect their morphology and the internal architecture and the mechanics of the cellular components. Using single, mesenchymal stem cells, we have developed a cell-specific three-dimensional finite-element model; cell models were developed from phase-contrast microscopy images. This allowed us to evaluate the mechanostructural response of the naturally occurring variety of cell morphologies to increase sustained compressive loading. We focus on the morphology of the cytoplasm and the nucleus, as the main mechanically responsive elements, and evaluate formation of tensional strains and area changes in cells undergoing increasing uniaxial compressions. Here, we study mesenchymal stem cells as a model, due to their important role in tissue engineering and regenerative medicine; the method and findings are, however, applicable to any cell type. We observe variability in the cell responses to compression, which correlate directly with the morphology of the cells. Specifically, in cells with or without elongated protrusions (i.e., lamellipodia) tensional strains were, respectively, distributed mostly in the thin extensions or concentrated around the stiff nucleus. Thus, through cell-specific computational modeling of mechanical loading we have identified an underlying cause for stiffening (by actin recruitment) along the length of lamellipodia as well as a role for cell morphology in inducing cell-to-cell variability in mechanostructural response to loading.

Planar cell polarity signaling regulates polarized second heart field morphogenesis to promote both arterial and venous pole septation

The second heart field (SHF) harbors progenitors that are important for heart formation, but little is known about its morphogenesis. We show that SHF population in the mouse splanchnic mesoderm (SpM-SHF) undergoes polarized morphogenesis to preferentially elongate anteroposteriorly. Loss of Wnt5, a putative ligand of the planar cell polarity (PCP) pathway, causes the SpM-SHF to expand isotropically. Temporal tracking reveals that the Wnt5a lineage is a unique subpopulation specified as early as E7.5, and undergoes bi-directional deployment to form specifically the pulmonary trunk and the dorsal mesenchymal protrusion (DMP). In Wnt5a-/- mutants, Wnt5a lineage fails to extend into the arterial and venous poles, leading to both outflow tract and atrial septation defects that can be rescued by an activated form of PCP effector Daam1. We identify oriented actomyosin cables in the medial SpM-SHF as a potential Wnt5a-mediated mechanism that promotes SpM-SHF lengthening and restricts its widening. Finally, the Wnt5a lineage also contributes to the pulmonary mesenchyme, suggesting that Wnt5a/PCP is a molecular circuit recruited by the recently identified cardiopulmonary progenitors to coordinate morphogenesis of the pulmonary airways and the cardiac septations necessary for pulmonary circulation.This article has an associated 'The people behind the papers' interview.

Molecular signatures identify immature mesenchymal progenitors in early mouse limb buds that respond differentially to morphogen signaling

The key molecular interactions governing vertebrate limb bud development are a paradigm for studying the mechanisms controlling progenitor cell proliferation and specification during vertebrate organogenesis. However, little is known about the cellular heterogeneity of the mesenchymal progenitors in early limb buds that ultimately contribute to the chondrogenic condensations prefiguring the skeleton. We combined flow cytometric and transcriptome analyses to identify the molecular signatures of several distinct mesenchymal progenitor cell populations present in early mouse forelimb buds. In particular, jagged 1 (JAG1)-positive cells located in the posterior-distal mesenchyme were identified as the most immature limb bud mesenchymal progenitors (LMPs), which crucially depend on SHH and FGF signaling in culture. The analysis of gremlin 1 (Grem1)-deficient forelimb buds showed that JAG1-expressing LMPs are protected from apoptosis by GREM1-mediated BMP antagonism. At the same stage, the osteo-chondrogenic progenitors (OCPs) located in the core mesenchyme are already actively responding to BMP signaling. This analysis sheds light on the cellular heterogeneity of the early mouse limb bud mesenchyme and on the distinct response of LMPs and OCPs to morphogen signaling.

Oscillatory cortical forces promote three dimensional cell intercalations that shape the murine mandibular arch

Multiple vertebrate embryonic structures such as organ primordia are composed of confluent cells. Although mechanisms that shape tissue sheets are increasingly understood, those which shape a volume of cells remain obscure. Here we show that 3D mesenchymal cell intercalations are essential to shape the mandibular arch of the mouse embryo. Using a genetically encoded vinculin tension sensor that we knock-in to the mouse genome, we show that cortical force oscillations promote these intercalations. Genetic loss- and gain-of-function approaches show that Wnt5a functions as a spatial cue to coordinate cell polarity and cytoskeletal oscillation. These processes diminish tissue rigidity and help cells to overcome the energy barrier to intercalation. YAP/TAZ and PIEZO1 serve as downstream effectors of Wnt5a-mediated actomyosin polarity and cytosolic calcium transients that orient and drive mesenchymal cell intercalations. These findings advance our understanding of how developmental pathways regulate biophysical properties and forces to shape a solid organ primordium.

Physical control of tissue morphogenesis across scales

During embryogenesis, tissues and organs are progressively shaped into their functional morphologies. While the information about tissue and organ shape is encoded genetically, the sculpting of embryonic structures in the 3D space is ultimately a physical process. The control of physical quantities involved in tissue morphogenesis originates at cellular and subcellular scales, but it is their emergent behavior at supracellular scales that guides morphogenetic events. In this review, we highlight the physical quantities that can be spatiotemporally tuned at supracellular scales to sculpt tissues and organs during embryonic development of animal species, and connect them to the cellular and molecular mechanisms controlling them.

Strain-triggered mechanical feedback in self-organizing optic-cup morphogenesis

Organogenesis is a self-organizing process of multiple cells in three-dimensional (3D) space, where macroscopic tissue deformations are robustly regulated by multicellular autonomy. It is clear that this robust regulation requires cells to sense and modulate 3D tissue formation across different scales, but its underlying mechanisms are still unclear. To address this question, we developed a versatile computational model of 3D multicellular dynamics at single-cell resolution and combined it with the 3D culture system of pluripotent stem cell-derived optic-cup organoid. The complementary approach enabled quantitative prediction of morphogenesis and its corresponding verification and elucidated that the macroscopic 3D tissue deformation is fed back to individual cellular force generations via mechanosensing. We hereby conclude that mechanical force plays a key role as a feedback regulator to establish the robustness of organogenesis.

A quantitative analysis of 3D-cell distribution in regenerating muscle-skeletal system with synchrotron X-ray computed microtomography

One of the greatest enigmas of modern biology is how the geometry of muscular and skeletal structures are created and how their development is controlled during growth and regeneration. Scaling and shaping of vertebrate muscles and skeletal elements has always been enigmatic and required an advanced technical level in order to analyse the cell distribution in 3D. In this work, synchrotron X-ray computed microtomography (µCT) and chemical contrasting has been exploited for a quantitative analysis of the 3D-cell distribution in tissues of a developing salamander (Pleurodeles waltl) limb – a key model organism for vertebrate regeneration studies. We mapped the limb muscles, their size and shape as well as the number and density of cells within the extracellular matrix of the developing cartilage. By using tomographic approach, we explored the polarity of the cells in 3D, in relation to the structure of developing joints. We found that the polarity of chondrocytes correlates with the planes in joint surfaces and also changes along the length of the cartilaginous elements. Our approach generates data for the precise computer simulations of muscle-skeletal regeneration using cell dynamics models, which is necessary for the understanding how anisotropic growth results in the precise shapes of skeletal structures.