Effect of network architecture on synchronization and entrainment properties of the circadian oscillations in the suprachiasmatic nucleus

The Functional Connectome Mediating Circadian Synchrony in the Suprachiasmatic Nucleus

Circadian rhythms in mammals arise from the spatiotemporal synchronization of ~20,000 neuronal clocks in the Suprachiasmatic Nucleus (SCN). While anatomical, molecular, and genetic approaches have revealed diverse cell types and signaling mechanisms, the network wiring that enables SCN cells to communicate and synchronize remains unclear. To overcome the challenges of revealing functional connectivity from fixed tissue, we developed MITE (Mutual Information & Transfer Entropy), an information theory approach that infers directed cell-cell connections with high fidelity. By analyzing 3447 hours of continuously recorded clock gene expression from 9011 cells in 17 mice, we found that the functional connectome of SCN was highly conserved bilaterally and across mice, sparse, and organized into a dorsomedial and a ventrolateral module. While most connections were local, we discovered long-range connections from ventral cells to cells in both the ventral and dorsal SCN. Based on their functional connectivity, SCN cells can be characterized as circadian signal generators, broadcasters, sinks, or bridges. For example, a subset of VIP neurons acts as hubs that generate circadian signals critical to synchronize daily rhythms across the SCN neural network. Simulations of the experimentally inferred SCN networks recapitulated the stereotypical dorsal-to-ventral wave of daily PER2 expression and ability to spontaneously synchronize, revealing that SCN emergent dynamics are sculpted by cell-cell connectivity. We conclude that MITE provides a powerful method to infer functional connectomes, and that the conserved architecture of cell-cell connections mediates circadian synchrony across space and time in the mammalian SCN.

Impact of Network Topology on Neural Synchrony in a Model of the Subthalamic Nucleus-Globus Pallidus Circuit

Synchronous neural oscillations within the beta frequency range are observed across the parkinsonian basal ganglia network, including within the subthalamic nucleus (STN) - globus pallidus (GPe) subcircuit. The emergence of pathological synchrony in Parkinson's disease is often attributed to changes in neural properties or connection strength, and less often to the network topology, i.e. the structural arrangement of connections between neurons. This study investigates the relationship between network structure and neural synchrony in a model of the STN-GPe circuit comprised of conductance-based spiking neurons. Changes in net synaptic input were controlled for through a synaptic scaling rule, which facilitated separation of the effects of network structure from net synaptic input. Five topologies were examined as structures for the STN-GPe circuit: Watts-Strogatz, preferential attachment, spatial, stochastic block, k-regular random. Beta band synchrony generally increased as the number of connections increased, however the exact relationship was topology specific. Varying the wiring pattern while maintaining a constant number of connections caused network synchrony to be enhanced or suppressed, demonstrating the ability of purely structural changes to alter synchrony. This relationship was well-captured by the algebraic connectivity of the network, the second smallest eigenvalue of the network's Laplacian matrix. The structure-synchrony relationship was further investigated in a network model designed to emulate the action selection role of the STN-GPe circuit. It was found that increasing the number of connections and/or the overlap of action selection channels could lead to a rapid transition to synchrony, which was also predicted by the algebraic connectivity.

A minimal model of peripheral clocks reveals differential circadian re-entrainment in aging

The mammalian circadian system comprises a network of endogenous oscillators, spanning from the central clock in the brain to peripheral clocks in other organs. These clocks are tightly coordinated to orchestrate rhythmic physiological and behavioral functions. Dysregulation of these rhythms is a hallmark of aging, yet it remains unclear how age-related changes lead to more easily disrupted circadian rhythms. Using a two-population model of coupled oscillators that integrates the central clock and the peripheral clocks, we derive simple mean-field equations that can capture many aspects of the rich behavior found in the mammalian circadian system. We focus on three age-associated effects that have been posited to contribute to circadian misalignment: attenuated input from the sympathetic pathway, reduced responsiveness to light, and a decline in the expression of neurotransmitters. We find that the first two factors can significantly impede re-entrainment of the clocks following perturbation, while a weaker coupling within the central clock does not affect the recovery rate. Moreover, using our minimal model, we demonstrate the potential of using the feed-fast cycle as an effective intervention to accelerate circadian re-entrainment. These results highlight the importance of peripheral clocks in regulating the circadian rhythm and provide fresh insights into the complex interplay between aging and the resilience of the circadian system.

Systematic analysis of negative and positive feedback loops for robustness and temperature compensation in circadian rhythms

Temperature compensation and robustness to biological noise are two key characteristics of the circadian clock. These features allow the circadian pacemaker to maintain a steady oscillation in a wide range of environmental conditions. The presence of a time-delayed negative feedback loop in the regulatory network generates autonomous circadian oscillations in eukaryotic systems. In comparison, the circadian clock of cyanobacteria is controlled by a strong positive feedback loop. Positive feedback loops with substrate depletion can also generate oscillations, inspiring other circadian clock models. What makes a circadian oscillatory network robust to extrinsic noise is unclear. We investigated four basic circadian oscillators with negative, positive, and combinations of positive and negative feedback loops to explore network features necessary for circadian clock resilience. We discovered that the negative feedback loop system performs the best in compensating temperature changes. We also show that a positive feedback loop can reduce extrinsic noise in periods of circadian oscillators, while intrinsic noise is reduced by negative feedback loops.

Generation and Disruption of Circadian Rhythms in the Suprachiasmatic Nucleus: A Core-Shell Model

We focus our research on how the core-shell organization controls behavior of the suprachiasmatic nucleus (SCN), how the core and shell are synchronized to the environment, what impact they have on the behavior of the SCN under different lighting conditions, and what mechanisms disrupt synchronization. To this end, we use a reduced Kuramoto model, with parameters inferred from experimental observations and calibrated for mice, and perform a detailed comparison between the model and experimental data under light-dark (LD), dark-dark (DD), and light-light (LL) conditions. The operating limits of free-running and entrained SCN activity under symmetric LD cycles are analyzed, with particular focus on the phenomena of anticipation and dissociation. Results reveal that the core-shell organization of the SCN enables anticipation of future events over circadian cycles. The model predicts the emergence of a second (dissociated) rhythm for large and small LD periods. Our results are in good qualitative and quantitative agreement with experimental observations of circadian dissociation. We further describe SCN activity under LL conditions and show that our model satisfies Aschoff’s first rule, according to which the endogenous free-running circadian period observed under complete darkness will shorten in diurnal animals and lengthen in nocturnal animals under constant light. Our results strongly suggest that the Kuramoto model captures essential features of synchronization and entrainment in the SCN. Moreover, our approach is easily extendible to an arbitrary number of groups, with dynamics described by explicit equations for the group phase and synchronization index. Viewed together, the reduced Kuramoto model presents itself as a useful tool for exploring open problems in the study of circadian rhythms, one that can account for evolving views of the circadian system’s organization, including peripheral clocks and inter-hemispheric interaction, and can be translated to other nocturnal and diurnal animals, including humans.

Deciphering clock cell network morphology within the biological master clock, suprachiasmatic nucleus: From the perspective of circadian wave dynamics

The biological master clock, suprachiasmatic nucleus (of rat and mouse), is composed of ~10,000 clock cells which are heterogeneous with respect to their circadian periods. Despite this inhomogeneity, an intact SCN maintains a very good degree of circadian phase (time) coherence which is vital for sustaining various circadian rhythmic activities, and it is supposedly achieved by not just one but a few different cell-to-cell coupling mechanisms, among which action potential (AP)-mediated connectivity is known to be essential. But, due to technical difficulties and limitations in experiments, so far very little information is available about the morphology of the connectivity at a cellular scale. Building upon this limited amount of information, here we exhaustively and systematically explore a large pool (~25,000) of various network morphologies to come up with some plausible network features of SCN networks. All candidates under consideration reflect an experimentally obtained 'indegree distribution' as well as a 'physical range distribution of afferent clock cells.' Then, importantly, with a set of multitude criteria based on the properties of SCN circadian phase waves in extrinsically perturbed as well as in their natural states, we select out appropriate model networks: Some important measures are, 1) level of phase dispersal and direction of wave propagation, 2) phase-resetting ability of the model networks subject to external circadian forcing, and 3) decay rate of perturbation induced "phase-singularities." The successful, realistic networks have several common features: 1) "indegree" and "outdegree" should have a positive correlation; 2) the cells in the SCN ventrolateral region (core) have a much larger total degree than that of the dorsal medial region (shell); 3) The number of intra-core edges is about 7.5 times that of intra-shell edges; and 4) the distance probability density function for the afferent connections fits well to a beta function. We believe that these newly identified network features would be a useful guide for future explorations on the very much unknown AP-mediated clock cell connectome within the SCN.

Time-keeping and decision-making in living cells: Part I

To survive and reproduce, a cell must process information from its environment and its own internal state and respond accordingly, in terms of metabolic activity, gene expression, movement, growth, division and differentiation. These signal–response decisions are made by complex networks of interacting genes and proteins, which function as biochemical switches and clocks, and other recognizable information-processing circuitry. This theme issue of Interface Focus (in two parts) brings together articles on time-keeping and decision-making in living cells—work that uses precise mathematical modelling of underlying molecular regulatory networks to understand important features of cell physiology. Part I focuses on time-keeping: mechanisms and dynamics of biological oscillators and modes of synchronization and entrainment of oscillators, with special attention to circadian clocks.

Network rewiring and plasticity promotes synchronization of suprachiasmatic nucleus neurons

In mammals, circadian rhythms throughout the body are orchestrated by the master clock in the hypothalamic suprachiasmatic nucleus (SCN), where SCN neurons are coupled with neurotransmitters to generate a uniform circadian rhythm. How the SCN circadian rhythm is so robust and flexible is, however, unclear. In this paper, we propose a temporal SCN network model and investigate the effects of dynamical rewiring and flexible coupling due to synaptic plasticity on the synchronization of the neural network in SCN. In networks consisting of simple Poincaré oscillators and complex circadian clocks, we found that dynamical rewiring and coupling plasticity enhance the synchronization in inhomogeneous networks. We verified the effect of enhanced synchronization in different architectures of random, scale-free, and small-world networks. A simple mean-field analysis for synchronization in plastic networks is proposed. Intuitively, the synchronization is greatly enhanced because both the random rewiring and coupling plasticity in the heterogeneous network have effectively increased the coupling strength in the whole network. Our results suggest that a proper network model for the master SCN circadian rhythm needs to take into account the effects of dynamical changes in topology and plasticity in neuron interactions that could help the brain to generate a robust circadian rhythm for the whole body.

Light-induced synchronization of the SCN coupled oscillators and implications for entraining the HPA axis

The suprachiasmatic nucleus (SCN) synchronizes the physiological rhythms to the external light-dark cycle and tunes the dynamics of circadian rhythms to photoperiod fluctuations. Changes in the neuronal network topologies are suggested to cause adaptation of the SCN in different photoperiods, resulting in the broader phase distribution of neuron activities in long photoperiods (LP) compared to short photoperiods (SP). Regulated by the SCN output, the level of glucocorticoids is elevated in short photoperiod, which is associated with peak disease incidence. The underlying coupling mechanisms of the SCN and the interplay between the SCN and the HPA axis have yet to be fully elucidated. In this work, we propose a mathematical model including a multiple-cellular SCN compartment and the HPA axis to investigate the properties of the circadian timing system under photoperiod changes. Our model predicts that the probability-dependent network is more energy-efficient than the distance-dependent network. Coupling the SCN network by intra-subpopulation and inter-subpopulation forces, we identified the negative correlation between robustness and plasticity of the oscillatory network. The HPA rhythms were predicted to be strongly entrained to the SCN rhythms with a pro-inflammatory high-amplitude glucocorticoid profile under SP. The fast temporal topology switch of the SCN network was predicted to enhance synchronization when the synchronization is not complete. These synchronization and circadian dynamics alterations might govern the seasonal variation of disease incidence and its symptom severity.

Cut the noise or couple up: Coordinating circadian and synthetic clocks

Circadian clocks are important to much of life on Earth and are of inherent interest to humanity, implicated in fields ranging from agriculture and ecology to developmental biology and medicine. New techniques show that it is not simply the presence of clocks, but coordination between them that is critical for complex physiological processes across the kingdoms of life. Recent years have also seen impressive advances in synthetic biology to the point where parallels can be drawn between synthetic biological and circadian oscillators. This review will emphasize theoretical and experimental studies that have revealed a fascinating dichotomy of coupling and heterogeneity among circadian clocks. We will also consolidate the fields of chronobiology and synthetic biology, discussing key design principles of their respective oscillators.

Light entrainment of the SCN circadian clock and implications for personalized alterations of corticosterone rhythms in shift work and jet lag

The suprachiasmatic nucleus (SCN) functions as the central pacemaker aligning physiological and behavioral oscillations to day/night (activity/inactivity) transitions. The light signal entrains the molecular clock of the photo-sensitive ventrolateral (VL) core of the SCN which in turn entrains the dorsomedial (DM) shell via the neurotransmitter vasoactive intestinal polypeptide (VIP). The shell converts the VIP rhythmic signals to circadian oscillations of arginine vasopressin (AVP), which eventually act as a neurotransmitter signal entraining the hypothalamic–pituitary–adrenal (HPA) axis, leading to robust circadian secretion of glucocorticoids. In this work, we discuss a semi-mechanistic mathematical model that reflects the essential hierarchical structure of the photic signal transduction from the SCN to the HPA axis. By incorporating the interactions across the core, the shell, and the HPA axis, we investigate how these coupled systems synchronize leading to robust circadian oscillations. Our model predicts the existence of personalized synchronization strategies that enable the maintenance of homeostatic rhythms while allowing for differential responses to transient and permanent light schedule changes. We simulated different behavioral situations leading to perturbed rhythmicity, performed a detailed computational analysis of the dynamic response of the system under varying light schedules, and determined that (1) significant interindividual diversity and flexibility characterize adaptation to varying light schedules; (2) an individual’s tolerances to jet lag and alternating shift work are positively correlated, while the tolerances to jet lag and transient shift work are negatively correlated, which indicates trade-offs in an individual’s ability to maintain physiological rhythmicity; (3) weak light sensitivity leads to the reduction of circadian flexibility, implying that light therapy can be a potential approach to address shift work and jet lag related disorders. Finally, we developed a map of the impact of the synchronization within the SCN and between the SCN and the HPA axis as it relates to the emergence of circadian flexibility.

Network Structure of the Master Clock Is Important for Its Primary Function

A master clock located in the suprachiasmatic nucleus (SCN) regulates the circadian rhythm of physiological and behavioral activities in mammals. The SCN has two main functions in the regulation: an endogenous clock produces the endogenous rhythmic signal in body rhythms, and a calibrator synchronizes the body rhythms to the external light-dark cycle. These two functions have been determined to depend on either the dynamic behaviors of individual neurons or the whole SCN neuronal network. In this review, we first introduce possible network structures for the SCN, as revealed by time series analysis from real experimental data. It was found that the SCN network is heterogeneous and sparse, that is, the average shortest path length is very short, some nodes are hubs with large node degrees but most nodes have small node degrees, and the average node degree of the network is small. Secondly, the effects of the SCN network structure on the SCN function are reviewed based on mathematical models of the SCN network. It was found that robust rhythms with large amplitudes, a high synchronization between SCN neurons and a large entrainment ability exists mainly in small-world and scale-free type networks, but not other types. We conclude that the SCN most probably is an efficient small-world type or scale-free type network, which drives SCN function.

Modelling the functional roles of synaptic and extra-synaptic γ-aminobutyric acid receptor dynamics in circadian timekeeping

In the suprachiasmatic nucleus (SCN), γ-aminobutyric acid (GABA) is a primary neurotransmitter. GABA can signal through two types of GABAA receptor subunits, often referred to as synaptic GABAA (gamma subunit) and extra-synaptic GABAA (delta subunit). To test the functional roles of these distinct GABAA in regulating circadian rhythms, we developed a multicellular SCN model where we could separately compare the effects of manipulating GABA neurotransmitter or receptor dynamics. Our model predicted that blocking GABA signalling modestly increased synchrony among circadian cells, consistent with published SCN pharmacology. Conversely, the model predicted that lowering GABAA receptor density reduced firing rate, circadian cell fraction, amplitude and synchrony among individual neurons. When we tested these predictions, we found that the knockdown of delta GABAA reduced the amplitude and synchrony of clock gene expression among cells in SCN explants. The model further predicted that increasing gamma GABAA densities could enhance synchrony, as opposed to increasing delta GABAA densities. Overall, our model reveals how blocking GABAA receptors can modestly increase synchrony, while increasing the relative density of gamma over delta subunits can dramatically increase synchrony. We hypothesize that increased gamma GABAA density in the winter could underlie the tighter phase relationships among SCN cells.

Analysis of input-induced oscillations using the isostable coordinate framework

Many reduced order modeling techniques for oscillatory dynamical systems are only applicable when the underlying system admits a stable periodic orbit in the absence of input. By contrast, very few reduction frameworks can be applied when the oscillations themselves are induced by coupling or other exogenous inputs. In this work, the behavior of such input-induced oscillations is considered. By leveraging the isostable coordinate framework, a high-accuracy reduced set of equations can be identified and used to predict coupling-induced bifurcations that precipitate stable oscillations. Subsequent analysis is performed to predict the steady state phase-locking relationships. Input-induced oscillations are considered for two classes of coupled dynamical systems. For the first, stable fixed points of systems with parameters near Hopf bifurcations are considered so that the salient dynamical features can be captured using an asymptotic expansion of the isostable coordinate dynamics. For the second, an adaptive phase-amplitude reduction framework is used to analyze input-induced oscillations that emerge in excitable systems. Examples with relevance to circadian and neural physiology are provided that highlight the utility of the proposed techniques.

Evolutionary Constraints on Connectivity Patterns in the Mammalian Suprachiasmatic Nucleus

The mammalian suprachiasmatic nucleus (SCN) comprises about 20,000 interconnected oscillatory neurons that create and maintain a robust circadian signal which matches to external light cues. Here, we use an evolutionary game theoretic framework to explore how evolutionary constraints can influence the synchronization of the system under various assumptions on the connection topology, contributing to the understanding of the structure of interneuron connectivity. Our basic model represents the SCN as a network of agents each with two properties-a phase and a flag that determines if it communicates with its neighbors or not. Communication comes at a cost to the agent, but synchronization of phases with its neighbors bears a benefit. Earlier work shows that when we have "all-to-all" connectivity, where every agent potentially communicates with every other agent, there is often a simple trade-off that leads to complete communication and synchronization of the system: the benefit must be greater than twice the cost. This trade-off for all-to-all connectivity gives us a baseline to compare to when looking at other topologies. Using simulations, we compare three plausible topologies to the all-to-all case, finding that convergence to synchronous dynamics occurs in all considered topologies under similar benefit and cost trade-offs. Consequently, sparser, less biologically costly topologies are reasonable evolutionary outcomes for organisms that develop a synchronizable oscillatory network. Our simulations also shed light on constraints imposed by the time scale on which we observe the SCN to arise in mammals. We find two conditions that allow for a synchronizable system to arise in relatively few generations. First, the benefits of connectivity must outweigh the cost of facilitating the connectivity in the network. Second, the game at the core of the model needs to be more cooperative than antagonistic games such as the Prisoner's Dilemma. These results again imply that evolutionary pressure may have driven the system towards sparser topologies, as they are less costly to create and maintain. Last, our simulations indicate that models based on the mutualism game fare the best in uptake of communication and synchronization compared to more antagonistic games such as the Prisoner's Dilemma.

Topology-induced dynamics in a network of synthetic oscillators with piecewise affine approximation

In synthetic biology approaches, minimal systems are used to reproduce complex molecular mechanisms that appear in the core functioning of multi-cellular organisms. In this paper, we study a piecewise affine model of a synthetic two-gene oscillator and prove existence and stability of a periodic solution for all parameters in a given region. Motivated by the synchronization of circadian clocks in a cluster of cells, we next consider a network of N identical oscillators under diffusive coupling to investigate the effect of the topology of interactions in the network's dynamics. Our results show that both all-to-all and one-to-all coupling topologies may introduce new stable steady states in addition to the expected periodic orbit. Both topologies admit an upper bound on the coupling parameter that prevents the generation of new steady states. However, this upper bound is independent of the number of oscillators in the network and less conservative for the one-to-all topology.

Astrocytic Modulation of Neuronal Activity in the Suprachiasmatic Nucleus: Insights from Mathematical Modeling

The suprachiasmatic nucleus (SCN) of the hypothalamus consists of a highly heterogeneous neuronal population networked together to allow precise and robust circadian timekeeping in mammals. While the critical importance of SCN neurons in regulating circadian rhythms has been extensively studied, the roles of SCN astrocytes in circadian system function are not well understood. Recent experiments have demonstrated that SCN astrocytes are circadian oscillators with the same functional clock genes as SCN neurons. Astrocytes generate rhythmic outputs that are thought to modulate neuronal activity through pre- and postsynaptic interactions. In this study, we developed an in silico multicellular model of the SCN clock to investigate the impact of astrocytes in modulating neuronal activity and affecting key clock properties such as circadian rhythmicity, period, and synchronization. The model predicted that astrocytes could alter the rhythmic activity of neurons via bidirectional interactions at tripartite synapses. Specifically, astrocyte-regulated extracellular glutamate was predicted to increase neuropeptide signaling from neurons. Consistent with experimental results, we found that astrocytes could increase the circadian period and enhance neural synchronization according to their endogenous circadian period. The impact of astrocytic modulation of circadian rhythm amplitude, period, and synchronization was predicted to be strongest when astrocytes had periods between 0 and 2 h longer than neurons. Increasing the number of neurons coupled to the astrocyte also increased its impact on period modulation and synchrony. These computational results suggest that signals that modulate astrocytic rhythms or signaling (e.g., as a function of season, age, or treatment) could cause disruptions in circadian rhythm or serve as putative therapeutic targets.

The Concept of Coupling in the Mammalian Circadian Clock Network

The circadian clock network regulates daily rhythms in mammalian physiology and behavior to optimally adapt the organism to the 24-h day/night cycle. A central pacemaker, the hypothalamic suprachiasmatic nucleus (SCN), coordinates subordinate cellular oscillators in the brain, as well as in peripheral organs to align with each other and external time. Stability and coordination of this vast network of cellular oscillators is achieved through different levels of coupling. Although coupling at the molecular level and across the SCN is well established and believed to define its function as pacemaker structure, the notion of coupling in other tissues and across the whole system is less well understood. In this review, we describe the different levels of coupling in the mammalian circadian clock system - from molecules to the whole organism. We highlight recent advances in gaining knowledge of the complex organization and function of circadian network regulation and its significance for the generation of stable but plastic intrinsic 24-h rhythms.

Disassortative Network Structure Improves the Synchronization between Neurons in the Suprachiasmatic Nucleus

In mammals, an endogenous clock located in the suprachiasmatic nucleus (SCN) of the brain regulates the circadian rhythms of physiological and behavioral activities. The SCN is composed of about 20,000 neurons that are autonomous oscillators with nonidentical intrinsic periods ranging from 22 h to 28 h. These neurons are coupled through neurotransmitters and synchronized to form a network, which produces a robust circadian rhythm of a uniform period. The neurons, which are the nodes in the network, are known to be heterogeneous in their characteristics, which is reflected in different phenotypes and different functionality. This heterogeneous nature of the nodes of the network leads to the question as to whether the structure of the SCN network is assortative or disassortative. Thus far, the disassortativity of the SCN network has not been assessed and neither have its effects on the collective behaviors of the SCN neurons. In the present study, we build a directed SCN network composed of hundreds of neurons for a single slice using the method of transfer entropy, based on the experimental data. Then, we measured the synchronization degree as well as the disassortativity coefficient of the network structure (calculated by either the out-degrees or the in-degrees of the nodes) and found that the network of the SCN is a disassortative network. Furthermore, a positive relationship is observed between the synchronization degree and disassortativity of the network, which is confirmed by simulations of our modeling. Our finding suggests that the disassortativity of the network structure plays a role in the synchronization between SCN neurons; that is, the synchronization degree increases with the increase of the disassortativity, which implies that a more heterogeneous coupling in the network of the SCN is important for proper function of the SCN.

Molecular and Cellular Networks in The Suprachiasmatic Nuclei

Why do we experience the ailments of jetlag when we travel across time zones? Why is working night-shifts so detrimental to our health? In other words, why can’t we readily choose and stick to non-24 h rhythms? Actually, our daily behavior and physiology do not simply result from the passive reaction of our organism to the external cycle of days and nights. Instead, an internal clock drives the variations in our bodily functions with a period close to 24 h, which is supposed to enhance fitness to regular and predictable changes of our natural environment. This so-called circadian clock relies on a molecular mechanism that generates rhythmicity in virtually all of our cells. However, the robustness of the circadian clock and its resilience to phase shifts emerge from the interaction between cell-autonomous oscillators within the suprachiasmatic nuclei (SCN) of the hypothalamus. Thus, managing jetlag and other circadian disorders will undoubtedly require extensive knowledge of the functional organization of SCN cell networks. Here, we review the molecular and cellular principles of circadian timekeeping, and their integration in the multi-cellular complexity of the SCN. We propose that new, in vivo imaging techniques now enable to address these questions directly in freely moving animals.

Chronobiology of limbic seizures: Potential mechanisms and prospects of chronotherapy for mesial temporal lobe epilepsy

Mesial Temporal Lobe Epilepsy (mTLE) characterized by progressive development of complex partial seizures originating from the hippocampus is the most prevalent and refractory type of epilepsy. One of the remarkable features of mTLE is the rhythmic pattern of occurrence of spontaneous seizures, implying a dependence on the endogenous clock system for seizure threshold. Conversely, circadian rhythms are affected by epilepsy too. Comprehending how the circadian system and seizures interact with each other is essential for understanding the pathophysiology of epilepsy as well as for developing innovative therapies that are efficacious for better seizure control. In this review, we confer how the temporal dysregulation of the circadian clock in the hippocampus combined with multiple uncoupled oscillators could lead to periodic seizure occurrences and comorbidities. Unraveling these associations with additional research would help in developing chronotherapy for mTLE, based on the chronobiology of spontaneous seizures. Notably, differential dosing of antiepileptic drugs over the circadian period and/or strategies that resynchronize biological rhythms may substantially improve the management of seizures in mTLE patients.

Coherency of circadian rhythms in the SCN is governed by the interplay of two coupling factors

Circadian clocks are autonomous oscillators driving daily rhythms in physiology and behavior. In mammals, a network of coupled neurons in the suprachiasmatic nucleus (SCN) is entrained to environmental light-dark cycles and orchestrates the timing of peripheral organs. In each neuron, transcriptional feedbacks generate noisy oscillations. Coupling mediated by neuropeptides such as VIP and AVP lends precision and robustness to circadian rhythms. The detailed coupling mechanisms between SCN neurons are debated. We analyze organotypic SCN slices from neonatal and adult mice in wild-type and multiple knockout conditions. Different degrees of rhythmicity are quantified by pixel-level analysis of bioluminescence data. We use empirical orthogonal functions (EOFs) to characterize spatio-temporal patterns. Simulations of coupled stochastic single cell oscillators can reproduce the diversity of observed patterns. Our combination of data analysis and modeling provides deeper insight into the enormous complexity of the data: (1) Neonatal slices are typically stronger oscillators than adult slices pointing to developmental changes of coupling. (2) Wild-type slices are completely synchronized and exhibit specific spatio-temporal patterns of phases. (3) Some slices of Cry double knockouts obey impaired synchrony that can lead to co–existing rhythms (“splitting”). (4) The loss of VIP-coupling leads to desynchronized rhythms with few residual local clusters. Additional information was extracted from co–culturing slices with rhythmic neonatal wild-type SCNs. These co–culturing experiments were simulated using external forcing terms representing VIP and AVP signaling. The rescue of rhythmicity via co–culturing lead to surprising results, since a cocktail of AVP-antagonists improved synchrony. Our modeling suggests that these counter-intuitive observations are pointing to an antagonistic action of VIP and AVP coupling. Our systematic theoretical and experimental study shows that dual coupling mechanisms can explain the astonishing complexity of spatio-temporal patterns in SCN slices.

The mammalian circadian clock is orchestrated by a network of coupled neurons. Brain slice preparations allow the analysis of coupling mechanisms mediated by neuropeptides. From bioluminescence recordings, we extract single cell characteristics such as period, amplitude and damping rate. Our data-based stochastic network model involves local coupling between cells and additional external forcing. Available experimental data guide our simulations with two distinct coupling and forcing mechanisms representing the neuropeptides VIP and AVP. We compare our simulations with experiments from neonatal and adult wild-type brain slices and multiple knockouts. Furthermore, we study co–culturing of slices with synchronized neonatal wild-type slices. The extreme complexity of the spatio-temporal patterns is quantified using empirical orthogonal functions (EOFs). The experimental reduction of AVP coupling leads to surprising observations. In double knockouts, inhibition of AVP signaling can improve synchrony, whereas, in triple knockouts, coherency is reduced. Our network modeling shows that these counter-intuitive observations can be explained by an antagonistic action of VIP and AVP signaling. The agreement of experiments and simulations suggests that quite complex spatio-temporal patterns can appear as emergent properties of oscillator networks with dual coupling mechanisms.

Protein sequestration versus Hill-type repression in circadian clock models

Circadian (∼24 h) clocks are self-sustained endogenous oscillators with which organisms keep track of daily and seasonal time. Circadian clocks frequently rely on interlocked transcriptional-translational feedback loops to generate rhythms that are robust against intrinsic and extrinsic perturbations. To investigate the dynamics and mechanisms of the intracellular feedback loops in circadian clocks, a number of mathematical models have been developed. The majority of the models use Hill functions to describe transcriptional repression in a way that is similar to the Goodwin model. Recently, a new class of models with protein sequestration-based repression has been introduced. Here, the author discusses how this new class of models differs dramatically from those based on Hill-type repression in several fundamental aspects: conditions for rhythm generation, robust network designs and the periods of coupled oscillators. Consistently, these fundamental properties of circadian clocks also differ among Neurospora, Drosophila, and mammals depending on their key transcriptional repression mechanisms (Hill-type repression or protein sequestration). Based on both theoretical and experimental studies, this review highlights the importance of careful modelling of transcriptional repression mechanisms in molecular circadian clocks.

Stabilization of dynamics of oscillatory systems by nonautonomous perturbation

Synchronization and stability under periodic oscillatory driving are well understood, but little is known about the effects of aperiodic driving, despite its abundance in nature. Here, we consider oscillators subject to driving with slowly varying frequency, and investigate both short-term and long-term stability properties. For a phase oscillator, we find that, counterintuitively, such variation is guaranteed to enlarge the Arnold tongue in parameter space. Using analytical and numerical methods that provide information on time-variable dynamical properties, we find that the growth of the Arnold tongue is specifically due to the growth of a region of intermittent synchronization where trajectories alternate between short-term stability and short-term neutral stability, giving rise to stability on average. We also present examples of higher-dimensional nonlinear oscillators where a similar stabilization phenomenon is numerically observed. Our findings help support the case that in general, deterministic nonautonomous perturbation is a very good candidate for stabilizing complex dynamics.

A Meta-Analysis Characterizing Stem-Like Gene Expression in the Suprachiasmatic Nucleus and Its Circadian Clock

Cells expressing proteins characteristic of stem cells and progenitor cells are present in the suprachiasmatic nucleus (SCN) of the adult mammalian hypothalamus. Any relationship between this distinctive feature and the master circadian clock of the SCN is unclear. Considering the lack of obvious neurogenesis in the adult SCN relative to the hippocampus and other structures that provide neurons and glia, it is possible that the SCN has partially differentiated cells that can provide neural circuit plasticity rather than ongoing neurogenesis. To test this possibility, available databases and publications were explored to identify highly expressed genes in the mouse SCN that also have known or suspected roles in cell differentiation, maintenance of stem-like states, or cell-cell interactions found in adult and embryonic stem cells and cancer stem cells. The SCN was found to have numerous genes associated with stem cell maintenance and increased motility from which we selected 25 of the most relevant genes. Over ninety percent of these stem-like genes were expressed at higher levels in the SCN than in other brain areas. Further analysis of this gene set could provide a greater understanding of how adjustments in cell contacts alter period and phase relationships of circadian rhythms. Circadian timing and its role in cancer, sleep, and metabolic disorders are likely influenced by genes selected in this study.

Dissociation between two subgroups of the suprachiasmatic nucleus affected by the number of damped oscillated neurons

In mammals, the main clock located in the suprachiasmatic nucleus (SCN) of the brain synchronizes the body rhythms to the environmental light-dark cycle. The SCN is composed of about 2×10^{4} neurons which can be classified into three oscillatory phenotypes: self-sustained oscillators, damped oscillators, and arrhythmic neurons. Exposed to an artificial external light-dark cycle with a period of 22h instead of 24h, two subgroups of the SCN can become desynchronized (dissociated). The ventrolateral (VL) subgroup receives photic input and is entrained to the external cycle and a dorsomedial (DM) subgroup oscillates with its endogenous (i.e., free running) period and is synchronized to the external light-dark cycle through coupling from the VL. In the present study, we examined the effects of damped oscillatory neurons on the dissociation between VL and DM under an external 22h cycle. We found that, with increasing numbers of damped oscillatory neurons located in the VL, the dissociation between the VL and DM emerges, but if these neurons are increasingly present in the DM the dissociation disappears. Hence, the damped oscillatory neurons in different subregions of the SCN play distinct roles in the dissociation between the two subregions of the SCN. This shows that synchrony between SCN subregions is affected by the number of damped oscillatory neurons and the location of these cells. We suggest that more knowledge on the number and the location of these cells may explain why some species do show a dissociation between the subregions and others do not, as the distribution of oscillatory types of neurons offers a plausible and novel candidate mechanism to explain heterogeneity.

The autoregulatory loop: A common mechanism of regulation of key sex determining genes in insects

Sex determination in most insects is structured as a gene cascade, wherein a primary signal is passed through a series of sex-determining genes, culminating in a downstream double-switch known as doublesex that decides the sexual fate of the embryo. From the literature available on sex determination cascades, it becomes apparent that sex determination mechanisms have evolved rapidly. The primary signal that provides the cue to determine the sex of the embryo varies remarkably, not only among taxa, but also within taxa. Furthermore, the upstream key gene in the cascade also varies between species and even among closely related species. The order Insecta alone provides examples of astoundingly complex diversity of upstream key genes in sex determination mechanisms. Besides, unlike key upstream genes, the downstream double-switch gene is alternatively spliced to form functional sex-specific isoforms. This sex-specific splicing is conserved across insect taxa. The genes involved in the sex determination cascade such as Sex-lethal (Sxl) in Drosophila melanogaster, transformer (tra) in many other dipterans, coleopterans and hymenopterans, Feminizer (fem) in Apis mellifera, and IGF-II mRNA-binding protein (Bmimp) in Bombyx mori are reported to be regulated by an autoregulatory positive feedback loop. In this review, by taking examples from various insects, we propose the hypothesis that autoregulatory loop mechanisms of sex determination might be a general strategy. We also discuss the possible reasons for the evolution of autoregulatory loops in sex determination cascades and their impact on binary developmental choices.

Regulating the Suprachiasmatic Nucleus (SCN) Circadian Clockwork: Interplay between Cell-Autonomous and Circuit-Level Mechanisms

The suprachiasmatic nucleus (SCN) is the principal circadian clock of the brain, directing daily cycles of behavior and physiology. SCN neurons contain a cell-autonomous transcription-based clockwork but, in turn, circuit-level interactions synchronize the 20,000 or so SCN neurons into a robust and coherent daily timer. Synchronization requires neuropeptide signaling, regulated by a reciprocal interdependence between the molecular clockwork and rhythmic electrical activity, which in turn depends on a daytime Na⁺ drive and nighttime K⁺ drag. Recent studies exploiting intersectional genetics have started to identify the pacemaking roles of particular neuronal groups in the SCN. They support the idea that timekeeping involves nonlinear and hierarchical computations that create and incorporate timing information through the interactions between key groups of neurons within the SCN circuit. The field is now poised to elucidate these computations, their underlying cellular mechanisms, and how the SCN clock interacts with subordinate circadian clocks across the brain to determine the timing and efficiency of the sleep-wake cycle, and how perturbations of this coherence contribute to neurological and psychiatric illness.

Mechanists Must be Holists Too! Perspectives from Circadian Biology

The pursuit of mechanistic explanations in biology has produced a great deal of knowledge about the parts, operations, and organization of mechanisms taken to be responsible for biological phenomena. Holist critics have often raised important criticisms of proposed mechanistic explanations, but until recently holists have not had alternative research strategies through which to advance explanations. This paper argues both that the results of mechanistic strategies has forced mechanists to confront ways in which whole systems affect their components and that new representational and modeling strategies are providing tools for understanding these effects of whole systems upon components. Drawing from research on the mechanism responsible for circadian rhythms in mammals, I develop two examples in which mechanistic analysis is being integrated into a more holist perspective: research revealing intercellular integration of circadian mechanisms with those involved in cell metabolism and research revealing that stable␣rhythms are dependent on how individual cells in the suprachiasmatic nucleus synchronize with each other to generate regular rhythms. Tools such as network diagramming and computational modeling are providing means to integrate mechanistic models into accounts of whole systems.

Coupling Controls the Synchrony of Clock Cells in Development and Knockouts

In mammals, a network of coupled neurons within the hypothalamus coordinates physiological rhythms with daily changes in the environment. In each neuron, delayed negative transcriptional feedbacks generate oscillations, albeit noisy and unreliable ones. Coupling mediated by diffusible neuropeptides lends precision and robustness to circadian rhythms. The double knockout of Cryptochrome Cry turns adult mice arrhythmic. But, remarkably, double knockout neonates continue to show robust oscillation much like wild-type neonates and appear to lose rhythmicity with development. We study quantitatively dispersed neurons and brain slices from wild-type and Cry double knockout mice to understand the links between single cell rhythmicity and intercellular coupling. We quantify oscillator properties of dispersed cells using nonlinear regression and study bifurcations diagrams of network models. We find that varying just three parameters-oscillator strength, strength of coupling, and timing of coupling-can reproduce experimentally observed features. In particular, modeling reveals that minor changes in timing of coupling can destroy synchronization as observed in adult slices from knockout mice.

Resynchronization of circadian oscillators and the east-west asymmetry of jet-lag

Cells in the brain's Suprachiasmatic Nucleus (SCN) are known to regulate circadian rhythms in mammals. We model synchronization of SCN cells using the forced Kuramoto model, which consists of a large population of coupled phase oscillators (modeling individual SCN cells) with heterogeneous intrinsic frequencies and external periodic forcing. Here, the periodic forcing models diurnally varying external inputs such as sunrise, sunset, and alarm clocks. We reduce the dimensionality of the system using the ansatz of Ott and Antonsen and then study the effect of a sudden change of clock phase to simulate cross-time-zone travel. We estimate model parameters from previous biological experiments. By examining the phase space dynamics of the model, we study the mechanism leading to the difference typically experienced in the severity of jet-lag resulting from eastward and westward travel.

The synchronization of neuronal oscillators determined by the directed network structure of the suprachiasmatic nucleus under different photoperiods

The main function of the principal clock located in the suprachiasmatic nucleus (SCN) of mammals is synchronizing the body rhythms to the 24 h light-dark cycle. Additionally, the SCN is able to adapt to the photoperiod of the cycle which varies among seasons. Under the long photoperiod (LP), the synchronization degree of the SCN neurons is lower than that under the photoperiod (SP). In the present study, a potential explanation is given for this phenomenon. We propose that the asymmetrical coupling between the light-signal-sensitive part (the ventralateral part, abbreviation: VL) and the light-signal-insensitive part (the dorsalmedial part, abbreviation: DM) of the SCN plays a role in the synchronization degree, which is reflected by the ratio of the number of the directed links from the VL neurons to the DM neurons to the total links of both directions between the VL and the DM. The ratio is assumed to characterize the directed network structure under different photoperiods, which is larger under the SP and smaller under the LP. We found that with the larger ratio in the situation of the SP, the synchronization degree is higher. Our finding may shed new light on the asymmetrical coupling between the VL and the DM, and the network structure of the SCN.

The circadian rhythm induced by the heterogeneous network structure of the suprachiasmatic nucleus

In mammals, the master clock is located in the suprachiasmatic nucleus (SCN), which is composed of about 20 000 nonidentical neuronal oscillators expressing different intrinsic periods. These neurons are coupled through neurotransmitters to form a network consisting of two subgroups, i.e., a ventrolateral (VL) subgroup and a dorsomedial (DM) subgroup. The VL contains about 25% SCN neurons that receive photic input from the retina, and the DM comprises the remaining 75% SCN neurons which are coupled to the VL. The synapses from the VL to the DM are evidently denser than that from the DM to the VL, in which the VL dominates the DM. Therefore, the SCN is a heterogeneous network where the neurons of the VL are linked with a large number of SCN neurons. In the present study, we mimicked the SCN network based on Goodwin model considering four types of networks including an all-to-all network, a Newman-Watts (NW) small world network, an Erdös-Rényi (ER) random network, and a Barabási-Albert (BA) scale free network. We found that the circadian rhythm was induced in the BA, ER, and NW networks, while the circadian rhythm was absent in the all-to-all network with weak cellular coupling, where the amplitude of the circadian rhythm is largest in the BA network which is most heterogeneous in the network structure. Our finding provides an alternative explanation for the induction or enhancement of circadian rhythm by the heterogeneity of the network structure.

Functional network inference of the suprachiasmatic nucleus

In the mammalian suprachiasmatic nucleus (SCN), noisy cellular oscillators communicate within a neuronal network to generate precise system-wide circadian rhythms. Although the intracellular genetic oscillator and intercellular biochemical coupling mechanisms have been examined previously, the network topology driving synchronization of the SCN has not been elucidated. This network has been particularly challenging to probe, due to its oscillatory components and slow coupling timescale. In this work, we investigated the SCN network at a single-cell resolution through a chemically induced desynchronization. We then inferred functional connections in the SCN by applying the maximal information coefficient statistic to bioluminescence reporter data from individual neurons while they resynchronized their circadian cycling. Our results demonstrate that the functional network of circadian cells associated with resynchronization has small-world characteristics, with a node degree distribution that is exponential. We show that hubs of this small-world network are preferentially located in the central SCN, with sparsely connected shells surrounding these cores. Finally, we used two computational models of circadian neurons to validate our predictions of network structure.

Heterogeneity induces rhythms of weakly coupled circadian neurons

The main clock located in the suprachiasmatic nucleus (SCN) regulates circadian rhythms in mammals. The SCN is composed of approximately twenty thousand heterogeneous self-oscillating neurons, that have intrinsic periods varying from 22 h to 28 h. They are coupled through neurotransmitters and neuropeptides to form a network and output a uniform periodic rhythm. Previous studies found that the heterogeneity of the neurons leads to attenuation of the circadian rhythm with strong cellular coupling. In the present study, we investigate the heterogeneity of the neurons and of the network in the condition of constant darkness. Interestingly, we found that the heterogeneity of weakly coupled neurons enables them to oscillate and strengthen the circadian rhythm. In addition, we found that the period of the SCN network increases with the increase of the degree of heterogeneity. As the network heterogeneity does not change the dynamics of the rhythm, our study shows that the heterogeneity of the neurons is vitally important for rhythm generation in weakly coupled systems, such as the SCN, and it provides a new method to strengthen the circadian rhythm, as well as an alternative explanation for differences in free running periods between species in the absence of the daily cycle.

Effects of different per translational kinetics on the dynamics of a core circadian clock model

Living beings display self-sustained daily rhythms in multiple biological processes, which persist in the absence of external cues since they are generated by endogenous circadian clocks. The period (per) gene is a central player within the core molecular mechanism for keeping circadian time in most animals. Recently, the modulation PER translation has been reported, both in mammals and flies, suggesting that translational regulation of clock components is important for the proper clock gene expression and molecular clock performance. Because translational regulation ultimately implies changes in the kinetics of translation and, therefore, in the circadian clock dynamics, we sought to study how and to what extent the molecular clock dynamics is affected by the kinetics of PER translation. With this objective, we used a minimal mathematical model of the molecular circadian clock to qualitatively characterize the dynamical changes derived from kinetically different PER translational mechanisms. We found that the emergence of self-sustained oscillations with characteristic period, amplitude, and phase lag (time delays) between per mRNA and protein expression depends on the kinetic parameters related to PER translation. Interestingly, under certain conditions, a PER translation mechanism with saturable kinetics introduces longer time delays than a mechanism ruled by a first-order kinetics. In addition, the kinetic laws of PER translation significantly changed the sensitivity of our model to parameters related to the synthesis and degradation of per mRNA and PER degradation. Lastly, we found a set of parameters, with realistic values, for which our model reproduces some experimental results reported recently for Drosophila melanogaster and we present some predictions derived from our analysis.

Brain clock driven by neuropeptides and second messengers

The master circadian pacemaker in mammals is localized in a small portion of the brain called the suprachiasmatic nucleus (SCN). It is unclear how the SCN produces circadian rhythms. A common interpretation is that the SCN produces oscillations through the coupling of genetic oscillators in the neurons. The coupling is effected by a network of neuropeptides and second messengers. This network is crucial for the correct function of the SCN. However, models that study a possible oscillatory behavior of the network itself have received little attention. Here we propose and analyze a model to examine this oscillatory potential. We show that an intercellular oscillator emerges in the SCN as a result of the neuropeptide and second messenger dynamics. We find that this intercellular clock can produce circadian rhythms by itself with and without genetic clocks. We also found that the model is robust to perturbation of parameters and can be entrained by light-dark cycles.

Molecular mechanisms that regulate the coupled period of the mammalian circadian clock

In mammals, most cells in the brain and peripheral tissues generate circadian (∼24 h) rhythms autonomously. These self-sustained rhythms are coordinated and entrained by a master circadian clock in the suprachiasmatic nucleus (SCN). Within the SCN, the individual rhythms of each neuron are synchronized through intercellular signaling. One important feature of SCN is that the synchronized period is close to the population mean of cells' intrinsic periods. In this way, the synchronized period of the SCN stays close to the periods of cells in peripheral tissues. This is important because the SCN must entrain cells throughout the body. However, the mechanism that drives the period of the coupled SCN cells to the population mean is not known. We use mathematical modeling and analysis to show that the mechanism of transcription repression in the intracellular feedback loop plays a pivotal role in regulating the coupled period. Specifically, we use phase response curve analysis to show that the coupled period within the SCN stays near the population mean if transcriptional repression occurs via protein sequestration. In contrast, the coupled period is far from the mean if repression occurs through highly nonlinear Hill-type regulation (e.g., oligomer- or phosphorylation-based repression), as widely assumed in previous mathematical models. Furthermore, we find that the timescale of intercellular coupling needs to be fast compared to that of intracellular feedback to maintain the mean period. These findings reveal the important relationship between the intracellular transcriptional feedback loop and intercellular coupling. This relationship explains why transcriptional repression appears to occur via protein sequestration in multicellular organisms, mammals, and Drosophila, in contrast with the phosphorylation-based repression in unicellular organisms and syncytia. That is, transition to protein sequestration is essential for synchronizing multiple cells with a period close to the population mean (∼24 h).

Timing of neuropeptide coupling determines synchrony and entrainment in the mammalian circadian clock

Robust synchronization is a critical feature of several systems including the mammalian circadian clock. The master circadian clock in mammals consists of about 20000 ‘sloppy’ neuronal oscillators within the hypothalamus that keep robust time by synchronization driven by inter-neuronal coupling. The complete understanding of this synchronization in the mammalian circadian clock and the mechanisms underlying it remain an open question. Experiments and computational studies have shown that coupling individual oscillators can achieve robust synchrony, despite heterogeneity and different network topologies. But, much less is known regarding the mechanisms and circuits involved in achieving this coupling, due to both system complexity and experimental limitations. Here, we computationally study the coupling mediated by the primary coupling neuropeptide, vasoactive intestinal peptide (VIP) and its canonical receptor, VPAC2R, using the transcriptional elements and generic mode of VIP-VPAC2R signaling. We find that synchrony is only possible if VIP (an inducer of Per expression) is released in-phase with activators of Per expression. Moreover, anti-phasic VIP release suppresses coherent rhythms by moving the network into a desynchronous state. Importantly, experimentally observed rhythms in VPAC2R have little effect on network synchronization, but can improve the amplitude of the SCN network rhythms while narrowing the network entrainment range. We further show that these findings are valid across several computational network models. Thus, we identified a general design principle to achieve robust synchronization: An activating coupling agent, such as VIP, must act in-phase with the activity of core-clock promoters. More generally, the phase of coupling is as critical as the strength of coupling from the viewpoint of synchrony and entrainment.

Synchronization among multiple oscillators is a common theme in many biological and engineered systems. Here, we look at its use by the mammalian biological clock to keep accurate time. Through biochemical interactions among a network of inaccurate neuron clocks, a strong precise clock is produced. Although we are gradually learning more about these biochemical interactions, the details still remain largely unclear. Studies, both computational and experimental, have shown that the strength of the rhythmic interaction critically decides if a system can synchronize, i.e., the interactions must be strong enough. In this work, we show that the rhythmic interaction between these neuronal clocks must be timed correctly (in the right phase) in addition to being strong enough to synchronize the network. Activating (repressing) interactions must coincide with other activators (repressors) in each neuron to achieve synchrony. Since this principle imposes certain properties on synchronizing interactions, these properties can be used to identify and understand novel interaction mechanisms. Further, these principles are applicable to interactions between cellular oscillators in other tissues and organisms.

GABA networks destabilize genetic oscillations in the circadian pacemaker

Systems of coupled oscillators abound in nature. How they establish stable phase relationships under diverse conditions is fundamentally important. The mammalian suprachiasmatic nucleus (SCN) is a self-sustained, synchronized network of circadian oscillators that coordinates daily rhythms in physiology and behavior. To elucidate the underlying topology and signaling mechanisms that modulate circadian synchrony, we discriminated the firing of hundreds of SCN neurons continuously over days. Using an analysis method to identify functional interactions between neurons based on changes in their firing, we characterized a GABAergic network comprised of fast, excitatory, and inhibitory connections that is both stable over days and changes in strength with time of day. By monitoring PERIOD2 protein expression, we provide the first evidence that these millisecond-level interactions actively oppose circadian synchrony and inject jitter into daily rhythms. These results provide a mechanism by which circadian oscillators can tune their phase relationships under different environmental conditions.

Entrainment range of nonidentical circadian oscillators by a light-dark cycle

The suprachiasmatic nucleus (SCN) is a principal circadian clock in mammals, which controls physiological and behavioral daily rhythms. The SCN has two main features: Maintaining a rhythmic cycle of approximately 24 h in the absence of a light-dark cycle (free-running period) and the ability to entrain to external light-dark cycles. Both free-running period and range of entrainment vary from one species to another. To understand this phenomenon, we investigated the diversity of a free-running period by the distribution of coupling strengths in our previous work [Phys. Rev. E 80, 030904(R) (2009)]. In this paper we numerically found that the dispersion of intrinsic periods among SCN neurons influence the entrainment range of the SCN, but has little influence on the free-running periods under constant darkness. This indicates that the dispersion of coupling strengths determines the diversity in free-running periods, while the dispersion of intrinsic periods determines the diversity in the entrainment range. A theoretical analysis based on two coupled neurons is presented to explain the results of numerical simulations.

Weakly circadian cells improve resynchrony

The mammalian suprachiasmatic nuclei (SCN) contain thousands of neurons capable of generating near 24-h rhythms. When isolated from their network, SCN neurons exhibit a range of oscillatory phenotypes: sustained or damping oscillations, or arrhythmic patterns. The implications of this variability are unknown. Experimentally, we found that cells within SCN explants recover from pharmacologically-induced desynchrony by re-establishing rhythmicity and synchrony in waves, independent of their intrinsic circadian period We therefore hypothesized that a cell's location within the network may also critically determine its resynchronization. To test this, we employed a deterministic, mechanistic model of circadian oscillators where we could independently control cell-intrinsic and network-connectivity parameters. We found that small changes in key parameters produced the full range of oscillatory phenotypes seen in biological cells, including similar distributions of period, amplitude and ability to cycle. The model also predicted that weaker oscillators could adjust their phase more readily than stronger oscillators. Using these model cells we explored potential biological consequences of their number and placement within the network. We found that the population synchronized to a higher degree when weak oscillators were at highly connected nodes within the network. A mathematically independent phase-amplitude model reproduced these findings. Thus, small differences in cell-intrinsic parameters contribute to large changes in the oscillatory ability of a cell, but the location of weak oscillators within the network also critically shapes the degree of synchronization for the population.

Circadian rhythms are daily, near 24-h oscillations in biological processes that nearly all organisms on Earth experience. Single cells contain a molecular clock that drives circadian rhythms in physiology and, when many cells synchronize in a population, daily behaviors. We hypothesized that small differences in intrinsic cellular properties allow for a diversity of circadian periods and amplitudes across cells. We observed circadian cells and their synchrony before, during, and after limiting communication between cells and then compared their intrinsic properties to their resynchronization behavior. We found that arrhythmic, weakly oscillating, and self-sustained circadian cells rejoined the rhythmic population independent of their cell-intrinsic oscillations. Using a mechanistic computational model of circadian cells, we found that resynchronization could be enhanced by including more weak oscillators or by placing weak oscillators at more connected nodes in the network. We conclude that intrinsic properties (e.g. oscillator weakness and responsiveness) and network structure (e.g. positions of weak oscillators) can independently buffer tissue rhythms from perturbations. This reveals how cellular and network properties impose rules on systems of circadian cells that must achieve synchrony from a desynchronized state, for example during perinatal development or when forced to overcome societal constraints on sleep-wake behavior, such as working early or late shifts.

Modeling the seasonal adaptation of circadian clocks by changes in the network structure of the suprachiasmatic nucleus

The dynamics of circadian rhythms needs to be adapted to day length changes between summer and winter. It has been observed experimentally, however, that the dynamics of individual neurons of the suprachiasmatic nucleus (SCN) does not change as the seasons change. Rather, the seasonal adaptation of the circadian clock is hypothesized to be a consequence of changes in the intercellular dynamics, which leads to a phase distribution of electrical activity of SCN neurons that is narrower in winter and broader during summer. Yet to understand this complex intercellular dynamics, a more thorough understanding of the impact of the network structure formed by the SCN neurons is needed. To that effect, we propose a mathematical model for the dynamics of the SCN neuronal architecture in which the structure of the network plays a pivotal role. Using our model we show that the fraction of long-range cell-to-cell connections and the seasonal changes in the daily rhythms may be tightly related. In particular, simulations of the proposed mathematical model indicate that the fraction of long-range connections between the cells adjusts the phase distribution and consequently the length of the behavioral activity as follows: dense long-range connections during winter lead to a narrow activity phase, while rare long-range connections during summer lead to a broad activity phase. Our model is also able to account for the experimental observations indicating a larger light-induced phase-shift of the circadian clock during winter, which we show to be a consequence of higher synchronization between neurons. Our model thus provides evidence that the variations in the seasonal dynamics of circadian clocks can in part also be understood and regulated by the plasticity of the SCN network structure.

Circadian clocks drive the temporal coordination of internal biological processes, which in turn determine daily rhythms in physiology and behavior in the most diverse organisms. In mammals, the 24-hour timing clock resides in the suprachiasmatic nucleus (SCN) of the hypothalamus. The SCN is a network of interconnected neurons that serves as a robust self-sustained circadian pacemaker. The electrical activity of these neurons and their synchronization with the 24-hour cycle is established via the environmental day and night cycles. Apart from daily luminance changes, mammals are exposed to seasonal day length changes as well. Remarkably, it has been shown experimentally that the seasonal adaptations to different photoperiods are related to the modifications of the neuronal activity of the SCN due to the plasticity of the network. In our paper, by developing a mathematical model of the SCN architecture, we explore in depth the role of the structure of this important neuronal network. We show that the redistribution of the neuronal activity during winter and summer can in part be explained by structural changes of the network. Interestingly, the alterations of the electrical activity patterns can be related with small-world properties of our proposed SCN network.