The unprecedented light curves of the Kepler space telescope document how the brightness of some stars pulsates at primary and secondary frequencies whose ratios are near the golden mean, the most irrational number. A nonlinear dynamical system driven by an irrational ratio of frequencies generically exhibits a strange but nonchaotic attractor. For Kepler 's "golden" stars, we present evidence of the first observation of strange nonchaotic dynamics in nature outside the laboratory. This discovery could aid the classification and detailed modeling of variable stars.

Stochasticity in gene expression impacts the dynamics and functions of gene regulatory circuits. Intrinsic noises, including those that are caused by low copy number of molecules and transcriptional bursting, are usually studied by stochastic simulations. However, the role of extrinsic factors, such as cell-to-cell variability and heterogeneity in the microenvironment, is still elusive. To evaluate the effects of both the intrinsic and extrinsic noises, we develop a method, named sRACIPE, by integrating stochastic analysis with random circuit perturbation (RACIPE) method. RACIPE uniquely generates and analyzes an ensemble of models with random kinetic parameters. Previously, we have shown that the gene expression from random models form robust and functionally related clusters. In sRACIPE we further develop two stochastic simulation schemes, aiming to reduce the computational cost without sacrificing the convergence of statistics. One scheme uses constant noise to capture the basins of attraction, and the other one uses simulated annealing to detect the stability of states. By testing the methods on several synthetic gene regulatory circuits and an epithelial–mesenchymal transition network in squamous cell carcinoma, we demonstrate that sRACIPE can interpret the experimental observations from single-cell gene expression data. We observe that parametric variation (the spread of parameters around a median value) increases the spread of the gene expression clusters, whereas high noise merges the states. Our approach quantifies the robustness of a gene circuit in the presence of noise and sheds light on a new mechanism of noise-induced hybrid states. We have implemented sRACIPE as an R package.

We study the stability of the synchronized state in time-varying complex networks using the concept of basin stability, which is a nonlocal and nonlinear measure of stability that can be easily applied to high-dimensional systems [P. J. Menck, J. Heitzig, N. Marwan, and J. Kurths, Nature Phys. 9, 89 (2013)]. The time-varying character is included by stochastically rewiring each link with the average frequency f. We find that the time taken to reach synchronization is lowered and the stability range of the synchronized state increases considerably in dynamic networks. Further we uncover that small-world networks are much more sensitive to link changes than random ones, with the time-varying character of the network having a significant effect at much lower rewiring frequencies. At very high rewiring frequencies, random networks perform better than small-world networks and the synchronized state is stable over a much wider window of coupling strengths. Lastly we show that the stability range of the synchronized state may be quite different for small and large perturbations, and so the linear stability analysis and the basin stability criterion provide complementary indicators of stability.

Stochasticity in gene expression impacts the dynamics and functions of gene regulatory circuits. Intrinsic noises, including those that are caused by low copy number of molecules and transcriptional bursting, are usually studied by stochastic analysis methods, such as Gillespie algorithm and Langevin simulation. However, the role of extrinsic factors, such as cell-to-cell variability and heterogeneity in microenvironment, is still elusive. To evaluate the effects of both intrinsic and extrinsic noises, we develop a new method, named sRACIPE, by integrating stochastic analysis with random circuit perturbation (RACIPE) method. Unlike traditional methods, RACIPE generates and analyzes an ensemble of mathematical models with random kinetic parameters. Previously, we have shown that the gene expression from random models form robust and functionally related clusters. Under the framework of this randomization-based approach, here we develop two stochastic simulation schemes, aiming to reduce the computational cost without sacrificing the convergence of statistics. One scheme uses constant noise to capture the basins of attraction, and the other one uses simulated annealing to detect the stability of states. By testing the methods on several gene regulatory circuits, we found that high noise, but not large parameter variation, merges clusters together. Our approach quantifies the robustness of a gene circuit in the presence of noise and sheds light on a new mechanism of noise induced hybrid states. We have implemented sRACIPE into a freely available R package.