Diurnal influences on electrophysiological oscillations and coupling in the dorsal striatum and cerebellar cortex of the anesthetized rat

Circadian rhythms modulate behavioral processes over a 24 h period through clock gene expression. What is largely unknown is how these molecular influences shape neural activity in different brain areas. The clock gene Per2 is rhythmically expressed in the striatum and the cerebellum and its expression is linked with daily fluctuations in extracellular dopamine levels and D2 receptor activity. Electrophysiologically, dopamine depletion enhances striatal local field potential (LFP) oscillations. We investigated if LFP oscillations and synchrony were influenced by time of day, potentially via dopamine mechanisms. To assess the presence of a diurnal effect, oscillatory power and coherence were examined in the striatum and cerebellum of rats under urethane anesthesia at four different times of day zeitgeber time (ZT1, 7, 13 and 19—indicating number of hours after lights turned on in a 12:12 h light-dark cycle). We also investigated the diurnal response to systemic raclopride, a D2 receptor antagonist. Time of day affected the proportion of LFP oscillations within the 0–3 Hz band and the 3–8 Hz band. In both the striatum and the cerebellum, slow oscillations were strongest at ZT1 and weakest at ZT13. A 3–8 Hz oscillation was present when the slow oscillation was lowest, with peak 3–8 Hz activity occurring at ZT13. Raclopride enhanced the slow oscillations, and had the greatest effect at ZT13. Within the striatum and with the cerebellum, 0–3 Hz coherence was greatest at ZT1, when the slow oscillations were strongest. Coherence was also affected the most by raclopride at ZT13. Our results suggest that neural oscillations in the cerebellum and striatum, and the synchrony between these areas, are modulated by time of day, and that these changes are influenced by dopamine manipulation. This may provide insight into how circadian gene transcription patterns influence network electrophysiology. Future experiments will address how these network alterations are linked with behavior.


INTRODUCTION
Local field potential (LFP) oscillations are recorded when rhythmic fluctuations in membrane potentials synchronize within a nearby group of cells. These oscillations can provide a mechanism for the storage and transfer of information by coordinating the neuronal activity between distant networks (Başar et al., 2001;Buehlmann and Deco, 2010;Uhlhaas et al., 2010). The many ion channels involved in maintaining the membrane potential and in generating post-synaptic potentials, along with neurotransmitter modulation, influence the dominant oscillatory parameters within neural networks (Hutcheon and Yarom, 2000;Buzsáki and Draguhn, 2004). In the suprachiasmatic nucleus (SCN), a strong influence on many of these cellular processes is the circadian expression of proteins that regulate neuron function and local neural network connections (Herzog, 2007;Colwell, 2011). Circadian protein expression shows distinct phase relationships throughout widespread areas of the brain, including the striatum and cerebellum (Namihira et al., 1999;Rath et al., 2012;Harbour et al., 2013), however, little is known about the role of circadian protein expression in regulating neuron function in these areas. In this paper, we examine how striatal and cerebellar networks may be under the influence of circadian rhythms. The role of the cerebellum in the context of circadian rhythms has scarcely been studied, however it has been shown to share a role with the dorsal striatum in regulating locomotor activity patterns under restricted feeding, a behavior that is modified by circadian genes (Mendoza et al., 2010;Verwey and Amir, 2012). Furthermore, there is an interest in the subcortical connectivity between these two areas and their cooperative contributions to motor control (Middleton and Strick, 2000;. In order to understand the neurophysiological changes associated with the circadian modulation influencing the striatum and cerebellum, we simultaneously recorded LFPs in these two areas and examined concurrent oscillations and their coherence. Within the striatal networks, and the granule cell layer of the posterior lobe of the cerebellum of rats and primates, LFP oscillations around 5-30 Hz are present during a variety of tasks (Courtemanche et al., 2002Berke et al., 2004;DeCoteau et al., 2007b). These oscillations, and their synchrony with more distant brain areas, are modulated under behavioral conditions such as sensorimotor associations, learning, reward expectancy and sleep (Destexhe et al., 1999;Steriade, 2003;Buzsáki, 2006;Thorn and Graybiel, 2014). Anesthetics, such as urethane, shape the LFPs of cortical and subcortical structures into predominant slow wave activity, inducing a sleep-like state in the LFP that shows widespread synchrony across many brain areas (Magill et al., 2004;Clement et al., 2008;Ros et al., 2009;Sharma et al., 2010). In this study, we used this oscillation-permissive state to evaluate the circadian expression of oscillatory network activity.
Both under anesthesia and in awake behaving animals, dopamine changes can induce important modifications to striatal oscillations (Walters et al., 2007;Lemaire et al., 2012). Dopamine depletion causes an increase in the power of oscillations below 30 Hz and an increase in coherence of these oscillations with the cerebral cortex and throughout the basal ganglia (Brown et al., 2001;Sharott et al., 2005). This is accompanied by a decrease in power of the oscillations above 70 Hz. Results from animal models that use dopamine antagonists or midbrain dopaminergic lesions are consistent with the findings in Parkinson's disease in humans, where abnormal oscillations can be normalized with dopamine agonists or dopamine-related treatments such as L-DOPA (Brown et al., 2001;Rivlin-Etzion et al., 2006;Ballion et al., 2009), further demonstrating the specific role of dopamine in modulating striatal oscillations.
The dopaminergic system is under the influence of circadian rhythmicity. Endogenous extracellular dopamine levels show cyclic fluctuations throughout the day and peak during the middle of the dark phase in rodents (Owasoyo et al., 1979;Hood et al., 2010;Ferris et al., 2014). Furthermore, activation of D2 dopamine receptors has been linked with normal expression of Per2, one of the core circadian genes, in the rat dorsal striatum (Hood et al., 2010;Gravotta et al., 2011). Per2 is involved in regulating daily physiological and behavioral cycles and is found in most tissues throughout the body (Albrecht et al., 2007). However, it is unclear how its rhythmic expression leads to such behavioral changes. Vesicle transporter activity throughout the brain (Darna et al., 2009) and changes in hippocampal synaptic plasticity (Chaudhury et al., 2005) exhibit diurnal variations; similarly, modulation of neurotransmitter function by the circadian clock provide potential ways through which circadian molecular mechanisms may be translated into diurnal electrophysiological and behavioral changes.
A limited number of studies have assessed the electrophysiological correlates of the circadian cycle downstream of the SCN, the master circadian clock, or its direct connections (Guilding and Piggins, 2007). Mordel et al. (2013) assessed diurnal variations of in vitro Purkinje cell firing and spontaneous inhibitory postsynaptic potentials, but they found that these measures were not modulated in a circadian manner. Restrictive feeding influences the rhythm of Per2 expression in the Purkinje and granule cell layers (Mendoza et al., 2010). There is limited evidence that dopamine levels are also rhythmic in the cerebellum (Owasoyo et al., 1979) but little more is known about circadian modulation of cerebellar physiology; including the functional role of clock genes or their regulatory mechanisms. Clock gene expression may help to organize the activity of striatal and cerebellar networks through mechanisms that affect the generation and synchronization of LFP oscillations. Such oscillations, recorded under urethane anesthesia, can provide information about the local activity and the neural synchronization between distant brain structures.
The purpose of this study was to determine if there are diurnal changes in LFP oscillations and coherence in the rat dorsal striatum and cerebellar cortex. Furthermore, since Per2 rhythms in the striatum are responsive to D2 receptor antagonists and not D1 antagonists (Hood et al., 2010), we also explored the effects of raclopride, a D2 antagonist, at different times of day. Given the endogenous rhythm of dopamine in the striatum, we hypothesized that circadian electrophysiological correlates would be expressed through diurnal changes. These changes could be seen in the potency of the LFP oscillations in each area, and in their synchrony. The basal ganglia and cerebellum show a cooperative role in sensorimotor control (Doya, 2000), and in their regulation of locomotor activity patterns related to circadian timing. These areas can interact through the cerebral cortex (Middleton and Strick, 2000) and via subcortical pathways (Ichinohe et al., 2000;Hoshi et al., 2005;. In order to identify their network interactions, we thus measured the LFP coherence between distant electrodes (DeCoteau et al., 2007a;Thorn et al., 2010;Frederick et al., 2013). The variables thus sought were the diurnal effects in (1) the power of LFP oscillations; (2) the LFP coherence within and across the striatum and cerebellar cortex; and also (3) if raclopride would induce an alteration in the circadian pattern of oscillations or coherence.

ANIMALS
Sixteen Sprague-Dawley male rats (Charles River, St-Constant, QC) 300-400 g in size were used for the experiment. Rats were housed separately in cages equipped with running wheels and kept in individual lightproof and sound-attenuated cabinets. Prior to electrophysiological recordings, rats were maintained on a 12:12 h light-dark cycle for a minimum of 2 weeks. Rats were assigned randomly to groups designated by four different target experimental times, corresponding to zeitgeber time (ZT) 1, 7, 13 and 19, where ZT0 represents when lights were turned on and ZT12 when lights were turned off. Light cycles were Frontiers in Systems Neuroscience www.frontiersin.org September 2014 | Volume 8 | Article 145 | 2 adjusted for each animal so that the target experimental time occurred at 1:30 pm for all animals. Environmental conditions were maintained at a constant temperature and rats were provided ad libitum access to tap water and laboratory chow throughout the entire entrainment period. To confirm that rats had adjusted to their light cycles, wheel-running activity was monitored continuously using the VitalView software (Mini Mitter, Bend, OR). All experimental procedures followed the guidelines of the Canadian Council on Animal Care and were approved by the Concordia University Animal Research Ethics Committee.

DRUG TREATMENT AND RECORDINGS
Surgical preparation for each animal was scheduled such that the animal's targeted circadian time occurred mid-way through the electrophysiological recording session. LFP recordings were obtained both before and after administration of raclopride. At the start of the recording procedures (ZT time − ∼4 h) the rats were anesthetized with a 5% isoflurane and 95% oxygen mixture until placement of a catheter into the jugular vein. They were then transferred to urethane anesthesia (0.8 g/ml), introduced through the jugular catheter in 1.0 ml boluses as needed. Verification of depth of anesthesia was monitored throughout the experiment by lack of a foot-pinch reflex. In 11 animals, transcardial heartbeat signals were also recorded along with LFPs using subdermal electrodes inserted at the level of the axillar forelimb fossa. Rats were placed in a stereotaxic apparatus and body temperature was maintained at 37 • C, monitored by a rectal thermometer. Throughout the experiment, the eyes of the animal were covered to prevent light exposure from influencing changes in the SCN. Four pairs of tungsten microelectrodes (1.0-1.2 MΩ; FHC, Bowdoin, ME) were used, and the difference in depth within each pair of electrodes was set to 0.5 mm. Electrodes were attached to microdrives, held by 2 modified head stages one for the striatal electrodes and one for the cerebellar electrodes, which were anchored to the stereotaxic apparatus by custom built arms. Two pairs of striatal electrodes were used, one placed in the medial dorsal striatum [0.48 mm AP, 1.40 mm ML, 4.70 mm DV; (Paxinos and Watson, 1998)] and another placed in the lateral dorsal striatum (0.48 mm AP, 5.40 mm ML, 5.00 mm DV). Two pairs of electrodes were also inserted into the cerebellar cortex, separated by approximately 2 mm, in the contralateral Crus2 or Paramedian lobule of the cerebellar cortex. Access to Crus2/Paramedian lobule was gained through a craniotomy roughly 2.5 mm in diameter in the occipital bone, located 2-3 mm lateral from the midline. Crus2 was then identified visually and electrodes were inserted at a 45 • angle to a depth of about 0.5 mm. The electrode position was adjusted until the granule cell layer was identified by characteristic multiunit activity and LFP oscillatory activity was identified on an oscilloscope. The 0.5 mm distance between electrodes within a pair suggests that both electrodes were not necessarily in the granule cell layer; this distance allowed for a comparison of nearby circuit activity with a strong granule cell layer influence, given the density of the granule cell layer and the electrode impedance. The nearby comparison could then be contrasted with the 2 mm distance of the second pair. All signals were grounded to stainless steel screws anchored to the skull on the contralateral side and referenced to a 16 gauge circular steel wire placed either in the cerebral cortex approximately 5 mm posterior to the striatal electrodes or on the unrecorded side of the cerebellum.
Electrophysiological signals were collected with a Neuralynx Cheetah system (Neuralynx, Bozeman, MT) at a sampling rate of 30,300 Hz and band-pass filtered between 0.1 and 9000 Hz. LFPs were then extracted from this signal with a 300-Hz lowpass filter. A typical recording session began with a first set of recordings, including 5 samples of 120 s separated by an interval of 90 s without recording, following which raclopride tartrate (Sigma-Aldrich, Oakville, ON) was administered via intraperitoneal injection (1 mg/kg dissolved in distilled water at 1 µg/ml). After a 15 min waiting period, 5 additional 120 s samples were recorded. Raclopride was administered at a dosage shown to induce akinesia and prevent operant conditioning, and to have effects on neural activity in the striatum as early as 5 min following intraperitoneal injection (Fowler and Liou, 1998;Dejean et al., 2011). Even though the animals were heated between recordings, a drop in body temperature sometimes occurred in the later recording periods; therefore, the second recording periods both before and after raclopride administration were used for analysis, since these showed stability across all subjects.

TISSUE STAINING
Following recordings, the brain was lesioned using a lesion maker (Grass Products, Warwick, RI) to mark electrode placement (300 µA for 30 s). Rats were then perfused with 300 ml cold 0.9% saline followed by 300 ml cold 4% paraformaldehyde in 0.1M phosphate buffer. Brains were extracted and kept overnight at 4 • C in a 4% paraformaldehyde solution then dehydrated in phosphate buffered saline (PBS) with 30% sucrose. The cerebrum was sliced coronally and the cerebellum was sliced horizontally into 50 µm sections on a cryostat and mounted directly onto slides. Sections were stained using cresyl violet and lesions were located using light microscopy.

ELECTROPHYSIOLOGICAL DATA PROCESSING
Raw microelectrode data was downsampled to 2020 Hz using Neuralynx software and off-line analysis of recorded data was conducted using custom Matlab (MathWorks, Natick, MA) functions and scripts. Once in Matlab, the data was low-pass filtered at 300 Hz and structured into 2-s windows. To analyze the rhythmicity content of the low-passed signal, Fast Fourier Transform (FFT) was calculated on each 2-s window. The following bands were used for analysis: (1) (0.1-3 Hz, delta); (2) (3-8 Hz, theta); (3) (8-13 Hz, alpha); (4) (13-30 Hz, beta); and (5) (30-55 Hz, gamma). The power spectrum was integrated within those bands, and then averaged for the 60 windows, providing one value per recording period. Coherence between LFP pairs was also calculated on each 2-s window and values within the same frequency bands were also averaged per recording period. The average FFT and coherence of each period was entered into the statistical analyses.
In order to establish the percentage of time where slow (0-3 Hz) vs. 3-8 Hz oscillations predominate, the integrated FFT power for each band was compared for each 2-s time window.
Plotting those values in time for each recording period, we identified the power level around which the dominant power Frontiers in Systems Neuroscience www.frontiersin.org September 2014 | Volume 8 | Article 145 | 3 between 0-3 Hz and 3-8 Hz was alternating. To compute this, we calculated the difference between the 0-3 and 3-8 Hz integrated power and took the mean of the values where this difference equaled 0 ± 10 of the total percent FFT power. This difference maintained a narrow range while still including enough data points to normalize, thus providing a baseline for each recording session, from which a threshold could be calculated. This mean +20% value determined if a window had strong 0-3 Hz activity, while the mean -15% determined if a window had strong 3-8 Hz activity. These thresholds were chosen empirically as the most conservative levels where the 0-3 Hz or 3-8 Hz oscillations could be detected across multiple recording sessions. The number of windows that fell above the threshold was calculated for each band, and then divided by the total number of windows to get a value expressed as a percentage.
To avoid duplication of data, spindles were analyzed for one electrode in each of the MStr and LStr pairs of electrodes and for one electrode placed in the cerebellum. The average of the number of spindles detected from the two striatal electrodes was then used for statistical analysis. A detectable spindle was defined as an oscillation within 8-20 Hz that occurred at the peak of the slow wave. Spindles were identified by first band-pass filtering the LFP signal between 0-3 Hz to identify the slow oscillation phase. The LFP was also filtered between 8-20 Hz, and this filtered signal was rectified, and convolved with a 14-Hz (71 ms) window. These 8-20 Hz filtered and rectified product-values were Z-scored and a threshold was set at 2 standard deviations. Spindles were detected when the 8-20 Hz signal reached above this threshold value within ±0.25 s of the peak of a slow wave oscillation. An example is shown in Figure 5A.

NEOCORTICAL RECORDINGS
In order to verify if coherence between the striatum and cerebellar cortex under urethane anesthesia involves neocortical connections, four additional recording sessions were done with neocortical electrodes added, one at each time of day (ZT1, 7, 13 and 19). The same electrophysiological set-up was used as described above, but instead of placing two pairs of electrodes in the striatum, one electrode from each pair was stopped in the overlying neocortex, corresponding to the sensorimotor cortex. These recordings were pooled together with the rest of the dataset to provide an overall description of synchrony between the neocortex, the striatum and the cerebellum. Coherence within each frequency band was analyzed as above, however the trend was similar across all bands, so we report below the 0.1-8 Hz and overall coherence from 0.1 to 55 Hz.

STATISTICAL ANALYSIS
Statistical tests were performed using Statistica 10 (StatSoft, Tulsa, OK) or in Matlab. All results are presented as the mean ± the standard error of the mean (SEM). The second recording period before and after raclopride showed the most stability with body temperature; so all further analysis was performed on these two periods. The dependent variables of oscillatory power and coherence (for each frequency band) were statistically assessed for two different analyses. Circadian effect on baseline recordings was assessed by one-way repeated measures ANOVA, with time-of-day as the independent variable. The circadian effects of raclopride assessed the interaction between pre-post raclopride conditions and time-of day using two-way repeated-measures ANOVAs. Results were considered significant at p < 0.05. Post-hoc comparisons used Tukey tests. For identifying the relationship between 0-3 Hz and 3-8 Hz power across windows, linear correlation (p < 0.05) was used.

IN VIVO RECORDINGS
LFP recordings were obtained from 16 urethane-anesthetized animals (4 animals × 4 ZT). There was no effect of either time-of-day or pre-post raclopride on the heart rate (ZT, F (3,7) = 0.14, p = 0.93, pre-post, F (1,7) = 0.14, p = 0.72). This confirms that at least from the heart rate marker, the depth of anesthesia was consistent across ZT and pre-post conditions. One additional (17th) animal from the ZT13 group failed to maintain normal body temperature during the experiment and was excluded from the analysis. Each of the remaining recordings included simultaneous medial and lateral striatal (MStr and LStr) recordings and cerebellar Crus2/Paramedian lobule (CB) recordings. As explained above, in 4 of these animals, one electrode from each pair aimed at the MStr and LStr was instead placed in the overlying neocortex. Electrode locations determined from electrolytic lesion sites confirmed the location of microelectrode tips in the striatum and granule cell layer of the cerebellum or close vicinity. Placement of 17 of the 32 cerebellar pairs, 10 of the 16 MStr and LStr pairs and all of the neocortical electrode tips are shown in Figure 1, providing a representative spread of recording location for each brain area.
Both the striatal and cerebellar LFPs showed predominant large amplitude slow wave activity between 0.8-2 Hz with intermittent smaller amplitude oscillations around 4-5 Hz (Figure 2A). The slow wave activity is consistent with other findings in which urethane anesthesia induces widespread cortical synchrony resulting in oscillations at ∼1 Hz that are intermittently interrupted by a desynchronized state (Destexhe et al., 1999;Clement et al., 2008). Slow wave activity tended to co-occur in the striatum and cerebellum, consistent with findings that anesthesia induces coherent slow-wave oscillations between the neocortex and striatum (Tseng et al., 2001;Mallet et al., 2005) and between the neocortex and cerebellum (Ros et al., 2009). We also observed transient spindle-like 8-20 Hz oscillations that were nested within the slow-wave activity. These spindles were most evident in the striatal LFPs, but similar electrophysiological events also occurred in the cerebellar LFPs ( Figure 2B).

DIURNAL MODULATION OF LFP OSCILLATIONS
To determine the effects of time of day on both slow-wave and 3-8 Hz oscillations, the power of the FFT within each band recorded at baseline, was compared across the four times of day ( Figure 3A). The MStr and LStr recordings showed the same variations by time of day and were thus combined for statistical analysis. In both striatal and cerebellar LFPs, slow-wave power showed a large diurnal rhythm in which power was the greatest at ZT1 and the lowest at ZT13 (Str, F (3,52) = 11.80, p < 0.001; CB, F (3,60) = 14.75, p < 0.001). There was an inverse diurnal pattern in the power of 3-8 Hz oscillations, in which power was highest at ZT13 (Str, F (3,52) = 5.78, p < 0.01; CB, F (3,60) = 11.56, p < 0.001).
No diurnal effect was found in frequency bands above 8 Hz for either region (not shown). Because the overall FFT power at 0-3 Hz and 3-8 Hz were inversely related to each other by time of day, we also sought to determine the relationship between the two frequency bands within each 2-s window, and to quantify differences in time spent in slow oscillations vs. 3-8 Hz oscillations. The relationship between the power of slow waves vs. 3-8 Hz oscillations was assessed by plotting the percent power of the FFT within the 0-3 Hz band vs. the 3-8 Hz band for each of the 180 2-s windows sampled in baseline recordings ( Figure 3B). There was a strong negative correlation between power in the two bands for all recording sessions ( Figure 3B shows two examples of this negative correlation). The inverse relationship of the power at 0-3 Hz vs. 3-8 Hz within each 2-s window indicated that the two oscillatory modes did not co-occur. This was then used to create the threshold to detect the duration of 0-3 Hz or 3-8 Hz activity within each recording session. Since no diurnal effect was observed in bands above 8 Hz, these bands were excluded from this detection analysis.

EFFECTS OF RACLOPRIDE ON DIURNAL VARIATIONS IN THE OCCURRENCE OF LFP OSCILLATIONS
Diurnal rhythms in the power of the LFP were highly associated to the proportion of time spent within each oscillatory mode, and raclopride influenced this distribution in a time-of-daydependent manner (Figure 4). The percentage of time spent preferentially in slow wave oscillations vs. 3-8 Hz oscillations was calculated for each recording period, and both cerebellar and striatal LFPs showed robust, inverse diurnal rhythms for the proportion of time spent in each band (see Figures 4C,D, closed symbols; 0-3 Hz: Str, F (3,52) = 6.63, p < 0.001, CB, F (3,60) = 11.85, p < 0.001; 3-8 Hz: Str, F (3,52) = 7.78, p < 0.001, CB, F (3,60) = 13.86, p < 0.001). At ZT1, a large proportion of time was spent in slow wave oscillations (86.9 ± 5.9% of the time for cerebellum and 94.2 ± 5.7% for striatum) and at ZT13 there was a much greater prevalence of 3-8 Hz oscillations in both the striatum and cerebellar cortex (57.7 ± 5.2% of the time for cerebellum and 40.0 ± 4.3% for striatum). Overall, both cerebellar and striatal LFPs showed similar diurnal effects in the proportion of time spent in each of the two frequency bands.
Since circadian PER2 rhythms in the striatum are dependent upon dopamine and D2 receptor activation (Hood et al., 2010), we assessed the effects of raclopride on LFP oscillations as a function of time of day. For both striatal and cerebellar LFPs, systemic raclopride administration resulted in an increase in slow oscillations a decrease in in 3-8 Hz oscillations (Figures 4A,B). This effect was only significant at ZT13, except the 3-8 Hz band in the striatum, where only a trend was observed (Figures 4C,D; 0-3 Hz: Str, F (3,52) = 3.00, p < 0.05; CB, F (3,60) = 4.67, p < 0.01; 3-8 Hz: Str, F (3,52) = 2.37, p = 0.08, CB, F (3,60) = 11.91, p < 0.001). The raclopride-induced increase in the occurrence of 0-3 Hz oscillations at ZT13 was also mirrored in the overall power of the FFT (not shown). Typically, the changes in oscillatory activity were complementary between the bands, such that raclopride induced an approximate 15% increase in time spent in the 0-3 Hz band and an approximate 15% decrease in time spent in the 3-8 Hz band for both the striatal and cerebellar LFPs.

DIURNAL VARIATION IN THE INCIDENCE OF SPINDLES
Spindles occur during slow wave oscillations during sleep and under urethane anesthesia, and can be modulated by distinct thalamic mechanisms (Valencia et al., 2013). We sought to determine if the incidence of spindles varies in a diurnal manner. The incidence of spindles was determined for striatal and cerebellar recording sites for each animal, and then compared between ZT and before and after raclopride administration. The occurrence of spindles in the LFP showed a similar pattern to that of slow waves, with a significant diurnal effect for striatal (F (3,28) = 4.93, p < 0.01), but not cerebellar recordings (F (3,12) = 2.12, p = 0.15) ( Figure 5B). The diurnal pattern of spindles was dependent upon the occurrence of slow waves, such that when normalized to the amount of time spent in 0-3 Hz oscillations there was no diurnal variation in spindle incidence (Str, F (3,28) = 1.16, p = 0.34; CB, F (3,12) = 1.54, p = 0.25) ( Figure 5C). In addition, raclopride did not have an effect on the incidence of spindles (Str, F (3,28) = 1.17, p = 0.34; CB, F (3,12) = 1.88, p = 0.20; data not shown). Therefore, diurnal variations in the incidence of spindles are dependent on the presence of slow wave oscillations under our urethane anesthesia preparation, and are not significantly affected by D2 receptor antagonism.

EFFECTS OF RACLOPRIDE ON DIURNAL VARIATIONS IN COHERENCE
To assess the overall interregional differences in the synchrony of LFP oscillations recorded in the striatum and cerebellum, coherence within and between these structures was evaluated. Within the 0-3 Hz and 3-8 Hz bands, coherence showed similar levels and were therefore combined. There was also no difference between the following comparisons, so this data was combined: MStr and LStr within pair, the MStr-CB and LStr-CB, the two cortical locations with the CB, and the two off-set cortico-striatal comparisons. 0-8 Hz coherence followed a similar pattern as 8-55 Hz coherence, but was consistently higher (see Figure 6A, 0-8 Hz vs. 0-55 Hz). Coherence was higher for recordings obtained from adjacent tips within each of the striatal electrode pairs (Within Pair) as compared to recordings obtained between electrode pairs (Btwn Pair) in the MStr vs. LStr (p < 0.01) or across electrode pairs in the cerebellum (p < 0.001). This demonstrates enhanced coherence of LFPs with spatial proximity in both the striatum and the cerebellar cortex ( Figure 6A, Within Pair vs. Btwn Pair comparisons). Four additional recording sessions were performed to evaluate the contribution of the neocortex to the striato-cerebellar coherence and included simultaneous recordings from two striatal electrodes, two electrodes in the overlying neocortex, and four electrodes in the cerebellar cortex. Corticostriatal coherence was greatest between electrodes in the same sagittal plane and was comparable to the striatal-within pair comparisons ( Figure 6A, 0-8 Hz: Str-within pair = 0.92 ± 0.024; Ctx-Str vertical = 0.90 ± 0.042). Cortico-striatal coherence was reduced for comparisons in different medial vs. lateral positions, showing a similar coherence as the striatal-between pair comparisons (Str-between pair = 0.82 ± 0.016; Ctx-Str off-set = 0.79 ± 0.042). The levels of Str-CB or Ctx-CB coherence (Str-CB = 0.53 ± 0.008; Ctx-CB = 0.54 ± 0.021) were substantially lower than levels of coherence obtained from more nearby electrodes pairs within the striatum or cerebellum. However, levels of both Str-CB and Ctx-CB coherence had similar values. This is consistent with the neocortex playing a role in coordinating oscillations in the striatum and cerebellum under urethane anesthesia. The presence of diurnal effects on mean coherence for each recording session was assessed for the 0-3 Hz and 3-8 Hz bands within striatal, cerebellar and striato-cerebellar comparisons. CB within pair and between pair along with CB-MStr and CB-LStr comparisons showed the same modulation by time-of-day and pre-post raclopride, so these data were combined. Coherence of striatal LFPs showed a significant main effect for time-ofday for slow oscillations between medial and lateral recording sites (Str-between pair F (3,48) = 7.24, p < 0.001) but not with closer electrode pairs (Str-within pair,  Figure 4C). These comparisons also showed a secondary increase at ZT1, suggesting a biphasic pattern of coherence in the 3-8 Hz band.
To assess the effect of D2 receptor antagonism on striatal and cerebellar coherence, the mean coherence before and after raclopride administration was compared for each of the four times of day. Raclopride had the effect of increasing coherence of slow oscillations at ZT13 between the MStr and LStr only (Str-between pair, F (3,44) = 7.52, p < 0.001). Raclopride minimally affected 3-8 Hz coherence, where it induced a decrease only in the striato-cerebellar comparisons at ZT1 and ZT13 (F (3,220) = 11.291, p < 0.001, Figures 6F-I). In general, changes in coherence in 0-3 Hz oscillations followed a diurnal pattern that is similar to the overall FFT power within these bands (compare Figures 4C  and 6). This implies that D2 receptor antagonism has the effect of enhancing both the power and coherence of slow oscillations, in a manner that is more potent at ZT13.

DISCUSSION
This study was aimed at determining if the neural network activity of the basal ganglia and the cerebellum, and their coupling, are influenced by time-of-day, and how raclopride, a D2 dopamine receptor antagonist, affects those relationships. A marked diurnal change in the rhythmicity of the striatal and cerebellar cortex networks was observed, with slow oscillatory activity decreasing at the beginning of the dark phase (when the rat would usually be more active), and 3-8 Hz activity increasing in the same period. The diurnal patterns of oscillations in the slow vs. 3-8 Hz bands were generally inversed, particularly in the cerebellum where 3-8 Hz oscillations were more prominent. This diurnal variation was also seen in the incidence of spindle activity, which was strongly influenced by the presence of the slow wave component. In addition, the shifting of activity from slow to 3-8 Hz at ZT13 was decreased by the administration of raclopride, which promoted slow band activity. Lastly, LFP coherence showed that cerebellar, striatal between-pair and striato-cerebellar slow wave coherence showed a strong diurnal pattern of modulation. Raclopride administration also influenced striatal coherence in a time-of-day dependent manner that was consistent with the shifts in oscillatory activity occurring in this area, such that coherence was enhanced when slow wave activity increased. These results show a clear diurnal modulation of striatal and cerebellar network activities and their coupling, and a time-specific effect of dopamine transmission that may serve to sculpt the activity of these networks.

OSCILLATORY ACTIVITY AND COUPLING UNDER URETHANE ANESTHESIA
Our experimental preparation, using urethane anesthesia, was chosen in order to minimize the chance of LFP signal contamination by behavioral artifacts such as neural movement-related LFP signal changes. Using the conventional parameters of heart rate and foot-pinch reflex, depth of anesthesia was stable across conditions. Urethane anesthesia is only minimally depressive, and thus permissive to neural oscillations (Maggi and Meli, 1986;Steriade, 2003;Clement et al., 2008). As such, urethane could be used for studying circadian effects on oscillations and communication between networks. This preparation favors slow wave oscillations and widespread synchrony (circa 1 Hz), while permitting transient oscillations in the 3-8 Hz range. In the context of sleep and urethane anesthesia, these slow oscillations have sculpting features in the neocortex (Steriade, 2003;Clement et al., 2008) whereby cortical "up" and "down" states greatly influence larger scale networks (Wilson and Kawaguchi, 1996;Stern et al., 1998;Ros et al., 2009). This cortical influence helps local circuits maintain synchronized activity during sleep, or sleep-like states (Sanchez-Vives and McCormick, 2000). We observed that cortical LFPs had a close relationship with striatal LFPs, and were more dissimilar with the cerebellar LFPs (Figure 6). Still, Ctx-CB coherence between 0-3 Hz could reach above 0.7, potentially indicating a strong cortical influence. In both structures, cortical slow oscillations can influence the cell membrane potential and firing pattern: slow oscillations are indeed related to granule and Golgi firing, and Purkinje cell complex spikes (Ros et al., 2009). Slow oscillations also have a strong influence on neuronal excitability in the striatum (Wilson and Kawaguchi, 1996;Sharott et al., 2012). The slow oscillations that we observed, synchronized in the striatum and cerebellar cortex, are therefore likely to arise from a common contribution of the neocortex, which may drive activity in both regions. Spindles recorded under urethane anesthesia are under the influence of thalamic pacemakers, particularly the reticular thalamic nucleus (Steriade, 2003). They occur during cortical "up" states and are dependent on the corticothalamic inputs. From our results, it appears that the diurnal modulation was more driven by the presence of the slow oscillations, further confirming the potential extrinsic influence on the local striatal and cerebellar cortex networks. Oscillatory interactions between the striatum and cerebellum could, however, involve a variety of reverberating circuits. The faster 3-8 Hz oscillations appear less uniform in their expression across recording sites and may be generated by more local processes that could be differentially organized in the cerebellum and striatum (Schnitzler and Gross, 2005). In the cerebellar cortex, the granule and Golgi cells seem to follow 5-30 Hz oscillations in a steady manner (Courtemanche et al., 2002;Dugué et al., 2009). However, faster 10-25 Hz oscillations occur in the awake primate and are stopped by movement (Pellerin and Lamarre, 1997;Courtemanche et al., 2002); they can also be synchronized with the neocortex under specific behavioral conditions (Courtemanche and Lamarre, 2005). These oscillations follow a parasagittal modulatory pattern (Courtemanche et al., 2009 adapting to the demands of the movement being performed. As for the striatum, both 4-12 Hz oscillations in rodents, and 10-30 Hz oscillations in primates also adapt to the task demands (Courtemanche et al., 2003;DeCoteau et al., 2007a;Thorn et al., 2010). This provides evidence of a capacity of the faster oscillations, and their synchrony, to adapt to behavioral conditions, organizing the networks in a task-dependent manner with a potentially more local organization than the slow oscillations.

DIURNAL AND DOPAMINE MODULATION OF LFP OSCILLATIONS
Many biological processes and daily changes in behavior are influenced by circadian gene expression patterns, yet little is known about the neural mechanisms that may be driving these changes. The striatum and the cerebellum are important subcortical structures for a number of functions subserving behavior, including motivation, adaptability of motor control, sequencing, and elaboration of cognitive and mental programs (Ito, 2006;Graybiel, 2008). To the best of our knowledge, this study is the first time subcortical LFP network activity has been analyzed in a circadian context. Our findings are consistent with the diurnal modulation of EEG oscillations observed in other brain sites (Grasing and Szeto, 1992), even though in this latter study, both sleep and awake states were included in their analysis. When states of sleep vs. wakefulness are compared, an increase in slow wave activity and inter-regional synchrony is expected early in the light phase (ZT1) when the rat would be falling asleep and, conversely, the slow oscillations should give way to faster oscillations in the transition into the dark phase (ZT13), when the animals become more active (Buzsáki and Draguhn, 2004). We found here that, even under urethane anesthesia, in both the striatum and the cerebellum, slow wave oscillations showed a clear diurnal variation that was consistent with this sleep-wake cycle. This indicates that the circadian mechanisms are robust enough to control network oscillations under anesthesia, suggesting that circadian clocks may influence how signal transmission occurs between local networks and play a role in establishing neural conditions that allow behavioral states to change throughout the day. The acceleration of the LFP oscillations and the shifts in coherence at the beginning of the dark period, at a time when the animal would normally be transitioning into its active phase, could be establishing optimal neural conditions that create a permissive state for networks to perform their awake activities. Changes in neurotransmitter availability influence neuronal excitability and can contribute to the switching between oscillatory modes (McCormick, 2004), where different neuromodulators can promote shifts in frequencies in a manner that is anatomically specific (Roopun et al., 2010). Dopamine is known to promote wakefulness, but has scarcely been studied in its effect on cortical excitability during slow wave activity and our results seem similar to the neuromodulatory effects of dopamine in the  (Mayne et al., 2013;Zeitzer, 2013). We found here that blocking D2 receptors could increase the slow oscillatory content of the LFP signal. Raclopride injections enhanced slow striatal and cerebellar oscillations, particularly at ZT13. This is the time of day that has the greatest reduction on the slow oscillations. This result points to a role of dopamine as a neuromodulator of the slow oscillations. Dopamine itself has a strong influence on striatal oscillatory properties (Costa et al., 2006). Under urethane anesthesia, dopamine depletion increases basal ganglia entrainment to cortical slow oscillations (Tseng et al., 2001;Walters et al., 2007). Lesion to substantia nigra dopamine neurons causes depolarization of striatal medium spiny neurons and increases their spontaneous burst firing at frequencies coherent with the neocortex (Tseng et al., 2001;Sharott et al., 2012). This results in a decrease in the striatal filtering of the slow components from the neocortex, consequently increasing slow wave activity throughout the basal ganglia. Therefore, in our experiments, the decrease in slow wave activity in the striatum at ZT13 could be due to an increase in dopamine transmission. This is supported by raclopride having the largest effects at this time. In contrast, peak extracellular dopamine concentrations occur in the striatum somewhat later in the dark phase, around ZT18 (Owasoyo et al., 1979;Hood et al., 2010;Ferris et al., 2014). This suggests that oscillatory activity is being affected by factors in addition to dopamine availability, such as membrane receptor activation. Diurnal variations in the striatal LFP are likely dependent, in part, on alteration of dopamine function, which may help to provide an oscillatory state in the striatum which is more "disconnected" from the neocortex during awake conditions, thus enhancing the occurrence of higher frequency oscillations.

NETWORK INTERACTIONS
In the cerebellar LFPs, we found slow oscillations that are thought to have a neocortical origin (Ros et al., 2009), and 3-8 Hz oscillations, which are likely dependent on reverberating properties in the cerebellar granule cell layer (Dugué et al., 2009;Courtemanche et al., 2013). The cerebellum followed a similar diurnal pattern in oscillations as the striatum (mostly slow oscillations at ZT1; fewer slow oscillations and more 3-8 Hz at ZT1) and its oscillatory profile was similarly affected by raclopride. While the oscillations recorded could in part come from local cerebellar circuit resonance (Dugué et al., 2009), it is quite possible that there are contributions from larger circuit interactions, including the subcortical connections between the basal ganglia and the cerebellum and the neocortico-ponto-cerebellar circuits.
Addressing the local circuits, it is unclear exactly how dopamine can affect cerebellar cortex circuitry. Even if not traditionally considered as a prominent cerebellar neuromodulator, there is evidence for dopaminergic transmission in the cerebellum (Takada et al., 1993;Hurley et al., 2003;Delis et al., 2004;Schweighofer et al., 2004;Giompres and Delis, 2005) There is also some evidence that dopamine levels in the cerebellum follow a circadian cycle (Owasoyo et al., 1979). However, the role of dopamine in the cerebellum remains poorly understood (Schweighofer et al., 2004). D2 and D3 type dopamine receptors have primarily been reported in the molecular cell layer of the vermis (Bouthenet et al., 1987), however raclopride does show low levels of binding to the cerebellar hemispheres in vivo (Kiss et al., 2011). Dopamine projections from the ventral tegmental area (VTA) to the cerebellar cortex have been reported to terminate in the granule and Purkinje cell layers in the posterior lobe (Ikai et al., 1992) and alterations of activity in the VTA both acutely and chronically induce changes in cFOS expression in the granule cell layer in multiple cerebellar regions, including Crus1 and 2 (Herrera-Meza et al., 2014). It may be that these VTA projections were affected by the systemic raclopride injections, altering neural activity in the areas we recorded, however the effects of dopamine on neural activity in the cerebellar cortex is also not well described. It is unknown if dopamine affects the resonance aspect in the cerebellar circuits, enhancing the 3-8 Hz resonance that appears optimal in those circuits, reducing the capacity to entrain local circuits to the 0-3 Hz inputs arriving from the cortex via the pontine nuclei. We have seen interactions between the GCL and the Purkinje cell layers in monkey ∼15 Hz resonance, so faster oscillations could be partially affected by this local connectivity . A potential mechanism could be the Lugaro connection affecting the Golgi cells, but at this point, since sufficient studies have not addressed this, these effects remains speculative.
A more probable action of dopamine in the cerebellar cortex however is its indirect effects through cerebellar interactions with other brain areas. The cerebellum is likely to interact with the striatum through the neocortex and subcortical connections. Interconnections between the neocortex and striatum as well as those between the neocortex and cerebellum are strong and numerous. Cortico-striatal connections are important inputs to the striatum (Graybiel, 2010), and these interactions are affected in disease models known to afflict the basal ganglia (Crittenden and Graybiel, 2011). In addition, corticopontocerebellar connections have an important effect on cerebellar processing (Morissette and Bower, 1996). This pathway constitutes one of the fastest central pathway in the CNS (Allen and Tsukahara, 1974). Overall, the potential avenue of communication between the striatum and cerebellar cortex taking advantage of the neocortical loops is a definite possibility. In addition, this interconnection is supplemented by the more recent finding of a plurisynaptic subcortical route between the cerebellum and basal ganglia (Hoshi et al., 2005;. In these papers and in their review,  identify that these pathways permit bidirectional information transfer from the basal ganglia to the cerebellum (via a connection from the subthalamic nucleus to the pontine nuclei), or vice-versa (via a connection from the dentate nucleus to the thalamusspecifically, the thalamic central lateral nucleus (Ichinohe et al., 2000)), without requiring the passage through cerebral cortex connections. This anatomical connectivity is supported both by molecular and behavioral studies reporting alterations in either the basal ganglia or the cerebellum inducing rapid changes in the other area (Koch et al., 2009;Calderon et al., 2011;Moers-Hornikx et al., 2011). Since functional changes appear to be bidirectional, one would also assume changes in neural activity to have a similar bidirectional influence. Electrical stimulation of the dentate nucleus can increase firing in the basal ganglia  (Li and Parker, 1969;Ratcheson and Li, 1969), suggesting that altered firing in the cerebellar cortex could also be produced by modifying the neural activity in the striatum. Dopamine may also affect circadian changes in oscillatory activity by altering more widespread circuit properties, such as thalamocortical circuits, which could then cascade down an effect on both the cerebellar and striatal oscillations. This is consistent with findings that systemic D2 antagonism reduces wakefulness and increases slow wave sleep, where it is believed to act through the nuclei responsible for regulating sleep and wakefulness in the brainstem, hypothalamus, and basal forebrain (Monti and Monti, 2007). With our particular LFP patterns showing widespread synchrony ∼1 Hz, and this oscillation having a strong link with the neocortex in the cerebellum (Ros et al., 2009) and striatum (Stern et al., 1998), it seems most likely that the coordination of the signals we recorded implicates neocortical activity. This would also support the appearance of spindles in the cerebellar LFPs. Spindles originate in the reticular thalamic and neocortical loops under conditions of slow wave oscillations, once appearing in the neocortex they could be transmit to the cerebellar cortex via the cortico-pontocerebellar pathway. It is possible that the neocortex circuits are also involved in the transmission of the 4-12 Hz activity between sites; many cortical areas do exhibit this 4-12 Hz activity, including the hippocampus (Buzsáki, 2006;Buzsáki and Moser, 2013), the cerebral somatosensory cortex (Nicolelis et al., 1995;Ahissar et al., 1997), and also the prefrontal cortex (Watson et al., 2014). However, this activity is better controlled locally, and less driven by large cortico-thalamic interactions that synchronize highly converging neocortical structures, as seems to be the case for slow rhythms during urethane anesthesia or sleep (Steriade, 2003).
The diurnal alterations between slow and 3-8 Hz oscillations are therefore likely influenced by the balance between widespread synchrony from the thalamocortical networks driving the slow oscillations and mechanisms underlying more locally generated oscillations. This is partly evidenced by the fact that when placing the electrodes, the cerebellum consistently showed weaker slow oscillations and stronger 3-8 Hz oscillations as compared to the striatum, reflective of a potentially more distant connection to the neocortex and stronger local resonance of this 3-8 Hz frequency range. This regional aspect would provide some rationale for the favoring of subcortical connections in the synchronization of theta-like activity patterns. A limitation of this study however, is that cortical activity can only be inferred in its role for driving these oscillations in both areas and to decipher regional connectivity will necessitate a more direct approach on the connectivity.

FUNCTIONAL AND CLINICAL IMPLICATIONS
This study uncovered a diurnal modulation that emphasizes the capacity of the striatal and cerebellar circuits to reorganize during the course of the day. The electrophysiological evidence presented here, with oscillations and synchrony alternating between ∼1 Hz or ∼4 Hz rhythms, shows that neural properties are being modulated in order to switch the underlying network states. This could allow the networks to optimally engage in various behaviors throughout the day. The drive in the diurnal modulation of these signals likely comes from the SCN, the central circadian pacemaker, and the subsidiary clocks that exist in a hierarchical network throughout the brain, including the striatum and the cerebellum (Namihira et al., 1999;Shieh, 2003;Imbesi et al., 2009;Rath et al., 2012;Harbour et al., 2013). Both latter areas are involved in locomotor activities driven by circadian gene expression (Masubuchi et al., 2000;Hood et al., 2010;Mendoza et al., 2010) and the diurnal modulation of LFP activity described here is likely to reflect network processes that promote these locomotor/movement activities and may further be important for fine-tuning sensorimotor performance and learning throughout the day. Future directions will include making the link between daily changes in neural network activity and behavioral processes that are driven by circadian genes.
There are a few neurological clinical conditions that affect both the basal ganglia and the cerebellum, and which display diurnal variations in symptoms or in response to pharmacological treatments. For example, Parkinson's disease patients with motor symptom fluctuations generally experience the least amount of impairment in the morning, with symptoms increasing throughout the day (Nutt et al., 1997;Bruguerolle and Simon, 2002). In addition, Parkinson's patients receiving dopamine treatments such as levodopa also experience a diurnal worsening of symptoms. Similarly, spinocerebellar type 3 ataxic symptoms also appear to show diurnal variations in symptoms (Wilder-Smith et al., 2003), though this seems to interact with a dyskinesia that is also levodopa responsive. The link could be direct with core function in these networks, as the ATX2 (ATAXIN-2) protein, which is implicated in the expression of spinocerebellar ataxia type 2 (or an increased risk of amyotrophic lateral sclerosis and Parkinsonism), is involved in the activation of the rate-limiting circadian clock component PERIOD in Drosophila (Lim and Allada, 2013) and their daily locomotor behavior (Zhang et al., 2013). In terms of circuit mechanisms, Nutt et al. (1997) report that circadian fluctuations in motor symptoms are independent of plasma levels of levodopa. Even when administering constant-rate infusions, symptoms still worsen despite plasma levels increasing throughout the day. This suggests that dopamine availability is not the only factor that is determining variations in motor behavior, but other factors such as receptor activation and interactions with other neuromodulators also likely play a role. Our results suggest that membrane targets responsible for changes in circuit properties change throughout the day, rather than simply availability of dopamine at the synapse. These daily alterations in circuit properties, both dopamine-related and otherwise, could underlie the diurnal symptomatology in basal ganglia and cerebellar patients.