Cell Fate Specification Based on Tristability in the Inner Cell Mass of Mouse Blastocysts

During development, interactions between transcription factors control the specification of different cell fates. The regulatory networks of genetic interactions often exhibit multiple stable steady states; such multistability provides a common dynamical basis for differentiation. During early murine embryogenesis, cells from the inner cell mass (ICM) can be specified in epiblast (Epi) or primitive endoderm (PrE). Besides the intracellular gene regulatory network, specification is also controlled by intercellular interactions involving Erk signaling through extracellular Fgf4. We previously proposed a model that describes the gene regulatory network and its interaction with Erk signaling in ICM cells. The model displays tristability in a range of Fgf4 concentrations and accounts for the self-organized specification process observed in vivo. Here, we further investigate the origin of tristability in the model and analyze in more detail the specification process by resorting to a simplified two-cell model. We also carry out simulations of a population of 25 cells under various experimental conditions to compare their outcome with that of mutant embryos or of embryos submitted to exogenous treatments that interfere with Fgf signaling. The results are analyzed by means of bifurcation diagrams. Finally, the model predicts that heterogeneities in extracellular Fgf4 concentration play a primary role in the spatial arrangement of the Epi/PrE cells in a salt-and-pepper pattern. If, instead of heterogeneities in extra- cellular Fgf4 concentration, internal fluctuations in the levels of expression of the transcription factors are considered as a source of randomness, simulations predict the occurrence of unrealistic switches between the Epi and the PrE cell fates, as well as the evolution of some cells toward one of these states without passing through the previous ICM state, in contrast to what is observed in vivo.