RESUMEN
Understanding cell fate selection remains a central challenge in developmental biology. We present a class of simple yet biologically motivated mathematical models for cell differentiation that generically generate oscillations and hence suggest alternatives to the standard framework based on Waddington's epigenetic landscape. The models allow us to suggest two generic dynamical scenarios that describe the differentiation process. In the first scenario, gradual variation of a single control parameter is responsible for both entering and exiting the oscillatory regime. In the second scenario, two control parameters vary: one responsible for entering, and the other for exiting the oscillatory regime. We analyse the standard repressilator and four variants of it and show the dynamical behaviours associated with each scenario. We present a thorough analysis of the associated bifurcations and argue that gene regulatory networks with these repressilator-like characteristics are promising candidates to describe cell fate selection through an oscillatory process.
Asunto(s)
Epigénesis Genética , Células Madre , Diferenciación Celular , Redes Reguladoras de Genes , Modelos Genéticos , Modelos TeóricosRESUMEN
Neural crest cells are crucial in development, not least because of their remarkable multipotency. Early findings stimulated two hypotheses for how fate specification and commitment from fully multipotent neural crest cells might occur, progressive fate restriction (PFR) and direct fate restriction, differing in whether partially restricted intermediates were involved. Initially hotly debated, they remain unreconciled, although PFR has become favoured. However, testing of a PFR hypothesis of zebrafish pigment cell development refutes this view. We propose a novel 'cyclical fate restriction' hypothesis, based upon a more dynamic view of transcriptional states, reconciling the experimental evidence underpinning the traditional hypotheses.
Asunto(s)
Diferenciación Celular/genética , Linaje de la Célula/genética , Cresta Neural/crecimiento & desarrollo , Pez Cebra/crecimiento & desarrollo , Animales , Linaje de la Célula/fisiología , Transición Epitelial-Mesenquimal/genética , Regulación del Desarrollo de la Expresión Génica/genética , Melanocitos/metabolismo , Pigmentación/genética , Pez Cebra/genética , Proteínas de Pez Cebra/genéticaRESUMEN
Cerebellar stellate cells (CSCs) are spontaneously active, tonically firing (5-30 Hz), inhibitory interneurons that synapse onto Purkinje cells. We previously analyzed the excitability properties of CSCs, focusing on four key features: type I excitability, non-monotonic first-spike latency, switching in responsiveness and runup (i.e., temporal increase in excitability during whole-cell configuration). In this study, we extend this analysis by using whole-cell configuration to show that these neurons can also burst when treated with certain pharmacological agents separately or jointly. Indeed, treatment with 4-Aminopyridine (4-AP), a partial blocker of delayed rectifier and A-type K+ channels, at low doses induces a bursting profile in CSCs significantly different than that produced at high doses or when it is applied at low doses but with cadmium (Cd2+), a blocker of high voltage-activated (HVA) Ca2+ channels. By expanding a previously revised Hodgkin-Huxley type model, through the inclusion of Ca2+-activated K+ (K(Ca)) and HVA currents, we explain how these bursts are generated and what their underlying dynamics are. Specifically, we demonstrate that the expanded model preserves the four excitability features of CSCs, as well as captures their bursting patterns induced by 4-AP and Cd2+. Model investigation reveals that 4-AP is potentiating HVA, inducing square-wave bursting at low doses and pseudo-plateau bursting at high doses, whereas Cd2+ is potentiating K(Ca), inducing pseudo-plateau bursting when applied in combination with low doses of 4-AP. Using bifurcation analysis, we show that spike adding in square-wave bursts is non-sequential when gradually changing HVA and K(Ca) maximum conductances, delayed Hopf is responsible for generating the plateau segment within the active phase of pseudo-plateau bursts, and bursting can become "chaotic" when HVA and K(Ca) maximum conductances are made low and high, respectively. These results highlight the secondary effects of the drugs applied and suggest that CSCs have all the ingredients needed for bursting.
Asunto(s)
4-Aminopiridina/farmacología , Potenciales de Acción/efectos de los fármacos , Cadmio/farmacología , Bloqueadores de los Canales de Calcio/farmacología , Cerebelo/efectos de los fármacos , Bloqueadores de los Canales de Potasio/farmacología , Células de Purkinje/efectos de los fármacos , 4-Aminopiridina/administración & dosificación , Animales , Cadmio/administración & dosificación , Cerebelo/citología , Cerebelo/fisiología , Relación Dosis-Respuesta a Droga , Ratones , Ratones Endogámicos C57BL , Modelos Biológicos , Técnicas de Placa-Clamp , Células de Purkinje/fisiologíaRESUMEN
Cerebellar stellate cells form inhibitory synapses with Purkinje cells, the sole output of the cerebellum. Upon stimulation by a pair of varying inhibitory and fixed excitatory presynaptic inputs, these cells do not respond to excitation (i.e., do not generate an action potential) when the magnitude of the inhibition is within a given range, but they do respond outside this range. We previously used a revised Hodgkin-Huxley type of model to study the nonmonotonic first-spike latency of these cells and their temporal increase in excitability in whole cell configuration (termed run-up). Here, we recompute these latency profiles using the same model by adapting an efficient computational technique, the two-point boundary value problem, that is combined with the continuation method. We then extend the study to investigate how switching in responsiveness, upon stimulation with presynaptic inputs, manifests itself in the context of run-up. A three-dimensional reduced model is initially derived from the original six-dimensional model and then analyzed to demonstrate that both models exhibit type 1 excitability possessing a saddle-node on an invariant cycle (SNIC) bifurcation when varying the amplitude of Iapp. Using slow-fast analysis, we show that the original model possesses three equilibria lying at the intersection of the critical manifold of the fast subsystem and the nullcline of the slow variable hA (the inactivation of the A-type K+ channel), the middle equilibrium is of saddle type with two-dimensional stable manifold (computed from the reduced model) acting as a boundary between the responsive and non-responsive regimes, and the (ghost of) SNIC is formed when the hA-nullcline is (nearly) tangential to the critical manifold. We also show that the slow dynamics associated with (the ghost of) the SNIC and the lower stable branch of the critical manifold are responsible for generating the nonmonotonic first-spike latency. These results thus provide important insight into the complex dynamics of stellate cells.