pubmed.ncbi.nlm.nih.gov

Modelling coupled oscillations in the Notch, Wnt, and FGF signaling pathways during somitogenesis: a comprehensive mathematical model - PubMed

doi: 10.1155/2015/387409. Epub 2015 Mar 17.

Affiliations

Modelling coupled oscillations in the Notch, Wnt, and FGF signaling pathways during somitogenesis: a comprehensive mathematical model

Hong-yan Wang et al. Comput Intell Neurosci. 2015.

Abstract

Somite formation in the early stage of vertebrate embryonic development is controlled by a complicated gene network named segmentation clock, which is defined by the periodic expression of genes related to the Notch, Wnt, and the fibroblast growth factor (FGF) pathways. Although in recent years some findings about crosstalk among the Notch, Wnt, and FGF pathways in somitogenesis have been reported, the investigation of their crosstalk mechanisms from a systematic point of view is still lacking. In this study, a more comprehensive mathematical model was proposed to simulate the dynamics of the Notch, Wnt, and FGF pathways in the segmentation clock. Simulations and bifurcation analyses of this model suggested that the concentration gradients of both Wnt, and FGF signals along the presomitic mesoderm (PSM) are corresponding to the whole process from start to stop of the segmentation clock. A number of highly sensitive parameters to the segmentation clock's oscillatory pattern were identified. By further bifurcation analyses for these sensitive parameters, and several complementary mechanisms in respect of the maintenance of the stable oscillation of the segmentation clock were revealed.

PubMed Disclaimer

Figures

Figure 1
Figure 1

The schematic diagram of the model. The diagram is created using CellDesigner [18]. In the diagram, the outer yellow ellipse represents cytomembrane; the inner yellow ellipse represents nuclear membrane; the light green rectangles represent protein; the bottle green parallelograms represent mRNA; the white rectangles which contain light green rectangles represent complex; Φ represents the resultants of degradation reactions or reactants of combination reactions. The arrow lines represent biochemical reactions.

Figure 2
Figure 2

The simulation results of the model. (a) The oscillatory expression of the target genes under the condition of a constant extracellular Wnt and FGF signal. (b) The expression patterns of the target genes after the induction of the Wnt signal to Msgn1 was removed at time point 240 minutes. (c) The expression patterns of the target genes after the induction of the Wnt signal to the Hes7 gene was removed at time point 240 minutes. (d) The expression patterns of the target genes after the FGF protein was removed at time point 240 minutes. (e) The expression patterns of the target genes after the repression of the Hes7 protein to the Dusp4 gene was removed at time point 240 minutes. (f) The expression patterns of the target genes after the repression of the Hes7 protein to the Dll1 gene was removed at time points 240 minutes.

Figure 3
Figure 3

The results of the bifurcation scanning to the Wnt signal. The X-axis represents time. The Y-axis represents the concentration of the Wnt signal. The Z-axis and the colour bar both represent the concentration of the target genes' expression. (a) The expression level of the Hes7 gene in two periods during scanning the Wnt signals. (b) The expression level of the Lfng gene in two periods during scanning the Wnt signals. (c) The expression level of the Axin2 gene in two periods during scanning the Wnt signals. (d) The expression level of the Dusp4 gene in two periods during scanning the Wnt signals. (e) The expression level of the Dll1 gene in two periods during scanning the Wnt signals. (f) The expression level of the Nkd1 gene in two periods during scanning the Wnt signals.

Figure 4
Figure 4

The results of one-parameter bifurcation analysis for the parameters in the feedback loops. The X-axis represents time. The Y-axis represents the values of the analyzed parameter. The Z-axis and the colour bar both represent the concentration of the target genes' expression. (a) The expression level of the Hes7 gene in respect to k mHes7 with the length of time 600 minutes. (b) The expression level of the Hes7 gene in respect to EX Hes7 with the length of time 600 minutes. (c) The expression level of the Hes7 gene in respect to K aDll1 with the length of time 600 minutes. (d) The expression level of the Hes7 gene in respect to k mLfng with the length of time 600 minutes. (e) The expression level of the Hes7 gene in respect to k NICD with the length of time 600 minutes. (f) The expression level of the Hes7 gene in respect to K ILfng with the length of time 600 minutes. (g) The expression level of the Hes7 gene in respect to k mDusp4 with the length of time 600 minutes. (h) The expression level of the Hes7 gene in respect to IM Erka with the length of time 600 minutes. (i) The expression level of the Hes7 gene in respect to EX Erkn with the length of 600 minutes.

Figure 5
Figure 5

The results of one-parameter bifurcation analysis to the parameters in the reactions of the Hes7 gene's expression. The X-axis represents time. The Y-axis represents the values of the analyzed parameter. The Z-axis and the colour bar both represent the concentration of the target genes' expression. (a) The expression level of the Hes7 gene in respect to K NsmHes7 with the length of time 600 minutes. (b) The expression level of the Hes7 gene in respect to K BsmHes7 with the length of time 600 minutes. (c) The expression level of the Hes7 gene in respect to K EsmHes7 with the length of time 600 minutes. (d) The expression level of the Hes7 gene in respect to K HsmHes7 with the length of time 600 minutes.

Figure 6
Figure 6

The results of two-parameter bifurcation analysis of the parameter pairs from the feedback loops centred on Hes7 and Lfng. The X-axis and Y-axis represent two parameters analyzed, respectively. The colour bar represents the normalized stable degree (see Equation (1)) of the Hes7 gene's oscillation at some a parameter value pair and it has no unit. (a) k mHes7-K aDll1. (b) EX Hes7-K aDll1. (c) k mHes7-k mLfng. (d) EX Hes7-k mLfng. (e) k mHes7-k NICD. (f) EX Hes7-k NICD. (g) k mHes7-K ILfng. (h) EX Hes7-K ILfng.

Figure 7
Figure 7

The expression of the Hes7 gene under the different conditions of k mHes7, k mLfng, k mDusp4 and K aDll1. (a) k mHes7 = 0.1 min−1. (b) k mLfng = 0.01 min−1. (c) k mHes7 = 0.1 min−1 and k mLfng = 0.01 min−1. (d) k mHes7 = 0.27 min−1. (e) k mLfng = 0.03 min−1. (f) k mHes7 = 0.1 min−1 and k mLfng = 0.03 min−1. (g) k mHes7 = 0.15 min−1. (h) k mDusp4 = 0.05 min−1. (i) k mHes7 = 0.15 min−1 and k mDusp4 = 0.05 min−1. (j) k mHes7 = 0.27 min−1. (k) k mDusp4 = 0.11 min−1. (l) k mHes7 = 0.27 min−1 and k mDusp4 = 0.11 min−1. (m) k mHes7 = 0.1 min−1. (n) K aDll1 = 0.5 nM. (o) k mHes7 = 0.1 min−1 and K aDll1 = 0.5 nM. (p) k mHes7 = 0.27 min−1. (q) K aDll1 = 1.1 nM. (r) k mHes7 = 0.27 min−1 and K aDll1 = 1.1 nM.

Similar articles

Cited by

References

    1. Gilbert S. F. Developmental Biology. 6th. Sunderland, Mass, USA: Sinauer Associates; 2000.
    1. Dequéant M.-L., Pourquié O. Segmental patterning of the vertebrate embryonic axis. Nature Reviews Genetics. 2008;9(5):370–382. doi: 10.1038/nrg2320. - DOI - PubMed
    1. Aulehla A., Pourquié O. Oscillating signaling pathways during embryonic development. Current Opinion in Cell Biology. 2008;20(6):632–637. doi: 10.1016/j.ceb.2008.09.002. - DOI - PubMed
    1. Dequéant M. L., Pourquié O. Segmental patterning of the vertebrate embryonic axis. Nature Reviews Genetics. 2008;9(5):370–382. doi: 10.1038/nrg2320. - DOI - PubMed
    1. Cooke J. A gene that resuscitates a theory—somitogenesis and a molecular oscillator. Trends in Genetics. 1998;14(3):85–88. doi: 10.1016/s0168-9525(98)01396-1. - DOI - PubMed

Publication types

MeSH terms

Substances

LinkOut - more resources