Morphogen gradients, positional information, and Xenopus: interplay of theory and experiment

Evaluation of temperature- and ethanol-related developmental degree variations by a new scoring system (FETAX-score) applicable to Frog Embryo Teratogenicity Assay: Xenopus

The aim of the present work is to propose a new quantitative assessment method (FETAX-score) for determining the degree of Xenopus laevis embryo development intended for use in embryotoxicity studies. Inspired by a similar scoring system used to evaluate developmental delays (young-for-age phenotypes) in rat embryos cultured in vitro, the FETAX-score was established by considering seven morphological features (head, naris, mouth, lower jaw, tentacles, intestine, anus) that are easily evaluable in tadpoles during the late stages of development at the conclusion of the test. Given that X. laevis development is temperature-dependent and that temperatures below 14°C and above 26°C are teratogenic, the FETAX-score was tested in embryos maintained at 17, 20, 23 and 26°C. No abnormalities were observed in any group, while the total score was temperature-related, suggesting that the FETAX-score is sensitive to moderate distress that does not influence general morphology. Intestine and anus were the least sensitive structures to temperature variations. To assess the applicability of the FETAX-score in developmental toxicological studies, we evaluated FETAX-score in tadpoles exposed during the morphogenetic period to Ethanol (Eth) at concentrations of 0, 0.25, 0.5, 1, 1.5, and 2 % v/v. Gross malformations were observed only in tadpoles from the Eth 2 % group. By contrast, data analysis of the other Eth groups showed dose-related reductions in the FETAX-score. Tentacles were the most sensitive structures to Eth-related delays. These results support the use of the FETAX-score to quantitatively assess developmental deviations in FETAX embryotoxicity studies.

Coupling and uncoupling of midline morphogenesis and cell flow in amniote gastrulation

Large-scale cell flow characterizes gastrulation in animal development. In amniote gastrulation, particularly in avian gastrula, a bilateral vortex-like counter-rotating cell flow, called 'polonaise movements', appears along the midline. Here, through experimental manipulations, we addressed relationships between the polonaise movements and morphogenesis of the primitive streak, the earliest midline structure in amniotes. Suppression of the Wnt/planar cell polarity (PCP) signaling pathway maintains the polonaise movements along a deformed primitive streak. Mitotic arrest leads to diminished extension and development of the primitive streak and maintains the early phase of the polonaise movements. Ectopically induced Vg1, an axis-inducing morphogen, generates the polonaise movements, aligned to the induced midline, but disturbs the stereotypical cell flow pattern at the authentic midline. Despite the altered cell flow, induction and extension of the primitive streak are preserved along both authentic and induced midlines. Finally, we show that ectopic axis-inducing morphogen, Vg1, is capable of initiating the polonaise movements without concomitant PS extension under mitotic arrest conditions. These results are consistent with a model wherein primitive streak morphogenesis is required for the maintenance of the polonaise movements, but the polonaise movements are not necessarily responsible for primitive streak morphogenesis. Our data describe a previously undefined relationship between the large-scale cell flow and midline morphogenesis in gastrulation.

Coupling and uncoupling of midline morphogenesis and cell flow in amniote gastrulation

Large-scale cell flow characterizes gastrulation in animal development. In amniote gastrulation, particularly in avian gastrula, a bilateral vortex-like counter-rotating cell flow, called 'polonaise movements', appears along the midline. Here, through experimental manipulations, we addressed relationships between the polonaise movements and morphogenesis of the primitive streak, the earliest midline structure in amniotes. Suppression of the Wnt/planar cell polarity (PCP) signaling pathway maintains the polonaise movements along a deformed primitive streak. Mitotic arrest leads to diminished extension and development of the primitive streak and maintains the early phase of the polonaise movements. Ectopically induced Vg1, an axis-inducing morphogen, generates the polonaise movements, aligned to the induced midline, but disturbs the stereotypical cell flow pattern at the authentic midline. Despite the altered cell flow, induction and extension of the primitive streak are preserved along both authentic and induced midlines. Finally, we show that ectopic axis-inducing morphogen, Vg1, is capable of initiating the polonaise movements without concomitant PS extension under mitotic arrest conditions. These results are consistent with a model wherein primitive streak morphogenesis is required for the maintenance of the polonaise movements, but the polonaise movements are not necessarily responsible for primitive streak morphogenesis. Our data describe a previously undefined relationship between the large-scale cell flow and midline morphogenesis in gastrulation.

Reconstructing gene regulatory networks of biological function using differential equations of multilayer perceptrons

BACKGROUND

Building biological networks with a certain function is a challenge in systems biology. For the functionality of small (less than ten nodes) biological networks, most methods are implemented by exhausting all possible network topological spaces. This exhaustive approach is difficult to scale to large-scale biological networks. And regulatory relationships are complex and often nonlinear or non-monotonic, which makes inference using linear models challenging.

RESULTS

In this paper, we propose a multi-layer perceptron-based differential equation method, which operates by training a fully connected neural network (NN) to simulate the transcription rate of genes in traditional differential equations. We verify whether the regulatory network constructed by the NN method can continue to achieve the expected biological function by verifying the degree of overlap between the regulatory network discovered by NN and the regulatory network constructed by the Hill function. And we validate our approach by adapting to noise signals, regulator knockout, and constructing large-scale gene regulatory networks using link-knockout techniques. We apply a real dataset (the mesoderm inducer Xenopus Brachyury expression) to construct the core topology of the gene regulatory network and find that Xbra is only strongly expressed at moderate levels of activin signaling.

CONCLUSION

We have demonstrated from the results that this method has the ability to identify the underlying network topology and functional mechanisms, and can also be applied to larger and more complex gene network topologies.

Patterning principles of morphogen gradients

Metazoan embryos develop from a single cell into three-dimensional structured organisms while groups of genetically identical cells attain specialized identities. Cells of the developing embryo both create and accurately interpret morphogen gradients to determine their positions and make specific decisions in response. Here, we first cover intellectual roots of morphogen and positional information concepts. Focusing on animal embryos, we then provide a review of current understanding on how morphogen gradients are established and how their spans are controlled. Lastly, we cover how gradients evolve in time and space during development, and how they encode information to control patterning. In sum, we provide a list of patterning principles for morphogen gradients and review recent advances in quantitative methodologies elucidating information provided by morphogens.

The developing optic tectum: An asymmetrically organized system and the need for a redefinition of the notion of sensitive period

The present article summarizes the main events involved in the isthmic organizer and optic tectum determination and analyses how optic tectum patterning is translated, by the organized operation of several specific cell behaviors, into the terminally differentiated optic tectum. The paper proposes that this assembling of temporally/spatially organized cell behaviors could be incorporated into a wider notion of patterning and that, given the asymmetric organization of the developing optic tectum, the notion of "sensitive period" does not capture the whole complexity of midbrain development and the pathogenesis of congenital disorders. The cell behaviors involved in the optic tectum development are organized in time and space by the isthmic organizer. A comprehensive description of the normal optic tectum development, and also its alterations, should consider both domains. Significantly, the identity of each neuronal cohort depends critically on its "time and place of birth". Both parameters must be considered at once to explain how the structural and functional organization of the optic tectum is elaborated. The notion of "patterning" applies only to the early events of the optic tectum development. Besides, the notion of "sensitive period" considers only a temporal domain and disregards the asymmetric organization of the developing optic tectum. The present paper proposes that these notions might be re-defined: (a) a wider meaning of the term patterning and (b) a replacement of the term "sensitive period" by a more precise concept of "sensitive temporal/spatial window".

Extracellular interactions and ligand degradation shape the nodal morphogen gradient

The correct distribution and activity of secreted signaling proteins called morphogens is required for many developmental processes. Nodal morphogens play critical roles in embryonic axis formation in many organisms. Models proposed to generate the Nodal gradient include diffusivity, ligand processing, and a temporal activation window. But how the Nodal morphogen gradient forms in vivo remains unclear. Here, we have measured in vivo for the first time, the binding affinity of Nodal ligands to their major cell surface receptor, Acvr2b, and to the Nodal inhibitor, Lefty, by fluorescence cross-correlation spectroscopy. We examined the diffusion coefficient of Nodal ligands and Lefty inhibitors in live zebrafish embryos by fluorescence correlation spectroscopy. We also investigated the contribution of ligand degradation to the Nodal gradient. We show that ligand clearance via degradation shapes the Nodal gradient and correlates with its signaling range. By computational simulations of gradient formation, we demonstrate that diffusivity, extra-cellular interactions, and selective ligand destruction collectively shape the Nodal morphogen gradient.

DOI:http://dx.doi.org/10.7554/eLife.13879.001

Animals develop from a single fertilized egg cell into multicellular organisms. This process requires chemical signals called “morphogens” that instruct the cells how to behave during development. The morphogens move across cells and tissues to form gradients of the signal. Cells then respond in different ways depending on how much of the signal they receive. This, in turn, depends on several factors: first, how quickly or slowly the signal moves; second, how well the morphogen binds to responding cells and other molecules in its path; and third, how much signal is lost or destroyed during the movement.

Many researchers study morphogen gradients in the transparent zebrafish, since it grows quickly and it is easy to see developmental changes. However, until now it was not fully clear how the well-known morphogen called Nodal moves in live zebrafish as they develop.

Wang, Wang et al. have now investigated how well Nodal signals bind to the surface of cells that receive the signal and to a molecule called “Lefty”, which is present in the same path and interferes with Nodal signals. Advanced techniques called fluorescence correlation and cross-correlation spectroscopy were used to measure Nodal signals at the level of single molecules in growing zebrafish. The experiments gave insights into how far Nodal signals move and remain active. The results showed that, in addition to Nodal diffusing and binding to receiving cells, one of the most important factors determining how far and quickly Nodal moves is its inactivation and destruction. Lastly, Wang, Wang et al. built computational models to test their observations from live zebrafish.

The current work was based on forcing zebrafish to produce molecules including Nodal at locations within the fish that normally do not make them. Therefore future experiments will aim to examine these molecules and their interactions when they are produced at their normal locations in the animal over time.

DOI:http://dx.doi.org/10.7554/eLife.13879.002

Multicellular Mathematical Modelling of Mesendoderm Formation in Amphibians

The earliest cell fate decisions in a developing embryo are those associated with establishing the germ layers. The specification of the mesoderm and endoderm is of particular interest as the mesoderm is induced from the endoderm, potentially from an underlying bipotential group of cells, the mesendoderm. Mesendoderm formation has been well studied in an amphibian model frog, Xenopus laevis, and its formation is driven by a gene regulatory network (GRN) induced by maternal factors deposited in the egg. We have recently demonstrated that the axolotl, a urodele amphibian, utilises a different topology in its GRN to specify the mesendoderm. In this paper, we develop spatially structured mathematical models of the GRNs governing mesendoderm formation in a line of cells. We explore several versions of the model of mesendoderm formation in both Xenopus and the axolotl, incorporating the key differences between these two systems. Model simulations are able to reproduce known experimental data, such as Nodal expression domains in Xenopus, and also make predictions about how the positional information derived from maternal factors may be interpreted to drive cell fate decisions. We find that whilst cell-cell signalling plays a minor role in Xenopus, it is crucial for correct patterning domains in axolotl.

Patterning of the angiosperm female gametophyte through the prism of theoretical paradigms

The FG (female gametophyte) of flowering plants (angiosperms) is a simple highly polar structure composed of only a few cell types. The FG develops from a single cell through mitotic divisions to generate, depending on the species, four to 16 nuclei in a syncytium. These nuclei are then partitioned into three or four distinct cell types. The mechanisms underlying the specification of the nuclei in the FG has been a focus of research over the last decade. Nevertheless, we are far from understanding the patterning mechanisms that govern cell specification. Although some results were previously interpreted in terms of static positional information, several lines of evidence now show that local interactions are important. In the present article, we revisit the available data on developmental mutants and cell fate markers in the light of theoretical frameworks for biological patterning. We argue that a further dissection of the mechanisms may be impeded by the combinatorial and dynamical nature of developmental cues. However, accounting for these properties of developing systems is necessary to disentangle the diversity of the phenotypic manifestations of the underlying molecular interactions.

Spatiotemporal proliferation of human stromal cells adjusts to nutrient availability and leads to stanniocalcin-1 expression in vitro and in vivo

Cells and tissues are intrinsically adapted to molecular gradients and use them to maintain or change their activity. The effect of such gradients is particularly important for cell populations that have an intrinsic capacity to differentiate into multiple cell lineages, such as bone marrow derived mesenchymal stromal cells (MSCs). Our results showed that nutrient gradients prompt the spatiotemporal organization of MSCs in 3D culture. Cells adapted to their 3D environment without significant cell death or cell differentiation. Kinetics data and whole-genome gene expression analysis suggest that a low proliferation activity phenotype predominates in stromal cells cultured in 3D, likely due to increasing nutrient limitation. These differences implied that despite similar surface areas available for cell attachment, higher cell concentrations in 3D reduced MSCs proliferation, while activating hypoxia related-pathways. To further understand the in vivo effects of both proliferation and cell concentrations, we increased cell concentrations in small (1.8 μl) implantable wells. We found that MSCs accumulation and conditioning by nutrient competition in small volumes leads to an ideal threshold of cell-concentration for the induction of blood vessel formation, possibly signaled by the hypoxia-related stanniocalcin-1 gene.

Dynamics of gene circuits shapes evolvability

To what extent does the dynamical mechanism producing a specific biological phenotype bias the ability to evolve into novel phenotypes? We use the interpretation of a morphogen gradient into a single stripe of gene expression as a model phenotype. Although there are thousands of three-gene circuit topologies that can robustly develop a stripe of gene expression, the vast majority of these circuits use one of just six fundamentally different dynamical mechanisms. Here we explore the potential for gene circuits that use each of these six mechanisms to evolve novel phenotypes such as multiple stripes, inverted stripes, and gradients of gene expression. Through a comprehensive and systematic analysis, we find that circuits that use alternative mechanisms differ in the likelihood of reaching novel phenotypes through mutation. We characterize the phenotypic transitions and identify key ingredients of the evolutionary potential, such as sensitive interactions and phenotypic hubs. Finally, we provide an intuitive understanding on how the modular design of a particular mechanism favors the access to novel phenotypes. Our work illustrates how the dynamical mechanism by which an organism develops constrains how it can evolve. It is striking that these dynamical mechanisms and their impact on evolvability can be observed even for such an apparently simple patterning task, performed by just three-node circuits.

Brief report: Rx1 defines retinal precursor identity by repressing alternative fates through the activation of TLE2 and Hes4

The molecular mechanisms underlying the acquisition of retinal precursor identity are scarcely defined. Although the homeobox gene Rx1 (also known as Rax) plays a major role in specifying retinal precursors and maintaining their multipotent state, the involved mechanisms remain to be largely deciphered. Here, following a highthroughput screen for genes regulated by Rx1, we found that this transcription factor specifies the fate of retinal progenitors by repressing genes normally activated in adjacent ectodermal territories. Unexpectedly, we also observed that Rx1, mainly through the activation of the transcriptional repressors TLE2 and Hes4, is necessary and sufficient to inhibit endomesodermal gene expression in retinal precursors of the eye field. In particular, Rx1 knockdown leads retinogenic blastomeres to adopt an endomesodermal fate, indicating a previously undescribed function for Rx1 in preventing the expression of endomesoderm determinants known to inhibit retinal fate. Altogether these data suggest that an essential requirement to establish a retinal precursor identity is the active inhibition of pathways leading to alternative fates.

Design principles of stripe-forming motifs: the role of positive feedback

Interpreting a morphogen gradient into a single stripe of gene-expression is a fundamental unit of patterning in early embryogenesis. From both experimental data and computational studies the feed-forward motifs stand out as minimal networks capable of this patterning function. Positive feedback within gene networks has been hypothesised to enhance the sharpness and precision of gene-expression borders, however a systematic analysis has not yet been reported. Here we set out to assess this hypothesis, and find an unexpected result. The addition of positive-feedback can have different effects on two different designs of feed-forward motif– it increases the parametric robustness of one design, while being neutral or detrimental to the other. These results shed light on the abundance of the former motif and especially of mutual-inhibition positive feedback in developmental networks.

Quantitative imaging of morphogen gradients in Drosophila imaginal discs

Cells at different positions in a developing tissue receive different concentrations of signaling molecules, called morphogens, and this influences their cell fate. Morphogen concentration gradients have been proposed to control patterning as well as growth in many developing tissues. Some outstanding questions about tissue patterning by morphogen gradients are the following: What are the mechanisms that regulate gradient formation and shape? Is the positional information encoded in the gradient sufficiently precise to determine the positions of target gene domain boundaries? What are the temporal dynamics of gradients and how do they relate to patterning and growth? These questions are inherently quantitative in nature and addressing them requires measuring morphogen concentrations in cells, levels of downstream signaling activity, and kinetics of morphogen transport. Here we first present methods for quantifying morphogen gradient shape in which the measurements can be calibrated to reflect actual morphogen concentrations. We then discuss using fluorescence recovery after photobleaching to study the kinetics of morphogen transport at the tissue level. Finally, we present particle tracking as a method to study morphogen intracellular trafficking.

Whole-genome microRNA screening identifies let-7 and mir-18 as regulators of germ layer formation during early embryogenesis

Tight control over the segregation of endoderm, mesoderm, and ectoderm is essential for normal embryonic development of all species, yet how neighboring embryonic blastomeres can contribute to different germ layers has never been fully explained. We postulated that microRNAs, which fine-tune many biological processes, might modulate the response of embryonic blastomeres to growth factors and other signals that govern germ layer fate. A systematic screen of a whole-genome microRNA library revealed that the let-7 and miR-18 families increase mesoderm at the expense of endoderm in mouse embryonic stem cells. Both families are expressed in ectoderm and mesoderm, but not endoderm, as these tissues become distinct during mouse and frog embryogenesis. Blocking let-7 function in vivo dramatically affected cell fate, diverting presumptive mesoderm and ectoderm into endoderm. siRNA knockdown of computationally predicted targets followed by mutational analyses revealed that let-7 and miR-18 down-regulate Acvr1b and Smad2, respectively, to attenuate Nodal responsiveness and bias blastomeres to ectoderm and mesoderm fates. These findings suggest a crucial role for the let-7 and miR-18 families in germ layer specification and reveal a remarkable conservation of function from amphibians to mammals.

Myotubularin-related protein 4 (MTMR4) attenuates BMP/Dpp signaling by dephosphorylation of Smad proteins

Bone morphogenetic proteins (BMPs) signaling essentially regulates a wide range of biological responses. Although multiple regulators at different layers of the receptor-effectors axis have been identified, the mechanisms of homeostatic BMP signaling remain vague. Herein we demonstrated that myotubularin-related protein 4 (MTMR4), a FYVE domain-containing dual-specificity protein phosphatase (DUSP), preferentially associated with and dephosphorylated the activated R-Smads in cytoplasm, which is a critical checkpoint in BMP signal transduction. Therefore, transcriptional activation by BMPs was tightly controlled by the expression level and the intrinsic phosphatase activity of MTMR4. More profoundly, ectopic expression of MTMR4 or its Drosophila homolog CG3632 genetically interacted with BMP/Dpp signaling axis in regulation of the vein development of Drosophila wings. By doing so, MTMR4 could interact with and dephosphorylate Mothers against Decapentaplegic (Mad), the sole R-Smad in Drosophila BMP pathway, and hence affected the target genes expression of Mad. In conclusion, this study has suggested that MTMR4 is a necessary negative modulator for the homeostasis of BMP/Dpp signaling.

Morphogengineering roots: comparing mechanisms of morphogen gradient formation

BACKGROUND

In developmental biology, there has been a recent focus on the robustness of morphogen gradients as possible providers of positional information. It was shown that functional morphogen gradients present strong biophysical constraints and lack of robustness to noise. Here we explore how the details of the mechanism which underlies the generation of a morphogen gradient can influence those properties.

RESULTS

We contrast three gradient-generating mechanisms, (i) a source-decay mechanism; and (ii) a unidirectional transport mechanism; and (iii) a so-called reflux-loop mechanism. Focusing on the dynamics of the phytohormone auxin in the root, we show that only the reflux-loop mechanism can generate a gradient that would be adequate to supply functional positional information for the Arabidopsis root, for biophysically reasonable kinetic parameters.

CONCLUSIONS

We argue that traits that differ in spatial and temporal time-scales can impose complex selective pressures on the mechanism of morphogen gradient formation used for the development of the particular organism.

On the evolution of morphogenetic models: mechano-chemical interactions and an integrated view of cell differentiation, growth, pattern formation and morphogenesis

In the 1950s, embryology was conceptualized as four relatively independent problems: cell differentiation, growth, pattern formation and morphogenesis. The mechanisms underlying the first three traditionally have been viewed as being chemical in nature, whereas those underlying morphogenesis have usually been discussed in terms of mechanics. Often, morphogenesis and its mechanical processes have been regarded as subordinate to chemical ones. However, a growing body of evidence indicates that the biomechanics of cells and tissues affect in striking ways those phenomena often thought of as mainly under the control of cell-cell signalling. This accumulation of data has led to a revival of the mechano-transduction concept in particular, and of complexity in general, causing us now to consider whether we should retain the traditional conceptualization of development. The researchers' semantic preferences for the terms 'patterning', 'pattern formation' or 'morphogenesis' can be used to describe three main 'schools of thought' which emerged in the late 1970s. In the 'molecular school', the term patterning is deeply tied to the positional information concept. In the 'chemical school', the term 'pattern formation' regularly implies reaction-diffusion models. In the 'mechanical school', the term 'morphogenesis' is more frequently used in relation to mechanical instabilities. Major differences among these three schools pertain to the concept of self-organization, and models can be classified as morphostatic or morphodynamic. Various examples illustrate the distorted picture that arises from the distinction among differentiation, growth, pattern formation and morphogenesis, based on the idea that the underlying mechanisms are respectively chemical or mechanical. Emerging quantitative approaches integrate the concepts and methods of complex sciences and emphasize the interplay between hierarchical levels of organization via mechano-chemical interactions. They draw upon recent improvements in mathematical and numerical morphogenetic models and upon considerable progress in collecting new quantitative data. This review highlights a variety of such models, which exhibit important advances, such as hybrid, stochastic and multiscale simulations.

Quantitative analysis of transient and sustained transforming growth factor-β signaling dynamics

Mathematical modeling and experimental analyses reveal that TGF-β ligand depletion has an important role in converting short-term graded signaling responses to long-term switch-like responses.

Cells respond in real time to the absolute number of TGF-β molecules in their environment.

A single pulse of TGF-β stimulation results in transient SMAD activation whereas repeated short pulses of stimulation result in sustained SMAD activation.

Ligand-induced short-term TGF-β/SMAD signaling activation is graded while long-term signaling response is switch-like or ultrasensitive.

TGF-β ligand depletion is a major cause of conversion from graded short-term responses to ultrasensitive long-term responses.

The transforming growth factor-β (TGF-β) pathway is a prominent signaling pathway that regulates diverse aspects of cellular homeostasis, including proliferation, differentiation, migration, and death (Massague, 1998). Remarkably, the pleiotropic biological effects of TGF-β are mediated by a relatively simple signaling module (Clarke and Liu, 2008). An interesting question is how such an apparently straightforward and simple cascade can generate a wide array of biological responses depending on the cellular context.

Members of the TGF-β superfamily are frequently used as morphogens in early embryo development (Green, 2002). The best-studied examples include Dpp in Drosophila and Activin in Xenopus (Gurdon and Bourillot, 2001; Lander, 2007). In the developmental context, cells can respond to a graded ligand concentration and produce discrete biological responses (e.g., transcription of certain genes, proliferation, or differentiation; Green, 2002). To convert continuous morphogen stimulation into discrete responses, mechanisms must exist to provide a threshold for the cellular response. How variable TGF-β ligand doses quantitatively control intracellular signaling dynamics and how continuous ligand doses are translated into discontinuous cellular fate decisions remains poorly understood.

We have previously reported that ligand molecules per cell is the input variable to which the cells respond, and ligand number per cell is the best predictor of signaling responses (Zi and Klipp, 2007a; Clarke et al, 2009). Here, we developed an improved mathematical model to predict TGF-β signaling responses by calibrating the model with various experimental data sets from different TGF-β stimulations. Using a combined experimental and mathematical modeling approach, we showed that TGF-β pulse stimulation results in transient activation of the pathway while repeated short pulses at short time intervals lead to a sustained activation similar to persistent ligand exposure.

We next investigate the system response to variable doses of TGF-β in HaCaT cells. Our mathematical model predicts that the short-term Smad2 phosphorylation (after 45 min of TGF-β stimulation) is a graded response, while long-term Smad2 activation (after 24 h of TGF-β stimulation) is a switch-like response (Figure 5A and B). As shown in Figure 5A–D, both short- and long-term Smad2 phosphorylation can be saturated but doses of TGF-β that cause maximum response are quite different. Additionally, the shapes of response curves were different. The short-term Smad2 activation was a graded (Michaelis–Menten-like) response with a very low apparent Hill coefficient of about 0.8 (Figure 5A and C) while the long-term Smad2 activation (P-Smad2 at 24 h) yielded a switch-like response with an apparent Hill coefficient of about 4.5 (Figure 5B and D). Thus, the Smad2 response is initially graded and sharpens over time to become ultrasensitive. To address whether TGF-β-inducible gene expression responses are graded or switch-like in the short and long term, we measured mRNA levels of Smad7, an early responsive gene of TGF-β and protein levels of p21 and PAI-1 whose inductions are delayed and late, respectively. The experimental data show that Smad7 induction exhibits a graded response with corresponding Hill coefficients of about 1.3 (Figure 5E), which is consistent with the graded P-Smad2 response at 45 min (Figure 5A and C). PAI-1 induction in response to variable doses of TGF-β for 24 h is highly ultrasensitive with an apparent Hill coefficient of ∼5.3. Compared with Smad7 and PAI-1, p21 induction is only modest ultrasensitive (n_Hill≈2) (Figure 5G). These results suggest short-term gene induction by TGF-β appears to be graded while long-term targets are more switch-like. Finally, we measured the growth inhibitory response of HaCaT cells to variable doses of TGF-β. The level of BrdU incorporation is also ultrasensitive with an apparent Hill coefficient of about 4.3 (Figure 5H). Therefore, the long-term TGF-β growth inhibitory response also shows a switch-like behavior. Finally, we show that TGF-β depletion affects long-term Smad phosphorylation and switch-like response of TGF-β signaling system. These findings shed new light on how continuous ligand doses are translated into discontinuous cell fate decisions in biological systems.

In summary, we have shown that the dose and time course of TGF-β stimulation have profound effects on Smad signaling dynamics. The rate of ligand depletion controls the duration of Smad2 phosphorylation. Cells can respond to a short pulse of TGF-β stimulation, and periodic short ligand exposures are sufficient to generate long-term signaling responses. Short-term TGF-β stimulation causes only transient pathway activation and can be terminated by ligand depletion. TGF-β-induced Smad2 phosphorylation is graded in the short-term but ultrasensitive (switch-like) in the long-term (Figure 7). Additionally, cell growth arrest in response to TGF-β shows switch-like rather than graded behavior. Our modeling and experimental analyses suggest that ligand depletion is likely to be involved in sharpening a graded response into a switch-like response.

Mammalian cells can decode the concentration of extracellular transforming growth factor-β (TGF-β) and transduce this cue into appropriate cell fate decisions. How variable TGF-β ligand doses quantitatively control intracellular signaling dynamics and how continuous ligand doses are translated into discontinuous cellular fate decisions remain poorly understood. Using a combined experimental and mathematical modeling approach, we discovered that cells respond differently to continuous and pulsating TGF-β stimulation. The TGF-β pathway elicits a transient signaling response to a single pulse of TGF-β stimulation, whereas it is capable of integrating repeated pulses of ligand stimulation at short time interval, resulting in sustained phospho-Smad2 and transcriptional responses. Additionally, the TGF-β pathway displays different sensitivities to ligand doses at different time scales. While ligand-induced short-term Smad2 phosphorylation is graded, long-term Smad2 phosphorylation is switch-like to a small change in TGF-β levels. Correspondingly, the short-term Smad7 gene expression is graded, while long-term PAI-1 gene expression is switch-like, as is the long-term growth inhibitory response. Our results suggest that long-term switch-like signaling responses in the TGF-β pathway might be critical for cell fate determination.

Patterned femtosecond-laser ablation of Xenopus laevis melanocytes for studies of cell migration, wound repair, and developmental processes

Ultrafast (femtosecond) lasers have become an important tool to investigate biological phenomena because of their ability to effect highly localized tissue removal in surgical applications. Here we describe programmable, microscale, femtosecond-laser ablation of melanocytes found on Xenopus laevis tadpoles, a technique that is applicable to biological studies in development, regeneration, and cancer research. We illustrate laser marking of individual melanocytes, and the drawing of patterns on melanocyte clusters to help track their migration and/or regeneration. We also demonstrate that this system can upgrade scratch tests, a technique used widely with cultured cells to study cell migration and wound healing, to the more realistic in vivo realm, by clearing a region of melanocytes and monitoring their return over time. In addition, we show how melanocyte ablation can be used for loss-of-function experiments by damaging neighboring tissue, using the example of abnormal tail regeneration following localized spinal cord damage. Since the size, shape, and depth of melanocytes vary as a function of tadpole age and melanocyte location (head or tail), an ablation threshold chart is given. Mechanisms of laser ablation are also discussed.

An atlas of gene regulatory networks reveals multiple three-gene mechanisms for interpreting morphogen gradients

Although >450 different topologies can achieve the same multicellular patterning function, they can be grouped into six main classes, which operate using different underlying dynamics.

Alternative designs for the same functions can therefore split into two types: (a) topology alterations that retain the same underlying dynamics and (b) alterations that utilize a completely different underlying dynamical mechanism.

This segregation of networks into distinct dynamical mechanisms can be revealed by the shape of the topology atlas itself.

Cell–cell communication is not usually part of the causal mechanism underlying a band-pass response during morphogen interpretation, but it can tune the result or increase robustness.

Understanding how gene regulatory networks (GRNs) achieve particular biological functions is a central question in systems biology. Systems biology promises to go beyond a case-by-case understanding of individual networks to map out the complete design space of mechanistic possibilities that underlie biological functions. Can such maps serve as useful theoretical frameworks in which to explore the general design principles for these functions? Towards addressing these questions, we created the first design space for a morphogen interpretation function.

In order to generate a design space for such a function, we enumerated all possible wiring designs of GRNs consisting of three genes and tested their ability to perform one particular morphogen interpretation function; stripe formation, as it represents a simplified form of the much studied French flag problem and is a commonly found gene expression pattern (Figure 1A). We found that only 5% of GRNs had the ability to generate a single stripe of gene expression when simulated with a fixed morphogen input in a one-dimensional model.

We hypothesized that the core mechanisms for producing the stripe of gene expression should be represented by topologies that contain only the necessary and sufficient gene–gene interactions for that function. Hence, we utilized the notions of complexity and neighborhood to generate a complexity atlas. GRNs of such an atlas (represented by nodes) are considered neighbors if they differ by a single gene–gene interaction (neighboring GRN nodes are connected by edges). Such a metagraph (graph of graphs) can then be reorganized using complexity (number of gene–gene interactions) to determine a GRNs position in the y axis, whereas GRNs are spaced in the x axis with the aim of reducing edge crossing (Figure 5A). This reorganization reveals a striking structure, where ‘stalactites' of complexity can be seen protruding from the bottom of the atlas. Each of these stalactites converges on a single ‘core' topology that by extensive analysis we find represents a distinct mechanism.

The mechanisms employ a diverse range of distinct space–time behaviors, and the underlying core topologies display design features such as modularity and feed-forward. We mapped the mechanisms to the complexity atlas by analyzing how each particular GRN of the atlas was working. The GRNs functioning via the different mechanisms are highlighted by the different colors in Figure 5A. Mechanisms thus occupy large regions of separated topology space, suggesting them to be discrete. Analyzing transitions between mechanisms through parameter space confirms this to be the case.

We find that three of the mechanisms are employed in real patterning systems, including both blastoderm patterning in Drosophila and mesoderm specification in Xenopus (Figure 5B). The remaining three mechanisms are thus candidates for employment in other patterning systems. We explored the performance features of these mechanisms, which suggest that some have features such as robustness to parameter variation that make them highly likely to be employed in particular patterning contexts.

Only one of the six-core mechanisms absolutely requires cell–cell communication for functionality, prompting us to predict that cell–cell communication will rarely be responsible for the basic dose response of morphogen interpretation networks. However, we show how cell–cell communication has an important role in robust stripe generation in the face of a noisy morphogen input and in fine tuning the quantitative details of stripe patterning.

In summary, the complexity atlas approach is an amendable approach to any system with a clear genotype–function relationship. We demonstrate how certain functions such as morphogen interpretation may have a range of potential solutions in contrast to previous studies that analyzed more constrained functions. Furthermore, we demonstrate how such an approach can be utilized to define a ‘design space' for a given biological function that describes the different mechanistic possibilities and how they relate to one another (Figure 5). Such a design space can be used practically as a guide to discern which patterning mechanisms are likely be at work in a particular context throwing up less intuitive possibilities with powerful performance features.

The interpretation of morphogen gradients is a pivotal concept in developmental biology, and several mechanisms have been proposed to explain how gene regulatory networks (GRNs) achieve concentration-dependent responses. However, the number of different mechanisms that may exist for cells to interpret morphogens, and the importance of design features such as feedback or local cell–cell communication, is unclear. A complete understanding of such systems will require going beyond a case-by-case analysis of real morphogen interpretation mechanisms and mapping out a complete GRN ‘design space.' Here, we generate a first atlas of design space for GRNs capable of patterning a homogeneous field of cells into discrete gene expression domains by interpreting a fixed morphogen gradient. We uncover multiple very distinct mechanisms distributed discretely across the atlas, thereby expanding the repertoire of morphogen interpretation network motifs. Analyzing this diverse collection of mechanisms also allows us to predict that local cell–cell communication will rarely be responsible for the basic dose-dependent response of morphogen interpretation networks.

Establishing positional information through gradient dynamics: a lesson from the Hedgehog signaling pathway

A long standing question in developmental biology is how morphogen gradients establish positional information during development. Although the existence of gradients and their role in developmental patterning is no longer in doubt, the ability of cells to respond to different morphogen concentrations has been controversial. In the Drosophila wing disc, Hedgehog (Hh) forms a concentration gradient along the anterior-posterior axis and establishes at least three different gene expression patterns. In a recent study, we challenged the prevailing idea that Hh establishes positional information in a dose-dependent manner and proposed a model in which dynamics of the gradient, resulting from the Hh gene network architecture, determines pattern formation in the wing disc. In this Extra View, we discuss further the methodology used in this study, highlight differences between this and other models of developmental patterning, and also present some questions that remain to be answered in this system.

APOBEC2, a selective inhibitor of TGFβ signaling, regulates left-right axis specification during early embryogenesis

The specification of left-right asymmetry is an evolutionarily conserved developmental process in vertebrates. The interplay between two TGFβ ligands, Derrière/GDF1 and Xnr1/Nodal, together with inhibitors such as Lefty and Coco/Cerl2, have been shown to provide the signals that lead to the establishment of laterality. However, molecular events leading to and following these signals remain mostly unknown. We find that APOBEC2, a member of the cytidine deaminase family of DNA/RNA editing enzymes, is induced by TGFβ signaling, and that its activity is necessary to specify the left-right axis in Xenopus and zebrafish embryos. Surprisingly, we find that APOBEC2 selectively inhibits Derrière, but not Xnr1, signaling. The inhibitory effect is conserved, as APOBEC2 blocks TGFβ signaling, and promotes muscle differentiation, in a mammalian myoblastic cell line. This demonstrates for the first time that a putative RNA/DNA editing enzyme regulates TGFβ signaling and plays a major role in development.

FGF8 acts as a classic diffusible morphogen to pattern the neocortex

Gain- and loss-of-function experiments have demonstrated that a source of fibroblast growth factor (FGF) 8 regulates anterior to posterior (A/P) patterning in the neocortical area map. Whether FGF8 controls patterning as a classic diffusible morphogen has not been directly tested. We report evidence that FGF8 diffuses through the mouse neocortical primordium from a discrete source in the anterior telencephalon, forms a protein gradient across the entire A/P extent of the primordium, and acts directly at a distance from its source to determine area identity. FGF8 immunofluorescence revealed FGF8 protein distributed in an A/P gradient. Fate-mapping experiments showed that outside the most anterior telencephalon, neocortical progenitor cells did not express Fgf8, nor were they derived from Fgf8-expressing cells, suggesting that graded distribution of FGF8 results from protein diffusion from the anterior source. Supporting this conclusion, a dominant-negative high-affinity FGF8 receptor captured endogenous FGF8 at a distance from the FGF8 source. New FGF8 sources introduced by electroporation showed haloes of FGF8 immunofluorescence indicative of FGF8 diffusion, and surrounding cells reacted to a new source of FGF8 by upregulating different FGF8-responsive genes in concentric domains around the source. Reducing endogenous FGF8 with the dominant-negative receptor in the central neocortical primordium induced cells to adopt a more posterior area identity, demonstrating long-range area patterning by FGF8. These observations support FGF8 as a classic diffusible morphogen in neocortex, thereby guiding future studies of neocortical pattern formation.

Predicting embryonic patterning using mutual entropy fitness and in silico evolution

During vertebrate embryogenesis, the expression of Hox genes that define anterior-posterior identity follows general rules: temporal colinearity and posterior prevalence. A mathematical measure for the quality or fitness of the embryonic pattern produced by a gene regulatory network is derived. Using this measure and in silico evolution we derive gene interaction networks for anterior-posterior (AP) patterning under two developmental paradigms. For patterning during growth (paradigm I), which is appropriate for vertebrates and short germ-band insects, the algorithm creates gene expression patterns reminiscent of Hox gene expression. The networks operate through a timer gene, the level of which measures developmental progression (a candidate is the widely conserved posterior morphogen Caudal). The timer gene provides a simple mechanism to coordinate patterning with growth rate. The timer, when expressed as a static spatial gradient, functions as a classical morphogen (paradigm II), providing a natural way to derive the AP patterning, as seen in long germ-band insects that express their Hox genes simultaneously, from the ancestral short germ-band system. Although the biochemistry of Hox regulation in higher vertebrates is complex, the actual spatiotemporal expression phenotype is not, and simple activation and repression by Hill functions suffices in our model. In silico evolution provides a quantitative demonstration that continuous positive selection can generate complex phenotypes from simple components by incremental evolution, as Darwin proposed.

Cell fate determination during tooth development and regeneration

Teeth arise from sequential and reciprocal interactions between the oral epithelium and the underlying cranial neural crest-derived mesenchyme. Their formation involves a precisely orchestrated series of molecular and morphogenetic events, and gives us the opportunity to discover and understand the nature of the signals that direct cell fates and patterning. For that reason, it is important to elucidate how signaling factors work together in a defined number of cells to generate the diverse and precise patterned structures of the mature functional teeth. Over the last decade, substantial research efforts have been directed toward elucidating the molecular mechanisms that control cell fate decisions during tooth development. These efforts have contributed toward the increased knowledge on dental stem cells, and observation of the molecular similarities that exist between tooth development and regeneration.

Particle tracking model of electrophoretic morphogen movement reveals stochastic dynamics of embryonic gradient

Some developmental events rely on an electrophoretic force to produce morphogenetic gradients. To quantitatively explore the dynamics of this process, we constructed a stochastic model of an early phase of left-right patterning: serotonin movement through the gap junction-coupled blastomeres of the Xenopus embryo. Particle-tracking simulations showed that a left-right gradient is formed rapidly, quickly reaching a final stable level. The voltage difference was critical for producing a morphogen gradient of the right steepness; gap junctional connectivity and morphogen mass determined the timing of the gradient. Endogenous electrophoresis drives approximately 50% of the particles across more than one cell width, and approximately 20% can travel across half the embryo. The stochastic behavior of the resulting gradients exhibited unexpected complexity among blastomeres' morphogen content, and showed how spatiotemporal variability within individual cells resulted in robust and consistent gradients across the embryonic left-right axis. Analysis of the distribution profile of gradient gain values made quantitative predictions about the conditions that result in the observed background level of laterality defects in unperturbed frog embryos. This work provides a general model that can be used to quantitatively analyze the unexpectedly complex dynamics of morphogens in a wide variety of systems.

Back to the future: evolution of computational models in plant morphogenesis

There has been a recent surge of studies in plant biology that combine experimental data with computational modeling. Here, we categorize a diversity of theoretical models and emphasize the need to tailor modeling approaches to the questions at hand. Models can start from biophysical or purely heuristic basic principles, and can focus at several levels of biological organization. Recent examples illustrate that this entire spectrum can be useful to understand plant development, and point to a future direction where more approaches are combined in fruitful ways--either by proving the same result with different basic principles or by exploring interactions across levels, in the so-called multilevel models.

Morphogen gradient formation

How morphogen gradients are formed in target tissues is a key question for understanding the mechanisms of morphological patterning. Here, we review different mechanisms of morphogen gradient formation from theoretical and experimental points of view. First, a simple, comprehensive overview of the underlying biophysical principles of several mechanisms of gradient formation is provided. We then discuss the advantages and limitations of different experimental approaches to gradient formation analysis.

Rab5-mediated endocytosis of activin is not required for gene activation or long-range signalling in Xenopus

Morphogen gradients provide positional cues for cell fate specification and tissue patterning during embryonic development. One important aspect of morphogen function, the mechanism by which long-range signalling occurs, is still poorly understood. In Xenopus, members of the TGF-beta family such as the nodal-related proteins and activin act as morphogens to induce mesoderm and endoderm. In an effort to understand the mechanisms and dynamics of morphogen gradient formation, we have used fluorescently labelled activin to study ligand distribution and Smad2/Smad4 bimolecular fluorescence complementation (BiFC) to analyse, in a quantitative manner, the cellular response to induction. Our results indicate that labelled activin travels exclusively through the extracellular space and that its range is influenced by numbers of type II activin receptors on responding cells. Inhibition of endocytosis, by means of a dominant-negative form of Rab5, blocks internalisation of labelled activin, but does not affect the ability of cells to respond to activin and does not significantly influence signalling range. Together, our data indicate that long-range signalling in the early Xenopus embryo, in contrast to some other developmental systems, occurs through extracellular movement of ligand. Signalling range is not regulated by endocytosis, but is influenced by numbers of cognate receptors on the surfaces of responding cells.

Bistability in a model of mesoderm and anterior mesendoderm specification in Xenopus laevis

In this paper we develop a model of mesendoderm specification in Xenopus laevis based on an existing gene regulation network. The mesendoderm is a population of cells that may contribute to either the mesoderm or endoderm. The model that we develop encompasses the time evolution of transcription factor concentrations in a single cell and is shown to have stable steady states that correspond to mesoderm and anterior mesendodermal cell types, but not endoderm (except in cells where Goosecoid expression is inhibited). Both in vitro and in vivo versions of the model are developed and analysed, the former indicating how cell fate is determined in large part by the concentration of Activin administered to a cell, with the model results comparing favourably with current quantitative experimental data. A numerical investigation of the in vivo model suggests that cell fate is determined largely by a VegT and beta-Catenin pre-pattern, subsequently being reinforced by Nodal. We argue that this sensitivity of the model to a VegT and beta-Catenin pre-pattern indicates that a key VegT self-limiting mechanism (for which there is experimental evidence) is absent from the model. Furthermore, we find that the lack of a steady state corresponding to endoderm is entirely consistent with current in vivo data, and that the in vivo model corresponds to mesendoderm specification on the dorsal, but not the ventral, side of the embryo.

A mathematical model of the stoichiometric control of Smad complex formation in TGF-beta signal transduction pathway

Cell fate in multicellular organism is regulated by the diffusible factor from surrounding cells in concentration-dependent manner. TGF-beta is a large protein family of the diffusible proteins secreted from a localized source. The signal of TGF-beta is transduced by Smad family transcription factor. Though it is well known that the stoichiometry of Smads in the transcriptional complex regulates the specificity of target genes of TGF-beta signal, little is known what the stoichiometry of Smads in the transcriptional complex is determined in TGF-beta signal transduction in concentration-dependent manner. To investigate the dynamics of Smad complex formation, we construct a two-compartment model for Smad complex formation in TGF-beta signal transduction. A simplified one-way oligomerization model, which ignores dissociation and well appropriate the full model under high expression levels of R- and Co-Smad, is constructed to analytically investigate the effect of the oligomerization of Smad. Our one-way model reveals that not only shuttling of the Smad from the cytoplasm to the nucleus but also the preferential accumulation of the heteromeric complex in oligomerization can contribute to the predominant production of the heteromeric complex of Smad including both R- and Co-Smad. It is also shown that oligomerization of Smad can contribute to the specificity of signal transduction. In endothelial cells, both Smad-1/5/8 and -2/3 pathways are activated by TGF-beta. The difference of the activity between the two pathways is amplified by trimerization but not by dimerization, suggesting possible importance of trimerization in maintaining the specificity of signal transduction.

Regulative feedback in pattern formation: towards a general relativistic theory of positional information

Positional specification by morphogen gradients is traditionally viewed as a two-step process. A gradient is formed and then interpreted, providing a spatial metric independent of the target tissue, similar to the concept of space in classical mechanics. However, the formation and interpretation of gradients are coupled, dynamic processes. We introduce a conceptual framework for positional specification in which cellular activity feeds back on positional information encoded by gradients, analogous to the feedback between mass-energy distribution and the geometry of space-time in Einstein's general theory of relativity. We discuss how such general relativistic positional information (GRPI) can guide systems-level approaches to pattern formation.

Decoding the quantitative nature of TGF-beta/Smad signaling

How transforming growth factor-beta (TGF-beta) signaling elicits diverse cell responses remains elusive, despite the major molecular components of the pathway being known. We contend that understanding TGF-beta biology requires mathematical models to decipher the quantitative nature of TGF-beta/Smad signaling and to account for its complexity. Here, we review mathematical models of TGF-beta superfamily signaling that predict how robustness is achieved in bone-morphogenetic-protein signaling in the Drosophila embryo, how changes in receptor-trafficking dynamics can be exploited by cancer cells and how the basic mechanisms of TGF-beta/Smad signaling conspire to promote Smad accumulation in the nucleus. These studies demonstrate the power of mathematical modeling for understanding TGF-beta biology.

Theoretical and experimental approaches to understand morphogen gradients

Morphogen gradients, which specify different fates for cells in a direct concentration-dependent manner, are a highly influential framework in which pattern formation processes in developmental biology can be characterized. A common analysis approach is combining experimental and theoretical strategies, thereby fostering relevant data on the dynamics and transduction of gradients. The mechanisms of morphogen transport and conversion from graded information to binary responses are some of the topics on which these combined strategies have shed light. Herein, we review these data, emphasizing, on the one hand, how theoretical approaches have been helpful and, on the other hand, how these have been combined with experimental strategies. In addition, we discuss those cases in which gradient formation and gradient interpretation at the molecular and/or cellular level may influence each other within a mutual feedback loop. To understand this interplay and the features it yields, it becomes essential to take system-level approaches that combine experimental and theoretical strategies.

Deriving structure from evolution: metazoan segmentation

Segmentation is a common feature of disparate clades of metazoans, and its evolution is a central problem of evolutionary developmental biology. We evolved in silico regulatory networks by a mutation/selection process that just rewards the number of segment boundaries. For segmentation controlled by a static gradient, as in long-germ band insects, a cascade of adjacent repressors reminiscent of gap genes evolves. For sequential segmentation controlled by a moving gradient, similar to vertebrate somitogenesis, we invariably observe a very constrained evolutionary path or funnel. The evolved state is a cell autonomous 'clock and wavefront' model, with the new attribute of a separate bistable system driven by an autonomous clock. Early stages in the evolution of both modes of segmentation are functionally similar, and simulations suggest a possible path for their interconversion. Our computation illustrates how complex traits can evolve by the incremental addition of new functions on top of pre-existing traits.

Wnt and TGF-beta expression in the sponge Amphimedon queenslandica and the origin of metazoan embryonic patterning

BACKGROUND

The origin of metazoan development and differentiation was contingent upon the evolution of cell adhesion, communication and cooperation mechanisms. While components of many of the major cell signalling pathways have been identified in a range of sponges (phylum Porifera), their roles in development have not been investigated and remain largely unknown. Here, we take the first steps toward reconstructing the developmental signalling systems used in the last common ancestor to living sponges and eumetazoans by studying the expression of genes encoding Wnt and TGF-beta signalling ligands during the embryonic development of a sponge.

METHODOLOGY/PRINCIPAL FINDINGS

Using resources generated in the recent sponge Amphimedon queenslandica (Demospongiae) genome project, we have recovered genes encoding Wnt and TGF-beta signalling ligands that are critical in patterning metazoan embryos. Both genes are expressed from the earliest stages of Amphimedon embryonic development in highly dynamic patterns. At the time when the Amphimedon embryos begin to display anterior-posterior polarity, Wnt expression becomes localised to the posterior pole and this expression continues until the swimming larva stage. In contrast, TGF-beta expression is highest at the anterior pole. As in complex animals, sponge Wnt and TGF-beta expression patterns intersect later in development during the patterning of a sub-community of cells that form a simple tissue-like structure, the pigment ring. Throughout development, Wnt and TGF-beta are expressed radially along the anterior-posterior axis.

CONCLUSIONS/SIGNIFICANCE

We infer from the expression of Wnt and TGF-beta in Amphimedon that the ancestor that gave rise to sponges, cnidarians and bilaterians had already evolved the capacity to direct the formation of relatively sophisticated body plans, with axes and tissues. The radially symmetrical expression patterns of Wnt and TGF-beta along the anterior-posterior axis of sponge embryos and larvae suggest that these signalling pathways contributed to establishing axial polarity in the very first metazoans.

Tsukushi modulates Xnr2, FGF and BMP signaling: regulation of Xenopus germ layer formation

Background

Cell-cell communication is essential in tissue patterning. In early amphibian development, mesoderm is formed in the blastula-stage embryo through inductive interactions in which vegetal cells act on overlying equatorial cells. Members of the TGF-β family such as activin B, Vg1, derrière and Xenopus nodal-related proteins (Xnrs) are candidate mesoderm inducing factors, with further activity to induce endoderm of the vegetal region. TGF-β-like ligands, including BMP, are also responsible for patterning of germ layers. In addition, FGF signaling is essential for mesoderm formation whereas FGF signal inhibition has been implicated in endoderm induction. Clearly, several signaling pathways are coordinated to produce an appropriate developmental output; although intracellular crosstalk is known to integrate multiple pathways, relatively little is known about extracellular coordination.

Methodology/Principal Findings

Here, we show that Xenopus Tsukushi (X-TSK), a member of the secreted small leucine rich repeat proteoglycan (SLRP) family, is expressed in ectoderm, endoderm, and the organizer during early development. We have previously reported that X-TSK binds to and inhibits BMP signaling in cooperation with chordin. We now demonstrate two novel interactions: X-TSK binds to and inhibits signaling by FGF8b, in addition to binding to and enhancement of Xnr2 signaling. This signal integration by X-TSK at the extracellular level has an important role in germ layer formation and patterning. Vegetally localized X-TSK potentiates endoderm formation through coordination of BMP, FGF and Xnr2 signaling. In contrast, X-TSK inhibition of FGF-MAPK signaling blocks ventrolateral mesoderm formation, while BMP inhibition enhances organizer formation. These actions of X-TSK are reliant upon its expression in endoderm and dorsal mesoderm, with relative exclusion from ventrolateral mesoderm, in a pattern shaped by FGF signals.

Conclusions/Significance

Based on our observations, we propose a novel mechanism by which X-TSK refines the field of positional information by integration of multiple pathways in the extracellular space.

Computational analysis of BMP gradients in dorsal-ventral patterning of the zebrafish embryo

The genetic network controlling early dorsal-ventral (DV) patterning has been extensively studied and modeled in the fruit fly Drosophila. This patterning is driven by signals coming from bone morphogenetic proteins (BMPs), and regulated by interactions of BMPs with secreted factors such as the antagonist short gastrulation (Sog). Experimental studies suggest that the DV patterning of vertebrates is controlled by a similar network of BMPs and antagonists (such as Chordin, a homologue of Sog), but differences exist in how the two systems are organized, and a quantitative comparison of pattern formation in them has not been made. Here, we develop a computational model in three dimensions of the zebrafish embryo and use it to study molecular interactions in the formation of BMP morphogen gradients in early DV patterning. Simulation results are presented on the dynamics BMP gradient formation, the cooperative action of two feedback loops from BMP signaling to BMP and Chordin synthesis, and pattern sensitivity with respect to BMP and Chordin dosage. Computational analysis shows that, unlike the case in Drosophila, synergy of the two feedback loops in the zygotic control of BMP and Chordin expression, along with early initiation of localized Chordin expression, is critical for establishment and maintenance of a stable and appropriate BMP gradient in the zebrafish embryo.

The competence of Xenopus blastomeres to produce neural and retinal progeny is repressed by two endo-mesoderm promoting pathways

Only a subset of cleavage stage blastomeres in the Xenopus embryo is competent to contribute cells to the retina; ventral vegetal blastomeres do not form retina even when provided with neuralizing factors or transplanted to the most retinogenic position of the embryo. These results suggest that endogenous maternal factors in the vegetal region repress the ability of blastomeres to form retina. Herein we provide three lines of evidence that two vegetal-enriched maternal factors (VegT, Vg1), which are known to promote endo-mesodermal fates, negatively regulate which cells are competent to express anterior neural and retinal fates. First, both molecules can repress the ability of dorsal-animal retinogenic blastomeres to form retina, converting the lineage from neural/retinal to non-neural ectodermal and endo-mesodermal fates. Second, reducing the endogenous levels of either factor in dorsal-animal retinogenic blastomeres expands expression of neural/retinal genes and enlarges the retina. The dorsal-animal repression of neural/retinal fates by VegT and Vg1 is likely mediated by Sox17alpha and Derriere but not by XNr1. VegT and Vg1 likely exert their effects on neural/retinal fates through at least partially independent pathways because Notch1 can reverse the effects of VegT and Derriere but not those of Vg1 or XNr1. Third, reduction of endogenous VegT and/or Vg1 in ventral vegetal blastomeres can induce a neural fate, but only allows expression of a retinal fate when both BMP and Wnt signaling pathways are concomitantly repressed.

CHD4/Mi-2beta activity is required for the positioning of the mesoderm/neuroectoderm boundary in Xenopus

Experiments in Xenopus have illustrated the importance of extracellular morphogens for embryonic gene regulation in vertebrates. Much less is known about how induction leads to the correct positioning of boundaries; for example, between germ layers. Here we report that the neuroectoderm/mesoderm boundary is controlled by the chromatin remodeling ATPase CHD4/Mi-2beta. Gain and loss of CHD4 function experiments shifted this boundary along the animal-vegetal axis at gastrulation, leading to excess mesoderm formation at the expense of neuroectoderm, or vice versa. This phenotype results from specific alterations in gene transcription, notably of the neural-promoting gene Sip1 and the mesodermal regulatory gene Xbra. We show that CHD4 suppresses Sip1 transcription by direct binding to the 5' end of the Sip1 gene body. Furthermore, we demonstrate that CHD4 and Sip1 expression levels determine the "ON" threshold for Nodal-dependent but not for eFGF-dependent induction of Xbra transcription. The CHD4/Sip1 epistasis thus constitutes a regulatory module, which balances mesoderm and neuroectoderm formation.

Two-dimensional morphogen gradient in Xenopus: boundary formation and real-time transduction response

Morphogen gradients play an important role in pattern formation in embryo. However, the interpretation of position in a morphogen gradient is not well understood. Because it is hard to analyze morphogen gradients especially in opaque embryos such as those of Xenopus, it is necessary to fix and section the embryo, thereby eliminating the possibility of real-time observation, and making more difficult the interpretation of events that take place in three dimensions. We describe here a two-dimensional preparation of cells from a Xenopus blastula animal cap, in which an activin concentration gradient appears to be formed and interpreted at the same rate and in the same way as in normal embryos. We use two-dimensional preparations of this kind to contribute the following new information about gradient formation and interpretation in embryo. We determine the dynamics of formation of an activin activity gradient in real time. We demonstrate that this gradient is established by diffusion of activin through intercellular space and does not require internalization of receptor or ligand. We also show that the generation of a boundary of gene expression depends on the interpretation, rather than a change of composition, of the concentration gradient.

Morpheus Unbound: Reimagining the Morphogen Gradient

The theory that the spatial organization of cell fate is orchestrated by gradients of diffusing molecules was a major contribution to 20th century developmental biology. Although the existence of morphogens is no longer in doubt, studies on the formation and function of their gradients have yielded far more puzzles than answers. On close inspection, every morphogen gradient seems to use a rich array of regulatory mechanisms, suggesting that the tasks carried out by such systems are far more extensive than previously thought.

21st century neontology and the comparative development of the vertebrate skull

Classic neontology (comparative embryology and anatomy), through the application of the concept of homology, has demonstrated that the development of the gnathostome (jawed vertebrate) skull is characterized both by a fidelity to the gnathostome bauplan and the exquisite elaboration of final structural design. Just as homology is an old concept amended for modern purposes, so are many of the questions regarding the development of the skull. With due deference to Geoffroy-St. Hilaire, Cuvier, Owen, Lankester et al., we are still asking: How are bauplan fidelity and elaboration of design maintained, coordinated, and modified to generate the amazing diversity seen in cranial morphologies? What establishes and maintains pattern in the skull? Are there universal developmental mechanisms underlying gnathostome autapomorphic structural traits? Can we detect and identify the etiologies of heterotopic (change in the topology of a developmental event), heterochronic (change in the timing of a developmental event), and heterofacient (change in the active capacetence, or the elaboration of capacity, of a developmental event) changes in craniofacial development within and between taxa? To address whether jaws are all made in a like manner (and if not, then how not), one needs a starting point for the sake of comparison. To this end, we present here a "hinge and caps" model that places the articulation, and subsequently the polarity and modularity, of the upper and lower jaws in the context of cranial neural crest competence to respond to positionally located epithelial signals. This model expands on an evolving model of polarity within the mandibular arch and seeks to explain a developmental patterning system that apparently keeps gnathostome jaws in functional registration yet tractable to potential changes in functional demands over time. It relies upon a system for the establishment of positional information where pattern and placement of the "hinge" is driven by factors common to the junction of the maxillary and mandibular branches of the first arch and of the "caps" by the signals emanating from the distal-most first arch midline and the lamboidal junction (where the maxillary branch meets the frontonasal processes). In this particular model, the functional registration of jaws is achieved by the integration of "hinge" and "caps" signaling, with the "caps" sharing at some critical level a developmental history that potentiates their own coordination. We examine the evidential foundation for this model in mice, examine the robustness with which it can be applied to other taxa, and examine potential proximate sources of the signaling centers. Lastly, as developmental biologists have long held that the anterior-most mesendoderm (anterior archenteron roof or prechordal plate) is in some way integral to the normal formation of the head, including the cranial skeletal midlines, we review evidence that the seminal patterning influences on the early anterior ectoderm extend well beyond the neural plate and are just as important to establishing pattern within the cephalic ectoderm, in particular for the "caps" that will yield medial signaling centers known to coordinate jaw development.

Apoptosis is required during early stages of tail regeneration in Xenopus laevis

The Xenopus tadpole is able to regenerate its tail, including skin, muscle, notochord, spinal cord and neurons and blood vessels. This process requires rapid tissue growth and morphogenesis. Here we show that a focus of apoptotic cells appears in the regeneration bud within 12 h of amputation. Surprisingly, when caspase-3 activity is specifically inhibited, regeneration is abolished. This is true of tails both before and after the refractory period. Programmed cell death is only required during the first 24 h after amputation, as later inhibition has no effect on regeneration. Inhibition of caspase-dependent apoptosis results in a failure to induce proliferation in the growth zone, a mispatterning of axons in the regenerate, and the appearance of ectopic otoliths in the neural tube, in the context of otherwise normal continued development of the larva. Larvae amputated during the refractory stage exhibit a much broader domain of caspase-3-positive cells, suggesting a window for the amount of apoptosis that is compatible with normal regeneration. These data reveal novel roles for apoptosis in development and indicate that a degree of apoptosis is an early and obligate component of normal tail regeneration, suggesting the possibility of the existence of endogenous inhibitory cells that must be destroyed by programmed cell death for regeneration to occur.

Cell lineage transport: a mechanism for molecular gradient formation

Gradient formation is a fundamental patterning mechanism during embryo development, commonly related to secreted proteins that move along an existing field of cells. Here, we mathematically address the feasibility of gradients of mRNAs and non-secreted proteins. We show that these gradients can arise in growing tissues whereby cells dilute and transport their molecular content as they divide and grow, a mechanism we termed ‘cell lineage transport.' We provide an experimental test by unveiling a distal-to-proximal gradient of Hoxd13 in the vertebrate developing limb bud driven by cell lineage transport, corroborating our model. Our study indicates that gradients of non-secreted molecules exhibit a power-law profile and can arise for a wide range of biologically relevant parameter values. Dilution and nonlinear growth confer robustness to the spatial gradient under changes in the cell cycle period, but at the expense of sensitivity in the timing of gradient formation. We expect that gradient formation driven by cell lineage transport will provide future insights into understanding the coordination between growth and patterning during embryonic development.

Dual function of the Drosophila Alk1/Alk2 ortholog Saxophone shapes the Bmp activity gradient in the wing imaginal disc

Wing patterning in Drosophila requires a Bmp activity gradient created by two Bmp ligands, Gbb and Dpp, and two Bmp type I receptors, Sax and Tkv. Gbb provides long-range signaling, while Dpp signals preferentially to cells near its source along the anteroposterior (AP) boundary of the wing disc. How each receptor contributes to the signaling activity of each ligand is not well understood. Here, we show that while Tkv mediates signals from both Dpp and Gbb, Sax exhibits a novel function for a Bmp type I receptor: the ability to both promote and antagonize signaling. Given its high affinity for Gbb, this dual function of Sax impacts the function of Gbb in the Bmp activity gradient more profoundly than does Dpp. We propose that this dual function of Sax is dependent on its receptor partner. When complexed with Tkv, Sax facilitates Bmp signaling, but when alone, Sax fails to signal effectively and sequesters Gbb. Overall, our model proposes that the balance between antagonizing and promoting Bmp signaling varies across the wing pouch, modulating the level and effective range, and, thus, shaping the Bmp activity gradient. This previously unknown mechanism for modulating ligand availability and range raises important questions regarding the function of vertebrate Sax orthologs.