- Open Access
Mathematical models for the Notch and Wnt signaling pathways and the crosstalk between them during somitogenesis
© Wang et al.; licensee BioMed Central Ltd. 2013
- Received: 5 February 2013
- Accepted: 15 April 2013
- Published: 20 April 2013
Somitogenesis is a fundamental characteristic feature of development in various animal embryos. Molecular evidence has proved that the Notch and Wnt pathways play important roles in regulating the process of somitogenesis and there is crosstalk between these two pathways. However, it is difficult to investigate the detailed mechanism of these two pathways and their interactions in somitogenesis through biological experiments. In recent years some mathematical models have been proposed for the purpose of studying the dynamics of the Notch and Wnt pathways in somitogenesis. Unfortunately, only a few of these models have explored the interactions between them.
In this study, we have proposed three mathematical models for the Notch signalling pathway alone, the Wnt signalling pathway alone, and the interactions between them. These models can simulate the dynamics of the Notch and Wnt pathways in somitogenesis, and are capable of reproducing the observations derived from wet experiments. They were used to investigate the molecular mechanisms of the Notch and Wnt pathways and their crosstalk in somitogenesis through the model simulations.
Three mathematical models are proposed for the Notch and Wnt pathways and their interaction during somitogenesis. The simulations demonstrate that the extracellular Notch and Wnt signals are essential for the oscillating expressions of both Notch and Wnt target genes. Moreover, the internal negative feedback loops and the three levels of crosstalk between these pathways play important but distinct roles in maintaining the system oscillation. In addition, the results of the parameter sensitivity analysis of the models indicate that the Notch pathway is more sensitive to perturbation in somitogenesis.
- Negative Feedback Loop
- Notch Pathway
- Parameter Sensitivity Analysis
- Notch Target Gene
- Hes7 Gene
Animals have a segmented aspect of the body axis and somitogenesis has long been thought to be a key aspect of the basic design of animals. In the early developing animal embryo the body is organized in a series of embryonic tissue masses called somites . Somites are progressively pinched off in pairs from the anterior end of two rods of mesenchymal tissue called the presomitic mesoderm (PSM) . It is accepted that somite formation is controlled by a complicated gene network named the segmentation clock. For nearly a decade it has been known that the Notch pathway, the Wnt pathway and the fibroblast growth factor (FGF) pathway are the important components of the segmentation clock . In particular, the Notch and Wnt pathways regulate the oscillating expressions of their target genes, which play major roles in controlling somite formation [4–6]. In recent years mathematical models have been proposed to reveal the mechanisms of the two pathways and their crosstalk in the process of somitogenesis. In 2003 Julian Lewis et al. proposed a simple mathematical model of the Notch pathway in zebrafish somitogenesis . They modeled the oscillating expressions of Notch pathway target genes by introducing two feedback loops. In 2009, Smita Agrawal et al. proposed a model of the Notch pathway during somitogenesis to elucidate the mechanisms of context-dependent signaling of the Notch pathway . They modeled bistability in Notch signaling. In 2010 Alan J. Terry et al. proposed a spatio-temporal model of Notch signaling in the zebrafish segmentation clock . They adopted a spatially-explicit modeling approach that can display intracellular protein diffusion graphically. In the same year, Peter B. Jensen et al. proposed a mathematical model to capture the oscillation of the Wnt pathway in somitogenesis . The core of their model was a negative feedback loop centered on Axin2. Now, more and more evidence supports the view that somite formation relies on complex cooperation among multiple signaling pathways. In 2007, J.G. Rodríguez-González et al. proposed a mathematical model to investigate the interaction between the Notch and the Wnt pathways in the segmentation clock in mice . In 2008, Albert Goldbeter et al. proposed a theoretical model for understanding the mechanism of interactions among the Notch, Wnt and FGF pathways . Moirés Santillán et al. also proposed a mathematical model for the gene regulatory network of the mouse embryo to elucidate somite formation . In 2009, A. Kazama et al. proposed a mathematical model to reveal the interaction of the Notch and Wnt pathways in the segmentation clock . Although the simulation results of these models agree well with some results of biological experiments, they only considered simple interaction relationships between the two pathways. More accurate and complicated mathematical models are still needed to further our understanding of the detailed mechanism of the Notch and Wnt pathways and their crosstalk in somitogenesis.
In this study, we have proposed three more complicated mathematical models for the Notch and Wnt signaling pathways and their crosstalk in somitogenesis, taking the mouse as example. In modeling the Notch and Wnt signaling pathways in isolation, three core negative feedback loops centered on Lfng, Hes7 and Axin2, respectively, were considered, while in the combined model of the two pathways, three levels of cross-regulation were modeled. These models not only simulate the periodic expressions of the Notch and Wnt target genes in somitogenesis, but also reproduce the wet experimental results in the literature. The simulations demonstrate that the extracellular Notch and Wnt signals are essential for the oscillating expressions of both Notch and Wnt target genes. Moreover, the internal negative feedback loops and the three levels of crosstalk between the Notch and Wnt pathways play important but distinct roles in maintaining the system oscillation. In addition, the results of the parameter sensitivity analysis of the models indicate that the Notch pathway is more sensitive to perturbation in somitogenesis.
The model for the Notch pathway in isolation
Notch-mediated signaling is initiated via the binding of the delta-like 1 (Dll1) ligand to the Notch receptor. Then the intracellular domain of Notch (NICD) is cleaved from the membrane tether. It is transported into the nucleus and associates with the recombining binding protein (RBP-j) to form a transcriptional activator that activates the transcription of a set of target genes, including the Lunatic Fringe (Lfng) and the hairy and enhancer of split 7 (Hes7) genes. The Lfng and Hes7 mRNAs are transported from nucleus and are translated into proteins in the cytoplasm. Lfng inhibits the cleavage of NICD from Notch leading to repression of the transcription factor NICD/ RBP-j, so a negative feedback loop in the Notch pathway, which is termed “the big feedback loop”, is formed . Hes7 inhibits the transcription of both itself and the Lfng gene, and thus another negative feedback loop of the Notch pathway, termed “the small feedback loop”, is formed . The periodic expressions of Lfng and Hes7 genes are essential for somite formation .
Model simulation for the Notch pathway in isolation
We used this model to perform simulations and tried to reveal the molecular mechanism behind the phenomena. First, we investigated the influence of the upstream Notch signals on the expressions of the target genes. When knocking out the Dll1 gene, the ligand of the Notch pathway, at time point 120 minutes, we found the expressions of the Notch target genes do not oscillate and the expression levels descend markedly (see Figure 2 (B)). This suggests that Dll1 is essential for the normal expressions of Notch target genes. When knocking out the Notch gene, receptor of the Notch pathway, at time point 120 minutes, the oscillating expressions of the target genes disappear as in the knockout of Dll1 (see Figure 2 (B)). Moreover, knockout of Dll1 or Notch made the expressions of NICD and the NICD-RBP-j transcriptional activator disappear (see Figure 2 (C)). NICD is the direct regulator of the Notch pathway target genes, so when it disappeared the oscillating expressions of the target genes were destroyed. The results demonstrate that the activity of upstream Notch signals is necessary for the oscillating expressions of the Notch pathway target genes.
Next, we investigated the influence of the feedback loops on the oscillating expressions of the Notch pathway target genes. After knocking out the Lfng gene at time point 120 minutes, we found that expression of Hes7 gene still oscillated, though its maximum expression level increased a little (see Figure 2 (D)). This suggests that the big negative feedback loop formed by Lfng is not essential for the oscillating expressions of the Notch pathway target genes. Similarly, we knocked out the Hes7 gene after one period. It was found that the expression of the Lfng gene increased markedly and its oscillating expression was destroyed (see Figure 2 (E)). This suggests that the small negative feedback loop formed by Hes7 is necessary for the oscillating expressions of the Notch pathway target genes. NICD in the big feedback loop induced the expressions of the target genes, whereas Hes7 in the small feedback loop inhibited them, so the oscillating pattern of the target genes was in phase with NICD but in antiphase with Hes7 (see Figure 2 (F)). The phase relationship of these two feedback loops in the Notch pathway is crucial for the correctly oscillating expressions of the target genes. After knocking out the Lfng gene, NICD increased quickly but the NICD-RBP-j transcriptional activator only increased a little owing to the constant concentration of RBP-j (see Figure 2 (G)). The expression of Hes7 then increased a little with the NICD-RBP-j transcriptional activator, but the concentration increase of NICD did not influence the oscillating expression of Hes7 (see Figure 2 (G)). On the other hand, after Hes7 is knocked out, NICD and the NICD-RBP-j transcriptional activator did not change immediately, but decreased to nearly zero after about one period (see Figure 2 (H)). Obviously, the increase of Lfng gene expression is not due to the concentration increase of its transcriptional activator (the NICD-RBP-j complex) but to the concentration decrease of its inhibitor (Hes7). The delayed concentration decrease of NICD and the NICD and RBP-j transcriptional activator was mainly due to a large increase in Lfng after Hes7 was knocked out, which precluded the cleavage of NICD from Notch. The influence of Hes7 on upstream Notch signals is not instantaneous but follows a time delay. Therefore, although the two negative feedback loops are all-important for somitogenesis, the small feedback loop is closely related to the periodic expressions of the Notch pathway target genes, while the big feedback loop is complementary to the periodic expressions of those genes. However, the big feedback loop is essential for the oscillating expressions of the Wnt pathway target genes in the combined model (see the section describing the combined model). Research of Ferjentsik et al. indicated that the big feedback loop is important for the formation of the somite in mouse embryo development . In their study, after the Lfng gene of the mouse embryo was knocked out, Notch activity was still dynamic but the somite was irregular in these embryos.
In summary, the simulation results demonstrate that the model for the Notch signaling pathway in isolation is capable of simulating the oscillating expressions of the Notch pathway target genes, and can also reproduce the wet experimental results. It has the potential for further use in research on the molecular mechanism of somitogenesis.
The model for the Wnt signaling pathway in isolation
A total of 13 ODEs for the Wnt signaling model and their biological explanations are given in Additional file 1. We assumed the activation of Dsh is reversible and obeys Michaelis-Menten kinetics. The catalysis of Dsh activation by Wnt was modeled using a Hill equation with Hill coefficient 1 (Eq 2.1 in Additional file 1). The transcription of the Dll1 gene in the nucleus when the β-catenin-Lef1 complex is activated was modeled using a Hill equation with Hill coefficient 2 (Eq 2.13 in Additional file 1), and the translation of the Dll1 mRNA was modeled using a mass action equation (Eq 2.7 in Additional file 1). The negative feedback loop centered on Axin2 was modeled using six major reactions. A mass action equation was used to model the reversible binding of GSK3 to Axin2 to form the degradation complex that phosphorylates β-catenin (Eq 2.3 in Additional file 1). The reversible phosphorylation of β-catenin when the GSK3-Axin2 complex acts as catalyst was modeled using a Michaelis-Menten equation with the catalysis rate proportional to the GSK3 of the GSK3-Axin2 complex in the total GSK3 (Eq 2.6 in Additional file 1). The two above reactions guarantee a very low concentration of the unphosphorylated β-catenin in the cytoplasm so as not to initiate the expressions of the Wnt pathway target genes, including the Axin2 gene. The binding of active Dsh to Axin2 and degradation of Axin2 were modeled using a mass action equation (Eq 2.2 in Additional file 1). This reaction can destroy the degradation complex by degrading Axin2, so that the unphosphorylated β-catenin protein can enter the nucleus to activate target genes. The binding of unphosphorylated β-catenin to Lef1 in the nucleus to form the transcription activator was modeled using a mass action law (Eq 2.11 in the Additional file 1). The transcription of the Axin2 gene when the β-catenin-Lef1 complex is activated was modeled using a Hill equation with Hill coefficient 2 because the complex binds to the Axin2 gene at 2 sites (Eq 2.12 in Additional file 1). The sixth reaction is the translation of the Axin2 mRNA, which was modeled using a mass action equation on the assumption that the translation rate is proportional to the Axin2 mRNA concentration (Eq 2.4 in Additional file 1).
Model simulation for the Wnt pathway in isolation
In summary, the simulation results demonstrate that the model for the Wnt signaling pathway in isolation is capable of simulating the oscillating expressions of the Wnt pathway target genes, and also reproduces the wet experimental results. It has the potential for use in further research on the molecular mechanism of somitogenesis.
The combined model for crosstalk between the Notch and Wnt pathways
The crosstalk between the Notch and Wnt pathways in somitogenesis is very complicated. Some biological experiments have confirmed several different levels of crosstalk. The first level is through Dll1, which is the ligand of Notch pathway and is also a downstream target gene of the Wnt pathway. Research by Juan Galceran et al.  confirmed that Lef1 binds to β-catenin in the nucleus to form the Lef1-β-catenin complex that activates the expression of endogenous Dll1 in fibroblasts. Research by Michael Hofmann et al.  indicated that the Lef/Tcf factors, which are regulators of the Wnt pathway, activate transcription of Dll1 cooperating with TBX6. The second level is through NICD in the Notch pathway and Dsh in the Wnt pathway, respectively. Studies by Alexander Aulehla et al.  indicated that NICD binds to active Dsh and thus inhibits the Wnt pathway; meanwhile, Notch signals are also inhibited. Moreover, the oscillating expressions of the Notch pathway target genes are in antiphase with the Wnt pathway target genes, so it is conjectured that the two pathways inhibit each other alternately. There is molecular evidence supporting this viewpoint. In Drosophila, it was found that NICD binds to the PDZ domain of Dsh, while Dsh interacts antagonistically with NICD . It is likely that Dsh blocks Notch signaling directly through binding of NICD . JG Rodríguez-González et al.  established a mathematical model to simulate this level of crosstalk between the two pathways. Their model simulates the results of Alexander Aulehla et al. . The third level is through naked cuticle 1 (Nkd1) protein, which is encoded by a Wnt downstream gene. Research by Aki Ishikawa et al.  indicated that the transcription of Nkd1 was extremely downregulated in the PSM of vestigial tail (vt/vt), a hypomorphic mutant of Wnt3a, so it can be speculated that the Nkd1 gene is regulated by Wnt signals. Moreover, Nkd1 oscillation had a similar phase to Lfng transcription and they were arrested in the Hes7 mutation embryo, hence it is likely that Nkd1 is regulated by Hes7. Furthermore, Nkd1 inhibits Wnt signals by binding to Dsh.
Model simulation for crosstalk between the Notch and Wnt pathways
From our simulation experiments, we are able to draw some interesting conclusions: The big feedback loop centered on Lfng in the Notch pathway is complementary to the periodic expressions of the Notch target genes, but is essential for the oscillating expressions of the Wnt target genes in the combined model. In contrast, the small feedback loop centered on Hes7 in the Notch pathway is essential for the oscillating expressions of the Notch target genes, but is not related to the oscillating expressions of the Wnt target genes in the combined model.
Lastly, we investigated the regulatory effect of Nkd1 on the crosstalk between the Notch and Wnt pathways. Aki Ishikawa et al.  indicated that Nkd1 is a downstream gene of the Wnt pathway, but it is also regulated by Hes7, a downstream gene of the Notch pathway. Nkd1 was expressed in the same phase as Lfng in the middle PSM . Yan et al.  and KA Wharton et al.  found that Nkd1 inhibits the Wnt pathway by binding with Dsh. From Figure 6 (A) we can see that expression of the Nkd1 gene was oscillating in the same phase as Lfng and Hes7 under normal oscillation conditions. This is in accord with the biological facts. After Wnt signals were knocked out at time point 240 minutes, the expression level of the Nkd1 gene decreased, but still oscillated (see Figure 7 (E)). This suggests that the Nkd1 gene is enhanced by Wnt signals. Wnt signals regulate the expression of the Dll1 gene, which is the ligand of the Notch pathway, so when Wnt signals are removed the Notch pathway target genes will be influenced. In order to see whether or not the regulation of the Nkd1 gene by Wnt signals is through the Notch pathway, we knocked out the Dll1 gene at time point 240 minutes. It was found that the expression level of the Nkd1gene increased rather than decreased (see Figure 7 (F)). This is different from the simulation result of Wnt knockout. So we conclude that the regulation of Nkd1 by Wnt signals is not through the Notch pathway. These observations are supported by the research of Aki Ishikawa et al. . Next, we knocked out the Hes7 gene at time point 120 minutes. It was found that the Nkd1gene was upregulated and ceased to oscillate (see Figure 7 (G)). After the Nkd1 gene was knocked out at time point 120 minutes, the expression levels of Wnt target genes (i.e. amplitudes) increased with a time delay about 60 minutes, whereas the expression levels of the Notch target genes decreased immediately and then began to increased about one period later (see Figure 7 (H)). Therefore, Hes7 is essential for the oscillating expression of Nkd1, and Nkd1 inhibits the Wnt pathway, though there is a time delay. These observations are supported by the research of Aki Ishikawa et al. , Yan et al.  and KA Wharton et al. . By the mathematical model, we can give a reasonable explanation of these observations. Nkd1 competes with NICD for Dsh. When Nkd1 is knocked out, more NICD binds to Dsh, so that the unbound NICD concentration decreases. As a result, the expression levels of Notch target genes decrease. But after a while, the expression of the Dll1 gene begin to rise owing to the knocking-out of Nkd1, so that the expression levels of Notch target genes rise with it.
In summary, the simulation results demonstrate that the combined model of the Notch and Wnt pathways is capable of simulating the oscillation of the whole network, and can also reproduce wet experimental results. The model has the potential for further use in research on the molecular mechanisms of somitogenesis.
Parameters sensitivity analysis of the models
The parameter sensitivity analysis is described in the methods section. Here we define a parameter as having significant sensitivity if its sensitivity value is equal to or greater than 1. In simulation experiments, different perturbation levels, low (perturbation 1%), medium (perturbation 10%) and high (perturbation 50%), were used for the parameter sensitivity analysis. All the analysis results are presented in Figures S1 to S12 in Additional file 2. The parameters with significant sensitivity in oscillating periods or amplitudes were picked out and presented in Table S7 to Table S9 in Additional file 2. From the tables we can summarize the characteristics of the isolated pathway models’ sensitivities as follows. (1) The isolated Wnt pathway model is robust to nearly all the parameters in the period and amplitude at all three levels. (2) The periods of target genes in the single Notch pathway model are robust to most of the parameters, and are significantly sensitive to no more than five parameters relevant to the Lfng and Hes7 genes. (3) The target genes in the single Notch pathway model show a diverse spectrum of amplitude sensitivities at the same or different perturbation levels. For the combined model, the periods of the Notch and Wnt target genes are sensitive to the parameters relevant to the Notch pathway, especially the Lfng and Hes7 genes, while they are robust to the parameters relevant to the Wnt pathway. This suggests that the period of somitogenesis is more prone to abnormalities when the Notch pathway is disturbed. This conclusion is in accord with the findings of Leah Herrgen et al. . The amplitudes of the Wnt target genes in the combined model tend to be more sensitive than in the isolated Wnt pathway model.
In this study, we established three mathematical models concerning the Notch and Wnt pathways and their crosstalk during somitogenesis. These models can be used to simulate the oscillating expressions of the target genes of the Notch and/or Wnt pathways and to explore the molecular mechanism of network oscillation. Because the experimental data are limited, this research mainly focuses on the oscillation characteristics of target genes and trends of concentration change of molecules in the pathways, rather than precise experimental results. So the proposed models are qualitative rather than quantitative. For more precise modeling, more experimental data on somitogenesis are still desirable. In addition, it is known that the Notch and Wnt pathways in a single cell are regulated by ligands from adjacent cells in the PSM. In this study, we assumed for simplicity that the pathways in a single cell are regulated by the ligands from that cell. The somitogenesis-related pathways of the cells in the PSM are oscillation-synchronous in somitogenesis, so this assumption has a rational basis. We just modeled the dynamics of the Notch and Wnt pathways in a single cell here. We will consider the communication between cells and establish mathematical models among a group of cells in the PSM in future research.
Although some models for the Notch and Wnt pathways in somitogenesis have been proposed in previous research, these models are too simple for use in researching the molecular mechanism of somitogenesis. Based on recent literature, we proposed two more complicated models for the isolated Notch and Wnt pathways by improving structure design and the representation of chemical reactions. On this basis, we proposed a new structure of the combined model by considering multiple levels of crosstalk between the two pathways. These models are able to simulate the periodic expressions of the Notch or/and Wnt target genes correctly, and can also reproduce biological experimental results from many different publications. They have the potential for use in exploring the molecular mechanisms of the Notch and Wnt signaling pathways and their crosstalk during somitogenesis. The simulation experiments demonstrate that the extracellular signals of the Notch and Wnt signaling pathways are essential for the oscillating expressions not only of their own but also each other’s target genes. Moreover, the negative feedback loops in the Notch and Wnt pathways play important but distinct roles in maintaining the system oscillation. The parameter sensitivity analysis of the models suggests that compared to the Wnt pathway, the Notch pathway is more sensitive to perturbation, so it is prone to abnormalities in somitogenesis.
The approaches to modeling
All models in this study were established using ODEs. The structure and the equations of the models were created using CellDesigner , and were exported as SBML format files. The SBML files were imported into MATLAB and the format converted, so that they can be read by the SBToolBox . In the SBToolBox framework, we accomplished the tasks of parameter learning and parameter sensitivity analysis. Ultimately, the final models were imported into COPASI . All of the model simulations were performed under the COPASI environment. The equations of the models are presented in Additional file 1. The models established using CellDesigner, COPASI and SBToolBox are presented in Additional file 3.
Parameter learning algorithm
We first assigned each species a random initial value, and then trained the model parameters. In general, the system reached semi-steady oscillation after a few abnormal periods. When the semi-steady oscillation was reached, the values of species at that time point were chosen as their initial values. Some of the parameters in the models were derived from , but most of them were estimated by using the parameter learning algorithm. The training set for parameter learning was a set of microarray data derived from  available in Additional file 3. The initial values and parameters of the three models are presented in Tables S1 to S6 in Additional file 1.
Where N is the number of target genes; is the simulating expression profile of the i th target gene; V i is the real expression profile of the i th target gene; M is the time point number of the gene expression profile; x j is the expression value of gene expression profile X at time point j; is the average expression value of gene expression profile X at all time points; y j is the expression value of gene expression profile Y at time point j; is the average expression value of gene expression profile Y at all time points; corr (X, Y) is the correlation coefficient between gene expression profiles X and Y. A PSO-embedded SA algorithm was proposed to explore the optimal parameter set. The algorithm pseudocode is presented in Additional file 4. Its MATLAB code is available in Additional file 5.
Parameter sensitivity analysis
Where S τ is the period sensitivity; τ norm is the period of a gene’s oscillation at normal parameter value; τ pert is the period of a gene’s oscillation after perturbation; p norm is the normal parameter value; p pert is the perturbed parameter value; S A is the amplitude sensitivity; A norm is the amplitude of a gene’s oscillation at normal parameter value; A pert is the amplitude of a gene’s oscillation after perturbation.
This work was supported by the National Natural Science Foundation of China (No. 61172183), the Natural Science Foundation of Jilin Province (No. 20101503), the Research Foundation of Jilin Provincial Science & Technology Development (No. 201201065), the Fundamental Research Funds for the Central Universities (Grant No. 11QNJJ021), and a Grant from Jilin Province Science & Technology Committee (No. 20110711).
- Gilbert SF: Developmental Biology, 6th edition. 2000, Sunderland (MA): Sinauer AssociatesGoogle Scholar
- Dequéant M-L, Pourquié O: Segmental patterning of the vertebrate embryonic axis. Nat Rev Genet. 2008, 9: 370-382.View ArticlePubMedGoogle Scholar
- Aulehla A, Pourquie O: Oscillating signaling pathways during embryonic development. Curr Opin Cell Biol. 2008, 20: 632-637. 10.1016/j.ceb.2008.09.002.View ArticlePubMedGoogle Scholar
- Gibb S, Zagorska A, Melton K, Tenin G, Vacca I, Trainor P, Maroto M, Dale JK: Interfering with Wnt signalling alters the periodicity of the segmentation clock. Dev Biol. 2009, 330 (1): 21-31. 10.1016/j.ydbio.2009.02.035.PubMed CentralView ArticlePubMedGoogle Scholar
- Gibb S, Maroto M, Dale JK: The segmentation clock mechanism moves up a notch. Trends Cell Biol. 2010, 20 (10–2): 593-600.PubMed CentralView ArticlePubMedGoogle Scholar
- Aulehla A, Wehrle C, Brand-Saberi B, Kemler R, Gossler A, Kanzler B, Herrmann BG: Wnt3a Plays a Major Role in the Segmentation Clock Controlling Somitogenesis. Dev Cell. 2003, 4: 395-406. 10.1016/S1534-5807(03)00055-8.View ArticlePubMedGoogle Scholar
- Lewis J: Autoinhibition with Transcriptional Delay: A Simple Mechanism for the Zebrafish Somitogenesis Oscillator. Curr Biol. 2003, 13: 1398-1408. 10.1016/S0960-9822(03)00534-7.View ArticlePubMedGoogle Scholar
- Agrawal S, Archer C, Schaffer DV: Computational Models of the Notch Network Elucidate Mechanisms of Context-dependent Signaling. PLoS Comput Biol. 2009, 5 (5): e1000390-10.1371/journal.pcbi.1000390.PubMed CentralView ArticlePubMedGoogle Scholar
- Terry AJ, Sturrock M, Dale JK, Maroto M, Chaplain MAJ: A Spatio-Temporal Model of Notch Signalling in the Zebrafish Segmentation Clock: Conditions for Synchronised Oscillatory Dynamics. PLoS One. 2010, 6 (2): e16980-View ArticleGoogle Scholar
- Jensen PB, Pedersen L, Krishna S, Jensen MH: A Wnt Oscillator Model for Somitogenesis. Biophys J. 2010, 98: 943-950. 10.1016/j.bpj.2009.11.039.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodríguez-González JG, Santillán M, Fowler AC, Mackey MC: The segmentation clock in mice: Interaction between the Wnt and Notch signalling pathways. J Theor Biol. 2007, 248: 37-47. 10.1016/j.jtbi.2007.05.003.View ArticlePubMedGoogle Scholar
- Goldbeter A, Pourquie O: Modeling the segmentation clock as a network of coupled oscillations in the Notch, Wnt and FGF signaling pathways. J Theor Biol. 2008, 252: 574-585. 10.1016/j.jtbi.2008.01.006.View ArticlePubMedGoogle Scholar
- Santilla´n M, Mackey MC: A Proposed Mechanism for the Interaction of the Segmentation Clock and the Determination Front in Somitogenesis. PLoS One. 2008, 3 (2): e1561-10.1371/journal.pone.0001561.View ArticleGoogle Scholar
- Kazama A, Karashima A, Katayama N, Nakao M: Modeling of Segmentation Clock Mechanism in Presomitic Mesoderm. 31st Annual International Conference of the IEEE EMBS. 2009, Hilton Minneapolis, Minnesota, 3267-3270.Google Scholar
- Monk NAM: Oscillatory Expression of Hes1, p53, and NF-κB Driven by Transcriptional Time Delays. Curr Biol. 2003, 13 (16): 1409-1413. 10.1016/S0960-9822(03)00494-9.View ArticlePubMedGoogle Scholar
- Jouve C, Palmeirim I, Henrique D, Beckers J, Gossler A, Ish-Horowicz D, Pourquié O: Notch signalling is required for cyclic expression of the hairy-like gene HES1 in the presomitic mesoderm. Development. 2000, 127: 1421-1429.PubMedGoogle Scholar
- Dequ M-L, Pourqui O: Segmental patterning of the vertebrate embryonic axis. Nat Rev Genet. 2008, 9: 370-382.View ArticleGoogle Scholar
- Ferjentsik Z, Hayashi S, Dale JK, Bessho Y, An Herreman BDS, Monte G, Pompa JL, Maroto M: Notch Is a Critical Component of the Mouse Somitogenesis Oscillator and Is Essential for the Formation of the Somites. PLoS Genet. 2009, 5 (9): e1000662-10.1371/journal.pgen.1000662.PubMed CentralView ArticlePubMedGoogle Scholar
- Clevers H: Wnt/β -Catenin Signaling in Development and Disease. Cell Mol Life Sci. 2006, 127: 469-480.Google Scholar
- Huelsken J, Behrens J: The Wnt signalling pathway. J Cell Sci. 2002, 115: 3977-3978. 10.1242/jcs.00089.View ArticlePubMedGoogle Scholar
- Boras-Granic K, Chang H, Grosschedl R, Hamel PA: Lef1 is require d for the transition of Wnt signaling from mesenchymal to epithelial cells in the mouse embryonic mammary gland. Dev Biol. 2006, 295: 219-231. 10.1016/j.ydbio.2006.03.030.View ArticlePubMedGoogle Scholar
- Galceran J, Sustmann C, Hsu S-C, Folberth S, Grosschedl R: LEF1-mediated regulation of Delta-like1links Wnt and Notch signaling in somitogenesis. RESEARCH COMMUNICATION. 2004, 18: 2718-2723.Google Scholar
- Hofmann M, Schuster-Gossler K, Watabe-Rudolph M, Aulehla A, Herrmann BG, Gossler A: WNT signaling, in synergy with T/TBX6, controls Notch signaling by regulating Dll1 expression in the presomitic mesoderm of mouse embryos. RESEARCH COMMUNICATION. 2004, 18: 2712-2717.Google Scholar
- Axelrod JD, Matsuno K, Artavanis-Tsankonas S, Perrimon N: Interaction between Wingless and Notch signaling pathways mediated by dishevelled. Science. 1996, 271: 1826-1832. 10.1126/science.271.5257.1826.View ArticlePubMedGoogle Scholar
- Ishikawa A, Kitajima S, Takahashi Y, Kokubo H, Kanno J, Inoue T, Saga Y: Mouse Nkd1, a Wnt antagonist, exhibits oscillatory gene expression in the PSM under the control of Notch signaling. Mech Dev. 2004, 121: 1443-1453. 10.1016/j.mod.2004.08.003.View ArticlePubMedGoogle Scholar
- Yan D, Wallingford JB, Sun T-Q, Nelson AM, Sakanaka C, Reinhard C, Harland RM, Fantl WJ, Williams LT: Cell autonomous regulation of multiple Dishevelled-dependent pathways by mammalian Nkd. PNAS. 2001, 98 (7): 3802-3807. 10.1073/pnas.071041898.PubMed CentralView ArticlePubMedGoogle Scholar
- Wharton KA, Zimmermann G, Rousset R, Scott MP: Vertebrate Proteins Related to Drosophila Naked Cuticle Bind Dishevelled and Antagonize Wnt Signaling. Dev Biol. 2001, 234: 93-106. 10.1006/dbio.2001.0238.View ArticlePubMedGoogle Scholar
- Herrgen L, Ares S, Morelli LG, Schröter C, Jülicher F, Oates AC: Intercellular coupling regulates the period of the segmentation clock. Curr Biol. 2010, 20: 1244-1253. 10.1016/j.cub.2010.06.034.View ArticlePubMedGoogle Scholar
- Funahashi A, Morohashi M, Kitano H, Tanimura N: Cell Designer: a process diagram editor for gene-regulatory and biochemical networks. Biosilico. 2003, 1 (5): 159-162. 10.1016/S1478-5382(03)02370-9.View ArticleGoogle Scholar
- Schmidt H, Jirstrand M: Systems Biology Toolbox for MATLAB: A computational platform for research in Systems Biology. Bioinformatics Applications Note. 2005, 22 (6): 514-515.Google Scholar
- Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U: COPASI — a COmplex PAthway SImulator. Bioinformatics. 2006, 22 (24): 3067-3074. 10.1093/bioinformatics/btl485.View ArticlePubMedGoogle Scholar
- Dequéant M-L, Glynn E, Gaudenz K, Wahl M, Chen J, Mushegian A, Pourquié O: A Complex Oscillating Network of Signaling Genes Underlies the Mouse Segmentation Clock. Science. 2006, 314: 1595-1598. 10.1126/science.1133141.View ArticlePubMedGoogle Scholar
- Varma A, Morbidelli M, Wu H: Parametric Sensitivity in Chemical Systems. 1999, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Stelling J, Gilles ED, III FJD: Robustness properties of circadian clock architectures. PNAS. 2004, 101: 13210-13215. 10.1073/pnas.0401463101.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.