Abstract
Modulation of oscillatory activity through basal ganglia–cortical loops in specific frequency bands is thought to reflect specific functional states of neural networks. A specific negative correlation between beta and gamma sub-bands has been demonstrated in human basal ganglia and may be key for normal basal ganglia function. However, these studies were limited to Parkinson's disease patients. To confirm that this interaction is a feature of normal basal ganglia, we recorded local field potential (LFP) from electrodes in globus pallidus (GP) of intact rats. We found significant negative correlation between specific frequencies within gamma (≈60 Hz) and beta (≈14 Hz) bands. Furthermore, we show that fluctuations in power at these frequencies are differentially nested within slow (≈3 Hz) oscillations in the delta band, showing maximum power at distinct and different phases of delta. These results suggest a hierarchical organization of LFP frequencies in the rat GP, in which a low-frequency signal in the basal ganglia can predict the timing and interaction of power fluctuations across higher frequencies. Finally, we found that dopamine D1 and D2 receptor antagonists differentially affected power in gamma and beta bands and also had different effects on correlation between them and the nesting within delta, indicating an important role for endogenous dopamine acting on direct and indirect pathway neurons in the maintenance of the hierarchical organization of frequency bands. Disruption of this hierarchical organization and subsequent disordered beta–gamma balance in basal ganglia disorders such as Parkinson's disease may be important in the pathogenesis of their symptoms.
Introduction
Oscillatory synchronization of neuronal assemblies is thought to enable periodic functional linkages between brain areas for specific information processing (Buzsáki et al., 1994; Singer, 2001; Fries et al., 2002, 2007). In neural circuits of the basal ganglia, modulation of oscillatory activity at different frequencies is particularly important for normal processing (Brown et al., 2002; Courtemanche et al., 2003; Berke et al., 2004).
Data from recordings in human subthalamic nucleus (STN) have shown that the power in gamma and beta sub-bands is negatively correlated, in the absence of movement (Fogelson et al., 2005). This suggests the possibility that a regulatory mechanism controls the power relationship between these bands within a single behavioral state, over and above any alternation in power that may occur across different movement states. Furthermore, in other brain regions not directly within the basal ganglia, evidence is accumulating for phase–amplitude coupling between low- and high-frequency oscillations. This has been interpreted as reflecting a fundamentally hierarchical organization of oscillatory frequencies (Bragin et al., 1995; Lakatos et al., 2005; Sirota and Buzsáki, 2005; Canolty et al., 2006; Isomura et al., 2006; Mena-Segovia et al., 2008; Sirota et al., 2008). However, the human STN data were obtained from patients with Parkinson's disease, leaving the possibility that the results reflected this pathology. The first aim of the present study was therefore to determine whether such hierarchical organization of frequencies occurs in the normal basal ganglia.
If such a hierarchical organization exists, an additional question of interest is the role of dopamine D1 and D2 receptors in regulating it. The striatum, STN, and the globus pallidus (GP) all receive dopaminergic input (Smith et al., 1998; Gauthier et al., 1999). Dopamine D1 and D2 receptors are very highly expressed in the basal ganglia but differentially in the direct and indirect pathways (Surmeier et al., 2007). The D1 receptor is found on “direct” pathway striatal neurons that project to the internal globus pallidus. In contrast, the D2 receptor dominates in the “indirect” pathway, which passes through the external GP and STN (Smith et al., 1998). Interaction between the tightly interconnected GP and STN (Bevan et al., 2002) has been proposed as an engine for basal ganglia oscillations (Plenz and Kital, 1999; Humphries et al., 2006), raising the possibility of a special role for D2 receptor in the regulation of basal ganglia rhythmic activity by dopamine. We therefore further investigated whether D1 and D2 receptors have differential roles in modulating basal ganglia oscillatory activity, using selective receptor antagonists.
Materials and Methods
All procedures were approved by the Animal Ethics Committee at the University of Otago (approval 59/04) and conducted in accordance with the New Zealand Animal Welfare Act. Eleven male Wistar rats weighing 250–400 g at the time of surgery were used for the main study. Food and water were available ad libitum in the home cage, and experiments were performed during the dark portion of a reversed 12 h light/dark cycle.
Recording electrodes were implanted after at least 7 d of daily handling habituation, under full general anesthesia [70 mg/kg ketamine and 0.5 mg/kg Domitor (medetomidine HCl), s.c.] with local anesthetic (0.5% Marcaine) at the incision site and under prophylactic antibiotic [Procaine Penicillin G and dihydrostreptomycin sulfate (Strepcin), 1 ml/kg, s.c.]. The recording electrodes [eight 0.001-inch coated diameter Formvar-insulated nichrome microwires (A-M Systems), fixed in a 30-gauge stainless-steel cannula] were lowered to position through a small craniotomy. The electrode tips were targeted to the dorsal margin of the left GP (1.1 mm posterior to bregma, 3.1 mm lateral to midline, and 5.0 mm below the dura surface) and advanced by small steps (≈50 μm) each day using an on-head microdrive. Seven stainless-steel skull screws were inserted, and the implant was affixed to the skull with dental acrylic. A contralateral skull screw positioned over the cerebellum ∼14 mm caudal of bregma was used for animal ground connection. Animals were injected with the sedation-reversal agent Antisedan (atipamezole HCl, 1.25 mg/kg, s.c.) and allowed to recover before being returned to the animal housing facilities.
Beginning 5 d after surgery, the electrodes were advanced daily in small increments for 2–3 weeks, until the tips reached the depth coordinate for the center of the GP. Electrode positions were marked at the end of the experiments by electrolytic lesions (10–50 mA direct current for 20–30 s), and only those animals with tips confirmed within the GP were included in analysis.
Recordings were made within a clear acrylic box in a dimly lit room. After connection of the recording cable (incorporating impedance-matching unity-gain amplifiers), animals were allowed to settle for 15–20 min. Recordings of local field potentials (LFPs) were then made for 10 min before (predrug) systemic injection of raclopride, SCH23390 [R(+)-7-chloro-8-hydroxy-3-methyl-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine hydrochloride], or vehicle and again for 10 min, beginning 5 min after the injection (postdrug). LFP recordings were made from the single electrode with the best signal-to-noise ratio in single-ended mode (referenced to amplifier signal ground), bandpass filtered (0.1–400 Hz; Axon CyberAmp 380 signal conditioner; Molecular Devices), digitized at 1 kHz (Power1401 data-acquisition system; Cambridge Electronic Design), and recorded using Spike2 software (Cambridge Electronic Design). The bandpass filters for other wires in the bundle were set at 0.5–10 kHz, in case single neurons were present.
In a separate group of three animals, 0.001-foot nichrome wire tetrodes were used instead of bundled single wires to improve yield of single unit spikes. Raw data were initially recorded as wide band (0.1 Hz to 20 kHz) signals and split into two copies; one was downsampled to 1 kHz for LFP analysis and the other bandpass filtered (100 Hz to 20 kHz) for spike waveform detection (supplemental data, available at www.jneurosci.org as supplemental material). This allowed spiking and LFP to be analyzed from the same electrodes. Data from this group were used to check for phase relationship of cell spiking activity with locally recorded LFP.
Drugs were dissolved in 0.9% saline and administered intraperitoneally. Raclopride (Sigma-Aldrich), a selective dopamine D2 antagonist, was dissolved at 1.0 mg/ml and the selective D1 antagonist SCH23390 (Sigma-Aldrich) at 0.5 mg/ml with both drugs administered in a volume of 1.0 ml/kg for final doses of 1.0 and 0.5 mg/kg, respectively. It was important for the present study that the behavioral effect of the two drugs was matched to allow comparisons and that the behavioral state was stable for the period of the recordings. Doses of the two drugs were chosen on the basis of preliminary experiments as sufficient to produce a comparable cataleptic response [similar levels of akinesia measured by the bar test (Marin et al., 1993)] and are higher than required to affect operant behavior. At this dose, SCH23390 remains selective for D1 over D2 receptors (Bischoff et al., 1986). Saline injected at 1.0 ml/kg served as a negative control. Drug tests were made at least 2 weeks apart.
The animals were continuously monitored during data acquisition, and activity periods were marked on an event channel of the recorded data file using keystroke input. Example traces are shown in supplemental Figure S1 (available at www.jneurosci.org as supplemental material). During postdrug catalepsy, animals were primarily akinetic. Movement can be a powerful modulator of power in gamma and beta bands, as is illustrated in supplemental Figure S2 (available at www.jneurosci.org as supplemental material). Therefore, to ensure that predrug–postdrug comparisons were not confounded by changes in the level of movement, analysis was limited to periods of immobility (no locomotion, rearing, or grooming) of at least 3 min duration, from within the 10 min predrug and postdrug recording periods. Thus, all postdrug data were obtained from within a time window extending from 5 to 15 min after injection. Periods with high voltage spindles were excluded in the recordings. Animals remained upright, and, although rapid eye movement sleep episodes cannot be ruled out, large-amplitude slow waves typical of slow-wave sleep were not observed. For each animal, the selected predrug and postdrug periods were of equal lengths.
The level of akinesia achieved by drug administration was quantified at the end of the postdrug recording period (i.e., 15 min after injection) using the bar test (Marin et al., 1993). For this, the animal's front paws were placed on a horizontal bar 10 cm above the bench, and the time until a single paw moved from the bar was recorded, to a maximum of 120 s.
Frequency spectrum analysis.
A second-order Butterworth band-stop filter with settings 49 and 51 Hz (Matlab; MathWorks) was used to remove 50 Hz mains interference. The power spectral density (PSD) of the LFP was then computed using fast Fourier transform (FFT), with the Welch spectral estimator and a Hamming window of 2048 points (2 s) with 50% overlap and 1 Hz resolution. A significant peak in the power spectrum was defined as four consecutive increasing power values.
Power correlation analysis.
We analyzed the temporal correlation between spectral power at different frequencies f as described previously (Masimore et al., 2004), by computing the time-resolved PSD of the LFP. The appropriate time resolution for evolving PSD to capture power fluctuations of beta and gamma was determined from inspection of the raw signal, and was set at 50 ms. We used FFT with a 75% overlapping Hamming window of 200 points (0.2 s) in the 0.1–1000 Hz range with 1 Hz resolution, with each window representing a single time point in the time-resolved PSD. This window size does not allow analysis of oscillations with period <200 ms, so the present analysis was limited to frequencies between 5 and 100 Hz. This generated time series of power values [Pf1,I, Pf1,i + 1, Pf1,i + n …; Pf2,I, Pf2,i + n; etc.] corresponding to each 50 ms epoch I, for each frequency f. Because power values are not Gaussian, we calculated the nonparametric Spearman's rank correlation (ρ(f1, f2)) for all possible pairs of power time series, generating two-dimensional matrices of negative and positive spectral correlations for each animal. The significance of the negative and positive correlation values were calculated using the nonparametric Fisher's test.
Gabor analysis.
We used Gabor functions to further investigate temporal fluctuations in power in LFP bands. Gabor functions are commonly used to fit autocorrelation (AC) and cross-correlation histograms of nonstationary rhythmic biological time series such as neuronal spiking activity (Gray et al., 1989; Engel et al., 1990; Young et al., 1992). For our LFP analysis, they are damped sine waves, with two terms, first the sine wave frequency (fo), the amplitude of which is modulated by the second term, here the damping frequency (fd). The Gabor functions served as predicted AC (pAC), used to fit the actual AC of the LFP frequencies of interest.
We constructed a set of Gabor functions as follows (Young et al., 1992; Nargeot et al., 2007): with fo and fd both ranging from 1 to 100 Hz, hence creating 100 × 100 predicted AC. The quality of the fit of each predicted AC was then assessed by its correlation (Spearman's ρ) with the actual AC of specific LFP bands (calculated for lags t of 0–500 ms), and this correlation score was plotted for each f0, fd pair. Points showing the highest correlation thus represent candidate f0, fd pairs capable of predicting fluctuations in LFP power.
Perievent time power spectral analysis.
To determine whether fluctuations in LFP power were linked to a specific phase of candidate fd frequencies revealed by the Gabor analysis, we generated averaged perievent time power spectra (PETPS). In these PETPS, the power of each frequency was calculated over lags of ±0.5 s with respect to the troughs in the candidate fd from 100 s of data for each animal. First, PSD (1 Hz resolution) were computed for 75% overlapping Hamming windows of 200 points (0.2 s), for 1 s periods centered on each trough. These individual PSD were then averaged to generate the PETPS for each animal. Finally, to compute unbiased group average data, the power values for each animal were normalized by converting to z-scores (zPETPSf), as follows: where f denotes the frequency component being analyzed, and t is time relative to the trigger. These zPETPS values were then averaged across animals.
To quantify the relationships suggested by the PETPS analysis, we computed the phase of maximum power with respect to the candidate fd using the Hilbert transform in Matlab, with LFP troughs set to 0° phase. We then determined the local power maxima within 0° and 360° of the trough. The mean phase value for the maximum was calculated for each rat, and the distribution across all rats was analyzed. The presence of a significant preferred phase was tested using the Rayleigh's test for directional data (Fisher, 1993). For this, the 0–360° circular space was divided into 12 bins of equal size, giving an angular resolution of 30°. The Rayleigh's test is based on the hypothesis that phases of power maxima are uniformly distributed along the circular space. The presence of a significant preferred phase of a given frequency was calculated for data before injection of saline (control) and after saline, SCH23390, and raclopride injections. We assessed the effect of each drug on the preferred phase using the Watson's U2 test (Fisher, 1993).
Results
We recorded LFP in the GP of nine Wistar rats with confirmed electrode positions (Fig. 1A). Spectral analysis of LFP revealed two distinct peaks (Fig. 1B). One was located between 13 and 17 Hz across animals, with a mean ± SD of 14.8 ± 1.3 (low beta). We refer to this frequency component as β14. The second peak was found within the gamma band at 54–61 Hz across animals, with a mean of 59.3 ± 2.4 Hz, and is referred to here as γ60. The frequency spectrum was also characterized by the presence of considerable power density at slower oscillatory frequencies, peaking in different animals at frequencies between 0.1 and 10 Hz but with no particular preferred frequency in this range.
These recordings were made with fixed bandpass optimized for LFP, so it was not possible to record action potential spikes on the same wires. In some cases (n = 6), cells were present on other wires in the bundle. For all these neurons, spike-triggered averaging of LFP suggested that they fired at a preferential phase of beta and/or gamma band oscillations (supplemental Fig. S3, available at www.jneurosci.org as supplemental material). However, although the electrodes were all in a single bundle, the exact spatial relationship between the LFP and neuronal recording locations was not known. Therefore, to further assess the extent to which neurons in GP fire in relation to locally recorded LFP signals, we implanted tetrodes into GP of another group of three rats and analyzed spike and LFP data acquired from wide-band recordings obtained from the same tetrodes. Example phase histograms are shown in supplemental Figure S4 (available at www.jneurosci.org as supplemental material). From 19 neurons, 14 (74%) showed significant entrainment to LFP (p < 0.05, Rayleigh's test for homogeneity), with five linked to both gamma and beta, five to gamma only, and four to beta frequency oscillations.
Interplay between frequency bands within GP LFP
Figure 1C shows a short section of typical recording from one rat filtered at 0.1–8, 12–30, and 55–65 Hz bands, during which the amplitude of oscillations (the power) in the separate bands can be seen to fluctuate in specific relation to each other. In particular, when β14 is high amplitude, γ60 tends to be low and vice versa, i.e., power in the two bands appears to be negatively correlated. Quantitative analysis confirmed that this was a significant feature in all rats. For this, we applied a power correlation analysis (Masimore et al., 2004; Fogelson et al., 2005), which returns the correlation coefficient between every possible pair of frequencies. An example analysis for one rat is shown in Figure 1D, which reveals a clear, significant negative correlation between narrow frequency bands centered at ≈14 and ≈60 Hz. All rats (n = 9) revealed similar significant interactions. No other frequency pairs were significantly positively or negatively correlated in >50% of animals, apart from the expected high positive correlation on and close to the diagonal (i.e., correlation of a frequency with itself and the neighboring frequencies).
Switching between beta and gamma dominance is predicted by delta phase
If a hierarchical organization of LFP frequency bands exists in basal ganglia, then the power fluctuations observed in γ60 and β14 may be related to oscillations at lower frequencies. To explore this, we first tested whether fluctuations in γ60 and β14 power, as expressed in the shape of their AC, can be predicted by Gabor functions derived from different damping frequency terms (see Materials and Methods). Our hypothesis was that, if the predicted AC derived from a damped sinusoid can fit the actual AC, then the damping frequency or frequencies that provided the best fit with the actual data would be candidates for playing a part in such a hierarchy. The 100 × 100 matrices in Figure 2A show the Spearman's correlation of the actual AC for γ60 and β14 bands, with all the Gabor-predicted AC derived for each combination of fo and fd terms, averaged across animals (n = 9 rats). The data clearly show maximum correlation values for fo terms close to 14 and 60 Hz, as would be expected. More interestingly, for both γ60 and β14, the damping term (fd) with the highest correlation occurred in the delta range, with maxima centered on fd values of ≈2–3 Hz. We defined this frequency as δ3. An illustrative example showing time course of the damping of the AC for γ60 and β14 compared with the duration of δ3 AC, the overlap between the actual and predicted AC, and the associated scatter plots for a single animal is shown in Figure 2B, using fd = 3 Hz. Both predicted and actual AC for both frequencies can be seen to decrement within the temporal envelope defined by the δ3 AC function.
The Gabor analysis indicates a potential relationship between δ3 and each of the two frequencies that were negatively correlated with each other. Another aspect of interest is whether the negative correlation itself could also be related to δ3. For instance, in the example data trace in Figure 1C, periods of high-amplitude β14 tended to be at high amplitude within delta band peaks, whereas γ60 amplitude appeared to increase during delta troughs. To determine whether such nesting of power fluctuations in each band within different and specific phases of delta rhythm was indeed a significant characteristic of the data, we computed PETPS (Fig. 3). This analysis reveals the timing of fluctuations in β14 and γ60 with respect to the time of delta troughs (see Materials and Methods).
An example PETPS analysis from one animal is shown in Figure 3A. In this example, it can once again be seen that power in β14 and γ60 fluctuates in anti-phase and that β14 is at maximum power during delta peaks, whereas γ60 power maxima are synchronous with the triggering δ3 trough. These timing relationships were consistent across animals, as shown by the normalized group data from the PETPS analysis shown in Figure 3B and confirmed by the statistical analysis illustrated in Figure 3C. For this, we computed the delta phase at which maximum power was found for β14 and γ60. Across all animals, the distribution of maximum power for both these frequencies was significantly clustered at specific phases (p < 0.05, Rayleigh's test), and the phase position of the clusters for γ60 and β14 were significantly different from each other (Watson's U2 test). As a negative control, we also analyzed a third frequency, 40 Hz, which showed no peak in the power spectral density (Fig. 1B) and was not involved in any significant interfrequency interactions (Fig. 1D). In contrast to γ60 and β14, the trough-triggered average for this frequency was relatively flat (Fig. 3B), and the distribution of 40 Hz power maxima was not significantly different from random (Fig. 3C).
Oscillation power and interfrequency interplay is under control of dopamine through differential actions at D1 and D2 receptors
We investigated the role of endogenous dopamine in regulating the interfrequency dependencies described in the previous sections, by administering the selective dopamine D1 and D2 receptor antagonists SCH23390 and raclopride. Six animals were tested with raclopride, SCH23390, and vehicle and an additional three with raclopride only. Doses were selected such that the drugs generated similar levels of mild catalepsy, confirmed by elevated bar test times (mean ± SEM bar time: raclopride, 18 ± 6 s; SCH23390, 25 ± 9 s) compared with saline control injections, in which bar times were <1 s and so scored as 0 s in all cases (p < 0.05, paired t test, n = 6).
We first determined the effect of dopamine antagonists on LFP PSD. Example traces are shown in supplemental Figure S1 (available at www.jneurosci.org as supplemental material). Analysis of group data, shown in Figure 4A, revealed that drug treatments altered the amplitude of PSD peaks, without altering their characteristic frequency. The specific effects depended on the drug and the frequency band. Both drugs elevated the spectral density over the beta range. In contrast, only raclopride influenced the low gamma peak, and the effect was to flatten it compared with the predrug control. These observations were confirmed in the quantitative analysis shown in Figure 4B. For this, beta was divided into low and high beta sub-bands, and quantification of gamma was restricted to the low gamma sub-band that exhibited the peak in the control PSD. Other classically defined delta, theta, and alpha bands were included for comparison. There was a significant effect of treatment (Kruskal–Wallis test, p < 0.05), and post hoc tests confirmed that this was attributable to significant elevation in low beta power by both drugs compared with saline and significant depression of low gamma power by the D2 receptor antagonist raclopride but not SCH23390 (p < 0.05, Wilcoxon's signed rank test).
The box plots in Figure 4B also suggested the possibility of a trend toward increases in alpha power after all treatments, including saline. However, one-sample (Wilcoxon's signed rank) tests showed that the only point that was significantly different from no change was alpha band power after SCH23390 injection (p < 0.05), and this value was not significantly different from the post-injection values for saline or raclopride. These results suggest the possibility of a small nonspecific effect on alpha power from some aspect of the drug administration procedure itself.
Finally, we asked whether constitutive dopamine has a role in modulating interplay between gamma and beta and nesting within delta frequency bands. First, we compared power correlation matrices across drug treatments. Figure 5A shows example post-injection data from one animal. The post-saline matrix shows the same pattern as seen in Figure 1D for the predrug condition, with clear negative correlation in power between γ60 and β14. In contrast, the negative correlation almost completely disappeared after administration of either dopamine antagonist. This was confirmed in the quantitative analysis of group data shown in Figure 5B. Both raclopride and SCH23390 significantly reduced the mean value of the Spearman's correlation coefficient between power at γ60 and β14 (p < 0.05).
Second, to determine whether SCH23390 and raclopride differentially affected the nesting of β14 and γ60 within δ3, we analyzed PETPS under the different pharmacological conditions, as shown in Figure 6. This analysis revealed that the D1 and D2 antagonists affected power modulation and correlation patterns in a complex manner, which was differential both across frequencies, and between drugs. For the γ60 band, inspection of the plots in Figure 6A reveals that administration of the D2 antagonist SCH23390 completely reversed the phase relationship, confirmed by a significant shift in the position of the delta phase of the maximum power (Fig. 6B). As illustrated in Figure 6B, before injection the phase of γ60 maximum power was significantly clustered at ≈0° (i.e., at the time of the δ3 trough), but after SCH23390, it was significantly clustered at ≈180°, i.e., the time of the δ3 peaks surrounding the troughs (p < 0.05, Rayleigh's test). This new cluster position was furthermore significantly different from predrug (p < 0.05, Watson's U2 test). In contrast, the D2 antagonist raclopride had no effect on the relationship of γ60 power with δ3. The data in Figure 6 also reveal that drug effects were quite different for the β14 band. Here, both drugs simply eliminated the relationship of β14 power fluctuations with γ60 phase, so that normalized data plots (Fig. 6A) became flat and the position of power maxima became randomly scattered with respect to δ3 phase (Fig. 6B). It is important to recall that this loss of temporal relationship was in the face of a significant increase in mean β14 power. Saline administration did not affect the temporal patterning of power for either band.
Discussion
The main findings of the present study are summarized in Figure 7. First, we found that power in beta and gamma LFP bands is negatively correlated in rat GP, in resting animals. Second, we show that the modulation of beta and gamma power is nested within slow oscillations in the delta frequency band, with peak power in β14 and γ60 at different phases of δ3. Third, endogenous dopamine is a critical determinant of not only the total power in each of the bands but also the negative correlation between beta and gamma power, and the relationship power in these bands to delta phase, with differential roles for D2 and D1 receptors.
Although volume conduction of signals from distant sources cannot be ruled out, particularly for structures such as basal ganglia nuclei that lack a layered structure (Berke, 2009), we found that the majority of neurons recorded on the GP electrodes showed phase locking of activity to LFP oscillations. This entrainment of spiking to the LFP indicates that the rhythms recorded in the LFP waves at least partially reflect local electrophysiological phenomena. Some previous data showing enhanced coherence at gamma frequencies between GP and other structures in anesthetized rats (but during periods of cortical activation considered physiologically similar to the awake state) also support a local component to LFP signals recorded in GP (Magill et al., 2004; Sharott et al., 2005). There is also considerable evidence for expression of anesthetic-induced slow-wave activity in rat GP (Magill et al., 2004; Walters et al., 2007; Zold et al., 2007). How the local rhythms are generated remains uncertain, but cocultures of GP and STN neurons form networks that show intrinsic slow oscillations (Plenz and Kital, 1999), suggesting that basal ganglia networks including GP have the necessary machinery to generate such rhythms, and modeling studies indicate that the properties of the intact GP–STN network would support oscillations in both the beta (Holgado et al., 2010) and gamma (Humphries et al., 2006) ranges.
Nesting of negative beta–gamma power correlation within delta rhythm
A main finding of the present study is that moment-to-moment alternation (negative correlation) of power in beta and gamma bands is a feature of activity in basal ganglia of intact rat. This is consistent with findings in l-3,4-dihydroxyphenylalanine-treated Parkinson's patients (Fogelson et al., 2005) and suggests that there may be a fundamental linkage of the neural generators underlying beta and gamma frequency oscillations.
Furthermore, we found the periods of negatively correlated power fluctuations in beta and gamma bands to exist over short time frames that were defined by (nested within) delta phase. These data from unanesthetized rats therefore support and extend previous evidence for nesting of basal ganglia gamma within slower oscillations, in striatum and pedunculopontine tegmental nucleus of anesthetized rats (Stern et al., 1997; Mena-Segovia et al., 2008) and striatum of unanesthetized rats after amphetamine administration (Berke, 2009) or during navigation tasks (Tort et al., 2008). Nesting of gamma within lower-frequency waves is well known in cortico-hippocampal networks and is thought to arise from the interaction of inputs oscillating at low frequency with local circuit mechanisms for generation of higher-frequency oscillations (Bragin et al., 1995). It is plausible to imagine that similar processes operate through the basal ganglia network.
It is well known that beta and gamma in the cortex–basal ganglia network are modulated according to movement state (Bergman et al., 1994; Donoghue et al., 1998; Brown and Marsden, 1999; Aoki et al., 2001; Brown et al., 2001, 2002; Cassidy et al., 2002; Levy et al., 2002; Lee, 2003; Masimore et al., 2005; Sharott et al., 2005; van der Meer and Redish, 2009). In contrast, the present results, and data showing negative power correlation between beta and gamma bands in humans (Fogelson et al., 2005), were obtained at rest. Thus, beta–gamma alternation may be a default baseline state in the basal ganglia. However, movement may bring out additional complexity; for instance, it has been shown recently that specific movement tasks enhance activity in a higher-frequency band within gamma (Berke, 2009).
Whereas we found the nesting frequency to center in the delta band, studies of nesting of gamma in hippocampus and cortical networks emphasize the role played by slightly higher (theta) frequencies (Bragin et al., 1995; Chrobak and Buzsáki, 1998; Lakatos et al., 2005; Canolty et al., 2006; Sirota et al., 2008). However, similar to the present finding, activity at lower frequencies in the delta range has also been shown to form the base of a complex “oscillatory hierarchy” in auditory cortex of conscious monkeys, in which delta oscillations modulate theta, which in turn modulates gamma band activity (Lakatos et al., 2005).
It has been suggested that the function of interfrequency nesting recorded elsewhere in the brain may be to allow spatially distributed subpopulations to form temporally defined ensembles through the action of propagating slow oscillations (Chrobak and Buzsáki, 1998) or to modulate excitability to optimize processing of rhythmic inputs (Lakatos et al., 2005). In basal ganglia, an additional possible functional role for nesting might be to organize interfrequency phase–amplitude relationships. We found that beta and gamma fluctuated in amplitude in anti-phase relative to specific phases of delta waves so that beta and gamma amplitude were locked in a delicate temporal counterpoint. Such conjoint nesting of dual frequencies does not appear to have been reported previously and presumably reflects the special importance of beta–gamma interactions in basal ganglia functions (Brown, 2003; Gatev et al., 2006; van der Meer and Redish, 2009). Propagating delta waves therefore appear to lie at the base of an oscillatory hierarchy that may be critical for organizing functional states within the basal ganglia–cortex loop.
Dopamine dependence of frequency power interactions: different roles of D1 and D2 receptors
Activity through basal ganglia pathways is strongly modulated by dopamine acting within the striatum and also to some extent directly in STN and GP (Smith et al., 1998; Gauthier et al., 1999). During hypokinesia induced by blockade of endogenous dopamine action at either D1 or D2 receptors with SCH23390 or raclopride, respectively, we found changes in multiple LFP measures, including power in different frequency bands, the correlation in power between bands, and the nesting of this correlation according to delta phase. Changes in neural circuit activity could in principle be secondary to, rather than causative of, changes in motility. However, this is unlikely to account for the findings here, because the two drugs had differential effects on power and interfrequency correlations, despite having similar effects on movement.
Given the difference in relative abundance of the two receptor subtypes in the direct (dominated by D1) and indirect (dominated by D2) pathways, differential effects of blockade of these receptors might indicate differing functional roles of the two pathways in the modulations seen in LFP recorded from GP. Effects on beta frequency were relatively straightforward. Both D1 and D2 antagonists acutely increased low-beta power in the GP. This is consistent with some previous acute studies (Sebban et al., 1999; Dejean et al., 2009) and with the idea that beta power is linked to akinetic states, although increases in beta power appear to be much more robust in chronic lesion rat models of Parkinson's disease (Sharott et al., 2005).
Furthermore, both antagonists abolished the entrainment of beta power fluctuation to the phase of δ3. Given the fact that the negative correlation observed between beta and gamma power appears to reflect an inverse relationship of the two frequency bands to delta, it is not surprising that the abolition of nesting of beta was accompanied by elimination of the negative correlation between beta and gamma frequencies that is observed in control conditions.
Although the underlying mechanisms relating beta power to delta phase remain uncertain, these data raise the possibility that the loss of entrainment to delta phase could be causally linked to the increase in beta power that was seen. Decoupling beta power from a temporal constraint provided by delta rhythms could enable beta power to remain at abnormally high levels for protracted periods and thus contribute to the exaggeration of beta oscillations seen in virtually all structures of the basal ganglia–cortex loop in parkinsonism (Hammond et al., 2007). Whatever the case, because beta power and beta–delta nesting in GP are equally affected by both dopamine receptor subtypes, these parameters are unlikely to be specifically modulated by indirect pathway activity limited to GP but may be more generally related to changes in capacity for movement, generated elsewhere in the basal ganglia–cortex circuit after disruption of dopamine transmission.
In contrast to the nonspecific effects of dopamine antagonists on beta power, we found that power in the gamma band was only affected significantly after administration of the D2 antagonist raclopride, with no change after SCH23390, despite both drugs resulting in similar levels of akinesia. Gamma recorded from the STN, which is strongly interconnected with GP in the indirect pathway (Smith et al., 1998; Bevan et al., 2007), is also modulated by D2 receptor (Brown et al., 2002). Together, these data support the idea that gamma in GP reflects a specific function of local circuit interactions within the indirect pathway (Bevan et al., 2002), as suggested by basal ganglia simulations (Humphries et al., 2006). Constitutive dopamine acting on the D2 receptor tends to decrease neural excitability (Surmeier et al., 2007), so raclopride is expected to increase the excitability of indirect pathway neurons projecting to GP. The resulting increased inhibitory input to GP could be responsible for lowering gamma power, for instance, by modulating local circuit interactions between GP and STN involved in generating local gamma (Humphries et al., 2006).
Although raclopride reduced average gamma power, our analysis revealed that mechanisms linking gamma power fluctuations with delta phase remained intact, because the remaining gamma activity was still nested within delta. These data suggest that the nesting mechanism does not depend on dopamine acting via D2-receptor-bearing elements of the indirect pathway. Instead, D1 antagonism had a dramatic effect, reversing the phase of delta at which gamma was at maximum, indicating that D1-receptor-expressing direct pathway neurons may be more involved in regulating the nesting phenomenon. Direct pathway neurons could modulate the gamma in GP indirectly, by altering thalamocortical loop activity that in turn feeds back into striatal and STN networks. More directly, some D1-expressing striatal neurons have axon collaterals that innervate the GP, which could play a more local role in this switching (Kawaguchi et al., 1990; Wu et al., 2000; Deng et al., 2006).
Footnotes
This work was supported by the New Zealand Neurological Foundation and Health Research Council. We thank Aaron Sheerin and Roseanna Smither for technical assistance.
- Correspondence should be addressed to Assoc. Prof. Brian Hyland, Department of Physiology, School of Medical Sciences, University of Otago, P.O. Box 913, Dunedin 9054, New Zealand. brian.hyland{at}otago.ac.nz