|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
INTEGRATIVE |
1 Medical Research Council Anatomical Neuropharmacology Unit, University of Oxford, Oxford, UK
2 Sobell Department of Motor Neuroscience and Movement Disorders, Institute of Neurology, London, UK
| Abstract |
|---|
|
|
|---|
(Received 29 March 2006;
accepted after revision 12 May 2006;
first published online 18 May 2006)
Corresponding author P. J. Magill: Medical Research Council Anatomical Neuropharmacology Unit, University of Oxford, Mansfield Road, Oxford OX1 3TH, UK. Email: peter.magill{at}pharm.ox.ac.uk
| Introduction |
|---|
|
|
|---|
10 Hz) may be mechanistically important for information processing across different levels of the sensory and motor systems (MacKay, 1997; Engel & Singer, 2001; Engel et al. 2001; Buzsáki & Draguhn, 2004). Several studies have shown that oscillatory activity in the cortex and BG may be temporally coupled (Boraud et al. 2005). However, although the temporal coupling of oscillations suggests a synchronization of activity in cortex and BG, it provides little insight into the nature of any transformation of information that might be achieved, and whether synchronization per se is an integral part of this transformation, as has been suggested in the visual cortex (Singer, 1999). The observations that deranged synchronization in cortex and BG is common in disorders with prominent motor dysfunction, such as Parkinson's disease (PD) and its animal models (Bergman et al. 1998; Levy et al. 2000, 2001, 2002; Brown et al. 2001; Bevan et al. 2002b; Boraud et al. 2005; Sharott et al. 2005b; Fogelson et al. 2006), and that synchronization varies in a task-dependent fashion (Cassidy et al. 2002; Williams et al. 2003; Kühn et al. 2004), provide circumstantial evidence linking synchronization of activity within cortico-basal ganglia-thalamocortical circuits to information processing. One of the major patterns of synchronization in PD, that in the beta-frequency (1530 Hz) band, is driven from cortex, at least in the parkinsonian state at rest (Marsden et al. 2001; Williams et al. 2002; Fogelson et al. 2006). As such, an alternative interpretation is that the temporal coupling of cortical and BG activity is an epiphenomenon of cortical synchronization in the same frequency range, similar to the manifestation of the beta rhythm in motor unit discharges that probably arises simply as the product of corticospinal connectivity (Brown, 2000). Moreover, the temporal coupling of activity that is evident in BG local field potentials and cortical electroencephalogram recorded in PD patients must involve synchronization across large assemblies of neurons in both BG and cortex (Brown, 2003; Fogelson et al. 2006), and is therefore unlikely to contain signals related to specific channels of information processing. Nevertheless, some synchronization of neuronal activity does occur, particularly in the beta-frequency band, in the BG of nonparkinsonian animals and humans, raising the possibility that these oscillations are linked to information processing under physiological rather than pathological conditions (Courtemanche et al. 2003; Sochurkova & Rektor, 2003; Berke et al. 2004; Magill et al. 2004b; Sharott et al. 2005a). If relevant for cortico-basal ganglia-thalamocortical function, one might expect patterns of synchronization to change according to local input, thereby retaining information specific to connectivity, and that, when appropriate, the basal ganglia would lead cortical activity rather than be driven by cortex.
One way to address this issue is to examine the effects of brief cortical stimulation on the relative timing (i.e. phase relationships) and topography of temporally coupled oscillations in subthalamic nucleus (STN) and cortex. Studies employing electrical stimulation of cortex have provided valuable insight into the integration of corticosubthalamic and corticostriatal inputs in the BG (Kita, 1994; Nambu et al. 2002). Stimulation of ipsilateral frontal and prefrontal cortices often leads to a short-latency multiphasic response in STN neurons, which is initiated by excitation that is driven by the corticosubthalamic pathway (Kitai & Deniau, 1981; Ryan & Clark, 1992; Fujimoto & Kita, 1993; Maurice et al. 1998; Nambu et al. 2000; Kolomiets et al. 2001; Magill et al. 2004a). The first three phases of the response, which evolve rapidly (within 50 ms), are usually followed by a sustained reduction in firing of
200 ms in duration, which has been proposed to be due to disfacilitation of cortical inputs (Fujimoto & Kita, 1993; Nambu et al. 2000). The characteristics and functional relevance of the delayed activity that inevitably marks the end of this prolonged reduction in firing are yet to be determined. As such, it is important to determine whether the delayed responses of STN neurons contain oscillatory components that are related to ongoing rhythmic activity in cortex, or how any such relationship is in turn governed by the topography of the underlying circuits.
To address these issues and further elucidate the temporal coupling of oscillatory activity in STN and cortex, we compared the delayed (>200 ms) responses of STN units and cortical population activity, as recorded in the frontal electrocorticogram (ECoG), to cortical stimulation in anaesthetized rats. Delayed responses were analysed in the frequency domain and activity relationships were defined in terms of coherence and phase (Halliday et al. 1995). The dependence of activity relationships on topography was characterized by recording unit activity along the rostrocaudal axis of STN (Magill et al. 2004a), and by varying the stimulation site between (sensorimotor) frontal cortex, which is known to project directly to STN, and (auditory) temporal cortex, which does not contribute to the corticosubthalamic pathways (Canteras et al. 1990; Kolomiets et al. 2001).
| Methods |
|---|
|
|
|---|
Experimental procedures were carried out on adult male Sprague-Dawley rats (Charles River, Margate, UK) and were conducted in accordance with the UK Animals (Scientific Procedures) Act 1986.
Electrophysiological recordings were made in seven rats (230350 g). Anaesthesia was induced with isoflurane (Isoflo; Schering-Plough Ltd, Welwyn Garden City, UK) and maintained with urethane (1.3 g kg1
I.P., ethyl carbamate; Sigma, Poole, UK), and supplemental doses of ketamine (30 mg kg1
I.P., Ketaset; Willows Francis, Crawley, UK) and xylazine (3 mg kg1
I.P., Rompun; Bayer, Germany), as previously described (Magill et al. 2001). All wound margins were infiltrated with the local anaesthetic bupivacaine (0.75% w/v; Astra, Kings Langley, UK), and corneal dehydration was prevented with application of Hypromellose eye drops (Norton Pharmaceuticals Ltd, Harlow, UK). Animals were then placed in a stereotaxic frame. Body temperature was maintained at 37 ± 0.5°C with the use of a homeothermic heating device (Harvard Apparatus Ltd, Edenbridge, UK). Anaesthesia levels were assessed by testing reflexes to a cutaneous pinch or gentle corneal stimulation, and by examination of the frontal ECoG (see following text). In terms of using the ECoG to ascertain anaesthesia levels, the expression levels of a slow oscillation at
1 Hz and spindle oscillations at
10 Hz, which together constitute the major components of the slow-wave activity in this preparation (see Magill et al. 2004b), were assessed online; deeper anaesthesia is accompanied by higher levels of slow-wave activity. Electrocardiographic (ECG) activity and respiration rate were also monitored constantly to ensure the animals' well being (see following text). Mineral oil or saline solution (0.9% w/v NaCl) was applied to all areas of exposed cortex to prevent dehydration.
Parallel bipolar stimulating electrodes (constructed from nylon-coated stainless steel wires; California Fine Wire, Grover City, CA, USA), with tip diameters of
100 µm, a tip separation of
150 µm, and an impedance of
10 k
, were implanted into the right frontal cortex and right temporal cortex, and then affixed to the skull with dental acrylic cement (Associated Dental Products Ltd, Swindon, UK). The coordinates of the frontal stimulation site (AP, +4.2 mm; ML, 3.5 mm (bregma reference); and a depth of 2.3 mm below the dura; Paxinos & Watson, 1986) correspond to the border region between the lateral and medial agranular fields of the somatic sensorimotor cortex (Donoghue & Wise, 1982; Donoghue & Parham, 1983). The coordinates of the temporal stimulation site (AP, 5.2 mm; ML, 6.9 mm; and a depth of 2.3 mm below the dura) approximately correspond to the primary auditory cortex (see Kolomiets et al. 2001).
The ECoG was recorded via a 1 mm diameter steel screw juxtaposed to the dura mater above the right frontal cortex (AP, +4.5 mm; ML, 2.0 mm; which corresponds to the medial agranular field of the somatic sensorimotor cortex; Donoghue & Wise, 1982) and referenced against an indifferent electrode placed adjacent to the temporal musculature. This and adjacent regions of cortex project to both striatum and STN, as demonstrated in anatomical and electrophysiological studies (Fujimoto & Kita, 1993; Smith et al. 1998; Kolomiets et al. 2001; Magill et al. 2004a), and, thus, activity in these regions is of direct functional relevance. Raw ECoG was band-pass filtered (0.12000 Hz, 3 dB limits) and amplified (x2000, NL104 preamplifier; Digitimer Ltd, Welwyn Garden City, UK) before acquisition. The ECG was differentially recorded via two silver wires that were inserted subcutaneously into the ipsilateral forelimb and hindlimb. Raw ECG was band-pass filtered (10100 Hz) and amplified (x5000, NL104; Digitimer) before acquisition. The chest movements accompanying respiration were recorded using a miniature accelerometer (AP19; Bay Systems Ltd, Somerset, UK) and charge amplifier (Type 5007; Kistler Instrumente AG, Winterthur, Switzerland). The signal from the accelerometer allowed the depth and rate of respiration to be accurately assessed on- and off-line (Magill et al. 2004b).
Extracellular recordings of action potentials in the ipsilateral STN were made with glass electrodes (612 M
measured at 10 Hz in situ, tip diameters of 2.03.0 µm) that were filled with a 0.5 M NaCl solution containing 1.5% w/v Neurobiotin (Vector Laboratories, Peterborough, UK). Electrodes were lowered into the brain using a computer-controlled stepper motor (Burleigh IW-711; Scientifica Ltd, Harpenden, UK), which allowed the electrode depth to be determined with a resolution of 0.5 µm. Extracellular signals from the electrode were amplified (x10) through the active bridge circuit of an Axoprobe 1A amplifier (Axon Instruments, Union City, CA, USA), AC-coupled and amplified a further x100 (NL106; Digitimer), before band-pass filtering (between 0.4 and 4 kHz; NL125, Digitimer). Signals were digitized on-line (see below) and simultaneously displayed on an oscilloscope (DSO 610; Gould Instruments, Ilford, Essex, UK). Action potentials were typically between 0.4 and 1.2 mV in amplitude and always exhibited an initial positive deflection (Magill et al. 2000, 2001).
The responses of STN neurons to cortical inputs were determined by focal electrical stimulation of the cortex (Maurice et al. 1998; Kolomiets et al. 2001; Magill et al. 2004a). Note that the first three phases of the short-latency responses to cortical stimulation have been previously reported in the same animals (Magill et al. 2004a). Electrical stimuli, which consisted of single square-wave current pulses of 0.3 ms duration and variable amplitude (300600 µA to ensure maximal responses; see Magill et al. 2004a), were delivered to the ipsilateral frontal or temporal cortices at a frequency of 0.67 Hz using a constant current isolator (A360D; World Precision Instruments Ltd, Stevenage, UK) that was gated by a programmable pulse generator (Master-8; AMPI, Jerusalem, Israel). Between 210 and 230 cortical stimuli were delivered for each STN neuron tested. Stimulation of frontal cortex at this and similar frequencies (0.41.4 Hz) evokes robust and stable single-cell and population responses in STN (Maurice et al. 1998; Nambu et al. 2000; Kolomiets et al. 2001; Magill et al. 2004a). However, we cannot rule out the possibility that stimulation at 0.67 Hz may have induced synaptic plasticity in the cortex and/or STN (or related neural circuits), and thus, it would be of interest to test in the future whether the responses described hereafter vary with stimulation parameters.
Between two and five electrode tracks were performed in each animal. The final recording location in each experiment was marked by a discrete extracellular deposit of Neurobiotin (100 nA anodal current; 1 s (50%) duty cycle for 60 min; Magill et al. 2004b). Following a period of 12 h for the uptake and transport of the Neurobiotin by neurons and glia at the recording site, animals were given a lethal dose of ketamine anaesthetic (200 mg kg1, I.P.) and perfused via the ascending aorta with 100 ml of 0.01 M phosphate-buffered saline (PBS) at pH 7.4, followed by 300 ml of 4% w/v paraformaldehyde and 0.1% w/v glutaraldehyde in 0.1 M phosphate buffer, pH 7.4, and then by 150 ml of the same solution without glutaraldehyde. Brains were then postfixed in the latter solution at 4°C for at least 12 h before sectioning.
Histochemistry
Standard techniques were used to visualize the Neurobiotin deposits (see Magill et al. 2001). Briefly, the fixed brain was cut into 60 µm thick sections in the coronal plane on a vibrating blade microtome (VT1000S; Leica Microsystems, Milton Keynes, UK). Sections were washed in PBS and incubated overnight in avidin-biotin-peroxidase complex (1:100; Vector) in PBS containing 0.2% v/v Triton X-100 and 1% w/v bovine serum albumin (Sigma). After washing, the sections were incubated in hydrogen peroxide (0.002% w/v; Sigma) and diaminobenzidine tetrahydrochloride (0.025% w/v; Sigma) in the presence of nickel ammonium sulphate (0.5% w/v; Sigma) dissolved in Tris buffer (0.05 M, pH 8.0) for 1530 min. Neurobiotin-filled cells were intensely labelled with an insoluble, black/blue precipitate (Fig. 1A and B; see also Magill et al. 2004b). Finally, sections were dehydrated, cleared and mounted for light microscopy using standard techniques (Bolam, 1992).
|
Data acquisition
Unit activity was sampled at 10 kHz. The ECoG signal was sampled at 5 kHz. The ECG and respiration signals were each sampled at 400 Hz. All biopotentials were digitized on-line with a PC running Spike2 acquisition and analysis software (version 4; Cambridge Electronic Design Ltd, Cambridge, UK).
Peristimulus time histograms and cusums
The short-latency (<200 ms) responses of units to cortical stimulation were analysed by peristimulus time histograms (PSTHs). These were constructed for each STN unit from 210 consecutive stimulation trials (Magill et al. 2004a). The timings of changes in single unit activity were determined by change-point analysis, using commercial software (Change-Point Analyser 2.0 shareware program, http://www.variation.com/cpa/tech/changepoint.html), and techniques previously described (Cassidy et al. 2002; Williams et al. 2003). Change-point analysis iteratively uses a combination of time-varying cumulative sum charts (cusums; see Ellaway, 1978) and random re-ordering of the original data (bootstrapping) to set confidence limits for the cusum of the original data set (Change-Point Analyser 2.0 shareware program). PSTHs of unit activity from 200 ms before, to 400 ms after cortical stimuli were first plotted with 1 ms bins. Cusums were then determined by plotting the sequentially summed deviations from the average discharge count per millisecond determined from the whole record segment (total of 600 ms), as this affords a more conservative estimate of significant changes than the use of any a priori defined control period. One thousand bootstraps were performed in each test and only changes with probabilities of
0.05 were considered significant. Upward gradients in cusums denote increases in discharge rate, whereas downward gradients denote decreases in discharge rate. The mean discharge count per millisecond over each identified interval of significant change was calculated and expressed as a percentage of the mean prestimulus (background) discharge count per millisecond, which was estimated from 200 to 0 ms before cortical stimulation. A reduction in firing was considered complete if mean activity over the period of reduced discharge identified by change point analysis was <5% of background activity.
Time-domain representations of coupling between STN unit activity and ECoG were determined by histograms using only the first action potential fired in the series of delayed oscillatory STN unit discharges that followed the prolonged period of reduced activity, or a long pause in firing, while at the same time reducing the ECoG signal to a digital representation of consecutive positive peaks. The latter was achieved by automatically selecting all ECoG peaks of an amplitude greater than 0.1 mV, with a minimum interval of 3 ms, using the peak-detection function of Spike2 after the ECoG had been first smoothed using a moving-average filter that averaged five consecutive data points before shifting one data point and repeating the procedure. The ECoG data were searched for all peaks that were followed by a fall in voltage of at least 0.1 mV. The criterion for detection was thus a peak-to-trough amplitude of 0.1 mV. Only positive peaks were considered for the analysis. Spike-triggered averages (STAs) of ECoGs were also calculated, after first discriminating STN unit discharges to give a digital signal.
Coherence and phase analyses
The long-latency, delayed (>200 ms) effects of cortical stimulation were characterized by frequency analysis of point process data (STN unit discharge) and time series data (cortical activity). For example applications of this technique to similar signals in neurophysiological studies see Brown et al. (1999) and Halliday et al. (1999). Point process unit activity in STN was downsampled to 5 kHz to match the ECoG. Two periods of unit activity and simultaneously recorded ECoG of 410 ms duration were exported for frequency analysis: 205615 ms and 10241434 ms after cortical stimulation. These 410 ms windows were the minimum that afforded reasonable frequency resolution (2.4 Hz) in the coherence and phase analyses. The first window, with an onset of 205 ms after stimulation, approximated, on the whole, to the timing and duration of robust oscillatory responses in the cortex and subthalamic nucleus. The latter window, with onset 1024 ms after stimulation, was chosen to test whether oscillatory activity and corresponding coherence was short-lived and recovered prior to the next electrical stimulus in the train. Both periods so chosen involved response elements with local stationarity in the STN and ECoG signals. Unit activity and ECoG were examined after 210 cortical stimuli, so this gave 86.1 s of data for each of these periods. Time-series data, describing STN unit activity and ECoG activity during the 205615 ms and 10241434 ms epochs, were separately concatenated and analysed in the frequency domain using blocks of 2048 data points (410 ms), which exactly corresponded to the exported sections. Data were passed through a Hanning window filter so that discontinuities between successive concatenated blocks were suppressed. Contiguous data sections of 1028 s (mean duration 16 s) were also drawn from recordings of spontaneous activity made immediately prior to the start of the train of 210 cortical stimuli in the 15 STN units that exhibited responses to corticosubthalamic input and had >10 s of spontaneous activity recorded before stimulation. As the duration of recording of spontaneous activity varied between units and never reached 86.1 s, spontaneous activity was separately compared to stimulation periods of matching duration for each unit.
The ECoG, denoted by A, was assumed to be a realization of a stationary, zero-mean time series. Unit activity in STN was assumed to be a realization of a stationary, stochastic point process and was denoted by B (Perkel et al. 1967). Both processes were further assumed to satisfy a mixing condition, whereby sample values widely separated in time were independent (Brillinger, 1981). The principal statistical tool used for the data analyses was the discrete Fourier transform and parameters derived from it, all of which were estimated by dividing the records into a number of disjointed sections of equal duration (2048 data points, affording a frequency resolution of 2.4 Hz). Spectra were estimated by averaging across these discrete sections (Halliday et al. 1995). In the frequency domain, estimates of the power spectra of the ECoG, fAA(
), and unit activity, fBB(
), were constructed, along with estimates of coherence, |RAB(
)|2, between the ECoG and unit activity. The coherence was calculated according to:
|
|
denotes the frequency, and fAB(
) is the cross spectrum between A and B. Confidence levels (95% CL) for coherence spectra were determined according to methods described in Halliday et al. (1995). A peak of significant coherence in the spectrum was defined as >2 contiguous frequency bins (i.e. >4.8 Hz in duration) above the 95% CL, provided that the onset frequency of the peak was >4.8 Hz. The latter criterion was included so as to distinguish coherence due to slow shifts in the DC potential from discrete peaks of coherence at higher frequencies.
Phase,
AB, defined as the argument of the cross spectrum, is:
|
|
The phase estimated from a single point (e.g. frequency) is ambiguous. Measuring phase relationships that are linear over a band of frequencies reduces this ambiguity (Gotman, 1983). Under these circumstances, the temporal delay between the signals can be calculated from delay (ms) =
g
x (1/2
), where g is the gradient of the line in plots of phase (radians) against frequency (Hz). A negative gradient indicates that the input/reference signal leads. The constant time lag between the STN unit discharge and ECoG could be calculated from the slope of the phase estimate after a line had been fitted by linear regression in seven spectra from seven units, from which the corresponding coherence estimates had a total of eight significant peaks (one spectrum had two peaks; coherence determined as significant where it exceeded the 95% CL as determined above). These peaks had >4 contiguous frequency bins and linear regression of the phase estimates established a r2 > 0.65 and a P < 0.05.
Statistics
Linear correlation was tested by Pearson's correlation coefficient (unless otherwise stated). All t tests, Fisher Exact tests,
2 tests, and tests of the significance of correlation coefficients were two tailed. Significance was considered as P < 0.05.
| Results |
|---|
|
|
|---|
10 Hz, but with a variety of firing patterns. The spontaneous activity of STN neurons and their responses to cortical stimulation were recorded during the whole gamut of brain states commonly observed in the ECoG of urethane-anaesthetized rats, including slow-wave activity, spontaneous (global) activation, and mixtures thereof (see Magill et al. 2000, 2004b). Accordingly, STN neurons tended to fire bursts of spikes in a periodic manner (
1 Hz) during slow-wave activity, whereas they tended to discharge with tonic, irregular patterns during the activated brain state (Magill et al. 2000). As reported previously, the responses of STN neurons to cortical stimulation are often multiphasic (Kitai & Deniau, 1981; Ryan & Clark, 1992; Fujimoto & Kita, 1993; Maurice et al. 1998; Nambu et al. 2000; Kolomiets et al. 2001; Magill et al. 2004a). Such multiphasic responses can include, in order of temporal evolution, a short-latency excitation, a brief reduction in firing, a secondary excitation (a combination hereafter referred to as a short-latency triphasic response), a prolonged reduction in firing and, finally, a delayed recovery of activity. Because the short-latency triphasic responses of STN units to electrical stimulation of the ipsilateral frontal cortex, as described here, have been previously reported in the same animals (Magill et al. 2004a), and given that these early response phases have also been described in detail by others (see above), they will only be briefly summarized below.
Responses of units recorded in rostral STN to stimulation of frontal cortex
Peristimulus time histograms of the responses of 24 units in rostral STN (Figs 1C and 2A) to electrical stimulation of the ipsilateral frontal cortex demonstrated a short-latency triphasic pattern, beginning with an excitation of a mean latency of 4.3 ms (S.D. 2.0 ms) that was followed by a brief period of suppressed firing and then another excitation. This short-latency triphasic response was in turn followed at 33 ms (S.D. 19 ms) by a prolonged silencing of unit discharge. Thereafter, upon recovery from this prolonged quiescence at around 200 ms, unit discharges had a tendency to occur at periodic intervals for up to 650 ms after stimulation, evident over the 200650 ms period as a series of peaks and troughs in PSTHs (Fig. 2A). These were mirrored by undulations in the frontal cortical ECoG over the same period but tended to dampen thereafter (Fig. 2), such that ECoG peaks (see Methods) were only very rarely detected between 650 and 1450 ms after stimulation. At the range of stimulation current amplitudes used here, which ensured maximal responses (see Methods), these delayed unit responses and ECoG undulations did not systematically differ with gross fluctuations in brain state (as assessed from the ECoG), in agreement with our previous finding that short-latency triphasic responses evoked during slow-wave activity did not differ substantially from those evoked during global activation (see Magill et al. 2004a). However, although these delayed responses were observed in recordings made during either slow-wave activity or global activation, the two extremes of brain state observed in this preparation (Magill et al. 2004b), we were unable to make direct comparisons of the evoked responses of any one rostral STN neuron, together with those of a single cortical locus, during both brain states.
|
|
|
|
|
To investigate the topography of frequency relationships between cortex and STN, we also recorded the responses of neurons in the caudal half of STN (n
= 13; Fig. 1D). The responses of neurons located in rostral and caudal STN to stimulation of the frontal cortex were different (see Magill et al. 2004a). No short-latency triphasic responses were seen in the PSTHs of caudal STN units (Fig. 2B). A delayed and prolonged reduction in firing was present in all but one caudal STN neuron, but the reduction was less than that displayed by rostral STN units (mean discharge frequency of caudal STN units fell to 73% (S.D. 18%) of baseline, as compared to complete quiescence for rostral units). The mean onset latency of this reduced firing, as determined from cusums, was 37 ms (S.D. 11 ms), which was not significantly different from the mean onset latency of the delayed pause in activity exhibited by rostral STN units (unpaired t test). The mean durations of the sustained periods of reduced firing, as determined from cusums, were also similar (unpaired t test; 236 ms (S.D. 86 ms) for caudal units and 232 ms (S.D. 54 ms) for rostral units). The appearance of the mean ECoG waveform after frontal cortical stimulation did not differ in experiments in which units were recorded in rostral and caudal STN (Fig. 2A and B). In particular, the same undulations were evident in the frontal ECoG at 200650 ms after cortical stimuli (Fig. 2B). Despite these apparent similarities, the tendency for the discharges of caudal STN units to occur at periodic intervals immediately following recovery from reduced discharge was less marked as compared to rostral STN units (cf. Figure 2A and B). Accordingly, spectra of the coherence between unit discharges in caudal STN and frontal ECoG over the period of 205615 ms after stimulation demonstrated fewer significant peaks, particularly at those frequencies above 15 Hz (Fig. 5C). Comparison of the coherence spectra of rostral STN units and caudal STN units revealed that there was no difference in the incidence of significant coherence below 15 Hz, but there was a significant difference over the 1530 Hz range (
2 test, P < 0.00001), such that rostral STN units exhibited a proportionally greater incidence of coherence above 15 Hz.
Responses of units recorded in rostral STN to stimulation of temporal cortex
To determine whether the coherence profiles over the period of 205615 ms after stimulation were relatively specific to stimulation of the frontal cortex, we also recorded the responses of some of the same rostral STN units (n
= 12 of 24) to stimulation of the ipsilateral temporal cortex. Stimulation of temporal cortex did not produce short-latency triphasic responses in PSTHs of rostral STN units, and only one neuron demonstrated the delayed and prolonged reduction in firing. The grouping of unit discharges at periodic intervals over the 200650 ms period after stimulation were diminished as compared to responses to frontal cortical stimulation, and the delayed undulations were absent from the frontal ECoG. Accordingly, there was a significant reduction (
2 test, P < 0.00001) in coherence between STN unit activity and frontal ECoG over 530 Hz after temporal cortical stimulation as compared to frontal stimulation (Fig. 5D). Furthermore, phase relationships of coherence over 530 Hz were mixed after temporal cortex stimulation, so that phase spectra could not be fitted by a linear model. Thus, electrical stimulation of the temporal cortex was relatively ineffective in eliciting coherence between unit activity in the rostral STN and frontal ECoG.
Responses of units recorded in caudal STN to stimulation of temporal cortex
To further define the functional specificity of evoked coherence, we recorded the responses of some of the same units in caudal STN (n
= 9 of 13) to stimulation of ipsilateral temporal cortex. Stimulation of temporal cortex did not elicit short-latency triphasic responses in the PSTHs of caudal STN units. However, all but one neuron demonstrated a delayed reduction in firing in response to stimulation. The mean onset latency of this reduction, as determined from cusums, was 44 ms (S.D. 24 ms), which was not different from the onset latency of the delayed reduction in firing in caudal STN units following frontal cortex stimulation (unpaired t test). The mean duration of the periods of reduced firing was also similar (279 ms (S.D. 69 ms); unpaired t test). The mean discharge frequency of caudal STN units fell to 66% (S.D. 12%) of baseline levels during these delayed reductions in firing. The tendency for unit discharges in the caudal STN to occur at periodic intervals over the 200650 ms period after stimulation of the temporal cortex was not as clear as in the rostral STN following stimulation of the frontal cortex. Undulations were virtually absent in the frontal ECoG following stimulation of temporal cortex. Coherence between unit activity in caudal STN and ECoG in the period of 205615 ms after stimulation did not differ to that seen following frontal cortex stimulation over the 530 Hz band (
2 test, P > 0.05), as illustrated in Fig. 5E. Thus, electrical stimulation of the temporal and frontal cortices was equally (in)effective in eliciting coherence between frontal ECoG and unit activity in the caudal STN. However, coherence at 530 Hz was still significantly less than that recorded in rostral STN following frontal cortex stimulation (
2 test, P < 0.0001).
Correlation of the coherence between STN units and frontal ECoG at 515 Hz with the delayed reductions in firing rate
Coherence between unit activity and ECoG at frequencies above 15 Hz was thus largely restricted to those STN neurons exhibiting a typical, short-latency triphasic response to frontal cortical stimulation, including a brief short-latency excitation response indicative of direct corticosubthalamic input arising from the region of stimulated cortex. Coherence between STN units and ECoG below 15 Hz was more widely distributed in STN, as was the case for the delayed and sustained reductions in firing. Coherence at 515 Hz was elicited following stimulation of either the frontal or temporal cortices, although the strongest coherence tended to occur in the rostral STN after frontal cortical stimulation. Accordingly, we correlated the degree of significant coherence in the 515 Hz band with the magnitude of the delayed reductions in firing elicited by frontal and temporal cortical stimulation in all 58 recordings (including those in which the same unit was studied after stimulation of the frontal and temporal cortex). Coherence was greater when the reduction in unit activity within STN was greater (Spearman's rho 0.481, P < 0.001). This correlation was not significant when only units in rostral STN were considered (because most units were recorded following stimulation of the frontal cortex, when the reduction in firing was complete and neurons were silent). However, when only units in caudal STN were considered, coherence was again greater when the reduction in unit activity within STN was greater in magnitude (Spearman's rho 0.478, P < 0.01).
| Discussion |
|---|
|
|
|---|
Synaptic inputs and intrinsic mechanisms underlying the short-latency and delayed responses of STN neurons to stimulation
The STN receives monosynaptic inputs from select areas of the ipsilateral cerebral cortex, including prefrontal, premotor, primary motor and somatosensory areas (see review by Smith et al. 1998). Neurons in the rostral half of STN typically responded to stimulation of the ipsilateral frontal (sensorimotor) cortex with a short-latency triphasic pattern of activity, as previously described in the same paradigm (Magill et al. 2004a). The first phase of the short-latency triphasic response was a brief, but powerful, excitation, which appears to be driven by the direct, excitatory corticosubthalamic projection (Kitai & Deniau, 1981; Ryan & Clark, 1992; Fujimoto & Kita, 1993; Maurice et al. 1998; Nambu et al. 2000; Kolomiets et al. 2001). The second and third phases of the typical unit response develop from polysynaptic interactions (see Fujimoto & Kita, 1993; Maurice et al. 1998; Bevan et al. 2002b). The accompanying short-latency responses of the ECoG to frontal cortical stimulation were almost certainly due to direct orthodromic and antidromic excitation of intracortical and corticopetal neuronal elements (Oldfield & Castro-Alamancos, 2003).
Units recorded in the rostral half of STN did not display a short-latency triphasic response to stimulation of ipsilateral temporal cortex, which is in keeping with studies demonstrating that the (auditory) temporal cortex does not project directly to STN (Canteras et al. 1990; Kolomiets et al. 2001). The topographical organization of unit responses was not only evident with respect to cortical areas, but also within the STN itself. Neurons located in the caudal half of STN did not respond at short latency to stimulation of either frontal or temporal cortex. These findings are in agreement with previous anatomical (Canteras et al. 1990; Kolomiets et al. 2001) and physiological (Kolomiets et al. 2001; Magill et al. 2004a) studies, which have shown that the frontal cortical areas stimulated in the present study project to the rostral two-thirds of STN only. Importantly, response differences were unlikely to be due to insufficient current flow at the site of stimulation because responses should have been maximal, with respect to pattern, at the stimulus intensities used (Magill et al. 2004a).
A prolonged reduction in activity, albeit of varying intensity, frequently preceded the delayed oscillatory responses of neurons in both rostral and caudal aspects of STN, which were, in turn, associated with increased coherence with cortex in the 530 Hz range. With one exception, neurons in rostral STN did not demonstrate a sustained reduction in activity after stimulation of temporal cortex, and their delayed oscillatory activity in response to the same stimulus was only weak (and weakly coherent with cortex at 515 Hz). The most robust unit activity oscillations and strongest coherence with cortex were characteristic of the rostral STN neurons that were consistently and completely silenced following frontal cortical stimulation. However, a prolonged reduction in the activity of STN units alone was not a requisite for the generation of subsequent coherent oscillations at 1530 Hz. Although many STN units displayed coherent oscillations after stimulation, coherence at 1530 Hz was almost exclusive to rostral STN in response to inputs from frontal cortex.
The cellular and synaptic mechanisms underlying the prolonged reduction or pause in firing are unknown, as are the factors responsible for generating the delayed oscillatory activity of STN neurons. A disfacilitation of cortical inputs has been hypothesized to be the cause of the prolonged reduction in firing (Fujimoto & Kita, 1993; Kita, 1994; Nambu et al. 2002). Yet, this interpretation does not fit well with extracellular recording data showing that many STN neurons are not silenced by acute ablation of much of the ipsilateral frontal and parietal cortices (Magill et al. 2001). The cumulative activation of Ca2+-dependent K+ channels during the third (excitation) phase of the short-latency triphasic response, eventually leading to a prolonged afterhyperpolarization, might also have contributed in some neurons (Beurrier et al. 1999; Bevan & Wilson, 1999; Hallworth et al. 2003). More recent in vitro data suggest that the delayed pause in firing could also be a reflection of an intrinsic depolarized plateau potential, during which neurons are quiescent (Kass & Mintz, 2006). Whatever the underlying mechanism(s), delayed and sustained reductions in activity were seen throughout the rostro-caudal axis of STN. Furthermore, these responses could be evoked by stimulation of either frontal or temporal cortices. Taken together, these findings suggest that the underlying circuits are less strictly organized than the highly topographic corticosubthalamic projections.
The fact that the discharge rates of most STN neurons were significantly reduced for about 200 ms after cortical stimulation indicates that they underwent prolonged membrane hyperpolarization from voltages associated with prestimulation (spontaneous) activity (but see Kass & Mintz, 2006). Although this scheme is corroborated by in vivo intracellular data (Kitai & Deniau, 1981), it should be noted that such a hyperpolarization might not be caused by inhibitory synaptic input (Fujimoto & Kita, 1993). Subthalamic nucleus neurons exhibit autonomous oscillatory activity, nonlinear inputoutput properties, and may express up to four self-sustained membrane states (Bevan & Wilson, 1999; Beurrier et al. 2000; Hallworth et al. 2003; Kass & Mintz, 2006). Hyperpolarization of the somatodendritic domains of STN neurons can reset the pacemaker oscillation of STN neurons such that, upon repeated stimulation, subsequent discharges appear to be phase-locked to the end of the hyperpolarization (Bevan et al. 2000, 2002a). Furthermore, prolonged membrane hyperpolarization may be sufficient to de-inactivate low-threshold (T-type) Ca2+ channels and/or gait hyperpolarization-activated mixed-cation (Ih) channels, both of which are non-operational during sustained single-spike activity (Beurrier et al. 1999; Bevan & Wilson, 1999; Bevan et al. 2000, 2002a). The cessation of membrane hyperpolarization may be marked by a period of rebound activity that is supported by these channels. Such rebound activity can take the form of rhythmic, single-spike firing or bursting activity, depending on the duration and magnitude of the preceding hyperpolarization (Bevan et al. 2000, 2002a). Thus, the delayed oscillatory discharges of STN neurons observed in the present study may be a function, at least in part, of the recovery from hyperpolarization and the subsequent resetting of activity. Of particular note in this regard is the possibility that rebound burst activity in STN neurons might constitute either a parallel or supplementary mechanism to the thalamocortical drive that is an important mediator of the delayed cortical response that occurs after the prolonged inhibition evoked by electrical stimulation (Grenier et al. 1998). That said, the delayed oscillations could also arise during recovery from a depolarized plateau potential (Kass & Mintz, 2006). It should be remembered, however, that although STN neurons can oscillate in the absence of synaptic inputs in vitro, it is likely that the activity of STN neurons in vivo is subject to afferent control, at least to some degree. Many globus pallidus neurons will also be active 200600 ms after cortical stimulation and the ensuing barrage of pallidal inputs may play a role in sculpting the oscillatory discharge of STN neurons (Ryan & Clark, 1991; Nambu et al. 2000; Bevan et al. 2002a,b). We also cannot rule out the possibility that the delayed oscillatory activity of STN neurons is driven, at least in part, by renewed excitatory input from, for example, the thalamus, the pedunculopontine tegmental nucleus and/or the cortex (Smith et al. 1998; Mena-Segovia et al. 2004).
Significance of temporal coupling of STN unit activity to frontal ECoG
Coherence between the discharges of single neurons in STN and the ECoG demonstrates that temporal coupling can occur, irrespective of confounding volume-conducted effects from cortex. The systematic phase differences between some STN and ECoG signals also present a strong argument against volume conduction between these sites, while arguing in favour of a true functional association. In addition, although the firing rates of STN neurons could overlap with the band of coherence (see Fig. 4), there was no linear relation between firing rates and coherence across spectra. Thus, it seems unlikely that the delayed effect of cortical stimulation was merely a resetting of periodic neuronal discharge after a sustained period of reduced firing because this mechanism would not explain why the resultant delayed activity was coherent with ECoG. Could coherence between the ECoG and STN arise without functional coupling of the two signals, if, for example, physiologically independent rhythms were reset or contemporaneously evoked at the cortical and STN levels by cortical stimulation? This too seems unlikely as it would not easily explain why one oscillatory signal (the STN unit discharge) was delayed in onset and it would require persistent amplitude- and phase-locking of the signals for up to 600 ms after cortical stimulation. Instead, the delayed and prolonged nature of the coherence between cortex and STN is highly suggestive of an ongoing functional coupling between these structures that maintains amplitude- and phase-locking over time.
It should be remembered that the present results were derived from responses to low-frequency stimulation of cortex in anaesthetized animals and thus, there is at least a possibility that our data might not be wholly representative of the single neuron and population activities present during behaviour. Nevertheless, the spontaneous activity of STN neurons recorded under urethane anaesthesia is similar to that of STN neurons recorded in the unanaesthetized state (Urbain et al. 2000). Furthermore, the short-latency responses and prolonged reductions in firing of STN neurons are not strongly dependent on ongoing brain state (Magill et al. 2004a) or anaesthetic type (cf. with Kolomiets et al. 2001). In the absence of comparative (stimulation) data from unanaesthetized animals however, we cannot rule out the possibility that the urethane anaesthesia used here influenced the delayed oscillatory responses of STN neurons. On the other hand, the significance of our findings is supported by studies suggesting that neurons in cortical areas projecting to the STN exhibit a wide range of synchronous network activity during sensory-motor integration and other complex behaviours that are likely to involve the basal ganglia as a whole (Murthy & Fetz, 1996; MacKay, 1997; Roelfsema et al. 1997; Donoghue et al. 1998; Brown, 2000; Marsden et al. 2000, 2001; Baker et al. 2001; Engel & Singer, 2001; Engel et al. 2001; Lee, 2003; Buzsáki & Draguhn, 2004). Under physiological conditions, the cortical output precipitated by electrical stimulation might be paralleled by the movement-related corollary discharge from cortex to STN (Kitai & Deniau, 1981; Giuffrida et al. 1985; Marsden et al. 2001). Rhythmic ECoG responses with frequency components similar to those described here have been observed in humans following stimulation of the motor cortex (Paus et al. 2001). Furthermore, long-latency oscillations are not a response feature of cortex that is limited to electrical stimulation. Indeed, such oscillations are known to occur in response to natural sensory inputs (Oldfield & Castro-Alamancos, 2003).
While keeping these caveats in mind, we have shown that the exact pattern of synchronization between population activity in the cortex, as reflected in the ECoG, and single-cell activity in STN depends on the type of cortical input received. For a given stimulated cortical area, the response patterns were similar in the same neuron across trials, and across different neurons recorded in the same half of the STN at different times. This suggests that changes in the temporal coupling of oscillations are a stereotyped and topographically organized response to synchronized cortical outputs. Direct excitatory input to STN from the frontal cortex via the so-called hyperdirect corticosubthalamic pathway (Nambu et al. 2002), and the ensuing cycle of excitations and inhibitions, are associated with a delayed, and temporary, periodic discharge of STN neurons. This delayed periodic discharge is coherent with oscillations in the frontal cortex over a wide band from 5 to 30 Hz. In contrast, the oscillatory responses of units in caudal STN, which did not receive corticosubthalamic input, but which did display the sustained firing reductions, are associated with a delayed and temporary coherence with ECoG activity in the 515 Hz band only. Similarly, stimulation of the temporal cortical areas that do not project to STN only resulted in minor coherence with cortex at 515 Hz. The fact that the delayed periodic discharge of STN neurons tended to lead oscillations in the frontal cortex suggests that at least some STN neurons may indirectly feed-back to, and thence influence, the synchronized activity of cortical neurons (as indexed by ECoG oscillations) following frontal cortical stimulation. Taken together, the stereotyped and organized modifications in temporal coupling between cortex and STN following synchronized cortical output suggest that these changes may be a mechanistically important feature of the response of the STN to cortical efferent activity.
| References |
|---|
|
|
|---|
Bergman H, Feingold A, Nini A, Raz A, Slovin H, Abeles M & Vaadia E (1998). Physiological aspects of information processing in the basal ganglia of normal and parkinsonian primates. Trends Neurosci 21, 3238.[CrossRef][Medline]
Berke JD, Okatan M, Skurski J & Eichenbaum HB (2004). Oscillatory entrainment of striatal neurons in freely moving rats. Neuron 43, 883896.[CrossRef][Medline]
Beurrier
C, Bioulac
B
&
Hammond
C (2000). Slowly inactivating sodium current (I(NaP) underlies single-spike activity in rat subthalamic neurons. J Neurophysiol
83, 19511957.
Beurrier
C, Congar
P, Bioulac
B
&
Hammond
C (1999). Subthalamic nucleus neurons switch from single-spike activity to burst-firing mode. J Neurosci
19, 599609.
Bevan
MD, Magill
PJ, Hallworth
NE, Bolam
JP
&
Wilson
CJ (2002a). Regulation of the timing and pattern of action potential generation in rat subthalamic neurons in vitro by GABA-A IPSPs. J Neurophysiol
87, 13481362.
Bevan MD, Magill PJ, Terman D, Bolam JP & Wilson CJ (2002b). Move to the rhythm: oscillations in the subthalamic nucleus-external globus pallidus network. Trends Neurosci 25, 525531.[CrossRef][Medline]
Bevan
MD
&
Wilson
CJ (1999). Mechanisms underlying spontaneous oscillation and rhythmic firing in rat subthalamic neurons. J Neurosci
19, 76177628.
Bevan
MD, Wilson
CJ, Bolam
JP
&
Magill
PJ (2000). Equilibrium potential of GABAA current and implications for rebound burst firing in rat subthalamic neurons in vitro. J Neurophysiol
83, 31693172.
Bolam JP (1992). Experimental Neuroanatomy: A Practical Approach. Oxford University Press, Oxford, UK.
Boraud T, Brown P, Goldberg JA, Graybiel AM & Magill PJ (2005). Oscillations in the basal ganglia: the good, the bad and the unexpected. In The Basal Ganglia VIII, ed. Bolam JP, Ingham CA & Magill PJ, pp. 324. Springer, New York.
Brillinger DR (1981). Time Series Data Analysis and Theory, 2nd edn. Holden Day, San Francisco.
Brown P (2000). Cortical drives to human muscle: the Piper and related rhythms. Prog Neurobiol 60, 97108.[CrossRef][Medline]
Brown P (2003). Oscillatory nature of human basal ganglia activity: relationship to the pathophysiology of Parkinson's disease. Mov Disord 18, 357363.[CrossRef][Medline]
Brown
P, Farmer
SF, Halliday
DM, Marsden
J
&
Rosenberg
JR (1999). Coherent cortical and muscle discharge in cortical myoclonus. Brain
122, 461472.
Brown
P, Oliviero
A, Mazzone
P, Insola
A, Tonali
P
&
Di Lazzaro
V (2001). Dopamine dependency of oscillations between subthalamic nucleus and pallidum in Parkinson's disease. J Neurosci
21, 10331038.
Buzsáki
G
&
Draguhn
A (2004). Neuronal oscillations in cortical networks. Science
304, 19261929.
Canteras NS, Shammah-Lagnado SJ, Silva BA & Ricardo JA (1990). Afferent connections of the subthalamic nucleus: a combined retrograde and anterograde horseradish peroxidase study in the rat. Brain Res 513, 4359.[CrossRef][Medline]
Cassidy
M, Mazzone
P, Oliviero
A, Insola
A, Tonali
P, Di Lazarro
V
&
Brown
P (2002). Movement-related changes in synchronization in the human basal ganglia. Brain
125, 12351246.
Courtemanche
R, Fujii
N
&
Graybiel
AM (2003). Synchronous, focally modulated beta-band oscillations characterize local field potential activity in the striatum of awake behaving monkeys. J Neurosci
23, 1174111752.
Donoghue JP & Parham C (1983). Afferent connections of the lateral agranular field of the rat motor cortex. J Comp Neurol 217, 390404.[CrossRef][Medline]
Donoghue
JP, Sanes
JN, Hastopoulos
NG
&
Gaál
G (1998). Neural discharge and local field potential oscillations in primate motor cortex during voluntary movements. J Neurophysiol
79, 159173.
Donoghue JP & Wise SP (1982). The motor cortex of the rat: cytoarchitecture and microstimulation mapping. J Comp Neurol 212, 7688.