Circadian rhythms are daily cycles in behavior and physiology which repeat approximately every 24 hours. The master regulator of these rhythms is located in a small part of the brain called the supra-chiasmatic nucleus. This brain structure regulates the timing of sleep and wakefulness and is also thought to control the daily rhythms of cells throughout the body on a molecular level. It does this by synchronizing the activity of a set of genes called clock genes.

Under normal conditions, the levels of proteins coded for by clock genes change throughout the day following a rhythm that matches sleep-wake patterns. However, keeping animals and humans awake at their preferred sleeping times affects the protein levels of clock genes in many tissues of the body. This suggests that, in addition to the supra-chiasmatic nucleus, sleep-wake cycles may also influence clock-gene rhythms throughout the body.

To test this theory, Hoekstra, Jan et al. measured the levels of PERIOD-2, a protein coded for by the clock gene Period-2, while tracking sleep-wake states in mice. They did this by imaging a bioluminescent version of the PERIOD-2 protein in the brain and the kidneys, at the same time as they recorded the brain activity, movement and muscle response of animals. Results showed that PERIOD-2 increased on waking and decreased when mice fell asleep. Additionally, in mice lacking a circadian rhythm in sleep-wake behavior – whose changes in PERIOD-2 levels with respect to time were greatly reduced – imposing a regular sleep-wake cycle restored normal PERIOD-2 rhythmicity.

Next, Hoekstra, Jan et al. developed a mathematical model to understand how sleep-wake cycles together with circadian rhythms affect clock-gene activity in the brain and kidneys. Computer simulations suggested that sleep-wake cycles and circadian factors act as forces of comparable strength driving clock-gene dynamics. Both need to act in concert to keep clock-genes rhythmic. The model also predicted the large and immediate effects of sleep deprivation on PERIOD-2 levels, giving further credence to the idea that waking accelerated clock-gene rhythms while sleeping slowed them down. Modelling also suggested that having regular clock-gene rhythms protects against sleep disturbances.

In summary, this work shows how sleep patterns contribute to the daily rhythms in clock genes in the brain and body. The findings support the idea that well-timed sleep-wake schedules could help people to adjust to new time zones. It might also be useful to inform other strategies to reduce the health impacts of shift work.

The sleep-wake distribution is coordinated by the interaction of a circadian and a homeostatic process (Daan et al., 1984). The biological substrates underlying the circadian process are relatively well understood: circadian rhythms in overt behavior of mammals are generated by the suprachiasmatic nuclei (SCN) located in the hypothalamus (Hastings et al., 2018). At the molecular level, so-called ‘clock genes’ interact through negative transcriptional/translational feedback loops (TTFLs), where the CLOCK/NPAS2:ARNTL (BMAL1) heterodimers drive the transcription of their target genes, among them the Period (Per1,2) and Cryptochrome (Cry1,2) genes. Subsequently, PER and CRY proteins assemble into repressor complexes that inhibit CLOCK/NPAS2:ARNTL-mediated transcription, including their own. The resulting reduction of the repressor complex allows a new cycle to start. This feedback loop, present in almost each cell of the mammalian body, interacts with other molecular pathways, together ensuring a period of ~24 hr (Hastings et al., 2018). The SCN synchronizes peripheral clock gene expression rhythms through its rhythmic behavioral, electrical, and humoral output generated across the day (Schibler et al., 2015).

Accumulating evidence suggests that, perhaps surprisingly, clock genes are also involved in the homeostatic aspect of sleep regulation (Franken, 2013). This is illustrated by the sleep deprivation (SD)-induced increase in the expression of the clock gene Per2 in tissues peripheral to the SCN, including the cerebral cortex, liver, and kidney (Curie et al., 2013; Curie et al., 2015; Franken et al., 2007; Maret et al., 2007). Moreover, the highest level of peripheral Per2 expression is reached after the time of day mice were awake most, suggesting that also during spontaneous periods of waking Per2 expression accumulates. Accordingly, lesioning of the SCN, which eliminates the circadian sleep-wake distribution, attenuates the circadian amplitude of clock gene transcripts and proteins in peripheral tissues (Akhtar et al., 2002; Curie et al., 2015; Tahara et al., 2012; Saini et al., 2013; Sinturel et al., 2021). Together, these studies suggest that sleeping and waking are important contributors to clock gene expression, but dissecting the contribution of the sleep-wake distribution and circadian time is challenging because the two change in parallel.

By simultaneously recording electroencephalogram (EEG), electromyogram (EMG), locomotor activity (LMA), and PER2-dependent bioluminescence signals from cortex and kidney in freely behaving mice, we established that the circadian sleep-wake distribution importantly contributes to the daily rhythmic changes in central and peripheral PER2 levels. To further test this hypothesis, we predicted that (i) in behaviorally arrhythmic SCN-lesioned (SCNx) mice, daily recurring SDs mimicking a circadian sleep-wake distribution will temporarily reinstate high-amplitude PER2 bioluminescence rhythms, and (ii) in intact rhythmic animals, reducing the amplitude of the circadian sleep-wake distribution will result in a reduced amplitude of PER2 rhythms. While daily SDs indeed enhanced the amplitude of peripheral PER2 rhythms in SCNx mice, the protocol used to reduce the amplitude of the sleep-wake distribution did not reduce PER2 amplitude in all mice. To reconcile the sleep-wake-driven and circadian aspects of PER2 dynamics, we implemented a mathematical model in which waking represents a force that sets in motion a harmonic oscillator describing PER2 dynamics and found that the sleep-wake distribution, also under undisturbed conditions, is an important contributor to the daily changes in PER2 bioluminescence. Moreover, we discovered a second, SCN and sleep-wake independent force with a circadian period that underlay the residual circadian PER2 rhythms in SCNx mice, and that the phase relationship between these two forces is important for predicting the amplitude response in PER2 rhythms to sleep-wake perturbations.

To quantify PER2 levels, we used mice expressing a knock-in (KI) construct encoding a fused PER2::LUCIFERASE (PER2::LUC) protein and in which changes in emitted bioluminescence can be used as proxy for changes in PER2 protein levels (Yoo et al., 2004). Per2^Luc KI mice have been used to follow clock gene expression in vivo (Curie et al., 2015; Ohnishi et al., 2014; Tahara et al., 2012; van der Vinne et al., 2018). However, in these studies, mice had to be anesthetized for each measurement, while in the set-up used in our study (RT-Biolumicorder; Saini et al., 2013), we assessed PER2 bioluminescence continuously in freely moving mice, in central and peripheral tissues. For the central quantification of bioluminescence, we used mice in which the Per^Luc construct was back-crossed onto a C57BL/6J (B6) background (see Materials and methods). For the experiments that quantified bioluminescence in the periphery, we used hairless SKH1 mice carrying the Per2^Luc KI construct because lack of fur allows for the unobstructed measurement of emitted photons (see also Figure 1—figure supplement 1 for experimental design and Figure 1—figure supplement 2A for imaging of bioluminescence in SKH1 mice). Under standard LD12:12 conditions, SKH1 mice exhibited sleep-wake patterns characteristic of mice, that is, during the light phase they spent more time in both non-rapid eye movement (NREM) sleep and REM sleep relative to the dark phase, the latter being their habitual active phase. Moreover, they showed the typical sleep homeostatic response to a 6 hr SD during the first 18 hr of recovery, both in sleep time and EEG delta power, although the increase in REM sleep did not reach significance levels (Figure 1—figure supplement 3).

In three pilot experiments, we optimized our experimental set-up. We established that the most important source contributing to the peripheral bioluminescence signals in the SKH1 mice are the kidneys (Figure 1—figure supplement 2A). We have previously shown that the central bioluminescence signal obtained in B6 mice is of cortical origin (Curie et al., 2015). To accomplish this, luciferin was infused directly into the brain, the skull locally thinned and equipped with a glass cone, and animals were not shaved, thereby preventing passage of photons from the periphery (see Materials and methods). We confirmed that in B6 mice thus prepared only photons emitted by the brain could be detected (Figure 1—figure supplement 2A). Finally, by using mice expressing luciferase under the control of synthetic CAG promotor (CAG-Luc mice; Cao et al., 2004), we determined that its substrate luciferin is best delivered through implantable osmotic mini-pumps compared to administration through the drinking water. Under the latter condition, strong daily rhythms in bioluminescence were observed, likely as a result of rhythms in drinking behavior, thereby driving luciferin availability (Figure 1—figure supplement 2B).

Sleep-wake state affects PER2 bioluminescence

It now has been well documented that enforced wakefulness affects Per2 mRNA and protein levels in various tissues and mammalian species (Franken, 2013; Hoekstra et al., 2019; Möller-Levet et al., 2013; Vassalli and Franken, 2017), but it is not known whether circadian rhythms in spontaneous sleep-wake behavior contribute to the daily changes in PER2 levels. To address this question, we equipped Per2^Luc KI B6 and SKH1 mice (n = 6 and 5, respectively) with wireless EEG/EMG recorders (i.e., NeuroLoggers) and monitored simultaneously sleep-wake state, PER2 bioluminescence, and LMA under constant darkness (DD; see Figure 1—figure supplement 1). Time asleep in mice kept in the RT-Biolumicorder quantified with the NeuroLoggers under DD was similar to that quantified under LD conditions using a tethered EEG acquisition system (48 hr baseline in DD vs. LD in C57BL/6J mice: NREM sleep: 42.2% ± 1.6%, REM sleep 6.1% ± 0.3% of total recording time; n = 6; compared to 41.3 ± 1.1 and 5.2% ± 0.2%, respectively; n = 12, data taken from Diessler et al., 2018). Sleep and EEG have not been recorded in SKH1 mice previously. Also in these mice, similar sleep durations were obtained under the two conditions (48 hr baseline in DD vs. LD: NREM sleep: 40.4% ± 2.1%, REM sleep 5.8% ± 0.9%; n = 5; compared to 40.0% ± 1.0% and 6.9% ± 0.2%, respectively; n = 8; see Figure 1—figure supplement 3).

Next, we reproduced the circadian changes in PER2 bioluminescence as well as the response to a 6 hr SD (Figure 1A) as described previously (Curie et al., 2015) in both peripheral bioluminescence, mainly of renal origin, and central bioluminescence, emitted by the cortex (see Materials and methods for explanation of the SD procedure). Similar to what was observed in that publication, SD elicited a tissue-specific response, with an immediate increase in central PER2 bioluminescence after SD, while in the periphery the response was delayed and PER2-bioluminescence increases were observed in the second and fifth hour of recovery (Figure 1A). In addition to these direct effects on PER2 bioluminescence within the first 6 hr of recovery, the SD also caused a long-term reduction of rhythm amplitude during both recovery days; that is, the first (REC1) and last (REC2) 24 hr period of recovery. In the periphery, this decrease amounted to ca. 30% on both days (REC1: 17.5–41.1%; REC2: 18.7–42.4%, 95% confidence intervals [95% CI]). In the central recordings, PER2-bioluminescence amplitude decreased significantly during REC1 (32% [10.2–54.5%]), whereas the decrease during REC2 (22% [0.06–44.8%]) no longer reached significance levels (linear mixed model with fixed conditional effect [‘BSL,’ ‘REC1,’ ‘REC2’] and random intercept effect [‘Mouse’]; periphery: BSL vs. REC1 p=0.0014; vs. REC2 p=0.0011; central: BSL vs. REC1 p=0.022; vs. REC2 p=0.082; sinewave-fitted baseline amplitudes: periphery 0.258 [0.15–0.36]; central 0.193 [0.08–0.31 a.u.]). This long-term reduction of PER2-biolumnescence is reminiscent of the long-term SD effects on rhythm amplitude we observed for Per2 expression and for other clock genes in the cortex (Hor et al., 2019).

Figure 1 with 7 supplements see all

Download asset Open asset

Figure 1—source data 1

This file contains the numerical values on which the graphs in Figure 1A–C are based.

https://cdn.elifesciences.org/articles/69773/elife-69773-fig1-data1-v2.xlsx

Download elife-69773-fig1-data1-v2.xlsx

Although a circadian modulation of PER2 bioluminescence of both peripheral and central origin is evident (see Figure 1B and Figure 1—figure supplement 4), we observed additional changes in bioluminescence that occurred simultaneously with changes in sleep-wake state, indicating that PER2 bioluminescence increases during wake-dominated periods and decreases during sleep-dominated periods in both tissues. To quantify this observation, transitions from sleep (irrespective of sleep state) to wake and from wake to sleep were selected (see Materials and methods for selection criteria). Examples of the selected transitions are indicated in Figure 1B as a hypnogram. A similar number of sleep-to-wake and wake-to-sleep transitions passed selection criteria during the two-and-a-half baseline days (periphery: 31.2 ± 3.9 and 29.0 ± 3.9; central: 25.6 ± 1.4 and 24.6 ± 1.4, respectively; mean ± SEM). Transitions obtained in mice in which central bioluminescence was recorded were shorter (longest common wake period after sleep-to-wake transitions: 39 vs. 57 min; longest common sleep period after wake-to-sleep transitions: 63 vs. 75 min, for central and periphery, respectively). Although PER2 bioluminescence increased during wakefulness and decreased during sleep in both tissues, tissue-specific differences were observed (Figure 1C). At sleep-to-wake transitions, tissue differences concerned an initial decrease in peripheral PER2 bioluminescence after wake onset, followed by a steep increase saturating at 125%. In contrast, central levels of PER2 bioluminescence increased from the start and followed a linear time course throughout the waking period. Despite these different dynamics, similar bioluminescence levels were reached in both tissues. After wake-to-sleep transitions, peripheral PER2 bioluminescence initially increased before quickly decreasing and then leveling out at around 84% (Figure 1C). In contrast, the central signal decreased linearly throughout the sleep period reaching levels of 91% at the end.

Changes in sleep-wake state are accompanied by physiological alterations, such as changes in body and brain temperature (Sela et al., 2020; Vishwakarma et al., 2021), which could influence bioluminescence by changing substrate levels and/or the rate of the enzymatic reaction. Although the circadian rhythms of subcutaneous temperature and PER2 bioluminescence are ca. 4 hr out of phase (Figure 1—figure supplement 6), this does not exclude the possibility that fast changes in physiology associated with sleep-wake transitions contribute to luciferin availability. Therefore, we assessed sleep-wake-related changes in mice carrying two other luciferase reporter constructs. Besides the CAG-Luc mice we used to decide on the route of luciferin administration (Figure 1—figure supplement 2B), we also had access to Pkg1-Luc mice in which bioluminescence is under the control of the promotor of the housekeeping gene Pkg1 (see Materials and methods). As expected from a housekeeping gene, bioluminescence in Pkg1-Luc mice was relatively constant and did not increase at sleep-wake transitions (Figure 1—figure supplement 5). In contrast, in CAG-Luc mice increases at transitions were larger than in Per2^Luc mice, consistent with the rapidly activated CMV-immediate-early enhancer contained within the CAG synthetic promotor (Brightwell et al., 1997; Collaco and Geusz, 2003). These results show that the bioluminescence changes at sleep-wake transitions are specific to the luciferase reporter construct and therefore not an artifact of sleep-wake-related changes in physiology.

Taken together, three observations indicate that the sleep-wake distribution importantly contributes to both central and peripheral PER2-bioluminescence dynamics: (1) enforced wakefulness has both acute and long-term effects on PER2 bioluminescence; (2) the changes in PER2 bioluminescence at sleep-wake transitions demonstrate that also spontaneous waking is associated with an increase, and sleep with a decrease, in PER2 bioluminescence; and (3) that PER2 bioluminescence is high when mice are spontaneously awake more.

The central PER bioluminescence signals tended to be noisier than peripheral signals (signal-to-noise ratio: –1.15 ± 1.0 and 0.37 ± 1.0 dB, respectively; estimated on 3 min values in baseline of the 6 hr SD experiments according to Leise et al., 2012). Moreover, the peripheral signal was easier to obtain and minimally invasive to the animal. We therefore decided for the next experiments to focus on peripheral PER2 bioluminescence using the Per2^Luc SKH1 mice. Because waking correlates highly with LMA (Figure 1—figure supplement 7), we recorded LMA as proxy for wakefulness to avoid invasive EEG/EMG surgery and mice having to adapt to a ca. 2.7 g recorder mounted on the head.

Repeated sleep deprivations in SCNx mice temporarily reinstate a circadian rhythm in PER2 bioluminescence

Lesioning the SCN rendered LMA arrhythmic under DD conditions (Figure 2—figure supplement 1). In arrhythmic SCNx mice, we confirmed that rhythms in PER2 bioluminescence were strongly reduced, but not completely abolished (Figure 2A–C). During the second baseline measurement after SCNx (BSL2 in Figure 2B), we observed in two of the four mice erratic high values (see also Figure 1—figure supplement 2) that we cannot readily explain especially because in the subsequent four SD days the variance in the PER2 bioluminescence signal among mice was again as small as in the earlier recordings (Figure 2A vs. B). The erratic values in these two mice led to a >10-fold increase in the residual sum of squares (RSS), indicating a poorer fit, and to a higher amplitude of the fitted sinewaves in the second compared to the first baseline recording (Figure 2C, BSL1 vs. BSL2). The repeated SDs induced a robust oscillation in the PER2-bioluminescence signal, restoring amplitude to the levels observed prior to SCN lesioning and higher than those observed in the baseline recordings (Figure 2C). The latter observation shows that the increased amplitude during the SD did not result from an SD-mediated alignment of different individual phases in baseline. Importantly, the oscillation continued after the end of the SDs while decreasing in amplitude (Figure 2C).

Figure 2 with 2 supplements see all

Download asset Open asset

Figure 2—source data 1

This file contains the numerical values on which the graphs in Figure 2A–C are based.

https://cdn.elifesciences.org/articles/69773/elife-69773-fig2-data1-v2.xlsx

Download elife-69773-fig2-data1-v2.xlsx

The 2hOnOff protocol reduces the circadian sleep-wake amplitude without consistently modulating PER2 dynamics

In the next experiment, we aimed at reducing the amplitude of the circadian sleep-wake distribution in intact mice using the 2-day 2hOnOff protocol. We measured sleep-wake state and PER2 bioluminescence during baseline conditions, the 2hOnOff procedure, and the two recovery days in two cohorts of mice. In the first cohort, PER2 bioluminescence was monitored, and mice were taken out of the RT-Biolumicorder for the 2 hr SDs. Bioluminescence was therefore quantified only during the SOWs, thus biasing the read-out of PER2 bioluminescence during the 2hOnOff protocol towards levels reached during sleep. A second cohort was implanted with tethered EEG/EMG electrodes to determine the efficacy of the 2hOnOff protocol in reducing the sleep-wake distribution amplitude.

As expected, mice exhibited a circadian PER2-bioluminescence rhythm under baseline conditions (Figure 3A). Contrary to expectation, the 2hOnOff protocol did not significantly decrease the ongoing circadian PER2 oscillation when analyzed at the group level (Figure 3B). When inspecting the individual responses, four out of the six mice did show a consistent 27% reduction (range 23–31%) in PER2-biolumininescence amplitude during the 2-day 2hOnOff protocol, while in the remaining two mice amplitude increased by 40% (Figure 3C). In all six mice, PER2-bioluminescence amplitude reverted to baseline values during recovery, irrespective of whether it was increased or decreased during the 2hOnOff protocol. The distinct bimodal, opposing response in amplitude observed among individual mice and the subsequent reverting back to baseline suggest that the 2hOnOff protocol did affect the ongoing PER2-bioluminescence rhythm. This phenomenon might relate to factors not accounted for in our experimental design. In the modeling section at the end of the Results section, we explore this observation further and find that variation in circadian phase could underlie the individual differential response (for details, see ‘Modeling PER2 dynamics’).

Figure 3

Download asset Open asset

Figure 3—source data 1

This file contains the numerical values on which the graphs in Figure 3A–E are based.

https://cdn.elifesciences.org/articles/69773/elife-69773-fig3-data1-v2.xlsx

Download elife-69773-fig3-data1-v2.xlsx

We used the second cohort of mice to determine the efficacy with which the 2hOnOff protocol reduced the circadian sleep-wake distribution. In contrast to bioluminescence, the amplitude of the circadian sleep-wake distribution did decrease consistently in all mice (Figure 3D and E). The circadian rhythm in the sleep-wake distribution was, however, not eliminated because the sleep obtained during the 2 hr SDs as well as during the 2 hr SOWs both varied as a function of time of day (sleep during SDs: average: 12.4%, min-max: 4.9–23.0%; during SOWs: average: 54.3%, min-max: 35.7–76.4%). This was especially evident at the beginning of the subjective light phase of day 2 when the average time spent asleep during the SD reached 23%. The sleep obtained during the SDs at this time could be due to a substantial sleep pressure because of lost sleep during day 1 of the 2hOnOff protocol (total time spent asleep/day, mean ± SEM; baseline: 11.1 ± 0.1 hr; 2hOnOff: day 1: 7.3 ± 0.3 hr, day 2: 8.6 ± 0.7 hr, paired two-tailed t-test, BSL vs. 2hOnOff-day 1: t(7) = 8.86, p<0.001; BSL vs. 2hOnOff-day 2: t(7) = 4.18, p=0.004; 2hOnOff-day 1 vs. -day 2, t(7) = 1.97, p=0.09), combined with the difficulty for the experimenters to visually detect and prevent sleep in pinkish hairless mice under dim-red light conditions. During recovery from the 2hOnOff protocol, baseline amplitudes were re-established (Figure 3E).

Modeling circadian PER2 dynamics

The results above demonstrate that both circadian and sleep-wake-dependent factors need to be considered when studying PER2 dynamics. To understand how these two factors collectively generate the variance in PER2 bioluminescence, we put our experimental data into a theoretical framework and modeled changes in peripheral PER2 bioluminescence according to a driven harmonic oscillator (Curie et al., 2013). This type of oscillator is not self-sustained but depends on rhythmic forces to set in motion and maintain the oscillation. In the model, we assume the sleep-wake distribution to represent one such force. In the absence of rhythmic forces, the oscillator loses energy, resulting in a gradual reduction of its amplitude according to the rate set by its damping constant (γ). Besides amplitude and timing (i.e., phase) of the recurring forces and the damping constant, the oscillatory system is further defined by a string constant (${\omega }_0^2$) with the natural angular frequency (${\omega }_0$) defining the intrinsic period. Because waking was not quantified alongside bioluminescence in the SCNx and the 2hOnOff experiments, we used LMA as a proxy for the driving force provided by wakefulness (${\displaystyle {\overrightarrow{F}}_{WAKE}}$). The 6 hr SD experiment showed that LMA and time spent awake are strongly correlated and that their hourly dynamics closely changed in parallel (Figure 1—figure supplement 7). Moreover, applying the model to this data set (see below) using either wakefulness or LMA as ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ yielded similar fits (Figure 4—figure supplement 1), demonstrating that LMA is an appropriate proxy for wakefulness for the purpose of our model. We performed the analyses at the group level; that is, the mean LMA level of all mice was used to reflect ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ , and the free parameters in the model were optimized by fitting the motion of the oscillator to the mean PER2-bioluminescence levels (see Materials and methods). Besides ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ , we implemented a circadian force to account for the residual PER2-bioluminescence rhythm observed in the SCNx mice (Figure 2). We assumed that this additional SCN and sleep-wake-independent force reflects a peripheral circadian process (${\overrightarrow{F}}_{PERI}$) present in both intact and SCNx mice (Sinturel et al., 2021). This force was modeled as a sinewave with amplitude and phase as free parameters in both experiments while period was set to the period estimated from the baseline PER2-bioluminescence rhythm observed in the SCNx mice (23.7 ± 0.6 hr; n = 4). A schematic overview of the influence of ${\overrightarrow{F}}_{PERI}$ and ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ on PER2 bioluminescence is presented in Figure 4A.

Figure 4 with 2 supplements see all

Download asset Open asset

Figure 4—source data 1

This file contains the numerical values on which the graphs in Figure 4B–D are based.

https://cdn.elifesciences.org/articles/69773/elife-69773-fig4-data1-v2.xlsx

Download elife-69773-fig4-data1-v2.xlsx

We first optimized the parameters of the model describing PER2-bioluminescence dynamics in the SCNx experiment (data from Figure 2), and then predicted PER2 bioluminescence under the 2hOnOff experiment. ${\overrightarrow{F}}_{PERI}$ ’s amplitude and phase again required optimization as both differed between the two experiments (nonoverlapping 95% CI for both in Table 1). With the parameters listed in Table 1, the model captured the dynamic changes in PER2 bioluminescence in the SCNx experiment with high precision including the residual rhythmicity in baseline, the reinstated pronounced rhythmicity during the four SDs, and its subsequent dampening thereafter (Figure 4B, black line, Table 2). The model also accurately captured average bioluminescence dynamics in the 2hOnOff experiment (Figure 4C, black line, Table 2). It furthermore predicted an 18% reduction of PER2-bioluminescence amplitude during the 2hOnOff protocol compared to baseline (relative amplitudes estimated by the model: 0.31 vs. 0.38 [a.u.]; Figure 4C), consistent with our hypotheses and the 27% reduction in PER2-bioluminescence amplitude observed in the four mice in which amplitude did decrease (see Figure 3C). To further test the performance of the model, we predicted the peripheral PER2-bioluminescence data obtained in the 6 hr SD experiment using the parameters optimized for the 2hOnOff experiment (Figure 1A, Table 1). Also, the results of this experiment could be accurately predicted with the model, including the higher PER2-bioluminescence levels reached over the initial recovery hours after SD and the subsequent longer-term reduction in amplitude (Figures 1A and 4D, Table 2). The model could not reliably predict the central PER2-bioluminescence dynamics mainly due to an earlier phase of the central compared to the peripheral signal (Figure 4—figure supplement 2D). Moreover, the increase in PER2 levels immediately following the SD (Figure 1A) was missed by the model, further underscoring the tissue-specific relationship between time-spent-awake and PER2 requiring the model to be optimized according to the tissue under study.

Table 1

Parameters	SCNx	2hOnOff
${\displaystyle \gamma }$(damping constant)	0.0155 (0.0103–0.0223) (hr–1)
${\omega }_0$(natural angular frequency)	0.288 (0.285–0.291) (rad * hr–1)
$\beta$(coefficient ${\overrightarrow{F}}_{WAKE}$)	6.81e-5 (4.48e-5–9.01e-5)
Model intercept	0.92 (0.90–0.95)	0.91 (0.89–0.92)
$A$(amplitude ${\overrightarrow{F}}_{PERI}$)	2.92e-3 (2.58e-3–3.32e-3)	3.87e-3 (3.80e-3–4.00e-3)
$\phi$(phase ${\overrightarrow{F}}_{PERI}$)	4.73 (4.6–4.85)	2.78 (2.78–2.82)

Table 2

		BIC	RSS	Bayes factor(model vs. flat)
SCNx	$\overrightarrow{F}$=${\overrightarrow{F}}_{WAKE}$	–287.0	11.94	2.47e32
	$\overrightarrow{F}$=${\overrightarrow{F}}_{PERI}$	–224.5	13.60	6.48e18
	$\overrightarrow{F}$=${\overrightarrow{F}}_{WAKE}+{\overrightarrow{F}}_{PERI}$	–427.5	8.36	7.87e62
2hOnOff	$\overrightarrow{F}$=${\overrightarrow{F}}_{WAKE}$	–28.5	11.96	1.82e24
	$\overrightarrow{F}$=${\overrightarrow{F}}_{PERI}$	–536.8	1.29	4.41e86
	$\overrightarrow{F}$=${\overrightarrow{F}}_{WAKE}+{\overrightarrow{F}}_{PERI}$	–540.1	1.27	2.27e87
6 hr SD	$\overrightarrow{F}$=${\overrightarrow{F}}_{WAKE}$	18.2	92.50	9.15e-17
	$\overrightarrow{F}$=${\overrightarrow{F}}_{PERI}$	–205.9	4.16	5.80e48
	$\overrightarrow{F}$=${\overrightarrow{F}}_{WAKE}+{\overrightarrow{F}}_{PERI}$	–311.9	2.46	6.23e71

The model accurately captured peripheral PER2-bioluminescence dynamics under three different experimental conditions. The model’s robust performance prompted us to assess in silico whether (1) both forces are required to predict PER2-bioluminescence dynamics, (2) lesioning the SCN affects the forces exerted on PER2 bioluminescence under undisturbed, DD conditions, and (3) differences in the circadian phase could predict the opposing response in PER2-bioluminescence amplitude among mice in the 2hOnOff experiment.

To evaluate if ${\overrightarrow{F}}_{PERI}$ and ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ combined are required to predict PER2 dynamics, we compared the performance of the model by dropping these two modeling terms; that is, removing either ${\overrightarrow{F}}_{PERI}$ or ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ . Removing ${\overrightarrow{F}}_{PERI}$ from the model did not change the coefficient obtained for ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ (8.25e-5, which is still within the 95% CI estimated in the full model; see Table 1). However, this simpler model (only four free parameters to estimate compared to six in the full model) could not reliably predict the bioluminescence data of the SCNx and 2hOnOff experiments (Figure 4B and C, blue lines) and fit statistics indicated a poorer fit (i.e., a higher Bayesian information criterion [BIC] score, see Table 2). For instance, without ${\overrightarrow{F}}_{PERI}$ the model is unable to capture the residual PER2 rhythmicity in the baseline of the SCNx experiment (Figure 4B, ‘model 1’ until 72 hr), and it took longer for the PER2 bioluminescence to reach high-amplitude levels compared to the full model because of ${\overrightarrow{F}}_{PERI}$ being in phase with the timing of the SDs in the full model. With the same strategy, we evaluated the model’s performance when dropping ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ (Figure 4B and C, green lines). As the sleep-wake distribution no longer affects the model, ${\overrightarrow{F}}_{PERI}$ dynamics are unperturbed throughout the experiment and can therefore be solved as a sinewave in the differential equation. Fit statistics for the SCNx experiment were poor compared to the full model as the effects of SDs on PER2 amplitude could not be captured (Figure 4B, Table 2). Removing ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ when modeling the 2hOnOff experiment resulted in BIC scores similar to those obtained with the full model (Table 2), although the amplitude reduction during the 2hOnOff protocol could not be captured (Figure 4C). Also, the PER2-bioluminescence dynamics in the 6 hr SD experiment were captured best when using both ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ and ${\displaystyle {\overrightarrow{F}}_{PERI}}$ (Figure 4D, Table 2). Together, these results demonstrate that both forces are required to accurately predict PER2 dynamics under the tested experimental conditions.

To answer the second question, we compared the forces driving PER2 bioluminescence in the baseline recordings of the 2hOnOff (intact mice) and the SCNx experiments. The model estimated that the average absolute force (${\displaystyle {\overrightarrow{F}}_{WAKE}}$ , ${\displaystyle {\overrightarrow{F}}_{PERI}}$ , $\overrightarrow{F}\gamma$, and $\overrightarrow{F}{\omega }_0^2$ combined) exerted on the relative PER2-bioluminescense levels was more than twice as high in intact mice compared to SCNx mice (2.11e-2 vs. 0.90e-2 h-2). This large difference was not due to a difference in ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ , which was somewhat lower in intact mice (2hOnOff vs. SCNx; 2.01e-3 vs. 2.22e-3 h-2). Although ${\overrightarrow{F}}_{PERI}$ was higher by 32% (2.45e-3 vs. 1.85e-3 h-2), this difference did not substantially contribute to the higher amplitude of the PER2 rhythm in intact mice. To illustrate this, we substituted the value of ${\overrightarrow{F}}_{PERI}$ obtained in the 2hOnOff baseline with the lower value obtained in the SCNx baseline, which resulted in a 17% amplitude reduction of the PER2 rhythm. The presence of a circadian sleep-wake distribution impacted PER2 amplitude to a larger extent: running the simulation with the parameters estimated for the 2hOnOff experiment but with the sleep-wake distribution of the SCNx mice resulted in a 32% reduction. Therefore, the circadian sleep-wake organization is an important contributor to high-amplitude oscillations in PER2. In the 2hOnOff experiment, the larger PER2 momentum resulted in a four times higher string and damping forces (2hOnOff vs. SCNx; $\overrightarrow{F}{\omega }_0^2$ : 19.9e-3 vs. 4.4e-3 h-2; $\overrightarrow{F}\gamma$: 9.74e-4 vs. 2.52e-4 h-2) that together with the larger ${\overrightarrow{F}}_{PERI}$ underlie the larger average absolute force. Moreover, this analysis demonstrated that the direct effects of ${\displaystyle {\overrightarrow{F}}_{WAKE}}$ and ${\overrightarrow{F}}_{PERI}$ on PER2 bioluminescence do not depend on an intact SCN and that the magnitude of the two forces is comparable.

Although the model predicted the expected decrease in PER2-bioluminescense rhythm amplitude in the 2hOnOff experiment, this decrease was not observed in the mean bioluminescence data. The individual data showed, however, that in four mice this intervention did decrease PER2-bioluminescence amplitude, while in the remaining two amplitude increased (Figure 3C). With the model, we addressed the third question: Do circadian phase differences predict the opposite effects of the 2hOnOff intervention on amplitude? Surprisingly, when systematically varying the phase of ${\overrightarrow{F}}_{PERI}$ during the baseline prior to the 2hOnOff intervention (Figure 4—figure supplement 2C), we found that at phase advances larger than 2.5 hr, the model predicted an increase of PER2 amplitude during the subsequent 2hOnOff protocol instead of a decrease (illustrated in Figure 4—figure supplement 2A for a 6 hr phase advance). Circadian phase differences among animals might relate to individual differences in period length of their free-running rhythms that accumulate over time. We did not find evidence for a difference in PER2 bioluminescence or LMA phase at the end of the baseline recording (2.5 days under DD) between animals that showed a decrease in PER2-bioluminescence amplitude compared to animals that showed an increase. However, during the first day of recovery following the 2hOnOff protocol (5.5 days under DD), we observed a ca. 2 hr phase advance in bioluminescence and 1 hr phase advance in LMA in the two mice that increased their amplitude during the preceding 2hOnOff protocol (Figure 4—figure supplement 2B, left) compared to the four mice for which we observed the anticipated decrease in PER2 amplitude (Figure 4—figure supplement 2B, right). Whether these differences in phase contributed to the opposite response cannot be answered with the current data set. Nevertheless, the model yielded a perhaps counterintuitive but testable hypothesis by demonstrating that phase angle between the sleep-wake distribution and peripheral circadian clock-gene rhythms is an important variable in predicting outcome and emphasizes the importance of carefully controlling the initial conditions.

In this study, we assessed the contribution of the sleep-wake distribution to circadian peripheral PER2 rhythmicity. We presented four key findings, supporting the notion that sleep-wake state is indeed an important factor in determining the circadian amplitude of peripheral changes in PER2: (1) spontaneous sustained bouts of waking and sleep were associated with increased and decreased PER2 bioluminescence, respectively; (2) a single SD acutely increased PER2 bioluminescence and during subsequent recovery dampened its circadian amplitude; (3) repeated SDs temporarily reinstated robust rhythmicity in PER2 bioluminescence in behaviorally arrhythmic mice; and (4) mathematical modeling suggests that PER2 dynamics is best understood as a harmonic oscillator driven by two forces: a sleep-waking-dependent force (${\displaystyle {\overrightarrow{F}}_{WAKE}}$) and an SCN-independent and sleep-wake-independent, circadian peripheral force (${\displaystyle {\overrightarrow{F}}_{PERI}}$).

Data underlying the experimental figures is available through the source files.

1. 1. Akhtar RA
   2. Reddy AB
   3. Maywood ES
   4. Clayton JD
   5. King VM
   6. Smith AG
   7. Gant TW
   8. Hastings MH
   9. Kyriacou CP

   (2002) Circadian cycling of the mouse liver transcriptome, as revealed by cDNA microarray, is driven by the suprachiasmatic nucleus

   Current Biology 12:540–550.

   https://doi.org/10.1016/s0960-9822(02)00759-5

   • PubMed
   • Google Scholar
2. 1. Andrews RV

   (1968) Temporal secretory responses of cultured hamster adrenals

   Comparative Biochemistry and Physiology 26:179–193.

   https://doi.org/10.1016/0010-406x(68)90324-1

   • PubMed
   • Google Scholar
3. 1. Bainier C
   2. Mateo M
   3. Felder-Schmittbuhl M-P
   4. Mendoza J

   (2017) Circadian rhythms of hedonic drinking behavior in mice

   Neuroscience 349:229–238.

   https://doi.org/10.1016/j.neuroscience.2017.03.002

   • PubMed
   • Google Scholar
4. 1. Baker FC
   2. Angara C
   3. Szymusiak R
   4. McGinty D

   (2005) Persistence of sleep-temperature coupling after suprachiasmatic nuclei lesions in rats

   American Journal of Physiology. Regulatory, Integrative and Comparative Physiology 289:R827–R838.

   https://doi.org/10.1152/ajpregu.00093.2005

   • PubMed
   • Google Scholar
5. 1. Balsalobre A
   2. Brown SA
   3. Marcacci L
   4. Tronche F
   5. Kellendonk C
   6. Reichardt HM
   7. Schütz G
   8. Schibler U

   (2000) Resetting of circadian time in peripheral tissues by glucocorticoid signaling

   Science 289:2344–2347.

   https://doi.org/10.1126/science.289.5488.2344

   • PubMed
   • Google Scholar
6. 1. Bass J
   2. Takahashi JS

   (2010) Circadian integration of metabolism and energetics

   Science 330:1349–1354.

   https://doi.org/10.1126/science.1195027

   • PubMed
   • Google Scholar
7. 1. Brightwell G
   2. Poirier V
   3. Cole E
   4. Ivins S
   5. Brown KW

   (1997) Serum-dependent and cell cycle-dependent expression from a cytomegalovirus-based mammalian expression vector

   Gene 194:115–123.

   https://doi.org/10.1016/s0378-1119(97)00178-9

   • PubMed
   • Google Scholar
8. 1. Cao YA
   2. Wagers AJ
   3. Beilhack A
   4. Dusich J
   5. Bachmann MH
   6. Negrin RS
   7. Weissman IL
   8. Contag CH

   (2004) Shifting foci of hematopoiesis during reconstitution from single stem cells

   PNAS 101:221–226.

   https://doi.org/10.1073/pnas.2637010100

   • PubMed
   • Google Scholar
9. 1. Cao R
   2. Gkogkas CG
   3. de Zavalia N
   4. Blum ID
   5. Yanagiya A
   6. Tsukumo Y
   7. Xu H
   8. Lee C
   9. Storch K-F
   10. Liu AC
   11. Amir S
   12. Sonenberg N

   (2015) Light-regulated translational control of circadian behavior by eIF4E phosphorylation

   Nature Neuroscience 18:855–862.

   https://doi.org/10.1038/nn.4010

   • PubMed
   • Google Scholar
10. 1. Cederroth CR
    2. Park JS
    3. Basinou V
    4. Weger BD
    5. Tserga E
    6. Sarlus H
    7. Magnusson AK
    8. Kadri N
    9. Gachon F
    10. Canlon B

    (2019) Circadian Regulation of Cochlear Sensitivity to Noise by Circulating Glucocorticoids

    Current Biology 29:2477–2487.

    https://doi.org/10.1016/j.cub.2019.06.057

    • PubMed
    • Google Scholar
11. 1. Cheon S
    2. Park N
    3. Cho S
    4. Kim K

    (2013) Glucocorticoid-mediated Period2 induction delays the phase of circadian rhythm

    Nucleic Acids Research 41:6161–6174.

    https://doi.org/10.1093/nar/gkt307

    • PubMed
    • Google Scholar
12. 1. Chong SYC
    2. Ptáček LJ
    3. Fu YH

    (2012) Genetic insights on sleep schedules: this time, it’s PERsonal

    Trends in Genetics 28:598–605.

    https://doi.org/10.1016/j.tig.2012.08.002

    • PubMed
    • Google Scholar
13. 1. Cirelli C
    2. Tononi G

    (2000) Differential expression of plasticity-related genes in waking and sleep and their regulation by the noradrenergic system

    The Journal of Neuroscience 20:9187–9194.

    https://doi.org/10.1523/JNEUROSCI.20-24-09187.2000

    • PubMed
    • Google Scholar
14. 1. Collaco AM
    2. Geusz ME

    (2003) Monitoring immediate-early gene expression through firefly luciferase imaging of HRS/J hairless mice

    BMC Physiology 3:8.

    https://doi.org/10.1186/1472-6793-3-8

    • Google Scholar
15. 1. Cuesta M
    2. Cermakian N
    3. Boivin DB

    (2015) Glucocorticoids entrain molecular clock components in human peripheral cells

    FASEB Journal 29:1360–1370.

    https://doi.org/10.1096/fj.14-265686

    • PubMed
    • Google Scholar
16. 1. Curie T
    2. Mongrain V
    3. Dorsaz S
    4. Mang GM
    5. Emmenegger Y
    6. Franken P

    (2013) Homeostatic and circadian contribution to EEG and molecular state variables of sleep regulation

    Sleep 36:311–323.

    https://doi.org/10.5665/sleep.2440

    • PubMed
    • Google Scholar
17. 1. Curie T
    2. Maret S
    3. Emmenegger Y
    4. Franken P

    (2015) In Vivo Imaging of the Central and Peripheral Effects of Sleep Deprivation and Suprachiasmatic Nuclei Lesion on PERIOD-2 Protein in Mice

    Sleep 38:1381–1394.

    https://doi.org/10.5665/sleep.4974

    • PubMed
    • Google Scholar
18. 1. Daan S
    2. Beersma DG
    3. Borbély AA

    (1984) Timing of human sleep: recovery process gated by a circadian pacemaker

    The American Journal of Physiology 246:R161–R183.

    https://doi.org/10.1152/ajpregu.1984.246.2.R161

    • PubMed
    • Google Scholar
19. 1. Damiola F
    2. Le Minh N
    3. Preitner N
    4. Kornmann B
    5. Fleury-Olela F
    6. Schibler U

    (2000) Restricted feeding uncouples circadian oscillators in peripheral tissues from the central pacemaker in the suprachiasmatic nucleus

    Genes & Development 14:2950–2961.

    https://doi.org/10.1101/gad.183500

    • PubMed
    • Google Scholar
20. 1. Dickmeis T

    (2009) Glucocorticoids and the circadian clock

    The Journal of Endocrinology 200:3–22.

    https://doi.org/10.1677/JOE-08-0415

    • PubMed
    • Google Scholar
21. 1. Diessler S
    2. Jan M
    3. Emmenegger Y
    4. Guex N
    5. Middleton B
    6. Skene DJ
    7. Ibberson M
    8. Burdet F
    9. Götz L
    10. Pagni M
    11. Sankar M
    12. Liechti R
    13. Hor CN
    14. Xenarios I
    15. Franken P

    (2018) A systems genetics resource and analysis of sleep regulation in the mouse

    PLOS Biology 16:e2005750.

    https://doi.org/10.1371/journal.pbio.2005750

    • PubMed
    • Google Scholar
22. 1. Du NH
    2. Arpat AB
    3. De Matos M
    4. Gatfield D

    (2014) MicroRNAs shape circadian hepatic gene expression on a transcriptome-wide scale

    eLife 3:e02510.

    https://doi.org/10.7554/eLife.02510

    • Google Scholar
23. 1. D’Alessandro M
    2. Beesley S
    3. Kim JK
    4. Jones Z
    5. Chen R
    6. Wi J
    7. Kyle K
    8. Vera D
    9. Pagano M
    10. Nowakowski R
    11. Lee C

    (2017) Stability of Wake-Sleep Cycles Requires Robust Degradation of the PERIOD Protein

    Current Biology 27:3454–3467.

    https://doi.org/10.1016/j.cub.2017.10.014

    • PubMed
    • Google Scholar
24. 1. Edgar DM
    2. Dement WC
    3. Fuller CA

    (1993) Effect of SCN lesions on sleep in squirrel monkeys: evidence for opponent processes in sleep-wake regulation

    The Journal of Neuroscience 13:1065–1079.

    https://doi.org/10.1523/JNEUROSCI.13-03-01065.1993

    • PubMed
    • Google Scholar
25. 1. Eide EJ
    2. Woolf MF
    3. Kang H
    4. Woolf P
    5. Hurst W
    6. Camacho F
    7. Vielhaber EL
    8. Giovanni A
    9. Virshup DM

    (2005) Control of mammalian circadian rhythm by CKIepsilon-regulated proteasome-mediated PER2 degradation

    Molecular and Cellular Biology 25:2795–2807.

    https://doi.org/10.1128/MCB.25.7.2795-2807.2005

    • PubMed
    • Google Scholar
26. 1. Engeland WC
    2. Massman L
    3. Mishra S
    4. Yoder JM
    5. Leng S
    6. Pignatti E
    7. Piper ME
    8. Carlone DL
    9. Breault DT
    10. Kofuji P

    (2018) The Adrenal Clock Prevents Aberrant Light-Induced Alterations in Circadian Glucocorticoid Rhythms

    Endocrinology 159:3950–3964.

    https://doi.org/10.1210/en.2018-00769

    • PubMed
    • Google Scholar
27. 1. Franken P
    2. Dudley CA
    3. Estill SJ
    4. Barakat M
    5. Thomason R
    6. O’Hara BF
    7. McKnight SL

    (2006) NPAS2 as a transcriptional regulator of non-rapid eye movement sleep: genotype and sex interactions

    PNAS 103:7118–7123.

    https://doi.org/10.1073/pnas.0602006103

    • PubMed
    • Google Scholar
28. 1. Franken P
    2. Thomason R
    3. Heller HC
    4. O’Hara BF

    (2007) A non-circadian role for clock-genes in sleep homeostasis: a strain comparison

    BMC Neuroscience 8:87.

    https://doi.org/10.1186/1471-2202-8-87

    • PubMed
    • Google Scholar
29. 1. Franken P

    (2013) A role for clock genes in sleep homeostasis

    Current Opinion in Neurobiology 23:864–872.

    https://doi.org/10.1016/j.conb.2013.05.002

    • PubMed
    • Google Scholar
30. 1. Gerber A
    2. Esnault C
    3. Aubert G
    4. Treisman R
    5. Pralong F
    6. Schibler U

    (2013) Blood-borne circadian signal stimulates daily oscillations in actin dynamics and SRF activity

    Cell 152:492–503.

    https://doi.org/10.1016/j.cell.2012.12.027

    • PubMed
    • Google Scholar
31. 1. Hall MP
    2. Woodroofe CC
    3. Wood MG
    4. Que I
    5. Van’t Root M
    6. Ridwan Y
    7. Shi C
    8. Kirkland TA
    9. Encell LP
    10. Wood KV
    11. Löwik C
    12. Mezzanotte L

    (2018) Click beetle luciferase mutant and near infrared naphthyl-luciferins for improved bioluminescence imaging

    Nature Communications 9:132.

    https://doi.org/10.1038/s41467-017-02542-9

    • PubMed
    • Google Scholar
32. 1. Hastings MH
    2. Maywood ES
    3. Brancaccio M

    (2018) Generation of circadian rhythms in the suprachiasmatic nucleus

    Nature Reviews. Neuroscience 19:453–469.

    https://doi.org/10.1038/s41583-018-0026-z

    • PubMed
    • Google Scholar
33. 1. Hayasaka N
    2. Hirano A
    3. Miyoshi Y
    4. Tokuda IT
    5. Yoshitane H
    6. Matsuda J
    7. Fukada Y

    (2017) Salt-inducible kinase 3 regulates the mammalian circadian clock by destabilizing PER2 protein

    eLife 6:e24779.

    https://doi.org/10.7554/eLife.24779

    • PubMed
    • Google Scholar
34. 1. Hoekstra MM
    2. Emmenegger Y
    3. Hubbard J
    4. Franken P

    (2019) Cold-inducible RNA-binding protein (CIRBP) adjusts clock-gene expression and REM-sleep recovery following sleep deprivation

    eLife 8:e43400.

    https://doi.org/10.7554/eLife.43400

    • PubMed
    • Google Scholar
35. 1. Hor CN
    2. Yeung J
    3. Jan M
    4. Emmenegger Y
    5. Hubbard J
    6. Xenarios I
    7. Naef F
    8. Franken P

    (2019) Sleep-wake-driven and circadian contributions to daily rhythms in gene expression and chromatin accessibility in the murine cortex

    PNAS 116:25773–25783.

    https://doi.org/10.1073/pnas.1910590116

    • PubMed
    • Google Scholar
36. 1. Iwano S
    2. Sugiyama M
    3. Hama H
    4. Watakabe A
    5. Hasegawa N
    6. Kuchimaru T
    7. Tanaka KZ
    8. Takahashi M
    9. Ishida Y
    10. Hata J
    11. Shimozono S
    12. Namiki K
    13. Fukano T
    14. Kiyama M
    15. Okano H
    16. Kizaka-Kondoh S
    17. McHugh TJ
    18. Yamamori T
    19. Hioki H
    20. Maki S
    21. Miyawaki A

    (2018) Single-cell bioluminescence imaging of deep tissue in freely moving animals

    Science 359:935–939.

    https://doi.org/10.1126/science.aaq1067

    • PubMed
    • Google Scholar
37. 1. Kofuji P
    2. Mure LS
    3. Massman LJ
    4. Purrier N
    5. Panda S
    6. Engeland WC

    (2016) Intrinsically Photosensitive Retinal Ganglion Cells (ipRGCs) Are Necessary for Light Entrainment of Peripheral Clocks

    PLOS ONE 11:e0168651.

    https://doi.org/10.1371/journal.pone.0168651

    • PubMed
    • Google Scholar
38. 1. Kornmann B
    2. Schaad O
    3. Bujard H
    4. Takahashi JS
    5. Schibler U

    (2007) System-driven and oscillator-dependent circadian transcription in mice with a conditionally active liver clock

    PLOS Biology 5:e34.

    https://doi.org/10.1371/journal.pbio.0050034

    • PubMed
    • Google Scholar
39. 1. Le Minh N
    2. Damiola F
    3. Tronche F
    4. Schütz G
    5. Schibler U

    (2001) Glucocorticoid hormones inhibit food-induced phase-shifting of peripheral circadian oscillators

    The EMBO Journal 20:7128–7136.

    https://doi.org/10.1093/emboj/20.24.7128

    • PubMed
    • Google Scholar
40. 1. Leise TL
    2. Wang CW
    3. Gitis PJ
    4. Welsh DK

    (2012) Persistent Cell-Autonomous Circadian Oscillations in Fibroblasts Revealed by Six-Week Single-Cell Imaging of PER2::LUC Bioluminescence

    PLOS ONE 7:e33334.

    https://doi.org/10.1371/journal.pone.0033334

    • Google Scholar
41. 1. Liu J
    2. Zou X
    3. Gotoh T
    4. Brown AM
    5. Jiang L
    6. Wisdom EL
    7. Kim JK
    8. Finkielstein CV

    (2018) Distinct control of PERIOD2 degradation and circadian rhythms by the oncoprotein and ubiquitin ligase MDM2

    Science Signaling 11:eaau0715.

    https://doi.org/10.1126/scisignal.aau0715

    • Google Scholar
42. 1. Mang GM
    2. Franken P

    (2012) Sleep and EEG Phenotyping in Mice

    Current Protocols in Mouse Biology 2:55–74.

    https://doi.org/10.1002/9780470942390.mo110126

    • PubMed
    • Google Scholar
43. 1. Mang GM
    2. Franken P

    (2015) Genetic dissection of sleep homeostasis

    Current Topics in Behavioral Neurosciences 25:25–63.

    https://doi.org/10.1007/7854_2013_270

    • PubMed
    • Google Scholar
44. 1. Maret S
    2. Dorsaz S
    3. Gurcel L
    4. Pradervand S
    5. Petit B
    6. Pfister C
    7. Hagenbuchle O
    8. O’Hara BF
    9. Franken P
    10. Tafti M

    (2007) Homer1a is a core brain molecular correlate of sleep loss

    PNAS 104:20090–20095.

    https://doi.org/10.1073/pnas.0710131104

    • PubMed
    • Google Scholar
45. 1. Masuda S
    2. Narasimamurthy R
    3. Yoshitane H
    4. Kim JK
    5. Fukada Y
    6. Virshup DM

    (2020) Mutation of a PER2 phosphodegron perturbs the circadian phosphoswitch

    PNAS 117:10888–10896.

    https://doi.org/10.1073/pnas.2000266117

    • PubMed
    • Google Scholar
46. 1. Möller-Levet CS
    2. Archer SN
    3. Bucca G
    4. Laing EE
    5. Slak A
    6. Kabiljo R
    7. Lo JCY
    8. Santhi N
    9. von Schantz M
    10. Smith CP
    11. Dijk DJ

    (2013) Effects of insufficient sleep on circadian rhythmicity and expression amplitude of the human blood transcriptome

    PNAS 110:E1132–E1141.

    https://doi.org/10.1073/pnas.1217154110

    • PubMed
    • Google Scholar
47. 1. Mongrain V
    2. Hernandez SA
    3. Pradervand S
    4. Dorsaz S
    5. Curie T
    6. Hagiwara G
    7. Gip P
    8. Heller HC
    9. Franken P

    (2010) Separating the contribution of glucocorticoids and wakefulness to the molecular and electrophysiological correlates of sleep homeostasis

    Sleep 33:1147–1157.

    https://doi.org/10.1093/sleep/33.9.1147

    • PubMed
    • Google Scholar
48. 1. Moore RY
    2. Eichler VB

    (1972) Loss of a circadian adrenal corticosterone rhythm following suprachiasmatic lesions in the rat

    Brain Research 42:201–206.

    https://doi.org/10.1016/0006-8993(72)90054-6

    • PubMed
    • Google Scholar
49. 1. Nash JC
    2. Varadhan R

    (2011) Unifying Optimization Algorithms to Aid Software System Users: optimx for R

    Journal of Statistical Software 43:i09.

    https://doi.org/10.18637/jss.v043.i09

    • Google Scholar
50. 1. Ohnishi N
    2. Tahara Y
    3. Kuriki D
    4. Haraguchi A
    5. Shibata S
    6. Ebihara S

    (2014) Warm Water Bath Stimulates Phase-Shifts of the Peripheral Circadian Clocks in PER2::LUCIFERASE Mouse

    PLOS ONE 9:e100272.

    https://doi.org/10.1371/journal.pone.0100272

    • Google Scholar
51. 1. Ohsaki K
    2. Oishi K
    3. Kozono Y
    4. Nakayama K
    5. Nakayama KI
    6. Ishida N

    (2008) The role of {beta}-TrCP1 and {beta}-TrCP2 in circadian rhythm generation by mediating degradation of clock protein PER2

    Journal of Biochemistry 144:609–618.

    https://doi.org/10.1093/jb/mvn112

    • PubMed
    • Google Scholar
52. 1. Ono D
    2. Honma K
    3. Honma S

    (2015) Circadian and ultradian rhythms of clock gene expression in the suprachiasmatic nucleus of freely moving mice

    Scientific Reports 5:12310.

    https://doi.org/10.1038/srep12310

    • Google Scholar
53. 1. Reischl S
    2. Vanselow K
    3. Westermark PO
    4. Thierfelder N
    5. Maier B
    6. Herzel H
    7. Kramer A

    (2007) Beta-TrCP1-mediated degradation of PERIOD2 is essential for circadian dynamics

    Journal of Biological Rhythms 22:375–386.

    https://doi.org/10.1177/0748730407303926

    • PubMed
    • Google Scholar
54. 1. Saini C
    2. Morf J
    3. Stratmann M
    4. Gos P
    5. Schibler U

    (2012) Simulated body temperature rhythms reveal the phase-shifting behavior and plasticity of mammalian circadian oscillators

    Genes & Development 26:567–580.

    https://doi.org/10.1101/gad.183251.111

    • PubMed
    • Google Scholar
55. 1. Saini C
    2. Liani A
    3. Curie T
    4. Gos P
    5. Kreppel F
    6. Emmenegger Y
    7. Bonacina L
    8. Wolf JP
    9. Poget YA
    10. Franken P
    11. Schibler U

    (2013) Real-time recording of circadian liver gene expression in freely moving mice reveals the phase-setting behavior of hepatocyte clocks

    Genes & Development 27:1526–1536.

    https://doi.org/10.1101/gad.221374.113

    • PubMed
    • Google Scholar
56. 1. Satinoff E
    2. Prosser RA

    (1988) Suprachiasmatic nuclear lesions eliminate circadian rhythms of drinking and activity, but not of body temperature, in male rats

    Journal of Biological Rhythms 3:1–22.

    https://doi.org/10.1177/074873048800300101

    • PubMed
    • Google Scholar
57. Conference

    1. Schibler U
    2. Gotic I
    3. Saini C
    4. Gos P
    5. Curie T
    6. Emmenegger Y
    7. Sinturel F
    8. Gosselin P
    9. Gerber A
    10. Fleury-Olela F
    11. Rando G
    12. Demarque M
    13. Franken P

    (2015) Cold Spring Harbor Symposia on Quantitative Biology

    Clock-Talk: Interactions between Central and Peripheral Circadian Oscillators in Mammals. pp. 223–232.

    https://doi.org/10.1101/sqb.2015.80.027490

    • Google Scholar
58. 1. Schiedner G
    2. Hertel S
    3. Kochanek S

    (2000) Efficient transformation of primary human amniocytes by E1 functions of Ad5: Generation of new cell lines for adenoviral vector production

    Human Gene Therapy 11:2105–2116.

    https://doi.org/10.1089/104303400750001417

    • Google Scholar
59. 1. Segall LA
    2. Amir S

    (2010) Exogenous corticosterone induces the expression of the clock protein, PERIOD2, in the oval nucleus of the bed nucleus of the stria terminalis and the central nucleus of the amygdala of adrenalectomized and intact rats

    Journal of Molecular Neuroscience 42:176–182.

    https://doi.org/10.1007/s12031-010-9375-4

    • PubMed
    • Google Scholar
60. Preprint

    1. Sela Y
    2. Hoekstra MMB
    3. Franken P

    (2020) A Sub-Minute Resolution Prediction of Brain Temperature Based on Sleep-Wake State in the Mouse

    bioRxiv.

    https://doi.org/10.1101/2020.08.12.246405

    • Google Scholar
61. 1. Sinturel F
    2. Gos P
    3. Petrenko V
    4. Hagedorn C
    5. Kreppel F
    6. Storch KF
    7. Knutti D
    8. Liani A
    9. Weitz C
    10. Emmenegger Y
    11. Franken P
    12. Bonacina L
    13. Dibner C
    14. Schibler U

    (2021) Circadian hepatocyte clocks keep synchrony in the absence of a master pacemaker in the suprachiasmatic nucleus or other extrahepatic clocks

    Genes & Development 35:329–334.

    https://doi.org/10.1101/gad.346460.120

    • Google Scholar
62. 1. So AYL
    2. Bernal TU
    3. Pillsbury ML
    4. Yamamoto KR
    5. Feldman BJ

    (2009) Glucocorticoid regulation of the circadian clock modulates glucose homeostasis

    PNAS 106:17582–17587.

    https://doi.org/10.1073/pnas.0909733106

    • PubMed
    • Google Scholar
63. 1. Tahara Y
    2. Kuroda H
    3. Saito K
    4. Nakajima Y
    5. Kubo Y
    6. Ohnishi N
    7. Seo Y
    8. Otsuka M
    9. Fuse Y
    10. Ohura Y
    11. Komatsu T
    12. Moriya Y
    13. Okada S
    14. Furutani N
    15. Hirao A
    16. Horikawa K
    17. Kudo T
    18. Shibata S

    (2012) In vivo monitoring of peripheral circadian clocks in the mouse

    Current Biology 22:1029–1034.

    https://doi.org/10.1016/j.cub.2012.04.009

    • PubMed
    • Google Scholar
64. 1. Tamaru T
    2. Hattori M
    3. Honda K
    4. Benjamin I
    5. Ozawa T
    6. Takamatsu K

    (2011) Synchronization of circadian Per2 rhythms and HSF1-BMAL1:CLOCK interaction in mouse fibroblasts after short-term heat shock pulse

    PLOS ONE 6:e24521.

    https://doi.org/10.1371/journal.pone.0024521

    • PubMed
    • Google Scholar
65. 1. Travnickova-Bendova Z
    2. Cermakian N
    3. Reppert SM
    4. Sassone-Corsi P

    (2002) Bimodal regulation of mPeriod promoters by CREB-dependent signaling and CLOCK/BMAL1 activity

    PNAS 99:7728–7733.

    https://doi.org/10.1073/pnas.102075599

    • PubMed
    • Google Scholar
66. 1. van der Vinne V
    2. Swoap SJ
    3. Vajtay TJ
    4. Weaver DR

    (2018) Desynchrony between brain and peripheral clocks caused by CK1δ/ε disruption in GABA neurons does not lead to adverse metabolic outcomes

    PNAS 115:E2437–E2446.

    https://doi.org/10.1073/pnas.1712324115

    • PubMed
    • Google Scholar
67. 1. Vassalli A
    2. Franken P

    (2017) Hypocretin (orexin) is critical in sustaining theta/gamma-rich waking behaviors that drive sleep need

    PNAS 114:E5464–E5473.

    https://doi.org/10.1073/pnas.1700983114

    • PubMed
    • Google Scholar
68. 1. Vishwakarma LC
    2. Sharma B
    3. Singh V
    4. Jaryal AK
    5. Mallick HN

    (2021) Acute sleep deprivation elevates brain and body temperature in rats

    Journal of Sleep Research 30:e13030.

    https://doi.org/10.1111/jsr.13030

    • Google Scholar
69. 1. Wisor JP
    2. O’Hara BF
    3. Terao A
    4. Selby CP
    5. Kilduff TS
    6. Sancar A
    7. Edgar DM
    8. Franken P

    (2002) A role for cryptochromes in sleep regulation

    BMC Neuroscience 3:20.

    https://doi.org/10.1186/1471-2202-3-20

    • Google Scholar
70. 1. Wu T
    2. Ni Y
    3. Kato H
    4. Fu Z

    (2010) Feeding-induced rapid resetting of the hepatic circadian clock is associated with acute induction of Per2 and Dec1 transcription in rats

    Chronobiology International 27:1–18.

    https://doi.org/10.3109/07420520903398625

    • PubMed
    • Google Scholar
71. 1. Yan L
    2. Silver R

    (2002) Differential induction and localization of mPer1 and mPer2 during advancing and delaying phase shifts

    The European Journal of Neuroscience 16:1531–1540.

    https://doi.org/10.1046/j.1460-9568.2002.02224.x

    • PubMed
    • Google Scholar
72. 1. Yasenkov R
    2. Deboer T

    (2010) Circadian regulation of sleep and the sleep EEG under constant sleep pressure in the rat

    Sleep 33:631–641.

    https://doi.org/10.1093/sleep/33.5.631

    • PubMed
    • Google Scholar
73. 1. Yoo SH
    2. Yamazaki S
    3. Lowrey PL
    4. Shimomura K
    5. Ko CH
    6. Buhr ED
    7. Siepka SM
    8. Hong HK
    9. Oh WJ
    10. Yoo OJ
    11. Menaker M
    12. Takahashi JS

    (2004) PERIOD2::LUCIFERASE real-time reporting of circadian dynamics reveals persistent circadian oscillations in mouse peripheral tissues

    PNAS 101:5339–5346.

    https://doi.org/10.1073/pnas.0308709101

    • Google Scholar
74. 1. Zhang B
    2. Gao Y
    3. Li Y
    4. Yang J
    5. Zhao H

    (2016) Sleep Deprivation Influences Circadian Gene Expression in the Lateral Habenula

    Behavioural Neurology 2016:7919534.

    https://doi.org/10.1155/2016/7919534

    • Google Scholar

Decision letter after peer review:

[Editors’ note: the authors submitted for reconsideration following the decision after peer review. What follows is the decision letter after the first round of review.]

Thank you for submitting your work entitled "The sleep-wake distribution contributes to the peripheral rhythms in PERIOD-2" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and a Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Hiroki Ueda (Reviewer #1).

Our decision has been reached after consultation between the reviewers. Based on these discussions and the individual reviews below, we regret to inform you that your work will not be considered further for publication in eLife.

As you can see from the feedback provided by the reviewers they found the data interesting and thought-provoking, yet they also raised substantial concerns about different aspects of the work. Importantly, the reviewers noted that some of the main conclusions (based on the presented data) are not fully supported. A considerable amount of work that is beyond the scope of a revision at eLife would be needed to address these concerns.

Reviewer #1:

In this paper entitled "The sleep-wake distribution contributes to the peripheral rhythms in PERIOD-2", the authors observed PER2 expression in peripheral tissues of freely behaving animals and found that the sleep-wake distribution contributes to the peripheral PER2 dynamics. In addition, implementing the mathematical model, they showed that the peripheral PER2 rhythms can be generated by two factors: a sleep-wake state dependent one and a peripheral (not master) circadian dependent one. Despite the difficulty to dissect the effect of these two factors in vivo, the authors, dealing with this problem by the unique experimental designs, provided direct evidence that sleep-wake state itself affects peripheral daily rhythms. However, this paper has a couple of points to be improved, which are summarized in the comments below.

1. As for Figure 1B, 1C and 1D, the authors showed PER2 bioluminescence changes according to circadian rhythm and state transitions in only one individual mouse (ID: 2MP82), which might be insufficient for concluding that the PER2 bioluminescence changes depending on the circadian and sleep-wake cycle. Plus, there might be a mistake in the time label of Figure 1B (is this plot started from time = 6 or time = 0?). Thus, the reviewer suggests that the authors should add data about other mice and perform the statistical analyses to clarify the validity of the circadian- and sleep-wake-dependent change of PER2 bioluminescence. With Figure 4, the reviewer suggests the same, i.e., the authors should analyze not only one representative mouse but a group of mice.

2. In lines 202 to 205, the authors concluded that the behavioral rhythm reduction by 2hOnOff protocol appeared to be insufficient to consistently reduce the amplitude of the PER2-bioluminescence oscillation. However, it is difficult to follow its logic toward the conclusion. Thus, the reviewer recommends that the authors should describe their logic in more detail.

3. As for Figure 4 and Table 2, the authors compared two mathematical models and concluded that the two forces model has better prediction accuracy. However, it is unclear whether the circadian force itself or the combination of the two forces is important to predict PER2-bioluminescence dynamics. The reviewer suggests that the authors should add a comparison with a model that using a single circadian force (F = FPERI).

4. In lines 277 to 280, the authors described the increase in PER2 amplitude during the 2hOnOff experiment by advancing the phase of FPERI by 6h. The reviewer suggests that the authors should clearly explain why they apply their mathematical model to the 6h advancing, although they observed 1-2h phase advance in their experiments.

Reviewer #2:

By measuring EEG/EMG and PER2::LUC bioluminescence (most signal comes from the kidney) from freely behaving mice simultaneously or from separate cohort mice with similar conditions, the authors observed that sleep state acutely changed PER2::LUC bioluminescence, repeated daily sleep deprivation reinstated circadian PER2::LUC rhythm in SCN ablated behaviorally arrhythmic mice. The authors also did a mathematical simulation modeling peripheral PER2::LUC bioluminescence influenced by a sleep dependent force and a SCN independent force.

I believe the data presented in the current manuscript, however, I don't think the main conclusions are supported by the data. The acute changes (within 30 min) of luminescence they observed is most likely caused by in vivo factors (Not changed in PER2::LUC protein, artifact). It is very problematic to assess the amplitude change from averaged luminescence record from multiple behaviorally arrhythmic mice (analysis flaw). Although the above mentioned two major conclusions must be removed from the manuscript, the study validates previously published results done by ex vivo tissue and in vivo anesthetized animal with luminescence measure from freely behaving mice. So I encourage the authors to publish the data to another journal (please also see my specific points below), but I don't think discovery reported in the work is novel for the publication in eLife.

1. It is unlikely that the acute change of PER2::LUC luminescence reflects PER2 protein translation and degradation. As the authors refer to correctly, luminescence from PER2::LUC knockin mice reflects PER2::LUC fusion protein (not Per2 transcription). The authors didn't correctly refer to the studies measuring light induction of Per2 mRNA and protein in the SCN. The light induction of Per2 mRNA in the SCN can only be seen 60 min after onset of the light pulse (no induction was observed at 30 min after the light pulse. Takumi et al. 1998 and many other publications). Increased PER2 protein in the SCN can only be seen at 4 h after the light pulse (no change was observed 2 h after the light pulse. Yan and Sliver 2004, Antoun et al., 2012). Therefore, the acute luminescence change observed at sleep wake transitions (Figure 1D) is unlikely reflecting changes in the PER2::LUC fusion protein. Temperature, pH, and Ca/Mg and ATP levels all strongly affect luminescence produced by firefly luciferase and luciferin. It was shown that the ultradian rhythm of the core body temperature closely matched the ultradian behavior rhythm (Blessing 2018 PMID: 30454601). The authors referred to their own published study that luminescence rhythm change was vilified by qPCR. However, the time scale of their published work showed that changes occurred after 6 h of sleep deprivation, not in the 30 min range.

2. The advantage of continuous luminescence measurement from the same mouse is that the presence and absence of the rhythm can be determined from a mouse that has no phase marker (e.g. SCN lesioned mice). Arrhythmicity (or lowered amplitude of the rhythm) observed by conventional population sampling (from multiple mice) cannot be determined if individual mice lose rhythmicity, or if individual mice retained intact rhythm but they are out of sync. From the data presented in Figure 2, the authors claimed repeated sleep deprivations reinstated the rhythm, but they didn't analyze amplitude change before, during, and after sleep deprivation in individual mice (even they didn't provide individual data anywhere else, they provided individual data in other experiments). It is well known that the SCN-lesioned mice do not show a daily rhythm under LD cycle (Light cannot reinstate behavior rhythm). The authors stated that they performed daily sleep deprivation to the mice 5 weeks after SCN-lesioned. From average luminescence rhythms of 4 SCN-X mice, the authors cannot conclude their main claim that sleep deprivation reinstated the peripheral circadian rhythm. It is possible there is no change in the amplitude of individual peripheral circadian oscillators, but repeated sleep deprivation entrained those rhythms (the authors averaged individual luminescence relative to the phase of daily sleep deprivation).

3. In the complementary experiment presented in Figure 3, the authors analyzed the amplitude changes from individual mice, but didn't analyze from the averaged data. The authors observed a reduced amplitude by 2hOnOff procedure in 4 out of 6 mice, and an increased amplitude in 2 out of 6 mice. But the authors concluded that the 2hOnOff procedure did not "consistently" affect the rhythm (in another place, they wrote "2hOnOff protocol didn't noticeably affect the ongoing circadian PER2 oscillation").

4. Because of the issues pointed out above, the parameters used for the mathematical simulation are not valid. Use of locomotor activity as a wake parameter is also not appropriate. Mice can be wake without moving.

5. In the method section, the authors described each procedure. However, it was impossible for me to know how exactly each individual mouse went through the procedures. It seems some sleep deprivation was done under dim red light and other was done under light. I suggest the authors provide a timeline (e.g. flowchart figure) of procedures in each experiment.

6. Since in vivo bioluminescence recording became popular technology, I strongly suggest the authors move the data presented in supplementary Figure 2 to the main text as an actual figure. Were CAG mice also crossed with hairless strain? Does tissue explanted from this mouse in vitro not show a circadian rhythm? It even looks like luminescence from the mice with luciferin supplemented by osmic pump shows a weak rhythm. I wonder if CAG promoter may have circadian rhythm as CMV promoter does (Collaco et al., 2005). Amplitude of activity rhythm is lowered in mice with pump. Is this due to surgery (pump implantation)? If so, lowered amplitude of luminescence rhythm might also be affected by surgery. Having a control (saline in pump and luciferin supplied through drinking water) is necessary.

Reviewer #3:

This is an interesting study. As far as I know, there are no previous studies performing simultaneous measurement of sleep-wake state and bioluminescence (clock gene expression) in freely moving mice. In addition, their measurement method was well validated through several pilot studies, and measurements obtained using this method are therefore reliable. In fact, the presented data are clear and convincing. However, although I have no doubt that the results are useful and valuable for a wide range of circadian scientists, I feel the present version of their work may not be a strong candidate for publication in this journal because of insufficient novelty and impact.

1) The results clearly revealed the temporal relationship between sleep-wake state and PER2 bioluminescence through simultaneous monitoring in freely moving mice, with high time resolution. However, because previous studies have already found that sleep deprivation results in a change in PER2 levels, I feel that the present study does not make obvious progress toward further understanding the mechanism. The simultaneous monitoring experiments are interesting but each of the techniques was originally established by other groups. For example, some findings related to the identification of systemic or molecular pathways controlling PER2 expression during sleep deprivation may be required.

2) Although a detailed description of the method for sleep deprivation is important, I was not able to find it in the main text. The authors should clarify whether sleep deprivation itself or other factors such as stress signals triggered by sleep deprivation affect PER2 expression. At this time, I feel not sleep deprivation per se but some stress signals affect PER2 expression. Many previous studies have shown that glucocorticoids, well-known stress hormones, strongly affect peripheral clocks. Sleep deprivation, one of the stimuli activating stress signals, may therefore affect PER2 expression in peripheral clocks via glucocorticoid signaling pathways. The contribution of glucocorticoids to PER2 expression during sleep deprivation can be examined simply through surgical removal of the adrenal glands.

3) Their conclusions may be specific to the kidney. All bioluminescence data were obtained through monitoring only the kidney. The authors should justify why they can conclude that sleep deprivation affects "peripheral clocks" only judging from kidney data. To generalize kidney data to other peripheral clocks, they should repeat similar experiments using at least another organ or tissue. Real-time monitoring of bioluminescence from the liver is technically simple because injected adenoviruses (carrying a reporter gene driven by a clock gene promoter) specifically infect the liver.

[Editors’ note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "The sleep-wake distribution contributes to the peripheral rhythms in PERIOD-2" for further consideration by eLife. Your revised article has been evaluated by Kate Wassum (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

As you will be able to see in the accompanying comments, while reviewers valued the changes included in this new version of the MS, they still raised important questions. Some may be addressed by reanalyzinging the data already present in the paper (non-circadian reporter mice, CAG-luc). It is claimed that PER2LUC luminescence spikes in Figure 1B are caused by changes in sleep-wake status (analysis is shown in Figure 1C). Currently, for the CAG-LUC animals, there is only group average data (Supplementary Figure 3B). A valid suggestion, hinted by reviewers, is to show individual records: if spikes are not seen in these animals it helps excluding the possibility that the in vivo luminescence spikes are an artifact (e.g. caused by core body temperature changes). Although it doesn't seem to be recorded EEG/EMG from the non-circadian reporter mice, a similar analysis (based on defining sleep as sustained immobility) could be conducted.

Yet, what would really provide strong support to the claims, eliminating doubts about potential artifacts, is to validate the data with conventional methods (immunoblot). As it is mentioned in the MS, during previous work (Curie 2015, PMID: 25581923) the authors have partially validated by Western blot analyses peripheral Per2luc bioluminescence data (Figure S4, liver). Yet, the temporal resolution, and experimental condition of such analyses may not be fully comparable to what is present here,

It is also important to consider that most of the conclusions drawn from the data are based on the temporal correlation between kidney-derived bioluminescence and sleep-wake states. Therefore, it is difficult to generalize the temporal correlation between single organ-derived bioluminescence and sleep-wake states to other peripheral and central clocks. Such limitation of the work (in the absence of additional data) should be clearly acknowledged.

Reviewer #1:

In this paper entitled "The sleep-wake distribution contributes to the peripheral rhythms in PERIOD-2", the authors attempted to find the factor(s) directly contributing to clock gene expression in peripheral tissues. To achieve that, they observed PERIOD-2 (PER2) expression in peripheral tissues of freely moving mice and developed new mathematical models to analyze them. They found that the sleep-wake distribution correlates to the peripheral PER2 dynamics, and that the peripheral PER2 rhythms can be predicted by sleep-wake-state- and peripheral-circadian-dependent factors. Despite the difficulty to dissect the effects of circadian and the sleep-wake cycle on gene expression in vivo, they provided direct evidence that the sleep-wake cycle itself affects peripheral daily rhythms by the unique experimental designs. However, there are still a few points to be addressed before publication including a detailed validation of the PER2 bioluminescence system, which would further support their conclusion.

All things considered, their findings will have an impact on the chronobiology field, as giving direct evidence on sleep-wake-state-dependent regulation of peripheral clock gene expression. In addition, those methods can be useful to investigate the gene expression not only in the kidney or cerebral cortex but also in other various tissues of freely moving mice.

In the revised paper, the authors addressed almost all our concerns with proper correction and additional analyses, but it has some points to be improved before publication, which are summarized in the comments below.

1. The authors assumed that PER2 bioluminescence directly reflected expression of PER2 protein, but it has not been fully confirmed yet especially in the fast (within an hour) dynamics. For further progress of the method, the combination with other techniques to quantify fast and small (< 10%) protein changes would be required, as the authors mentioned in the paper.

2. This reviewer suggests that the authors should discuss about the controversial results regarding PER2-bioluminescene amplitude and sleep-wake distribution in Figure 3 and Suppl. Figure 7. When the amplitude of the sleep-wake distribution was reduced, PER2-bioluminescence amplitude was reduced in 4/6 mice as the model expected, while it was increased in 2/6 mice.

3. Regarding Figure 3E, the authors showed that the 2hOnOff protocol decreased the amplitude of sleep-wake distribution. However, the sleep-wake distribution during the recovery phase was not assessed. Thus, the reviewer recommends that the authors should add data of recovery day to Figure 3E and compare it to the PER2 bioluminescence trends (Figure 3C).

4. In lines 290 to 294, the authors concluded that the prediction from the model, the behavioral rhythm reduction reduces PER2 bioluminescence amplitude, which is consistent with their observation in the Figure 3. However, this conclusion could be difficult to understand because the 2hOnOff protocol increases the PER2 bioluminescence amplitude in one-third of mice (there may be no significant change between BSL and the 2hOnOff protocol; see Figure 3C). The authors provided an additional model (Suppl. Figure 7A-C) to explain the individual difference, and thus the reviewer suggests the authors to describe that conclusion after the modeling part to make the logic clearer.

5. Regarding Suppl. Figure 7A-C, the authors described the increase in PER2 amplitude during the 2hOnOff experiment by advancing the phase of FPERI more than 2.5h. The reviewer suggests that the authors could mention these results in Results or Discussion section.

Reviewer #2:

Hoekstra et al. investigated the relationship between sleep state and PER2::LUC luminescence in freely behaving mice and showed that luminescence acutely changed at sleep/wake transitions. The authors also showed that daily sleep deprivation enhanced the PER2::LUC luminescence rhythm amplitude in SCN ablated mice.

A strength of this study is the simultaneous recording of EEG/EMG/activity and PER2::LUC luminescence in freely moving mice.

Weakness: One obvious concern is that acute PER2::LUC luminescence changes occurring at sleep/wake transitions are artifacts (e.g. due to acute change in core body temperature), but not caused by synthesis/degradation of PER2 protein. Because those luminescence changes are relatively large and fast (30-40% change within 30 min), a conventional method may be used for validation. The authors recorded activity and bioluminescence from non-circadian luminescence reporter mice (CAG -luc). Alternatively, the authors can analyze that data and exclude artifact possibility if luminescence change doesn't occur at active/inactive transitions.

Another weakness is the PER2::LUC luminescence measurement from the brain. The authors must show IVIS imaging and show the signal coming from the brain and not from other tissues.

Abstract: Line 18, 19: instead of saying PER2 increases or PER2 rhythmicity, I suggest using PER2::LUC luminescence increases, PER2::LUC luminescence rhythmicity.

Line 89, Suppl. Figure 1: should be Suppl. Figure 3.

Figure 1B: The hypnogram only includes sleeping and waking, please include REM and NREM. I understand that the authors only wanted to indicate sleep/wake transition used for the luminescence analysis indicated in C, but I think the authors should plot standard hypnogram.

Figure 1C: To avoid confusion, it is better to use the same colors as in A and B (pink: peripheral, purple: central).

Figure 2: The authors should show individual records in the supplemental information. I know the authors discussed two possibilities that can explain changes in amplitude in group average data. The amplitude of an individual mouse was shown in Figure 2C. But I think, like me, the readers will still want to see the individual record. The amplitude of 3 cycles in REC should be analyzed separately (Plot as REC day 1, REC day 2 and REC day 3). Technically REC day 1 is the same as SD, but it is good to show how the amplitude in the group average and the individual mouse declines.

Line 573-584: This is an important discussion addressing the possible temperature artifact. However, all statements are based on analysis, "data not shown", unfortunately. The authors must show the data.

Line 763: I think the activity data was binned, not averaged.

Supplementary Figure 1: I know luminescence recording was done in DD, but please indicate this, as the lighting condition in the other treatments was stated.

Supplementary Figure 3B: Please indicate individual records as well, so the relationship between luminescence and activity can be evaluated.

Reviewer #3:

Response to the comment #1

It is my understanding that almost all of the conclusions made by the authors are based on the temporal correlation between kidney-derived bioluminescence and sleep-wake states. Therefore, I do not feel convinced about these conclusions without additional experimental evidence suggesting the causal link between spontaneous wakefulness and PER2 expression.

A few previous studies, although not cited in the present version, clearly demonstrated that arousal levels affect the circadian clock via the SCN (for example, PNAS 2016 Yamakawa el al). I think that the present data may not make obvious progress from these previous achievements.

Response to the comment #3

The authors claim that spontaneous wakefulness-driven PER2 increase is one of the major findings. Therefore, under not sleep deprivation but spontaneous wakefulness conditions, they need to confirm PER2 expression in other organs and apply expression data to their mathematical model.

Many previous studies have indicated that circadian phase in clock gene expression varies among peripheral organs and tissues, and it is therefore difficult to generalize the temporal correlation between single organ-derived bioluminescence and sleep-wake states to other peripheral and central clocks. In particular, circadian phase of the SCN is largely phase-advanced in comparison with that of peripheral clocks, raising a question on the author's thought that spontaneous wakefulness is a similar trigger of PER2 in both peripheral and central clocks.

[Editors’ note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "The sleep-wake distribution contributes to the peripheral rhythms in PERIOD-2" for further consideration by eLife. Your revised article has been evaluated by Kate Wassum (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

There are only two final pending issues that we anticipate should be easy to incorporate as part of the results/discussion and be included in the final rebuttal letter, without the need for any additional experiments.

Essential revisions:

1) Although the new data has clearly enriched the paper, this still raises the concern of one of the reviewers regarding the presence of low amplitude baseline rhythm in luminescence (particularly as CAG-luc appears to be activated by wake: the reviewer expects to see a basal circadian change of luminescence in CAG-luc mice).

2) Without showing an absolute luminescence value emitted from the brain, liver, and kidney in the same mouse, it cannot be utterly said that the luminescence signal from the entire mouse is just coming from the cortex. Part of this concern is based on the fact that the Curie et al., 2015 work not only shows signals (in B6 mice) coming from the cortex, but also clear signals coming from the kidney/liver (as visualized with the IVIS signals, whereas the RT-Biolumicorder does not allow sorting those signals).

Reviewer #2:

Hoekstra et al. investigated the relationship between sleep state and PER2::LUC luminescence in freely behaving mice and showed that luminescence acutely changed at sleep/wake transitions. The authors also showed that repeated daily sleep deprivation reinitiated the PER2::LUC luminescence rhythm in SCN ablated behaviorally arrhythmic mice. This re-initiated rhythm sustained a few cycles in constant conditions, indicating that sleep can synchronize (or coordinate) peripheral circadian oscillators.

A strength of this study is the simultaneous recording of EEG/EMG/activity and PER2::LUC luminescence in freely moving mice.

One obvious concern is that acute PER2::LUC luminescence changes occurring at sleep/wake transitions are artifacts (e.g. due to acute change in core body temperature, Calcium/Magnesium, pH, ATP), but not caused by the synthesis/degradation of PER2 protein. After two major revisions, the authors now provided the data showing that no luminescence change occurs at sleep/wake transitions in Pkg1-luc activity (housekeeping gene promoter activity measured by luminescence reporter). Although it is still difficult for this reviewer to believe 30-40% changes in PER2 protein can occur within 30 min in the brain and peripheral organs, the new Pkg1-luc data rejects the above-mentioned concern. The authors also provided an interesting new analysis that acute luminescence changes at sleep/wake transitions in CAG-luc mice are much faster compared to changes in PER2::LUC mice. The authors reasoned this by saying the CAG synthetic promoter is activated by transcription factors (e.g. NF-κB, CREB, AP-1). This is an interesting result but raises a new concern. The authors used this CAG-luc mouse to optimize the luciferin administration method, as this is a non-circadian reporter. However, if being awake activates CAG promoter, this reviewer expects to see a basal circadian change of luminescence in CAG-luc mice.

Another relatively minor concern is that the authors used PER2::LUC in pigmented C57BL/6 mice with a cranial window for "the central quantification of bioluminescence". This reviewer suggested the authors to perform IVIS imaging and show the signal coming from the brain and not from other tissues. Instead, the authors referred to their previous publication (Curie et al., 2015) that the luminescence signal is from the cortex. However, in this paper, the authors actually detected a signal from the liver and kidney as well. Because the authors only showed the hold changes in those tissues (brain, liver kidney), it is not known how strong the signal from cranial windows compare to that from the liver and kidney. Without showing an absolute luminescence value emitted from the brain, liver, and kidney in the same mouse, the authors cannot say that the luminescence signal from the entire mouse is from the cortex.

This reviewer appreciates the efforts of the authors to address the concerns raised by the reviewer.

This reviewer can now recommend this work to be published in this journal. However, those two concerns need to be addressed.

One additional suggestion is that it may not be clear to the reader that mice with a cranial window is used for measuring signal from the cortex. This is only mentioned in the methods section. It may be better to state this in the main text and figure caption.

[Editors’ note: the authors resubmitted a revised version of the paper for consideration. What follows is the authors’ response to the first round of review.]

Reviewer #1:

In this paper entitled "The sleep-wake distribution contributes to the peripheral rhythms in PERIOD-2", the authors observed PER2 expression in peripheral tissues of freely behaving animals and found that the sleep-wake distribution contributes to the peripheral PER2 dynamics. In addition, implementing the mathematical model, they showed that the peripheral PER2 rhythms can be generated by two factors: a sleep-wake state dependent one and a peripheral (not master) circadian dependent one. Despite the difficulty to dissect the effect of these two factors in vivo, the authors, dealing with this problem by the unique experimental designs, provided direct evidence that sleep-wake state itself affects peripheral daily rhythms. However, this paper has a couple of points to be improved, which are summarized in the comments below.

We thank the reviewer for his positive evaluation and particularly for stating that the results “provided direct evidence that sleep-wake state itself affects peripheral daily rhythms”, which is the central claim of our work. We also thank the reviewer for the suggestions to help improve our manuscript.

1. As for Figure 1B, 1C and 1D, the authors showed PER2 bioluminescence changes according to circadian rhythm and state transitions in only one individual mouse (ID: 2MP82), which might be insufficient for concluding that the PER2 bioluminescence changes depending on the circadian and sleep-wake cycle. Thus, the reviewer suggests that the authors should add data about other mice and perform the statistical analyses to clarify the validity of the circadian- and sleep-wake-dependent change of PER2 bioluminescence. With Figure 4, the reviewer suggests the same, i.e., the authors should analyze not only one representative mouse but a group of mice.

Indeed, Figure 1B and -C (now Figure 1B) represented an example of one mouse. Group means based on 5 mice were, however, depicted in Figure 1A and -D (now Figure 1A and C), as was mentioned in the figure legend, illustrated by the individual recordings of the 4 remaining mice in Suppl. Figure 3 (now Suppl. Figure 4), stated in the ‘transparent reporting’ form in our submission, and can be deduced from the raw data we provided for this and all other figures. The same holds true for Figure 4 where the model was based on recordings in 4 and 6 mice, respectively. In response to a suggestion from Reviewer 2, we have included a supplementary figure (Suppl. Figure 1) with an overview of all experiments explicitly stating the number of animals used in each.

Plus, there might be a mistake in the time label of Figure 1B (is this plot started from time = 6 or time = 0?).

The plot starts at 12h, at a time of day that under the previous LD cycles is referred to as ZT12; i.e., the beginning of the dark phase. The tick and label indeed disappeared. This has been corrected in the revised version.

2. In lines 202 to 205, the authors concluded that the behavioral rhythm reduction by 2hOnOff protocol appeared to be insufficient to consistently reduce the amplitude of the PER2-bioluminescence oscillation. However, it is difficult to follow its logic toward the conclusion. Thus, the reviewer recommends that the authors should describe their logic in more detail.

We thank the reviewer for this suggestion. With the use of the word ‘appeared’ we hoped it was clear that this was an interpretation. Our logic was as follows: If the sleep-wake distribution contributes to the diurnal changes in PER2 bioluminescence then we can expect that after experimentally reducing the amplitude of the sleep-wake distribution, the amplitude of the PER2 oscillation will be reduced as a result. With the 2hOnOff protocol we reduced, but did not eliminate, the amplitude of the sleep-wake distribution. This reduction was however accompanied by a reduction of PER2 amplitude in only 4 of the 6 mice and we therefore concluded that the reduction in sleep-wake amplitude appeared to be insufficient to elicit the anticipated PER2 effects in all mice. We now have eliminated this interpretative sentence from the Result section.

3. As for Figure 4 and Table 2, the authors compared two mathematical models and concluded that the two forces model has better prediction accuracy. However, it is unclear whether the circadian force itself or the combination of the two forces is important to predict PER2-bioluminescence dynamics. The reviewer suggests that the authors should add a comparison with a model that using a single circadian force (F = FPERI).

Again, we thank the reviewer for this valuable suggestion. We have now performed the suggested additional analysis in the revised manuscript, illustrating better that both forces are needed and their combination represents an important statistical gain in predicting PER2. The results have been added to Figure 4B and C and are presented in the text (l. 326-332). We now extended the modeling to predict results of the 6h SD experiments using the parameters obtained in the other two experiments (Figure 4D, Suppl. Figure 7D, Table 2, l. 297-302).

4. In lines 277 to 280, the authors described the increase in PER2 amplitude during the 2hOnOff experiment by advancing the phase of FPERI by 6h. The reviewer suggests that the authors should clearly explain why they apply their mathematical model to the 6h advancing, although they observed 1-2h phase advance in their experiments.

This analysis was meant as thought-experiment to illustrate on the one hand the power of our dynamic model to perform in silico experiments and, on the other hand, that the initial conditions could give rise to diametrically opposing results. We made this clearer in the manuscript. We have now regrouped this analysis together with the analysis mentioned under the preceding point 3, as one of three questions addressed by the model as in silico experiments (l. 356-379). In a re-analysis for the revised manuscript, we found that the same effect (i.e., an increase of PER2 bioluminescence amplitude instead of the expected decrease) could already be observed with a 2.5h phase advance. We have added this finding to the revised manuscript (l. 360-365) and illustrated this analysis in Suppl. Figure 7C.

Reviewer #2:

By measuring EEG/EMG and PER2::LUC bioluminescence (most signal comes from the kidney) from freely behaving mice simultaneously or from separate cohort mice with similar conditions, the authors observed that sleep state acutely changed PER2::LUC bioluminescence, repeated daily sleep deprivation reinstated circadian PER2::LUC rhythm in SCN ablated behaviorally arrhythmic mice. The authors also did a mathematical simulation modeling peripheral PER2::LUC bioluminescence influenced by a sleep dependent force and a SCN independent force.

I believe the data presented in the current manuscript, however, I don't think the main conclusions are supported by the data. The acute changes (within 30 min) of luminescence they observed is most likely caused by in vivo factors (Not changed in PER2::LUC protein, artifact). It is very problematic to assess the amplitude change from averaged luminescence record from multiple behaviorally arrhythmic mice (analysis flaw). Although the above mentioned two major conclusions must be removed from the manuscript, the study validates previously published results done by ex vivo tissue and in vivo anesthetized animal with luminescence measure from freely behaving mice. So I encourage the authors to publish the data to another journal (please also see my specific points below), but I don't think discovery reported in the work is novel for the publication in eLife.

We have addressed each of the issues mentioned in detail below. What the reviewer refers to as “very problematic” and “analysis flaw” could be easily addressed. Moreover, the reasons given to deem our bioluminescence signal as artifact are flawed.

1. It is unlikely that the acute change of PER2::LUC luminescence reflects PER2 protein translation and degradation. As the authors refer to correctly, luminescence from PER2::LUC knockin mice reflects PER2::LUC fusion protein (not Per2 transcription). The authors didn't correctly refer to the studies measuring light induction of Per2 mRNA and protein in the SCN. The light induction of Per2 mRNA in the SCN can only be seen 60 min after onset of the light pulse (no induction was observed at 30 min after the light pulse. Takumi et al. 1998 and many other publications). Increased PER2 protein in the SCN can only be seen at 4 h after the light pulse (no change was observed 2 h after the light pulse. Yan and Sliver 2004, Antoun et al., 2012). Therefore, the acute luminescence change observed at sleep wake transitions (Figure 1D) is unlikely reflecting changes in the PER2::LUC fusion protein. Temperature, pH, and Ca/Mg and ATP levels all strongly affect luminescence produced by firefly luciferase and luciferin. It was shown that the ultradian rhythm of the core body temperature closely matched the ultradian behavior rhythm (Blessing 2018 PMID: 30454601). The authors referred to their own published study that luminescence rhythm change was vilified by qPCR. However, the time scale of their published work showed that changes occurred after 6 h of sleep deprivation, not in the 30 min range.

First, we have to correct two of the reviewer’s statements: (1) We did not claim that the literature showed that 30min light-pulses induce PER2 in the SCN and (2) we previously used western blot, and not qPCR, to validate the sleep-deprivation induced increases in PER2 bioluminescence in brain and liver.

We are aware that no study has reported on such quick induction of PER2 after a light pulse. We did, however, not study light effects in the SCN, and the effects of wakefulness on PER2 in peripheral tissues are vastly under-explored and might differ from those of light in the SCN as a variety of molecular pathways activated by waking are known to affect PER2 (reviewed in the Discussion). The fact that in the literature there is no precedent of wakefulness rapidly inducing small but significant changes in PER2 using more conventional techniques, cannot be used as argument that we reported on a mere artefact. We already cited the study by Cao et al. (Nat Neurosci 2015), demonstrating that PER2 protein can indeed be rapidly induced through other means than light; i.e., within 30 min after a serum shock (see Discussion paragraph l. 400-414). The wake-related rapid induction of PER2 is an important novel aspect of our study. Whether the signal obtained with the PER2::Luciferase construct indeed reflects PER2 protein was already discussed in detail in the Discussion (“Do changes in bioluminescence reflect changes in PER2 levels?”).

Although the variables listed by the reviewer do affect the bioluminescence reaction, they seem not to play a large role (see Discussion section: ‘Do changes in bioluminescence reflect changes in PER2 levels?’). The reviewer’s example of a potential temperature effect can illustrate this. Body temperature changes in a circadian fashion, yet the bioluminescence time courses reporting on clock genes known to be regulated in anti-phase, remain anti-phasic and do not simply follow the temperature dynamics. Also the lack of a noticeable rhythm in the CAGbioluminescence signal (Suppl. Figure 3B) confirms that the circadian rhythm in activity and associated temperature changes have little or no effect (or even decreases when activity is highest, see our response to point 6 below). We have added this observation to the Discussion (l. 577-579). While temperature increases quickly (within minutes) after awakening and stays high while awake (see the Blessing et al. reference mentioned by the reviewer and our own work: Hoekstra et al. eLife 2019), PER2-bioluminescence follows a much slower time course and initially continues to decrease before subsequently increasing (see Figure 1C). For this rebuttal we followed both central (see also comment 3 of Reviewer 3) and peripheral PER2-bioluminescence for as long as waking prevailed as opposed to the 1st 30min after wake onset only. The new analysis shows that peripheral PER2 increases and then again decreases with increasing wake duration, while central levels monotonically keep increasing consistent with the western-blot verified tissue-specific effects of experimentally prolonging wakefulness (Figure 1C and brain vs. liver in Curie et al. 2015). These results further underscore that changes in PER2 bioluminescence do not follow wake related changes in temperature. We have replaced the transition analyses in Figure 1C with this new analysis.

2. The advantage of continuous luminescence measurement from the same mouse is that the presence and absence of the rhythm can be determined from a mouse that has no phase marker (e.g. SCN lesioned mice). Arrhythmicity (or lowered amplitude of the rhythm) observed by conventional population sampling (from multiple mice) cannot be determined if individual mice lose rhythmicity, or if individual mice retained intact rhythm but they are out of sync. From the data presented in Figure 2, the authors claimed repeated sleep deprivations reinstated the rhythm, but they didn't analyze amplitude change before, during, and after sleep deprivation in individual mice (even they didn't provide individual data anywhere else, they provided individual data in other experiments).

The authors stated that they performed daily sleep deprivation to the mice 5 weeks after SCN-lesioned. From average luminescence rhythms of 4 SCN-X mice, the authors cannot conclude their main claim that sleep deprivation reinstated the peripheral circadian rhythm. It is possible there is no change in the amplitude of individual peripheral circadian oscillators, but repeated sleep deprivation entrained those rhythms (the authors averaged individual luminescence relative to the phase of daily sleep deprivation).

We thank the reviewer for raising this issue. We had provided the individual amplitude information for the baseline bioluminescence rhythms of the SCNx mice in the first version of our manuscript (Suppl. Figure 5C and Figure 2B), but indeed did not include this information during and after the sleep deprivations. As mentioned in our reply to minor issue 1 from Reviewer 1, we have added these data to main Figure 2C. The results clarified that increased amplitude during the sleep deprivations does not simply result from a phase alignment as it is significantly larger than in the baseline recordings. This analysis and conclusion has now been added in the Results (l. 206-210).

It is well known that the SCN-lesioned mice do not show a daily rhythm under LD cycle (Light cannot reinstate behavior rhythm).

We did not make this claim in the manuscript and therefore fail to see the relevance of this comment. Moreover, the literature indicates that this claim is less well-known as the reviewer makes us believe (specifically addressed in Redlin and Mrosovsky J Comp Physiol A 1999, reviewed by Redlin Chronobiol Int 2001). Light is an important driver of rhythmic sleep-wake behavior (see e.g. Hubbard et al. PNAS 2021 in press; bioRxiv).

3. In the complementary experiment presented in Figure 3, the authors analyzed the amplitude changes from individual mice, but didn't analyze from the averaged data. The authors observed a reduced amplitude by 2hOnOff procedure in 4 out of 6 mice, and an increased amplitude in 2 out of 6 mice. But the authors concluded that the 2hOnOff procedure did not "consistently" affect the rhythm (in another place, they wrote "2hOnOff protocol didn't noticeably affect the ongoing circadian PER2 oscillation").

We did analyze the amplitude changes at the group level by comparing PER2 amplitude under baseline with the 2hOnOff regime, and mentioned the results in the figure legend (“resulting in overall lack of an effect of the 2hOnOff protocol”) along with the statistics. In the Results section, the lack of an overall effect was reiterated twice. We now have added the mean values and the SEMs to the figure legend for additional clarity.

4. Because of the issues pointed out above, the parameters used for the mathematical simulation are not valid.

The reviewer claims the parameters for the simulation are not valid, but does not explain how he or she arrives at this conclusion. Many issues were raised in the preceding three points, none of which make reference to the model. We will not try and guess what the reviewer might have meant especially since we have now addressed the issues pointed out above.

Use of locomotor activity as a wake parameter is also not appropriate. Mice can be wake without moving.

Indeed, locomotor activity can importantly deviate from time-spent-awake in certain situations. This was however not the case in the current experiments as could be deduced from the very strong positive correlation between waking and LMA (R²=0.79) which we presented in Suppl. Figure 5A (now Suppl. Figure 6B), confirming that the two variables strongly covary. Moreover, from a modeling perspective the important point is not whether locomotor activity one-to-one relates to waking but that the time courses of the two variables match. We have added another panel to this supplementary figure illustrating the time course of LMA and waking to underscore their similarity (Suppl. Figure 6A). Finally, we now have modelled PER2 dynamics in the 6h SD experiments based on either wakefulness or LMA. Both resulted in similar tight fits. These findings were added to the results and Suppl. Figure 6C. To conclude, in the context of the model, locomotor activity is a highly appropriate proxy for the sleep-wake distribution.

5. In the method section, the authors described each procedure. However, it was impossible for me to know how exactly each individual mouse went through the procedures. It seems some sleep deprivation was done under dim red light and other was done under light. I suggest the authors provide a timeline (e.g. flowchart figure) of procedures in each experiment.

We thank the reviewer for this useful suggestion and we have included the suggested figure as Suppl. Figure 1.

6. Since in vivo bioluminescence recording became popular technology, I strongly suggest the authors move the data presented in supplementary Figure 2 to the main text as an actual figure. Were CAG mice also crossed with hairless strain? Does tissue explanted from this mouse in vitro not show a circadian rhythm? It even looks like luminescence from the mice with luciferin supplemented by osmic pump shows a weak rhythm. I wonder if CAG promoter may have circadian rhythm as CMV promoter does (Collaco et al., 2005). Amplitude of activity rhythm is lowered in mice with pump. Is this due to surgery (pump implantation)? If so, lowered amplitude of luminescence rhythm might also be affected by surgery. Having a control (saline in pump and luciferin supplied through drinking water) is necessary.

Although the many suggested follow-up experiments related to the results obtained in the CAGmice have their merit, addressing them experimentally is clearly far removed from the aims of our manuscript and would be more appropriate for a separate, technical paper. With a straightforward and unique approach, we show that administering the luciferase substrate in the drinking water gives rise to high-amplitude bioluminescence rhythms that all but disappear when using mini-pumps. The by far most plausible explanation is that the large-amplitude rhythms are an artefact associated with the circadian timing of water intake (and thus luciferin availability) and that therefore mini-pumps are to be preferred to administer the substrate.

Indeed, a small yet significant bioluminescence rhythm is present in CAG mice implanted with mini-pumps which does, however, not challenge our choice of route of administration. Moreover, this residual bioluminescence rhythm peaks 6.3 circadian hours prior to when the activity rhythm peaks and bioluminescence thus decreases when animals are awake most, again demonstrating that changes in bioluminescence do not simply follow changes in e.g. body temperature associated with the circadian sleep-wake distribution.

The reviewer claims that the amplitude of the activity rhythm in mice implanted with minipumps is lower. Statistics fail to support the reviewer’s claim (P=0.67, t=0.47, paired t-test, n=4 on amplitudes derived from sine-wave fits of individual time-course data).

Reviewer #3:

This is an interesting study. As far as I know, there are no previous studies performing simultaneous measurement of sleep-wake state and bioluminescence (clock gene expression) in freely moving mice. In addition, their measurement method was well validated through several pilot studies, and measurements obtained using this method are therefore reliable. In fact, the presented data are clear and convincing. However, although I have no doubt that the results are useful and valuable for a wide range of circadian scientists, I feel the present version of their work may not be a strong candidate for publication in this journal because of insufficient novelty and impact.

We thank the reviewer for judging our work well validated and our data reliable, clear, and convincing. Our study is indeed the first to simultaneously record EEG/EMG and bioluminescence, as well as locomotor activity allowing to show that undisturbed sleep-wake behavior contributed to PER2 dynamics.

1) The results clearly revealed the temporal relationship between sleep-wake state and PER2 bioluminescence through simultaneous monitoring in freely moving mice, with high time resolution. However, because previous studies have already found that sleep deprivation results in a change in PER2 levels, I feel that the present study does not make obvious progress toward further understanding the mechanism. The simultaneous monitoring experiments are interesting but each of the techniques was originally established by other groups. For example, some findings related to the identification of systemic or molecular pathways controlling PER2 expression during sleep deprivation may be required.

The reviewer states our work lacks the novelty that would justify publication in eLife by referring to previous sleep-deprivation studies. In doing so, the reviewer overlooked an important novel aspect of our current study: we are the first to have quantified the relationship between spontaneous sleep-wake behavior and PER2-bioluminescence, whereas our previous work only assessed the effects of sleep deprivation. We discovered that like sleep deprivation, also spontaneous wakefulness, raises PER2 levels. The finding that PER2’s increase is associated with wakefulness per se and does not result from confounding physiological changes accompanying the sleep deprivation (e.g. stress, see next issue), is a core assumption of the mathematical model which, with high accuracy, could predict PER2’s behavior under three challenging experimental protocols. Thus, the circadian distribution of sleep-wake behavior under undisturbed conditions, considered by most a mere output of the clock, modulates the circadian time-keeping circuitry and is required to maintain high-amplitude PER2 oscillations. We believe this to be a major advance in understanding circadian-rhythm generation in tissues peripheral to the SCN as it challenges current dogma of the clockgene circuitry being responsible for all rhythmic gene expression, including clock genes themselves. In the Discussion we had already reviewed the molecular pathways by which sleep-wake state could affect the clock-gene circuitry. We recently published on one of these pathways (temperature – CIRBP – clock genes; Hoekstra et al. eLife 2019) through which the sleep-wake distribution could modulate clock-gene expression outside the SCN. We now have put more emphasis on the importance of spontaneous sleep-wake behavior on PER2 bioluminescence rhythms and added an analysis documenting this further (l. 344-353).

2) Although a detailed description of the method for sleep deprivation is important, I was not able to find it in the main text.

We indeed did not describe the sleep-deprivation method in detail but referred to our sleep- phenotyping methods paper instead [“…mice were sleep deprived for six hours at a time they were expected to rest (….) by gentle handling as described (Mang and Franken, 2012)”]. We have expanded on this method in the revised manuscript (l. 683-689).

The authors should clarify whether sleep deprivation itself or other factors such as stress signals triggered by sleep deprivation affect PER2 expression. At this time, I feel not sleep deprivation per se but some stress signals affect PER2 expression. Many previous studies have shown that glucocorticoids, well-known stress hormones, strongly affect peripheral clocks. Sleep deprivation, one of the stimuli activating stress signals, may therefore affect PER2 expression in peripheral clocks via glucocorticoid signaling pathways. The contribution of glucocorticoids to PER2 expression during sleep deprivation can be examined simply through surgical removal of the adrenal glands.

Short-term sleep deprivation by gentle handling is indeed considered a mild stressor and corticosterone does increase as we had already clarified in the Discussion (“How does wakefulness increase PER2?”). We published the experiment the reviewer suggested a while ago (Mongrain et al. Sleep 2010). This study shows that the increase in Per2 expression during sleep deprivation is not solely due to the sleep-deprivation associated corticosterone increase as in adrenalectomized mice sleep deprivation still significantly increases Per2 levels. These results, together with the results obtained after spontaneous periods of waking (which are not accompanied by an increase in corticosterone) presented in the current study, strongly argue against the simple scenario proposed by the reviewer. These issues are addressed in the Discussion (l. 499-511).

3) Their conclusions may be specific to the kidney. All bioluminescence data were obtained through monitoring only the kidney. The authors should justify why they can conclude that sleep deprivation affects "peripheral clocks" only judging from kidney data. To generalize kidney data to other peripheral clocks, they should repeat similar experiments using at least another organ or tissue. Real-time monitoring of bioluminescence from the liver is technically simple because injected adenoviruses (carrying a reporter gene driven by a clock gene promoter) specifically infect the liver.

In an earlier publication (Curie et al. Sleep 2015), using an in vivo imaging approach, we already demonstrated that sleep deprivation increases PER2 in cortex, liver, and kidney, which, for cortex and liver, were confirmed by western blot. The baseline and recovery bioluminescence dynamics obtained in central and peripheral tissue with the RT-Biolumicorder depicted in Figure 1A match remarkably well those reported for cortex and liver and kidney, respectively, with the in vivo imaging method. We are therefore confident that the effects we report can be generalized to other organs peripheral to the SCN.

When setting up the experiments, we also recorded central (cortical) PER2bioluminescence together with EEG/EMG in 6 C57BL/6J PER2::LUC mice. Because detected bioluminescence levels were somewhat noisier (l. 167-169), for the original manuscript we decided to focus on the better accessible (hairless SKH1 mice) and more robust peripheral (kidney) signals. Nevertheless, also the cortical bioluminescence signals demonstrate that the results we reported are not specific to kidney, albeit with tissue-specific dynamics during spontaneous wake episodes. We have included these central data in Figure 1 and addressed the tissue-specificity in the Discussion (l. 531-551).

[Editors’ note: what follows is the authors’ response to the second round of review.]

Reviewer #1:

In this paper entitled "The sleep-wake distribution contributes to the peripheral rhythms in PERIOD-2", the authors attempted to find the factor(s) directly contributing to clock gene expression in peripheral tissues. To achieve that, they observed PERIOD-2 (PER2) expression in peripheral tissues of freely moving mice and developed new mathematical models to analyze them. They found that the sleep-wake distribution correlates to the peripheral PER2 dynamics, and that the peripheral PER2 rhythms can be predicted by sleep-wake-state- and peripheral-circadian-dependent factors. Despite the difficulty to dissect the effects of circadian and the sleep-wake cycle on gene expression in vivo, they provided direct evidence that the sleep-wake cycle itself affects peripheral daily rhythms by the unique experimental designs. However, there are still a few points to be addressed before publication including a detailed validation of the PER2 bioluminescence system, which would further support their conclusion.

All things considered, their findings will have an impact on the chronobiology field, as giving direct evidence on sleep-wake-state-dependent regulation of peripheral clock gene expression. In addition, those methods can be useful to investigate the gene expression not only in the kidney or cerebral cortex but also in other various tissues of freely moving mice.

In the revised paper, the authors addressed almost all our concerns with proper correction and additional analyses, but it has some points to be improved before publication, which are summarized in the comments below.

1. The authors assumed that PER2 bioluminescence directly reflected expression of PER2 protein, but it has not been fully confirmed yet especially in the fast (within an hour) dynamics. For further progress of the method, the combination with other techniques to quantify fast and small (< 10%) protein changes would be required, as the authors mentioned in the paper.

We fully agree with the reviewer and if techniques were available that could reliable quantify the small but significant changes in protein levels that we reported on, we would have readily used them. However, for several reasons this is less straightforward as it might seem: (i) Western blots are not sensitive enough to pick up such small differences as this method is semi-quantitative at best and (ii) the design can no longer be longitudinally (i.e. within the same individual) in contrast to our current in vivo technology. This latter point is important when assessing the effects of spontaneous wakefulness on PER2 as this would require each animal to be closely monitored (using EEG/EMG signals) for periods of either sleep or wake of similar duration to appear at approximately the same circadian phase, before rapidly killing the mouse for tissue extraction such that waking up (if asleep) and stressing the animal does not affect the measurement. Moreover, because of the semi-quantitative nature of Western blots, small effect size, and group instead of within-subject comparisons, we would need to kill a large number of mice for quantifying PER2 for each time point and a time course analyses, as presented here, would be practically impossible. For these reasons, we decided to analyze bioluminescence dynamics in mice with other luciferase reporter constructs to rule out nonspecific sleep-wake related effects on the bioluminescence reaction. See our response to the first comment by Reviewer 2 below.

2. This reviewer suggests that the authors should discuss about the controversial results regarding PER2-bioluminescene amplitude and sleep-wake distribution in Figure 3 and Suppl. Figure 7. When the amplitude of the sleep-wake distribution was reduced, PER2-bioluminescence amplitude was reduced in 4/6 mice as the model expected, while it was increased in 2/6 mice.

We believe that, given the data, we went as far as we could in trying to explain these counterintuitive results. Our model gave us a tool to investigate this issue and the simulation results suggest that the phase of the circadian force relative to the phase of the wake force could be a contributing factor. This scenario was presented and discussed in detail in the Results section (Lines 386-410). However, this explanation follows from an in silico analysis and, as we pointed out in that section, it is clear that this hypothesis must be tested by follow-up experiments before making any firmer claims.

3. Regarding Figure 3E, the authors showed that the 2hOnOff protocol decreased the amplitude of sleep-wake distribution. However, the sleep-wake distribution during the recovery phase was not assessed. Thus, the reviewer recommends that the authors should add data of recovery day to Figure 3E and compare it to the PER2 bioluminescence trends (Figure 3C).

We recorded recovery after the 2hOnOff protocol, but did not further analyze the sleep-wake distribution during this period. Following the reviewer’s recommendation, we now have analyzed sleep-wake states for 48h of recovery and added the estimated amplitude of the sleep-wake distribution during these two recovery days to Figure 3E. After the reduction of the daily sleep-wake amplitude during the 48h 2hOnOff experiment, during recovery amplitude was restored to baseline levels.

4. In lines 290 to 294, the authors concluded that the prediction from the model, the behavioral rhythm reduction reduces PER2 bioluminescence amplitude, which is consistent with their observation in the Figure 3. However, this conclusion could be difficult to understand because the 2hOnOff protocol increases the PER2 bioluminescence amplitude in one-third of mice (there may be no significant change between BSL and the 2hOnOff protocol; see Figure 3C). The authors provided an additional model (Suppl. Figure 7A-C) to explain the individual difference, and thus the reviewer suggests the authors to describe that conclusion after the modeling part to make the logic clearer.

5. Regarding Suppl. Figure 7A-C, the authors described the increase in PER2 amplitude during the 2hOnOff experiment by advancing the phase of FPERI more than 2.5h. The reviewer suggests that the authors could mention these results in Results or Discussion section.

At the group level the reduction in PER2 amplitude during the 2hOnOff protocol did not reach significance, as was clearly stated in the Results section and figure legend. Yet, when inspecting the individual responses, we found an interesting pattern suggesting that the 2hOnOff protocol did affect PER2 amplitude albeit in opposite direction for individual mice (Lines 251-261). With the analysis presented in Suppl. Figure 7A-C (now Figure 4 —figure supplement 2A-C), we did not introduce an additional model to explain this phenomenon. Instead, we used the same model as proposed earlier, but instead of letting the model determine a single (optimal) phase of the circadian force, we systematically varied it to predict the resulting effects on PER2 amplitude. We described the rational and results of this in silico experiment in the Results section after the modeling part (Lines 386-410). To enhance clarity, we have now commented after line 256 (now Lines 260-263) on the results of the modelling section, stating that variation in circadian phase could underlie the differential individual response.

Reviewer #2:

Hoekstra et al. investigated the relationship between sleep state and PER2::LUC luminescence in freely behaving mice and showed that luminescence acutely changed at sleep/wake transitions. The authors also showed that daily sleep deprivation enhanced the PER2::LUC luminescence rhythm amplitude in SCN ablated mice.

A strength of this study is the simultaneous recording of EEG/EMG/activity and PER2::LUC luminescence in freely moving mice.

Weakness: One obvious concern is that acute PER2::LUC luminescence changes occurring at sleep/wake transitions are artifacts (e.g. due to acute change in core body temperature), but not caused by synthesis/degradation of PER2 protein. Because those luminescence changes are relatively large and fast (30-40% change within 30 min), a conventional method may be used for validation. The authors recorded activity and bioluminescence from non-circadian luminescence reporter mice (CAG -luc). Alternatively, the authors can analyze that data and exclude artifact possibility if luminescence change doesn't occur at active/inactive transitions.

Thank you for this great suggestion. With the reviewer we would have preferred to use a conventional method for confirmation but as explained in our answer to the first comment of Reviewer 1 this is unfortunately less straightforward as it might seem.

Following your suggestion, we analyzed bioluminescence dynamics at the onset of activity bouts in the CAG-Luc mice. As you can see in and the new Figure 1 —figure supplement 5, activity bouts in CAG-Luc mice are accompanied by robust increases in bioluminescence, even larger in amplitude than those observed in PER2^Luc mice. This result is consistent with the CMV major-immediate early enhancer/promotor, which is part of the CAG synthetic promotor driving luciferase, being rapidly activated by large number of transcription factors such as NF-κB, CREB/ATF, and AP-1 (see review by Stinski and Isomura Med Microbiol Immunol 2008; doi: 10.1007/s00430-007-0069-7), a response comparable to that of an immediate early response such as c-FOS (see e.g. Brightwell et al. Gene 1997; doi: 10.1016/S0378-1119(97)00178-9 and Collaco and Geusz BMC Physiol. 2003; doi: 10.1186/1472-6793-3-8).

We used the CAG-Luc mice to verify the effect of the route of luciferin administration on circadian bioluminescence dynamics because we expected this synthetic promotor not to be regulated in a circadian fashion (which we confirmed). We did, however, not chose this mouse line to rule out that PER2bioluminescence changes at sleep-wake transitions are artefacts, an issue that was raised during the review process. After discussing this issue with Prof. David Gatfield, a colleague at my institute studying circadian aspects of translation, he shared some of his recordings with other luciferase-reporter constructs using the same set-up (i.e., osmotic mini-pumps delivering luciferin and the RT-biolumicorder system recording bioluminescence). We recorded bioluminescence in mice expressing luciferase under the control of the promotor of the housekeeping gene Pgk1 (Pkg1-Luc mice). As expected of a housekeeping gene, overall bioluminescence levels do not consistently vary during the recording and no activity-bout related changes were observed (Figure 1 —figure supplement 5). These recordings demonstrate that the increases in PER2 bioluminescence we observed at sleep-wake transitions are not an artifact of the sleep-wake transition as similar activity-related increases should have been observed for all 3 constructs. We have added the Pgk1-Luc results together with the CAG-Luc results to the Results section (Lines 177-182) and added a new supplementary figure (Figure 1 —figure supplement 5). For details on the Pgk1-Luc mice see Methods.

Another weakness is the PER2::LUC luminescence measurement from the brain. The authors must show IVIS imaging and show the signal coming from the brain and not from other tissues.

We performed the experiment the reviewer suggested in a previous publication using IVIS imaging (Curie et al. Sleep; 2015 doi: 10.5665/sleep.4974). As described in the Results section these experiments were performed in C57BL6/J mice with pigmented pelage absorbing emitted photos. In addition, the substrate luciferin was administrated only centrally and the by far largest (i.e., only) source of photons in this setup is the glass cone directly mounted on top of the skull over the left cortical hemisphere. Moreover, in the same publication in situ immune-histofluorescence indicated the cerebral cortex as the main source of the almost 2-fold increase in brain bioluminescence after a 6h sleep deprivation. This work was cited when we first stated that brain signal represents cortical activity but now have stated this more explicit (Lines 126-130).

Abstract: Line 18, 19: instead of saying PER2 increases or PER2 rhythmicity, I suggest using PER2::LUC luminescence increases, PER2::LUC luminescence rhythmicity.

We changed the abstract to “PER2 bioluminescence” increases and rhythmicity, consistent with the terminology used in the rest of the manuscript.

Line 89, Suppl. Figure 1: should be Suppl. Figure 3.

We intended referring to the protocol figure (previously Suppl. Figure 1) here but now have added the reference to the other figure as suggested.

Figure 1B: The hypnogram only includes sleeping and waking, please include REM and NREM. I understand that the authors only wanted to indicate sleep/wake transition used for the luminescence analysis indicated in C, but I think the authors should plot standard hypnogram.

As the reviewer noted, this is not a classical hypnogram but our way of indicating which sleep-wake and wake-sleep transitions met our criteria and contributed to the analyses. In the manuscript we made no statements on the two substates of sleep (REM sleep and NREM sleep), and their respective possible effects on PER2 bioluminescence. Moreover, adding separate tracks for REM sleep and NREM sleep will clutter the figure and might confuse the reader. We therefore prefer to keep the figure as is but we now give values of REM sleep and NREM sleep obtained in the recordings on which Figure 1B is based in the text (Lines 122-125) to assure the reader that NREM and REM sleep were expressed at expected quantities under the different recording conditions.

Figure 1C: To avoid confusion, it is better to use the same colors as in A and B (pink: peripheral, purple: central).

In this figure panel we now have used the same red and purple hues for peripheral and central transition data, respectively.

Figure 2: The authors should show individual records in the supplemental information. I know the authors discussed two possibilities that can explain changes in amplitude in group average data. The amplitude of an individual mouse was shown in Figure 2C. But I think, like me, the readers will still want to see the individual record. The amplitude of 3 cycles in REC should be analyzed separately (Plot as REC day 1, REC day 2 and REC day 3). Technically REC day 1 is the same as SD, but it is good to show how the amplitude in the group average and the individual mouse declines.

We have added these data to Figure 2 —figure supplement 2.

Line 573-584: This is an important discussion addressing the possible temperature artifact. However, all statements are based on analysis, "data not shown", unfortunately. The authors must show the data.

We have added these data to the new Figure 1 —figure supplement 6.

Line 763: I think the activity data was binned, not averaged.

We indeed binned the data into (circadian) hourly intervals. However, values in corresponding circadian hours were then averaged across individuals. We have now specified this (see Lines 818-820).

Supplementary Figure 1: I know luminescence recording was done in DD, but please indicate this, as the lighting condition in the other treatments was stated.

We have added this information to the figure.

Supplementary Figure 3B: Please indicate individual records as well, so the relationship between luminescence and activity can be evaluated.

We have now provided an example of individual recordings of the CAG-Luc bioluminescence in Figure 1 —figure supplement 5, as well as an assessment of the association between changes in activity and bioluminescence in the CAG-Luc mice.

Reviewer #3:

Response to the comment #1

It is my understanding that almost all of the conclusions made by the authors are based on the temporal correlation between kidney-derived bioluminescence and sleep-wake states. Therefore, I do not feel convinced about these conclusions without additional experimental evidence suggesting the causal link between spontaneous wakefulness and PER2 expression.

Response to the comment #3

The authors claim that spontaneous wakefulness-driven PER2 increase is one of the major findings. Therefore, under not sleep deprivation but spontaneous wakefulness conditions, they need to confirm PER2 expression in other organs and apply expression data to their mathematical model.

We had already included brain (cortex) derived PER2-bioluminescence dynamics at sleep-wake transitions in our revised manuscript submitted in May 2021 (see Figure 1). As similar dynamics were observed at sleep-wake transitions in both cortex and kidney, we demonstrate that our observations go beyond the kidney and also beyond a mere correlation as PER2 bioluminescence consistently follows spontaneous as well as experimentally induced wakefulness in both tissues. Moreover, we were able to predict the brain dynamics of PER2 with the model parameters set for kidney (Lines 329-331), confirming that at least mathematically sleep-wake driven PER2 dynamics in kidney and cortex underlie similar (but not identical) rules. See next our reply concerning the phase of the SCN.

A few previous studies, although not cited in the present version, clearly demonstrated that arousal levels affect the circadian clock via the SCN (for example, PNAS 2016 Yamakawa el al). I think that the present data may not make obvious progress from these previous achievements.

The reference mentioned does not take position on clock-gene rhythms because their expression was not quantified. Our work highlights the importance of the sleep-wake distribution when evaluating clock-gene rhythms in tissues peripheral to the SCN (including the cortex) in the mouse. Moreover, altered ‘arousal levels’ such as associated with a sleep deprivation (Vassalli and Franken, PNAS 2017; doi: 10.1073/pnas.1700983114) or with food-restriction conditions, which redistributes activity patterns (Damiola et al., Gen Dev 2000; doi: 10.1101/gad.183500), does not affect clock-gene rhythms in the SCN in the mouse (see Discussion Lines 520-521). Therefore, at least under entrained conditions, the phase of clock-gene rhythms in the SCN remains locked to the light-dark cycle, while phase in tissues peripheral to the SCN (including the cortex) can be dissociated from that in SCN by e.g. an altered sleep-wake distribution.

Many previous studies have indicated that circadian phase in clock gene expression varies among peripheral organs and tissues, and it is therefore difficult to generalize the temporal correlation between single organ-derived bioluminescence and sleep-wake states to other peripheral and central clocks. In particular, circadian phase of the SCN is largely phase-advanced in comparison with that of peripheral clocks, raising a question on the author's thought that spontaneous wakefulness is a similar trigger of PER2 in both peripheral and central clocks.

[Editors’ note: what follows is the authors’ response to the third round of review.]

There are only two final pending issues that we anticipate should be easy to incorporate as part of the results/discussion and be included in the final rebuttal letter, without the need for any additional experiments.

Essential revisions:

1) Although the new data has clearly enriched the paper, this still raises the concern of one of the reviewers regarding the presence of low amplitude baseline rhythm in luminescence (particularly as CAG-luc appears to be activated by wake: the reviewer expects to see a basal circadian change of luminescence in CAG-luc mice).

This might indeed be expected provided that the CMV activation upon waking up (i) does not vary with circadian time and (ii) does not again decreases as wakefulness progresses. From the example in Figure 1 —figure supplement 5 it seems that the wake-related activation is strongest after consolidated bouts of sleep (i.e. rest phase), and weakest during the active phase, when the fewer and shorter periods of sleep do not appreciatively decrease the CAG-luc signal. Therefore, although CAG-luc bioluminescence increases after sleep-wake transitions, a circadian modulation in this response might explain the lack of a robust modulation of the average hourly bioluminescence levels. If examining the precise dynamics linking the bioluminescence signal of this synthetic promotor to the sleep-wake distribution is considered of interest by the scientific community, this should be addressed experimentally – however this clearly falls outside the scope of our PER2 work. We only invested in recording CAG-luc mice to assess which route of luciferin administration (drinking water versus osmotic mini-pumps) was best.

2) Without showing an absolute luminescence value emitted from the brain, liver, and kidney in the same mouse, it cannot be utterly said that the luminescence signal from the entire mouse is just coming from the cortex. Part of this concern is based on the fact that the Curie et al., 2015 work not only shows signals (in B6 mice) coming from the cortex, but also clear signals coming from the kidney/liver (as visualized with the IVIS signals, whereas the RT-Biolumicorder does not allow sorting those signals).

In the Curie et al. 2015 publication we indeed recorded bioluminescence in the brain, liver, and kidney of B6 mice. However, recording of peripheral and central signals has to be done in separate animals for two reasons. First, where luciferin is delivered (central versus peripheral) is key to where bioluminescence is detected, as luciferin does not readily cross the blood-brain barrier (Bakhsheshian et al. PNAS 2013; PMID: 24297888). To further facilitate central detection of bioluminescence, a glass cone was directly mounted on a thinned skull allowing passage of photons generated by the underlying cortex. Second, the body of the B6 mice, the strain we used specifically to detect central PER2 bioluminescence, was not shaved, thereby preventing the release of photons. In B6 mice thus prepared (not shaved, a cannula to infuse luciferin centrally, and a glass cone placed over the cortex), the PER2-bioluminescence signal by and large represents photons emitted by the brain (see Figure 1—figure supplement 2). The tissue-specific PER2-bioluminescence dynamics observed after sleep deprivation, which were confirmed by western blot analyses and brain immune-histofluorescence, confirm that the bioluminescence originated from different sources. We now have further emphasized in the Results section the specific preparation of the animals used to record central PER2 bioluminescence (Lines 102-106), detailed the IVIS imaging experiment in the Methods (Lines: 681-686), and added to Figure 1 —figure supplement 2 the image of central bioluminescence detected with the IVIS setup, clearly showing that besides the head no other above-background signal could be detected from the animal.

Author details

1. Marieke MB Hoekstra

   Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland

   Present address

   UK Dementia Research Institute at Imperial College London, Division of Brain Sciences, London, United Kingdom

   Contribution

   Conceptualization, Data curation, Formal analysis, Investigation, Project administration, Visualization, Writing – original draft, Writing – review and editing

   Contributed equally with

   Maxime Jan

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-0723-2026

2. Maxime Jan

   Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland

   Contribution

   Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – review and editing

   Contributed equally with

   Marieke MB Hoekstra

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-6483-7430

3. Georgia Katsioudi

   Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland

   Contribution

   Investigation, Methodology, Writing – review and editing

   Competing interests

   No competing interests declared

4. Yann Emmenegger

   Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland

   Contribution

   Investigation, Methodology

   Competing interests

   No competing interests declared

5. Paul Franken

   Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland

   Contribution

   Conceptualization, Data curation, Formal analysis, Funding acquisition, Project administration, Supervision, Visualization, Writing – review and editing

   For correspondence

   paul.franken@unil.ch

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-2500-2921

• 1,138

  Page views

• 184

  Downloads

• 5

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.