# A highly tunable dopaminergic oscillator generates ultradian rhythms of behavioral arousal

Ian D. Blum, Lei Zhu, Luc Moquin, Maia V. Kokoeva, Alain Gratton, Bruno

Giros, Kai-Florian Storch

► To cite this version:

Ian D. Blum, Lei Zhu, Luc Moquin, Maia V. Kokoeva, Alain Gratton, et al.

HAL Id: hal-01316664

https://hal.sorbonne-universite.fr/hal-01316664

Submitted on 17 May 2016

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# A highly tunable dopaminergic oscillator generates ultradian rhythms of behavioral arousal

Ian D Blum<sup>1,2,4</sup>, Lei Zhu<sup>1,2</sup>, Luc Moquin<sup>2</sup>, Maia V Kokoeva<sup>3</sup>, Alain Gratton<sup>1,2</sup>, Bruno Giros<sup>1,2,5</sup>, Kai-Florian Storch<sup>1,2*</sup>

<sup>1</sup>Department of Psychiatry, McGill University, Montreal, Canada; <sup>2</sup>Douglas Mental Health University Institute, Montreal, Canada; <sup>3</sup>Department of Medicine, McGill University, Montreal, Canada; <sup>4</sup>Integrated Program in Neuroscience, McGill University, Montreal, Canada; <sup>5</sup>INSERM UMR S1130 CNRS, UMR8246, Sorbonne University, Paris, France

**Abstract** Ultradian (~4 hr) rhythms in locomotor activity that do not depend on the master circadian pacemaker in the suprachiasmatic nucleus have been observed across mammalian species, however, the underlying mechanisms driving these rhythms are unknown. We show that disruption of the dopamine transporter gene lengthens the period of ultradian locomotor rhythms in mice. Period lengthening also results from chemogenetic activation of midbrain dopamine neurons and psychostimulant treatment, while the antipsychotic haloperidol has the opposite effect. We further reveal that striatal dopamine levels fluctuate in synchrony with ultradian activity cycles and that dopaminergic tone strongly predicts ultradian period. Our data indicate that an arousal regulating, dopaminergic ultradian oscillator (DUO) operates in the mammalian brain, which normally cycles in harmony with the circadian clock, but can desynchronize when dopamine tone is elevated, thereby producing aberrant patterns of arousal which are strikingly similar to perturbed sleep-wake cycles comorbid with psychopathology.

DOI: 10.7554/eLife.05105.001

## Introduction

Ultradian rhythms with periods ranging from one to several hours have been linked to various aspects of mammalian physiology. Usually superimposed on the 24-hr diurnal or circadian rhythm, ultradian oscillations have been observed in the context of locomotion, sleep, feeding, body temperature, serum hormones, and brain monoamines in species ranging from fruit flies to humans (Tannenbaum and Martin, 1976; Ibuka et al., 1977; Rusak, 1977; Daan and Slopsema, 1978; Honma and Hiroshige, 1978; Dowse et al., 1987; Rivkees, 2003; van Oort et al., 2007; Dowse et al., 2010; Seki and Tanimura, 2014). These physiological cycles most frequently exhibit periods of 2–6 hr, adopting harmonics of the 24-hr daily light–dark cycle or the endogenous circadian rhythm, when external timing cues are absent. However, the generation of such ultradian rhythms does not depend on a functional circadian system nor a light:dark cycle. Ultradian locomotor oscillations persist in rodents housed in constant darkness even upon ablation of the suprachiasmatic nucleus (SCN), the site of the master circadian pacemaker (Ibuka et al., 1977; Rusak, 1977), or genetic disruption of the circadian clock (Vitaterna et al., 1994; Bunger et al., 2000) (Figure 1A,B). These ultradian activity cycles may not simply be driven by metabolic demand since the 2- to 3-hr rhythm in foraging activity observed in the common vole persists even in the absence of food (Gerkema and van der Leest, 1991). Studies in this species further indicate that one adaptive value of ultradian activity rhythms may lie in the facilitation of social synchrony, which is suggested to reduce predator risk in this species

*For correspondence: florian.storch@mcgill.ca

**Competing interests:** The authors declare that no competing interests exist.

**Funding:** See page 19

**Received:** 10 October 2014

**Accepted:** 28 December 2014

**Published:** 29 December 2014

**Reviewing editor:** Richard D Palmiter, Howard Hughes Medical Institute, University of Washington, United States

© Copyright Blum et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
**eLife digest** The sleep-wake cycle of mammals is controlled by a 'circadian clock' within the brain, which is synchronized to the day-night cycle. However, other aspects of mammalian physiology including alertness and activity levels, as well as appetite and body temperature—fluctuate in cycles that repeat every few hours. These cycles are known as ultradian rhythms, and they may offer survival benefits by enabling potentially risky behaviors, such as foraging, to be coordinated between members of a group.

Despite their widespread nature and the fact that they appear to be conserved in evolution, virtually nothing is known about the molecular basis of ultradian rhythms. Blum et al. have now identified a second internal clock within the brain, which they name 'the DUO', and shown that this clock normally works in concert with the circadian clock to regulate daily patterns of activity and alertness.

Experiments in mice revealed that the DUO uses the brain chemical dopamine to generate bursts of activity roughly every four hours. Moreover, it continues to work when the circadian clock has been destroyed. Measurements of dopamine in freely moving mice showed that levels of the chemical fluctuate in synchrony with the animals' activity levels. Moreover, drugs that flood the brain with dopamine, such as methamphetamine, disrupt the 4-hour cycle by lengthening the period between bursts of activity, whereas drugs that block dopamine receptors have the opposite effect.

As well as revealing a mechanism by which the brain coordinates processes that repeat several times per day, the identification of the DUO could also provide insights into the biological basis of psychiatric disorders. Conditions such as schizophrenia and bipolar disorder are often accompanied by disturbances in patterns of activity and rest. While these have previously been attributed to the disruption of circadian rhythms, there is little direct evidence for this, which raises the possibility that these changes might instead reflect the disruption of ultradian rhythms.

DOI: 10.7554/eLife.05105.002

(Gerkema and Verhulst, 1990). However, despite their prevalence and hypothesized biological significance, ultradian locomotor rhythms have received little research attention, likely owing to their frequently masked expression and unstable periodicity in contrast with circadian activity rhythms (Ruis et al., 1987; Schibler, 2008) (Figure 1C,D).

In addition to ultradian oscillations, methamphetamine-induced rhythms of locomotor activity occur in the absence of a functional circadian clock. When provided in the drinking water, methamphetamines (Meth) induce a daily activity bout in addition to the expected circadian rhythm (Honma et al., 1986). This Meth-dependent component—which typically adopts a period in the circadian range—is not abolished upon SCN lesion nor by genetic disruption of circadian clock function (Honma et al., 1987; Mohawk et al., 2009). It was thus concluded that a methamphetamine-sensitive circadian oscillator (MASCO) outside the SCN exists which is capable of driving daily cycles of locomotor activity (Tataroglu et al., 2006). Despite the longstanding recognition of ultradian and Meth-dependent rhythms, the underlying cellular and molecular identity of the oscillator(s) driving them is unknown.

Here, we provide evidence for a highly tunable dopaminergic ultradian oscillator (DUO) which is continuously operative in the mammalian brain and which, together with the circadian clock, orchestrates the daily pattern of arousal. Our data suggest that dopamine acts as both the principal oscillator output as well as an integral component of the DUO, determining oscillator period. Our findings further indicate that the previously described MASCO represents a long-period manifestation of the DUO resulting from elevated dopamine tone. Importantly, our data support an intriguing proposition: that DUO, rather than circadian clock, dysregulation critically contributes to the sleep-wake abnormalities associated with psychopathology.

## Results

### Dopamine transporter deficiency results in ultradian locomotor period lengthening

To gain insights into the mechanistic basis of ultradian locomotor rhythm generation, we considered that locomotor activity is associated with an awakened state (Welsh et al., 1988) and consequently,
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Inter- and intra-animal variability of ultradian activity rhythms across time in circadian incompetent mice. (A) Representative double-plotted actograms of running wheel activity in SCN-lesioned (SCNx) and *Bmal1*−/− mice in DD. (B) Dot plot of locomotor period length in DD based on Lomb-Scargle periodogram analysis of seven consecutive days of activity recording (N = 65 for *Bmal1*−/− and N = 48 for SCNx; t₁₁₁ = 0.2785, p = 0.78, unpaired t-test). (C) Intra-animal period variability expressed as mean ± SD for each animal, ranked according to mean period length derived from continuous wavelet transforms (CWT) for the 1–12-hr frequency range (same animals and timespans for calculation as in B). (D) CWT-heatmaps showing decibel scaled and normalized amplitude of oscillations according to frequency and time with black traces indicating the ridge of local amplitude maxima.
DOI: 10.7554/eLife.05105.003</figcaption></figure>
the ultradian locomotor rhythms observed in mice that lack circadian clock function (Figure 1) could be interpreted as rhythms of heightened wakefulness or arousal. In mammals, a key role in arousal promotion has been attributed to distinct monoaminergic neuronal populations located in the upper brainstem and midbrain (Brown et al., 2012). While altering extracellular levels of the arousal-associated monoamines serotonin, norepinephrine, or histamine by genetic manipulation has only limited effects on locomotion (Thomas and Palmiter, 1997; Bengel et al., 1998; Xu et al., 2000; Parmentier et al., 2002; Zhao et al., 2006), depleting the brain of dopamine (DA) profoundly abrogates locomotor activity (Zhou and Palmiter, 1995). Moreover, increasing extracellular DA levels induces hyperlocomotion (Giros et al., 1996) and lengthens the time spent awake (Wisor et al., 2001). We therefore speculated that altering DA tone may affect ultradian rhythm generation. To test this, we examined running wheel activity in mice carrying a disruption in the Slc6a3 gene, which encodes the dopamine transporter (DAT). Slc6a3⁻/⁻ mice exhibit hyperdopaminergia due to the lack of DAT-mediated DA reuptake into dopaminergic neurons, leading to a hyperactivity phenotype (Giros et al., 1996; Gainetdinov et al., 1998). As the presence of the circadian clock and/or a light:dark cycle frequently masks ultradian activity rhythms (Schibler, 2008), we assessed the locomotor behavioral consequences of DAT elimination in the absence of the master SCN circadian pacemaker. To do so, we electrolytically lesioned the SCN of Slc6a3⁻/⁻ mice and their wildtype littermates and monitored their running wheel behavior in constant darkness (DD). While control mice (SCNx-Slc6a3⁺/⁺) exhibited ultradian activity rhythms with the expected ~4-hr period, SCNx-Slc6a3⁻/⁻ mice showed rhythms whose periods were three times longer (Figure 2A,B). Analysis of mice that were deficient for both DAT and the essential clock component BMAL1 (Bmal1⁻/⁻, Slc6a3⁻/⁻) corroborated this finding. Bmal1⁻/⁻, Slc6a3⁻/⁻ mice exhibited ~12–14-hr rhythms in locomotor activity, largely phenocopying the SCNx-Slc6a3⁻/⁻ mice, while their Bmal1⁻/⁻, Slc6a3⁺/⁺ littermates showed ~4-hr periods as expected for isodopaminergic mice lacking circadian clock function (Figure 2C,D). Together, these results suggest that DAT removal markedly increases ultradian cycle length. Alternatively, the ~12-hr rhythms observed in SCNx-Slc6a3⁻/⁻ and Bmal1⁻/⁻, Slc6a3⁻/⁻ animals may originate from an independent oscillator, one that is activated by DAT elimination, while the short period ultradian oscillator that operates in DAT intact, SCNx or Bmal1⁻/⁻ animals is disengaged or otherwise obscured.

## The psychostimulant methamphetamine lengthens ultradian locomotor period

In order to corroborate that DAT removal indeed lengthens the period of the ultradian activity cycles, we took into account that DAT-mediated DA uptake can be reversed by the selective action of the psychostimulant methamphetamine (Meth) (Howell and Kimmel, 2008). Because Meth leads to increased extracellular DA concentrations as in the case of Slc6a3 gene disruption, we speculated that Meth treatment would result in a, possibly gradual, period lengthening of the ultradian locomotor rhythms. Indeed when we treated Bmal1⁻/⁻ animals with increasing concentrations of Meth via drinking water in DD, we observed a gradual lengthening of the initial ~4-hr locomotor oscillations and this was accompanied by a corresponding increase in activity bout length (Figure 3A,B). Of note, the period increase in response to elevating Meth concentrations did not halt in the circadian range, rather, oscillations continued to lengthen with individual animals reaching periods of 100 hr or more (Figure 3C). Gradual period lengthening of ultradian locomotor rhythms was also observed in Bmal1⁻/⁻ animals exposed to amphetamine (Figure 3—figure supplement 1A,B), a drug similarly targeting DAT (Howell and Kimmel, 2008), but with lower efficacy than Meth (Goodwin et al., 2009). These results argue that DAT targeting psychostimulants affect an endogenous ultradian rhythm generator by increasing period length. However, due to the mode of delivery (drinking water) and the rhythmic Meth uptake that may in turn result from it, it is conceivable that the Meth-dependent, long-period oscillations are ‘driven’ by rhythmic drug intake rather than being generated endogenously. To address this possibility, we subcutaneously implanted Bmal1⁻/⁻ mice with osmotic minipumps that continuously infused Meth over a period of 2 weeks. Running wheel analysis in DD of Meth-infused Bmal1⁻/⁻ animals demonstrated a significant ultradian period lengthening upon drug infusion (Figure 3—figure supplement 1C,D) suggesting that rhythmic uptake is not required for Meth to exert its period lengthening effect, in line with a previous study performed in rats (Honma et al., 1987). The relatively limited change in period observed in this experimental paradigm could be due to the short, 2-week infusion timespan, which is perhaps insufficient to robustly lengthen periods beyond 12 hr (Figure 3—figure supplement 1D). Together, these findings support the notion that the long-period oscillations observed
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 2.** Dopamine transporter knockout alters periodicity of ultradian locomotor rhythms in circadian incompetent mice under constant darkness. (A) Representative, double-plotted actograms demonstrating marked lengthening of ultradian locomotor periods in *SCNx-Slc6a3*<sup>+/+</sup> mice as compared to *SCNx-Slc6a3*<sup>-/-</sup> littermates. *Tau(τ)* indicates individual period. (B) Period length averages of ultradian activity (N = 7; *F*<sub>1,6</sub> = 253.8, ***p* < 0.0001, ANOVA) from Lomb-Scargle periodogram analysis of 7-day time-spans. (C) Representative, double-plotted actograms demonstrating markedly increased ultradian period lengths in *Bmal1*<sup>-/-</sup>, *Slc6a3*<sup>+/+</sup> mice as compared to *Bmal1*<sup>-/-</sup>, *Slc6a3*<sup>+/+</sup> mice. (D) Period length averages of ultradian activity (N = 4; *F*<sub>1,3</sub> = 194.2, ***p* < 0.001, ANOVA) from Lomb-Scargle periodogram analysis of 7-day time-spans.</figcaption></figure>

DOI: 10.7554/eLife.05105.004

in Meth-treated animals, and likewise in *Slc6a3*<sup>-/-</sup> mice, are due to period expansion of an endogenous ultradian rhythm of arousal.

We next aimed to confirm these findings in intact mice. While activity rhythms in the ultradian range are easily discernable in voles and hamsters (Gerkema et al., 1990; van der Veen et al., 2006; Prendergast et al., 2012), the presence of such rhythms in mice or rats is much less obvious and this
<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 3.** Pharmacological tuning of DAT activity by incrementally increasing methamphetamine concentrations lengthens ultradian locomotor period into the infradian range. (A) Representative actograms of Meth-treated *Bmal1*⁻/⁻ mice in DD. Treatment intervals are highlighted with corresponding concentrations indicated in (B). (B) Mean periods from the last 7 days at a given Meth concentration. Repeated measures ANOVA revealed a significant main effect of concentration/time (*F*<sub>5,40</sub> = 34.30, *p* < 0.001) and significant period lengthening between consecutive concentrations (mean ± SEM; N = 9; **p** ≤ 0.005, planned comparisons; 30 day, 30 day exposure to 100 mg/l). (C) Modulo 110-hr actogram of a *Bmal1*⁻/⁻ animal after extended exposure to Meth revealing an ultra-long activity rhythm.</figcaption></figure>

DOI: 10.7554/eLife.05105.005

The following figure supplement is available for figure 3:

**Figure supplement 1.** Amphetamine treatment lengthens ultradian locomotor rhythms in *Bmal1*⁻/⁻ mice.

DOI: 10.7554/eLife.05105.006

has been attributed to masking by the circadian clock (Schibler, 2008). Indeed, running wheel activity data from intact mice often do not provide strong indications of ultradian activity rhythms. However, when we recorded ambulatory behavior using telemetry, we frequently detected three, evenly spaced activity bouts during the active (night) phase (Figure 4A). The observed ~4 hr peak-to-peak spacing is in line with a previous study on several circadian competent mouse strains where ultradian rhythms with similar period length were detected (Dowse et al., 2010). The absence of a clearly discernable ~4-hr rhythmic component during the light portion of the daily LD cycle (Figure 4A) is likely due to masking by the SCN and/or light. Unlike WT mice, *Slc6a3*⁻/⁻ mice never showed a triple-peak activity pattern at night (Figure 4B), and spectral density analysis in the ultradian (2–8 hr) range (Figure 4C) underscores that rhythm generation is significantly altered in these animals while total daily activity is not (Figure 4D). When we provided C57BL/6 mice carrying a telemetry implant with Meth in their
<figure><img src="image_14.png" /><figcaption>**Figure 4.** Ultradian activity in SCN-intact *Slc6a3*⁻/⁻ mice and their wildtype littermates. (A and B) Ambulatory activity recorded by telemetry implants across multiple days (left) and averaged over 24 hr (right). Traces represent 2-hr recursive smoothing (black) of the underlying raw DATa (dark grey; SEM envelope, light grey). Areas in yellow indicate lights on. (C) Amplitude spectral density in the ultradian range (2–8 hr) is significantly different between *Slc6a3*⁻/⁻ mice and their wildtype littermates revealing a loss of the ultradian component (mean ± SEM; N = 10; *F*<sub>1,18</sub> = 26.40, **p** < 0.0001, ANOVA). (D) Although the temporal pattern of locomotion differs between genotypes there is no significant difference in daily activity averaged over multiple days (mean ± SEM; N = 10, *F*<sub>1,18</sub> = 0.1793, ANOVA). (E–G) Addition of Meth to the drinking water of C57BL/6 mice lengthens the night-time activity bouts in a concentration-dependent manner. Averaged daily locomotor activity of individual mice at different Meth concentrations (G) derived from the time-span indicated by colored bars next to the representative actogram (F) and subjected to Butterworth filtering. The three night-time activity peaks before treatment (white triangles), transform to 2 peaks after exposure to the highest concentration (black triangles). The night-time bout lengthening is also reflected in the reduction of ultradian amplitude spectral density in the 1 to 5-hr range (E, mean ± SEM, N = 7; *F*<sub>2,20</sub> = 8.08, p < 0.005, repeated measures ANOVA; \*p < 0.05, \*\*p < 0.01, post-hoc Bonferroni). White asterisks (in E) indicate cage changes.
DOI: 10.7554/eLife.05105.007</figcaption></figure>

drinking water, we observed a gradual lengthening of the interval between the night-time activity peaks (Figure 4F,G), resulting in the transformation of the triple peak (Figure 4G, white trace and white triangles) into a dual or single night-time activity peak pattern at the highest (100 mg/l) Meth concentration (Figure 4G, black trace and black triangles). This was corroborated by spectral density analysis revealing a marked reduction in ultradian frequencies upon Meth treatment (Figure 4E). Together, these results indicate that blocking DA reuptake also lengthens ultradian activity bout intervals in circadian competent mice, arguing that an ultradian oscillator is not 'activated' by circadian disruption, but rather continuously operative alongside the circadian timer. As with *Bmal1*⁻/⁻ mice, prolonged Meth-treatment in intact WT mice can lead to profound period lengthening of ultradian rhythms, often into the circadian (48 hr) range (for examples see Figures 5E, 7F, 9B).

## The antipsychotic haloperidol shortens ultradian locomotor period

Given that both DAT disruption and Meth treatment elevate extracellular DA and concurrently lengthen ultradian rhythm period, we speculated that manipulations aimed at lowering DA tone
<figure><img src="image_15.png" /><figcaption>**Figure 5.** Haloperidol shortens circadian-clock-independent locomotor rhythms. Representative actograms with Hal treatment periods indicated by colored bars (left); bar graphs indicate corresponding locomotor period (right). (A) Increasing concentrations of Hal provided in the drinking water gradually shortens the endogenous ultradian activity rhythms of *Bmal1*⁻/⁻ mice in DD (B, mean ± SEM, N = 12; *F*<sub>2,10</sub> = 14.36, *p* < 0.0001, repeated measures ANOVA; \*p < 0.05 and \*\*p < 0.01, planned comparison ANOVA). (C and E) Hal shortens the infradian rhythms in Meth-treated *Bmal1*⁻/⁻ mice (D, mean ± SEM, N = 9; *F*<sub>1,8</sub> = 2.357, \*p < 0.05 ANOVA) and WT mice (F, mean ± SEM, N = 9; *F*<sub>1,8</sub> = 3.525, \*\*p < 0.01, ANOVA) in DD. (G) Hal treatment increases the frequency of temperature fluctuation in *Slc6a3*⁻/⁻ mice under LD measured by telemetric implants. (H) Changes of amplitude spectral density in the ultradian range (2–8 hr) in response to increasing Hal concentrations (mean ± SEM; N = 8; *F*<sub>2,7</sub> = 14.74, repeated measures ANOVA, \*p < 0.05, \*\*\*p < 0.0005, post-hoc Bonferroni). For all experiments, periods are calculated based on the running wheel activity during the last 7 days of treatment at the indicated Meth/Hal concentration.</figcaption></figure>

DOI: 10.7554/eLife.05105.008

should conversely lead to period shortening. We therefore provided *Bmal1*⁻/⁻ mice in DD with the antipsychotic drug haloperidol (Hal) in their drinking water. Hal selectively blocks the DA receptor D2 (DRD2) (*Strange, 2008*), which is expressed on postsynaptic target sites of DA neuronal projections
<figure><img src="image_16.png" /><figcaption>**Figure 6.** DA fluctuations correlate with ultradian locomotor behavior in *Bmal1*⁻/⁻ mice. (A) Two representative examples of in vivo striatal microdialysis in *Bmal1*⁻/⁻ mice. Upper: locomotor activity as measured by beam breaks. Lower: DA dialysate concentration measured at 20-min intervals (red trace) plotted alongside the corresponding locomotor activity (solid black trace). Cosinor (dotted line) was computed from the 20 hr of locomotor activity following dialysate sampling. (B) False discovery rate of the fit between the DA profiles and corresponding cosinors (mean ± SEM, N = 7). (C) Linear regression analysis of period length vs mean DA concentration, dots representing individual animals. (D) Tissue punches of CPu, NAcc, and midbrain were analyzed for DA and DOPAC content in animals sacrificed during their active (■) vs their non-active (□) phase (mean ± SEM; active phase, N = 7, inactive phase, N = 6; \*p < 0.05, \*\*p < 0.01, \~p = 0.066, ANOVA). (E) Linear regression for period vs DA content in animals sacrificed during the active phase revealed significant correlations for the CPu and midbrain. (F) Representative actograms displaying activity rhythms of Meth-treated *Bmal1*⁻/⁻ mice used for tissue collection. Triangles indicate collection time points (in °), with locomotor activity bout onset set to 180°. (G) Illustrations of coronal mouse brain sections (*Paxinos and Franklin, 2001*) indicating position of the active membrane (black bar) and tissue punch placements (circles, colors correspond to labels in **D** and **E**).</figcaption></figure>

DOI: 10.7554/eLife.05105.012

The following figure supplement is available for figure 6:

**Figure supplement 1.** Rhythms of extracellular DA in the striatum of freely-moving *Bmal1*⁻/⁻ mice.

DOI: 10.7554/eLife.05105.013

but also presynaptically as an autoreceptor (Schmitz et al., 2003). Importantly, when given chronically as in our case, Hal has been reported to electrically silence midbrain DA neurons (White and Wang, 1983) and to markedly lower extracellular DA levels in the striatum/nucleus accumbens regions in rats (Lane and Blaha, 1987; Ichikawa and Meltzer, 1991). As predicted, *Bmal1*⁻/⁻ mice responded with successive locomotor period shortening to increasing concentrations of Hal (Figure 5A,B). We also observed Hal-mediated shortening of the long-period behavioral rhythms in *Bmal1*⁻/⁻ and wild-type mice treated concurrently with Meth (Figure 5C-F). A period shortening effect of Hal could also be discerned from the core body temperature fluctuations of *Slc6a3*⁻/⁻ mice, which increased in frequency (Figure 5G). This response, which mirrored the locomotor behavioral response (not shown), was concentration dependent and spectral density analysis confirmed a successive increase of the ultradian component upon Hal exposure (Figure 5H). Together, these data strongly support the hypothesis that both the ultradian rhythms observed in wildtype and *Bmal1*⁻/⁻ mice and the long-period rhythms previously attributed to the MASCO are driven by the exact same oscillatory mechanism.
<figure><img src="image_17.png" /><figcaption>**Figure 7.** *Slc6a3*^−/− mice show a second rhythmic locomotor activity component. Representative actograms displaying daily running wheel activity of *Slc6a3*^+/+ (*A*), *Slc6a3*^−/− mice (*D*) and C57BL/6 mice on Meth (*G*). (*B*, *E* and *H*) Lomb-Scargle periodograms generated from the time-span indicated by red bars, in (*A*, *C*, *G*), respectively. (*C*) There is no significant difference between genotypes in daily activity averaged over the time-span of analysis (mean ± SEM; N = 6, *F*<sub>1,10</sub> = 1.848, ANOVA) (*F* and *I*) Average periods of the highest 2 periodogram peaks for *Slc6a3*^−/− mice (*E*, mean ± SEM, N = 6)

*Figure 7. Continued on next page*</figcaption></figure>
Figure 7. Continued

and C57BL/6 mice on Meth (H, mean ± SEM, N = 9, Meth-treatment started on day 1 of the recordings). Areas in yellow indicate lights-on. Green line in the periodograms indicates the confidence threshold for rhythmicity (α = 0.01).

DOI: 10.7554/eLife.05105.009

The following figure supplement is available for figure 7:

**Figure supplement 1.** The SCN of *Slc6a3*⁻/⁻ mice is unperturbed.

DOI: 10.7554/eLife.05105.010

As altering extracellular DA is a common denominator of all three manipulations—DAT elimination, (meth) amphetamine, and haloperidol treatment—these findings collectively suggest that DA tone determines ultradian period. Given that DA is known to mediate arousal/wakefulness (*Brown et al.*, 2012), it appears plausible that DA also serves as principal oscillator output, driving rhythms in arousal. We thus hypothesized that extracellular DA must oscillate in synchrony with the observed activity rhythms.

## Extracellular DA fluctuates in synchrony with ultradian activity cycles

To test this hypothesis, *Bmal1*⁻/⁻ mice were unilaterally implanted with a microdialysis probe positioned along the ventro-dorsal extent of the right striatum (Figure 6G), a site heavily innervated by DA neurons. Mice were kept under constant dim red light (<5lux) and ultradian activity rhythms were monitored using infra-red beam breaking. The animals showed the expected short period rhythms in locomotor activity throughout the experiment (Figure 6A), and analysis of dialysates revealed that extracellular DA levels fluctuated concordantly with the activity cycles (Figure 6A and Figure 6—figure supplement 1, compare solid red and black traces). Importantly, we generated a model oscillation using the cosinor method from locomotor rhythms recorded for 20 hr after dialysate sampling with no investigator present and compared it to the DA fluctuations observed. The goodness of fit between the cosinor and the DA profiles indicated by the sum of squared errors (SSE = 0.050 ± 0.027, mean ± SEM, N = 7) suggests that the dialysate sampling procedure itself did not perturb the generation of the endogenous ultradian activity cycles. To determine the statistical significance of the observed agreement between the cosinor and the DA trace, we randomly permutated the individual time points of the DA trace and determined the percentage of permutated traces with an equal or better fit to the cosinor in comparison to the observed DA fluctuation. On average, only 3.2 ± 1.7% of 100,000 permutated DA profiles correlated as well or better than the experimentally observed DA measurements (Figure 6B), confirming that extracellular DA fluctuates in synchrony with the ultradian activity cycles. Of particular note, linear regression revealed a highly significant correlation between mean DA concentration in the dialysate and ultradian locomotor period in the *Bmal1*⁻/⁻ animals we tested (Figure 6C), which again supports a role of (extracellular) DA as a period determinant of the ultradian activity cycle.

As our findings suggest that the observed >24-hr activity rhythms are due to period lengthening of the ultradian oscillations, we similarly expected DA levels to fluctuate in synchrony with the activity cycles in Meth-treated animals. To test this, we measured DA content in tissue punches of Meth-treated *Bmal1*⁻/⁻ mice from the dorsal (caudate putamen, CPu) and ventral (nucleus accumbens, NAcc) striatum as well as the ventral midbrain, which includes both the substantia nigra (SN) and the ventral tegmental area (VTA) (Figure 6G). For punch collection, animals were sacrificed during their active and inactive phases, respectively (Figure 6F). We detected significantly higher levels of DA as well as its immediate metabolite, 3,4-dihydroxyphenylacetic acid (DOPAC), in tissue extracts of the CPu during the animal's active phase (Figure 6D). There was also a significant increase in DOPAC levels and a trend towards elevated DA (p = 0.066) during the active period in the NAcc, while extracts from the midbrain region, which contains the cell bodies of the striatal dopaminergic afferents, exhibited no detectable change (Figure 6D).

We also found a significant correlation between locomotor period and DA content in CPu and midbrain punches collected during the active phase (Figure 6E), again consistent with a role for DA as a period determinant of ultradian oscillations.

## SCN-intact *Slc6a3*⁻/⁻ mice show two components of rhythmic activity

As *Slc6a3*⁻/⁻ mice lack the triple peak night-time activity pattern that is characteristic of wild-type mice, we subjected *Slc6a3*⁻/⁻ mice to long-term running wheel activity monitoring to examine their
<figure><img src="image_18.png" /><figcaption>**Figure 8.** Chemogenetic activation of midbrain DA neurons lengthens ultradian locomotor period. (A) AAV-DIO-hM3Dq-mCherry was stereotaxically and bilaterally delivered into the VTA/SN region of DAT-Cre transgenic mice as indicated. (B) Representative immuno-fluorescent image of the ventral midbrain from a virus-injected and behaviorally responsive mouse showing extensive co-expression of the mCherry fusion-tag in TH-positive cells of the midbrain. (C) Locomotor response to CNO (red arrow, 1 mg/kg body weight i.p.) of a representative, AAV-hM3Dq transduced DAT-Cre mouse undergoing microdialysis. (D) 20min- binned locomotor activity and extracellular striatal DA content (20 min) of AAV-hM3Dq transduced mice 1 hr prior (Baseline) and 2 hr after CNO injection. Mice were implanted with a striatal microdialysis probe and DA content was measured as in **Figure 6A** (mean ± SEM; N = 3, *p* < 0.05, **p** < 0.01, paired t-test). (E–H) Representative actograms of AAV-hM3Dq transduced DAT-Cre, SCNx (E) and DAT-Cre, Bmal1~−/−~ mice (G). Switch to CNO-containing drinking water (7.5 mg/l) is indicated by red arrow. Blue arrow marks return to regular water (in E). Animal activity is plotted modulo according to the period measured during the last 7 days of CNO treatment. Periodogram analysis reveals CNO-dependent period lengthening in both DAT-Cre, SCNx (F, mean ± SEM; N = 6, F<sub>1,5</sub> = 3.68, *p* < 0.05, ANOVA) and DAT-Cre, Bmal1~−/−~ (H, mean ± SEM; N = 11; F<sub>1,10</sub> = 20.48, **p** < 0.0001, ANOVA) mice. VTA, ventral tegmental area; SNc, substantia nigra pars compacta.</figcaption></figure>

DOI: 10.7554/eLife.05105.014

The following figure supplement is available for figure 8:

**Figure supplement 1.** AAV-h3MDq targeting of DA neurons.

DOI: 10.7554/eLife.05105.015

locomotor behavior in more detail. Under LD conditions, we observed more fragmented daily activity in *Slc6a3*~−/−~ than their wildtype counterparts (Figure 7A,D). Upon release into constant darkness (DD), wild-type animals exhibited a daily rhythm with a period slightly shorter than 24 hr as expected for mice (Bunger et al., 2000) (Figure 7A). While this principal locomotor component was also observed
<figure><img src="image_19.png" /><figcaption>**Figure 9.** Proposed model of circadian and ultradian oscillator output integration to govern daily locomotor behavior. (A) Light input entrains the circadian (SCN) and (indirectly) the ultradian (DUO) oscillators creating a stable phase relationship. Their rhythmic outputs, upon integration at a common downstream effector, generate the daily pattern of locomotor activity. (B) Under conditions of high DA tone (e.g., DAT elimination or Meth-treatment), DUO period lengthens, leading to a second, separate activity rhythm. This rhythm either free-runs (see **Figure 7D, G**) or phase locks with the SCN pacemaker by adopting a subharmonic, that is, 48-hr period, as frequently observed upon Meth treatment. Representative output plots show average activity of individual mice computed from 8 days of ambulatory activity (A) or 14 days of running wheel activity under Meth treatment (B). DOI: 10.7554/eLife.05105.011</figcaption></figure>

in *Slc6a3*⁻/⁻ mice, they repeatedly exhibited an additional activity component with periods longer than 24 hr (Figure 7D-F). This component, which was more evident after the first few weeks in DD, extended from the main activity component and lasted for several cycles before disappearing late in the subjective day (Figure 7D). Overall, there was no significant difference in the total amount of daily activity even though the temporal pattern of locomotion in these mice is severely altered (Figure 7C). The observed dual component activity pattern shows striking similarities to the activity patterns of Meth-treated C57BL/6 mice (Figure 7G-I), suggesting that the second rhythmic component is generated by the dopaminergic oscillator.

However, such locomotor activity pattern, with two components oscillating at different frequencies, has also been observed in rats housed under a 22-hr LD cycle and have been attributed to a dissociation of molecular rhythms between the dorsal and ventral SCN (de la Iglesia et al., 2004). To rule out that the SCN is similarly affected in *Slc6a3*⁻/⁻ mice, we sacrificed mice at the peak and trough of the
main activity component when the second rhythm was maximally antiphasic (Figure 7—figure supplement 1A). In situ hybridization directed against Period1 transcripts did not indicate rhythmic desynchrony between or within SCN hemispheres. The riboprobe showed strong, homogenous staining throughout the SCN in the early subjective day, whereas little to no signal was detectable in the early subjective night, suggesting that SCN pacemaker function is unperturbed in Slc6a3⁻/⁻ mice (Figure 7—figure supplement 1B).

Together, these findings thus argue that the second rhythmic component of Slc6a3⁻/⁻ mice is indeed a product of the above described dopaminergic oscillator.

## Chemogenetic activation of DA neurons lengthens ultradian period

Since striatal and midbrain DA levels vary concordantly with activity rhythms in Bmal1⁻/⁻ mice, and DA neurons of the SN/VTA are by and large the exclusive source of striatal and midbrain DA, these neurons are a likely site of ultradian rhythm regulation, if not generation. To corroborate the importance of these neurons in the oscillator process, we employed a chemogenetic approach based upon DREADD (designer receptors exclusively activated by designer drugs) technology (Krashes et al., 2011). We stereotaxically delivered the adeno-associated virus AAV-DIO-hM3Dq-mCherry (Krashes et al., 2011) into the SN/VTA region of mice that carried an Slc6a3 promoter driven Cre transgene (DAT-Cre) (Figure 8A) (Turiault et al., 2007) and that additionally were either BMAL1-deficient (DAT-Cre, Bmal1⁻/⁻) or received an SCN lesion (DAT-Cre, SCNx). Upon Cre-mediated recombination, virally transfected cells express the stimulatory DREADD hM3Dq, which has been demonstrated to increase neuronal firing frequency upon binding of the compound clozapine-N-oxide (CNO), which is otherwise physiologically inert (Alexander et al., 2009). As augmenting firing frequency in DA neurons enhances DA release (Sulzer, 2011), and CNO stimulation of successfully transduced DAT-Cre mice evokes vigorous firing of dopamine neurons (Wang et al., 2013), CNO is expected to increase DA release. Consistently, i.p.-injection of CNO elevated extracellular DA in the striatum of those mice and led to a prolonged increase in locomotor activity (Figure 8C,D), which was absent in vehicle injected animals (Figure 8—figure supplement 1D). Responsive mice showed mCherry fusion tag expression selectively in tyrosine hydroxylase (TH) positive cells (Figure 8B, Figure 8—figure supplement 1A). Cell counting revealed that, 92.9 ± 5.1% and 87.8 ± 6.0% (mean ± SEM, N = 6 for both groups) of the midbrain TH⁺ neurons also expressed the mCherry fusion tag in responsive DAT-Cre, SCNX and DAT-Cre, Bmal1⁻/⁻ mice, respectively, confirming effective Cre-mediated hM3Dq expression in midbrain DA neurons. Running wheel activity monitoring of virus-injected mice showed the expected short period locomotor oscillations, however, upon switching to CNO-containing drinking water (red arrows), both, AAV-hM3Dq transduced DAT-Cre, SCNx and DAT-Cre, Bmal1⁻/⁻ mice responded with locomotor period lengthening (Figure 8E–H), an effect that was reversed when mice were returned to pure water (Figure 8G, blue arrows). As expected, a period lengthening was not observable in AAV-hM3Dq injected DAT-Cre<sup>0/0</sup>, Bmal1⁻/⁻ mice upon CNO exposure (Figure 8—figure supplement 1B,C). These results further corroborate that elevating DA tone, in this case by selective excitation of DA neurons, lengthens ultradian locomotor period and that the oscillator driving ultradian rhythmicity comprises DA neurons of the VTA/SN region.

## Discussion

Collectively, our results provide strong evidence that a dopaminergic ultradian oscillator (DUO) driving rhythms of behavioral arousal is continuously operative in the mammalian brain. We propose that under normal conditions, this DUO cycles in harmony with the circadian SCN pacemaker and that the rhythmic information of both the SCN and the DUO are integrated at a common downstream site to create the daily pattern in locomotor activity (Figure 9A). However, elevation of DA tone can lead to DUO period lengthening, which either results in DUO free-run (for example, Figure 7D,G) or reinstatement of oscillator synchrony albeit at a different harmonic (Figure 9B, e.g., 48 hr). The DUO appears as a highly tunable oscillator, able to adopt period lengths from a few hours to multiple days. This is in stark contrast to the circadian timer which cannot adopt periods that are more than a few hours off from 24 hr when its limits of entrainment are tested experimentally (Aschoff and Pohl, 1978).

Given that DA is known to stimulate locomotor activity (Zhou and Palmiter, 1995), our observation of a cyclical rise in extracellular striatal DA, which is in synchrony with ultradian activity rhythms, argues that DA acts as an output of the DUO. Our finding that manipulations affecting extracellular DA levels alter oscillator period and that extracellular DA tone shows a remarkably high correlation with activity
cycle length, strongly suggests that DA is a period determinant and therefore must be an integral
component of the oscillatory mechanism itself. As all period-altering manipulations directly impinge
upon DA neuronal physiology, this suggests that either (i) DA neurons are the site of ultradian rhythm
generation or (ii) they are a key cog in the oscillatory mechanism. The chemogenetic activation experi-
ments indicate that the relevant DA neuronal population is located in the midbrain as selective activa-
tion of DA neurons in this region had a period lengthening effect on the ultradian activity. Future
experiments will be aimed at delineating the precise DA population(s) required for the ultradian
rhythm generation process and how rhythmic synchrony between neurons is maintained, if the DUO is
indeed composed of a population of cellular oscillators.

Our data indicate that the circadian and ultradian locomotor rhythms are normally harmonized (for instance, Figure 9A), suggesting that the circadian pacemaker and the DUO interact. Of note, extracellular DA was reported to fluctuate diurnally in the rodent striatum (Hood et al., 2010), which was abrogated in Slc6a3⁻/⁻ mice (Gallardo et al., 2014). In these studies, only group averaged DA profiles were presented and DA was solely measured in circadian intact mice and rats, thus it is conceivable that any ultradian component (in WT) or infradian component (in Slc6a3⁻/⁻) escaped detection. Critically, the observation that the DA levels, on average, followed a diurnal pattern with a night-time peak, together with the observation that ultradian locomotor period in female hamsters is longer in the dark vs the daily light phase (Prendergast et al., 2012), is consistent with the LD cycle and/or the circadian pacemaker affecting DUO rhythmicity by altering (extracellular) DA. Note however, that dopamine content in whole brain extracts from circadian-intact rats kept under LD showed a strong ultradian variation with no obvious evidence for a diurnal rhythmic component (Scheving et al., 1968). Also, microdialysis did not reveal a day:night difference in extracellular DA levels in the striatal NAcc region when measured in DD (Chung et al., 2014). This same study reported elevated DA levels and hyperactivity in mice lacking clock gene Rev-erba, suggesting that the circadian clock, possibly intrinsic to DA neurons themselves, has a role in DA regulation. Indeed, knockdown of the core circadian clock component Clock in DA neurons of the murine VTA, increases electrical firing rate in VTA neurons and enhances locomotor response to novel objects (Mukherjee et al., 2010). Interestingly, we did not observe any systematic differences in ultradian locomotor period between SCNx and Bmal1⁻/⁻ mice in constant darkness conditions, regardless of DAT status (Figures 1, 2), suggesting that extra-SCN circadian clocks have no role in DUO-mediated ultradian locomotor rhythm generation.

Our data also provide evidence that the postulated MASCO (Tataroglu et al., 2006) represents a long-period manifestation of the DUO as a result of methamphetamine's action on the dopamine transporter, blocking dopamine reuptake thereby increasing extracellular DA levels. Interestingly, the dopamine system has been also implicated with another behavioral timing system: the food-entrainable oscillator (FEO). This circadian independent oscillator (Pitts et al., 2003; Pendergast et al., 2009; Storch and Weitz, 2009) drives food-anticipatory locomotor activity (FAA) that emerges when food access is restricted to a few hours each day (Mistlberger, 2011). D1 (DRD1) and D2 (DRD2) receptor antagonists attenuate FAA additively (Liu et al., 2012) while pharmacological DRD2, but not DRD1, activation altered the phase of FAA (Smit et al., 2013). Most recently, using knockout mice, it was revealed that DRD1 but not DRD2 is necessary for the appropriate expression of FAA and that rescuing dopamine signaling selectively within the dorsal striatum was sufficient to restore FAA in dopamine-deficient mice (Gallardo et al., 2014). Given the links to the dopaminergic system, it will be of interest to investigate whether the DUO has a role in the temporal regulation of FAA.

While these findings argue that food cues engage the dopamine system to alter daily locomotor activity patterns, it is clear that dopamine signaling also affects food intake as mice lacking dopamine are lethargic and do not actively consume food (Zhou and Palmiter, 1995). Considering that the DUO is a universal driver of ultradian behavioral rhythms in mammals, the finding that ultradian bouts of running-wheel and feeding activity are co-expressed in the common vole (van der Veen et al., 2006) suggests that dopamine can synchronously drive food seeking and general activity, which is in line with the view that dopamine acts as a general promoter of motivated arousal.

Slc6a3⁻/⁻ mice have been proposed as a model for schizophrenia (Gainetdinov et al., 2001) and the DA hypothesis of schizophrenia states that DA elevation is causal to the behavioral symptoms of this psychiatric condition. Intriguingly, circabidian (48 hr) or free-running rhythms in locomotor activity reminiscent of the behavioral patterns we detected in Slc6a3⁻/⁻ or Meth-treated mice (Figures 4F, 7C,F, 9B) have been observed in schizophrenic subjects (Wirz-Justice et al., 2001; Wulff et al., 2012),
suggesting that DUO dysregulation underlies the rest-activity aberrations associated with schizophrenia. Furthermore, actigraphy recordings in schizophrenic patients also revealed that Hal treatment reduces circadian/diurnal locomotor amplitude and leads to the emergence of ultradian activity bouts (Wirz-Justice et al., 1997, 2009), effects we likewise observed in *Slc6a3<sup>-/-</sup>* mice in response to Hal. Transient DUO period lengthening might equally account for the rest-activity pattern abnormalities observed during manic episodes in bipolar disorder, which has been also associated with altered DA tone (Berk et al., 2007). Bipolar subjects have been reported to show rapid, 48-hr cycling between mania and depression (Gann et al., 1993; Wilk and Hegerl, 2010), one to multiple 48-hr sleep-wake cycles when switching from depression to mania (Wehr et al., 1982), or a long-period 'free-running' rhythm of wakefulness (Wehr et al., 1998). Thus schizophrenic and bipolar subjects both appear to exhibit rest-activity cycle aberrations strikingly similar to those observed in *Slc6a3<sup>-/-</sup>* or Meth-treated mice, suggesting that DUO dysregulation is a common indicator for these psychopathologies and perhaps even a common disease cause.

A switch to sleep-wake cycles with a period much longer than 24 hr have also been observed in subjects that were studied in temporal isolation (Aschoff, 1965). Because other physiological parameters such as urine secretion and core body temperature showed phase-aligned, circadian fluctuations with a period much closer to 24 hr, the subjects were considered internally desynchronized. Notably, affected subjects tend to exhibit high scores of neuroticism (Wever, 1979). It is therefore conceivable that the observed internal desynchronization is also due to a dysregulation of the DUO.

# Materials and methods

## Animals

*Bmal1<sup>-/-</sup>* (Bunger et al., 2000) and *DAT-Cre* (Turiault et al., 2007) mice were on a C57BL/6J genetic background, while *Slc6a3<sup>-/-</sup>* mice (*Giros et al., 1996*) were maintained on a mixed C57BL/6JBA/2J background (*Morice et al., 2004*). Animals were housed under LD 12 hr:12 hr unless otherwise stated. *Slc6a3<sup>-/-</sup>* mice were found to exhibit high attrition rates (~80%) in experiments that involve long-term locomotor activity monitoring. To increase survival, *Slc6a3<sup>-/-</sup>* mice were first group-housed for 1 week after transfer into the light-tight, ventilated cabinets used for activity monitoring. This was followed by 1 week of individual housing in running wheel cages with the wheels locked before commencing baseline activity recordings. During this adaptation phase and throughout the subsequent experimental period, both *Slc6a3<sup>-/-</sup>* mice and their wild-type littermates were provided with ad libitum chocolate flavored chow (Supreme Mini-Treats, BioServ, Flemington, NJ) as well as cotton nestlets and ample shredded corrugated card stock. This regimen significantly reduced attrition rates to less than 20% on average. 'Material and methods' were performed in accordance with the Canadian Council on Animal Care guidelines and approved by the McGill University Animal Care Committee.

## Locomotor activity and core body temperature monitoring

*Running wheels*: Animals were individually housed in light-controlled cabinets and activity was recorded continuously (ClockLab, Actimetrics, Wilmette, IL). Actograms, displaying binned running wheel revolutions per 6 min (0.1 hr), and the associated Lomb-Scargle periodograms, displaying amplitude, were generated using ClockLab software. *Telemetry*: Animals were individually housed in standard cages placed atop energizer/receiver units (ER-4000, Starr Life Science Corp., Oakmont, PA). 1 week prior to data collection, electromagnetic induction powered telemetry probes (G2 E-mitter, Starr Life Science Corp.) were implanted intraperitoneally. Locomotion, measured in counts per minute, and core body temperature (BT, in °C), was collected in 6-min bins (0.1 hr) using Vitalview software (Starr Life Science Corp.). BT data was exported into Clocklab to generate actogram-style data displays with tick mark height corresponding to temperatures from 34-38°C; for waveform generation, averaging, and 2 hr recursive smoothing, locomotor activity and BT data was exported to Excel (Microsoft, Redmond, WA). Matlab (Mathworks, Natick, MA) was used for ribbon plot generation, and low pass Butterworth filtering (1 hr).

## Pharmacology

Stock solutions of (+)-Methamphetamine hydrochloride (100 mg/l, Sigma-Aldrich, St. Louis, MO), D-Amphetamine hemisulfate (100 mg/l, Sigma-Aldrich), Haloperidol (11.25 mg/ml, Sigma-Aldrich) and Clozapine-N-Oxide (15 mg/l, National Institute of Mental Health, Bethesda, MD) were prepared
using distilled water. Methamphetamine solutions were adjusted to pH 7 using sodium hydroxide and haloperidol was dissolved by stirring at 40°C.

## Chronic methamphetamine infusion

Mice were continuously infused with methamphetamine (0.6 mg/day in 0.9% saline) for 14 days using subcutaneously implanted osmotic minipumps (Alzet Model 1002, Durect Corp., Cupertino CA). Pumps were fitted with a 65-mm polyvinyl chloride catheter (0.69-mm inner diameter, Plastics One, Roanoke, VA) and backfilled with saline to allow for an approximately 4 day post-surgery recovery prior to drug exposure.

## In situ hybridization

In situ hybridization was performed as previously described (Kraves and Weitz, 2006). Briefly, following decapitation, brains were removed and quickly frozen at -80°C. Serial coronal hypothalamic brain slices (25 µm) were collected using a cryostat and stored at -80°C until hybridization. Sections were hybridized overnight at 60°C to a digoxigenin-labeled riboprobe targeting the coding regions of mouse *Per1* (nucleotides 579-1478 of the *Per1* mRNA, Genbank, NM_011065.4).

## SCN lesions

Electrolytic lesions of the suprachiasmatic nucleus were performed as described (Storch et al., 2007). Briefly, an electrode (RNE-300X, Rhodes Medical Instruments, Woodland Hills, CA) was lowered through a hole drilled in the skull at the mid-sagittal sinus according to stereotaxic coordinates (AP -0.25 mm, DV -6.00 mm from Bregma) and a constant current (2 mA, 10 s; D.C. Constant Lesion Maker, Grass Instruments, Quincy, MA) was applied. Only mice with behavioral circadian arrhythmia (in the 20-28-hr range assessed by Lomb-Scargle periodogram analysis) and subsequent post-mortem histological verification using DAPI staining mounting medium (Vectashield, Vector Labs) were included for analysis. This represented approximately 50% of lesioned mice when averaged across all studies.

## In vivo microdialysis

1 week prior to sampling, mice were stereotaxically implanted with a guide cannula (C312G/spc 2.5 mm below pedestal, Plastics One), targeting the striatum (AP +1.10 mm, DV +5.5 mm, ML +0.9 mm). 1 day prior to dialysis, mice were transferred to a beam break monitor (VersaMax, Omnitech Electronics, Columbus, OH) and tethered using a dummy probe assembly under constant, dim, red light (<5lux). On the day of sampling, in-house recording probes (Lupinsky et al., 2010), were connected to a central swivel (Model 72-0000, Harvard Apparatus, Holliston, MA). Artificial cerebrospinal fluid was delivered via syringe pump (Model 403, CMA Microdialysis, Kista, Sweden) at a flow rate of 1.0 µl/min. This procedure was similarly followed to test the effects of CNO in AAV-transduced *DAT-Cre* mice where intraperitoneal injections of 1 mg/kg CNO were given after approximately 1 hr of baseline sampling. Samples were collected every 20 min for 4 hr and immediately admixed with 1 µl of perchloric acid. DA and DOPAC content in the dialysate were determined by high-performance liquid chromatography with electrochemical detection (HPLC-EC) as described (Domenger et al., 2012). Chromatographic peak analysis was conducted using CoulArray software (ESA Inc., Chelmsford, MA). After data collection, brains were removed, sectioned, and stained with hematoxylin to verify probe placement. Locomotor activity data, recorded as the number of beam breaks per minute, were exported into Matlab (Mathworks) for waveform generation, 20-min binning, cosinor modeling, and false discovery rate analysis.

## DA/DOPAC content determination in tissue extracts

After decapitation, brains were quickly removed and frozen (-80°C). 320 µm coronal brain slices were obtained by cryosectioning and then microdissected (1- or 2-mm diameter sample corers, Fine Science Tools Inc., Foster City, CA) to obtain tissue of the CPu (2 mm), NAcc (1 mm), and SN/VTA midbrain (2 mm) regions (see Figure 6G for punch location). DA/DOPAC was quantified as previously described (Domenger et al., 2012). Briefly, individual punches from each region were homogenized in 45 µl perchloric acid (0.25 M) to which 15 µl of a 100 ng/ml solution of 3,4-dihydroxybenzylamine was added, which served as the internal standard. Concentrations were determined by HPLC-EC. After perchloric acid extraction, the protein containing pellets were reconstituted in 0.1 N sodium hydroxide for protein quantification (Pierce BCA Kit, ThermoFisher Scientific, Waltham, MA).
## Chemogenetics

Mice were anaesthetized with isofluorane and placed in a stereotaxic apparatus (David Kopf Instruments). Recombinant AAV8-DIO-h3MDq-mCherry (Krashes et al., 2011) (titer = 3 × 10¹² genomes copies per ml, UNC Vector Core Services, Chapel Hill, NC) was bilaterally injected into the VTA/SN area (coordinates from bregma: AP: -3.44 mm, DV: -4.40 mm, L: ±0.48 mm) (Tsai et al., 2009) through a cannula (33 gauge, Plastics One) at a flow rate of 0.05 µl/min for 10 min (0.5 µl total volume per side) using a syringe pump (Harvard Apparatus). Mice were subsequently maintained in individual housing for at least 2 weeks prior to CNO treatment.

## Immunohistochemistry

Immunostaining was performed as previously described (Chu et al., 2013) using cryoprotected (30% Sucrose), free-floating coronal sections from fixed brains collected after intra-cardiac perfusion (4% paraformaldehyde) and cut at 40 microns using a cryostat (Leica, Solms, Germany). For fluorescent labeling, sections were incubated overnight with primary antibodies for mCherry (rabbit anti-RFP, 1:1000, Rockland, Limerick, PA) to enhance detection of the mCherry expression and antibodies for tyrosine hydroxylase (mouse anti-TH, 1:1000, EMD Millipore, Etibicoke, Canada). This was followed by 2 hr incubation with secondary antibodies (anti-rabbit Alexa Fluor 567 and anti-mouse Alexa Fluor 647, 1:250, Life Technologies, Carlsbad, CA) after which sections were mounted on superfrost slides (VWR, Radnor, PA), coverslipped with Vectashield mounting medium (Vector Labs, Burlington, Canada) and imaged by fluorescence microscopy (AxioObserver Z1, Zeiss, Jena, Germany). Quantification of co-expression was performed on a single, medial, midbrain section from each behaviorally responsive mouse using the cell counter plugin of ImageJ (National Institute of Health, Bethesda, MD) in order to determine the percentage of all TH-positive cells that also expressed the hM3Dq--mCherry fusion protein (N = 6 per genotype).

## Statistical analysis

One-way and repeated measures ANOVA, t-tests, and linear regression analyses were performed using Prism 5 (GraphPad, La Jolla, CA). Planned comparisons for the repeated measures were carried out to determine significant period changes between subsequent measurements. Bonferroni correction was used for post-hoc testing of individual group differences.

## False discovery rate (FDR) and cosinor determination

To evaluate the probabilistic significance of the least-square fit between the DA trace and the cosinor model that was computed based on the animal's locomotor activity oscillations, a false discovery rate approach was employed (Storch et al., 2007). For cosinor computation, we first determined the locomotor period of *Bmal1*~/*~ mice by Lomb-Scargle periodogram analysis of 20 hr of locomotor recordings following dialysate sampling. The determined period was then used as an input parameter for a least-squares cosinor analysis (Nelson et al., 1979) using a custom script and function written for Matlab (https://github.com/storchlab/Cosinor-FDR.git). This procedure generated a best-fit (co)sine wave modeling the activity time-series. This model wave was then used to assess the concordance between the extracellular DA fluctuations and locomotor oscillations. To calculate FDR, we randomly permutated the temporal order of the measured DA concentrations 100,000 times, and assessed its fit to the cosinor. The permutation procedure will degrade the fit of any truly rhythmic signal but will have little to no effect on noisy or flat profiles. The probability of false discovery was calculated as the percentage of randomly permuted traces that show an equal or better fit to the cosinor than the measured DA fluctuation.

## Butterworth filtering

Low pass filtering (1-hr cut-off) of raw data was conducted using a Butterworth zero-phase filter (Butterworth, 1930) in Matlab. This allowed for better visualization of ultradian frequencies in the >1-hr range, in a manner similar to recursive smoothing but without phase distortion.

## Period determination

Ultradian rhythms show variability in amplitude and period. We thus used Lomb-Scargle periodogram analysis (Clocklab) to estimate ultradian period length as this method is relatively tolerant to noisy data and data gaps, which can result from intermittent ultradian rhythm expression (Press and Rybicki, 1989). Unless otherwise stated, period length was determined by identifying the highest peak above the significance threshold (α = 0.01) in the Lomb-Scargle periodogram. Subsequent
plotting of an actogram at the determined modulus was conducted in each case in order to visually confirm rhythmicity and rule out the false identification of harmonic frequencies or side-lobes due to leakage (Van Dongen et al., 1999).

## Continuous wavelet transforms

To visualize and quantify the dynamics of period and amplitude in circadian incompetent mice, continuous wavelet transforms using the Morlet wavelet (Goupillaud et al., 1984) were performed for the 1–12-hr range using custom algorithms written in Matlab (https://github.com/storchlab/CWT.git). Locomotor activity data of seven consecutive days per individual SCN-lesioned or *Bmal1*⁻/⁻ mouse was used. For each spectrogram, a ridge identifying the peak amplitude of the spectrum was generated as described (Leise, 2013) and the associated values were used to calculate the average and standard deviation of the dominant ultradian frequency for each animal.

## Amplitude spectral density (ASD)

To determine the prevalence of oscillations in a given period range we calculated the area under the curve of all significantly rhythmic periodicities (α = 0.01) in the Lomb-Scargle periodogram, that is, the spectral density. For normalization, the obtained spectral density was divided by the total significant spectral density in the 0–40-hr range and expressed as percentage.

## Acknowledgements

We thank A Villemain, E Vigneault, and M-E Desaulnier for mouse-colony maintenance and F Tronche for the *DAT-Cre* mice. This work was supported by a NSERC discovery grant (K-FS), CIHR operating grants (K-FS and AG), a CIHR New Investigator Award (K-FS), by the Graham Boeckh Foundation (BG) and a Canada Research Chair (BG).

## Additional information

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Canadian Institutes of Health Research</td><td>Operating Grants</td><td>Alain Gratton, Kai-Florian Storch</td></tr><tr><td>Natural Sciences and Engineering Research Council of Canada</td><td>Discovery Grant</td><td>Kai-Florian Storch</td></tr><tr><td>Graham Boeckh Foundation</td><td>Graham Boeckh Chair in Schizophrenia</td><td>Bruno Giros</td></tr><tr><td>Canadian Institutes of Health Research</td><td>New Investigator Award</td><td>Kai-Florian Storch</td></tr><tr><td>Canadian Institutes of Health Research</td><td>Canada Research Chair</td><td>Bruno Giros</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

### Author contributions

IDB, K-FS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article; LZ, LM, Acquisition of data, Analysis and interpretation of data; MVK, AG, BG, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents

### Author ORCIDs

Ian D Blum, http://orcid.org/0000-0001-6999-8364

### Ethics

Animal experimentation: All experimental procedures were performed in accordance with the Canadian Council on Animal Care guidelines and approved by the McGill University Animal Care Committee (animal use protocol #2010-5945).
## References

Alexander GM, Rogan SC, Abbas Al, Armbruster BN, Pei Y, Allen JA, Nonneman RJ, Hartmann J, Moy SS, Nicolelis MA, McNamara JO, Roth BL. 2009. Remote control of neuronal activity in transgenic mice expressing evolved G protein-coupled receptors. *Neuron* **63**:27–39. doi: 10.1016/j.neuron.2009.06.014.

Aschoff J. 1965. Circadian rhythms in man. *Science* **148**:1427–1432. doi: 10.1126/science.148.3676.1427.

Aschoff J, Pohl H. 1978. Phase relations between a circadian rhythm and its zeitgeber within the range of entrainment. *Die Naturwissenschaften* **65**:80–84. doi: 10.1007/BF00440545.

Bengel D, Murphy DL, Andrews AM, Wichems CH, Feltner D, Heils A, Mössner R, Westphal H, Lesch KP. 1998. Altered brain serotonin homeostasis and locomotor insensitivity to 3, 4-methylenedioxymethamphetamine ('Ecstasy') in serotonin transporter-deficient mice. *Molecular Pharmacology* **53**:649–655.

Berk M, Dodd S, Kauer-Sant'anna M, Malhi GS, Bourin M, Kapczinski F, Norman T. 2007. Dopamine dysregulation syndrome: implications for a dopamine hypothesis of bipolar disorder. *Acta Psychiatrica Scandinavica* *Supplementum*:41–49.

Brown RE, Basheer R, McKenna JT, Strecker RE, McCarley RW. 2012. Control of sleep and wakefulness. *Physiological Reviews* **92**:1087–1187. doi: 10.1152/physrev.00032.2011.

Bunger MK, Wilsbacher LD, Moran SM, Clendenin C, Radcliffe LA, Hogenesch JB, Simon MC, Takahashi JS, Bradfield CA. 2000. Mop3 is an essential component of the master circadian pacemaker in mammals. *Cell* **103**:1009–1017. doi: 10.1016/S0092-8674(00)00205-1.

Butterworth S. 1930. On the theory of filter amplifiers. *Experimental Wireless* 7:536–541.

Chu A, Zhu L, Blum ID, Mai O, Leliavski A, Fahrenkrug J, Oster H, Boehm U, Storch KF. 2013. Global but not gonadotrope-specific disruption of Bmal1 abolishes the luteinizing hormone surge without affecting ovulation. *Endocrinology* **154**:2924–2935. doi: 10.1210/en.2013-1080.

Chung S, Lee EJ, Yun S, Choe D, Park SB, Son HJ, Kim KS, Dluzen DE, Lee I, Hwang O, Son GH, Kim K. 2014. Impact of circadian nuclear receptor REV-ERBα on midbrain dopamine production and mood regulation. *Cell* **157**:858–868. doi: 10.1016/j.cell.2014.03.039.

Daan S, Slopsema S. 1978. Short-term rhythms in foraging behaviour of the common vole, Microtus arvalis. *Journal of Comparative Physiology* **127**:215–227. doi: 10.1007/BF01350112.

de la Iglesia HO, Cambras T, Schwartz WJ, Diez-Noguera A. 2004. Forced desynchronization of dual circadian oscillators within the rat suprachiasmatic nucleus. *Current Biology* **14**:796–800. doi: 10.1016/j.cub.2004.04.034.

Domenger D, Dea D, Theroux L, Moquin L, Gratton A, Poirier J. 2012. The MPTP neurotoxic lesion model of Parkinson's disease activates the apolipoprotein E cascade in the mouse brain. *Experimental Neurology* **233**:513–522. doi: 10.1016/j.expneurol.2011.11.031.

Dowse HB, Hall JC, Ringo JM. 1987. Circadian and ultradian rhythms in period mutants of *Drosophila melanogaster*. *Behavior Genetics* **17**:19–35. doi: 10.1007/BF01066008.

Dowse H, Umemori J, Koide T. 2010. Ultradian components in the locomotor activity rhythms of the genetically normal mouse, Mus musculus. *The Journal of Experimental Biology* **213**:1788–1795. doi: 10.1242/jeb.038877.

Gainetdinov RR, Jones SR, Fumagalli F, Wightman RM, Caron MG. 1998. Re-evaluation of the role of the dopamine transporter in dopamine system homeostasis. *Brain Research Brain Research Reviews* **26**:148–153. doi: 10.1016/S0165-0173(97)00063-5.

Gainetdinov RR, Mohn AR, Caron MG. 2001. Genetic animal models: focus on schizophrenia. *Trends in Neurosciences* **24**:527–533. doi: 10.1016/S0166-2236(00)01886-5.

Gallardo CM, Darvas M, Oviatt M, Chang CH, Michalik M, Huddy TF, Meyer EE, Shuster SA, Aguayo A, Hill EM, Kiani K, Ikpeazu J, Martinez JS, Purpura M, Smit AN, Patton DF, Mistlberger RE, Palmiter RD, Steele AD. 2014. Dopamine receptor 1 neurons in the dorsal striatum regulate food anticipatory circadian activity rhythms in mice. *eLife* **3:e03781**. doi: 10.7554/eLife.03781.

Gann H, Riemann D, Hohagen F, Strauss LG, Dressing H, Müller WE, Berger M. 1993. 48-hour rapid cycling: results of psychopathometric, polysomnographic, PET imaging and neuro-endocrine longitudinal investigations in a single case. *Journal of Affective Disorders* **28**:133–140. doi: 10.1016/S0165-0327(93)90042-I.

Gerkema MP, Groos GA, Daan S. 1990. Differential elimination of circadian and ultradian rhythmicity by hypothalamic lesions in the common vole, Microtus arvalis. *Journal of Biological Rhythms* **5**:81–95. doi: 10.1177/074873049000500201.

Gerkema MP, van der Leest F. 1991. Ongoing ultradian activity rhythms in the common vole, Microtus arvalis, during deprivations of food, water and rest. *Journal of Comparative Physiology A, Sensory, Neural, and Behavioral Physiology* **168**:591–597. doi: 10.1007/BF00215081.

Gerkema MP, Verhulst S. 1990. Warning against an unseen predator: a functional aspect of synchronous feeding in the common vole, Microtus arvalis. *Animal Behaviour* **40**:1169–1178. doi: 10.1016/S0003-3472(05)80183-6.

Giros B, Jaber M, Jones SR, Wightman RM, Caron MG. 1996. Hyperlocomotion and indifference to cocaine and amphetamine in mice lacking the dopamine transporter. *Nature* **379**:606–612. doi: 10.1038/379606a0.

Goodwin JS, Larson GA, Swant J, Sen N, Javitch JA, Zahniser NR, De Felice LJ, Khoshbouei H. 2009. Amphetamine and methamphetamine differentially affect dopamine transporters in vitro and in vivo. *The Journal of Biological Chemistry* **284**:2978–2989. doi: 10.1074/jbc.M805298200.

Goupillaud P, Grossmann A, Morlet J. 1984. Cycle-octave and related transforms in seismic signal analysis. *Geoexploration* **23**:85–102. doi: 10.1016/0016-7142(84)90025-5.
Honma KI, Hiroshige T. 1978. Endogenous ultradian rhythms in rats exposed to prolonged continuous light. *The American Journal of Physiology* **235**:R250-R256.

Honma K, Honma S, Hiroshige T. 1986. Disorganization of the rat activity rhythm by chronic treatment with methamphetamine. *Physiology & Behavior* **38**:687-695. doi: 10.1016/0031-9384(86)90265-9.

Honma K, Honma S, Hiroshige T. 1987. Activity rhythms in the circadian domain appear in suprachiasmatic nuclei lesioned rats given methamphetamine. *Physiology & Behavior* **40**:767-774. doi: 10.1016/0031-9384(87)90281-2.

Hood S, Cassidy P, Cossette MP, Weigl Y, Verwey M, Robinson B, Stewart J, Amir S. 2010. Endogenous dopamine regulates the rhythm of expression of the clock protein PER2 in the rat dorsal striatum via daily activation of D2 dopamine receptors. *The Journal of Neuroscience* **30**:14046-14058. doi: 10.1523/JNEUROSCI.2128-10.2010.

Howell LL, Kimmel HL. 2008. Monoamine transporters and psychostimulant addiction. *Biochemical Pharmacology* **75**:196-217. doi: 10.1016/j.bcp.2007.08.003.

Ibuka N, Inouye SI, Kawamura H. 1977. Analysis of sleep-wakefulness rhythms in male rats after suprachiasmatic nucleus lesions and ocular enucleation. *Brain Research* **122**:33-47. doi: 10.1016/0006-8993(77)90660-6.

Ichikawa J, Meltzer HY. 1991. Differential effects of repeated treatment with haloperidol and clozapine on dopamine release and metabolism in the striatum and the nucleus accumbens. *The Journal of Pharmacology and Experimental Therapeutics* **256**:348-357.

Krashes MJ, Koda S, Ye C, Rogan SC, Adams AC, Cusher DS, Maratos-Flier E, Roth BL, Lowell BB. 2011. Rapid, reversible activation of AgRP neurons drives feeding behavior in mice. *The Journal of Clinical Investigation* **121**:1424-1428. doi: 10.1172/JCI46229.

Kraves S, Weitz CJ. 2006. A role for cardiotrophin-like cytokine in the circadian control of mammalian locomotor activity. *Nature Neuroscience* **9**:212-219. doi: 10.1038/nn1633.

Lane RF, Blaha CD. 1987. Chronic haloperidol decreases dopamine release in striatum and nucleus accumbens in vivo: depolarization block as a possible mechanism of action. *Brain Research Bulletin* **18**:135-138. doi: 10.1016/0361-9230(87)90042-6.

Leise TL. 2013. Wavelet analysis of circadian and ultradian behavioral rhythms. *Journal of Circadian Rhythms* **11**:5. doi: 10.1186/1740-3391-11-5.

Liu YY, Liu TY, Qu WM, Hong ZY, Urade Y, Huang ZL. 2012. Dopamine is involved in food-anticipatory activity in mice. *Journal of Biological Rhythms* **27**:398-409. doi: 10.1177/0748730412455913.

Lupinsky D, Moquin L, Gratton A. 2010. Interhemispheric regulation of the medial prefrontal cortical glutamate stress response in rats. *The Journal of Neuroscience* **30**:7624-7633. doi: 10.1523/JNEUROSCI.1187-10.2010.

Mistlberger RE. 2011. Neurobiology of food anticipatory circadian rhythms. *Physiology & Behavior* **104**:535-545. doi: 10.1016/j.physbeh.2011.04.015.

Mohawk JA, Baer ML, Menaker M. 2009. The methamphetamine-sensitive circadian oscillator does not employ canonical clock genes. *Proceedings of the National Academy of Sciences of USA* **106**:3519-3524. doi: 10.1073/pnas.0813366106.

Morice E, Denis C, Giros B, Nosten-Bertrand M. 2004. Phenotypic expression of the targeted null-mutation in the dopamine transporter gene varies as a function of the genetic background. *The European Journal of Neuroscience* **20**:120-126. doi: 10.1111/j.1460-9568.2004.03465.x.

Mukherjee S, Coque L, Cao JL, Kumar J, Chakravarty S, Asaithamby A, Graham A, Gordon E, Enwright JF III, DiLeone RJ, Birnbaum SG, Cooper DC, McClung CA. 2010. Knockdown of Clock in the ventral tegmental area through RNA interference results in a mixed state of mania and depression-like behavior. *Biological Psychiatry* **68**:503-511. doi: 10.1016/j.biopsych.2010.04.031.

Nelson W, Tong YL, Lee JK, Halberg F. 1979. Methods for cosinor-rhythmometry. *Chronobiologia* **6**:305-323.

Parmentier R, Ohtsu H, Djebbara-Hannas Z, Valatx JL, Watanabe T, Lin JS. 2002. Anatomical, physiological, and pharmacological characteristics of histidine decarboxylase knock-out mice: evidence for the role of brain histamine in behavioral and sleep-wake control. *The Journal of Neuroscience* **22**:7695-7711.

Paxinos G, Franklin KBJ. 2001. *The Mouse Brain in Stereotaxic Coordinates*. New York, NY: Academic Press.

Pendergast JS, Nakamura W, Friday RC, Hatanaka F, Takumi T, Yamazaki S. 2009. Robust food anticipatory activity in BMAL1-deficient mice. *PLOS ONE* **4**:e4860. doi: 10.1371/journal.pone.0004860.

Pitts S, Perone E, Silver R. 2003. Food-entrained circadian rhythms are sustained in arrhythmic Clk/Clk mutant mice. *American Journal of Physiology Regulatory, Integrative and Comparative Physiology* **285**:R57-R67.

Prendergast BJ, Cisse YM, Cable EJ, Zucker I. 2012. Dissociation of ultradian and circadian phenotypes in female and male Siberian hamsters. *Journal of Biological Rhythms* **27**:287-298. doi: 10.1177/0748730412448618.

Press WH, Rybicki GB. 1989. Fast algorithm for spectral analysis of unevenly sampled data. *The Astrophysical Journal* **338**:277-280. doi: 10.1086/167197.

Rivkees SA. 2003. Developing circadian rhythmicity in infants. *Pediatrics* **112**:373-381. doi: 10.1542/peds.112.2.373.

Ruis JF, Rietveld WJ, Buys PJ. 1987. Effects of suprachiasmatic nuclei lesions on circadian and ultradian rhythms in body-temperature in ocular enucleated rats. *Journal of Interdisciplinary Cycle Research* **18**:259-273. doi: 10.1080/09291018709359952.

Rusak B. 1977. The role of the suprachiasmatic nuclei in the generation of circadian rhythms in the golden hamster, Mesocricetus auratus. *Journal of Comparative Physiology* **118**:145-164. doi: 10.1007/BF00611819.

Scheving LE, Harrison WH, Gordon P, Pauly JE. 1968. Daily fluctuation (circadian and ultradian) in biogenic amines of the rat brain. *The American Journal of Physiology* **214**:166-173.

Schibler U. 2008. The mammalian circadian timekeeping system. In: Lloyd D, Rossi E, editors. *Ultradian Rhythms From Molecules to Mind*. Netherlands: Springer. p. 261-279.
Schmitz Y, Benoit-Marand M, Gonon F, Sulzer D. 2003. Presynaptic regulation of dopaminergic neurotransmission. *Journal of Neurochemistry* **87**:273-289. doi: 10.1046/j.1471-4159.2003.02050.x.

Seki Y, Tanimura T. 2014. Ultradian rhythm unmasked in the Pdf clock mutant of *Drosophila*. *Journal of Biosciences* **39**:585-594. doi: 10.1007/s12038-014-9450-z.

Smit AN, Patton DF, Michalik M, Opiol H, Mistlberger RE. 2013. Dopaminergic regulation of circadian food anticipatory activity rhythms in the rat. *PLOS ONE* **8**:e82381. doi: 10.1371/journal.pone.0082381.

Storch KF, Weitz CJ. 2009. Daily rhythms of food-anticipatory behavioral activity do not require the known circadian clock. *Proceedings of the National Academy of Sciences of USA* **106**:6808-6813. doi: 10.1073/pnas.0902063106.

Storch KF, Paz C, Signorovitch J, Raviola E, Pawlyk B, Li T, Weitz CJ. 2007. Intrinsic circadian clock of the mammalian retina: importance for retinal processing of visual information. *Cell* **130**:730-741. doi: 10.1016/j.cell.2007.06.045.

Strange PG. 2008. Antipsychotic drug action: antagonism, inverse agonism or partial agonism. *Trends in Pharmacological Sciences* **29**:314-321. doi: 10.1016/j.tips.2008.03.009.

Sulzer D. 2011. How addictive drugs disrupt presynaptic dopamine neurotransmission. *Neuron* **69**:628-649. doi: 10.1016/j.neuron.2011.02.010.

Tannenbaum GS, Martin JB. 1976. Evidence for an endogenous ultradian rhythm governing growth hormone secretion in the rat. *Endocrinology* **98**:562-570. doi: 10.1210/endo-98-3-562.

Tataroglu O, Davidson AJ, Benvenuto LJ, Menaker M. 2006. The methamphetamine-sensitive circadian oscillator (MASCO) in mice. *Journal of Biological Rhythms* **21**:185-194. doi: 10.1177/0748730406287529.

Thomas SA, Palmiter RD. 1997. Disruption of the dopamine beta-hydroxylase gene in mice suggests roles for norepinephrine in motor function, learning, and memory. *Behavioral Neuroscience* **111**:579-589. doi: 10.1037/0735-7044.111.3.579.

Tsai HC, Zhang F, Adamantidis A, Stuber GD, Bonci A, de Lecea L, Deisseroth K. 2009. Phasic firing in dopaminergic neurons is sufficient for behavioral conditioning. *Science* **324**:1080-1084. doi: 10.1126/science.1168878.

Turiault M, Parnaudeau S, Milet A, Parlato R, Rouzeau JD, Lazar M, Tronche F. 2007. Analysis of dopamine transporter gene expression pattern – generation of DAT-iCre transgenic mice. *The FEBS Journal* **274**: 3568-3577. doi: 10.1111/j.1742-4658.2007.05886.x.

van der Veen DR, Minh NL, Gos P, Arneric M, Gerkema MP, Schibler U. 2006. Impact of behavior on central and peripheral circadian clocks in the common vole Microtus arvalis, a mammal with ultradian rhythms. *Proceedings of the National Academy of Sciences of USA* **103**:3393-3398. doi: 10.1073/pnas.0507825103.

Van Dongen HP, Olofsen E, Van Hartevelt JH, Kruyt EW. 1999. A procedure of multiple period searching in unequally spaced time-series with the Lomb-Scargle method. *Biological Rhythm Research* **30**:149-177. doi: 10.1076/brhm.30.2.149.1424.

van Oort BE, Tyler NJ, Gerkema MP, Folkow L, Stokkan KA. 2007. Where clocks are redundant: weak circadian mechanisms in reindeer living under polar photic conditions. *Die Naturwissenschaften* **94**:183-194. doi: 10.1007/s00114-006-0174-2.

Vitaterna MH, King DP, Chang AM, Kornhauser JM, Lowrey PL, McDonald JD, Dove WF, Pinto LH, Turek FW, Takahashi JS. 1994. Mutagenesis and mapping of a mouse gene, Clock, essential for circadian behavior. *Science* **264**:719-725. doi: 10.1126/science.8171325.

Wang S, Tan Y, Zhang JE, Luo M. 2013. Pharmacogenetic activation of midbrain dopaminergic neurons induces hyperactivity. *Neuroscience Bulletin* **29**:517-524. doi: 10.1007/s12264-013-1327-x.

Wehr TA, Goodwin FK, Wirz-Justice A, Breitmaier J, Craig C. 1982. 48-hour sleep-wake cycles in manic-depressive illness: naturalistic observations and sleep deprivation experiments. *Archives of General Psychiatry* **39**:559-565. doi: 10.1001/archpsyc.1982.04290050037008.

Wehr TA, Turner EH, Shimada JM, Lowe CH, Barker C, Leibenluft E. 1998. Treatment of a rapidly cycling bipolar patient by using extended bed rest and darkness to stabilize the timing and duration of sleep. *Biological Psychiatry* **43**:822-828. doi: 10.1016/S0006-3223(97)00542-8.

Welsh D, Richardson GS, Dement WC. 1988. Effect of running wheel availability on circadian patterns of sleep and wakefulness in mice. *Physiology & Behavior* **43**:771-777. doi: 10.1016/0031-9384(88)90375-7.

Wever RA. 1979. In: Schaefer KE, editor. *The Circadian System of Man: Results of Experiments Under Temporal Isolation*. New York: Springer-Verlag Berlin Heidelberg.

White FJ, Wang RY. 1983. Comparison of the effects of chronic haloperidol treatment on A9 and A10 dopamine neurons in the rat. *Life Sciences* **32**:983-993. doi: 10.1016/0024-3205(83)90929-3.

Wilk K, Hegerl U. 2010. Time of mood switches in ultra-rapid cycling disorder: a brief review. *Psychiatry Research* **180**:1-4. doi: 10.1016/j.psychres.2009.08.011.

Wirz-Justice A, Bromundt V, Cajochen C. 2009. Circadian disruption and psychiatric disorders: the importance of entrainment. *Sleep Medicine Clinics* **4**:273-284. doi: 10.1016/j.jsmc.2009.01.008.

Wirz-Justice A, Cajochen C, Nussbaum P. 1997. A schizophrenic patient with an arrhythmic circadian rest-activity cycle. *Psychiatry Research* **73**:83-90. doi: 10.1016/S0165-1781(97)00117-0.

Wirz-Justice A, Haug HJ, Cajochen C. 2001. Disturbed circadian rest-activity cycles in schizophrenia patients: an effect of drugs? *Schizophrenia Bulletin* **27**:497-502. doi: 10.1093/oxfordjournals.schbul.a006890.

Wisor JP, Nishino S, Sora I, Uhl GH, Mignot E, Edgar DM. 2001. Dopaminergic role in stimulant-induced wakefulness. *The Journal of Neuroscience* **21**:1787-1794.

Wulff K, Dijk DJ, Middleton B, Foster RG, Joyce EM. 2012. Sleep and circadian rhythm disruption in schizophrenia. *The British Journal of Psychiatry* **200**:308-316. doi: 10.1192/bjp.bp.111.096321.
Xu F, Gainetdinov RR, Wetsel WC, Jones SR, Bohn LM, Miller GW. Wang YM, Caron MG. 2000. Mice lacking the norepinephrine transporter are supersensitive to psychostimulants. *Nature Neuroscience* **3**:465-471.
doi: 10.1038/74839.

Zhao ZQ, Scott M, Chiechio S, Wang JS, Renner KJ, RW Gereau IV, Johnson RL, Deneris ES, Chen ZF. 2006. Lmx1b is required for maintenance of central serotonergic neurons and mice lacking central serotonergic system exhibit normal locomotor activity. *The Journal of Neuroscience* **26**:12781-12788. doi: 10.1523/JNEUROSCI.4143-06.2006.

Zhou QY, Palmiter RD. 1995. Dopamine-deficient mice are severely hypoactive, adipsic, and aphagic. *Cell* **83**:1197-1209. doi: 10.1016/0092-8674(95)90145-0.