PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Dev Biol. Author manuscript; available in PMC 2010 October 15.
Published in final edited form as:
PMCID: PMC2811698
NIHMSID: NIHMS145522

Spots and stripes: Pleomorphic patterning of stem cells via p-ERK-dependent cell-chemotaxis shown by feather morphogenesis and mathematical simulation

Summary

A key issue in stem cell biology is the differentiation of homogeneous stem cells towards different fates which are also organized into desired configurations. Little is known about the mechanisms underlying the process of periodic patterning. Feather explants offer a fundamental and testable model in which multi-potential cells are organized into hexagonally arranged primordia and the spacing between primordia. Previous work explored roles of a Turing reaction-diffusion mechanism in establishing chemical patterns. Here we show that a continuum of feather patterns, ranging from stripes to spots, can be obtained when the level of p-ERK activity is adjusted with chemical inhibitors. The patterns are dose-dependent, tissue stage-dependent, and irreversible. Analyses show ERK-activity-dependent mesenchymal cell chemotaxis is essential for converting micro-signaling centers into stable feather primordia. A mathematical model based on short-range activation, long-range inhibition, and cell chemotaxis is developed and shown to simulate observed experimental results. This generic cell behavior model can be applied to model stem cell patterning behavior at large.

Keywords: pattern formation, feather morphogenesis, stem cells, placode, ERK, mathematical modeling, chemotaxis

Introduction

By definition, stem cells are a population of multi-potential or pluri-potential cells. In response to environmental signals they differentiate and assemble into organized tissues and organs. While much has been learned about the molecular signals that induce cell differentiation (Keller, 2005), less is known about the rules that govern their morphogenesis. This is vividly illustrated in the case of pluri-potential embryonic stem cells, as seen both in the differentiated, yet unorganized, embryoid body and in teratoma. Attempts to engineer multi-potential stem cells into organs have found the requirement for architecture diminished when stem cell products are released into the bloodstream (which does not require structural organization), as seen in bone marrow stem cells (Weissman, 2000). However, it becomes an acute issue when the architecture of the stem cell product is critical to its function, as is the case for skeletal regeneration and engineered hair follicles (Watt, 2001; Widelitz et al., 2006). Progress in the tissue engineering of stem cells has pushed this issue beyond basic interest and into the realm of practical ramifications. While analytical approaches have provided valuable data, enabling us to identify, for example, essential molecules for morphogenetic processes, the rules that govern spatial patterning remain elusive. We now need to pursue a systematic approach to comprehend how these molecular pathways work together to build the architecture of a tissue / organ (Chuong et al., 2006).

One of the most fundamental processes in molding organ architecture is the ability of cell populations to form periodically arranged spots and stripes (Fig. 1A). These patterns are striking when they are observed in the integument pigment pattern of leopards, tigers, fishes, etc., (Wolpert, 1971; Ball, 1999; Kondo, 2002), the segmented vertebrae (Dequéant and Pourquié, 2008), and skin appendage arrangement (Jiang et al., 2004). Theoretical models have been proposed to explain these phenomena (Liu et al., 2006; Murray, 2003; Maini et al., 2006; Meinhardt and Gierer, 2000). Some developmental systems may be based on combinatorial molecular coding which can be interpreted at the enhancer / transcription factor level (e.g., Drosophila segments; Small and Levine, 1991). Others may be based on molecular oscillations and wave fronts (e.g., somitogenesis; Pourquie, 2003). It is likely these complex patterning processes may involve epigenetic, stochastic, and self-organizing processes (Newman and Comper 1990; Jiang et al., 2004; Newman et al., 2008). Mechanisms involving the reaction and diffusion of chemical substances and the chemotactic response of cells to such chemicals have been proposed for a number of biological phenomena including the pattern forming behavior of the slime mold, Dictyostelium discoidium (Vasiev, 1994) and the bacterium E. coli (Zorzano, 2005), patterns on the integument of snakes and fish (Myerscough, 1991; Painter, 1999), in the induction of hairs and feathers (Nagorcka and Mooney, 1992; Jung et al., 1998) and feather branching (Harris et al., 2005; Yu et al., 2005). These principles also have been used to explain skeletal patterning in the limb bud (Hentschel et al., 2004; Kiskowski et al., 2004; Miura and Maini, 2004). However, much remains to be learned about how molecular pathway and cellular events are coupled to the patterning process at the tissue / organ level.

Figure 1
Expression of p-ERK and related molecules

The periodically arranged feather arrays on embryonic chicken skin (Fig. 1B; Lin et al., 2006) provide an excellent paradigm for studying this issue. In the feather, at the outset, both epithelia and mesenchyma are homogenous, i.e., every cell is multi-potential and has an equal probability of becoming bud or interbud (equipotent), as demonstrated by DiI labeling in a reconstitution experiment (Chuong et al., 1996; Jiang et al., 1999). When feathers start to form, this homogeneity is disrupted, leading to the emergence of new arrangements or structures. In the chick, feathers arise in a sequential manner, from the dorsal midline out to the lateral regions (Figure 7). The exquisite two-dimensional layout makes it easier to see experimentally driven alterations in feather arrangements. Skin epithelium and dermis recombination experiments show that the initial patterning signal arises from the mesenchyme (Sengel, 1976). However, the molecular networks and mechanisms involved in establishing the periodic pattern remain elusive. Experimental work has shown that members of the FGF family, such as FGFs 1, 2, 4, as well as noggin and follistatin, function as activators by promoting feather bud formation (Jung et al., 1998; Song et al., 1996; Song et al., 2004; Patel et al., 1999; Widelitz et al., 1996). On the other hand, the BMPs and Delta-1 function as inhibitors by suppressing bud formation (Jung et al., 1998; Noramly and Morgan, 1998; Crowe et al., 1998; Viallet et al., 1998). Since gene networks are likely to drive the complex patterning process, tilting the balance of activator or inhibitor concentrations may result in different feather patterns (Jiang et al., 2004).

Figure 7
Schematic summary showing the periodic patterning process during feather morphogenesis

While activators and inhibitors are involved in the spatial patterning of skin precursor cells (Jung et al., 1998; Sick et al., 2006), exactly how molecular events combine to produce the feather primordia array is unclear. We have previously proposed that the interaction of activators and inhibitors could lead to a stable chemical concentration pattern. However, a reaction-diffusion mechanism alone is not sufficient to account for the changes in cell density which lead to the progression from a homogeneous field to small cell clusters and, finally, to committed stable dermal condensations (Jiang et al., 1999). In addition, the complexity of the patterning process is further highlighted by the dramatically different phenotypes obtained when the timing of signal pathway perturbation is subtly altered (Widelitz et al., 1996; Drew et al., 2006). Therefore, with our current, novel observations we propose that, along with the reaction and diffusion of activators and inhibitors, a cell-chemotaxis mechanism is required to achieve mesenchymal condensation and complete the periodic pattern-forming process.

As mentioned previously, FGFs are involved in feather pattern formation (Song et al., 1996; Song et al., 2004; Mandler and Neubuser, 2004; Widelitz et al., 1996). FGFs 1, 2 and 4 can induce many smaller feather buds from embryonic chicken skin explant cultures (Widelitz et al., 1996) and induce dermal condensations (Song et al., 2004). On the other hand, FGF10 induces epidermal thickening and enlarges feather primordia while decreasing the number of feather buds (Tao et al., 2002). The addition of dominant negative soluble FGFRs 1 and 2 at an early developmental stage inhibit feather bud formation (Mandler and Neubuser, 2004). FGFR1 is expressed initially beneath the feather placode and subsequently in the anterior feather bud mesenchyme. FGFR2 is expressed in the mesenchyme between the feather placode and the ectoderm of feather buds. On the other hand, FGFR3 expression is more ubiquitous (Noji et al., 1993). In scaleless mutant skin, beta-catenin is induced in the track field, but fails to form periodic patterns (Widelitz et al, 2000), while FGF can rescue feather bud formation in the skin of embryonic scaleless mutant chickens (Song et a., 1996; Viallet et al., 1998). These results imply an essential role for FGF in epithelial mesenchymal interaction during feather morphogenesis. However, how it works at the cellular level remains unknown.

FGFs exert their effects through the Raf-MEK-ERK pathway in a number of developmental systems (Cabernard and Affolter, 2005; Delfini et al., 2005; Fisher et al., 2001; Matsubayashi et al., 2004; Sawada et al., 2001). Here we show that FGF-ERK activity plays a role in controlling the migratory behavior of cells observed during the dermal condensation process, which helps to transform a homogeneous field of feather precursor cells into discrete cellular condensations. Adding inhibitors of ERK phosphorylation to pattern-forming explant cultures produces a spectrum of placode patterns ranging from broad stripes when added early in placode formation to anastomosing, segmented stripes and spots when added later. We propose that these stripes and spots represent a continuum of equilibrium states of morphing dermal condensations during the patterning process (Fig. 6A). To gain a global understanding, we develop a mathematical model with cell movement driven by chemotaxis towards areas containing high chemo-attractant concentrations. We show that chemotactic movement toward initial chemo-attractant concentration peaks is essential for stabilizing the bud pattern and that this process requires ERK signaling in the mesenchyme (Fig. 7; Supplemental movie 6Supplemental movie 8). The mathematical formulation of our model can reproduce the experimentally observed results and predict cell behavior. This study integrates theoretical and experimental approaches and opens new avenues of research in patterning stem cells.

Figure 6
Simulated model of periodic pattern formation for skin appendages

Results

ERK is specifically phosphorylated in feather buds during the periodic patterning process

In the spinal tract, feather formation begins with the development of a primary row of buds along the midline. The primary row initially forms as a stripe, as visualized by the molecular expression of several molecules, which contains an initiation point for feather formation. This serves as the point of convergence between stripes and spots (feather buds). After the primary row of feather buds form, a morphogenetic patterning wave propagates bilaterally to establish the lateral rows. FGFs and FGFRs are expressed in a restrictive mode (initially expressed all over, becoming restricted to the bud region; Widelitz et al., 2000) in the forming feather primordia (Jung et al., 1998; Song et al., 1996; Mandler and Neubuser, 2004; Tao et al., 2002; Widelitz et al., 1996). Here we focus on the intracellular signaling events downstream of the FGFs. To study the involvement of MAPK/ERK signaling, the phosphorylated form of ERK was examined in stages 28 and 35 skin using whole-mount immunostaining. Along the primary row, buds progressively form in a posterior to anterior direction. p-ERK is first homogeneously expressed in the morphogenetic zone (the periodic pattern forming region; Jiang et al., 1999). It then becomes restricted to the bud region, with a surrounding clear, lateral, halo-shaped inhibition zone, and finally it segregates into periodically arranged buds (Fig. 1C). A high power view shows that at the feather placode stage, there are medium levels of p-ERK present distributed in a diffuse way. At the short bud stage, medium levels of p-ERK are still present, but a central domain with high p-ERK expression starts to appear (Fig. 1C'). In the long bud stage, strong expression is restricted to the bud region, while the inter-bud region becomes deficient of p-ERK, establishing a distinct boundary between the two. This is typical of the restrictive mode of expression. Sections show p-ERK staining in the mesenchyme (Fig. 1D). Expression of Raf, an upstream ERK pathway member, was assessed by in situ hybridization. It also shows a restrictive mode expression pattern: first expressed throughout the morphogenetic field at stage 29, then becoming progressively restricted to the feather buds. Later, at stage 35, Raf expression becomes further restricted to the posterior feather bud (Fig. 1E). Raf is also expressed in both the epithelium and mesenchyme (Widelitz et al., 1996).

Suppressing p-ERK activity leads to a range of stripe patterns

To test the role of p-ERK, we first perturbed p-ERK pathways using siRNA targeted suppression of ERK transcripts. Skin explants were electroporated with siRNA-ERK or a scrambled siRNA (control) sequence. CMV-RFP was co-electroporated to enable identification of the affected regions. Results showed feather buds appear as short anastomosing stripes and merged plateaus (Fig. 2A).

Figure 2
Converting spots to stripes by suppression of ERK pathway

To manipulate the level of p-ERK more effectively, we used a chemical inhibitor of ERK phosphorylation. Western blot analysis showed that U0126 reduces p-ERK levels in skin explants (Fig. 2B). When U0126 was added to stage 31 skin explants, the feather patterns changed in a dose-dependent manner (Fig. 2C). At low dosage (2.5 µg/ml or 6 µM), feather bud elongation was inhibited. The diameter of the bud base expanded while the interbud spacing diminished. At a high dosage (10 µg/ml or 25 µM), the buds fused into stripes. At the highest dosage (20 µg/ml), feather bud formation was completely inhibited. The feather buds were all fused together to form broad stripe.

Inhibition of p-ERK in the presence of U0126 was verified by whole mount immunostaining for p-ERK using Rabbit anti-p-ERK antibodies. This showed that p-ERK was expressed in feather bud regions and was gradually reduced in explants treated with 2.5 and 10 µg/ml of U0126. In explants treated with 20 µg/ml of U0126 (or 50 µM), p-ERK expression was totally suppressed (Fig. 2D). Thus there is a U0126 dose-dependent reduction of ERK phosphorylation levels that correlates with the gradual morphological transformation from spots to stripes.

The feather explant culture system provides an excellent opportunity to observe the consequence of suppressing p-ERK at different time points. Since feather buds first form in the posterior midline, then anterior midline, and later the lateral regions of the explants (Fig. 6A. schematic bird drawing), different parts of stage 28 and 31 explants can actually represent different time periods from early to late development. Specifically, we consider that the lateral region of a stage 28 explant, mid-region of a stage 28 explant (Fig. 3C), lateral region of a stage 31 explant and mid-region of a stage 31 explant (Fig. 3B) represent different stages of feather development. If suppression takes place very early, as in the lateral region of a stage 28 explant, formation of the bud domain is entirely inhibited. The explant forms broad sheets which gradually break into long stripes towards the midline. Long stripes can be seen in the lateral region of a stage 31 explant. Towards the midline, these stripes further break into anastomosing stripes, and eventually become short stripe segments and some discrete buds. Thus, these stripes and spots represent a spectrum of patterns when the dermal condensation process is abolished at different times. Similar altered patterns can also be observed in the control (Fig. 4A and and5A).5A). Short “bud stripes” resemble those observed on skins treated with the DN FGFR (Mandler and Neubuser, 2004).

Figure 3
Molecular expression of the stripes caused by inhibition of ERK signaling
Figure 4
Involvement of p-ERK signaling in mesenchymal condensation
Figure 5
Involvement of FGF in p-ERK dependent chemotactic cell migration

To test the reversibility of the effect of U0126, we cultured skin explants in the presence of U0126 for two days. U0126 was then washed away, and the explant continued to be cultured for two or four more days (Fig. 3F). The temporal sequence of skin maturation from posterior midline, anterior midline to lateral flank regions was investigated (Fig. 3F, the red, green and blue regions, respectively). If the region has already formed feather primordia when it encounters U0126, the feather primordia cannot continue to attract dermal cells to form distinct dermal condensations, and therefore a broad area of stripes form (Fig. 3F, red region). If the region is in the middle of the periodic patterning process when it encounters U0126, many small buds form as the dermal condensation processes that are supposed to continue to build up a bigger buds are now blocked (Fig. 3F, green region). If the region has not started the periodic patterning process when U0126 was present, and started patterning after U0126 was removed, then the pattern is not affected (Fig. 3F, blue region). Thus, differences in temporal competence are revealed as different phenotypes in different explant regions. Furthermore, the phenotypes, once formed, are stable as seen in the explant after either 2 or 4 days of culture (Fig. 3F, right panel). Thus the effect of U0126 on the periodic patterning of feather buds is irreversible.

Molecular characterization of the stripes

What is the differentiation status of these stripes? Can they be considered as equivalent to an interbud region, are they a tract field which cannot progress, or do they represent giant, fused buds? We examined molecular expression patterns to help interpret these findings. We tested both restrictive and de novo mode signaling molecules, using L-fringe (L-fng) and Sonic hedgehog (Shh) as representatives, respectively. When stage 31 skin was cultured for one day, L-fng was present in the periodically arranged array, in either a circular or posterior configuration, representing its expression patterns at different developmental stages (Chen and Chuong, 2000). At day two, in the presence of U0126, the medial rows form buds which fuse laterally to form short “bud stripes”. The lateral rows are still in the early morphogenetic stage with unsegregated L-fng staining (Fig. 3A). Notch has similar expression patterns (not shown). Shh showed weak staining at day two. At day three, Shh appeared in the short horizontal bud stripes which form in the midline region, and in the two long longitudinal stripes which form at the lateral edges of the explants (Fig. 3B). An inter-stripe region is present between them. When stage 28 explants were used, a wider and more homogeneous stripe appeared with homogeneous Shh expression (Fig. 3C).

Tissue sections showed that the placode marker, Shh, is diffusely present at a low level in the stripe epithelia (Fig. 3D). β-catenin is expressed but also at lower levels than in control bud epithelium (not shown). The epidermis is thickened (Fig. S1D and Fig. 3D), but remains unpatterned. The dermal condensation marker neural cell adhesion molecule (NCAM) is expressed periodically at the base of individual buds (Jiang and Chung, 1992), but in a continuous and diffuse layer of dermal cells in the fused bud mesenchyme (Fig. 3D, H&E staining of these sections is shown in Fig. S1D).

Cell density in bud, interbud and striped regions was quantified. From homogeneity, cell density becomes mosaic-like: high (1.86 ± 0.13/100 µm2) in the bud and low (1.24 ± 0.04/100 µm2) in the interbud regions (Sengel., 1976; Fig. 3E). The striped region showed an average cell density of 1.47 ± 0.07/100 µm2 (Fig. 3E), suggestive of a plateau state somewhere between the bud and basal state (Fig. 7).

p-ERK activity is required for dermal condensation

E30 chicken skin was cultured with or without 25µM U0126 for 16 hours and the epithelium and mesenchyme then separated after 2x CMF +0.25% EDTA buffer incubation on ice for 10 min. The separated epithelium and mesenchyme were then recombined in different ways (Fig. 4A). As the effect of U0126 on bud morphogenesis is irreversible (Fig. 2F), we can now determine whether the action of U0126 is more dependent on epithelium or mesenchyme. U0126-treated or non-treated epithelia and mesenchyme were recombined, and the chimeric explants were cultured for two more days. U0126-treated epithelium recombined with non-treated mesenchyme resembled control feather arrays. In contrast, control epithelium recombined with U0126-treated mesenchyme caused the formation of short stripes, co-existing with discrete primordia. The phenotype was more severe when U0126-treated epithelium was recombined with U0126-treated mesenchyme.

To test the effects of FGF4-ERK signaling on cell migration, we used DiI to track cell motility (Fig. 4B). DiI was injected directly to the mesenchymal side of skin explants, labeling different mesenchymal cell regions. The original labeled spot is about 150–200 µm in diameter. The location of the DiI was determined at 0, 19, 43 and 67 hours (Fig. 4B). In a period of 67 hours, cells migrated 109.5 µm ± 11.5 µm in control skin and 231 µm ± 34.6 µm in U0126-treated skin (Figure 4D). Therefore, there is a tendency in the control skin for cell movement to be condensed toward each bud. When the MAPK pathway is suppressed, cells scatter over a much wider range, implying a failure of the dermal condensation process. An interesting observation is that cell movement in the treated skin tends to be along the anterior-posterior axis, consistent with the stripes formed near the midline.

Next we applied DiI to stage 30 embryonic skin, in both bud and interbud mesenchyme to determine whether migration differed between these two regions. We injected the mesenchymal regions of ten control and ten U0126-treated skin explants and cultured them for 68 hours. The location of DiI was photographed at 0 (green) and 68 hours (red). We then overlaid the red photograph over the green photograph to visualize how far the cells migrated over this period of time (Fig. 4B’). Mesenchymal cells did not migrate uniformly in all directions but in each example DiI labeled cells migrated further in U0126-treated skin. Both bud and interbud regions showed increased migration after U0126 treatment.

The degree of mesenchymal cell migration seen in B' implies unstable dermal condensation in early buds. As can be seen in the bottom panel (Fig. 4C): when some cells were labeled in the putative feather bud, they were able to migrate into the interbud region. We think this mixing is facilitated by inhibition of pERK. This cell mixing does not occur in later feather bud, as the boundary between bud and interbud become stabilized.

P-ERK activity is downstream to epithelial FGF4

In order to discover molecules that may be upstream to p-ERK, we tested several growth factors (FGFs and EGFs) which may activate receptor tyrosine kinase. For example, FGF4-coated beads were implanted on stage 31 skin explants. FGF4 was used because it is expressed in the bud and known to induce feather buds (Widelitz et al., 1996). Four hours later, a large zone (2 mm in diameter) of p-ERK was induced around the bead (Fig. 5A). In controls, normal p-ERK expression remains in each bud domain. When FGF10-coated (100 µg/ml) beads were placed on stage 31 skin explants, a zone expressing comparable levels of p-ERK, as seen in feather buds, is also induced (Fig. 5A). When FGF 4 and FGF 10 beads were implanted side by side, the p-ERK-expressing region was enlarged. However, regions under the influence of both FGFs do not appear to express more p-ERK, suggesting that the two FGFs share the same pathway. p-ERK expression is also reduced in the region immediately outside of this zone, presumably because the majority of dermal cells in this region, which normally would contribute to feather buds, were instead recruited to the bead source of high FGF/p-ERK activity. These data suggest that the FGF/Raf/p-ERK pathway is involved in the feather periodic patterning process.

We then examined the effects of an FGF4-coated bead on dermal condensation formation using visible light captured by time-lapse video-microscopy (5 mins per frame) for 10 hours (Fig. 5B, Movie 5). By 10 hours, a dermal condensation had formed around the FGF4 bead. The condensation continued to grow through 20 hours but remained at the same size through 30 hours (Fig. 5B and data not shown).

To explore the role of FGF4 as a chemoattractant, 1% of the dissociated mesenchymal cells were labeled with Vybrant CFDA SD (Invitrogen, Carlsbad, CA) and plated at a density of 1 × 105 cells/10 µl. Heparin beads were placed in 3 µl of an FGF4 solution (100 µg/ml) for 1 hour at 4°C. The FGF-coated beads then were placed on top of plated cells. The migration of cells towards the bead can be seen as the accumulation of fluorescent cells around the beads by 26 hours (Fig. 5C). Cells did not accumulate around control heparin beads. The graph summarizing 5 separate experiments shows that significantly more cells accumulated around the FGF4-coated beads (location 1 – red) than at distances progressively further from the bead (locations 2 – green and 3 – blue, respectively; Fig. 5C’).

p-ERK expression in the mesenchyme is epithelium-dependent. We tested whether FGF4-coated beads can substitute for the effect of the epithelium on denuded mesenchyme by removing part of the epithelium from a skin explant (Fig. 5D). An FGF4-coated bead placed on denuded mesenchyme was able to induce mesenchymal cells around the bead to form condensations. Thus, there is a partial rescue of condensation formation. In the presence of U0126, distinct condensations towards the bead are suppressed and the dense mesenchyme appears as anastomosing cords. These results suggest that FGF4 induces p-ERK activity in the mesenchyme, causing mesenchymal cells to migrate towards areas with high FGF4 concentration and accumulate there, forming distinct condensations.

Mathematical simulation of cell behavior in the chemotactic stage of periodic patterning

In earlier work, we discussed the role of a reaction-diffusion mechanism in establishing the feather pattern (Jung et al., 1998). Here we focus on the role of chemotaxis in expanding feather primordial signaling centers during early feather bud formation. To further understand our results described above we developed a simple mathematical model based on the chemotactic response of cells to FGF4 (Fig. 6). We hypothesize that, in line with experimental observations, spot or stripe patterns in cell density are formed from an initially uniform field of mesenchymal cells by directed cell migration. We assume that cell motility becomes directed in such a way that cells move up gradients in chemo-attractant (FGF4) concentration i.e. chemotaxis. We further assume that cells produce the chemo-attractant, which itself undergoes random motion and is removed via decay or by endogenous antagonists or inhibitors of the relevant signaling pathway. The model can be summarized using the following “word” equations. More details of the mathematical model are provided in the Supplementary Information.

Chemical equation

Rate of change of chemical concentration = random diffusion + chemical production − chemical decay/removal.

Cell equation

Rate of change of cell density = random motion + directed motion (chemotaxis).

Intuitively, the model works as follows: small random fluctuations in chemo-attractant concentration or cell density are amplified by a feedback loop of chemo-attractant production and cell movement. For example, a small peak in chemo-attractant concentration causes cells to move preferentially in the direction of the peak, where chemo-attractant production increases (due to increased cell density). “Positive feedback” competition between neighboring peaks results in a pattern of cell density. Some peaks are eliminated while others are later stabilized to become dermal condensations. We define three major states as a competent basal state, a bud state, and an interbud state (Fig. 6B). The competent state is homogeneous. “Spot patterns” represent the conversion from a competent state to an oscillating bud and interbud state. “Stripe patterns” represent the one-dimensional waning of the oscillations to reach a plateaued state that lies somewhere between the competent and interbud states.

Using numerical tools we simulated the patterns in cell densities that arise via the model. Fig. 6C shows a control case: the resulting pattern along the primary row is spots of high cell density, which represent the bud state, interspersed with regions of low cell density, which represent the interbud state. The effect of perturbing p-ERK activity, as described earlier, can be tested using this model by varying the rate at which the chemo-attractant is removed. Fig. 6D shows the result of increasing the rate of chemo-attractant removal: oscillations in the medio-lateral direction gradually disappear, leading to horizontal stripes with cell density lying in a state between the basal and bud states. Fig. 6E shows the result as the removal rate is increased even further: the feather field remains in the homogeneous state. Computer simulations of the processes are shown in the Supplementary Information as movie 6Supplementary Information as movie 8. The results of the mathematical model are consistent with those observed experimentally and support the hypothesis that p-ERK-dependent cell chemotaxis is involved in the periodic patterning of feather buds.

Discussion

In this work we use an experimentally manipulatable patterning model in which cells evolve from a homogeneous state (basal state) into become periodically array of elements arranged with different size, shape and spacing (Fig. 1A). We found that p-ERK-dependent chemotaxis is essential for the successful formation of feather patterns, suggesting this process may work alongside or in sequence with a reaction-diffusion mechanism (Fig. 7B). Another unique feature of this model is that there are two components, epithelia and mesenchyme, in this patterning process. We found that the successful patterning process requires cooperation of epithelium and mesenchyme, with molecular signals coming from one component or the other, thus demonstrating the essence of tissue interactions (Fig. 7A). An earlier study blocked FGF activity completely by delivering an FGF dominant receptor to early developing skin that suppressed bud formation (Mandler and Neubuser, 2004). While it demonstrates the essentialness of FGF activity, it is informative to study the time course of pattern formation as the whole process is inhibited. Here we adjusted the level of p-ERK activity during later stages of feather morphogenesis and by doing so, revealed the pleomorphic patterns generated at different developmental times when the dermal condensation process is blocked. Analyses of these results with time lapse movies and mathematical modeling provide new insights for understanding the process of feather pattern formation.

FGF4 / ERK activity regulates the chemotactic process of mesenchymal condensations

Previously, we and others explored the role of FGFs as activators of feather formation (Jung et al., 1998; Song et al., 1996, 2004; Widelitz et al., 1996). We proposed that FGF works in conjunction with inhibitors (ie, BMPs; Jung et al., 1998; Noramly and Morgan, 1998) through a reaction-diffusion mechanism to establish feather pattern formation (Jung et al., 1998). Here we further explore the role of FGFs functioning as a chemo-attractant at a later feather morphogenesis phase. The FGF pathway belongs to the receptor tyrosine kinase (RTK) signaling family. RTK signaling can activate different downstream effectors including the MAPK/ERK pathway, the phosphatidylinositol 3 kinase (PI3K) pathway, or the phospholipase Cr (PLCr) signaling pathway. The phospholipase C inhibitor, U73122, has no effect on feather pattern formation (our unpublished data). Inhibition of PI3K affects feather bud outgrowth but did not affect the expansion of feather primordia (Atit et al., 2003). Here we focus on the MAPK/ERK signaling pathway. We show that inhibiting this pathway at different stages produces stripes of intermediate cell density, rather than clear-cut bud or interbud fates. Analyses of these stripes showed that feather placode markers, such as Shh (Ting-Berreth and Chuong, 1996a), remain diffusely distributed. Bud growth is retarded and buds fail to express differentiation markers. Surprisingly, EGF, another growth factor which functions through a tyrosine kinase receptor, does not invoke this pathway in feathers (Attit et al., 2003).

Cell tracing and time lapse movies show that in early feather morphogenesis, the ERK pathway affects the process of mesenchymal dermal condensation. Further support comes from the diffuse presence of the dermal condensation marker NCAM (Jiang and Chuong, 1992) in ERK-suppressed specimens. This parallels the involvement of growth factors in regulating adhesion molecules during limb morphogenesis (Hentschel et al., 2004; Kiskowski et al., 2004). Using a localized bead source, we showed that FGF4 can specifically induce phosphorylation of ERK. In the absence of epithelium, mesenchymal cells move randomly and cannot form discrete, spot-like dermal condensations. FGF4-coated beads placed on denuded mesenchyme rescued dermal condensation formation around the FGF source. This rescue can be neutralized by an inhibitor of ERK phosphorylation. Epithelial-mesenchymal recombination showed that the ERK inhibitor effect is mediated mainly through the mesenchyme. Consistent with these findings, FGF4 is produced by the epithelium and FGFR1 is present in the feather bud mesenchyme (Jung et al., 1998; Noji et al., 1993). These data support the notion that the FGF/FGFR/RAF/ERK pathway is involved in mediating epithelial-mesenchymal interactions with FGF from the epithelial placode binding FGFR in the mesenchyme leading to downstream events. This does not rule out an effect of FGF activity on the epithelium, which then affect mesenchyme indirectly. Previously we showed that TGFβ 2 and Shh mediate the epithelial effect on dermal condensation formation (Ting-Berreth and Chuong, 1996a, b). The relationship between TGFβ 2, Shh and FGF4 in the induction of dermal condensation remains to be investigated.

We believe that short-range activation, long-range inhibition is involved in forming micro-aggregates during the early stages of feather patterning, and many different classes of molecules may be involved in this process as activators or inhibitors (growth factors and their antagonists, signaling molecules, adhesion molecules). However, what matters to the patterning is the summed activity of activators or inhibitors. In line with this view, recently, based on the counteracting roles of BMP7 and BMP2 in chicken skin, a model involving cell migration and cell adhesion was proposed as underlying dermal condensation formation during feather morphogenesis (Michon et al., 2008). Whereas BMP7 induced cell migration in these studies, BMP2 induced integrin α4 and altered the splicing of fibronectin to exclude the fibronectin EIIIA domain, decrease migration and foster intracellular adhesion within the forming placode. Based on these observations, Michon et al. presented a simulation of the process of feather bud formation. Whilst their model can explain a number of experimentally observed phenomena, such as the appearance of a “zone of inhibition” around an FGF4-coated bead, it cannot explain the feather bud fusions we observe in this work.

Pleomorphic patterning of stem cells

It appears that a moderate basal level of β-catenin marks the basal state of epithelia in the morphogenetic field, with competence to become either bud or interbud. Through the patterning process, this homogeneous field breaks into regions of high and low β-catenin expression levels. If a localized region with high level of β-catenin is achieved, then feather buds form (state A; Noramly et al., 1999; Widelitz et al., 2000). If a region becomes deficient in β-catenin, it assumes the inter-bud fate (state B). Do the stripes observed here represent the basal, bud or interbud states? The stripe epithelium can express the placode marker, Shh, so it is not in the interbud state. On the other hand, Shh expressing cells are diffusely distributed all over the stripes, and the stripe cannot grow in height (i.e., it fails to form a localized growth zone), nor express feather keratin (i.e., it fails to differentiate). Therefore, the stripe cannot be considered as in the bud state either. Mesenchyme under the stripe expresses NCAM diffusely. Proliferation normally occurs in periodically arranged localized growth zones in bud regions, but the stripe exhibits un-patterned proliferation. By tracing the dynamics of DiI labeled mesenchyme cells, cell motility was random and active at first. Gradually, cells move towards dermal condensations where they become more adhesive and stationery. In U0126-treated specimens, mesenchymal cell motility remains high and discrete dermal condensations fail to form. Cell density in the stripes is higher than that of the interbud regions, but lower than that of the bud regions. Thus the stripe we study here has a property between the original basal state and the bud state (Fig. 6A).

We propose that an organ field in the embryo represents a multi-potential status, competent to form different cell types in that organ. Through differentiation and patterning, these precursor cells are specified simultaneously or sequentially to build the topological order of the tissue. If the patterning process is altered, differentiated cells can be distributed in different configurations, in spots or stripes, consisting of a spectrum of differentiation states (Fig. 1A). Thus cells with a homogeneous stem cell status can give rise to pleomorphic patterns, depending on local signaling cues. Recently, it was shown that Wnt and DKK can also work as Turing activator- inhibitor pairs during the formation of hair germs (Sick et al., 2006; Maini et al., 2006). In accordance with the pleomorphic patterning concept here, forced expression of β-catenin resulted in a failure of the periodic patterning process, leading to the randomization of placode size and spacing, forming giant placode-like structures. These giant “placodes” exhibit differentiation markers intermediate to epidermis and hair follicle keratinocytes (Närhi et al., 2008; Zhang et al., 2008). Such results are consistent with the pleomorphic patterning concept here.

A generic cellular behavior model for a two-component periodic patterning process

Combining knowledge from our current and previous results (Jiang et al., 1999) and using FGF/p-ERK signaling as an example, we envision the following epithelial-mesenchymal interactions taking place during the patterning process (Fig. 6A, Fig. 7A). (a) Competence stage. Epithelia and mesenchyme gain competence when a feather tract field forms. (b) Reaction-diffusion stage. At this competent state, cells exhibit basal adhesiveness and random motility. Reactions occur and unstable microaggregates form. These microaggregates are only capable of eliciting a weak signal. (b') Dermal-epidermal signaling. Since FGFs are produced in the epithelium, and FGFR1 is present in the mesenchyme, it is possible that microaggregates may induce the epithelia above to express more FGF, helping to create a larger signaling center that would elicit a greater signal. The molecular nature of this signal is still undergoing investigation. p-ERK expression is homogenous (Figure 1C, flank regions in the left panels). (c) Dermal condensation stage. Mesenchymal p-ERK is induced in the mesenchyme beneath the epithelia with higher FGF, and its activity induces chemotaxis of mesenchymal cells, initiating the uneven distribution of cell density within the dermis. p-ERK positive primordia start to emerge at the junction of homogenously stained and distinct bud region (Figure 1C, left panels). (d) Putative primordia stage. The increasing size of dermal condensations further induces the overlying epithelia to express higher levels of FGF, which induces more p-ERK in the mesenchyme, and bigger dermal condensations. We called this “putative” because we know the placodes at this stage are unstable (Fig. 4D) and can be reset (Chuong et al., 1996; Jiang et al., 1999). p-ERK staining shows moderate levels in the putative primordia with a blurred margin. The center regions show higher expression levels that are expanding (Fig. 1C', left and central panels). (E) Stabilized bud stage. The positive feedback loops continue until each dermal condensation / placode is stabilized. Thus the fate of feather bud and interbud states are established in both epithelia and mesenchyme. p-ERK is now strongly expressed in the bud but excluded from the interbud region, forming a sharp boundary (Fig. 1C', right panel). Our model focuses on the mesenchyme. We think the mesenchymal message is received and interpreted by the overlying epithelium which causes epithelial cell rearrangements. Although we have focused on the role of ERK in patterning the mesenchyme, we cannot rule out a role for intraepithelial signaling which is the subject of current investigations ongoing in our laboratory.

Based on these results, we develop a generic model that will not only be applicable to skin appendage formation, but also can be applied to understand how stem cells are patterned in general. We propose a mathematical model to aid in understanding the periodic patterning process. In early stages, the short-range activation, long-range inhibition is the major mechanism patterning the mesenchyme. This process is epithelial-independent, i.e. it can occur in the absence of epithelia (Jung et al., 1998; Jiang et al., 1999). During later stages, chemotactic cell migration towards micro signaling centers is converted into larger stable dermal condensations. During this process, some original microaggregates will merge to be bigger or disappear, due to the phenomenon of peak competition inherent in chemotaxis. These two stages transit gradually without a sharp time boundary and will largely overlap in time (Fig. 7B) although further evidence would have to be demonstrated through future video time lapse cinematography experiments. The process is epithelial-dependent and here we show that this epithelial-mesenchymal interactive patterning process involves the expression of FGF by the epithelia and the phosphorylation of ERK in the mesenchyma.

Besides identifying the molecular members essential for feather formation, we aim to understand how the number, size and arrangement of feather primordia are determined. To this end, we developed a mathematical model based on the ability of cells to migrate toward chemo-attractants which promote bud formation as a way of testing our experimental hypothesis. Our model was able to predict the characteristic “spotted” patterns of cell density observed in control embryos and, further, to predict the effects of p-ERK perturbation. It can also be used to predict the effects of other possible chemo-attractants and their antagonists. Increasing the rate of chemo-attractant removal produced “striped” patterns with cell density lying between the bud and interbud states.

We acknowledge that there are other mechanisms that can lead to the formation of patterning in biology. The formation of somites (spot) from presomitic mesoderm (stripe) involves oscillating expression levels of FGF, p-ERK, HES1 and Wnt (Delfini et al., 2005). In limb bud micromass cultures, local activation coupled with lateral inhibition is proposed to be important for the formation of cartilage nodules (Newman and Bhat, 2007). It is tempting to propose some unifying fundamental mechanism may underlie these periodic patterning processes. More studies that combine experimental and theoretical approaches will be required for us to further understand the essence of biological pattern formation (Chuong and Richardson, in press).

Experimental Procedures

Embryos, skin explant culture and recombination

Specific pathogen-free (SPF) White Leghorn chick embryos were obtained from Charles River Laboratories (Preston, CT, USA) and staged according to Hamburger and Hamilton (H&H) (1951). Explant cultures, partial removal of epithelia and bead preparations were performed as described (Ting-Berreth and Chuong, 1996b). Application of the inhibitor, U0126 (Sigma, St. Louis, MO, USA), to chicken skin explant cultures was achieved by supplementation of the growth media. U0126 was dissolved in DMSO and added to a final concentration of 10 µg/ml. This concentration was optimized in preliminary studies within the laboratory. Control cultures were supplemented with an equivalent volume of DMSO.

In situ hybridization and immunohistochemistry

Procedures were performed as described (Jiang and Chuong, 1992; Ting-Berreth et al., 1996a). The automated Discovery in situ hybridization unit from Ventana was used for some specimens.

Electroporation, probe and siRNA preparation

The ERK siRNA sense sequence is AGA UCU UAC UGC GCU UCA GTT and antisense sequence is CUG AAG CGC AGU AAG AUC UTT. The target oligonucleotide was suspended in DEPC-water at a concentration of 50 µM. The ERK siRNA oligonucleotides were annealed with 5X annealing buffer (250 mM Tris-HCl, pH 7.5, 500 mM NaCl, DEPC-water) for 2 minutes at 95°C and allowed to cool to room temperature. siRNA-ERK was diluted in electroporation hypo-osmolar buffer (eppendorf) 5µM and RFP added as a marker to visualize affected sites by electroporation. E31 chicken skin was electroporated with an siRNA –ERK mixture in which the negative pole faces the epithelium. The electric current was delivered as 3 pulse of 5v/50ms. After electroporation, the skin was cultured on 60mm culture inserts (Falcon, San Jose, CA, USA) for 2 days. The electroporation efficiency was checked by florescence microscopy. The ERK-siRNA explants have reduced feather bud elongation. The feather buds were partially inhibited by ERK-siRNA.

Cell labeling

To observe cell migration towards an FGF4 coated bead, 1% of the dissociated mesenchymal cells isolated from stage 31 skin explants were labeled with Vybrant CFDA SE Cell Tracer (10 mM; Invitrogen, Carlsbad, CA, USA) for 15 mins at 37°C and washed with media for 30 mins at 37°C following the manufacturer’s recommended protocol. The labeled cells were then mixed with unlabeled cells and 1 × 105 cells were plated in 10 µl on a culture insert (Falcon, San Jose, CA, USA). For lineage tracing cells in skin explants were microinjected with 1–5 nl of 0.1% DiI in 70% ethanol (Invitrogen, Carlsbad, CA, USA) via glass capillaries followed by two washes with PBS. The labeled skin explants were cultured for 72 hours. The results were visualized by fluorescent microscopy.

Time lapse video microscopy

Skin explants were grown on 60 mm culture inserts (Falcon, San Jose, CA, USA) in a SmartSlide environmental chamber (Wafergen Biosystems, Fremont, CA, USA) which maintained temperature and CO2 conditions. Time lapse photographs were obtained every 5 minutes through an Olympus IMT-2 microscope.

Quantification of cell density

Mesenchymal cell density was calculated in feather buds and interbud regions in control specimens. The cells found within three feather buds and interbuds were counted and cell density calculated by dividing the total size of the bud and interbud regions (µm2). The mesenchymal cells from three fused skin samples were counted and cell density calculated by dividing by the size of the fused area (µm2).

Supplementary Material

02

03

04

05

06

07

09

10

Acknowledgements

We are grateful to grant support from NIAMS AR2177, AR47364 (CMC), and AR052397 (RW). REB would like to thank Microsoft Research for a European Postdoctoral Research Fellowship, Research Councils UK for an RCUK Academic Fellowship in Mathematical Biology and St Hugh's College, Oxford for a Junior Research Fellowship. REB also gratefully acknowledges support from the Astor Travel Fund at the University of Oxford in order to visit the University of Southern California. PKM was partially supported by a Royal Society Wolfson Merit Award. Some photos use microscopes from the microscopy sub core at the USC Center for liver diseases (NIH 1 P30 DK48522).

Footnotes

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Disclosure

Widelitz serves on the advisory board to Wafergen Biosystems.

Accession Numbers

Raf-1 accession number NM205307.

References

  • Atit R, Conlon RA, Niswander L. EGF signaling patterns the feather array by promoting the interbud fate. Dev. Cell. 2003;4:231–240. [PubMed]
  • Ball P. The self-made tapestry: pattern formation in nature. New York: Oxford University Press; 1999.
  • Cabernard C, Affolter M. Distinct roles for two receptor tyrosine kinases in epithelial branching morphogenesis in Drosophila. Dev. Cell. 2005;9:831–842. [PubMed]
  • Chang CH, Jiang TX, Lin CM, Burrus LW, Chuong CM, Widelitz R. Distinct Wnt members regulate the hierarchical morphogenesis of skin regions (spinal tract) and individual feathers. Mech. Dev. 2004;121:157–171. [PMC free article] [PubMed]
  • Chen CW, Chuong CM. Dynamic expression of lunatic fringe during feather morphogenesis: a switch from medial-lateral to anterior-posterior asymmetry. Mech. Dev. 2000;91:351–354. [PubMed]
  • Christley S, Alber MS, Newman SA. Patterns of mesenchymal condensation in a multiscale, discrete stochastic model. PLoS Comput. Biol. 2007;3:e76. [PMC free article] [PubMed]
  • Chuong CM, Richardson MR. Pattern formation. Int. J. Dev. Biol. 2009 In press. [PubMed]
  • Chuong C-M, Ting-Berreth S, Widelitz RB, Jiang T-X. Early events during the regeneration of skin appendages: order of molecular reappearance following epithelial-mesenchymal recombination with rotation. J. Invest. Dermatol. 1996;107:639–646. [PubMed]
  • Chuong C-M, Wu P, Plikus MV, Jiang TX, Widelitz RB. Engineering Stem cells into organs: Topobiological transformations demonstrated by beak, feather and other ectodermal organ morphogenesis. Curr. Top. Dev. Biol. 2006;72:237–274. [PMC free article] [PubMed]
  • Crowe R, Henrique D, Ish-Horowicz D, Niswander L. A new role for Notch and Delta in cell fate decisions: patterning the feather array. Development. 1998;125:767–775. [PubMed]
  • Delfini MC, Dubrulle J, Malapert P, Chal J, Pourquie O. Control of the segmentation process by graded MAPK/ERK activation in the chick embryo. Proc. Natl. Acad. Sci. USA. 2005;102:11343–11348. [PubMed]
  • Dequéant ML, Pourquié O. Segmental patterning of the vertebrate embryonic axis. Nat. Rev. Genet. 2008;9:370–382. [PubMed]
  • Drew CF, Lin CM, Jiang TX, Blunt G, Mou C, Chuong CM, Headon DJ. The Edar subfamily in feather placode formation. Dev. Biol. 2007;305:232–245. [PMC free article] [PubMed]
  • Fisher CE, Michael L, Barnett MW, Davies JA. Erk MAP kinase regulates branching morphogenesis in the developing mouse kidney. Development. 2001;128:4329–4338. [PubMed]
  • Hamburger V, Hamilton HL. A series of normal stages in the development of the chick embryo. J. Morphol. 1951;88:49–92. [PubMed]
  • Harris MP, Williamson S, Fallon JF, Meinhardt H, Prum RO. Molecular evidence for an activator-inhibitor mechanism in development of embryonic feather branching. Proc. Natl. Acad. Sci. USA. 2005;102:11734–11739. [PubMed]
  • Hentschel HG, Glimm T, Glazier JA, Newman SA. Dynamical mechanisms for skeletal pattern formation in the vertebrate limb. Proc. Biol. Sci. 2004;271:1713–1722. [PMC free article] [PubMed]
  • Jiang TX, Chuong CM. Mechanism of skin morphogenesis. I. Analyses with antibodies to adhesion molecules tenascin, N-CAM, and integrin. Dev. Biol. 1992;150:82–98. [PubMed]
  • Jiang T-X, Jung H-S, Widelitz RB, Chuong C-M. Self organization is the initial event in periodic feather patterning: roles of signaling molecules and adhesion molecules. Development. 1999;126:4997–5009. [PubMed]
  • Jiang TX, Widelitz RB, Shen WM, Will P, Wu DY, Lin CM, Jung HS, Chuong CM. Integument pattern formation involves genetic and epigenetic controls: feather arrays simulated by digital hormone models. Int. J. Dev. Biol. 2004;48:117–135. [PMC free article] [PubMed]
  • Jung H-S, Francis-West PH, Widelitz RB, Jiang T-X, Ting-Berreth S, Tickle C, Wolpert L, Choung C-M. Local inhibitory action of BMPs and their relationships with activators in feather formation: implications for periodic patterning. Dev. Biol. 1998;196:11–23. [PubMed]
  • Keller G. Embryonic stem cell differentiation: emergence of a new era in biology and medicine. Genes Dev. 2005;19:1129–1155. [PubMed]
  • Kiskowski MA, Alber MS, Thomas GL, Glazier JA, Bronstein NB, Pu J, Newman SA. Interplay between activator-inhibitor coupling and cell-matrix adhesion in a cellular automaton model for chondrogenic patterning. Dev. Biol. 2004;271:372–387. [PubMed]
  • Kondo S. The reaction-diffusion system: a mechanism for autonomous pattern formation in the animal skin. Genes Cells. 2002;7:535–541. [PubMed]
  • Lin CM, Jiang TX, Widelitz RB, Chuong CM. Molecular signaling in feather morphogenesis. Curr. Opin. Cell. Biol. 2006;18:730–741. [PMC free article] [PubMed]
  • Liu RT, Liaw SS, Maini PK. Two-stage Turing model for generating pigment patterns on the leopard and the jaguar. Phys. Rev. E. Stat. Nonlin. Soft Matter Phys. 2006;74:011914-1–011914-8. [PubMed]
  • Maini PK, Baker RE, Chuong CM. The Turing model comes of molecular age. Science. 2006;314:1397–1398. [PMC free article] [PubMed]
  • Mandler M, Neubuser A. FGF signaling is required for initiation of feather placode development. Development. 2004;131:3333–3343. [PubMed]
  • Matsubayashi Y, Ebisuya M, Honjoh S, Nishida E. ERK activation propagates in epithelial cell sheets and regulates their migration during wound healing. Curr. Biol. 2004;14:731–735. [PubMed]
  • Meinhardt H, Gierer A. Pattern formation by local self-activation and lateral inhibition. Bioessays. 2000;22:753–760. [PubMed]
  • Michon F, Forest L, Collomb E, Demongeot J, Dhouailly D. BMP2 and BMP7 play antagonistic roles in feather induction. Development. 2008;135:2797–2805. [PMC free article] [PubMed]
  • Miura T, Maini PK. Speed of pattern appearance in reaction-diffusion models: implications in the pattern formation of limb bud mesenchyme cells. Bull. Math. Biol. 2004;66:627–649. [PubMed]
  • Murray JD. Mathematical Biology II: Spatial Models and Biomedical Applications. New York: Springer-Verlag; 2003.
  • Myerscough MR, Murray JD. Pigmentation pattern formation on snakes. J. Theor. Biol. 1991;149:339–360. [PubMed]
  • Nagorcka BN, Mooney JR. From stripes to spots: prepatterns which can be produced in the skin by a reaction–diffusion system. IMA J. Math. Appl. Med. Biol. 1992;9:249–267. [PubMed]
  • Närhi K, Järvinen E, Birchmeier W, Taketo MM, Mikkola ML, Thesleff I. Sustained epithelial β-catenin activity induces precocious hair development but disrupts hair follicle down-growth and hair shaft formation. Development. 2008;135:1019–1028. [PubMed]
  • Newman SA, Bhat R. Activator-inhibitor dynamics of vertebrate limb pattern formation. Birth Defects Res. C Embryo Today. 2007;81:305–319. [PubMed]
  • Newman SA, Comper WD. 'Generic' physical mechanisms of morphogenesis and pattern formation. Development. 1990;110:1–18. [PubMed]
  • Newman SA, Christley S, Glimm T, Hentschel HG, Kazmierczak B, Zhang YT, Zhu J, Alber M. Multiscale models for vertebrate limb development. Curr. Top. Dev. Biol. 2008;81:311–340. [PubMed]
  • Noji S, Koyama E, Myokai F, Nohno T, Ohuchi H, Nishikawa K, Taniguchi S. Differential expression of three chick FGF receptor genes, FGFR1, FGFR2 and FGFR3, in limb and feather development. Prog. Clin. Biol. Res. 1993;383B:645–654. [PubMed]
  • Noramly S, Morgan BA. BMPs mediate lateral inhibition at successive stages in feather tract development. Development. 1998;125:3775–3787. [PubMed]
  • Noramly S, Freeman A, Morgan BA. β-catenin signaling can initiate feather bud development. Development. 1999;126:3509–3521. [PubMed]
  • Painter KJ, Maini PK, Othmer HG. Stripe formation in juvenile Pomacanthus explained by a generalized Turing mechanism with chemotaxis. Proc. Natl. Acad. Sci., USA. 1999;96:5549–5554. [PubMed]
  • Patel K, Makarenkova H, Jung H-S. The role of long range, local and direct signalling molecules during chick feather bud development involving the BMP’s, follistatin and the Eph receptor tyrosine kinase Eph-A4. Mech. Dev. 1999;86:51–62. [PubMed]
  • Pourquie O. Vertebrate somitogenesis: a novel paradigm for animal segmentation? Int. J. Dev. Biol. 2003;47:597–603. [PubMed]
  • Sawada A, Shinya M, Jiang YJ, Kawakami A, Kuroiwa A, Takeda H. Fgf/MAPK signalling is a crucial positional cue in somite boundary formation. Development. 2001;128:4873–4880. [PubMed]
  • Sengel P. Morphogenesis of skin. In: Abercrombie M, Newth DR, Torrey JG, editors. Developmental and Cell Biology Series. Cambridge: Cambridge Univ. Press; 1976.
  • Sick S, Reinker S, Timmer J, Schlake T. WNT and DKK determine hair follicle spacing through a reaction-diffusion mechanism. Science. 2006;314:1447–1450. [PubMed]
  • Small S, Levine M. The initiation of pair-rule stripes in the Drosophila blastoderm. Curr. Opin. Genet. Dev. 1991;1:255–260. [PubMed]
  • Song H, Wang Y, Goetinck PF. Fibroblast growth factor 2 can replace ectodermal signaling for feather development. Proc. Natl. Acad. Sci. USA. 1996;93:10246–10249. [PubMed]
  • Song HK, Lee SH, Goetinck PF. FGF-2 signaling is sufficient to induce dermal condensations during feather development. Dev. Dyn. 2004;231:741–749. [PubMed]
  • Tao H, Yoshimoto Y, Yoshioka H, Nohno T, Noji S, Ohuchi H. FGF10 is a mesenchymally derived stimulator for epidermal development in the chick embryonic skin. Mech. Dev. 2002;116:39–49. [PubMed]
  • Ting-Berreth SA, Chuong C-M. Sonic Hedgehog in feather morphogenesis: induction of mesenchymal condensation and association with cell death. Dev. Dyn. 1996a;207:157–170. [PubMed]
  • Ting-Berreth SA, Chuong C-M. Local delivery of TGF β2 can substitute for placode epithelium to induce mesenchymal condensation during skin appendage morphogenesis. Dev. Biol. 1996b;179:347–359. [PubMed]
  • Vasiev BN, Hogeweg P, Panfilov AV. Simulation of Dictyostelium discoideum aggregation via reaction-diffusion model. Phys. Rev. Lett. 1994;73:3173–3176. [PubMed]
  • Viallet JP, Prin F, Olivera-Martinez I, Hirsinger E, Pourquie O, Dhouailly D. Chick Delta-1 gene expression and the formation of the feather primordia. Mech. Dev. 1998;72:159–168. [PubMed]
  • Watt FM. Stem cell fate and patterning in mammalian epidermis. Curr. Opin. Genet. Dev. 2001;11:410–417. [PubMed]
  • Weissman IL. Translating stem and progenitor cell biology to the clinic: barriers and opportunities. Science. 2000;287:1442–1446. [PubMed]
  • Widelitz RB, Jiang T-X, Noveen A, Chen C-W, Chuong C-M. FGF induces new feather buds from developing avian skin. J. Invest. Dermatol. 1996;107:797–803. [PubMed]
  • Widelitz RB, Baker RE, Plikus M, Lin CM, Maini PK, Paus R, Chuong CM. Distinct mechanisms underlie pattern formation in the skin and skin appendages. Birth Defects Res. C Embryo Today. 2006;15:280–291. [PMC free article] [PubMed]
  • Wolpert L. Positional information and pattern formation. Curr. Top. Dev. Biol. 1971;6:183–224. [PubMed]
  • Yu M, Wu P, Widelitz RB, Chuong CM. The morphogenesis of feathers. Nature. 2002;420:308–312. [PMC free article] [PubMed]
  • Zhang Y, Andl T, Yang SH, Teta M, Liu F, Seykora JT, Tobias JW, Piccolo S, Schmidt-Ullrich R, Nagy A, Taketo MM, Dlugosz AA, Millar SE. Activation of β-catenin signaling programs embryonic epidermis to hair follicle fate. Development. 2008;135:2161–2172. [PMC free article] [PubMed]
  • Zorzano M-P, Hochberg D, Cuevas M-T, Gómez-Gómez J-M. Reaction-diffusion model for pattern formation in E. coli swarming colonies with slime. Phys. Rev. 2005 E 71, 031908. [PubMed]