PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(5): e10892.
Published online May 28, 2010. doi:  10.1371/journal.pone.0010892
PMCID: PMC2878345
Bare Bones Pattern Formation: A Core Regulatory Network in Varying Geometries Reproduces Major Features of Vertebrate Limb Development and Evolution
Jianfeng Zhu,1 Yong-Tao Zhang,1 Mark S. Alber,1,2 and Stuart A. Newman3*
1Department of Mathematics, University of Notre Dame, Notre Dame, Indiana, United States of America
2Center for the Study of Biocomplexity, University of Notre Dame, Notre Dame, Indiana, United States of America
3Department of Cell Biology and Anatomy, New York Medical College, Valhalla, New York, United States of America
Mark Isalan, Editor
Center for Genomic Regulation, Spain
* E-mail: newman/at/nymc.edu
Conceived and designed the experiments: JZ YTZ MSA SAN. Performed the experiments: JZ YTZ. Analyzed the data: MSA SAN. Wrote the paper: YTZ MSA SAN.
Received March 4, 2010; Accepted May 7, 2010.
Background
Major unresolved questions regarding vertebrate limb development concern how the numbers of skeletal elements along the proximodistal (P-D) and anteroposterior (A-P) axes are determined and how the shape of a growing limb affects skeletal element formation. There is currently no generally accepted model for these patterning processes, but recent work on cartilage development (chondrogenesis) indicates that precartilage tissue self-organizes into nodular patterns by cell-molecular circuitry with local auto-activating and lateral inhibitory (LALI) properties. This process is played out in the developing limb in the context of a gradient of fibroblast growth factor (FGF) emanating from the apical ectodermal ridge (AER).
Results
We have simulated the behavior of the core chondrogenic mechanism of the developing limb in the presence of an FGF gradient using a novel computational environment that permits simulation of LALI systems in domains of varying shape and size. The model predicts the normal proximodistal pattern of skeletogenesis as well as distal truncations resulting from AER removal. Modifications of the model's parameters corresponding to plausible effects of Hox proteins and formins, and of the reshaping of the model limb, bud yielded simulated phenotypes resembling mutational and experimental variants of the limb. Hypothetical developmental scenarios reproduce skeletal morphologies with features of fossil limbs.
Conclusions
The limb chondrogenic regulatory system operating in the presence of a gradient has an inherent, robust propensity to form limb-like skeletal structures. The bare bones framework can accommodate ancillary gene regulatory networks controlling limb bud shaping and establishment of Hox expression domains. This mechanism accounts for major features of the normal limb pattern and, under variant geometries and different parameter values, those of experimentally manipulated, genetically aberrant and evolutionary early forms, with no requirement for an independent system of positional information.
The limbs of vertebrate animals emerge from the embryonic flank as buds of somatopleure-derived mesenchymal cells covered by an epithelial layer, the ectoderm, flattening into paddle shapes as they grow. The most thoroughly studied aspect of limb development is the formation of the skeleton, an array of jointed bone or cartilage elements having a stereotypical pattern that has sustained only modest alterations over the course of evolution [1], [2]. The mechanism of limb skeletal pattern formation is incompletely understood. There is broad agreement, however, concerning the cellular and molecular-genetic interactions underlying the differentiation of the cartilage tissue that forms the embryonic primordia of the bony skeleton of tetrapod limbs and the endoskeleton of fish fins (see [3], [4] for reviews).
A major question regarding limb development concerns how the number and positioning of skeletal elements along the proximodistal (P-D) and anteroposterior (A-P) axes is determined [4]. There is a general P-D increase in the number of skeletal elements which occurs even in cases, such as the chicken forelimb (Fig. 1), in which the A-P width remains essentially constant while the skeletal pattern is being laid out. In fish, non-tetrapod vertebrates, the fin endoskeleton is a mixture of bars and nodules which have no discernable P-D numerical trend [5].
Figure 1
Figure 1
Relationship between core cartilage patterning network and “bare bones” framework for limb development.
Cartilage differentiation, or chondrogenesis, is preceded by “condensation” of the precartilage mesenchyme, in which cell density increases and the cells enter into broad, transient contact with one another [6]. The precartilage cells are embedded in a dilute extracellular matrix (ECM) and condensation is accompanied by and dependent on local accumulation of the ECM molecule fibronectin [7], with markers of prospective condensation appearing earlier than ECM and morphological changes [8]. Molecules secreted by the dorsal and ventral ectoderm, including FGFs and Wnt, inhibit chondrogenesis [8], [9], thereby confining the developing one-bar proximal cartilage primordium (stylopod, i.e., humerus, femur) to a central planar sector of the paddle-shaped limb bud [4]. As development proceeds, the skeleton remains confined to this plane but expands laterally in more distal regions as the stylopod gives way to the two-bar (zeugopod, i.e., radius and ulna, tibia and fibula) and multiple-bar (autopod, i.e., digits) primordia of the mid and terminal regions of the limb (Fig. 1C). This brings the developing skeletal elements increasingly closer to the anterior and posterior edges of the limb bud, reflecting attenuation of the peripheral inhibitory effect. Attenuation of inhibition can also be seen in the proximity of the more distal elements to the dorsal and ventral surfaces as the limb bud tapers towards its tip, and to the apical boundary, as the potency of the AER wanes [10].
There is a large literature that treats limb pattern formation as a question of establishing spatiotemporal informational fields that cells subsequently interpret according to their genetic programs [11], [12]. Other approaches to this problem (and those of other developing systems), take into account the conditional, self-organizational properties of interacting cells and tissues (see refs. [13] and [14] for reviews). This latter approach has the advantage over the “positional information” perspective in its ability to formulate hypotheses for why the limbs (or other embryonic structures) assume the forms they actually do.
A mechanism that has gained increasing attention as a generator of spot- and stripe-like patterns in developmental systems is related to the chemical reaction-diffusion process described by Turing in 1952 [15], [16]. This has been schematized in a biological “local autoactivation-lateral inhibition” (LALI) framework by Meinhardt and Gierer [17]. LALI systems, while formally similar to reaction-diffusion systems, are more suitable to biological applications, since they include cases where the activator and inhibitor terms are mediated by cellular “reactors” rather than simple chemical reactions [18], and spatial transport can be mediated by mechanisms in addition to simple diffusion [19]. LALI systems are self-organizing: that is, with appropriate choices of activator-inhibitor and transport parameters, in domains of appropriate size, shape and boundary conditions, a spatially uniform distribution of the morphogens becomes dynamically unstable, giving way to nonuniform distributions, typically periodic arrangements of stripes and spots in two dimensions, and their three-dimensional analogs (bars, nodules) [13], [20].
The developing limb bud contains all the components of a skeletogenic LALI mechanism [21], [22], [23]. The locus of activity of this process is a population of multipotent skeletogenic progenitor cells (termed “precartilage” cells in what follows) [24], [25], [26], which are maintained in a non-terminally differentiated state in the distal-most 0.3 mm of the avian and mammalian limb bud by the action of the apical ectodermal ridge (AER), a narrow A-P-oriented ectodermal thickening, along the limb bud tip [27]. Whereas the precartilage cells can give rise to cartilage, connective tissue, and bone, and can undergo programmed cell death, the limb musculature is formed by a separate population of cells [26]. The non-myogenic limb mesenchyme thus forms the skeletal pattern by being caused to choose among its potential fates of cartilage, soft connective tissue fibroblasts and apoptotic cells, in a spatiotemporally regulated fashion.
Precartilage condensation and subsequent chondrogenesis are promoted by morphogens of the TGF-β superfamily [28], [29] that form an “activator subnetwork” (Fig. 1B). In brief, one or more TGF-βs and activins set off a train of events whereby BMPs (one or more of BMP2, 4 and 7), acting via the receptor BMPR-Ia, induce regions of high BMP signaling activity, marked by phosphorylated Smads [30]. Other classes of molecules, such as galactose-binding lectins (“galectins”) [31], may also be involved in this initiation step. Signaling by the activator subnetwork induces extracellular matrix, matricellular and cell adhesion molecules, promoting mesenchymal condensation and subsequent chondrogenesis [32], [33] (Fig. 1B).
The spatial profile of BMP action is dictated not by its receptor distribution, which is uniform, but by a “prepattern” (i.e., a molecular pattern prefiguring the cellular one) of the diffusible morphogen itself, which in the autopod takes the form of narrow crescents at the tips of the forming digits [30]. This prepattern arises from the spatiotemporal dynamics of the activator network morphogens interacting with inhibitory factors. The latter include the extracellular antagonistic BMP binding partners Noggin, Gremlin and Ventroptin (CHL2), the BMP receptor antagonist BAMBI and the antagonistic intracellular co-receptor Smad6, the expression of which are all induced by the activator subnetwork [30], [34]. Consistent with this, the prospective digits and interdigits are uncommitted as to their fates late into autopod formation [30], and this plasticity likely obtains in the fields of precartilage cells that form the stylopod and zeugopod at earlier stages. In particular, exogenous TGF-β can cause noncondensing limb bud precartilage mesenchyme to condense and form ectopic cartilage in vitro [28] and in vivo [34].
The FGF pathway [8] and the Notch pathway [35], [36] are involved in restricting the expansion of precartilage condensations as they form. Along with the inhibitory factors described above, these perinodular chondrogenesis-restricting components constitute what we refer to as the “inhibitor subnetwork” (I) (Fig. 1B).
The activator (A) and inhibitor (I) subnetworks together constitute a core network for the regulation of chondrogenesis (Fig. 1B). All cells of the precartilage mesenchyme appear to be equally capable of producing all components of A and I. Subnetwork A functions in a paracrine and (in part via its TGF-β component [29]), a positively autoregulatory fashion. It induces the production of cell-to-cell transmissible components of subnetwork I, which in turn limit subnetwork A's radius of action (Fig. 1B).
Limb bud mesenchyme can in fact self-organize into cartilage patterns with limb-like features. Dissociated and reaggregated limb mesoderm cells, packed into an ectodermal hull, form limb-like arrays of rod-like and nodular skeletal elements [37]. In high-density culture these cells form spots and stripes of cartilage with dimensions and spacing similar to the in vivo patterns [38]. When the physical properties of the culture microenvironment [39] or putative LALI components are manipulated in vitro and vivo, the pattern changes are predictable from LALI models [8], [28], [29], [30].
In most amniote species, emergence of the cartilage primordia of the limb skeleton occurs first in proximal regions of the limb bud and then in increasingly more distal regions [40] (Figs. 1; ;2,2, left). This process is dependent on an intact AER, the removal of which leads to truncations of the skeleton with progressively more distal elements forming the longer the AER remains in place [40]). The growth of the skeletal structures (and portions thereof) that are spared after AER removal is unaffected, with these components attaining nearly their normal size and shapes, except for terminal deletions, as development proceeds [40]. The limbs of certain amphibians seem to represent partial exceptions to the proximal-to-distal rule [41].
Figure 2
Figure 2
Simulation of chicken wing development.
Although the idea that a LALI mechanism could account for major features of limb development in its normal, genetically and experimentally perturbed, and evolutionary aspects is several decades old [23], what has been lacking for most of the succeeding period has been a fast computational modeling environment that would enable simulation of a LALI system as a continuous (“real-time”) process, in domains with natural contours and changing geometries. There has been progress in this direction more recently, using finite element algorithms, producing modeling environments suitable for investigating time- and shape-dependent changes in patterns generated by reaction-diffusion systems of biological interest [42], [43]. Here we use one such framework to test the capability of a LALI limb model to generate realistic patterns under normal and biologically altered conditions.
Below, we analyze the dynamics of a mathematical representation of a limb LALI system in the context of changing domains of realistic proportions, under conditions that take into account uncertainties in the values of parameters such as biosynthetic and diffusion rates and interaction strengths. We show that the contextualized LALI mechanism accounts for major features of the normal limb pattern as well as those of experimentally manipulated, genetically aberrant and evolutionary transitional forms with no a priori requirement for an independent system of positional information.
Mathematical representations of the limb skeletogenic system
Simulations of simple LALI systems under geometric and growth constraints partially approximating those of developing limbs show that pattern formation of normal and mutant (e.g., polydactyly, human Apert syndrome, mouse Doublefoot) limbs is consistent with such mechanisms [22], [44]. These simulations, however, have variously employed highly schematic and ad hoc molecular and cellular representations, and stationary approximations of inherently temporal dynamics [21], [22]. Up till now, unavailability of adequate computational methods has prevented the continuous time-dependent simulation of a realistic LALI-type system for limb development.
The model of Hentschel and coworkers [18] is the most biologically detailed of the LALI mechanisms for limb development that have been considered. It consists of eight coupled partial differential equations that model all of the following: reaction-diffusion dynamics of a network based on a TGF-β-related activator and a FGF/FGF receptor 2-induced inhibitor, diffusion of FGF from fixed and moving sources, TGF-β-induced production of fibronectin, and the short-range movement of cells and their differentiation from uncondensed mesenchyme to condensed mesenchyme, and then into cartilage, including regulatory dynamics of the associated FGF receptor types. While it would be desirable to follow the time evolution of this system in an initially unpatterned 3D domain that mimics the growing limb bud, the complexity of this prospective simulation exceeds the capacity of available methods.
We have simplified this problem using both biologically motivated assumptions and mathematical arguments [45]. According to a recent classification of developmental mechanisms [46], the limb model of Hentschel et al. [18] is “morphodynamic,” since differentiation of new cell types occurs simultaneously with cell rearrangement. Many developmental mechanisms, in contrast, conform to a simpler “morphostatic” scenario, in which cell identity becomes established independently of cell rearrangement. We have assumed for the purposes of this study that the short-range movement involved in precartilage mesenchymal condensation occurs in populations of cells that have already been determined to do so by developmental signals such as the activating and inhibitory morphogens described above. While the morphostatic assumption is difficult to prove experimentally, it is biologically reasonable. Moreover, the mechanistic separation of cell movement and differentiation can be dispensed with in more elaborate modeling strategies as computational power and methods improve.
Employing the morphostatic limit has allowed us to extract a LALI activator-inhibitor subnetwork from the 8-equation system described in ref. [18]:
A mathematical equation, expression, or formula.
 Object name is pone.0010892.e001.jpg
(1)
Here An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e002.jpg represents the concentration of activating morphogen and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e003.jpg the concentration of inhibitor. The parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e004.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e005.jpg are the respective coefficients of (effective) diffusion, An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e006.jpg a decay or breakdown constant, and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e007.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e008.jpg functions that define the production of the morphogens based on the cross- and auto-feedback relationships. Two reaction kinetic parameters incorporated in the functional forms of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e009.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e010.jpg, An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e011.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e012.jpg (see ref. [45] and File S1) represent factors that are plausibly altered in value by modulatory transcription factors such as Hox and other gene regulatory proteins [47], and the effects of their variation are described below.
This biologically motivated mathematical simplification of the limb skeletogenic pattern forming system makes it amenable to set of new computational methods for simulation of its time-dependent behavior in regions of arbitrary shape ([42]; also see File S1). In the following sections we describe pattern development in a growing limb-bud-like domain with no imposed structure or pattern other than a simple P-D gradient representing the effect of the AER [27] (Fig. 1C). Our analysis demonstrates that major features of normal, experimentally manipulated, genetically aberrant and evolutionary transitional limb development emerge from the inherent self-organizing properties of the core skeletal patterning mechanism. The roles of factors important in establishing asymmetries in the limb axes and differences among the skeletal elements, such as the morphogen Sonic hedgehog and the Hox gene families (reviewed in ref. [3]), have natural interpretations in this computational framework as regulators of limb bud shape and interaction parameters of the activator-inhibitor network.
Simulation of normal development
To simulate the developing limb we allowed the LALI system (1) to operate in a geometric setting that represents the section of the growing limb bud equidistant from the dorsal and ventral surfaces (Fig. 1C). A proximodistal gradient of FGF, representing the AER, is initially put in place on this template, with its source at the limb tip. The morphogens are constrained not to leak from the borders of the domain (“no-flux” boundary conditions) and the model limb is programmed to grow uniformly in the P-D direction.
As noted above, the AER keeps precartilage mesenchyme in a labile, undifferentiated state, an effect that has been suggested to be related to its ability to promote outgrowth and distal progression [27], [48]. Moreover, it can be replaced in both these functions by one or more FGFs [49], which diffuse through the distal limb bud mesenchyme [50], [51] forming a proximodistal gradient [52] and regulate cell number and survival [53]. Because the P-D length of the uncondensed distal tip of the chicken limb bud decreases during the period of pattern formation [54] as the potency of the AER wanes [10] we have incorporated into our model the hypothesis that a temporally declining gradient of FGF with its source in the AER determines the P-D length of this distal region, which we term the “LALI zone.” Although the curved distal boundary allows us to simulate the reaction-diffusion dynamics of the LALI zone in more realistic shapes, we did not use the boundary to define a curved source of FGF. Instead the distribution of FGF is modeled as a 1D gradient varying purely along the PD axis, whose value is highest at the distal-most point of the simulation domain. This gradient is represented not by an explicit diffusion equation, but parametrically, by a function that decreases nonlinearly from a maximal value at the distal tip and in which the tip value itself decreases nonlinearly with time (see File S2).
Depending on the geometry of the LALI zone, the activator concentration profile (Fig. 2, right, color-scale) takes the form of one or more elongated spots, analogous to the crescents observed for BMP in vivo [30] (see also Fig. 1C). We define the morphogenetically active region of the LALI zone (the “active zone”; violet bands) as that portion which is at a sufficient distance from the AER to escape the suppressive effect of FGF (see also Fig. 1C). In the full pattern-forming model, cells in the active zone exposed to peak values of activator produce fibronectin and other condensation-promoting molecules [18]. Here we impose a simpler rule by which such cells directly and stably differentiate into cartilage (blue) once elongation of the limb brings them outside the LALI zone. The skeletal elements in this “frozen zone” (with respect to pattern formation) elongate by both growth and addition from newly emerging active zone tissue.
We caused the AER to wane in potency by decreasing the maximal value of the FGF gradient over time. This reduces the P-D length of the LALI zone. Based on the properties of the LALI system (1), this contraction by itself can lead to abrupt increases in the numbers of parallel elements along the A-P axis [45] (see also [18], [23]). We obtained simulation results with a more authentic time-course and proportions (Fig. 2, right; Movie S1) by also changing the parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e013.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e014.jpg (components of the functions An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e015.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e016.jpg) in a stage-dependent fashion (i.e., when the stylopod, zeugopod and autopod are forming) (see File S3). These parameter changes leave the network topology of the reaction-diffusion system unaltered but modify the biosynthetic responses of the morphogens to themselves and each other, and indirectly [45], of ECM production to the morphogens, similarly to demonstrated and proposed functions of the Hox family of transcription factors [47], [55].
Establishment of appropriate expression domains of Hox proteins, some of which change in abundance in the distal mesenchyme in a stage-dependent fashion [56], [57], although not essential for the development of arrays of skeletal elements, are necessary for refined, species-characteristic patterns [58]. Correspondingly, in the parameter range of the standard model the stage-specific variations of the parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e017.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e018.jpg (along with contraction of the LALI zone), are essential to a normal developmental outcome (File S3; see Discussion).
Robustness of the pattern forming mechanism
Varying key parameters in the LALI system (1) and the geometric properties of the limb simulation template (Fig. 1C) allowed us to explore their effects on the generated pattern and the robustness of the entire process. The reaction kinetic parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e019.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e020.jpg (see ref. [45] and File S1), and the A-P width and P-D length of the LALI zone, are key determinants of the character of the patterns formed. This can be illustrated by the effect of changes in these dimensionless parameters (Fig. 3) in the context of the standard developmental sequence (Fig. 2, right; Movie S1). During the first phase of simulated development (up to 1.4 arbitrary time units), as the LALI zone shrinks in the P-D direction to 2/3 its original length, one (the normal result), two, or three stripes will form when An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e021.jpg takes on a particular one (±7–11%) of three separated values over a two-fold range, coupled with low, middle or high values (±1%) of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e022.jpg in an overlapping range (Fig. 3).
Figure 3
Figure 3
Dependence of the number of stripes on the kinetic parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e023.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e024.jpg during the first phase of normal development.
Based on an extensive set of simulations using altered limb shapes (see below), we have found that the model is much more tolerant to variations in An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e031.jpg than in An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e032.jpg. This disparity could be anticipated from the critical role of a narrow range of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e033.jpg values in the capacity of the system to exhibit pattern formation [45].
The change of parameters in the different phases of the simulation is a plausible computational implementation of the changing distributions of Hox gene products in the apical zone at the different phases of limb development [57]. Hox factors have multiple targets, including genes for cell attachment and ECM molecules, as well as morphogen receptors [47]. Our simulations of the standard developmental progression are sensitive to large and combined changes in parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e034.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e035.jpg, but robust against smaller and individual changes (Fig. 3). For this parameter sweep a single P-D length for the LALI zone was chosen, corresponding to the initial stage of the standard simulation (Fig. 2), when a single stripe is normally specified.
Simulation of AER removal experiments
Simulations using the standard parameter set, but with the AER removed part-way during development, show the same result as when this manipulation is performed on developing limbs [40]. In silico, as in vivo, removal of the AER is equivalent to eliminating the distal source of the FGF gradient. After a time lag during which the suppressive gradient decays, the entire LALI zone up to the limb bud tip becomes susceptible to being organized by the activator-inhibitor system, leaving the limb distally truncated (Fig. 4; Movies S2 and S3).
Figure 4
Figure 4
Simulations of AER removal.
Simulation of anteroposteriorly expanded limb buds
In contrast to the chicken forelimb, many limb buds have an apical zone that expands along the A-P axis during skeletogenic patterning (Fig 5, left). This occurs in the avian hindlimb, in reptilian and mammalian limbs, and to an even greater extent in pathological limbs, such as those of chickens which have received ectopic grafts of the zone of polarizing activity (ZPA) [59], embryos which bear the talpid2 mutation [60], or mice in which the Sonic hedgehog (Shh) and Gli3 genes have been jointly knocked down [61]. In the dogfish, while there is no autopod (this being considered an innovation distinguishing fish from tetrapods), a limb bud that is expanded relative to those of birds and mammals produces an array of parallel elements [5], which while not true digits, assume a similar configuration (Fig. 5, left).
Figure 5
Figure 5
Simulations of effect of limb bud distal expansion.
The changing limb bud profiles were represented by heuristic equations chosen to emulate plausible reshaping trajectories. Values of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e036.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e037.jpg were chosen for 2–7 phases of the continuous simulations, by trial and error to provide a good fit to the final pattern (see File S4).
We found that the number of parallel stripes in the A-P direction in the model limb bud increased in a discrete fashion with expansion in the length of the LALI zone, consistent with the general expectation for reaction-diffusion and other LALI systems (Fig. 6). Simulations in which the limb bud or its tip was allowed to expand in a fashion similar to that seen with ZPA grafts into the chicken forelimb, Shh/Gli3 null mouse embryos, and the chicken mutant talpid2, show an increased number of parallel “digits” in the autopod (Fig. 5, right; File S4).
Figure 6
Figure 6
Relation between A-P width and number of stripes.
Simulation of fossil limb skeletons
Many fossil vertebrate limbs with endoskeletons have been characterized, including some from extinct fish-like ancestors and quasi-tetrapods with authentic autopods [62] (Fig. 7, left). Dinosaur limbs resemble those of modern reptiles and mammals, but ichthyosaurs, swimming dinosaurs, are of interest in that their paddles consist largely of nodular rather than rod-like skeletal elements [63] (Fig. 7, left). We performed simulations based on a wide range of hypothetical developmental scenarios (the one leading to the ichthyosaur-like pattern is shown in File S5) in which the network topology and rules of the standard progression that generated Fig. 2 were maintained, but the contours of the distal region of the limb bud were varied during the simulations (though without prior knowledge of the actual limb bud profiles available for the simulation in Fig. 5), and the parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e041.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e042.jpg were changed in a stage-specific fashion (see File S5). The simulation end-points shown in Fig. 7 (right) indicate that our model exhibits sufficient flexibility to reproduce the general features of a wide variety of fossil limb skeletons.
Figure 7
Figure 7
Simulation of fossil limb skeletons.
The results presented here confirm that a reaction-diffusion process operating in a reshaping domain in the context of a distal suppressive gradient can reproduce major aspects of skeletal patterning in vertebrate limbs. In particular, the proximodistal temporal sequence of development of skeletal elements (Fig. 2), the proximodistally increasing number of elements in fully developed limbs (Fig. 2) and the dependence of development on the suppressive FGF gradient (i.e., the AER), with distal truncation occurring upon its removal (Fig. 4), are all readily accounted for by this mechanism. The “bare bones” patterning mechanism also produces increased numbers of elements when the apical domain is expanded in an anteroposterior direction during development, as occurs with ectopic grafts of the ZPA, with double knockouts of Shh and Gli3, development of the chicken talpid2 mutant or the dogfish (Fig. 5).
All of these simulations employed model limbs with curved apical contours, which are truer to the living limb shapes than the simple geometries generally used in reaction-diffusion simulations. The ability to perform such simulations in domains with natural contours and continuously changing nonstandard shapes was enabled by our previous development of new finite element computational methods which are applicable to any biological system in which fields of chemicals or other mobile agents are produced by activator-inhibitor interactions on domains of changing size and shape [42].
By our main hypothesis, limb skeletal patterning is established by a LALI process within a reshaping zone of tissue in the presence of a distal suppressive gradient. While the skeletal patterns in all the limbs we have simulated are mediated by the same self-organizing system, the changing limb bud shapes within which this system operates were imposed arbitrarily according to the schemes described in Files S3, S4, S5. Since limb bud shape is controlled by a different set of molecular determinants from the ones regulating the initiation of chondrogenesis (reviewed in ref. [2]), a more complete model would incorporate an independent system of equations for limb bud shaping [64], [65] involving the growth [66] and viscoelastic [67] properties of the limb bud tissue regulated by molecules such as Wnt, Shh, Gli3, and FGFs [68], [69].
In our model, an increase in the number of parallel elements occurs when either the P-D length of the LALI zone decreases ([45]; first figure of File S3) or its A-P width increases (Fig. 6). While certain other conditions may lead to peak splitting in activator-inhibitor systems [22], [70], we have seen only occasional examples of intercalation of new peaks between preexisting ones with parameter change (e.g., second figure of File S3). Since such patterning modes characterize certain variant limb types [22], this may indicate a limitation of our representation that would improve with the introduction of additional modulatory parameters corresponding to the molecular complexity of the biological system [3], [4], [11].
Joints can form in our model as a result of discontinuities in the pattern of elements along the AP axis as different spatial solutions become stable as the LALI zone changes in size and shape (e.g., Fig. 2; Movie S1). Oscillations are well known to occur in reaction-diffusion system in certain parameter regimes [71], including our morphostatic system (H.G.E. Hentschel and S.A. Newman, unpublished data). Although such oscillations would not generally serve to segment stripe-like LALI pattern elements [70], the elongated elements in our system are typically formed by temporal persistence of spot-like patterns rather than as de novo stripes (e.g., Brachypterygius simulation in File S5), making oscillatory modes plausible factors in joint formation. In embryonic limbs the generation of joints depends on members of the BMP superfamily beyond those specified in our core mechanism, such as GDF5 [72]. These factors have activating and inhibitory effects on chondrogenesis that would result in more complex spatiotemporal waves than those seen in the basic model.
In all the simulations shown we have varied the values of the reaction kinetic parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e043.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e044.jpg in a P-D level-specific fashion; that is, different values were used for the stylopod, zeugopod and autopod (see Files S3, S4, S5 for details). Although we determined previously that the contraction of the P-D length of the LALI zone (a consequence of the attenuation of the AER suppressive gradient), was by itself sufficient to increase in the number of parallel skeletal elements [45], the variation of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e045.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e046.jpg (which does not change the network topology of the core patterning system nor impart any element-specific positional information to the model limb), fine tuned the auto- and cross-regulatory interactions between the morphogens and led, with appropriate choices to relatively authentic skeletal patterns.
The parameters An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e047.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e048.jpg were defined by a detailed mathematical analysis [45] in which factors relating to cell movement and extracellular matrix production were folded into functions (U and V) governing the activating and inhibitory morphogens operating within a system encompassing a fuller range of biological phenomena [18]. This simplification was necessitated by limitations on the complexity of systems that can be simulated by available finite element algorithms [42]. While attribution of molecular functions to the parameters in the original eight-equation system is more straightforward, it is nonetheless still possible to discern the roles of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e049.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e050.jpg in the morphostatic system (1) and thereby attempt to put a “molecular face” on these parameters.
Both parameters appear in the production rates of the activator and inhibitor morphogens, and An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e051.jpg is also the association rate constant of the activator and the inhibitor. The parameter An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e052.jpg denotes the activator morphogen concentration at which there is a transition from a linear response to a saturation phase. In particular, the Turing bifurcation (a transition in the behavioral characteristics of the system) that enables pattern formation only occurs in a narrow range of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e053.jpg values [45]. This constraint is in evidence in all the simulations shown here (with displacement of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e054.jpg±2% around its mean value significantly affecting the pattern), and can be considered a required relationship for pattern formation among the system's activator and inhibitory components under our biological assumptions and their mathematical representations.
The limb deformity (ld) locus in the mouse encodes several functionalities corresponding to the role of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e055.jpg in our model. Formin1, disruption of which leads to the absence of the fibula and reduction in digit number, probably regulates expression of BMPs [73], components of the activator subnetwork (Fig. 1b), whereas a transcriptional global control region (GCR) in this domain activates the limb-specific expression of the BMP antagonist Gremlin [74], a component of the inhibitor subnetwork. Not included in the present model, however, are the buffering systems accumulated over 400 million years of evolution that protect the biological equivalent of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e056.jpg from deviating from its prescribed range of values in present-day limbs.
In contrast to An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e057.jpg, the value of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e058.jpg can vary extensively without compromising the pattern-forming capacity of the system. Variations in An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e059.jpg lead to extensive changes in the number and arrangement of skeletal elements, particularly in the autopod, the most variable region of normal, mutant and fossil limbs. As suggested above, the behavior of the system under parameter variation may provide insight into the elusive functions of Hox transcription factors in the developing limb. Different members of the Hoxa family vary in abundance in the limb tip in a stage-dependent fashion [57], and like the parameter An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e060.jpg, influence activator dynamics. Hoxa13, for example, which is the most highly expressed member in the prospective autopod, alters the expression in the developing limb of two morphogens of the TGF-β superfamily, members of which comprise the activator subnetwork (Fig. 1b) [75]. Hoxa9, abundantly expressed in more proximal regions, modulates TGF-β superfamily signaling [76]. We can thus tentatively identify the role of An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e061.jpg in our scheme with certain functions of the Hoxa gene products.
A different set of Hox gene products, the Hoxd class, vary along the limb bud A-P axis during development [57]. A more detailed model could potentially account for the morphological differences between the radius and ulna, or the various fingers, by parameter variations in this direction as well. As with the dynamical network for limb bud shaping which is not yet included in our framework (see above), equations for setting parameter values in a biological fashion, ultimately based on the Hox gene regulatory network [58], [77], would fill an important gap in the model. It should be noted, however, that such shape and parameter regulatory systems presuppose a core chondrogenic mechanism (Fig. 1B) and a bare bones pattern generating scheme (Fig. 1C), as described above. They do not by themselves provide positional information to the cells of the developing limb.
Despite the extensive sampling of parameter space represented by the simulation results shown, and scores of additional simulations with other parameter combinations, some in which the parameter An external file that holds a picture, illustration, etc.
Object name is pone.0010892.e062.jpg varied by more than 10-fold (Files S3, S4, S5), it was difficult to obtain an outcome that did not resemble a limb skeleton. This in no way means that any arbitrary pattern could be produced: the simulated skeletons are all composed of spot- and stripe-like elements. The LALI system (1), in the context in which we have simulated it, is therefore inherently “skeletogenic.” It appears, moreover, that the topology of the core network, rather than the specific identity of the relevant gene products (of which our knowledge is incomplete, and for which there is variation between tetrapod species), may be decisive for this process.
A recent study has provided evidence that Tiktaalik [78] and other vertebrates known from the fossil record to have noncanonical bony limb skeletons previously thought to be transitional to definitive tetrapod limbs actually coexisted with early tetrapods [79]. The drastically narrowed the time span during which evolution of all known limb morphotypes must have occurred calls for a skeletogenic mechanism with a propensity to generate a profusion of patterns due to small genetic changes affecting limb bud shaping and the rates and strengths of core interactions. The mechanism we propose is the only empirically based one currently under consideration that has these properties.
While our model provides a plausible account of the general form of the limb skeletal pattern (a capability absent in competing models), its main role at present is as a framework for further experimental tests. In particular, it will be important to gain knowledge of the earliest acting activators and inhibitors of precartilage condensation, the details of their transmission and interaction, and the role of FGF or other signals from the AER in defining the morphogenetically active region of the developing limb.
Numerical simulations of the LALI system (1) were performed using a new class of finite element algorithms on moving and deforming domains [42], [80] based on the discontinuous Galerkin (DG) method [81], [82]. The DG method is a means for converting an ordinary or partial differential equation system into to a problem represented by a system of algebraic equations in a more restricted space than that of the original system. It employs “independent” polynomials on every element to approximate the system's behavior in the restricted space and provides more flexibility than the continuous Galerkin method. These novel methods enable numerical solutions of system (1) and other LALI and reaction-diffusion systems on deforming and moving grids in domains with complicated geometries. To approximate the irregular geometries of limb buds, we used the cubic spline interpolation technique of de Boor [83]. This method approximates the curved boundary of a limb bud by piece-wise cubic polynomials and maintains global smoothness of the obtained approximation curve.
File S1
Description of the partial differential equation model.
(0.07 MB PDF)
File S2
Representation and characteristics of the FGF gradient.
(0.06 MB PDF)
File S3
Dependence of pattern on P-D length of LALI zone and reaction-diffusion parameters.
(0.09 MB PDF)
File S4
Details of limb bud expansion experiments and mutants.
(0.03 MB PDF)
File S5
Details of fossil simulations.
(0.10 MB PDF)
Movie S1
Simulation video of normal development.
(1.77 MB AVI)
Movie S2
Simulation video of early AER removal.
(0.48 MB AVI)
Movie S3
Simulation video of late AER removal.
(1.06 MB AVI)
Footnotes
Competing Interests: The authors have declared that no competing interests exist.
Funding: S. A. Newman acknowledges support from the National Science Foundation (NSF) Frontiers in Integrative Biological Research program (grant EF-0526854). Y.-T. Zhang's research is partially supported by NSF grant DMS-0810413 and the Oak Ridge Associated Universities (ORAU) Ralph E. Powe Junior Faculty Enhancement Award. This research of M. S. Alber is partially supported by NSF grant DMS-0800612. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
1. Darwin C. On the origin of species by means of natural selection, or, The preservation of favoured races in the struggle for life. London: J. Murray; 1859.
2. Newman SA, Müller GB. Origination and innovation in the vertebrate limb skeleton: an epigenetic perspective. J Exp Zoolog B Mol Dev Evol. 2005;304:593–609. [PubMed]
3. Newman SA, Bhat R. Activator-inhibitor dynamics of vertebrate limb pattern formation. Birth Defects Res C Embryo Today. 2007;81:305–319. [PubMed]
4. Yang Y. Growth and patterning in the limb: signaling gradients make the decision. Sci Signal. 2009;2:pe3. [PubMed]
5. Yonei-Tamura S, Abe G, Tanaka Y, Anno H, Noro M, et al. Competent stripes for diverse positions of limbs/fins in gnathostome embryos. Evolution & Development. 2008;10:737–745. [PubMed]
6. Thorogood PV, Hinchliffe JR. An analysis of the condensation process during chondrogenesis in the embyronic hind limb. J Embryol Exp Morphol. 1975;33:581–606. [PubMed]
7. Frenz DA, Jaikaria NS, Newman SA. The mechanism of precartilage mesenchymal condensation: a major role for interaction of the cell surface with the amino-terminal heparin-binding domain of fibronectin. Dev Biol. 1989;136:97–103. [PubMed]
8. Moftah MZ, Downie SA, Bronstein NB, Mezentseva N, Pu J, et al. Ectodermal FGFs induce perinodular inhibition of limb chondrogenesis in vitro and in vivo via FGF receptor 2. Dev Biol. 2002;249:270–282. [PubMed]
9. ten Berge D, Brugmann SA, Helms JA, Nusse R. Wnt and FGF signals interact to coordinate growth with cell fate specification during limb development. Development. 2008;135:3247–3257. [PMC free article] [PubMed]
10. Rubin L, Saunders JW., Jr Ectodermal-mesodermal interactions in the growth of limb buds in the chick embryo: constancy and temporal limits of the ectodermal induction. Dev Biol. 1972;28:94–112. [PubMed]
11. Tabin C, Wolpert L. Rethinking the proximodistal axis of the vertebrate limb in the molecular era. Genes Dev. 2007;21:1433–1442. [PubMed]
12. Zeller R, Lopez-Rios J, Zuniga A. Vertebrate limb bud development: moving towards integrative analysis of organogenesis. Nat Rev Genet. 2009;10:845–858. [PubMed]
13. Forgacs G, Newman SA. Biological physics of the developing embryo. Cambridge: Cambridge Univ. Press; 2005.
14. Schnell S, Maini PK, Newman TJ, Newman SA, editors. Multiscale models of developmental systems. San Diego, CA: Elsevier; 2008.
15. Turing AM. The chemical basis of morphogenesis. Phil Trans Roy Soc Lond B. 1952;237:37–72.
16. Maini PK, Baker RE, Chuong CM. Developmental biology. The Turing model comes of molecular age. Science. 2006;314:1397–1398. [PubMed]
17. Meinhardt H. Models of biological pattern formation: from elementary steps to the organization of embryonic axes. Curr Top Dev Biol. 2008;81:1–63. [PubMed]
18. Hentschel HG, Glimm T, Glazier JA, Newman SA. Dynamical mechanisms for skeletal pattern formation in the vertebrate limb. Proc R Soc Lond B Biol Sci. 2004;271:1713–1722. [PMC free article] [PubMed]
19. Lander AD. Morpheus unbound: reimagining the morphogen gradient. Cell. 2007;128:245–256. [PubMed]
20. Alber M, Glimm T, Hentschel HGE, Kazmierczak B, Newman SA. Stability of n-dimensional patterns in a generalized Turing system: implications for biological pattern formation. Nonlinearity. 2005;18:125–138.
21. Newman SA, Christley S, Glimm T, Hentschel HG, Kazmierczak B, et al. Multiscale models for vertebrate limb development. Curr Top Dev Biol. 2008;81:311–340. [PubMed]
22. Miura T, Shiota K, Morriss-Kay G, Maini PK. Mixed-mode pattern in Doublefoot mutant mouse limb–Turing reaction-diffusion model on a growing domain during limb development. J Theor Biol. 2006;240:562–573. [PubMed]
23. Newman SA, Frisch HL. Dynamics of skeletal pattern formation in developing chick limb. Science. 1979;205:662–668. [PubMed]
24. Newman SA. Lineage and pattern in the developing wing bud. In: Ede DA, Hinchliffe JR, Balls M, editors. Vertebrate Limb and Somite Morphogenesis. Cambridge: Cambridge University Press; 1977. pp. 181–197.
25. Newman SA. Fibroblast progenitor cells of the embryonic chick limb. J Embryol Exp Morphol. 1980;56:191–200. [PubMed]
26. Pearse RV, 2nd, Scherz PJ, Campbell JK, Tabin CJ. A cellular lineage analysis of the chick limb bud. Dev Biol. 2007;310:388–400. [PMC free article] [PubMed]
27. Kosher RA, Savage MP, Chan SC. In vitro studies on the morphogenesis and differentiation of the mesoderm subjacent to the apical ectodermal ridge of the embryonic chick limb-bud. J Embryol Exp Morphol. 1979;50:75–97. [PubMed]
28. Leonard CM, Fuld HM, Frenz DA, Downie SA, Massagué J, et al. Role of transforming growth factor-β in chondrogenic pattern formation in the embryonic limb: stimulation of mesenchymal condensation and fibronectin gene expression by exogenous TGF-β and evidence for endogenous TGF-β-like activity. Dev Biol. 1991;145:99–109. [PubMed]
29. Miura T, Shiota K. TGFβ2 acts as an “activator” molecule in reaction-diffusion model and is involved in cell sorting phenomenon in mouse limb micromass culture. Dev Dyn. 2000;217:241–249. [PubMed]
30. Montero JA, Lorda-Diez CI, Ganan Y, Macias D, Hurle JM. Activin/TGFβ and BMP crosstalk determines digit chondrogenesis. Dev Biol. 2008;321:343–356. [PubMed]
31. Matsutani E, Yamagata T. Chick endogenous lectin enhances chondrogenesis of cultured chick limb bud cells. Dev Biol. 1982;92:544–548. [PubMed]
32. Pizette S, Niswander L. BMPs are required at two steps of limb chondrogenesis: formation of prechondrogenic condensations and their differentiation into chondrocytes. Dev Biol. 2000;219:237–249. [PubMed]
33. Yoon BS, Pogue R, Ovchinnikov DA, Yoshii I, Mishina Y, et al. BMPs regulate multiple aspects of growth-plate chondrogenesis through opposing actions on FGF pathways. Development. 2006;133:4667–4678. [PubMed]
34. Chimal-Monroy J, Rodriguez-Leon J, Montero JA, Gañan Y, Macias D, et al. Analysis of the molecular cascade responsible for mesodermal limb chondrogenesis: sox genes and BMP signaling. Dev Biol. 2003;257:292–301. [PubMed]
35. Fujimaki R, Toyama Y, Hozumi N, Tezuka K. Involvement of Notch signaling in initiation of prechondrogenic condensation and nodule formation in limb bud micromass cultures. J Bone Miner Metab. 2006;24:191–198. [PubMed]
36. Pan Y, Liu Z, Shen J, Kopan R. Notch1 and 2 cooperate in limb ectoderm to receive an early Jagged2 signal regulating interdigital apoptosis. Dev Biol. 2005;286:472–482. [PubMed]
37. Ros MA, Lyons GE, Mackem S, Fallon JF. Recombinant limbs as a model to study homeobox gene regulation during limb development. Dev Biol. 1994;166:59–72. [PubMed]
38. Christley S, Alber MS, Newman SA. Patterns of mesenchymal condensation in a multiscale, discrete stochastic model. PLoS Comput Biol. 2007;3: (e76):0743–0753. [PMC free article] [PubMed]
39. Miura T, Shiota K. Extracellular matrix environment influences chondrogenic pattern formation in limb bud micromass culture: Experimental verification of theoretical models. Anat Rec. 2000;258:100–107. [PubMed]
40. Saunders JW., Jr The proximo-distal sequence of origin of the parts of the chick wing and the role of the ectoderm. J Exp Zool. 1948;108:363–402. [PubMed]
41. Franssen RA, Marks S, Wake D, Shubin N. Limb chondrogenesis of the seepage salamander, Desmognathus aeneus (Amphibia: Plethodontidae). J Morphol. 2005;265:87–101. [PubMed]
42. Zhu J, Zhang Y-T, Newman S, Alber M. Application of discontinuous Galerkin methods for reaction-diffusion systems in developmental biology. Journal of Scientific Computing. 2009;40:391–418.
43. Madzvamuse A, Maini PK, Wathen AJ. A moving grid finite element method for the simulation of pattern generation by Turing models on growing domains. J Scientif Comput. 2005;24:247–262.
44. Chaturvedi R, Huang C, Kazmierczak B, Schneider T, Izaguirre JA, et al. On multiscale approaches to three-dimensional modelling of morphogenesis. J R Soc Interface. 2005;2:237–253. [PMC free article] [PubMed]
45. Alber M, Glimm T, Hentschel HG, Kazmierczak B, Zhang YT, et al. The morphostatic limit for a model of skeletal pattern formation in the vertebrate limb. Bull Math Biol. 2008;70:460–483. [PubMed]
46. Salazar-Ciudad I, Jernvall J, Newman SA. Mechanisms of pattern formation in development and evolution. Development. 2003;130:2027–2037. [PubMed]
47. Svingen T, Tonissen KF. Hox transcription factors and their elusive mammalian gene targets. Heredity. 2006;97:88–96. [PubMed]
48. Gregg BC, Rowe A, Brickell PM, Wolpert L. Ectodermal inhibition of cartilage differentiation in micromass culture of chick limb bud mesenchyme in relation to gene expression and cell shape. Development. 1989;105:769–777. [PubMed]
49. Niswander L, Tickle C, Vogel A, Booth I, Martin GR. FGF-4 replaces the apical ectodermal ridge and directs outgrowth and patterning of the limb. Cell. 1993;75:579–587. [PubMed]
50. Li S, Anderson R, Reginelli AD, Muneoka K. FGF-2 influences cell movements and gene expression during limb development. J Exp Zool. 1996;274:234–247. [PubMed]
51. Vargesson N, Kostakopoulou K, Drossopoulou G, Papageorgiou S, Tickle C. Characterisation of hoxa gene expression in the chick limb bud in response to FGF. Dev Dyn. 2001;220:87–90. [PubMed]
52. Nikbakht N, McLachlan JC. A proximo-distal gradient of FGF-like activity in the embryonic chick limb bud. Cell Mol Life Sci. 1997;53:447–451. [PubMed]
53. Sun X, Mariani FV, Martin GR. Functions of FGF signalling from the apical ectodermal ridge in limb development. Nature. 2002;418:501–508. [PubMed]
54. Summerbell D. A descriptive study of the rate of elongation and differentiation of skeleton of the developing chick wing. J Embryol Exp Morphol. 1976;35:241–260. [PubMed]
55. Newman SA. Sticky fingers: Hox genes and cell adhesion in vertebrate limb development. Bioessays. 1996;18:171–174. [PubMed]
56. Yokouchi Y, Nakazato S, Yamamoto M, Goto Y, Kameda T, et al. Misexpression of Hoxa-13 induces cartilage homeotic transformation and changes cell adhesiveness in chick limb buds. Genes and Development. 1995;9:2509–2522. [PubMed]
57. Nelson CE, Morgan BA, Burke AC, Laufer E, DiMambro E, et al. Analysis of Hox gene expression in the chick limb bud. Development. 1996;122:1449–1466. [PubMed]
58. Zakany J, Duboule D. The role of Hox genes during vertebrate limb development. Curr Opin Genet Dev. 2007;17:359–366. [PubMed]
59. Saunders JW, Jr, Gasseling MT. Ectodermal-mesenchymal interactions in the origin of limb symmetry. In: Fleischmajer R, Billingham RE, editors. Epithelial-Mesenchymal Interactions. Baltimore: Williams and Wilkins; 1968. pp. 78–97.
60. Goetinck PF, Abbott UK. Studies on limb morphogenesis. I. Experiments with the polydactylous mutant, talpid2. J Exp Zool. 1964;155:161–170. [PubMed]
61. Litingtung Y, Dahn RD, Li Y, Fallon JF, Chiang C. Shh and Gli3 are dispensable for limb skeleton formation but regulate digit number and identity. Nature. 2002;418:979–983. [PubMed]
62. Shubin N, Tabin C, Carroll S. Deep homology and the origins of evolutionary novelty. Nature. 2009;457:818–823. [PubMed]
63. Motani R. On the evolution and homologies of ichthyopterygian forefins. Journal of Vertebrate Paleontology. 1999;19:28–41.
64. Dillon R, Othmer HG. A mathematical model for outgrowth and spatial patterning of the vertebrate limb bud. J Theor Biol. 1999;197:295–330. [PubMed]
65. Morishita Y, Iwasa Y. Growth based morphogenesis of vertebrate limb bud. Bull Math Biol. 2008;70:1957–1978. [PMC free article] [PubMed]
66. Morishita Y, Iwasa Y. Estimating the spatiotemporal pattern of volumetric growth rate from fate maps in chick limb development. Dev Dyn. 2009;238:415–422. [PubMed]
67. Damon BJ, Mezentseva NV, Kumaratilake JS, Forgacs G, Newman SA. Limb bud and flank mesoderm have distinct “physical phenotypes” that may contribute to limb budding. Dev Biol. 2008;321:319–330. [PubMed]
68. Benazet JD, Bischofberger M, Tiecke E, Goncalves A, Martin JF, et al. A self-regulatory system of interlinked signaling feedback loops controls mouse limb patterning. Science. 2009;323:1050–1053. [PubMed]
69. Mackem S, Lewandoski M. Limb development takes a measured step toward systems analysis. Sci Signal. 2009;2:pe33. [PubMed]
70. Meinhardt H. The algorithmic beauty of sea shells. Berlin, Heidelberg: Springer-Verlag Berlin Heidelberg; 2009.
71. Lefever R, Nicolis G, Prigogine I. On the occurrence of oscillations around the steady state in systems of chemical reactions far from equilibrium. J Chem Phys. 1967;47:1045–1047.
72. Settle SH, Jr, Rountree RB, Sinha A, Thacker A, Higgins K, et al. Multiple joint and skeletal patterning defects caused by single and double mutations in the mouse Gdf6 and Gdf5 genes. Dev Biol. 2003;254:116–130. [PubMed]
73. Zhou F, Leder P, Zuniga A, Dettenhofer M. Formin1 disruption confers oligodactylism and alters Bmp signaling. Hum Mol Genet. 2009;18:2472–2482. [PMC free article] [PubMed]
74. Pavel E, Zhao W, Powell KA, Weinstein M, Kirschner LS. Analysis of a new allele of limb deformity (ld) reveals tissue- and age-specific transcriptional effects of the Ld Global Control Region. Int J Dev Biol. 2007;51:273–281. [PubMed]
75. Perez WD, Weller CR, Shou S, Stadler HS. Survival of Hoxa13 homozygous mutants reveals a novel role in digit patterning and appendicular skeletal development. Dev Dyn. 2010;239:446–457. [PMC free article] [PubMed]
76. Bai S, Shi X, Yang X, Cao X. Smad6 as a transcriptional corepressor. J Biol Chem. 2000;275:8267–8270. [PubMed]
77. Deschamps J. Developmental biology. Hox genes in the limb: a play in two acts. Science. 2004;304:1610–1611. [PubMed]
78. Shubin NH, Daeschler EB, Jenkins FA., Jr The pectoral fin of Tiktaalik roseae and the origin of the tetrapod limb. Nature. 2006;440:764–771. [PubMed]
79. Niedzwiedzki G, Szrek P, Narkiewicz K, Narkiewicz M, Ahlberg PE. Tetrapod trackways from the early Middle Devonian period of Poland. Nature. 2010;463:43–48. [PubMed]
80. Zhu J, Zhang Y-T, Newman SA, Alber MS. A discontinuous Galerkin finite element model on moving grids for vertebrate limb pattern formation. Mathematical Modelling of Natural Phenomena. 2009;4:131–148.
81. Johnson C. Numerical solution of partial differential equations by the finite element method. Cambridge-New York-New Rochelle-Melbourne-Sydney: Cambridge University Press; 1987.
82. Cheng Y, Shu C-W. A discontinuous Galerkin finite element method for time dependent partial differential equations with higher order derivatives. Mathematics of Computation. 2008;77:699–730.
83. De Boor C. A practical guide to splines. New York: Springer; 2001.
Articles from PLoS ONE are provided here courtesy of
Public Library of Science