Sage Journals: Discover world-class research

Abstract

Synchronizing the circadian timing system (CTS) to external light/dark cycles is crucial for homeostasis maintenance and environmental adaptation. The CTS is organized hierarchically, with the central pacemaker located in the suprachiasmatic nuclei (SCN) generating coherent oscillations that are entrained to light/dark cycles. These oscillations regulate the release of glucocorticoids by the hypothalamus-pituitary-adrenal (HPA) axis, which acts as a systemic entrainer of peripheral clocks throughout the body. The SCN adjusts its network plasticity in response to variations in photoperiod, leading to changes in the rhythmic release of glucocorticoids and ultimately impacting peripheral clocks. However, the effects of photoperiod-induced variations of glucocorticoids on the synchronization of peripheral clocks are not fully understood, and the interaction between jetlag adaption and photoperiod changes is unclear. This study presents a semi-mechanistic mathematical model to investigate how the CTS responds to changes in photoperiod. Specifically, the study focuses on the entrainment properties of a system composed of the SCN, HPA axis, and peripheral clocks. The results show that high-amplitude glucocorticoid rhythms lead to a more coherent phase distribution in the periphery. In addition, our study investigates the effect of photoperiod exposure on jetlag recovery time and phase shift, proposing different interventional strategies for eastward and westward jetlag. The findings suggest that decreasing photic exposure before jetlag during eastward traveling and after jetlag during westward traveling can accelerate jetlag readaptation. The study provides insights into the mechanisms of CTS organization and potential recovery strategies for transitions between time zones and lighting zones.

Keywords

circadian rhythm SCN HPA axis seasonality jetlag

Introduction

Adapting to predictable environmental changes is regulated by the circadian timing system (CTS), a multi-scale, hierarchical, and distributed mechanism. At the single-cell level, molecular clocks are driven by positive and negative feedback loops consisting of a group of “clock” genes and proteins that generate self-sustained, autonomous oscillations with a nearly 24-h period. The endogenous molecular clock machinery is adaptable and can be synchronized to external periodic stimuli (zeitgebers). Among different entraining factors (temperature, feeding, exercise, etc.), daily light-dark cycles serve as the primary synchronizer of the circadian rhythm.

In mammals, light is sensed by photoreceptors in the retina, and photic information is transduced to the hypothalamic suprachiasmatic nuclei (SCN) via the intrinsically photoreceptive retinohypothalamic tract (S.-K. Chen et al., 2011). The SCN is a bilateral structure that comprises approximately 20,000 neurons; each neuron contains an autonomous circadian oscillator. Those neurons are heterogeneous. They couple with each other and function as a network, generating coherent rhythmicity through spatial and temporal coordination (Herzog, 2007). As a result, the topographically organized neurons exhibit spontaneous cyclic fluctuations in circadian clock gene expression. This rhythm is endogenous, meaning that the oscillation will be preserved even in the absence of light/dark entrainment (Webb et al., 2009). By receiving direct photic input and employing organized coupling mechanisms to reinforce its intrinsic rhythmicity, the SCN plays a role as the central pacemaker in the CTS.

Circadian regulation is astonishingly widespread: nearly half of all mammalian genes are rhythmically expressed in one or more tissues (Yan et al., 2008; Acevedo et al., 2021). Besides the SCN neurons, clock genes and proteins have been observed in cells and tissues throughout the body (such as the lung, liver, skeletal muscle, and heart), indicating the existence of endogenous clocks among the peripheral cells (Welsh et al., 2004; Yoo et al., 2004; Zhang et al., 2014). The peripheral clocks exert broad-ranging control over many physiological processes (Nahmias and Androulakis, 2021; Schlesinger et al., 2022), including steroidogenesis (H. Chen et al., 2017), glucose homeostasis (Gachon et al., 2017), and immune function (Beam et al., 2020). The peripheral cellular oscillators are heterogeneous in their intrinsic periods and, thus, must be synchronized to establish coherent tissue rhythms. While SCN exhibits inherent coupling and can maintain synchronization ex vivo for extended periods, peripheral cells display ensemble oscillations that rapidly dampen without intercellular coupling signals (Abraham et al., 2010), indicating that systemic entrainment is necessary for peripheral clock synchronization.

The light-sensitive SCN conveys photic signals through its rhythmic output, entraining the neuroendocrine system, which eventually synchronizes the peripheral clocks (Kalsbeek et al., 2010). The hypothalamic-pituitary-adrenal (HPA) axis plays an important role in coupling the oscillation between the SCN and the periphery. In the paraventricular nucleus (PVN), corticotropin-releasing hormone (CRH) is rhythmically secreted by CRH neurons, regulated by the SCN-produced neurotransmitter AVP (Kalsbeek et al., 2012). The CRH travels to the anterior pituitary, stimulating the production of adrenocorticotropic hormone (ACTH), which activates the release of glucocorticoids (CORT) from the adrenal cortex. The daily rhythmicity of plasma glucocorticoid levels functions as a systemic synchronizer that coordinates the phase of peripheral clocks by binding to the glucocorticoid-response elements (GREs) regions, activating the transcription of a number of clock genes (Yamamoto et al., 2005; Reddy et al., 2007; So et al., 2009; Cuesta et al., 2015). Glucocorticoids also exert negative feedback on both PVN and the anterior pituitary by inhibiting the secretion of CRH and ACTH, therefore indirectly regulating its own secretion (Stanojević et al., 2018; Focke and Iremonger, 2020). Entrained to the SCN, glucocorticoids display a daily peak at the habitual sleep-wake transition, with levels peaking at the onset of the active phase and declining during the inactive phase (defined as day for diurnal species and night for nocturnal species) (Oster et al., 2017).

Apart from the light/dark cycle, the SCN neurons also adapt to photoperiod changes by adjusting their ensemble phase plasticity. During short photoperiod months (winter), the SCN neurons exhibit a narrower phase distribution compared to long photoperiod months (summer) (VanderLeest et al., 2007; Meijer et al., 2010). Governed by the substantial differences in the SCN phase distribution, glucocorticoid rhythms exhibit distinct variations under photoperiod changes, with amplitudes highest during short photoperiod (winter) and lowest during long photoperiods (summer) (Amirat et al., 1980; Otsuka et al., 2012). Studies have shown that the robustness of glucocorticoid rhythmicity is strongly associated with peripheral circadian alignment (Morris et al., 2016). Therefore, it has been speculated that photoperiod changes induced by SCN and HPA rhythm will differentially regulate the synchronization of peripheral oscillators while detailed mechanisms remain to be elucidated (Pierre et al., 2018).

Optimal physiological functions occur when the central and peripheral rhythms are aligned with the external light schedules. Conversely, the misalignment between the central clock, peripheral clock, and the environment may increase the risk of disease, with the immune system, metabolism, and cardiovascular system negatively impacted (Morris et al., 2016; Olejniczak et al., 2022). Considering the complex hierarchical structure of the CTS, mathematical modeling offers advantages in accessing the physiological system at a holistic level, uncovering potential mechanisms, and quantifying properties that are experimentally challenging to identify (Androulakis, 2021). While previous theoretical studies have investigated the effects of jetlag and seasonality using relatively simplified mathematical models, and have identified the impact of intrinsic period and amplitude on jetlag recovery (VanderLeest et al., 2009; Schmal et al., 2015; Ananthasubramaniam et al., 2020; Tokuda et al., 2020), a comprehensive model that can be used to study jetlag and seasonal-related mechanisms in detail is still lacking. Furthermore, the interaction between time zone transition and lighting zone transition has been overlooked in previous research.

The current study presents a semi-mechanistic mathematical model that represents the hierarchical structure of the CTS and examines its systemic response to light-dark entrainment, annual photoperiod changes, and photoperiod-dependent jetlag responses. Focusing on peripheral clock synchronization properties, our model identified differential peripheral phase distribution under different photoperiods. By investigating the systemic circadian alterations under lighting zone and time zone transitions, our model shows that reducing photic exposure before a phase-advanced jetlag can accelerate the transition rate of the system while reducing photic exposure after a phase-delayed jetlag can also speed up the recovery rate. Furthermore, our study revealed a positive correlation between the magnitude of phase shift during jetlag and the duration of jetlag recovery, suggesting that preserving the intrinsic phase shift of the system may lead to more effective jetlag recovery. These findings shed light on the mechanisms underlying photoperiod adaptation and non-pharmacological interventions for jetlag, providing a basis for further studies and potential clinical applications.

Materials and Methods

The “Model Construction” section of the Supplementary Material (Suppl.) provides an in-depth description of the developed mathematical model that characterizes the seasonal entrainment of the HPA axis induced by the SCN organization. In brief, the hierarchical organization of the CTS is modeled at both the central and peripheral levels. As shown in Figure 1, in the central compartment, the SCN is entrained to the light/dark cycles and can adjust its network plasticity according to photoperiod changes. The SCN synchronizes the autonomous oscillations in the HPA axis, which, in turn, regulates cells in the periphery. SCN neurons can be neuroanatomically divided into two subpopulations: the ventrolateral “core” (VL) and the dorsomedial “shell” (DM). The core neurons secret vasoactive intestinal peptide (VIP), neurotensin (NT), and gastrin-releasing peptide (GRP) and receive afferent signals from the retina. The shell neurons express arginine vasopressin (AVP), angiotensin II (AII), and met-enkephalin (mENK). Neuropeptides VIP and AVP are the main synchronizers of the core and shell neuron population (Ono et al., 2021). Based on the experimentally identified structures, the SCN model in the present work consists of two ensembles of neuronal oscillators, which represent the core and the shell, respectively. At the tissue level, the neurons communicate with each other via neurotransmitters VIP and AVP (Moore et al., 2002). Both neurotransmitters are released upon the activity of PER/CRY proteins and can activate the transcription of Per/Cry mRNA (Bernard et al., 2007; Li and Androulakis, 2021).

Figure 1.

The modeled circadian timing system (CTS) is structured in three levels: the suprachiasmatic nucleus (SCN), the hypothalamic-pituitary-adrenal (HPA) axis, and the periphery. The SCN comprises of a light-sensitive core region and a light-insensitive shell region, with each neuron subpopulation consisting of mutually coupled cellular oscillators. The SCN’s coupling plasticity can change under different photoperiods. The HPA axis exhibits a negative feedback loop consisting of corticotropin-releasing hormone (CRH), adrenocorticotropic hormone (ACTH), and corticosterone (CORT). The HPA axis is entrained by the SCN via the negative regulatory role of arginine vasopressin (AVP) on CRH, as the model is based on a nocturnal species. In the peripheral compartment, cells are composed of independent oscillators that are synchronized through glucocorticoid regulation. Abbreviations: VL = ventrolateral; DM = dorsomedial; LP = long photoperiod; SP = short photoperiod; mRNA = messenger ribonucleic acid.

The SCN synchronizes the oscillations of the HPA axis through the SCN-mediated secretion of AVP in the paraventricular hypothalamic nucleus (Kalsbeek et al., 2010, 2012). To simulate the endogenous HPA oscillatory property, we used a previously developed HPA axis model (Mavroudis et al., 2012, 2014; Rao and Androulakis, 2017, 2019; Androulakis, 2021) and incorporated the synchronization signals from the SCN (Li and Androulakis, 2021, 2022). After diffusing into the peripheral cytoplasm, the glucocorticoids exert their effects by binding to glucocorticoid receptors in the periphery. The hyperphosphorylation of the GR upon its activation is considered the rate-determining step. A glucocorticoid-mediated indirect activation mechanism was used to model the phosphorylation kinetics. Binding to glucocorticoid ligand, the glucocorticoid receptor complex (FGR) translocates to the nucleus where they interact with glucocorticoid-responsive elements (GREs) promoter regions to activate glucocorticoid-responsive clock genes.

The periphery consists of a heterogeneous ensemble of cells whose free-running periods and phases vary significantly (Honma et al., 2004). To account for the heterogeneity of the periphery, we simulated cell-cell variability using Sobol sampling ( $\pm 3 %$ ) of the parameters of peripheral circadian dynamics, which reflect the stochastic expression of intracellular activities. A population of 1000 peripheral cells was generated for all the simulation experiments performed in the present study. The detailed model of the peripheral clock genes contains a positive and negative feedback loop, which leads to an autonomous circadian oscillation (Geier et al., 2005). Studies have shown that GCs acutely induce Per2 gene expression in rodent peripheral clocks and human peripheral blood mononuclear cells (PBMCs) (Cuesta et al., 2015). In accordance, the transcriptional oscillation of clock genes in our model is regulated by the rhythmicity of glucocorticoids. The entrainment effects of GCs on PCG activity are accounted for by including the activation term of the nucleus glucocorticoid receptor complex. The coupling strength between the HPA axis and peripheral clocks is quantified by parameter $k_{v}$ . The coupling parameter was estimated to calibrate the peripheral phase of Per/Cry under L12/D12 with experimental data in rodents’ liver (the peaking time of Per/Cry peaks 6 h after the offset of light period) (Korenčič et al., 2012).

In the periphery, synchronizing heterogeneous oscillations through the entrainment of the systemic signal is essential for generating a consistent physiological rhythm (Aton and Herzog, 2005). To study the impact of systemic photoperiod-related seasonal plasticity in the synchronization state of the peripheral compartment, we evaluate the degree of synchronization among the peripheral clocks. The temporal profile synchronization $R_{s y n}$ was introduced for quantifying the coherency of Per/Cry mRNA concentration profiles, represented by the ratio between the mean-field variance and the mean-variance of each oscillator (Suppl. Eqs. 45-47) (Garcia-Ojalvo et al., 2004; Pierre et al., 2018).

Results

Characterizing the Hierarchical Structure of the CTS and the Entrainment Properties of the Peripheral Clocks

To explore the stochastic nature of the periphery in the absence of the systemic glucocorticoid entraining signal, we calculated the phase and period distributions of peripheral cells when the coupling coefficient $k_{v}$ (Suppl. Eq. 38) is zero. As shown in Figure 2b, the peripheral cell phases exhibit a uniform distribution from 0 to 2π. The periods of individuals cells, on the other hand, are normally distributed with an average period of 24.70 h (Figure 2a). Therefore, in the absence of the glucocorticoid entraining signal, the populations of peripheral cells exhibit damped oscillations at the ensemble level. Moreover, as we restore the value of $k_{v}$ , a quick synchronization among the peripheral cells was observed, indicating the critical role of glucocorticoids in entraining the peripheral cells (Suppl. Fig. S1).

Figure 2.

The systemic glucocorticoid signal transduced by the SCN and the HPA axis synchronizes the peripheral clocks. (a) The phase distribution and period distribution of peripheral cells in the absence of glucocorticoid entrainment. (b) The peripheral cells exhibit random phase and period distribution when not entrained by glucocorticoids. Abbreviations: SCN = suprachiasmatic nucleus; HPA = hypothalamic-pituitary-adrenal.

Regulated by the zeitgeber, the entrainment status of an oscillator is highly dependent on (a) period mismatch $τ - T$ ( $τ$ is the intrinsic period of the oscillator, whereas $T$ is the period of the zeitgeber) and (b) the intensity of the zeitgeber (Bordyugov et al., 2015). The entrainment range is defined as period mismatches ( $τ - T$ ) under which an oscillator can still be entrained and exhibit the same period with the zeitgeber. To explore the entrainment properties of the SCN—HPA—periphery system, we examined the “Arnold tongue” by identifying regions of entrainment on a zeitgeber period—zeitgeber intensity plane for compartments including the SCN, the HPA, and the periphery, as shown in Figure 3. Zeitgeber intensity is quantified by the parameter $l i s$ , which represents the light intensity received by the SCN neurons (Suppl. Eq. 1). To account for the fact that the SCN and the periphery in the model are composed of multiple cells rather than a single cell, we determined the synchronization status of each compartment by comparing their ensemble average periods with the zeitgeber period and computed the synchronization level $R_{s y n}$ for the entrained state. This approach extended the scope of the Arnold tongue characterization from a single oscillator to multiple oscillators. The level of synchronization in the SCN and the periphery was assessed based on the circadian dynamics of Per/Cry mRNA expression across all the compartments. In Figure 3, the synchronization levels of the SCN and peripheral ensembles are depicted by color maps, indicating the degree of synchronization when the period of the ensemble matches the period of the zeitgeber.

Figure 3.

Entrainment properties of the system were identified by the “Arnold tongue” method. The entrainment range was identified in a zeitgeber period-zeitgeber intensity plane for (a) the SCN, (b) the HPA axis, (c) the peripheral clocks. (d) three Arnold tongues derived from different compartments were plotted together for comparison. Gray dots denote the SCN, blue dots denote the HPA, and black dots denote the periphery. In the SCN, the Per/Cry mRNA concentration in the light-sensitive core was used as a representative component to determine the entrainment status. In the periphery, the Per/Cry mRNA concentration was used as a representative component to determine the entrainment status. The HPA axis’ entrainment status was represented by its output glucocorticoids. An oscillation is considered to be entrained when its period equals the zeitgeber period. The color in (a) and (c) denotes the synchronization level among the SCN and peripheral ensembles, respectively. The entrainment ranges diminish as the photic signal is transduced from upstream to downstream in the circadian timing system. Abbreviations: SCN = suprachiasmatic nucleus; HPA = hypothalamic-pituitary-adrenal; mRNA = messenger ribonucleic acid.

The measured entrainment range for the SCN compartment increases with increasing zeitgeber intensity, leading to a triangular region of entrainment, with higher synchronization levels concentrated in areas where the zeitgeber period is close to $24 h$ and when zeitgeber intensities are high (Figure 3a). In contrast, downstream HPA and peripheral compartments exhibit significantly diminished entrainment ranges, being entrained only when the zeitgeber period is close to their intrinsic period (Figure 3b and 3c). The comparison of entrainment ranges for the three compartments in Figure 3d reveals that the SCN exhibits a wider range of entrainment compared to the peripheral compartment. This is due to the direct entrainment of the SCN by light, whereas the periphery relies on the output of the HPA axis for entrainment. This design enables the CTS to function as a gatekeeper mechanism, ensuring that the periphery does not synchronize with zeitgebers that have periods that are either too large or too small in vivo.

Photoperiod-Induced Variations in HPA and Their Impact on Peripheral Clock Synchronization

The photoperiod-induced variations in the HPA axis circadian dynamics are represented in Figure 4. Increasing photoperiod from $8 h$ (PP 8) to $16 h$ (PP 16) resulted in a decrease in glucocorticoid amplitude. Figure 4a shows that a decreased photoperiod induces an increase in glucocorticoid activity. In Figure 4b, the simulation results predicting corticosterone amplitude over a year were compared with experimental data obtained from rodent studies (Amirat et al., 1980). Our model prediction displayed significant annual fluctuations in corticosterone amplitude, characterized by a decrease during spring, a nadir in early summer, and a resurgence in autumn, consistent with experimental observations. To compare phase and amplitude variations of glucocorticoids during different photoperiods, we plotted light onset and offset schedules with the length of the light period changing while the middle of the light period remains constant in Figure 4c. A compass plot was used to represent the phase and amplitude of the glucocorticoid profiles under the simulated light schedules, as shown in Figure 4d. Our model predicts that increasing photoperiod leads to a phase advance in the timing of peak glucocorticoid levels, which is consistent with experimental findings (Laakso et al., 1994; Vondrasová et al., 1997).

Figure 4.

Photoperiod-induced variations in the HPA dynamics were represented by the phase and amplitude changes in glucocorticoids. (a) Circadian profiles of the glucocorticoid under different photoperiods, all light schedule starts when zeitgeber time = 0 h (ZT 0). The colorful patterns denote the light onset and offset times, and the shadow denotes darkness. (b) Annual dynamics of the HPA axis compared to experimental data. The black solid line denotes the predicted corticosterone amplitude, while the black dashed line denotes the experimentally measured corticosterone amplitude. (c) The introduced photoperiod schedule at a real-time scale. The colorful patterns denote the light onset and offset times, and the shadow denotes darkness. (d) Phase and amplitude relations of glucocorticoids under different photoperiod schedules. As the photoperiod increases, the amplitude of glucocorticoids decreases. The phases of the circadian rhythms were advanced for the long photoperiod and delayed for the short photoperiod. Abbreviations: CORT = corticosterone; PP = photoperiod; HPA = hypothalamic-pituitary-adrenal.

Entrained by glucocorticoids, the peripheral clocks exhibit photoperiod-dependent variation in their synchronization status and phase distribution. To investigate this relationship, we modeled a winter, spring, summer, and autumn light schedule with photoperiods of $8 h$ , $12 h$ , $16 h$ , and $10 h$ , respectively, as shown in Figure 5a. Using Per/Cry mRNA as the representative peripheral clock component, we found that the seasonal distributions of peripheral clocks exhibited photoperiod-dependent variations. Our results showed that the peripheral clocks displayed the broadest phase distribution during the summer months, consistent with the observed behavior patterns of organisms exposed to long photoperiods (Saper et al., 2005). The compass plot in Figure 5b revealed a decrease in the peripheral ensemble amplitude as the photoperiod increased, and a $~ 2 h$ phase delay during short photoperiod winter compared to long photoperiod summer.

Figure 5.

The impact of photoperiod-induced variations in the hypothalamic-pituitary-adrenal (HPA) axis on peripheral clock synchronization. (a) The phase distribution of peripheral Per/Cry mRNA under different photoperiods. The widest phase distribution occurs during the summer, in agreement with the behavior patterns of organisms in long-photoperiod months. (b) Dynamics of the peripheral Per/Cry mRNA under different photoperiod schedules. Regulated by the photoperiod-responsive glucocorticoid rhythms, the peripheral cells exhibit a higher level of amplitude during short photoperiods (SP) and a lower amplitude during long photoperiods (LP). Together, the figure illustrates that the photoperiod-induced variations in the HPA axis play a crucial role in regulating the peripheral clock. Abbreviations: PP = photoperiod; mRNA = messenger ribonucleic acid.

Impact of Photoperiod and Time Zone Changes on Jetlag Readaptation Time

Light schedule disruptions due to long-distance flights in contemporary society can result in notable circadian misalignment among different compartments of the CTS, which is associated with travel fatigue, sleep loss, and decreased physical and mental performance (Bin et al., 2020). However, the impact of photoperiod zone transitions is often overlooked in research on travel fatigue and jetlag (Wieczorek et al., 2013, 2016). We aim to investigate the consequences of both lighting zone and time zone transitions using our model.

We vary the duration of light exposure in lighting zone transition schedules, starting from short photoperiod (SP) of 8 h to long photoperiod (LP) of 16 h and vice versa, to study the effects of lighting zone transition. As shown in Figure 6a, a quick synchronization process occurs when the system moves from LP to SP, along with a phase advance at the ensemble level. Conversely, when the system transitions from SP to LP, a gradual desynchronization takes place with a phase delay, suggesting that prolonged exposure to light can cause a phase delay in the circadian clock (Figure 6b). Similarly, we subject the system to a $6 h$ advance (east-bound) and $6 h$ delay (west-bound) jetlag to study the effect of time zone transition. The transition across time zones results in a gradual shift of the peripheral phase toward the new light schedule, with jetlag inducing misalignment among the peripheral clocks, especially in the phase advance schedule (Figure 6c and Suppl. Fig. S4).

Figure 6.

Peripheral clocks dynamics during light schedule changes. The transition of peripheral Per/Cry mRNA and protein concentrations were analyzed when (a) transitioning from long photoperiod (LP, PP = 16) to short photoperiod (SP, PP = 8) and when (b) transitioning from SP to LP. The phase transition of the peripheral ensemble upon lighting zone shifting was also examined. Results indicate that the system takes longer to transition from SP to LP compared to LP to SP, with a recovery time of 18 days for the transition from LP to SP and 27 days for the transition from SP to LP. Furthermore, a phase delay was observed when transitioning from SP to LP. (c) Desynchronization of the peripheral ensemble upon exposure to jetlag during a 6 h phase-advanced jetlag schedule and (d) a 6 h phase-delayed jetlag schedule. The photoperiod before and after jetlag remains constant at 8 h (PP = 8). The red line in panels represents the ensemble average phase calculated from Per/Cry mRNA, while the black dots represent the phase calculated from individual cells. In addition, the system exhibits a higher level of desynchronization and a longer recovery time during phase-advanced jetlag schedules. Abbreviations: ZT = zeitgeber time; mRNA = messenger ribonucleic acid.

To investigate the adaptability of the system under different photoperiod transitions across time zones, we defined the starting point A and the destination point B, each with its own photoperiods. The difference in light onset between A and B is defined as the time zone crossed. Jetlag is considered advanced if the light onset in B is earlier than that in A and delayed if occurs later. We varied factors such as photoperiods in point A ( $P P_{A}$ ) and B ( $P P_{B}$ ), the direction of jetlag, and time zone crossed during jet lag to create 175 different transition schedules. Detailed information on the conditions for each schedule can be found in the Supplementary Materials. To measure the transition time and phase, we recorded the peripheral ensemble peaking phase (represented by Per/Cry mRNA) when the system reaches steady states before and after jetlag. The system is considered to enter a steady state when the phase of the periphery reaches a constant value with a less than $1 h$ variation between different days. The time difference between the onset of jetlag and the establishment of a new steady state is defined as the “jetlag transition time.,” while the phase difference is the “phase shift.”

In Figure 7, we present an example to investigate the interaction between jetlag and photoperiod transitions. We plotted the jetlag transition time and phase shift under a $6 h$ time delay or advance with different photoperiods before and after time zone transition. Notably, we observed a significant discrepancy in phase shift and jetlag transition time under both delayed and advanced schedules, indicating that the differential photoperiods can efficiently enhance or impede jetlag transition speed. Figure 7a and 7b show an inverted phase-variation pattern between the advanced schedule and the delayed schedule. The augmentation of the photoperiod at point B notably reduces the phase shift during the phase-advanced schedule while concurrently exacerbating the phase shift during the phase-delayed schedule. The observation is consistent with the results presented in Figure 6a and 6b; as the transition to a prolonged photoperiod phase delays the system, it partially counteracts the effects of the externally imposed phase advance time zone transitions. Interestingly, a reduced phase shift appears to facilitate a quicker jetlag transition rate. We observed that under phase delay jetlag, shorter photoperiod schedules at point B accelerated the transition rate, while under phase advance schedules, shorter photoperiods at point A generally hastened the transition time. This same conclusion can be extrapolated from other time zone transition schedules. The corresponding figures illustrating jetlag transition and phase shift variation can be found in the Supplementary Figure S6.

Figure 7

The interplay between lighting zone transition and time zone transition as demonstrated through the measurement of Per/Cry mRNA in the peripheral ensemble. Panel (a) shows the phase shift under a 6-h delayed schedule, panel (b) shows the phase shift under a 6-h advanced schedule, panel (c) shows the jetlag transition time under a 6-h delayed schedule, and panel (d) shows the jetlag transition time under a 6-h advanced schedule. The system is considered fully recovered when the peripheral ensemble peaking phase lies within a 1-h range of the new steady-state phase. The time difference between the onset of jetlag and the establishment of a new steady state is recorded as the jetlag recovery time, and the phase difference between the old steady state and the new steady state phase in the peripheral ensemble is recorded as a phase shift. An increase in phase difference is observed with an increased photoperiod exposure after jetlag during phase delay schedules, and a decrease in phase difference is observed with an increased photoperiod exposure after jetlag during phase advance schedules. Jetlag transition time exhibits a positive correlation with the phase shift. Abbreviations: PP = photoperiod; mRNA = messenger ribonucleic acid.

To quantitatively assess the impact of various light schedule factors on the systemic transition response, we analyzed the peripheral response under various jetlag schedules and conducted multiple linear regression analyses separately for the advanced and delayed cases. Models with the form of $f_{w, b} (P P_{A}, P P_{B}, T Z) = w_{P P A} * P P_{A} + w_{P P B} * P P_{B} + w_{T Z} * T Z + b$ were employed for both phase-advanced and phase-delayed transition schedules. $f_{w, b}$ represents the responses of the linear regression model, including transition time and phase shift. $P P_{A}$ , $P P_{B}$ , $T Z$ represent three independent variables, including photoperiod at point A, photoperiod at point B, and time zone crossed during the jetlag, respectively. As demonstrated by the results obtained from the linear regression analysis in Figure 8, we plotted the correlations between the observed responses and predictors. A correlation was considered significant when the p-value was less than $0.05$ . Detailed information on the coefficient values and p-values of the independent variables in various scenarios are documented in Supplementary Table 3. During advanced phase jetlag, the intrinsic phase shift of the system exhibits a positive correlation with the photoperiod at point A and a negative correlation with the photoperiod at point B, whereas during phase-delayed jetlag, an inverse trend is observed. Such a correlation is consistent with the results presented in Figure 7 and strongly implies that the duration of light exposure before and after jetlag has a marked effect on the intrinsic phase shift of the body and should be taken into account when devising jetlag interventions.

Figure 8.

Linear regression analysis of the impact of photoperiod before and after jetlag, and time zone crossing on phase shift and jetlag transition time. Red dots represent data obtained from phase-advanced jetlag simulations, while blue dots represent data from phase-delayed schedules. Shadows represent the predicted range from the linear regression model. A linear regression model of $f_{w, b} (P P_{A}, P P_{B}, T Z) = w_{P P A} * P P_{A} + w_{P P B} * P P_{B} + w_{T Z} * T Z + b$ was used in both phase-advanced and phase-delayed transition schedules. $f_{w, b}$ denotes the responses of the linear regression model, including transition time and phase shift. Detailed information on the coefficient values and p-values of the independent variables in various scenarios are documented in Supplementary Table 3. A correlation was considered significant when the p-value was less than 0.05. The results show that phase shifts during jetlag are significantly correlated with photoperiod changes prior to jetlag, photoperiod changes following jetlag, and time zone crossing. Under phase delay schedules, the duration of the jetlag transition is substantially correlated to the photoperiod after jetlag and the time zone crossed. Under phase advance schedules, jetlag recovery time significantly depends on the photoperiod before jetlag and time zone crossed. Abbreviation: PP = photoperiod.

Light schedule factors exert different effects on the jetlag transition time depending on the direction of jetlag. In advanced schedules, both the time zone crossed and the photoperiod at point A significantly impact the jetlag transition time, while during delayed schedules, the time zone crossed and the photoperiod at point B are predicted to play a role. The jetlag transition time generally increases as the time zone crossed increases in both phase advance and phase delay schedules. In addition, our findings suggest that the jetlag transition time shows a significant decrease as the photoperiod at point A decreases in phase advance schedules, and when photoperiod at point B decreases in phase delay schedules.

Discussion

The CTS plays a crucial role in synchronizing an organism’s biological rhythms to the external light/dark cycles, facilitating successful adaptation to the environment. To gain a deeper understanding of the mechanisms underlying this hierarchically organized system, we developed a mathematical model that incorporates several key compartments of the CTS, including the SCN, the HPA axis, and the peripheral clocks. Building upon our previous work (Li and Androulakis, 2021, 2022), this model enabled us to conduct a range of novel investigations, such as characterizing hierarchical entrainment, exploring photoperiod-induced seasonal responses, and analyzing interactions between lighting zone and time zone transitions during jetlag.

Our study aimed to gain a deeper understanding of the mechanisms underlying the synchronization of peripheral clocks with external light/dark cycles. One crucial factor in this process is the circadian rhythm of glucocorticoids, which are regulated by the SCN and act as a synchronizer of peripheral clocks (Balsalobre et al., 2000; Cuesta et al., 2015). Our simulation results showed that the HPA axis, which serves as a “bridge” between the central and peripheral clocks, plays a critical role in synchronizing peripheral biological activities (Figure 2). Specifically, our analysis of the entrainment dynamics of the peripheral clocks revealed that when entrained in vivo, the downstream peripheral compartments exhibit an insensitive response when exposed to light entraining signals with extreme periods. This observation suggests a potential benefit to the organism’s homeostasis (Figure 3). It is worth noting that our findings do not contradict previous studies by Abraham et al. (2010). In their research, the entrainment ranges of both the SCN and peripheral compartments were tested under the direct entrainment of two identical entrainers. However, in our study, the SCN is directly entrained by light, whereas the downstream peripheral compartments are directly entrained by signals originating from upstream compartments. As a result, our results offer complementary insights into the entrainment dynamics of peripheral clocks in vivo, further enriching our understanding of the intrinsic homeostatic mechanisms within the CTS.

By incorporating the photoperiod-related SCN neuronal plasticity, our model exhibits responses to photoperiod variations, with glucocorticoids amplitude rising as photoperiod decreases (Figure 4). This phenomenon has also been observed experimentally in various species, where glucocorticoid levels are highest during winter months and lowest during summer months (Walker et al., 1997; Bilbo et al., 2002; Owen et al., 2005). In temperate regions with significant day length variation, photoperiod serves as a potent entrainer of immune function and physiological plasticity, which may explain some of the variability in disease incidence and severity over the year. For example, winter is generally associated with increased co-regulated proinflammatory gene expression and interleukin-6 (IL6) signaling (Dopico et al., 2015; Pierre et al., 2017). Elevated glucocorticoid levels during winter may be associated with the increased activity of numerous inflammatory immune responses (Pierre et al., 2017). Moreover, the onset of seasonal affective disorder (SAD) is often observed during winter, with symptoms reducing in summer and spring, which could lead to winter depression. Although the mechanisms of SAD are unclear, it is hypothesized that SAD is associated with delayed circadian rhythms during winter (Lewy et al., 2007). Consistent with this hypothesis, our model identified a short-photoperiod-dependent phase delay (Figures 4 and 5), suggesting that the peripheral phase delay during winter months may serve as a mechanism inducing SAD.

According to the model, varied glucocorticoid amplitude is predicted to exhibit differential synchronization effects on peripheral cells (Figure 5). Experimentally, the duration of daily activity in mammals is positively correlated with the photoperiod, with reduced activity during winter to reduce metabolic expenditure (Kenagy, 1973; Daan and Aschoff, 1975). In our model, the peripheral phase is narrowly distributed during winter while spread out during summer, which may lead to broader behavior patterns as a mechanism for adapting to longer photoperiods. It is important to acknowledge that the current model does not take into account all the factors associated with seasonal variations that may impact the periphery, such as fluctuations in temperature, variations in light intensity, changes in stress levels, and the regulatory role of the melatonin (Xu et al., 2018). However, our model highlights the potential limitations of lower amplitude glucocorticoids in synchronizing peripheral dynamics under prolonged photoperiods. In addition, it is hypothesized that high light intensity during summer months may significantly alter entrainment ability on the periphery, and further research incorporating various seasonal-related factors is expected to be conducted in the future to explore this possibility.

During jetlag, peripheral cells can become desynchronized due to the misalignment between biological clocks and social time. Our previous study reported that recovering from jetlag is correlated with an individual’s chronotype and light sensitivity, implying that bright light therapy might alleviate jetlag-associated fatigue (Li and Androulakis, 2021). However, few studies have investigated the impact of photoperiod on jetlag recovery. In this study, we sought to investigate the effect of light exposure length on jetlag recovery time and phase shift (Figures 6 –8). Although the time required for the oscillatory system to transition from one state to another in our model may be longer than what is observed in certain experimental settings, our model considers a more complex network in which different compartments are interconnected, which may result in slower resynchronization times. Our goal is to provide a theoretical framework for comprehending the overall dynamics of the system rather than predicting specific quantitative outcomes. In addition, the faster synchronization rate observed in vivo may be explained by the potential synergistic effects of multiple zeitgebers under physiological conditions. Our simulation results indicate that different interventional strategies are required for eastward (phase advance) and westward (phase delay) jetlag, as significant influencing factors differ between the two. Adjusting to a phase-advanced schedule is more challenging, with a more significant positive correlation between the number of time zones crossed and jetlag transition time compared to that of a phase-delayed schedule.

Notably, our study suggests that in phase-advanced schedules, reducing light exposure before eastward traveling can mitigate jetlag symptoms, while decreasing light exposure after jetlag can accelerate recovery during phase-delayed schedules. Our results also demonstrated a complex positive correlation between intrinsic phase shift and the jetlag transition time (Suppl. Fig. S7). When the phase shift of the host is moderate (below ~10 h), a greater phase shift is likely to result in a longer jetlag transition period. However, as the phase shift increases, jetlag recovery time loses its sensitivity to respond to phase shift changes, indicating a more severe disruption of the circadian rhythm. According to these predictions, to improve jetlag recovery, it is suggested to minimize the intrinsic phase shift by decreasing the light exposure before eastward traveling and after westward traveling. Our predictions have the potential to inspire more experimental designs focused on developing non-pharmacological interventions to alleviate jetlag through optimized light exposure before and after the transition. In addition, further investigation is anticipated to explore other methods, such as adjusting feeding time and exercise (as suggested by Bae and Androulakis, 2017; Bin et al., 2019), for preserving the intrinsic phase of the system.

Supplemental Material

sj-docx-1-jbr-10.1177_07487304231188541 – Supplemental material for The SCN-HPA-Periphery Circadian Timing System: Mathematical Modeling of Clock Synchronization and the Effects of Photoperiod on Jetlag Adaptation

Supplemental material, sj-docx-1-jbr-10.1177_07487304231188541 for The SCN-HPA-Periphery Circadian Timing System: Mathematical Modeling of Clock Synchronization and the Effects of Photoperiod on Jetlag Adaptation by Yannuo Li and Ioannis P. Androulakis in Journal of Biological Rhythms

Supplemental Material

sj-xlsx-2-jbr-10.1177_07487304231188541 – Supplemental material for The SCN-HPA-Periphery Circadian Timing System: Mathematical Modeling of Clock Synchronization and the Effects of Photoperiod on Jetlag Adaptation

Supplemental material, sj-xlsx-2-jbr-10.1177_07487304231188541 for The SCN-HPA-Periphery Circadian Timing System: Mathematical Modeling of Clock Synchronization and the Effects of Photoperiod on Jetlag Adaptation by Yannuo Li and Ioannis P. Androulakis in Journal of Biological Rhythms

Supplemental Material

sj-xlsx-3-jbr-10.1177_07487304231188541 – Supplemental material for The SCN-HPA-Periphery Circadian Timing System: Mathematical Modeling of Clock Synchronization and the Effects of Photoperiod on Jetlag Adaptation

Supplemental material, sj-xlsx-3-jbr-10.1177_07487304231188541 for The SCN-HPA-Periphery Circadian Timing System: Mathematical Modeling of Clock Synchronization and the Effects of Photoperiod on Jetlag Adaptation by Yannuo Li and Ioannis P. Androulakis in Journal of Biological Rhythms

Footnotes

Acknowledgements

The authors acknowledge support from NIH GM131800.

Conflict of Interest Statement

The authors have no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

ORCID iD

Yannuo Li

Data Availability

Simulations and analysis were performed in MATLAB_R2021a (The MathWorks, http://www.mathworks.com). Additional data and materials generated in this work are available in a public GitHub repository at .

Notes

References

Abraham

Granada

Westermark

Heine

Kramer

Herzel

(2010) Coupling governs entrainment range of circadian clocks. Mol Syst Biol 6:438.

Acevedo

Mavroudis

DuBois

Almon

Jusko

Androulakis

(2021) Pathway-level analysis of genome-wide circadian dynamics in diverse tissues in rat and mouse. J Pharmacokinet Pharmacodyn 48:361-374.

Amirat

Khammar

Brudieux

(1980) Seasonal changes in plasma and adrenal concentrations of cortisol, corticosterone, aldosterone, and electrolytes in the adult male sand rat (Psammomys obesus). Gen Comp Endocrinol 40:36-43.

Ananthasubramaniam

Schmal

Herzel

(2020) Amplitude effects allow short jet lags and large seasonal phase shifts in minimal clock models. J Mol Biol 432:3722-3737.

Androulakis

(2021) Circadian rhythms and the HPA axis: a systems view. Wires Mech Dis 13:e1518.

Aton

Herzog

(2005) Come together, right. . .now: synchronization of rhythms in a mammalian circadian clock. Neuron 48:531-534.

Bae

Androulakis

(2017) The synergistic role of light-feeding phase relations on entraining robust circadian rhythms in the periphery. Gene Regul Syst Bio 11:1-21.

Balsalobre

Brown

Marcacci

Tronche

Kellendonk

Reichardt

Schütz

Schibler

(2000) Resetting of circadian time in peripheral tissues by glucocorticoid signaling. Science 289:2344-2347.

Beam

Wasserfall

Woodwyk

Akers

Rauch

Blok

Mason

Vos

Perry

Brusko

(2020) Synchronization of the normal human peripheral immune system: a comprehensive circadian systems immunology analysis. Sci Rep 10:672.

10.

Bernard

Gonze

Čajavec

Herzel

Kramer

(2007) Synchronization-induced rhythmicity of circadian oscillators in the suprachiasmatic nucleus. PLoS Comput Biol 3:e68.

11.

Bilbo

Dhabhar

Viswanathan

Saul

Yellon

Nelson

(2002) Short day lengths augment stress-induced leukocyte trafficking and stress-induced enhancement of skin immune function. Proc Natl Acad Sci U S A 99:4067-4072.

12.

Bin

Ledger

Nour

Postnova

Stamatakis

Cistulli

deChazal

Allman-Farinelli

Caillaud

Bauman

, et al. (2020) How do travelers manage jetlag and travel fatigue? a survey of passengers on long-haul flights. Chronobiol Int 37:1621-1628.

13.

Bin

Postnova

Cistulli

(2019) What works for jetlag? a systematic review of non-pharmacological interventions. Sleep Med Rev 43:47-59.

14.

Bordyugov

Abraham

Granada

Rose

Imkeller

Kramer

Herzel

(2015) Tuning the phase of circadian entrainment. J R Soc Interface 12:20150282.

15.

Chen

Gao

Xiong

Yang

Wang

Jin

(2017) Circadian clock and steroidogenic-related gene expression profiles in mouse Leydig cells following dexamethasone stimulation. Biochem Biophys Res Commun 483:294-300.

16.

Chen

S-K

Badea

Hattar

(2011) Photoentrainment and pupillary light reflex are mediated by distinct populations of ipRGCs. Nature 476:92-95.

17.

Cuesta

Cermakian

Boivin

(2015) Glucocorticoids entrain molecular clock components in human peripheral cells. FASEB J 29:1360-1370.

18.

Daan

Aschoff

(1975) Circadian rhythms of locomotor activity in captive birds and mammals: their variations with season and latitude. Oecologia 18:269-316.

19.

Dopico

Evangelou

Ferreira

Guo

Pekalski

Smyth

Cooper

Burren

Fulford

Hennig

, et al. (2015) Widespread seasonal gene expression reveals annual differences in human immunity and physiology. Nat Commun 6:7000.

20.

Focke

CMB

Iremonger

(2020) Rhythmicity matters: Circadian and ultradian patterns of HPA axis activity. Mol Cell Endocrinol 501:110652.

21.

Gachon

Loizides-Mangold

Petrenko

Dibner

(2017) Glucose homeostasis: regulation by peripheral circadian clocks in rodents and humans. Endocrinology 158:1074-1084.

22.

Garcia-Ojalvo

Elowitz

Strogatz

(2004) Modeling a synthetic multicellular clock: repressilators coupled by quorum sensing. Proc Natl Acad Sci U S A 101:10955-10960.

23.

Geier

Becker-Weimann

Kramer

Herzel

(2005) Entrainment in a model of the mammalian circadian oscillator. J Biol Rhythms 20:83-93.

24.

Herzog

(2007) Neurons and networks in daily rhythms. Nat Rev Neurosci 8:790-802.

25.

Honma

Nakamura

Shirakawa

Honma

K-i

(2004) Diversity in the circadian periods of single neurons of the rat suprachiasmatic nucleus depends on nuclear structure and intrinsic period. Neurosci Lett 358:173-176.

26.

Kalsbeek

Fliers

Hofman

Swaab

Buijs

(2010) Vasopressin and the output of the hypothalamic biological clock. J Neuroendocrinol 22:362-372.

27.

Kalsbeek

van der Spek

Lei

Endert

Buijs

Fliers

(2012) Circadian rhythms in the hypothalamo-pituitary-adrenal (HPA) axis. Mol Cell Endocrinol 349:20-29.

28.

Kenagy

(1973) Daily and seasonal patterns of activity and energetics in a heteromyid rodent community. Ecology 54:1201-1219.

29.

Korenčič

Bordyugov

Košir

Rozman

Goličnik

Herzel

(2012) The interplay of cis-regulatory elements rules circadian rhythms in mouse liver. PLoS One 7:e46835.

30.

Laakso

Porkka-Heiskanen

Alila

Stenberg

Johansson

(1994) Twenty-four-hour rhythms in relation to the natural photoperiod: a field study in humans. J Biol Rhythms 9:283-293.

31.

Lewy

Rough

Songer

Mishra

Yuhas

Emens

(2007) The phase shift hypothesis for the circadian component of winter depression. Dialogues Clin Neurosci 9:291-300.

32.

Androulakis

(2021) Light entrainment of the SCN circadian clock and implications for personalized alterations of corticosterone rhythms in shift work and jet lag. Sci Rep 11:17929.

33.

Androulakis

(2022) Light-induced synchronization of the SCN coupled oscillators and implications for entraining the HPA axis. Front Endocrinol 13:960351.

34.

Mavroudis

Corbett

Calvano

Androulakis

(2014) Mathematical modeling of light-mediated HPA axis activity and downstream implications on the entrainment of peripheral clock genes. Physiol Genomics 46:766-778.

35.

Mavroudis

Scheff

Calvano

Lowry

Androulakis

(2012) Entrainment of peripheral clock genes by cortisol. Physiol Genomics 44:607-621.

36.

Meijer

Michel

VanderLeest

Rohling

(2010) Daily and seasonal adaptation of the circadian clock requires plasticity of the SCN neuronal network. Eur J Neurosci 32:2143-2151.

37.

Moore

Speh

Leak

(2002) Suprachiasmatic nucleus organization. Cell Tissue Res 309:89-98.

38.

Morris

Purvis

Scheer

(2016) Circadian misalignment increases cardiovascular disease risk factors in humans. Proc Natl Acad Sci U S A 113:E1402-E1411.

39.

Nahmias

Androulakis

(2021) Circadian effects of drug responses. Annu Rev Biomed Eng 23:203-224.

40.

Olejniczak

Oster

Ray

(2022) Glucocorticoid circadian rhythms in immune function. Semin Immunopathol 44:153-163.

41.

Ono

Honma

K-i

Honma

(2021) Roles of neuropeptides, VIP and AVP, in the mammalian central circadian clock. Front Neurosci 15:351.

42.

Oster

Challet

Ott

Arvat

de Kloet

Dijk

Lightman

Vgontzas

Van Cauter

(2017) The functional and clinical significance of the 24-hour rhythm of circulating glucocorticoids. Endocr Rev 38:3-45.

43.

Otsuka

Goto

Kawai

Togo

Sato

Katoh

Furuse

Yasuo

(2012) Photoperiod regulates corticosterone rhythms by altered adrenal sensitivity via melatonin-independent mechanisms in Fischer 344 rats and C57BL/6J mice. PLoS One 7:e39090.

44.

Owen

Czekala

Swaisgood

Steinman

Lindburg

(2005) Seasonal and diurnal dynamics of glucocorticoids and behavior in giant pandas. Ursus 16:208-221.

45.

Pierre

Rao

Hartmanshenn

Androulakis

(2018) Modeling the influence of seasonal differences in the HPA axis on synchronization of the circadian clock and cell cycle. Endocrinology 159:1808-1826.

46.

Pierre

Schlesinger

Androulakis

(2017) The hepato-hypothalamic-pituitary-adrenal-renal axis: mathematical modeling of cortisol’s production, metabolism, and seasonal variation. J Biol Rhythms 32:469-484.

47.

Rao

Androulakis

(2017) Modeling the sex differences and interindividual variability in the activity of the hypothalamic-pituitary-adrenal axis. Endocrinology 158:4017-4037.

48.

Rao

Androulakis

(2019) Allostatic adaptation and personalized physiological trade-offs in the circadian regulation of the HPA axis: a mathematical modeling approach. Sci Rep 9:11212.

49.

Reddy

Maywood

Karp

King

Inoue

Gonzalez

Lilley

Kyriacou

Hastings

(2007) Glucocorticoid signaling synchronizes the liver circadian transcriptome. Hepatology 45:1478-1488.

50.

Saper

Chou

Gooley

(2005) The hypothalamic integrator for circadian rhythms. Trends Neurosci 28:152-157.

51.

Schlesinger

Brunetti

Androulakis

(2022) Does seasonality of the microbiota contribute to the seasonality of acute gout flare? Clin Exp Rheumatol 40:1793-1800.

52.

Schmal

Myung

Herzel

Bordyugov

(2015) A theoretical study on seasonality. Front Neurol 6:94.

53.

Bernal

Pillsbury

Yamamoto

Feldman

(2009) Glucocorticoid regulation of the circadian clock modulates glucose homeostasis. Proc Natl Acad Sci U S A 106:17582-17587.

54.

Stanojević

Marković

Čupić Kolar-Anić

ŽL

Vukojević

(2018) Advances in mathematical modelling of the hypothalamic–pituitary–adrenal (HPA) axis dynamics and the neuroendocrine response to stress. Curr Opin Chem Eng 21:84-95.

55.

Tokuda

Schmal

Ananthasubramaniam

Herzel

(2020) Conceptual models of entrainment, jet lag, and seasonality. Front Physiol 11:334.

56.

VanderLeest

Houben

Michel

Deboer

Albus

Vansteensel

Block

Meijer

(2007) Seasonal encoding by the circadian pacemaker of the SCN. Curr Biol 17:468-473.

57.

VanderLeest

Rohling

Michel

Meijer

(2009) Phase shifting capacity of the circadian pacemaker determined by the SCN neuronal network organization. PLoS One 4:e4976.

58.

Vondrasová

Hájek

Illnerová

(1997) Exposure to long summer days affects the human melatonin and cortisol rhythms. Brain Res 759:166-170.

59.

Walker

Best

Noon

Watt

Webb

(1997) Seasonal variation in glucocorticoid activity in healthy men. J Clin Endocrinol Metab 82:4015-4019.

60.

Webb

Angelo

Huettner

Herzog

(2009) Intrinsic, nondeterministic circadian rhythm generation in identified mammalian neurons. Proc Natl Acad Sci U S A 106:16493-16498.

61.

Welsh

Yoo

S-H

Liu

Takahashi

Kay

(2004) Bioluminescence imaging of individual fibroblasts reveals persistent, independently phased circadian rhythms of clock gene expression. Curr Biol 14:2289-2295.

62.

Wieczorek

Blazejczyk

Bullet

Morita

(2013) Lighting characteristics during the polar day and their impact on changes in melatonin secretion. Geogr Pol 86:67-75.

63.

Wieczorek

Blazejczyk

Morita

(2016) Changes in melatonin secretion in tourists after rapid movement to another lighting zone without transition of time zone. Chronobiol Int 33:220-233.

64.

Liu

Guo

(2018) Association of melatonin production with seasonal changes, low temperature, and immuno-responses in hamsters. Molecules 23:703.

65.

Yamamoto

Nakahata

Tanaka

Yoshida

Soma

Shinohara

Yasuda

Mamine

Takumi

(2005) Acute physical stress elevates mouse period1 mRNA expression in mouse peripheral tissues via a glucocorticoid-responsive element. J Biol Chem 280:42036-42043.

66.

Yan

Wang

Liu

Shao

(2008) Analysis of gene regulatory networks in the mammalian circadian rhythm. PLoS Comput Biol 4:e1000193.

67.

Yoo

Yamazaki

Lowrey

Shimomura

Buhr

Siepka

Hong

Yoo

, et al. (2004) PERIOD2::LUCIFERASE real-time reporting of circadian dynamics reveals persistent circadian oscillations in mouse peripheral tissues. Proc Natl Acad Sci U S A 101:5339-5346.

68.

Zhang

Lahens

Ballance

Hughes

Hogenesch

(2014) A circadian gene expression atlas in mammals: implications for biology and medicine. Proc Natl Acad Sci U S A 111:16219-16224.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

1.70 MB

0.01 MB