A mechanism for the sharp transition of morphogen gradient interpretation in Xenopus
BMC Developmental Biology volume 7, Article number: 47 (2007)
One way in which positional information is established during embryonic development is through the graded distribution of diffusible morphogens. Unfortunately, little is known about how cells interpret different concentrations of morphogen to activate different genes or how thresholds are generated in a morphogen gradient.
Here we show that the concentration-dependent induction of the T-box transcription factor Brachyury (Xbra) and the homeobox-containing gene Goosecoid (Gsc) by activin in Xenopus can be explained by the dynamics of a simple network consisting of three elements with a mutual negative feedback motif that can function to convert a graded signal (activin) into a binary output (Xbra on and Gsc off, or vice versa). Importantly, such a system can display sharp thresholds. Consistent with the predictions of our model, Xenopus ectodermal cells display a binary response at the single cell level after treatment with activin.
This kind of simple network with mutual negative feedback might provide a general mechanism for selective gene activation in response to different levels of a single external signal. It provides a mechanism by which a sharp boundary might be created between domains of different cell types in response to a morphogen gradient.
One way in which positional information might be established during embryonic development is through the graded distribution of diffusible morphogens, including members of the TGF-β, FGF and Wnt families of growth factors [1–3]. Although progress is being made in coming to understand the ways in which morphogens can traverse fields of cells [4–6], rather little is known about how cells interpret different concentrations of morphogen to activate different genes or how thresholds are generated in a morphogen gradient. Recent studies indicate that morphogens frequently exert their effects through the post-translational activation of a single transcription factor, which in turn induces downstream target genes in a concentration-dependent manner. For example, dorso-ventral patterning in Drosophila embryo is controlled by the graded activation of the NF-kB-like transcription factor Dorsal .
In Xenopus, members of the TGF-β family such as activin and the nodal-related proteins act as morphogens and are essential for mesoderm formation [2, 8]. They activate downstream gene expression in a concentration-dependent manner, with low concentrations activating the T-box transcription factor Brachyury (Xbra) and high concentrations inducing the homeobox-containing gene Goosecoid (Gsc) . Changes in the extracellular activin concentration are reflected by differences in the concentration of nuclear Smad2, the effector of activin signalling . But how do differences in effector concentration cause the activation of different target genes? Several mechanisms to explain this phenomenon have been proposed , but it is not clear how they generate sharp thresholds, with small differences in morphogen concentration yielding qualitative differences in gene expression.
In this paper we show that the concentration-dependent induction of Xbra and Gsc by activin can be explained by the dynamics of a mutual negative feedback motif that can be adapted to function as a module to convert a graded signal (activin) into a binary output (Xbra on and Gsc off, or vice versa). We note that this system can display sharp thresholds, and it provides a plausible mechanism by which cells might distinguish between small concentration differences in a morphogen gradient and thereby create a boundary between two different cell types. Consistent with the predictions of our model, Xenopus ectodermal cells display a binary response at the single cell level after treatment with activin. We suggest that this mutual negative feedback gene network represents a general mechanism for selective gene activation in response to different levels of a single external signal.
A mutual negative feedback motif can explain the concentration-dependent induction of Xbra and Gsc
During Xenopus development, activin and the nodal-related proteins activate gene expression in a concentration-dependent manner, with low concentrations of activin inducing the expression of Xbra while high concentrations activate Gsc in isolated Xenopus animal pole regions. These two genes repress each other's expression, thereby creating a regulatory network based on mutual negative feedback: Gsc can repress transcription of Xbra by binding directly to its promoter [11, 12], while Xbra suppresses Gsc by inducing the expression of Xom, a repressor of Gsc (also known as Xvent2, Vox and Tbr-1), [13–17] (Fig. 1A, left panel).
We have attempted to predict the concentration-dependent effects of activin, and the dynamic aspects of its inductive activity, from this network structure. Our model consists of three interacting elements, M, A and B (Fig. 1A, right panel). M induces both A and B, while A suppresses the expression of B and vice versa. The dynamics of the network are described by two time-dependent ordinary differential equations:
M, A and B are the concentrations of each element, ka and kb are the synthesis rates (or maximum flux rates of A and B into a specific compartment of the system such as the nucleus), and kda and kdb are constants describing the decay (or irreversible inactivation such as protein degradation) of A and B, respectively. α and β are the cooperativities of repression by A of B and by B of A respectively, and μ is the cooperativity of induction by M. These cooperativities (or Hill coefficients) introduce non-linearity within the network and are important for its bistability [18, 19]. The cooperativity of induction (μ) is set to be the same for A and B, and kda and kdb are set to unity. Neither of these simplifications affects the overall behaviour of the system described in Additional file 1. M corresponds to the concentration of activin. The initial condition (at t = 0) for the simulations is (A, B) = (0, 0), reflecting the fact that expression of both Xbra and Gsc is induced by activin.
For the sake of simplicity we assume that M stays constant throughout the simulation. A few observations justify this assumption. First, cells respond by a 'ratchet mechanism' to the highest concentration of activin they are exposed to during their period of competence, and the timing of the response is related to the developmental stage of the cells, and not to the time of first exposure to activin . Second, the response of cells to activin is proportional to the absolute number of bound receptors that are internalised by the endocytic pathway and these remain active for several hours after a brief exposure of cells to activin . We note that the model described by equations (1) and (2) is a modified form of the dimensionless 'toggle model'  or 'repressor-repressor switch' .
Depending on the parameter values, the behaviour of the network illustrated in Fig. 1A can be classified as monostable or bistable (see Additional file 1). Bistability (i.e. the condition under which two stable states exist) requires the balanced rates of synthesis of A and B (ka and kb) (Fig. 2). Our result is consistent with previous studies showing that the product of cooperativity of the mutual repression between A and B must be greater than 1 (i.e. αβ >1) [18, 21], reflecting the fact that our model and theirs both have a mutual negative feedback motif at their cores. Most significantly, when certain conditions are met, the system reaches one of two opposing stable states (high A and low B, or vice versa) depending on the value of M (Fig. 1B), with bifurcation of the system occurring with a sharp threshold (Fig. 1C, see also Additional file 1). This behaviour resembles the concentration-dependent effects of a morphogen. The choice between the two opposing steady states with a given value of M depends on the parameter values and the initial conditions.
We have examined the conditions under which a threshold response can occur. Firstly, bistability is necessary but not sufficient for threshold formation (see Additional file 1). Secondly, the synthesis rates ka and kb must be well balanced for threshold formation to occur and the parameter choice here is more limited than for bistability (Fig. 3A, compare with Fig. 2). Thirdly, when the values of decay constants kda and kdb are different, the range of parameter choice for synthesis rates ka and kb becomes significantly larger and threshold values become more robust to parameter fluctuations (Fig. 3B). It also seems that, with a larger value of cooperativity of induction μ, threshold values are less sensitive to fluctuations of the other parameters (Fig. 3A–C). Interestingly, the establishment of the threshold proves to be rather insensitive to small changes of the parameters for cooperativity of repression α and β (Fig. 3C). Finally, the numerical simulations indicate that, if α = β and decay constants kda = kdb at the same time, a threshold cannot be established with any pair of values of synthesis rates ka and kb (Fig. 3C and data not shown). However, even if α = β, a threshold can be formed when appropriate parameter values of (ka, kb) and (kda, kdb) are chosen (Fig. 4).
The simulations also indicate that if α = β, the parameters must satisfy ka > kb and kda > kdb, or ka < kb and kda < kdb in order to create a threshold (Fig. 4 and data not shown). When ka < kb and kda < kdb, it seems always to be the case that A is on and B is off with low M at steady state, and vice versa with high M (Fig. 4 is an example). This principle of the system's behaviour might be explained as follows: at low M, the decay rate predominates over the synthesis rate, so that A, which has a lower decay rate, comes to be expressed at the expense of B. At high M, rates of synthesis dictate the behaviour of the system and B, with its higher synthesis rate, is eventually expressed at the expense of A.
The simulation in Fig. 1B shows the induction of A and B at two levels of M, one which results in expression of A at steady state and one which results in expression of B. Importantly, in both cases, the immediate response to M before steady state is to activate both A and B, and this behaviour recapitulates the expression dynamics of Xbra and Gsc in Xenopus animal pole regions after treatment with activin [22–24]. Incorporation of the transcription factor Xom in the simulation by adding an extra component in the scheme illustrated in Fig. 1A does not change the fundamental behaviour of the system and the dynamics of A and B can be categorised into three patterns: monostable, bistable without threshold, or bistable with threshold (see Additional file 1). This is presumably because the mutual negative relations between A and B stay the same. By choosing a set of appropriate parameters, it shows bistability with a threshold (Fig. 5).
Mutual exclusion of Xbra and Gscexpression at the single cell level
To ask whether the thresholds predicted by the simulation can exist at the single cell level in vivo, dissociated Xenopus animal pole cells were treated with different concentration of activin and cultured on fibronectin-coated glass. Cells were fixed after 7 hours and observed by indirect immunofluorescence microscopy using an anti-Xbra antibody. As described previously, Xbra was activated in a concentration-dependent manner, with maximal expression at 0.25 and 0.5 U/ml activin (Fig. 6). Interestingly, cells did not respond to activin uniformly; even at 0.5 U/ml some cells did not express Xbra (Fig. 6A, yellow arrowheads and 6B). The origin of this heterogeneity is unknown, although it may be due in part to the bifurcation properties of the system illustrated in Fig. 1 and also to the intrinsic stochastic nature of gene expression . In particular, the border (separatrix) that divides the phase plane is close to the point of the initial state (A = B = 0) when the value of M is near the threshold value (Fig. 1B). This would make the system sensitive to noise which, in a real in vivo system such as Xenopus animal cap cells, might influence the choice between the two stable states. Thus, two populations of cells (with high Xbra and low Gsc expression, or vice versa) would be produced when the concentration of activin is close to the threshold level. In fact, quantification of levels of fluorescence show that the expression level of Xbra is somewhat heterogeneous at both 0.25 and 0.5 U/ml activin, but much more so at 0.5 U/ml (Fig. 6B). It is likely that the higher of the two concentrations is close to the threshold at which Xbra expression is extinguished and Gsc is activated. The population of cells might therefore be a mixture of 'Xbra on' and 'Xbra off'. The lower concentration of 0.25 U/ml is likely to cause all cells to activate Xbra, and the lower level of heterogeneity observed at this concentration is likely to reflect the stochastic nature of gene expression.
The results described above suggest that the mutual exclusion of Xbra and Gsc expression occurs at the single cell level. To test this idea, we injected RNA encoding HA-tagged Gsc into Xenopus embryos at the one-cell stage and cultured these embryos to mid-blastula stage 9. Animal pole blastomeres derived from such embryos were mixed with animal pole cells from uninjected embryos, and they were then dissociated and cultured for 7 hours in the presence of 0.5 U/ml activin, after which Xbra and Gsc-HA were detected by indirect immunofluorescence. Xbra proved to be activated only in cells that expressed the lowest levels of Gsc-HA (Fig. 7A), a conclusion that was confirmed by measuring the levels of Xbra and Gsc-HA in each cell and plotting them against each other (Fig. 7B).
The results described in this paper show that the concentration-dependent induction of Xbra and Gsc by activin can be explained by the dynamics of a simple network with only three elements (Fig. 1A), in which a mutual negative feedback motif converts a graded signal (activin) into a binary output (Xbra on and Gsc off, or vice versa). The behaviour of mutual negative feedback motifs of this sort has been examined theoretically [18, 26] and experimentally , and it may provide a general mechanism for selective gene activation in response to different levels of a single external signal. Indeed, mutual negative feedback motifs play crucial roles in the behaviour of several biological processes, including cell type specification in the vertebrate neural tube [27, 28], embryonic segmentation  and photoreceptor cell fate decision in Drosophila . However, to our knowledge, the dynamic behaviour of such a mutual negative feedback motif to a single signal as illustrated in Fig. 1A has not previously been investigated. The results described in this paper demonstrate that such a system can display sharp thresholds. Although the now widely-accepted concept of a morphogen gradient requires threshold behaviours of this sort, rather few suggestions as to how thresholds are generated have so far been made [31, 32]. Our model provides a plausible mechanism by which a sharp boundary might be created between domains of different cell types in response to a morphogen gradient. It also explains how cells interpret a small concentration change in a morphogen gradient.
Our model builds on previous work [31, 32], and especially on experiments in Drosophila showing that different promoter affinities might underlie differential responses to the Bicoid morphogen gradient [33, 34]. In our model, promoter affinities and strengths are incorporated into the parameters of synthesis rates ka and kb. It seems that balanced but not necessarily equal rates of synthesis are required for threshold formation, but equally important are the stabilities of A and B, as represented by the decay rates kda and kdb. Indeed, the system becomes more robust when kda > kdb (or vice versa) because there are more choices of parameter pairs for synthesis (ka, kb) that permit threshold generation (compare Fig. 3B and 3C). It also means that thresholds become less sensitive to gene expression noise if the decay rates are different. Recent evidence suggests that pre-steady state interpretation of the Bicoid morphogen gradient contributes to the robustness of the system and ensures the accurate expression patterns of target genes . Although the morphogen level (M) stays the same in our model simulation, similar dynamic interpretations of changes in activin concentration might also provide robustness to threshold determination in our model. Our results also indicate that ka < kb and kda < kdb (or vice versa) must be satisfied in order to create a threshold if the other parameters such as the cooperativities of repression are the same. In principle these predictions might be tested by manipulating the promoter strengths of Xbra and Gsc together with the stabilities of their gene products, but such experiments would be technically challenging. Such attempts have also been hampered by the unavailability of anti-Gsc antisera.
Our model simulation recapitulates the induction and the co-expression of Xbra and Gsc immediately after exposure to activin, which cannot be deduced solely by the mutual repression between Xbra and Gsc. We note that the expression domains of Xbra and Gsc overlap at the early gastrula stage but then resolve by the end of gastrulation . And, also consistent with the model, the down regulation of Xbra in response to high concentrations of activin requires protein synthesis , indicating that mutual inhibition plays an essential role in this process. Although transient, the ubiquitous expression of Xbra in mesoderm is essential for morphogenetic movements during gastrulation [36, 37]. This illustrates the importance of the dynamic behaviour of a gene regulatory network as well as its steady state outcome, especially in the context of a dynamic process such as animal development.
Interestingly, we find that exposure to intermediate concentrations of activin causes different cells to make opposite binary decisions (Xbra on or off; Fig. 6). A similar heterogeneity in Xbra expression was found in single marginal zone cells of the Xenopus embryo at the early gastrula stage . This observation points to the importance of cell-cell communication in refining the borders between populations of different cell types created by a morphogen gradient. This communication may involve a community effect , perhaps requiring positive feedback between Xbra and FGF signalling . Lack of this positive feedback in dissociated Xenopus animal pole cells may also contribute to the observed heterogeneity in Xbra expression. It should be possible to incorporate such feedback to refine our model.
It is also possible, by combining additional feedback events, to create multiple thresholds in response to a morphogen gradient. For example, during neuronal cell fate specification in the ventral neural tube, 'class I' and 'class II' homeodomain transcription factors, which are regulated by the morphogen Sonic Hedgehog, inhibit each other's expression [27, 28]. This sort of combinatorial mechanism may be made more general (see Additional file 2).
Our simulation of a simple network of three elements with a mutual negative feedback motif recapitulates the dynamics of the concentration-dependent induction of Xbra and Gsc by activin in Xenopus. Our model provides a mechanism by which a sharp boundary might be created between domains of different cell types in response to a morphogen gradient.
Numerical simulations and calculations
Numerical simulations and calculations were performed using Mathematica (Ver 5.2, Wolfram Research). The package of functions and programs used is available upon request. Part of the algorithm was adopted from Murrell .
Xenopusembryo manipulation and microinjection
Fertilisation and culture of Xenopus embryos were performed as described  and recombinant human activin A was obtained as described . RNA injection was performed according to Smith . Animal pole regions were dissected from blastulae and cells were dissociated in calcium- and magnesium-free medium (Tris 75 mM pH7.5, NaCl 880 mM, KCl 10 mM, NaHCO3 24 mM) for 30–45 minutes at room temperature before treating with activin. Gsc tagged with an HA epitope at its C-terminus, cloned in the vector pCS2+, was transcribed in vitro using mMESSAGE mMACHINE (Ambion).
For the experiment shown in Fig. 7, fertilised Xenopus embryos were injected with 200 pg RNA encoding HA-tagged Gsc (Gsc-HA). Animal pole regions derived from these injected embryos and from uninjected embryos were mixed in a ratio of 5:1. Blastomeres derived from this mixed population of animal pole regions were dissociated and treated with 0.5 U/ml of activin. They were fixed after 7 hours of culture and stained with Hoechst 33342 and processed for anti-Xbra and anti-HA staining. Dissociated animal cap cells were cultured in 75% Normal Amphibian Medium  containing 0.1% bovine serum albumin (Sigma, Fraction V, A9647) in glass-bottomed microwell dishes (MatTek Co. P35GC-1.5-14-C) that had been covered overnight with a 0.002% solution of fibronectin (Sigma, F0895) and washed once with water before use.
Indirect immunofluorescence microscopy and quantification of fluorescence
Cells were fixed in MEMFA  containing 1 mM EDTA and incubated with anti-Xbra  (1:1000) and anti-HA (1:10,000, clone 3F10 high affinity, Roche) antibodies. Secondary antibodies were Cy3-conjugated goat anti-rabbit IgG (H+L) (1.5 μg/ml, Jackson ImmunoResearch Laboratories) and Alexa488-conjugated anti-mouse IgG (2 μg/ml, Molecular Probes). Hoechst33342 (0.5 μg/ml, Molecular Probes) was used to visualize nuclear DNA. Images were acquired with a digital CCD camera (Hamamatsu C4742-95 Orca100) mounted on Leica DMIRB microscope, using Openlab software (Improvision). Quantification of fluorescence was performed using ImageJ . Arbitrary fluorescence figures in Fig 7 were defined by the mean grey value of anti-Xbra (or anti-HA) fluorescence divided by mean grey value of DNA fluorescence in a region of interest (a nucleus).
Ashe HL, Briscoe J: The interpretation of morphogen gradients. Development. 2006, 133 (3): 385-394. 10.1242/dev.02238.
Gurdon JB, Bourillot PY: Morphogen gradient interpretation. Nature. 2001, 413 (6858): 797-803. 10.1038/35101500.
Wolpert L: Positional information and the spatial pattern of cellular differentiation. J Theor Biol. 1969, 25: 1-47. 10.1016/S0022-5193(69)80016-0.
Scholpp S, Brand M: Endocytosis controls spreading and effective signaling range of Fgf8 protein. Curr Biol. 2004, 14 (20): 1834-1841. 10.1016/j.cub.2004.09.084.
Kinoshita T, Jullien J, Gurdon JB: Two-dimensional morphogen gradient in Xenopus: boundary formation and real-time transduction response. Dev Dyn. 2006, 235 (12): 3189-3198. 10.1002/dvdy.20963.
Williams PH, Hagemann A, Gonzalez-Gaitan M, Smith JC: Visualizing long-range movement of the morphogen Xnr2 in the Xenopus embryo. Curr Biol. 2004, 14 (21): 1916-1923. 10.1016/j.cub.2004.10.020.
Stathopoulos A, Levine M: Dorsal gradient networks in the Drosophila embryo. Dev Biol. 2002, 246 (1): 57-67. 10.1006/dbio.2002.0652.
Piepenburg O, Grimmer D, Williams PH, Smith JC: Activin redux: specification of mesodermal pattern in Xenopus by graded concentrations of endogenous activin B. Development. 2004, 131 (20): 4977-4986. 10.1242/dev.01323.
Green JBA, New HV, Smith JC: Responses of embryonic Xenopus cells to activin and FGF are separated by multiple dose thresholds and correspond to distinct axes of the mesoderm. Cell. 1992, 71 (5): 731-739. 10.1016/0092-8674(92)90550-V.
Shimizu K, Gurdon JB: A quantitative analysis of signal transduction from activin receptor to nucleus and its relevance to morphogen gradient interpretation. Proc Natl Acad Sci U S A. 1999, 96 (12): 6791-6796. 10.1073/pnas.96.12.6791.
Artinger M, Blitz I, Inoue K, Tran U, Cho KW: Interaction of goosecoid and brachyury in Xenopus mesoderm patterning. Mech Dev. 1997, 65 (1-2): 187-196. 10.1016/S0925-4773(97)00073-7.
Latinkic BV, Umbhauer M, Neal KA, Lerchner W, Smith JC, Cunliffe V: The Xenopus Brachyury promoter is activated by FGF and low concentrations of activin and suppressed by high concentrations of activin and by paired-type homeodomain proteins. Genes Dev. 1997, 11 (23): 3265-3276.
Trindade M, Tada M, Smith JC: DNA-binding specificity and embryological function of Xom (Xvent-2). Dev Biol. 1999, 216 (2): 442-456. 10.1006/dbio.1999.9507.
Ladher R, Mohun TJ, Smith JC, Snape AM: Xom: a Xenopus homeobox gene that mediates the early effects of BMP-4. Development. 1996, 122 (8): 2385-2394.
Onichtchouk D, Gawantka V, Dosch R, Delius H, Hirschfeld K, Blumenstock C, Niehrs C: The Xvent-2 homeobox gene is part of the BMP-4 signalling pathway controlling [correction of controling] dorsoventral patterning of Xenopus mesoderm. Development. 1996, 122 (10): 3045-3053.
Papalopulu N, Kintner C: A Xenopus gene, Xbr-1, defines a novel class of homeobox genes and is expressed in the dorsal ciliary margin of the eye. Dev Biol. 1996, 174 (1): 104-114. 10.1006/dbio.1996.0055.
Schmidt JE, von Dassow G, Kimelman D: Regulation of dorsal-ventral patterning: the ventralizing effects of the novel Xenopus homeobox gene Vox. Development. 1996, 122 (6): 1711-1721.
Cherry JL, Adler FR: How to make a biological switch. J Theor Biol. 2000, 203 (2): 117-133. 10.1006/jtbi.2000.1068.
Ferrell JE: Self-perpetuating states in signal transduction: positive feedback, double-negative feedback and bistability. Curr Opin Cell Biol. 2002, 14 (2): 140-148. 10.1016/S0955-0674(02)00314-9.
Jullien J, Gurdon J: Morphogen gradient interpretation by a regulated trafficking step during ligand-receptor transduction. Genes Dev. 2005, 19 (22): 2682-2694. 10.1101/gad.341605.
Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli. Nature. 2000, 403 (6767): 339-342. 10.1038/35002131.
Green JB, Smith JC, Gerhart JC: Slow emergence of a multithreshold response to activin requires cell-contact-dependent sharpening but not prepattern. Development. 1994, 120 (8): 2271-2278.
Papin C, Smith JC: Gradual refinement of activin-induced thresholds requires protein synthesis. Dev Biol. 2000, 217 (1): 166-172. 10.1006/dbio.1999.9531.
Wardle FC, Smith JC: Refinement of gene expression patterns in the early Xenopus embryo. Development. 2004, 131 (19): 4687-4696. 10.1242/dev.01340.
Kaern M, Elston TC, Blake WJ, Collins JJ: Stochasticity in gene expression: from theories to phenotypes. Nat Rev Genet. 2005, 6 (6): 451-464. 10.1038/nrg1615.
Cinquin O, Page KM: Generalized, switch-like competitive heterodimerization networks. Bull Math Biol. 2007, 69 (2): 483-494. 10.1007/s11538-006-9114-8.
Briscoe J, Pierani A, Jessell TM, Ericson J: A homeodomain protein code specifies progenitor cell identity and neuronal fate in the ventral neural tube. Cell. 2000, 101 (4): 435-445. 10.1016/S0092-8674(00)80853-3.
Marquardt T, Pfaff SL: Cracking the transcriptional code for cell specification in the neural tube. Cell. 2001, 106 (6): 651-654. 10.1016/S0092-8674(01)00499-8.
Jaeger J, Blagov M, Kosman D, Kozlov KN, Manu, Myasnikova E, Surkova S, Vanario-Alonso CE, Samsonova M, Sharp DH, Reinitz J: Dynamical analysis of regulatory interactions in the gap gene system of Drosophila melanogaster. Genetics. 2004, 167 (4): 1721-1737. 10.1534/genetics.104.027334.
Mikeladze-Dvali T, Wernet MF, Pistillo D, Mazzoni EO, Teleman AA, Chen YW, Cohen S, Desplan C: The growth regulators warts/lats and melted interact in a bistable loop to specify opposite fates in Drosophila R8 photoreceptors. Cell. 2005, 122 (5): 775-787. 10.1016/j.cell.2005.07.026.
Lewis J, Slack JM, Wolpert L: Thresholds in development. J Theor Biol. 1977, 65 (3): 579-590. 10.1016/0022-5193(77)90216-8.
Goldbeter A, Wolpert L: Covalent modification of proteins as a threshold mechanism in development. J Theor Biol. 1990, 142 (2): 243-250. 10.1016/S0022-5193(05)80225-5.
Driever W, Thoma G, Nusslein-Volhard C: Determination of spatial domains of zygotic gene expression in the Drosophila embryo by the affinity of binding sites for the bicoid morphogen. Nature. 1989, 340 (6232): 363-367. 10.1038/340363a0.
Struhl G, Struhl K, Macdonald PM: The gradient morphogen bicoid is a concentration-dependent transcriptional activator. Cell. 1989, 57 (7): 1259-1273. 10.1016/0092-8674(89)90062-7.
Bergmann S, Sandler O, Sberro H, Shnider S, Schejter E, Shilo BZ, Barkai N: Pre-steady-state decoding of the Bicoid morphogen gradient. PLoS Biol. 2007, 5 (2): e46-10.1371/journal.pbio.0050046.
Tada M, Smith JC: Xwnt11 is a target of Xenopus Brachyury: regulation of gastrulation movements via Dishevelled, but not through the canonical Wnt pathway. Development. 2000, 127 (127): 2227-2238.
Heisenberg CP, Tada M, Rauch GJ, Saude L, Concha ML, Geisler R, Stemple DL, Smith JC, Wilson SW: Silberblick/Wnt11 mediates convergent extension movements during zebrafish gastrulation. Nature. 2000, 405 (6782): 76-81. 10.1038/35011068.
Gurdon JB: A community effect in animal development. Nature. 1988, 336 (6201): 772-774. 10.1038/336772a0.
Isaacs HV, Pownall ME, Slack JM: eFGF regulates Xbra expression during Xenopus gastrulation. Embo J. 1994, 13 (19): 4469-4481.
Murrell H: Planer phaseplots and bifurcation animations. Mathematica J. 1994, 4 (3):
Tada M, O'Reilly MA, Smith JC: Analysis of competence and of Brachyury autoinduction by use of hormone-inducible Xbra. Development. 1997, 124 (11): 2225-2234.
Smith JC: Purifying and assaying mesoderm-inducing factors from vertebrate embryos. Cellular Interactions in Development - a Practical Approach. Edited by: Hartley D. 1993, Oxford , Oxford University Press, 181-204.
Slack JMW: Regional biosynthetic markers in the early amphibian embryo. J Embryol exp Morph. 1984, 80: 289-319.
Cunliffe V, Smith JC: Specification of mesodermal pattern in Xenopus laevis by interactions between Brachyury, noggin and Xwnt-8. Embo J. 1994, 13 (2): 349-359.
Rasband W: ImageJ. [http://rsb.info.nih.gov/ij/index.html]
We thank Joanna Argasinska for the Gsc-HA plasmid, Kevin Dingwell and Fiona Wardle for comments on the manuscript, and John Gurdon and James Briscoe for discussion. This work was supported by the Wellcome Trust.
YS conceived and carried out the studies, and drafted the manuscript. JCS participated in the design of the study and helped to draft the manuscript. Both authors read and approved the final manuscript.
Electronic supplementary material
Additional File 1: Three possible behaviours for the system described in Figure 1A: monostable, bistable without a threshold, and bistable with a threshold. This system has two variables (A and B) that change over time according to equations (1) and (2) in the main text. The state of the system at a given time point is defined by the values (A, B) and can be plotted as a single point (i.e. phase point) on the two-dimensional phase plane illustrated in Fig. 1B. The phase point moves around the phase space with a speed defined by the equations (1) and (2) in the main text. This system like others with a mutual negative feedback motif undergoes a saddle-node bifurcation (Kobayashi H et al (2004) Proc Natl Acad Sci USA 101: 8414–8419). Nodes and saddles are the points on the phase plane where the speed of the phase point is zero, and therefore derive from the solution of a coupled equation where the left hand sides of equations (1) and (2) are replaced with 0. In our system, there are three fixed points when bifurcation occurs, two of which are stable nodes and the other is an unstable saddle point. A saddle point is analogous to a pointed summit where a ball never rests. At steady state, the phase point rests at one of the two stable nodes (see Fig. 1B). The figure shows how the positions of these nodes and the saddle change in the phase plane (A, B) as M increases from 0 to 8 as indicated by the arrows. Green curves indicate the nodes where the phase point reaches at steady state after starting from (0, 0). (A) Monostable. There is just one node irrespective of the value of M. (B) Bistable without threshold. When M is small there is only one node. As M increases, bifurcation occurs (red triangle) and another node appears, but the phase point always ends up at the original node at steady state. (C) Bistable with a threshold. When M reaches the threshold (yellow diamond), the position of the phase point at steady state shifts from the original node to the other node created at the bifurcation. Bistability is necessary but not sufficient for a threshold to be generated. (PDF 20 KB)
Additional File 2: Creation of multiple thresholds by combination of mutual negative feedback motifs. (A) A network structure with two mutual negative feedback motifs, which are controlled by a single morphogen (denoted by M) independently. (B) Multiple thresholds in a morphogen gradient. If the two mutual negative feedback motifs have different threshold values, the responding tissue is divided into three compartments (differently coloured) with two sharp boundaries. These compartments express different combinations of genes according to their position in the morphogen gradient. In theory, any number of mutual negative feedback motifs can be incorporated in the scheme, thereby generating multiple thresholds. (PDF 3 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Saka, Y., Smith, J.C. A mechanism for the sharp transition of morphogen gradient interpretation in Xenopus. BMC Dev Biol 7, 47 (2007). https://doi.org/10.1186/1471-213X-7-47