Autocrine FGF feedback can establish distinct states of Nanog expression in pluripotent stem cells: a computational analysis
© Lakatos et al.; licensee BioMed Central Ltd. 2014
Received: 12 March 2014
Accepted: 8 August 2014
Published: 27 September 2014
The maintenance of stem cell pluripotency is controlled by a core cluster of transcription factors, NANOG, OCT4 and SOX2 – genes that jointly regulate each other's expression. The expression of some of these genes, especially of Nanog, is heterogeneous in a population of undifferentiated stem cells in culture. Transient changes in expression levels, as well as heterogeneity of the population is not restricted to this core regulator, but involve a large number of other genes that include growth factors, transcription factors or signal transduction proteins.
As the molecular mechanisms behind NANOG expression heterogeneity is not yet understood, we explore by computational modeling the core transcriptional regulatory circuit and its input from autocrine FGF signals that act through the MAP kinase cascade. We argue that instead of negative feedbacks within the core NANOG-OCT4-SOX2 transcriptional regulatory circuit, autocrine signaling loops such as the Esrrb - FGF - ERK feedback considered here are likely to generate distinct sub-states within the “ON” state of the core Nanog switch. Thus, the experimentally observed fluctuations in Nanog transcription levels are best explained as noise-induced transitions between negative feedback-generated sub-states. We also demonstrate that ERK phosphorilation is altered and being anti-correlated with fluctuating Nanog expression - in accord with model simulations. Our modeling approach assigns an empirically testable function to the transcriptional regulators Klf4 and Esrrb, and predict differential regulation of FGF family members.
We argue that slow fluctuations in Nanog expression likely reflect individual cell-specific changes in parameters of an autocrine feedback loop, such as changes in ligand capture efficiency, receptor numbers or the presence of crosstalks within the MAPK signal transduction pathway. We proposed a model that operates with binding affinities of multiple transcriptional regulators of pluripotency, and the activity of an autocrine signaling pathway. The resulting model produces varied expression levels of several components of pluripotency regulation, largely consistent with empirical observations reported previously and in this present work.
KeywordsStem cell pluripotency Computational model Autocrine signaling Transcription factors Gene expression heterogeneity Fibroblast growth factors Mitogen-activated protein kinase Regulatory network
Embryonic stem (ES) cells are pluripotent cell populations that can be induced to differentiate into a variety of cell types. Mouse ES cells are derived from the inner cell mass of the blastocyst, and their capacity to either self-renew or differentiate into cells of the three germ layers; thus this cell type is a useful model to dissect the molecular regulation behind pluripotency and differentiation -. The maintenance of stem cell pluripotency is controlled by a core cluster of transcription factors, including NANOG, OCT4 and SOX2. The molecular mechanisms by which these factors act is complex and not completely characterized; however part of their critical activity includes the joint regulation of each other's expression ,. The pluripotent status of stem cells is maintained through high expression levels of these genes, and the downregulation of these factors accompanies and is required for cell differentiation. NANOG expression levels are crucial as its forced expression is sufficient to sustain pluripotency even in the absence of extracellular signaling factors such as Leukemia Inhibiting Factor (LIF) , which are otherwise required for stem cell maintenance.
In the last few years it also become clear that - despite their fundamental importance - the expression of many components of the stem cell self-renewal circuitry, including NANOG, is heterogeneous in a population of undifferentiated stem cells ,,. The heterogeneity is dynamically maintained, with individual cells exhibiting transient changes in expression levels. Genes with dynamic expression in mouse ES cells include growth factors, transcription factors, and signal transduction proteins ,,. Recently Galvin-Burgess et al. proposed that undifferentiated ES cells are in various distinct states, depending on the activity of various signal transduction pathways . The presence of states with distinct Nanog expression levels was also suggested based on statistical modeling of changes in flow-sorted populations . Conceptually, the pluripotent and differentiating states of these cells are thus not described well by a simple “ON/OFF” switch. Instead, a cell being in one of the various pluripotent states may be primed or biased in a way that influences its response to differentiation-inducing signals .
As we discuss in the Additional file 1, the equilibrium constants K as well as the cooperativity measures C are related to the binding energies between the transcription factors, the promoter and the RNA polymerase.
The magnitude of model parameters (Additional file 1: Tables S1 and S2) were set by the following considerations. The transcription and translation rates were chosen in such a way that the steady state transcription factor (protein) concentrations are in the nanomolar range (in the order of 100 copy of the TF is present in the cell) when the promoter is fully active ,. To get a functional transcriptional regulatory system, the nanomolar concentration range must be also characteristic for promoter binding affinities, which by Additional file 1: Eq. (S1) translates (at T=300 K) into binding energies around 12 kCal/mol. The transcription factors were assumed to work through stabilizing RNAP II binding - with protein-protein binding energies around 4 kCal/mol . This binding energy is increased for cooperative, multimolecular complexes. We assume that the probability of RNAP II binding in the absence of all the transcription factors considered is very low, [P] K g,P≈10-3.
The core network
First we consider various scenarios for the core NANOG circuit (Figure 2) and compare their behavior. Steady state system behavior was characterized by numerically obtaining intersections of nullcline planes (see Additional file 1: Figure S1).
Our starting point is the model A, which is symmetric in the roles of SOX2 and OCT4 (Figure 2a). This model exhibits bistability: there are two stable fixed points corresponding to the “ON” and “OFF” states of the system, separated by an unstable fixed point. Linear stability analysis reveals that the stable fixed points are stable nodes, thus, no oscillations are expected in their vicinity.
In model C (Figure 2c), OCT4 is independent of NANOG activation. In such a scenario the same parameters that were used for model A yield only the “OFF” fixed point as OCT4 never turns on. Considering the steady presence of OCT4 in embryonic stem cells, we argue that for this model it is reasonable to choose an increased probability for RNAP II binding to the Oct4 locus even in the absence of SOX2. This choice yields a bistable system similar to that of model A.
Model D (Figure 2d) is derived from model C by adding a NANOG autorepression feedback. As Additional file 1: Figure S1d demonstrates, this change does not alter substantially the systems dynamics as the Nanog promoter activity can be obtained by scaling Additional file 1: Eq. (S17) the activity of model C.
Finally, we combine model D with model B to see if NANOG autorepression can further promote the transformation of the “ON” fixed point into a spiral. Adding NANOG autorepression to model B reduces the equilibrium NANOG levels (as determined by the amount of OCT4 present) according to Additional file 1: Eq. (S17), which further stabilizes the fixed point. We also derived model E (Figure 2e) by adding OCT4 as a Nanog repressor to model D. In this scenario we found the fixed point still strongly attractive: oscillations decay fast and change both NANOG and OCT4 levels to a similar extent.
This analysis of the core NANOG circuit variants suggest that they are likely to behave as bistable systems. Hence, the experimentally seen heterogeneity, given the stability of OCT4 expression levels (, Figure 3), is an unlikely consequence of the core NANOG-OCT4-SOX2 dynamics. In the following we will focus on model D, as the simplest variant of the investigated networks, and functionally equivalent with the one proposed by the most recent experimental data . This choice, however, is somewhat arbitrary as all model networks function as a bistable switch.
Sensitivity of the core network to model parameters
Changes in model parameter values can gradually shift the nullclines and fixed points in the phase space. As changes in nullcline positions and shapes can create or remove intersections, the presence of both the “OFF” and “ON” states are parameter dependent.
As we demonstrated above, the most plausible assumptions do not suggest the presence of substantial oscillations within the core NANOG-OCT4-SOX2 system. To explain the observed broad distribution of NANOG expression within a population of mouse ES cells, we hypothesize that model parameters such as binding energies depend on the larger biochemical context. For example, ERK activity is a known potent negative regulator of NANOG transcription ,. In turn, FGF signaling is capable to create an autocrine feedback using cell surface receptors that feed into the MAPK pathway, and FGF activity is indeed a well established modulator of ES cell heterogeneity . Nodal signaling gives an alternative possibility for an autocrine extracellular regulation of Nanog, mediated through the Smad family of transcription factors .
To demonstrate that autocrine signaling loops can induce fluctuating Nanog expression levels, we consider a feedback through FGFs and ERK. As a particular example, we investigate a scenario in which the stability of the RNAP II complex containing the OCT4/SOX2 dimer depends on ERK activation. To explore the behavior of the FGF pathway during NANOG fluctuations, mouse ES cells were sorted based on their NANOG expression levels into four groups (low, medium, high, very high). For each group we determined the transcriptional activity of key genes, such as OCT4, ESRRB, and FGF family members expressed by mouse ES cells: FGF4, FGF5 and FGF8. As Figure 3 demonstrates, the range of Oct4 variability is less than half of that of Nanog. In contrast, the range of variability in the expression levels of FGFs and ESRRB is even greater than that of NANOG. High FGF4 expression is associated with high Nanog expression, whereas high FGF5 and FGF8 expression is characteristic for cells with low levels of Nanog. ERK activity (phosphorilation) assays revealed an inverse relationship between ERK activity and Nanog expression levels (Figure 3). This correlation between Nanog expression and ERK activity is seen in both sorted Nanog subpopulations (Figure 3b) and in a doxycycline-inducible Nanog ES cell line (Figure 3c).
KLF4 and ESRRB
Autocrine FGF signaling
where 0<p Fgf- i <1 is the probability of transcription at the Fgf4 and the Fgf5 locus, and δ FGF-i is the decay rate combined with the diffusive flux transporting the ligand off the cell surface. The α Fgf- i coefficient reflects both production and the conversion between autocrine ligand concentration and production. In the EGF system the relationship between steady state autocrine ligand concentration and its production rate was 0.05 pM/(molecules/cell/h) ≈ 0.05 pM/(30 pM/h) ≈ 0.1% h (see Figure six of ). Based on this result, we expect that less than 1% of the FGF molecules produced in an hour (our approximate time unit) will act as autocrine ligands at the cell surface. Therefore, providing autocrine ligand concentration in the nanomolar range requires higher production rates than the rates we assumed for transcription factors (see Additional file 1: Table S1).
where R t o t is the total amount of FGF receptors at the cell surface, and K is the binding constant between the receptors and their FGF ligands.
where the coupling factor a is chosen in such a way that for a typical simulation the magnitude of (14) is smaller than one.
NANOG expression determined by autocrine feedback
The mechanism behind Nanog fluctuations
Protein and mRNA levels can fluctuate due to the inherently stochastic kinetics of gene transcription and translation. Even an unregulated, constitutively expressed gene exhibits a 30% spread of expression levels over a population of identical cells . Clearly, the reported variations in the stem cell maintenance network, especially those genes that are downstream of Nanog, greatly exceed this baseline variability. Thus, we expect that the observed heterogeneity in Nanog expression levels is generated by the dynamics of the regulatory system as it greatly amplifies the molecular stochastic noise. To describe a potential mechanism generating the experimentally observed dynamics, we proposed a model that operates with binding affinities of multiple transcriptional regulators of pluripotency, the topology of the transcriptional regulatory networks, and activity of an autocrine signaling pathway. The resulting model produces varied expression levels of several components of pluripotency regulation, largely consistent with our (Figure 3) and previously reported empirical observations .
We argue that slow fluctuations in Nanog expression likely reflect individual cell-specific changes in parameters of an autocrine feedback loop, such as changes in ligand capture efficiency, receptor numbers or the presence of crosstalks within the MAPK signal transduction pathway. While high-frequency variability may be filtered out by the slow dynamics of transcription factor synthesis and accumulation, low frequency changes, such as a slow alteration in the cell’s microenvironment, are capable to push expression levels across substate boundaries. Given the complexities of a cell’s variable exposure to autocrine/juxtacrine signaling in culture, this model incorporates a plausible basis for a variable activity of an autocrine signaling pathway eliciting heterogeneous expression of intracellular components. In this view the fluctuations of Nanog reflect the response of a regulatory system with multiple feedbacks in a non-stationary environment.
This feedback mechanism, that does not involve changes in OCT4, is consistent with both the stability of Oct4 expression levels , (Figure 3) and the observation that alteration of OCT4 levels induces ES cells to differentiate rapidly; too much or too little OCT4 rapidly directs cells to differentiate ,. Very small increases in OCT4 expression causes differentiation to mesoderm and endoderm, and reduction of OCT4 levels induces loss of pluripotency and dedifferentiation to trophectoderm.
Comparison with experimental data
The proposed computational model needs to be compared with experimental data on at least three levels. (1) The dynamics of the whole system can be evaluated and compared with corresponding experimental data. (2) Our model makes explicit or implicit assumptions on the topology of the transcriptional regulatory network, and identifies functional modules. Finally, (3) most model parameters are expressed either as molecular binding affinities or as parameters effecting the stability of the transcription complex in addition to production and decay rates.
Our model calculations demonstrate that a biologically plausible autocrine feedback can create distinct sub-states within the “ON” state of the core Nanog switch. While stochasticity and feedback regulation has been proposed to explain Nanog fluctuations ,,, the previously proposed mechanisms involved noise induced transitions between the “ON” and “OFF” states of the core Nanog switch. Experimental evidence, however, suggests that Nanog expression in Nanog-low ES cells is still much higher than that in cells committed to differentiation . Furthermore, Nanog-low cells can still be maintained indefinitely without committing to differentiation - their Nanog expression level will, in fact, increase. These empirical observations clearly support a mechanism that operates with transitions between substates that do not involve the “OFF” state of the core Nanog switch. The existence of distinct substates of Nanog expression was also suggested by a recent study that analyzed experimentally observed changes in Nanog expression profiles in terms of a phenomenological model that made no explicit assumptions on the underlying molecular signaling network .
Our approach also demonstrates that the slow modulation of Nanog expression can reflect changes that are external to the core circuit - instead of the stochastic expression of the transcription factors . The topology of the regulatory network considered here is more elaborated than in previous studies ,, and thus able to represent and predict changes in expression of several downstream genes. Furthermore, our model does not include direct autoregulation of Nanog - a frequent assumption in simplified models which lacks empirical support.
FGF and NODAL signaling are clearly active in mouse embryonic stem cells, and function in an autocrine fashion in undifferentiated cells . Previous studies , have shown that these autocrine signaling pathways influence the dynamic heterogeneity of ES cells in culture, likely through specific molecular mechanisms that have yet to be elaborated. Our results show highly variable expression of FGF ligands and intracellular ERK activity in ES cells grown in serum-based media (Figure 3). The model behavior presented in Figure 9, is largely consistent with the anticorrelation found in ERK activity and Nanog expression (Figure 3) as well as existing information on FGF signaling .
Regulatory network topology
The regulatory network shown in Figure 5 contains motifs that are well established as well as hypothetical regulatory connections that we explore in this work. The full network can be broken up into four modules: the core NANOG switch, the downstream ESRRB amplifier, the autocrine FGF module and finally, the feedback through the MAPK cascade.
Our analysis indicates that each of the recently proposed regulatory architecture of the core NANOG-OCT4-SOX2 network (Figure 2) functions as a bistable switch. We argue, that irrespective of the underlying model details, the nonlinearity needed to create two stable fixed points is provided by the autoregulation and the dimerization of the OCT4/SOX2 transcription factors. The switch-like behavior is also maintained when OCT4 acts as a likely low affinity repressor - an assumption motivated by the finding that OCT4 overexpression was required for the manifestation of the repressive behavior .
The roles KLF4 and ESRRB play in stem cell maintenance are in the focus of recent scientific interest ,,. Based on chromatin immunoprecipitation data , both the Klf4 and Esrrb genes have several potential regulatory sites, which include transcription factors from the Nanog core as well as allow autoregulatory feedbacks. We demonstrated that a subset of the known likely regulatory interactions is capable to function as an amplifier, greatly expanding the variability of these factors. While the existence of (functional) autoregulation in these loci has not yet been established by targeted experiments, the proposed amplifier function of the KLF4-ESRRB system is in good agreement with the experimentally obtained large variations in Esrrb (Figure 3 and ) and downstream Fgf expression levels (Figure 3).
The Fgf genes which are known to be actively transcribed in mouse ES cells exhibit regulatory sites for the transcription factors considered in this work. Chromatin immunoprecipitation data reveals that Fgf4 has putative regulatory sites for KLF4 and ESRRB, Fgf5 has putative regulatory sites for NANOG and ESRRB, while Fgf8 exhibits a regulatory site for NANOG . While the nature of the regulatory connection is currently unknown, based on the strong correlation between Nanog and Fgf4 expression and the strong anticorrelation between Nanog and Fgf5 as well as Nanog and Fgf8 expression levels, we assume that Fgf5 and Fgf8 are repressed while Fgf4 expression is enhanced by the transcription factors considered in this work. The repressional regulation of Fgf genes is of key importance - it yields a negative feedback loop, and this will tested by targeted experiments in the future.
Finally, the MAPK module and the autocrine ligand concentration was treated phenomenologically by simplified input-output relationships reflecting previously established results. In particular, we postulated a linear relationship between autocrine ligand concentration and its rate of production ,. The MAPK pathway can exhibit transient and sustained activation ,, the latter lasting only for less than an hour . As changes in transcriptional regulation are operating on a slower time scale, here we considered a change in sustained ERK activation, which was also carefully studied and is expected to operate as an amplifier (ultrasensitive switch) linking input and output with a Hill coefficient much greater than one ,. This observation motivates our use of a coupling factor a in Eq. (14) that is able to substantially reduce Nanog expression, in agreement with specific observations in mouse ES cells ,.
A molecular regulatory network model is bound to have several parameters, for most of which very little empirical data is available. As one way to contain this problem, here we used a modular approach: first we identified the behavior of smaller units, like the core Nanog switch, or the Esrrb amplifier. Functional requirements can considerably constrain the possible parameter values within such a module. As a second effort to reduce the arbitrariness of model parameters we derived expression rates from binding energies ,. Finally, unless empirical evidence suggested otherwise, we kept the simplest (uniformly assigned) values of the parameters and did not try to match empirical data by fine tuning the parameters. Nevertheless, we do not consider our parameters to be predictive - many different combinations can give similar overall behavior. Yet, they are important to demonstrate that the proposed regulatory system can work with plausible parameter values.
The major focus of this model involves a regulatory network of Nanog, Oct4, and Sox2; these transcriptional regulators constitute a semi-independent regulatory module (‘core’ module ) in maintaining pluripotency. Clearly many other transcriptional regulators are involved in this complex module however. Our model represents an extensible platform for adding further components of the core module. The integration of other transcriptional regulators such as Tcf3 , Tbx3 , and signaling pathways such as Nodal and BMP signaling , will be an important future application of this model to study the inherent instability of the pluripotent phenotype in serum-based media. The model presented here does not take into account potential role of allelic expression of the Nanog locus  in directing heterogeneity in Nanog expression. Experimental data for this phenomenon are contradictory, however, . Further experimental data will be required to confirm this phenomenon, and if so, how this unusual mode of transcriptional regulation may be integrated into the model presented here.
The pluripotency of embryonic stem cells is sustained by a core cluster of co-regulated transcription factors, NANOG, OCT4 and SOX2. Surprisingly, Nanog as well as several other downstream transcription factors exhibit widely fluctuating expression in a population of undifferentiated stem cells. To explain the observed heterogeneity of expression levels, we propose a computational model that couples the transcriptional regulation of Nanog to autocrine extracellular signals. We argue that the likely source of fluctuations is not the core regulatory cluster, but stochasticity within the autocrine feedback loops. The model predicts fluctuating expression levels for several factors involved in pluripotency maintenance, largely consistent with empirical observations presented here or reported previously. Our model indicates the presence of distinct substates of pluripotency, each exhibiting various expression levels of Nanog and downstream transcription or signaling factors. Thus, we predict that subpopulations within undifferentiated embryonic stem cells can have non-uniform responses to extracellular stimuli. Finally, we assign an empirically testable function to the transcriptional regulators Klf4 and Esrrb, and predict differential regulation of FGF family members. In vivo, we expect autocrine feedbacks to be highly sensitive to alterations of the stem cell microenvironment, hence our model can guide future studies linking stem cell maintenance to physical and biochemical properties of the stem cell niche.
Mouse ES cell culture
Experiments used the BAC-Nanog-GFP (BNG) ES cells  in which Nanog-GFP bacterial artificial chromosome (BAC)  was introduced into Ainv15 ES cells . BNG ES cells were maintained as described previously , on gelatin-coated or fibroblast cocultured plates. Cells were grown in serum-based ES cell media: Dulbecco’s modified Eagle’s medium, 15% fetal bovine serum (FBS), penicillin-streptomycin, L-glutamine, nonessential amino acids, β-mercaptoethanol, and 103 units/ml leukemia inhibitory factor (LIF).
Fluorescent cell sorting and analysis
ES cells were trypsinized to a single cell suspension and analyzed by BD FACSAria (cell sorting) and BD LSRII flow cytometer (cell analysis) for GFP expression. A convention was established for sorting and analyzing subpopulations of BNG ES cells based on the profile of GFP expression of unsorted cells, as described in . The GFP-medium and GFP-high populations were determined by gating 30%–35% of the cells from the peak of the GFP distribution. Cells expressing higher levels of GFP than the GFP-high subpopulation were classified as GFP-very high; whereas cells with expression less than GFP-medium cells but above background were identified as GFP-low cells. Analysis of control E14 ES cells were used as a negative control for flow analysis and sorting. After cell sorting, subpopulations were analyzed to confirm purity of the population.
RNA was isolated (Qiagen, Hilden, Germany), and cDNA was synthesized (Invitrogen) following manufacturer’s instructions. TaqMan primer sets with the 7500 Real Time polymerase chain reaction (PCR) system (Applied BioSystems, Foster City, CA) were used for quantitative real-time PCR analysis. Upon request, specific ABI TaqMan Primer/Probe assay identification numbers are available.
Cells were pelleted and then lysed with radioimmune precipitation assay buffer supplemented with Halt Protease and Phosphatase Inhibitor Cocktails (Pierce, Rockford, IL). Protein samples were separated on BioRad Tris-HCl gels, and blots were probed with primary antibodies for GAPDH (Santa Cruz Biotechnology, Santa Cruz, CA) and pERK (Cell Signaling) and incubated overnight at 4 deg C with appropriate secondary antibodies. SuperSignal West Pico Chemiluminescence (Pierce) was used to detect the Western blots.
Dox-inducible Nanog ES cell line
The Nanog coding sequence was inserted into a dox-responsive locus via CRE-mediated insertion using previously described methods ,, into the AInv-15 ES cell line, to generate BNG-dox-Nanog cell line. These cells were cultured under feeder-free conditions in standard serum based ES cell media supplemented with doxycycline (1 ug/ml) for the indicated times. Cells were then harvested in RIPA buffer for western analysis.
DL programmed and conceptualized the model, performed the simulations and analyzed the data. EDT performed the experiments. KEP developed and analyzed the dox-Nanog cell line. JLV analyzed the experiments and conceptualized the model. AC conceptualized the model and analyzed the data. All authors read and approved the final manuscript.
This work was funded through grants from the Hungarian Development Agency (KTIA AIK 12-1-2013-0041 to AC), the Hungarian Science Fund (OTKA K72664 to AC) and the NIH (COBRE P20GM104936 to JLV).
- Evans MJ, Kaufman MH: Establishment in culture of pluripotential cells from mouse embryos. Nature. 1981, 292 (5819): 154-156. 10.1038/292154a0.View ArticlePubMedGoogle Scholar
- Smith AG: Embryo-derived stem cells: of mice and men. Annu Rev Cell Dev Biol. 2001, 17: 435-462. 10.1146/annurev.cellbio.17.1.435.View ArticlePubMedGoogle Scholar
- Ohtsuka S, Dalton S: Molecular and biological properties of pluripotent embryonic stem cells. Gene Ther. 2008, 15 (2): 74-81. 10.1038/sj.gt.3303065.View ArticlePubMedGoogle Scholar
- Boyer LA, Mathur D, Jaenisch R: Molecular control of pluripotency. Curr Opin Genet Dev. 2006, 16 (5): 455-462. 10.1016/j.gde.2006.08.009.View ArticlePubMedGoogle Scholar
- Pan G, Thomson JA: Nanog and transcriptional networks in embryonic stem cell pluripotency. Cell Res. 2007, 17 (1): 42-49. 10.1038/sj.cr.7310125.View ArticlePubMedGoogle Scholar
- Mitsui K, Tokuzawa Y, Itoh H, Segawa K, Murakami M, Takahashi K, Maruyama M, Maeda M, Yamanaka S: The homeoprotein nanog is required for maintenance of pluripotency in mouse epiblast and es cells. Cell. 2003, 113 (5): 631-642. 10.1016/S0092-8674(03)00393-3.View ArticlePubMedGoogle Scholar
- Chambers I, Silva J, Colby D, Nichols J, Nijmeijer B, Robertson M, Vrana J, Jones K, Grotewold L, Smith A: Nanog safeguards pluripotency and mediates germline development. Nature. 2007, 450 (7173): 1230-1234. 10.1038/nature06403.View ArticlePubMedGoogle Scholar
- Chickarmane V, Troein C, Nuber UA, Sauro HM, Peterson C: Transcriptional dynamics of the embryonic stem cell switch. PLoS Comput Biol. 2006, 2 (9): 123-10.1371/journal.pcbi.0020123.View ArticleGoogle Scholar
- Kalmar T, Lim C, Hayward P, Muñoz-Descalzo S, Nichols J, Garcia-Ojalvo J, Arias AM: Regulated fluctuations in nanog expression mediate cell fate decisions in embryonic stem cells. PLoS Biol. 2009, 7 (7): 1000149-10.1371/journal.pbio.1000149.View ArticleGoogle Scholar
- Navarro P, Festuccia N, Colby D, Gagliardi A, Mullin NP, Zhang W, Karwacki-Neisius V, Osorno R, Kelly D, Robertson M, Chambers I: Oct4/sox2-independent nanog autorepression modulates heterogeneous nanog gene expression in mouse es cells. EMBO J. 2012, 31 (24): 4547-4562. 10.1038/emboj.2012.321.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayashi K, de Sousa Lopes SMC, Tang F, Surani MA: Dynamic equilibrium and heterogeneity of mouse pluripotent stem cells with distinct functional and epigenetic states. Cell Stem Cell. 2008, 3 (4): 391-401. 10.1016/j.stem.2008.07.027.View ArticlePubMedGoogle Scholar
- Lanner F, Lee KL, Sohl M, Holmborn K, Yang H, Wilbertz J, Poellinger L, Rossant J, Farnebo F: Heparan sulfation-dependent fibroblast growth factor signaling maintains embryonic stem cells primed for differentiation in a heterogeneous state. Stem Cells. 2010, 28 (2): 191-200.PubMedGoogle Scholar
- Galvin-Burgess KE, Travis ED, Pierson KE, Vivian JL: Tgf-β-superfamily signaling regulates embryonic stem cell heterogeneity: self-renewal as a dynamic and regulated equilibrium. Stem Cells. 2013, 31 (1): 48-58. 10.1002/stem.1252.PubMed CentralView ArticlePubMedGoogle Scholar
- Luo Y, Lim CL, Nichols J, Martinez-Arias A, Wernisch L: Cell signalling regulates dynamics of nanog distribution in embryonic stem cell populations. J R Soc Interface. 2012, 10 (78): 20120525-10.1098/rsif.2012.0525.View ArticleGoogle Scholar
- Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, Jaenisch R, Wagschal A, Feil R, Schreiber SL, Lander ES: A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006, 125 (2): 315-326. 10.1016/j.cell.2006.02.041.View ArticlePubMedGoogle Scholar
- Buchler NE, Gerland U, Hwa T: On schemes of combinatorial transcription logic. Proc Natl Acad Sci U S A. 2003, 100 (9): 5136-5141. 10.1073/pnas.0930314100.PubMed CentralView ArticlePubMedGoogle Scholar
- Chickarmane V, Olariu V, Peterson C: Probing the role of stochasticity in a model of the embryonic stem cell: heterogeneous gene expression and reprogramming efficiency. BMC Syst Biol. 2012, 6: 98-10.1186/1752-0509-6-98.PubMed CentralView ArticlePubMedGoogle Scholar
- Lanner F, Rossant J: The role of fgf/erk signaling in pluripotent cells. Development. 2010, 137 (20): 3351-3360. 10.1242/dev.050146.View ArticlePubMedGoogle Scholar
- Flicek P, Ahmed I, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fairley S, Fitzgerald S, Gil L, García-Girón C, Gordon L, Hourlier T, Hunt S, Juettemann T, Kähäri AK, Keenan S, Komorowska M, Kulesha E, Longden I, Maurel T, McLaren WM, Muffato M, Nag R, Overduin B, Pignatelli M, Pritchard B, Pritchard E, et al: Ensembl 2013. Nucleic Acids Res 2013. 2003, 41 (Database issue): 48-55.Google Scholar
- Kunath T, Saba-El-Leil MK, Almousailleakh M, Wray J, Meloche S, Smith A: Fgf stimulation of the erk1/2 signalling cascade triggers transition of pluripotent embryonic stem cells from self-renewal to lineage commitment. Development. 2007, 134 (16): 2895-2902. 10.1242/dev.02880.View ArticlePubMedGoogle Scholar
- Festuccia N, Osorno R, Halbritter F, Karwacki-Neisius V, Navarro P, Colby D, Wong F, Yates A, Tomlinson SR, Chambers I: Esrrb is a direct nanog target gene that can substitute for nanog function in pluripotent cells. Cell Stem Cell. 2012, 11 (4): 477-490. 10.1016/j.stem.2012.08.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Jiang J, Chan Y-S, Loh Y-H, Cai J, Tong G-Q, Lim C-A, Robson P, Zhong S, Ng H-H: A core klf circuitry regulates self-renewal of embryonic stem cells. Nat Cell Biol. 2008, 10 (3): 353-360. 10.1038/ncb1698.View ArticlePubMedGoogle Scholar
- Chen X, Xu H, Yuan P, Fang F, Huss M, Vega VB, Wong E, Orlov YL, Zhang W, Jiang J, Loh Y-H, Yeo HC, Yeo ZX, Narang V, Govindarajan KR, Leong B, Shahab A, Ruan Y, Bourque G, Sung W-K, Clarke ND, Wei C-L, Ng H-H: Integration of external signaling pathways with the core transcriptional network in embryonic stem cells. Cell. 2008, 133 (6): 1106-1117. 10.1016/j.cell.2008.04.043.View ArticlePubMedGoogle Scholar
- Marson A, Levine SS, Cole MF, Frampton GM, Brambrink T, Johnstone S, Guenther MG, Johnston WK, Wernig M, Newman J, Calabrese JM, Dennis LM, Volkert TL, Gupta S, Love J, Hannett N, Sharp PA, Bartel DP, Jaenisch R, Young RA: Connecting microrna genes to the core transcriptional regulatory circuitry of embryonic stem cells. Cell. 2008, 134 (3): 521-533. 10.1016/j.cell.2008.07.020.PubMed CentralView ArticlePubMedGoogle Scholar
- Galvin KE, Travis ED, Yee D, Magnuson T, Vivian JL: Nodal signaling regulates the bone morphogenic protein pluripotency pathway in mouse embryonic stem cells. J Biol Chem. 2010, 285 (26): 19747-19756. 10.1074/jbc.M109.077347.PubMed CentralView ArticlePubMedGoogle Scholar
- DeWitt AE, Dong JY, Wiley HS, Lauffenburger DA: Quantitative analysis of the egf receptor autocrine system reveals cryptic regulation of cell response by ligand capture. J Cell Sci. 2001, 114 (Pt 12): 2301-2313.PubMedGoogle Scholar
- Huang CY, Ferrell JE: Ultrasensitivity in the mitogen-activated protein kinase cascade. Proc Natl Acad Sci U S A. 1996, 93 (19): 10078-10083. 10.1073/pnas.93.19.10078.PubMed CentralView ArticlePubMedGoogle Scholar
- Asthagiri AR, Nelson CM, Horwitz AF, Lauffenburger DA: Quantitative relationship among integrin-ligand binding, adhesion, and signaling via focal adhesion kinase and extracellular signal-regulated kinase 2. J Biol Chem. 1999, 274 (38): 27119-27127. 10.1074/jbc.274.38.27119.View ArticlePubMedGoogle Scholar
- Asthagiri AR, Lauffenburger DA: A computational study of feedback effects on signal dynamics in a mitogen-activated protein kinase (mapk) pathway model. Biotechnol Prog. 2001, 17 (2): 227-239. 10.1021/bp010009k.View ArticlePubMedGoogle Scholar
- Cohen AA, Kalisky T, Mayo A, Geva-Zatorsky N, Danon T, Issaeva I, Kopito RB, Perzov N, Milo R, Sigal A, Alon U: Protein dynamics in individual human cells: experiment and theory. PLoS One. 2009, 4 (4): 4901-10.1371/journal.pone.0004901.View ArticleGoogle Scholar
- Niwa H, Miyazaki J, Smith AG: Quantitative expression of oct-3/4 defines differentiation, dedifferentiation or self-renewal of es cells. Nat Genet. 2000, 24 (4): 372-376. 10.1038/74199.View ArticlePubMedGoogle Scholar
- Loh Y-H, Wu Q, Chew J-L, Vega VB, Zhang W, Chen X, Bourque G, George J, Leong B, Liu J, Wong K-Y, Sung KW, Lee CWH, Zhao X-D, Chiu K-P, Lipovich L, Kuznetsov VA, Robson P, Stanton LW, Wei C-L, Ruan Y, Lim B, Ng H-H: The oct4 and nanog transcription network regulates pluripotency in mouse embryonic stem cells. Nat Genet. 2006, 38 (4): 431-440. 10.1038/ng1760.View ArticlePubMedGoogle Scholar
- Glauche I, Herberg M, Roeder I: Nanog variability and pluripotency regulation of embryonic stem cells-insights from a mathematical model analysis. PLoS One. 2010, 5 (6): 11238-10.1371/journal.pone.0011238.View ArticleGoogle Scholar
- Pan G, Li J, Zhou Y, Zheng H, Pei D: A negative feedback loop of transcription factors that controls stem cell pluripotency and self-renewal. FASEB J. 2006, 20 (10): 1730-1732. 10.1096/fj.05-5543fje.View ArticlePubMedGoogle Scholar
- Zhang P, Andrianakos R, Yang Y, Liu C, Lu W: Kruppel-like factor 4 (klf4) prevents embryonic stem (es) cell differentiation by regulating nanog gene expression. J Biol Chem. 2010, 285 (12): 9180-9189. 10.1074/jbc.M109.077958.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim MO, Kim S-H, Cho Y-Y, Nadas J, Jeong C-H, Yao K, Kim DJ, Yu D-H, Keum Y-S, Lee K-Y, Huang Z, Bode AM, Dong Z: Erk1 and erk2 regulate embryonic stem cell self-renewal through phosphorylation of klf4. Nat Struct Mol Biol. 2012, 19 (3): 283-290. 10.1038/nsmb.2217.View ArticlePubMedGoogle Scholar
- van den Berg DLC, Zhang W, Yates A, Engelen E, Takacs K, Bezstarosti K, Demmers J, Chambers I, Poot RA: Estrogen-related receptor beta interacts with oct4 to positively regulate nanog gene expression. Mol Cell Biol. 2008, 28 (19): 5986-5995. 10.1128/MCB.00301-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Shvartsman SY, Wiley HS, Deen WM, Lauffenburger DA: Spatial range of autocrine signaling: modeling and computational analysis. Biophys J. 2001, 81 (4): 1854-1867. 10.1016/S0006-3495(01)75837-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Bhalla US, Iyengar R: Emergent properties of networks of biological signaling pathways. Science. 1999, 283 (5400): 381-387. 10.1126/science.283.5400.381.View ArticlePubMedGoogle Scholar
- Sasagawa S, Ozaki Y-I, Fujita K, Kuroda S: Prediction and validation of the distinct dynamics of transient and sustained erk activation. Nat Cell Biol. 2005, 7 (4): 365-373. 10.1038/ncb1233.View ArticlePubMedGoogle Scholar
- Aoki K, Yamada M, Kunida K, Yasuda S, Matsuda M: Processive phosphorylation of erk map kinase in mammalian cells. Proc Natl Acad Sci U S A. 2011, 108 (31): 12675-12680. 10.1073/pnas.1104030108.PubMed CentralView ArticlePubMedGoogle Scholar
- Shea MA, Ackers GK: The or control system of bacteriophage lambda. A physical-chemical model for gene regulation. J Mol Biol. 1985, 181 (2): 211-230. 10.1016/0022-2836(85)90086-5.View ArticlePubMedGoogle Scholar
- Kim J, Woo AJ, Chu J, Snow JW, Fujiwara Y, Kim CG, Cantor AB, Orkin SH: A myc network accounts for similarities between embryonic stem and cancer cell transcription programs. Cell. 2010, 143 (2): 313-324. 10.1016/j.cell.2010.09.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Wray J, Kalkan T, Gomez-Lopez S, Eckardt D, Cook A, Kemler R, Smith A: Inhibition of glycogen synthase kinase-3 alleviates tcf3 repression of the pluripotency network and increases embryonic stem cell resistance to differentiation. Nat Cell Biol. 2011, 13 (7): 838-845. 10.1038/ncb2267.PubMed CentralView ArticlePubMedGoogle Scholar
- Martello G, Sugimoto T, Diamanti E, Joshi A, Hannah R, Ohtsuka S, Göttgens B, Niwa H, Smith A: Esrrb is a pivotal target of the gsk3/tcf3 axis regulating embryonic stem cell self-renewal. Cell Stem Cell. 2012, 11 (4): 491-504. 10.1016/j.stem.2012.06.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu J, Tzanakakis ES: Distinct allelic patterns of nanog expression impart embryonic stem cell population heterogeneity. PLoS Comput Biol. 2013, 9 (7): 1003140-10.1371/journal.pcbi.1003140.View ArticleGoogle Scholar
- Faddah DA, Wang H, Cheng AW, Katz Y, Buganim Y, Jaenisch R: Single-cell analysis reveals that expression of nanog is biallelic and equally variable as that of other pluripotency factors in mouse escs. Cell Stem Cell. 2013, 13 (1): 23-29. 10.1016/j.stem.2013.04.019.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaniel C, Li F, Schafer XL, Moore T, Lemischka IR, Paddison PJ: Delivery of short hairpin rnas-triggers of gene silencing-into mouse embryonic stem cells. Nat Methods. 2006, 3 (5): 397-400. 10.1038/nmeth0506-397.View ArticlePubMedGoogle Scholar
- Kyba M, Perlingeiro RCR, Daley GQ: Hoxb4 confers definitive lymphoid-myeloid engraftment potential on embryonic stem cell and yolk sac hematopoietic progenitors. Cell. 2002, 109 (1): 29-37. 10.1016/S0092-8674(02)00680-3.View 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.