Cooperative differentiation through clustering in multicellular populations

Solitary death in coupled limit cycle oscillators with higher-order interactions

Coupled limit cycle oscillators with pairwise interactions are known to depict phase transitions from an oscillatory state to amplitude or oscillation death. This Research Letter introduces a scheme to incorporate higher-order interactions which cannot be decomposed into pairwise interactions and investigates the dynamical evolution of Stuart-Landau oscillators under the impression of such a coupling. We discover an oscillator death state through a first-order (explosive) phase transition in which a single, coupling-dependent stable death state away from the origin exists in isolation without being accompanied by any other stable state usually existing for pairwise couplings. We call such a state a solitary death state. Contrary to widespread subcritical Hopf bifurcation, here we report homoclinic bifurcation as an origin of the explosive death state. Moreover, this explosive transition to the death state is preceded by a surge in amplitude and followed by a revival of the oscillations. The analytical value of the critical coupling strength for the solitary death state agrees with the simulation results. Finally, we point out the resemblance of the results with different dynamical states associated with epileptic seizures.

Self-Regulated Symmetry Breaking Model for Stem Cell Differentiation

In conventional disorder-order phase transitions, a system shifts from a highly symmetric state, where all states are equally accessible (disorder) to a less symmetric state with a limited number of available states (order). This transition may occur by varying a control parameter that represents the intrinsic noise of the system. It has been suggested that stem cell differentiation can be considered as a sequence of such symmetry-breaking events. Pluripotent stem cells, with their capacity to develop into any specialized cell type, are considered highly symmetric systems. In contrast, differentiated cells have lower symmetry, as they can only carry out a limited number of functions. For this hypothesis to be valid, differentiation should emerge collectively in stem cell populations. Additionally, such populations must have the ability to self-regulate intrinsic noise and navigate through a critical point where spontaneous symmetry breaking (differentiation) occurs. This study presents a mean-field model for stem cell populations that considers the interplay of cell-cell cooperativity, cell-to-cell variability, and finite-size effects. By introducing a feedback mechanism to control intrinsic noise, the model can self-tune through different bifurcation points, facilitating spontaneous symmetry breaking. Standard stability analysis showed that the system can potentially differentiate into several cell types mathematically expressed as stable nodes and limit cycles. The existence of a Hopf bifurcation in our model is discussed in light of stem cell differentiation.

Transition from inhomogeneous limit cycles to oscillation death in nonlinear oscillators with similarity-dependent coupling

We consider a system of coupled nonlinear oscillators in which the interaction is modulated by a measure of the similarity between the oscillators. Such a coupling is common in treating spatially mobile dynamical systems where the interaction is distance dependent or in resonance-enhanced interactions, for instance. For a system of Stuart-Landau oscillators coupled in this manner, we observe a novel route to oscillation death via a Hopf bifurcation. The individual oscillators are confined to inhomogeneous limit cycles initially and are damped to different fixed points after the bifurcation. Analytical and numerical results are presented for this case, while numerical results are presented for coupled Rössler and Sprott oscillators.

Three faces of biofilms: a microbial lifestyle, a nascent multicellular organism, and an incubator for diversity

Biofilms are organised heterogeneous assemblages of microbial cells that are encased within a self-produced matrix. Current estimates suggest that up to 80% of bacterial and archaeal cells reside in biofilms. Since biofilms are the main mode of microbial life, understanding their biology and functions is critical, especially as controlling biofilm growth is essential in industrial, infrastructure and medical contexts. Here we discuss biofilms both as collections of individual cells, and as multicellular biological individuals, and introduce the concept of biofilms as unique incubators of diversity for the microbial world.

Effects of propagation delay in coupled oscillators under direct-indirect coupling: Theory and experiment

Propagation delay arises in a coupling channel due to the finite propagation speed of signals and the dispersive nature of the channel. In this paper, we study the effects of propagation delay that appears in the indirect coupling path of direct (diffusive)-indirect (environmental) coupled oscillators. In sharp contrast to the direct coupled oscillators where propagation delay induces amplitude death, we show that in the case of direct-indirect coupling, even a small propagation delay is conducive to an oscillatory behavior. It is well known that simultaneous application of direct and indirect coupling is the general mechanism for amplitude death. However, here we show that the presence of propagation delay hinders the death state and helps the revival of oscillation. We demonstrate our results by considering chaotic time-delayed oscillators and FitzHugh-Nagumo oscillators. We use linear stability analysis to derive the explicit conditions for the onset of oscillation from the death state. We also verify the robustness of our results in an electronic hardware level experiment. Our study reveals that the effect of time delay on the dynamics of coupled oscillators is coupling function dependent and, therefore, highly non-trivial.

Robustness and timing of cellular differentiation through population-based symmetry breaking

During mammalian development and homeostasis, cells often transition from a multilineage primed state to one of several differentiated cell types that are marked by the expression of mutually exclusive genetic markers. These observations have been classically explained by single-cell multistability as the dynamical basis of differentiation, where robust cell-type proportioning relies on pre-existing cell-to-cell differences. We propose a conceptually different dynamical mechanism in which cell types emerge and are maintained collectively by cell-cell communication as a novel inhomogeneous state of the coupled system. Differentiation can be triggered by cell number increase as the population grows in size, through organisation of the initial homogeneous population before the symmetry-breaking bifurcation point. Robust proportioning and reliable recovery of the differentiated cell types following a perturbation is an inherent feature of the inhomogeneous state that is collectively maintained. This dynamical mechanism is valid for systems with steady-state or oscillatory single-cell dynamics. Therefore, our results suggest that timing and subsequent differentiation in robust cell-type proportions can emerge from the cooperative behaviour of growing cell populations during development.

Emergence of multistability and strongly asymmetric collective modes in two quorum sensing coupled identical ring oscillators

The simplest ring oscillator is made from three strongly nonlinear elements repressing each other unidirectionally, resulting in the emergence of a limit cycle. A popular implementation of this scheme uses repressor genes in bacteria, creating the synthetic genetic oscillator known as the Repressilator. We consider the main collective modes produced when two identical Repressilators are mean-field-coupled via the quorum-sensing mechanism. In-phase and anti-phase oscillations of the coupled oscillators emerge from two Andronov-Hopf bifurcations of the homogeneous steady state. Using the rate of the repressor's production and the value of coupling strength as the bifurcation parameters, we performed one-parameter continuations of limit cycles and two-parameter continuations of their bifurcations to show how bifurcations of the in-phase and anti-phase oscillations influence the dynamical behaviors for this system. Pitchfork bifurcation of the unstable in-phase cycle leads to the creation of novel inhomogeneous limit cycles with very different amplitudes, in contrast to the well-known asymmetrical limit cycles arising from oscillation death. The Neimark-Sacker bifurcation of the anti-phase cycle determines the border of an island in two-parameter space containing almost all the interesting regimes including the set of resonant limit cycles, the area with stable inhomogeneous cycle, and very large areas with chaotic regimes resulting from torus destruction and period doubling of resonant cycles and inhomogeneous cycles. We discuss the structure of the chaos skeleton to show the role of inhomogeneous cycles in its formation. Many regions of multistability and transitions between regimes are presented. These results provide new insights into the coupling-dependent mechanisms of multistability and collective regime symmetry breaking in populations of identical multidimensional oscillators.

Revival and death of oscillation under mean-field coupling: Interplay of intrinsic and extrinsic filtering

Mean-field diffusive coupling was known to induce the phenomenon of quenching of oscillations even in identical systems, where the standard diffusive coupling (without mean-field) fails to do so [Phys. Rev. E 89, 052912 (2014)PLEEE81539-375510.1103/PhysRevE.89.052912]. In particular, the mean-field diffusive coupling facilitates the transition from amplitude to oscillation death states and the onset of a nontrivial amplitude death state via a subcritical pitchfork bifurcation. In this paper, we show that an adaptive coupling using a low-pass filter in both the intrinsic and extrinsic variables in the coupling is capable of inducing the counterintuitive phenomenon of reviving of oscillations from the death states induced by the mean-field coupling. In particular, even a weak filtering of the extrinsic (intrinsic) variable in the mean-field coupling facilitates the onset of revival (quenching) of oscillations, whereas a strong filtering of the extrinsic (intrinsic) variable results in quenching (revival) of oscillations. Our results reveal that the degree of filtering plays a predominant role in determining the effect of filtering in the extrinsic or intrinsic variables, thereby engineering the dynamics as desired. We also extend the analysis to networks of mean-field coupled limit-cycle and chaotic oscillators along with the low-pass filters to illustrate the generic nature of our results. Finally, we demonstrate the observed dynamical transition experimentally to elucidate the robustness of our results despite the presence of inherent parameter fluctuations and noise.

Coexistence of oscillation and quenching states: Effect of low-pass active filtering in coupled oscillators

Effects of a low-pass active filter (LPAF) on the transition processes from oscillation quenching to asymmetrical oscillation are explored for diffusively coupled oscillators. The low-pass filter part and the active part of LPAF exhibit different effects on the dynamics of these coupled oscillators. With the amplifying active part only, LPAF keeps the coupled oscillators staying in a nontrivial amplitude death (NTAD) and oscillation state. However, the additional filter is beneficial to induce a transition from a symmetrical oscillation death to an asymmetrical oscillation death and then to an asymmetrical oscillation state which is oscillating with different amplitudes for two oscillators. Asymmetrical oscillation state is coexisting with a synchronous oscillation state for properly presented parameters. With the attenuating active part only, LPAF keeps the coupled oscillators in rich oscillation quenching states such as amplitude death (AD), symmetrical oscillation death (OD), and NTAD. The additional filter tends to enlarge the AD domains but to shrink the symmetrical OD domains by increasing the areas of the coexistence of the oscillation state and the symmetrical OD state. The stronger filter effects enlarge the basin of the symmetrical OD state which is coexisting with the synchronous oscillation state. Moreover, the effects of the filter are general in globally coupled oscillators. Our results are important for understanding and controlling the multistability of coupled systems.

Integrated Information as a Measure of Cognitive Processes in Coupled Genetic Repressilators

Intercellular communication and its coordination allow cells to exhibit multistability as a form of adaptation. This conveys information processing from intracellular signaling networks enabling self-organization between other cells, typically involving mechanisms associated with cognitive systems. How information is integrated in a functional manner and its relationship with the different cell fates is still unclear. In parallel, drawn originally from studies on neuroscience, integrated information proposes an approach to quantify the balance between integration and differentiation in the causal dynamics among the elements in any interacting system. In this work, such an approach is considered to study the dynamical complexity in a genetic network of repressilators coupled by quorum sensing. Several attractors under different conditions are identified and related to proposed measures of integrated information to have an insight into the collective interaction and functional differentiation in cells. This research particularly accounts for the open question about the coding and information transmission in genetic systems.

Quenching and revival of oscillations induced by coupling through adaptive variables

An adaptive coupling based on a low-pass filter (LPF) is proposed to manipulate dynamic activity of diffusively coupled dynamical systems. A theoretical analysis shows that tracking either the external or internal signal in the coupling via a LPF gives rise to distinctly different ways of regulating the rhythmicity of the coupled systems. When the external signals of the coupling are attenuated by a LPF, the macroscopic oscillations of the coupled system are quenched due to the emergence of amplitude or oscillation death. If the internal signals of the coupling are further filtered by a LPF, amplitude and oscillation deaths are effectively revoked to restore dynamic behaviors. The applicability of this approach is demonstrated in laboratory experiments of coupled oscillatory electrochemical reactions by inducing coupling through LPFs. Our study provides additional insight into (ar)rhythmogenesis in diffusively coupled systems.

Networks of coupled oscillators: From phase to amplitude chimeras

We show that amplitude-mediated phase chimeras and amplitude chimeras can occur in the same network of nonlocally coupled identical oscillators. These are two different partial synchronization patterns, where spatially coherent domains coexist with incoherent domains and coherence/incoherence referring to both amplitude and phase or only the amplitude of the oscillators, respectively. By changing the coupling strength, the two types of chimera patterns can be induced. We find numerically that the amplitude chimeras are not short-living transients but can have a long lifetime. Also, we observe variants of the amplitude chimeras with quasiperiodic temporal oscillations. We provide a qualitative explanation of the observed phenomena in the light of symmetry breaking bifurcation scenarios. We believe that this study will shed light on the connection between two disparate chimera states having different symmetry-breaking properties.

Transition from homogeneous to inhomogeneous limit cycles: Effect of local filtering in coupled oscillators

We report an interesting symmetry-breaking transition in coupled identical oscillators, namely, the continuous transition from homogeneous to inhomogeneous limit cycle oscillations. The observed transition is the oscillatory analog of the Turing-type symmetry-breaking transition from amplitude death (i.e., stable homogeneous steady state) to oscillation death (i.e., stable inhomogeneous steady state). This novel transition occurs in the parametric zone of occurrence of rhythmogenesis and oscillation death as a consequence of the presence of local filtering in the coupling path. We consider paradigmatic oscillators, such as Stuart-Landau and van der Pol oscillators, under mean-field coupling with low-pass or all-pass filtered self-feedback and through a rigorous bifurcation analysis we explore the genesis of this transition. Further, we experimentally demonstrate the observed transition, which establishes its robustness in the presence of parameter fluctuations and noise.

The impact of propagation and processing delays on amplitude and oscillation deaths in the presence of symmetry-breaking coupling

We numerically investigate the impacts of both propagation and processing delays on the emergences of amplitude death (AD) and oscillation death (OD) in one system of two Stuart-Landau oscillators with symmetry-breaking coupling. In either the absence of or the presence of propagation delay, the processing delay destabilizes both AD and OD by revoking the stability of the stable homogenous and inhomogenous steady states. In the AD to OD transition, the processing delay destabilizes first OD from large values of coupling strength until its stable regime completely disappears and then AD from both the upper and lower bounds of the stable coupling interval. Our numerical study sheds new insight lights on the understanding of nontrivial effects of time delays on dynamic activity of coupled nonlinear systems.

A comparison of Monte Carlo-based Bayesian parameter estimation methods for stochastic models of genetic networks

We compare three state-of-the-art Bayesian inference methods for the estimation of the unknown parameters in a stochastic model of a genetic network. In particular, we introduce a stochastic version of the paradigmatic synthetic multicellular clock model proposed by Ullner et al., 2007. By introducing dynamical noise in the model and assuming that the partial observations of the system are contaminated by additive noise, we enable a principled mechanism to represent experimental uncertainties in the synthesis of the multicellular system and pave the way for the design of probabilistic methods for the estimation of any unknowns in the model. Within this setup, we tackle the Bayesian estimation of a subset of the model parameters. Specifically, we compare three Monte Carlo based numerical methods for the approximation of the posterior probability density function of the unknown parameters given a set of partial and noisy observations of the system. The schemes we assess are the particle Metropolis-Hastings (PMH) algorithm, the nonlinear population Monte Carlo (NPMC) method and the approximate Bayesian computation sequential Monte Carlo (ABC-SMC) scheme. We present an extensive numerical simulation study, which shows that while the three techniques can effectively solve the problem there are significant differences both in estimation accuracy and computational efficiency.

Explosive death induced by mean-field diffusion in identical oscillators

We report the occurrence of an explosive death transition for the first time in an ensemble of identical limit cycle and chaotic oscillators coupled via mean–field diffusion. In both systems, the variation of the normalized amplitude with the coupling strength exhibits an abrupt and irreversible transition to death state from an oscillatory state and this first order phase transition to death state is independent of the size of the system. This transition is quite general and has been found in all the coupled systems where in–phase oscillations co–exist with a coupling dependent homogeneous steady state. The backward transition point for this phase transition has been calculated using linear stability analysis which is in complete agreement with the numerics.

Basin stability measure of different steady states in coupled oscillators

In this report, we investigate the stabilization of saddle fixed points in coupled oscillators where individual oscillators exhibit the saddle fixed points. The coupled oscillators may have two structurally different types of suppressed states, namely amplitude death and oscillation death. The stabilization of saddle equilibrium point refers to the amplitude death state where oscillations are ceased and all the oscillators converge to the single stable steady state via inverse pitchfork bifurcation. Due to multistability features of oscillation death states, linear stability theory fails to analyze the stability of such states analytically, so we quantify all the states by basin stability measurement which is an universal nonlocal nonlinear concept and it interplays with the volume of basins of attractions. We also observe multi-clustered oscillation death states in a random network and measure them using basin stability framework. To explore such phenomena we choose a network of coupled Duffing-Holmes and Lorenz oscillators which are interacting through mean-field coupling. We investigate how basin stability for different steady states depends on mean-field density and coupling strength. We also analytically derive stability conditions for different steady states and confirm by rigorous bifurcation analysis.

Cell signaling as a cognitive process

Cellular identity as defined through morphology and function emerges from intracellular signaling networks that communicate between cells. Based on recursive interactions within and among these intracellular networks, dynamical solutions in terms of biochemical behavior are generated that can differ from those in isolated cells. In this way, cellular heterogeneity in tissues can be established, implying that cell identity is not intrinsically predetermined by the genetic code but is rather dynamically maintained in a cognitive manner. We address how to experimentally measure the flow of information in intracellular biochemical networks and demonstrate that even simple causality motifs can give rise to rich, context-dependent dynamic behavior. The concept how intercellular communication can result in novel dynamical solutions is applied to provide a contextual perspective on cell differentiation and tumorigenesis.

Flexible dynamics of two quorum-sensing coupled repressilators

Genetic oscillators play important roles in cell life regulation. The regulatory efficiency usually depends strongly on the emergence of stable collective dynamic modes, which requires designing the interactions between genetic networks. We investigate the dynamics of two identical synthetic genetic repressilators coupled by an additional plasmid which implements quorum sensing (QS) in each network thereby supporting global coupling. In a basic genetic ring oscillator network in which three genes inhibit each other in unidirectional manner, QS stimulates the transcriptional activity of chosen genes providing for competition between inhibitory and stimulatory activities localized in those genes. The "promoter strength", the Hill cooperativity coefficient of transcription repression, and the coupling strength, i.e., parameters controlling the basic rates of genetic reactions, were chosen for extensive bifurcation analysis. The results are presented as a map of dynamic regimes. We found that the remarkable multistability of the antiphase limit cycle and stable homogeneous and inhomogeneous steady states exists over broad ranges of control parameters. We studied the antiphase limit cycle stability and the evolution of irregular oscillatory regimes in the parameter areas where the antiphase cycle loses stability. In these regions we observed developing complex oscillations, collective chaos, and multistability between regular limit cycles and complex oscillations over uncommonly large intervals of coupling strength. QS coupling stimulates the appearance of intrachaotic periodic windows with spatially symmetric and asymmetric partial limit cycles which, in turn, change the type of chaos from a simple antiphase character into chaos composed of pieces of the trajectories having alternating polarity. The very rich dynamics discovered in the system of two identical simple ring oscillators may serve as a possible background for biological phenotypic diversification, as well as a stimulator to search for similar coupling in oscillator arrays in other areas of nature, e.g., in neurobiology, ecology, climatology, etc.

Generation and cessation of oscillations: Interplay of excitability and dispersal in a class of ecosystems

We investigate the complex spatiotemporal dynamics of an ecological network with species dispersal mediated via a mean-field coupling. The local dynamics of the network are governed by the Truscott-Brindley model, which is an important ecological model showing excitability. Our results focus on the interplay of excitability and dispersal by always considering that the individual nodes are in their (excitable) steady states. In contrast to the previous studies, we not only observe the dispersal induced generation of oscillation but also report two distinct mechanisms of cessation of oscillations, namely, amplitude and oscillation death. We show that the dispersal between the nodes influences the intrinsic dynamics of the system resulting in multiple oscillatory dynamics such as period-1 and period-2 limit cycles. We also show the existence of multi-cluster states, which has much relevance and importance in ecology.

Experimental demonstration of revival of oscillations from death in coupled nonlinear oscillators

We experimentally demonstrate that a processing delay, a finite response time, in the coupling can revoke the stability of the stable steady states, thereby facilitating the revival of oscillations in the same parameter space where the coupled oscillators suffered the quenching of oscillation. This phenomenon of reviving of oscillations is demonstrated using two different prototype electronic circuits. Further, the analytical critical curves corroborate that the spread of the parameter space with stable steady state is diminished continuously by increasing the processing delay. Finally, the death state is completely wiped off above a threshold value by switching the stability of the stable steady state to retrieve sustained oscillations in the same parameter space. The underlying dynamical mechanism responsible for the decrease in the spread of the stable steady states and the eventual reviving of oscillation as a function of the processing delay is explained using analytical results.

Multiscale modeling of tumor growth induced by circadian rhythm disruption in epithelial tissue

We propose a multiscale chemo-mechanical model of cancer tumor development in epithelial tissue. The model is based on the transformation of normal cells into a cancerous state triggered by a local failure of spatial synchronization of the circadian rhythm. The model includes mechanical interactions and a chemical signal exchange between neighboring cells, as well as a division of cells and intercalation that allows for modification of the respective parameters following transformation into the cancerous state. The numerical simulations reproduce different dephasing patterns--spiral waves and quasistationary clustering, with the latter being conducive to cancer formation. Modification of mechanical properties reproduces a distinct behavior of invasive and localized carcinoma.

Mixed-mode oscillation suppression states in coupled oscillators

We report a collective dynamical state, namely the mixed-mode oscillation suppression state where the steady states of the state variables of a system of coupled oscillators show heterogeneous behaviors. We identify two variants of it: The first one is a mixed-mode death (MMD) state, which is an interesting oscillation death state, where a set of variables show dissimilar values, while the rest arrive at a common value. In the second mixed death state, bistable and monostable nontrivial homogeneous steady states appear simultaneously to a different set of variables (we refer to it as the MNAD state). We find these states in the paradigmatic chaotic Lorenz system and Lorenz-like system under generic coupling schemes. We identify that while the reflection symmetry breaking is responsible for the MNAD state, the breaking of both the reflection and translational symmetries result in the MMD state. Using a rigorous bifurcation analysis we establish the occurrence of the MMD and MNAD states, and map their transition routes in parameter space. Moreover, we report experimental observation of the MMD and MNAD states that supports our theoretical results. We believe that this study will broaden our understanding of oscillation suppression states; subsequently, it may have applications in many real physical systems, such as laser and geomagnetic systems, whose mathematical models mimic the Lorenz system.

Revival of oscillation from mean-field-induced death: Theory and experiment

The revival of oscillation and maintaining rhythmicity in a network of coupled oscillators offer an open challenge to researchers as the cessation of oscillation often leads to a fatal system degradation and an irrecoverable malfunctioning in many physical, biological, and physiological systems. Recently a general technique of restoration of rhythmicity in diffusively coupled networks of nonlinear oscillators has been proposed in Zou et al. [Nat. Commun. 6, 7709 (2015)], where it is shown that a proper feedback parameter that controls the rate of diffusion can effectively revive oscillation from an oscillation suppressed state. In this paper we show that the mean-field diffusive coupling, which can suppress oscillation even in a network of identical oscillators, can be modified in order to revoke the cessation of oscillation induced by it. Using a rigorous bifurcation analysis we show that, unlike other diffusive coupling schemes, here one has two control parameters, namely the density of the mean-field and the feedback parameter that can be controlled to revive oscillation from a death state. We demonstrate that an appropriate choice of density of the mean field is capable of inducing rhythmicity even in the presence of complete diffusion, which is a unique feature of this mean-field coupling that is not available in other coupling schemes. Finally, we report the experimental observation of revival of oscillation from the mean-field-induced oscillation suppression state that supports our theoretical results.

Spatial coexistence of synchronized oscillation and death: A chimeralike state

We report an interesting spatiotemporal state, namely the chimeralike incongruous coexistence of synchronized oscillation and stable steady state (CSOD) in a network of nonlocally coupled oscillators. Unlike the chimera and chimera death state, in the CSOD state identical oscillators are self-organized into two coexisting spatially separated domains: In one domain neighboring oscillators show synchronized oscillation and in another domain the neighboring oscillators randomly populate either a synchronized oscillating state or a stable steady state (we call it a death state). We consider a realistic ecological network and show that the interplay of nonlocality and coupling strength results in two routes to the CSOD state: One is from a coexisting mixed state of amplitude chimera and death, and another one is from a globally synchronized state. We provide a qualitative explanation of the origin of this state. We further explore the importance of this study in ecology that gives insight into the relationship between spatial synchrony and global extinction of species. We believe this study will improve our understanding of chimera and chimeralike states.

Dispersal-induced synchrony, temporal stability, and clustering in a mean-field coupled Rosenzweig-MacArthur model

In spatial ecology, dispersal among a set of spatially separated habitats, named as metapopulation, preserves the diversity and persistence by interconnecting the local populations. Understanding the effects of several variants of dispersion in metapopulation dynamics and to identify the factors which promote population synchrony and population stability are important in ecology. In this paper, we consider the mean-field dispersion among the habitats in a network and study the collective dynamics of the spatially extended system. Using the Rosenzweig-MacArthur model for individual patches, we show that the population synchrony and temporal stability, which are believed to be of conflicting outcomes of dispersion, can be simultaneously achieved by oscillation quenching mechanisms. Particularly, we explore the more natural coupling configuration where the rates of dispersal of different habitats are disparate. We show that asymmetry in dispersal rate plays a crucial role in determining inhomogeneity in an otherwise homogeneous metapopulation. We further identify an unusual emergent state in the network, namely, a multi-branch clustered inhomogeneous steady state, which arises due to the intrinsic parameter mismatch among the patches. We believe that the present study will shed light on the cooperative behavior of spatially structured ecosystems.

The effect of mechanical stimulation on mineralization in differentiating osteoblasts in collagen-I scaffolds

Developing a viable and functional bone scaffold in vitro that is capable of surviving and bearing mechanical load in vivo requires an understanding of the cell biology of osteoprogenitor cells, particularly how they are influenced by mechanical stimulation during cell differentiation and maturation. In this study, mechanical load was applied using a modified FlexCell plate to impart confined compression to collagen-I scaffolds seeded with undifferentiated murine embryonic stem cells. The activity, presence, and expression of osteoblast-cadherin (OB-Cad) and connexin-43, as well as various pluripotent and osteogenic markers were examined at 5-30 days of differentiation as cells were stimulated to differentiate to osteoblasts with and without applied mechanical load. Fluorescence recovery after photobleaching, immunofluorescence, viability, von Kossa, and real-time polymerase chain reaction assessments revealed that mechanical prestimulation of this cell-seeded scaffold altered the expression of OB-Cad and connexin-43 and resulted in significant differences in the structure and organization of mineralization present in the collagen matrix. Specifically, cells in gels that were loaded for 40 h after 5 days of differentiation and then left to fully differentiate for 30 days produced a highly structured honeycomb-shaped mineralization in the matrix; an outcome that was previously shown to be indicative of late osteoblast/early osteocyte activity. This study highlights the potential of mechanical load to accelerate differentiation and enhance osteoblast communication and function during the differentiation process, and highlights a time point of cell differentiation within this scaffold to apply load in order to most effectively transduce a mechanical signal.

Pluripotency, Differentiation, and Reprogramming: A Gene Expression Dynamics Model with Epigenetic Feedback Regulation

Embryonic stem cells exhibit pluripotency: they can differentiate into all types of somatic cells. Pluripotent genes such as Oct4 and Nanog are activated in the pluripotent state, and their expression decreases during cell differentiation. Inversely, expression of differentiation genes such as Gata6 and Gata4 is promoted during differentiation. The gene regulatory network controlling the expression of these genes has been described, and slower-scale epigenetic modifications have been uncovered. Although the differentiation of pluripotent stem cells is normally irreversible, reprogramming of cells can be experimentally manipulated to regain pluripotency via overexpression of certain genes. Despite these experimental advances, the dynamics and mechanisms of differentiation and reprogramming are not yet fully understood. Based on recent experimental findings, we constructed a simple gene regulatory network including pluripotent and differentiation genes, and we demonstrated the existence of pluripotent and differentiated states from the resultant dynamical-systems model. Two differentiation mechanisms, interaction-induced switching from an expression oscillatory state and noise-assisted transition between bistable stationary states, were tested in the model. The former was found to be relevant to the differentiation process. We also introduced variables representing epigenetic modifications, which controlled the threshold for gene expression. By assuming positive feedback between expression levels and the epigenetic variables, we observed differentiation in expression dynamics. Additionally, with numerical reprogramming experiments for differentiated cells, we showed that pluripotency was recovered in cells by imposing overexpression of two pluripotent genes and external factors to control expression of differentiation genes. Interestingly, these factors were consistent with the four Yamanaka factors, Oct4, Sox2, Klf4, and Myc, which were necessary for the establishment of induced pluripotent stem cells. These results, based on a gene regulatory network and expression dynamics, contribute to our wider understanding of pluripotency, differentiation, and reprogramming of cells, and they provide a fresh viewpoint on robustness and control during development.

Characterization of pluripotent states, in which cells can both self-renew and differentiate, and the irreversible loss of pluripotency are important research areas in developmental biology. In particular, an understanding of these processes is essential to the reprogramming of cells for biomedical applications, i.e., the experimental recovery of pluripotency in differentiated cells. Based on recent advances in dynamical-systems theory for gene expression, we propose a gene-regulatory-network model consisting of several pluripotent and differentiation genes. Our results show that cellular-state transition to differentiated cell types occurs as the number of cells increases, beginning with the pluripotent state and oscillatory expression of pluripotent genes. Cell-cell signaling mediates the differentiation process with robustness to noise, while epigenetic modifications affecting gene expression dynamics fix the cellular state. These modifications ensure the cellular state to be protected against external perturbation, but they also work as an epigenetic barrier to recovery of pluripotency. We show that overexpression of several genes leads to the reprogramming of cells, consistent with the methods for establishing induced pluripotent stem cells. Our model, which involves the inter-relationship between gene expression dynamics and epigenetic modifications, improves our basic understanding of cell differentiation and reprogramming.

Synchronization and phase redistribution in self-replicating populations of coupled oscillators and excitable elements

We study the dynamics of phase synchronization in growing populations of discrete phase oscillatory systems when the division process is coupled to the distribution of oscillator phases. Using mean-field theory, linear-stability analysis, and numerical simulations, we demonstrate that coupling between population growth and synchrony can lead to a wide range of dynamical behavior, including extinction of synchronized oscillations, the emergence of asynchronous states with unequal state (phase) distributions, bistability between oscillatory and asynchronous states or between two asynchronous states, a switch between continuous (supercritical) and discontinuous (subcritical) transitions, and modulation of the frequency of bulk oscillations.

Oscillator death induced by amplitude-dependent coupling in repulsively coupled oscillators

The effects of amplitude-dependent coupling on oscillator death (OD) are investigated for two repulsively coupled Lorenz oscillators. Based on numerical simulations, it is shown that as constraint strengths on the amplitude-dependent coupling change, an oscillatory state may undergo a transition to an OD state. The parameter regimes of the OD domain are theoretically determined, which coincide well with the numerical results. An electronic circuit is set up to exhibit the transition process to the OD state with an amplitude-dependent coupling. These findings may have practical importance on chaos control and oscillation depression.

Mean-field dispersion-induced spatial synchrony, oscillation and amplitude death, and temporal stability in an ecological model

One of the most important issues in spatial ecology is to understand how spatial synchrony and dispersal-induced stability interact. In the existing studies it is shown that dispersion among identical patches results in spatial synchrony; on the other hand, the combination of spatial heterogeneity and dispersion is necessary for dispersal-induced stability (or temporal stability). Population synchrony and temporal stability are thus often thought of as conflicting outcomes of dispersion. In contrast to the general belief, in this present study we show that mean-field dispersion is conducive to both spatial synchrony and dispersal-induced stability even in identical patches. This simultaneous occurrence of rather conflicting phenomena is governed by the suppression of oscillation states, namely amplitude death (AD) and oscillation death (OD). These states emerge through spatial synchrony of the oscillating patches in the strong-coupling strength. We present an interpretation of the mean-field diffusive coupling in the context of ecology and identify that, with increasing mean-field density, an open ecosystem transforms into a closed ecosystem. We report on the occurrence of OD in an ecological model and explain its significance. Using a detailed bifurcation analysis we show that, depending on the mortality rate and carrying capacity, the system shows either AD or both AD and OD. We also show that the results remain qualitatively the same for a network of oscillators. We identify a new transition scenario between the same type of oscillation suppression states whose geneses differ. In the parameter-mismatched case, we further report on the direct transition from OD to AD through a transcritical bifurcation. We believe that this study will lead to a proper interpretation of AD and OD in ecology, which may be important for the conservation and management of several communities in ecosystems.

Transitions among the diverse oscillation quenching states induced by the interplay of direct and indirect coupling

We report the transitions among different oscillation quenching states induced by the interplay of diffusive (direct) coupling and environmental (indirect) coupling in coupled identical oscillators. This coupling scheme was introduced by Resmi et al. [Phys. Rev. E 84, 046212 (2011)] as a general scheme to induce amplitude death (AD) in nonlinear oscillators. Using a detailed bifurcation analysis we show that, in addition to AD, which actually occurs only in a small region of parameter space, this coupling scheme can induce other oscillation quenching states, namely oscillation death (OD) and a novel nontrvial AD (NAD) state, which is a nonzero bistable homogeneous steady state; more importantly, this coupling scheme mediates a transition from the AD state to the OD state and a new transition from the AD state to the NAD state. We identify diverse routes to the NAD state and map all the transition scenarios in the parameter space for periodic oscillators. Finally, we present the first experimental evidence of oscillation quenching states and their transitions induced by the interplay of direct and indirect coupling.

Transition from amplitude to oscillation death under mean-field diffusive coupling

We study the transition from the amplitude death (AD) to the oscillation death (OD) state in limit-cycle oscillators coupled through mean-field diffusion. We show that this coupling scheme can induce an important transition from AD to OD even in identical limit cycle oscillators. We identify a parameter region where OD and a nontrivial AD (NTAD) state coexist. This NTAD state is unique in comparison with AD owing to the fact that it is created by a subcritical pitchfork bifurcation and parameter mismatch does not support this state, but destroys it. We extend our study to a network of mean-field coupled oscillators to show that the transition scenario is preserved and the oscillators form a two-cluster state.

Diverse routes of transition from amplitude to oscillation death in coupled oscillators under additional repulsive links

We report the existence of diverse routes of transition from amplitude death to oscillation death in three different diffusively coupled systems, which are perturbed by a symmetry breaking repulsive coupling link. For limit-cycle systems the transition is through a pitchfork bifurcation, as has been noted before, but in chaotic systems it can be through a saddle-node or a transcritical bifurcation depending on the nature of the underlying dynamics of the individual systems. The diversity of the routes and their dependence on the complex dynamics of the coupled systems not only broadens our understanding of this important phenomenon but can lead to potentially new practical applications.

Parameter estimation methods for chaotic intercellular networks

We have investigated simulation-based techniques for parameter estimation in chaotic intercellular networks. The proposed methodology combines a synchronization–based framework for parameter estimation in coupled chaotic systems with some state–of–the–art computational inference methods borrowed from the field of computational statistics. The first method is a stochastic optimization algorithm, known as accelerated random search method, and the other two techniques are based on approximate Bayesian computation. The latter is a general methodology for non–parametric inference that can be applied to practically any system of interest. The first method based on approximate Bayesian computation is a Markov Chain Monte Carlo scheme that generates a series of random parameter realizations for which a low synchronization error is guaranteed. We show that accurate parameter estimates can be obtained by averaging over these realizations. The second ABC–based technique is a Sequential Monte Carlo scheme. The algorithm generates a sequence of “populations”, i.e., sets of randomly generated parameter values, where the members of a certain population attain a synchronization error that is lesser than the error attained by members of the previous population. Again, we show that accurate estimates can be obtained by averaging over the parameter values in the last population of the sequence. We have analysed how effective these methods are from a computational perspective. For the numerical simulations we have considered a network that consists of two modified repressilators with identical parameters, coupled by the fast diffusion of the autoinducer across the cell membranes.

Generalizing the transition from amplitude to oscillation death in coupled oscillators

Amplitude death (AD) and oscillation death (OD) are two structurally different oscillation quenching types in coupled nonlinear oscillators. The transition from AD to OD has been recently realized due to the interplay between heterogeneity and coupling strength [A. Koseska et al., Phys. Rev. Lett. 111, 024103 (2013)]. We identify here the transition from AD to OD in nonlinear oscillators with couplings of distinct natures. It is demonstrated that the presence of time delay in the coupling cannot induce such a transition in identical oscillators, but it can indeed facilitate its occurrence with a low degree of heterogeneity. Moreover, it is further shown that the AD to OD transition is reliably observed in identical oscillators with dynamic and conjugate couplings. The coexistence of AD and OD and rich stable OD configurations after the transition are revealed, which are of great significance for potential applications in physics, biology, and control studies.

Minimal model for stem-cell differentiation

To explain the differentiation of stem cells in terms of dynamical systems theory, models of interacting cells with intracellular protein expression dynamics are analyzed and simulated. Simulations were carried out for all possible protein expression networks consisting of two genes under cell-cell interactions mediated by the diffusion of a protein. Networks that show cell differentiation are extracted and two forms of symmetric differentiation based on Turing's mechanism and asymmetric differentiation are identified. In the latter network, the intracellular protein levels show oscillatory dynamics at a single-cell level, while cell-to-cell synchronicity of the oscillation is lost with an increase in the number of cells. Differentiation to a fixed-point-type behavior follows with a further increase in the number of cells. The cell type with oscillatory dynamics corresponds to a stem cell that can both proliferate and differentiate, while the latter fixed-point type only proliferates. This differentiation is analyzed as a saddle-node bifurcation on an invariant circle, while the number ratio of each cell type is shown to be robust against perturbations due to self-consistent determination of the effective bifurcation parameter as a result of the cell-cell interaction. Complex cell differentiation is designed by combing these simple two-gene networks. The generality of the present differentiation mechanism, as well as its biological relevance, is discussed.

Transition from amplitude to oscillation death via Turing bifurcation

Coupled oscillators are shown to experience two structurally different oscillation quenching types: amplitude death (AD) and oscillation death (OD). We demonstrate that both AD and OD can occur in one system and find that the transition between them underlies a classical, Turing-type bifurcation, providing a clear classification of these significantly different dynamical regimes. The implications of obtaining a homogeneous (AD) or inhomogeneous (OD) steady state, as well as their significance for physical and biological applications and control studies, are also pointed out.

Electronic implementation of a repressilator with quorum sensing feedback

We investigate the dynamics of a synthetic genetic repressilator with quorum sensing feedback. In a basic genetic ring oscillator network in which three genes inhibit each other in unidirectional manner, an additional quorum sensing feedback loop stimulates the activity of a chosen gene providing competition between inhibitory and stimulatory activities localized in that gene. Numerical simulations show several interesting dynamics, multi-stability of limit cycle with stable steady-state, multi-stability of different stable steady-states, limit cycle with period-doubling and reverse period-doubling, and infinite period bifurcation transitions for both increasing and decreasing strength of quorum sensing feedback. We design an electronic analog of the repressilator with quorum sensing feedback and reproduce, in experiment, the numerically predicted dynamical features of the system. Noise amplification near infinite period bifurcation is also observed. An important feature of the electronic design is the accessibility and control of the important system parameters.

Dimensionality reduction of bistable biological systems

Time hierarchies, arising as a result of interactions between system's components, represent a ubiquitous property of dynamical biological systems. In addition, biological systems have been attributed switch-like properties modulating the response to various stimuli across different organisms and environmental conditions. Therefore, establishing the interplay between these features of system dynamics renders itself a challenging question of practical interest in biology. Existing methods are suitable for systems with one stable steady state employed as a well-defined reference. In such systems, the characterization of the time hierarchies has already been used for determining the components that contribute to the dynamics of biological systems. However, the application of these methods to bistable nonlinear systems is impeded due to their inherent dependence on the reference state, which in this case is no longer unique. Here, we extend the applicability of the reference-state analysis by proposing, analyzing, and applying a novel method, which allows investigation of the time hierarchies in systems exhibiting bistability. The proposed method is in turn used in identifying the components, other than reactions, which determine the systemic dynamical properties. We demonstrate that in biological systems of varying levels of complexity and spanning different biological levels, the method can be effectively employed for model simplification while ensuring preservation of qualitative dynamical properties (i.e., bistability). Finally, by establishing a connection between techniques from nonlinear dynamics and multivariate statistics, the proposed approach provides the basis for extending reference-based analysis to bistable systems.

Characterizing an ensemble of interacting oscillators: the mean-field variability index

We introduce a way of characterizing an ensemble of interacting oscillators in terms of their mean-field variability index κ, a dimensionless parameter defined as the variance of the oscillators' mean field r divided by the mean square of r. Based on the assumption that the overall mean field is the sum of a very large number of oscillators, each giving a small contribution to the total signal, we show that κ depends on the mutual interactions between the oscillators, independently of their number or spectral properties. For purely random phasors, or a noninteracting ensemble of oscillators, κ converges on 0.215. Interactions push κ in different directions: lower where there is interoscillator phase coherence, tending to zero for complete phase synchronization, or higher for amplitude synchronization or intermittent synchronization. We calculate κ for several different cases to illustrate its utility, using both numerically simulated data and electroencephalograph signals from the brains of human subjects while awake, while anesthetized, and while undergoing an epileptic fit.

Inhomogeneous stationary and oscillatory regimes in coupled chaotic oscillators

The dynamics of linearly coupled identical Lorenz and Pikovsky-Rabinovich oscillators are explored numerically and theoretically. We concentrate on the study of inhomogeneous stable steady states ("oscillation death (OD)" phenomenon) and accompanying periodic and chaotic regimes that emerge at an appropriate choice of the coupling matrix. The parameters, for which OD occurs, are determined by stability analysis of the chosen steady state. Three model-specific types of transitions to and from OD are observed: (1) a sharp transition to OD from a nonsymmetric chaotic attractor containing random intervals of synchronous chaos; (2) transition to OD from the symmetry-breaking chaotic regime created by negative coupling; (3) supercritical bifurcation of OD into inhomogeneous limit cycles and further evolution of the system to inhomogeneous chaotic regimes that coexist with complete synchronous chaos. These results may fill a gap in the understanding of the mechanism of OD in coupled chaotic systems.

Inference of time-evolving coupled dynamical systems in the presence of noise

A new method is introduced for analysis of interactions between time-dependent coupled oscillators, based on the signals they generate. It distinguishes unsynchronized dynamics from noise-induced phase slips and enables the evolution of the coupling functions and other parameters to be followed. It is based on phase dynamics, with Bayesian inference of the time-evolving parameters achieved by shaping the prior densities to incorporate knowledge of previous samples. The method is tested numerically and applied to reveal and quantify the time-varying nature of cardiorespiratory interactions.

"Quorum sensing" generated multistability and chaos in a synthetic genetic oscillator

We model the dynamics of the synthetic genetic oscillator Repressilator equipped with quorum sensing. In addition to a circuit of 3 genes repressing each other in a unidirectional manner, the model includes a phase-repulsive type of the coupling module implemented as the production of a small diffusive molecule-autoinducer (AI). We show that the autoinducer (which stimulates the transcription of a target gene) is responsible for the disappearance of the limit cycle (LC) through the infinite period bifurcation and the formation of a stable steady state (SSS) for sufficiently large values of the transcription rate. We found conditions for hysteresis between the limit cycle and the stable steady state. The parameters' region of the hysteresis is determined by the mRNA to protein lifetime ratio and by the level of transcription-stimulating activity of the AI. In addition to hysteresis, increasing AI-dependent stimulation of transcription may lead to the complex dynamic behavior which is characterized by the appearance of several branches on the bifurcation continuation, containing different regular limit cycles, as well as a chaotic regime. The multistability which is manifested as the coexistence between the stable steady state, limit cycles, and chaos seems to be a novel type of the dynamics for the ring oscillator with the added quorum sensing positive feedback.

Excitatory and inhibitory coupling in a one-dimensional array of Belousov-Zhabotinsky micro-oscillators: theory

We study numerically the behavior of one-dimensional arrays of aqueous droplets containing the oscillatory Belousov-Zhabotinsky reaction. Droplets are separated by an oil phase that allows coupling between neighboring droplets via two species: an inhibitor, Br(2), and an activator, HBrO(2). Excitatory coupling alone (through the activator) generates in-phase oscillations and/or "waves," while inhibitory coupling alone (through Br(2)) gives rise to antiphase oscillations, Turing patterns, and their combinations. The simultaneous presence of excitatory and inhibitory coupling leads to a large number of new spatiotemporal patterns, including some that exhibit very complex behavior. Analysis of a simple model allows us to simulate patterns resembling those observed experimentally under similar conditions and to elucidate the contributions of droplet and gap sizes, activator and inhibitor partition coefficients, and malonic acid concentration to the coupling strengths.

Analysing dynamical behavior of cellular networks via stochastic bifurcations

The dynamical structure of genetic networks determines the occurrence of various biological mechanisms, such as cellular differentiation. However, the question of how cellular diversity evolves in relation to the inherent stochasticity and intercellular communication remains still to be understood. Here, we define a concept of stochastic bifurcations suitable to investigate the dynamical structure of genetic networks, and show that under stochastic influence, the expression of given proteins of interest is defined via the probability distribution of the phase variable, representing one of the genes constituting the system. Moreover, we show that under changing stochastic conditions, the probabilities of expressing certain concentration values are different, leading to different functionality of the cells, and thus to differentiation of the cells in the various types.

Communication-induced multistability and multirhythmicity in a synthetic multicellular system

Traditionally, the main role of cell-to-cell communication was thought of as synchronizing a population of cells, thereby coordinating cellular behavior. Here we show that cell density, which quantifies cellular communication, can induce multistability and multirhythmicity in a synthetic multicellular system, where individual oscillators are a combination of repressillator and hysteresis-based oscillators and are coupled through a quorum-sensing mechanism. Specifically, for moderately small cell densities, the coupled system can exhibit multistability including stable homogenous and inhomogeneous steady states. For moderately large cell densities, it has the potential to generate multirhythmicity including multimode oscillations such as an in-phase periodic solution, antiphase periodic solution, asymmetric periodic solution, mixed-mode oscillations, coexistence of periodic orbits of several different modes, and bursting oscillations such as periodic bursting, torus quasiperiodic bursting, and chaotic bursting. Such versatility of cell-to-cell communication would be beneficial for cells or organisms to live in diversely changeable environments.