J Physiol Wellcome Trust-funded researchers
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Physiol Volume 561, Number 1, 205-214, November 15, 2004 DOI: 10.1113/jphysiol.2004.075325
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
561/1/205    most recent
jphysiol.2004.075325v1
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Staba, R. J
Right arrow Articles by Barth, D. S
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Staba, R. J
Right arrow Articles by Barth, D. S

Dissociation of slow waves and fast oscillations above 200 Hz during GABA application in rat somatosensory cortex

Richard J Staba1, Peter C Bergmann1 and Daniel S Barth1

1 Department of Psychology, University of Colorado, Boulder, CO 80309-0345, USA


    Abstract
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Fast electrical oscillations (FOs; > 200 Hz), superimposed on vibrissa-evoked slow potentials, may support rapid sensory integration in neocortex. Yet, while it is well established that the positive/negative (P1/N1) slow wave of the somatosensory evoked potential primarily reflects sequential activation of supragranular and infragranular pyramidal cells mediated chiefly via excitatory chemical synaptic pathways, little is known about the generation of FOs. In this study, laminar current source–density analysis and principal component analysis indicated that FOs are generated by two dipolar current sources situated in the supra- and infragranular layers, similar in laminar location to the two current dipoles associated with the slow wave. However, exogenous GABA application reversibly abolished the N1 slow wave, leaving the P1 intact, while the FO was unaffected by GABA. Furthermore, reductions in both supra- and infragranular cortical unit discharge during application of GABA suggests that FO generation is not dependent on the same intracortical synaptic circuits that are associated with the N1 slow wave. These data suggest a marked functional dissociation between neural mechanisms underlying the slow and fast components of the vibrissa-evoked response.

(Received 7 September 2004; accepted after revision 1 October 2004; first published online 1 October 2004)
Corresponding author R. J. Staba: Department of Psychology, University of Colorado, Boulder, CO 80309-0345, USA. Email: staba{at}psych.colorado.edu


    Introduction
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Transient displacement of a rodent's facial vibrissa evokes a cascade of both slow and fast activity associated with the somatosensory evoked potential (SEP) that can be recorded within the posteromedial barrel subfield (PMBSF) of the rat. Superimposed on the initial slow wave components of the SEP (labelled P1/N1 to denote their polarity and sequence of occurrence) are fast electrical oscillations (FOs; > 200 Hz; Jones & Barth, 1999; Jones et al. 2000; Barth, 2003) that resemble fast electrical activity found in human somatosensory cortex during median nerve stimulation (Curio et al. 1994; Klostermann et al. 1998). Results from recent studies (Jones & Barth, 1999; Barth, 2003) suggest that FOs may uniquely reflect very rapid cellular interactions involved with somatosensory information processing.

While the somatosensory P1/N1 is thought to be produced by sequential activation of supra- and infragranular pyramidal cells via excitatory synaptic connections between cell groups (Simons, 1978; Simons & Woolsey, 1979; Di et al. 1990; Agmon & Connors, 1991), the neural mechanisms of FOs are unclear. Unit recordings indicate that FOs are associated with rapid and time-locked population discharge of cortical pyramidal cells (Kandel & Buzsaki, 1997; Jones et al. 2000; Barth, 2003). While pyramidal cell action potentials (APs) cannot follow FO frequencies, FOs appear to define preferred latencies for firing of individual cells within a population (Jones et al. 2000; Barth, 2003). The cellular mechanism responsible for synchronizing this population response remains a major question. One possibility is that FOs reflect synchronized inhibitory drive imposed by GABAergic interneurones, similar to the generator of high frequency ripples in hippocampus (Buzsaki et al. 1992). This possibility is consistent with intracellular recordings indicating that fast spiking (FS) cells, thought to represent inhibitory interneurones, present phase-locked APs at FO frequencies in response to physiological stimulation (Jones et al. 2000). However, the role of inhibitory interneurones has been called into question by evidence that the GABAA receptor antagonist, bicuculline methiodide, fails to attenuate FOs and instead increases their duration along with development of epileptiform slow waves (Jones & Barth, 2002). This raises the distinct possibility that, instead of inhibition, FOs may be synchronized by excitatory interactions between pyramidal cells.

In this light, FOs may rely on excitatory interactions within the pyramidal population mediated by chemical synaptic transmission similar to that responsible for the P1/N1 slow wave. Previous studies have shown that epicortical GABA application effectively abolishes the N1 slow wave, while leaving the P1 intact, suggesting that GABA differentially suppresses intracortical synaptic excitation of infragranular pyramidal cells responsible for the N1 (Brailowsky & Knight, 1984; Toldi & Feher, 1985; Kaneko & Hicks, 1988). In the present study, we tested the hypothesis that FOs rely on an excitatory pathway similar to that responsible for the slow wave complex by recording laminar field potentials and single unit responses before and during GABA application in the rat PMBSF.


    Methods
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Animals and surgery

All procedures were conducted within the guidelines established by the University of Colorado Institutional Animal Care and Use Committee. Adult male Sprague-Dawley rats (n = 28, 225–300 g) were anaesthetized using a mixture of ketamine HCl (64 mg kg–1), xylazine (13 mg kg–1) and acepromazine (2.1 mg kg–1), and placed on a regulated heating pad to maintain normal body temperature (37°C). A unilateral craniotomy was performed over the right hemisphere extending from bregma to lambda and from the midsagittal suture to the lateral aspect of the temporal bone, exposing a large area of parietotemporal cortex. The dura was reflected and the exposed cortex regularly wetted with Ringer solution containing (mM): NaCl 135, KCl 3, MgCl2 2 and CaCl2 2; pH 7.3 at 37°C. Anaesthesia levels were maintained throughout the experiment so that the eye blink and flexor withdrawal reflexes could barely be elicited. At the conclusion of the experiment, animals were killed with an overdose of anaesthesia without regaining consciousness.

Stimulation and recording

Vibrissa on the left mystacial pad were displaced ~300 µm in the ventrodorsal direction (0.3 ms duration at 1 Hz) using a laboratory-built solenoid (Jones & Barth, 1999; Barth, 2003). Field potential recordings were made using a laminar electrode consisting of a linear array of 24 platinum contacts (diameter, 40 µm; intercontact spacing, 100 µm) extending 2.3 mm beneath the cortical surface (Neurotrack KFT, Budapest, Hungary). The electrode was centrally positioned within the cortical barrel field and inserted perpendicular to the surface nearest areas generating maximum amplitude SEPs as determined from surface recordings using a 64-contact electrode array arranged in an 8 x 8 grid. For single unit recordings, a stabilizing pressure plate with central access hole for a microelectrode was positioned within the vibrissa/barrel field using the same techniques as used for laminar recordings. The pressure plate also had a single silver wire electrode mounted to the side of the access hole to monitor epicortical SEPs concurrently with unit responses. Fine tipped (1.5 µm; 20–30 M{Omega}) glass microelectrodes filled with Ringer solution were inserted through the pressure plate and advanced through the cortex in 2 µm increments until single units responding to stimulation of the vibrissa were isolated. Single unit isolation was verified visually by the presence of only one AP firing pattern with consistent amplitude and morphology. Laminar field potentials and unit recordings were referenced to a AgCl-coated electrode placed within the dorsal neck muscle. Field potentials were amplified (x 200), filtered (1–3000 Hz) and digitized at 10 kHz. Unit recordings were amplified (x 1000), filtered (300–3000 Hz) and also digitized at 10 kHz.

Data collection

For all experiments, single trial data (n = 100) was recorded during three conditions: before GABA application; during GABA exposure; and 30 min after washout. GABA (50 mM; Sigma, St Louis, MO, USA) was dissolved in Ringer solution buffered with 5 mM Hepes that had a final pH of 7.3. GABA was applied to the cortical barrel field using a 4-mm2 pledget that was notched to surround the laminar electrode (n = 10). The GABA-soaked pledget remained in place for ~5 min during the recording of evoked responses. To verify that GABA was penetrating the infragranular lamina, in five animals cortical depth injections were made using a microinjection system (WPI Inc., Sarasota, FL, USA). From the cortical surface to a depth of ~1 mm, and < 1 mm from the laminar electrode, GABA was injected at a rate of 23 nl s–1 for 5–10 s. In a separate series of animals used for single unit recordings (n = 13), the access hole in the pressure plate was filled with GABA with time allowed for GABA to diffuse throughout the cortical barrel field. Recovery trials were recorded 30 min after GABA washout using Ringer solution.

Laminar analysis

For each animal, laminar field potentials were averaged across trials and digitally band-pass filtered to separately examine the slow wave (1–100 Hz) and FO (250–600 Hz). To better resolve the locations of extracellular current sinks and sources along the axis of the electrode array, spatially filtered current source density (CSD) (low pass cut-off, 1 Hz mm–1; Nicholson & Freeman, 1975; Rappelsberger et al. 1981; Mitzdorf, 1985) was calculated using the equation described by Mitzdorf (1985). This equation computes CSD as the second spatial derivative of potential measured in 100 µm increments along an axis perpendicular to the cortical surface. Preliminary inspection of the 24 CSD waveforms recorded from the laminar electrode revealed that the eight most distal contacts located within superficial white matter did not contribute additional information to the CSD profile. Subsequent analysis and illustration used the remaining 16 contacts spanning the cortical laminae.

Principal component analysis (PCA) was used to separate the independent sources of covariance in the CSD across the laminar electrode (Ruchkin et al. 1964; Barth et al. 1989). In the context of this study, the principal components represented distinct neuronal populations generating unique spatial patterns of CSD, covarying across multiple contacts of the laminar electrode. Component loadings reflected the weighted contribution of each neuronal population to the CSD at each electrode contact. Preliminary results indicated that the laminar patterns of component loadings were approximately dipolar. To simplify the model of activity underlying the evoked response, principal components were rigidly rotated so that their laminar loading patterns best reflected dipolar sources. Component loadings were derived from the baseline CSD only. Multiple regression was then used to calculate a linear combination of these component loadings that best fitted the time course of the evoked response at baseline and during application of GABA. Reconstructions of the spatiotemporal patterns of activity contributed by each principal component were computed by multiplying the spatial loading pattern (xn) by their respective temporal series of regression weights (Reg Wgts) (tn). A composite model of activity was computed by summing the spatiotemporal patterns of activity from the first two principal components, such that:


{tjp_570_m1}

The P1/N1 and FO amplitude were computed from the regression weights as the maximum/minimum peak and root-mean-square values, respectively. FO centre frequency was defined as the largest peak appearing in a 512-point fast Fourier transform of the bandpass filtered data. Comparison between the composite model and CSD profile was calculated as follows:


{tjp_570_m2}

Single units were initially segregated into groups according to the phase relationship between their vibrissa-evoked firing pattern and concurrently recorded surface FOs. In order to eliminate slight intertrial variations in FO latency that could obscure phase relationships between FOs and unit discharges, FOs and APs were phase-aligned by a ± 1.5 ms shift determined from the cross-correlation function computed between single trial FOs and grand averaged FOs for a given unit. Distinct temporal patterns of a single unit's discharge were additionally segregated using a K-means cluster algorithm (Hartigan & Wong, 1979). The grouping variable was the poststimulus latencies of unit responses. In this application, if all trials consisted of unit responses at the same poststimulus latencies, then a single cluster containing all trials would yield the minimum variance around the mean response latency. However, if unit response latencies differed systematically in subgroups of trials (i.e. occurring at the latency of the second FO wave in some trials, first and fourth FO waves on other trials, etc.) then trials clustered within each subgroup would yield a smaller within-group variance than if all trials were treated as a single group. The maximum number of permitted clusters was six. Changes in unit discharge during application of GABA were then quantified by computing the ratio of mean APs during baseline to those during the application of GABA for each isolated unit. Analysis of variance (ANOVA) and paired t tests were used to compare evoked responses across conditions and changes in AP discharge, respectively. Significance for all comparisons was set at P = 0.05.


    Results
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Rapid, transient displacement of the contralateral vibrissae typically evoked a stereotyped SEP during the pre-GABA period that consisted of a large amplitude (mean peak-to-peak (p–p) amplitude, 0.93 ± 0.23 mV (S.D.)) biphasic slow wave (Fig. 1A), and small amplitude (56 ± 34 µV) FO superimposed on the P1 and rising phase of the N1 (Fig. 1A; FO indicated by asterisks; Fig. 1B). CSD profiles of pre-GABA P1/N1 complex (Fig. 1C; CSD) indicated that the P1 was associated with a source within superficial layers that reached peak amplitude 15 ± 2.1 ms following stimulus onset and complementary sink distributed across the middle/deep cortical layers. The N1 was associated with a sink near the surface that extended across the middle cortical layers, and a large complementary source within deeper cortical layers that reached maximum amplitude 22 ± 2.1 ms following stimulus onset.



View larger version (35K):
[in this window]
[in a new window]
 
Figure 1.  Slow and fast components of the somatosensory evoked potential (SEP)
A and B, SEP consists of slow wave (band pass filtered, 1–2000 Hz) and fast oscillation (FO; 250–600 Hz). P1 and N1 denote polarity of slow wave complex. Asterisks indicate superimposition of FO on P1 and rising phase of N1 components. Scale bar, 0.5 mV for slow wave, 50 µV for FO. C, averaged CSD profile (n = 15) of slow wave complex. Spatiotemporal pattern of sources (shaded grey) and sinks (shaded black) contributing to the P1/N1. Continuous lines denote maximum peak amplitude of supra- and infragranular sources contributing to P1 and N1, respectively. Reconstruction of supragranular (Comp 1) and infragranular (Comp 2) patterns of activity contributing to the P1/N1 complex computed by multiplying the spatial loading pattern (Load) by their respective temporal pattern of regression weights (Reg Wgts). Dashed lines surrounding Loads and Reg Wgts denote 95% confidence intervals calculated across all animals. Scale on the right represents approximate depth (in mm) and cortical lamina of recording electrode within the vibrissal/barrel field of the somatosensory cortex. D, same as C, but of vibrissa-evoked FO CSD profile recorded from a single animal. Superimposition of Reg Wgts reveal the variability in the initiation of FO following stimulus onset across all animals.

 
PCA accounted for 97 ± 2.0% of system variance in the CSD profile with only two components. Reconstructions of activity contributed by each component were computed by multiplying the component-loading pattern (Fig. 1C; ‘Load’) by their respective regression weight (Fig. 1C; ‘Reg Wgts’). The loading patterns reflect the relative amplitude and polarity of the components at each electrode contact, and the regression weights reflect the weighted contribution of the components over the time course of the evoked response. Reconstructions of the first two principal components produced patterns reflecting the asynchronous activity from two spatially distinct cellular populations. The first reconstruction depicted a vertically orientated current dipole reversing polarity ~300 µm beneath the cortical surface that attained maximum amplitude 15 ms following stimulus onset (Fig. 1C; Comp 1). Similarly, the second reconstruction revealed a current dipole reversing at ~1100 µm with a poststimulus latency of 22 ms (Fig. 1C; Comp 2).

Figure 1D displays similar laminar analysis of the band pass filtered (250–600 Hz) FO (mean, 338 ± 60 Hz) that coincided with the SEP slow wave. Inspection of the CSD profile revealed a complex pattern of spatiotemporal activity that appeared as a series of waves propagating between the cortical depth and surface with no clear reversal in polarity in the depth. Unlike the slow wave response, the poststimulus latency of FOs was variable across animals and precluded the calculation of a grand average CSD profile. The example of Fig. 1D is therefore from a single animal. Similar to the analysis of the slow wave, results from PCA showed that FOs could be modelled by the activation of two spatially distinct cellular populations. Reconstructions of the first two principal components again suggested contributions from supra- and infragranular populations. Dipolar loading patterns of the two components were consistent across animals and reversed polarity ~500 and 1000 µm below the cortical surface, respectively. Superimposition of regression weights show the variability in FO onset across animals. Yet, models of FOs, based on the patterns of activity from the first two principal components of the individual animals, accounted for 97 ± 1.6% of the variance in the CSD.

Application of exogenous GABA simplified the laminar profile of the vibrissa-evoked slow wave response (Fig. 2A), but had little effect on the generation of FOs (341 ± 58 Hz; Fig. 2B). GABA enhanced the P1 and largely suppressed the N1. Suppression of the N1 was consistently observed in all animals and occurred within 1 min of GABA exposure. During GABA application the evoked response began with a P1-like slow wave peaking at a slightly longer latency (18 ± 1.6 ms) and longer duration compared to the P1 before GABA application. Inspection of the CSD profile revealed that the evoked response was produced by a current sink distributed across the middle/deep cortical layers that was surrounded by a large source within superficial layers and a smaller source located within the deepest cortical layers (Fig. 2A). This pattern resembled the profile of activity that occurred during the time course of the P1 before GABA application (cf. Fig. 1C). However, the ensuing large polarity reversal of the sink-source pair within the depth, associated with the N1, was greatly diminished during GABA application.



View larger version (57K):
[in this window]
[in a new window]
 
Figure 2.  Evoked slow wave and FO during GABA application
A, during GABA application, the slow wave CSD profile reduces to a large superficial source and complementary deep sink that is accompanied by a smaller source within the deepest layers. Reconstructions show the evoked response results primarily from cellular activity within supragranular layers (Comp 1) with greatly reduced contributions from infragranular cells (Comp 2). B, in contrast to the slow wave, the FO CSD profile and cellular activity within supra- and infragranular layers (Comp 1 and 2, respectively) contributing to FOs were largely unchanged during GABA application compared to baseline.

 
Results of PCA showed that the GABA response resulted primarily from the activity of the supragranular cell population, indicated by the large current dipole reversing polarity in the supragranular layers (Fig. 2A; Comp 1). Activation of the infragranular cell population was reduced, reflected by greatly attenuated regression weights for the second component (Fig. 2A; Comp 2). Analysis of P1 and N1 amplitude revealed that the P1 was significantly larger during GABA application compared to before GABA (142 ± 39 versus 82 ± 36 µV mm–2, P < 0.0001), while the N1 amplitude was significantly smaller during GABA application (24 ± 10 versus 57 ± 19 µV mm–2, P < 0.0001). The linear combination of the first two principal components accounted for 96 ± 2.7% of the variance in the CSD profile.

In contrast to the pronounced effects that GABA had on the P1/N1 complex, no differences were observed in the CSD of FOs during GABA exposure (Fig. 2B). Using the component loadings derived from the baseline pre-GABA period, reconstructions of the first two principal components revealed that contributions of the supra- and infragranular current dipoles to the FOs were virtually unchanged by GABA. Analysis comparing amplitude before and during GABA application found no difference in FO amplitude generated from either supragranular (RMS amplitude, 0.77 ± 0.37 versus 0.74 ± 0.37 µV mm–2, P = 0.9) or infragranular populations (0.88 ± 0.61 versus 0.86 ± 0.53 µV mm–2, P = 0.9). The regression model, based on pre-GABA component loadings, accounted for 96 ± 3.2% of the system variance during exposure to GABA.

Based on distinct firing patterns and the presence or absence of phase-locking between APs and FOs, the 33 units studied were segregated into four groups (Fig. 3AD). The first group consisted of two cells that consistently responded with bursts of three to five APs at the same frequency and phase as concurrently recorded FOs (Fig. 3A). The phase relationship between temporal patterns of AP bursts and FOs varied across clusters of single trials (Fig. 3A; {kappa}1–{kappa}6) but each cluster contained multiple APs with interspike intervals matching the FOs. These two cells were also distinct in that GABA did not decrease their firing, and in this example, actually increased firing over baseline by 160%. It should be noted that these units also displayed APs with half amplitude widths of less than 0.5 ms, but this was not considered a reliable distinguishing feature as eight other cells showed the same characteristic, but without high frequency firing. The second group consisted of cells (n = 12) that typically fired one to two APs with an interspike interval greater than, but at a multiple of, the FO period. Single trials clustered into subgroups that showed AP discharge at preferred latencies of the FOs (Fig. 3B). Three cells comprised a third group that typically fired a single AP at only one phase of the FOs (Fig. 3C). Finally, approximately half the cells recorded (n = 16) fitted into a fourth group that typically fired two or more APs at interspike intervals exceeding the FO period and showing no locking to a consistent phase even after clustering of single trials (Fig. 3D). All units, except for the first group, exhibited a significant decrease in firing during GABA application compared to baseline (60 ± 6.5%, P = 0.00007, n = 31). Similar attenuation of unit firing during GABA exposure was observed when cells were separately analysed in the supragranular layers (60 ± 9.6%, P = 0.001, n = 10) and infragranular layers (60 ± 8.4%, P = 0.00005, n = 21), and differences between supra- and infragranular attenuation were not significant. Excluding the first group, there were no significant differences in GABA-dependent attenuation between cell groups.



View larger version (30K):
[in this window]
[in a new window]
 
Figure 3.  Raster plots of unit discharge before and during GABA application
A–D, examples of cortical single units grouped according to their distinct discharge patterns (n = number of cells within each group). Individual trials of evoked responses were segregated into clusters ({kappa}1–{kappa}6) based on their phase-locked AP discharge to concurrently recorded FOs (indicated by shaded bars). Note that most units demonstrated some phase-locked discharge associated with FOs and all units, with the exception of two in the first group (A), decreased their firing during exposure to GABA. E, patterns of thalamic MUA displayed no apparent phase-locked discharge with cortical FOs before or during GABA application.

 
The failure of GABA to attenuate FOs, despite its large effect on the N1 slow wave and modest effect on unit firing rates, prompted two additional experiments to further clarify this phenomenon. The first experiment examined whether phase-locked, high frequency thalamocortical input contributed to persistent FOs even after intracortical circuits were suppressed by GABA. Eighteen units in six additional animals were recorded from microelectrodes positioned sterotaxically within the ventrobasal nucleus (VB) of the thalamus (A–P, 3.6 mm; M–L, 2.7 mm; D–V, 6.0 mm), and concurrent FOs were recorded from a single macroelectrode centred in the PMBSF. Thalamic microelectrodes typically displayed lower impedance (< 2 M{Omega} at 1 kHz) than cortical microelectrodes and therefore were assumed to reflect multiunit activity (MUA) as opposed to single units. Firing patterns of thalamic units displayed no apparent relationship to cortical FOs either before or after GABA application (Fig. 3E). MUA typically consisted of one to three rapid discharges at 5–7 ms poststimulus, before cortical FOs appeared. This was often followed by activity during and after cortical FOs emerged but revealed no phase-locking with the cortical signal. The temporal pattern of MUA was not changed in any of the thalamic units by the epicortical application of GABA.

The second experiment examined whether gap junctions might be implicated in FO generation. Epipial slow waves and FOs were recorded from the PMBSF before and during the systemic administration of halothane, which is an anaesthetic known to close gap junctions (Terrar & Victory, 1988; Burt & Spray, 1989; Peinado et al. 1993). Halothane (1.5%) was administered to four additional animals for a period of 4–6 min following baseline recording. During this initial period, there was little effect on the amplitude of the slow wave (Fig. 4A; dark traces) compared to baseline recording (Fig. 4A; light traces). The RMS amplitude of the slow wave during halothane exposure was 91% that of baseline (P = 0.37). In contrast, concurrently recorded FO were attenuated by 70% (P = 0.0009) during halothane exposure (Fig. 4B; dark traces) compared to baseline (Fig. 4B; light traces). The differential effect of halothane on FO compared to the slow wave could be maintained for up to 1–3 min.



View larger version (25K):
[in this window]
[in a new window]
 
Figure 4.  Slow wave and FO during halothane
A, examples of surface recording within the PMBSF of the slow wave before (grey traces) and during halothane exposure (black traces). Overall, slow wave amplitude was nearly identical during halothane exposure compared to baseline. B, in contrast to the slow wave, FO amplitude was attenuated by 70% during halothane exposure compared to baseline. Example from a single animal.

 

    Discussion
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
The present data demonstrate a marked dissociation between the effects of GABA on slow wave and FO components of the SEP and suggest that they rely on distinct neural mechanisms. Thalamic unit recording indicates that the persistence of cortical FOs in the presence of GABA cannot be explained by high frequency thalamocortical drive as the response latency and firing pattern of VB cells is unrelated to concurrently recorded FOs. The selective suppression of FOs with halothane suggests that, similar to ripples in hippocampus, cortical FOs may rely on gap junctions for their generation and/or synchronization.

Laminar analysis of the slow wave closely resembles previous studies (Mitzdorf, 1987; Vaknin, 1989; Barth & Di, 1990; Di et al. 1990; Sukov & Barth, 1998), indicating thalamocortical activation of supra- and infragranular pyramidal cells producing the P1, followed by longer latency intracortical activation of infragranular pyramidal cells producing the N1. As previously shown in cat somatosensory cortex (Brailowsky & Knight, 1984), GABA suppresses the N1 with a reduced contribution from infragranular pyramidal cells. The simplest explanation for this effect is that GABA widely hyperpolarizes cortical neurones, resulting in an average reduction of firing to 60% of baseline. The resiliency of the P1 slow wave in the presence of GABA suggests that thalamocortical input is sufficient to evoke postsynaptic activity within supragranular cells, but that many of these cells fail to fire and thus fail to establish postsynaptic potentials in infragranular pyramidal cells. This possibility is consistent with evidence that suggests that thalamocortical transmission may be more reliable and efficient compared to intracortical synaptic transmission (Kharazia & Weinberg, 1994; Gil et al. 1999; see also Castro-Alamancos & Connors, 1997; Amitai, 2001).

Laminar analysis of FOs suggests that they too are produced by asynchronous activation of supra- and infragranular cells. This conclusion is supported by simultaneous field potential and multiunit recordings revealing neuronal discharge associated with individual cycles of FOs (Kandel & Buzsaki, 1997; Jones et al. 2000; Barth, 2003) and suggesting that cortical circuit interactions may contribute to FO generation. Data on higher frequency sigma-bursts observed in humans (Gobbele et al. 1998; Gobbele et al. 1999) and piglets (Ikeda et al. 2002) have suggested that an additional extra-cortical oscillatory component may also reside within the thalamus. However, rhythmic thalamic input does not appear to play a role in rat somatosensory cortex as FOs may be evoked during direct cortical stimulation with or without an intact thalamic VB (Staba et al. 2003), and are abolished during cortical cooling (Staba et al. 2003) or cortical application of kynurenic acid (Ikeda et al. 2002), despite intact, non-rhythmic thalamic input. This conclusion is also supported by the present results, indicating no rhythmic MUA in the VB corresponding to concurrently recorded cortical FOs.

While the laminar profile of FOs appears similar to that of the P1/N1 slow wave complex, FOs are unchanged by GABA, indicating that a distinct neural mechanism underlies their generation. It has been proposed that FOs may reflect excitatory interactions between cortical pyramidal cells (Jones et al. 2000; Jones & Barth, 2002; Barth, 2003) resulting in tightly synchronized population spikes. Recent evidence indicates that regular-spiking cells, presumably pyramidal cells (Connors & Gutnick, 1990), can discharge at preferred latencies that correspond with individual cycles of the oscillation (Jones et al. 2000; Barth, 2003). The responses of approximately half of the units studied here during baseline conditions were quite similar to those obtained with previous intracellular recording from regular-spiking cells. While no single pyramidal cell can fire at FO frequencies, coordinated AP firing within a population of pyramidal cells may generate coherent, small amplitude extracellular currents that could produce patterns in the laminar CSD similar to the patterns observed in the present study.

If FOs reflect synchronized population spikes in cortical pyramidal cells, the question remains how they persist in the presence of GABA, which substantially decreases cell firing and therefore intracortical transmission upon which the postsynaptic N1 slow wave relies. Simultaneous displacement of multiple vibrissae generates a synchronous thalamocortical volley, apparently of sufficient strength to overcome cellular hyperpolarization in a limited number of supra- and infragranular pyramidal cells, resulting in FOs without the accompanying large amplitude N1 slow wave. Yet, if FOs require circuit interactions between pyramidal cells that are based on excitatory synaptic transmission, it is unlikely that they would remain unchanged by GABA.

The present observation that the gap junction blocker, halothane, produces a marked attenuation of FOs compared to the slow wave, suggests the possibility that excitatory coupling giving rise to FOs is not based strongly on chemical synaptic transmission, and may occur instead through electrical coupling. Precedent for this suggestion does exist in studies of hippocampal ripple (~200 Hz) oscillations. Modelling studies predict that high frequency activity, in the absence of chemical synaptic transmission, may occur through a network of hippocampal pyramidal cells electrically coupled by their axons (Traub et al. 1999). Simulations of ripple activity can be produced with a sparse connectivity between axons and low frequency of ectopic spike generation (Traub et al. 1999). Indeed, evidence from in vitro hippocampal studies shows that only a small proportion of principal neurones participate in each ripple, and ripples can be extinguished with halothane, indicating that they may emerge through a network of electrically coupled principal cells (Ylinen et al. 1995; Draguhn et al. 1998; Schmitz et al. 2001). If a similar mechanism is responsible for FOs in cortex, it could explain why FOs persist during partial suppression of pyramidal cell firing sufficient to eliminate the N1, yet, unlike the N1, are nearly extinguished by halothane.

However, analogies between somatosensory cortex and hippocampus must be made tentatively. While evidence exists for cortical gap junctions (Gutnick & Prince, 1981), more evidence is needed concerning electrical coupling between axons of cortical pyramidal cells. In contrast to pyramidal cells, there is strong evidence for electrical coupling between fast-spiking inhibitory cells in neocortex (Galarreta & Hestrin, 1999), which could serve to synchronize their discharge at FO frequencies. We identified only two units with high frequency firing patterns and phase-locking to the FOs that were similar to intracortical recordings of inhibitory interneurones responding to vibrissa stimulation (Jones et al. 2000). Inhibitory interneurones do receive direct excitatory thalamocortical input that would be expected to persist in the absence of intracortical excitatory drive, and when synchronized in sufficient numbers, could produce small dipolar field potentials in dendrites that radiate in a predominately vertical orientation (Jones & Hendry, 1984). Yet, while the present study indicates that the mechanism of FO generation is clearly different than that of the slow wave-evoked potential and may rely on gap junctions, the effect of GABA and halothane on both principal cells and interneurones must now be examined and characterized with intracellular recording techniques to determine their specific contribution to FOs in somatosensory cortex.


    References
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Agmon A & Connors BW (1991). Thalamocortical responses of mouse somatosensory (barrel) cortex in vitro. Neuroscience 41, 365–379.[CrossRef][Medline]

Amitai Y (2001). Thalamocortical synaptic connections: efficacy, modulation, inhibition and plasticity. Rev Neurosci 12, 159–173.[Medline]

Barth DS (2003). Submillisecond synchronization of fast electrical oscillations in neocortex. J Neurosci 23, 2502–2510.[Abstract/Free Full Text]

Barth DS & Di S (1990). Three-dimensional analysis of auditory-evoked potentials in rat neocortex. J Neurophysiol 64, 1527–1536.[Abstract/Free Full Text]

Barth DS, Di S & Baumgartner C (1989). Laminar cortical interactions during epileptic spikes studied with principal component analysis and physiological modeling. Brain Res 484, 13–35.[CrossRef][Medline]

Brailowsky S & Knight RT (1984). Inhibitory modulation of cat somatosensory cortex: a pharmacological study. Brain Res 322, 310–315.[CrossRef][Medline]

Burt JM & Spray DC (1989). Volatile anesthetics block intercellular communication between neonatal rat myocardial cells. Circ Res 65, 829–837.[Abstract/Free Full Text]

Buzsaki G, Horvath Z, Urioste R, Hetke J & Wise K (1992). High-frequency network oscillation in the hippocampus. Science 256, 1025–1027.[Abstract/Free Full Text]

Castro-Alamancos MA & Connors BW (1997). Thalamocortical synapses. Prog Neurobiol 51, 581–606.[CrossRef][Medline]

Connors BW & Gutnick MJ (1990). Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci 13, 99–104.[CrossRef][Medline]

Curio G, Mackert BM, Burghoff M, Koetitz R, Abraham-Fuchs K & Harer W (1994). Localization of evoked neuromagnetic 600 Hz activity in the cerebral somatosensory system. Electroencephalogr Clin Neurophysiol 91, 483–487.[CrossRef][Medline]

Di S, Baumgartner C & Barth DS (1990). Laminar analysis of extracellular field potentials in rat vibrissa/barrel cortex. J Neurophysiol 63, 832–840.[Abstract/Free Full Text]

Draguhn A, Traub RD, Schmitz D & Jefferys JGR (1998). Electrical coupling underlies high-frequency oscillations in the hippocampus in vitro. Nature 394, 189–192.[CrossRef][Medline]

Galarreta M & Hestrin S (1999). A network of fast-spiking cells in the neocortex connected by electrical synapses. Nature 402, 72–75.[CrossRef][Medline]

Gil Z, Connors BW & Amitai Y (1999). Efficacy of thalamocortical and intracortical synaptic connections: quanta, innervation, and reliability. Neuron 23, 385–397.[CrossRef][Medline]

Gobbele R, Bucher H & Curio G (1998). High-frequency (600Hz) SEP activities originating in the subcortical and cortical human somatosensory system. Electroencephalogr Clin Neurophysiol 108, 182–189.[CrossRef][Medline]

Gobbele R, Buchner H, Scherg M & Curio G (1999). Stability of high-frequency (600Hz) components in human somatosensory evoked potentials under variation of stimulus rate – evidence for a thalamic origin. Clin Neurophysiol 110, 1659–1663.[CrossRef][Medline]

Gutnick MJ & Prince DA (1981). Dye coupling and possible electrotonic coupling in the guinea pig neocortical slice. Science 211, 67–70.[Abstract/Free Full Text]

Hartigan JA & Wong MA (1979). Algorithm AS 136: a K-means clustering algorithm. Appl Stat 28, 100–108.[CrossRef]

Ikeda H, Leyba L, Bartolo A, Wang Y & Okada C (2002). Synchronized spikes of thalamocortical axonal terminals and cortical neurons are detectable outside the pig brain with MEG. J Neurophysiol 87, 626–630.[Abstract/Free Full Text]

Jones MS & Barth DS (1999). Spatiotemporal organization of fast (>200 Hz) electrical oscillations in rat vibrissa/barrel cortex. J Neurophysiol 82, 1599–1609.[Abstract/Free Full Text]

Jones MS & Barth DS (2002). Effects of bicuculline methiodide on fast (>200Hz) electrical oscillations in rat somatosensory cortex. J Neurophysiol 88, 1016–1025.[Abstract/Free Full Text]

Jones EG & Hendry SHC (1984). Basket Cells. In Cellular Components of the Cerebral Cortex, vol. 1, ed. Peters A & Jones EG,pp. 309–336. Plenum Press, New York.

Jones MS, MacDonald KD, Choi B, Dudek FE & Barth DS (2000). Intracellular correlates of fast (>200 Hz) electrical oscillations in rat somatosensory cortex. J Neurophysiol 84, 1505–1518.[Abstract/Free Full Text]

Kandel A & Buzsaki G (1997). Cellular-synaptic generation of sleep spindles, spike-and-wave discharges, and evoked thalamocortical responses in the neocortex of the rat. J Neurosci 17, 6783–6797.[Abstract/Free Full Text]

Kaneko T & Hicks TP (1988). Baclofen and gamma-aminobutyric acid differentially suppress the cutaneous responsiveness of primary somatosensory cortical neurones. Brain Res 443, 360–366.[CrossRef][Medline]

Kharazia VN & Weinberg RJ (1994). Glutamate in thalamic fibers terminating in layer IV of primary sensory cortex. J Neurosci 14, 6021–6032.[Abstract]

Klostermann F, Nolte G, Losch F & Curio G (1998). Differential recruitment of high frequency wavelets (600Hz) and primary cortical response (N20) in human median nerve somatosensory evoked potentials. Neurosci Lett 256, 101–104.[CrossRef][Medline]

Mitzdorf U (1985). Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena. Physiol Rev 65, 37–100.[Free Full Text]

Mitzdorf U (1987). Properties of the evoked potential generators: current source-density analysis of visually evoked potentials in the cat cortex. Int J Neurosci 33, 33–59.[Medline]

Nicholson C & Freeman JA (1975). Theory of current source-density analysis and determination of conductivity tensor for anuran cerebellum. J Neurophysiol 38, 356–368.[Abstract/Free Full Text]

Peinado A, Yuste R & Katz LC (1993). Extensive dye coupling between rat neocortical neurons during the period of circuit formation. Neuron 10, 103–114.[CrossRef][Medline]

Rappelsberger P, Pockberger H & Petsche H (1981). Current source density analysis: methods and application to simultaneously recorded field potentials of the rabbit's visual cortex. Pflugers Arch 389, 159–170.[CrossRef][Medline]

Ruchkin DS, Villegas J & John ER (1964). An analysis of average evoked potentials making use of least mean square techniques. Ann N Y Acad Sci 115, 799–821.[Medline]

Schmitz D, Schuchmann S, Fisahn A, Draguhn A, Buhl EH, Petrasch-Parwez E et al. (2001). Axo-axonal coupling: a novel mechanism for ultrafast neuronal communication. Neuron 31, 831–840.[CrossRef][Medline]

Simons DJ (1978). Response properties of vibrissa units in rat SI somatosensory neocortex. J Neurophysiol 41, 798–820.[Abstract/Free Full Text]

Simons DJ & Woolsey TA (1979). Functional organization in mouse barrel cortex. Brain Res 165, 327–332.[CrossRef][Medline]

Staba RJ, Brett-Green B, Paulsen M & Barth DS (2003). Effects of ventrobasal lesion and cortical cooling on fast oscillations (>200Hz) in rat somatosensory cortex. J Neurophysiol 89, 2380–2388.[Abstract/Free Full Text]

Sukov W & Barth DS (1998). Three-dimensional analysis of spontaneous and thalamically evoked gamma oscillations in auditory cortex. J Neurophysiol 79, 2875–2884.[Abstract/Free Full Text]

Terrar DA & Victory JG (1988). Influence of halothane on electrical coupling in cell pairs isolated from guinea-pig ventricle. Br J Pharmacol 94, 509–514.[Medline]

Toldi J & Feher O (1985). The role of the gamma-aminobutyrate system in the interactions of cortical evoked potentials. Neuroscience 15, 977–981.[CrossRef][Medline]

Traub RD, Schmitz D, Jefferys JGR & Draguhn A (1999). High-frequency population oscillations are predicted to occur in hippocampal pyramidal neuronal networks interconnected by axoaxonal gap junctions. Neuroscience 92, 407–426.[CrossRef][Medline]

Vaknin G (1989). Major excitatory pathways in rat visual cortex. PhD Thesis, Kent State University, Kent, OH.

Ylinen A, Bragin A, Nadasdy Z, Jando G, Szabo I, Sik A et al. (1995). Sharp wave-associated high-frequency oscillation (200Hz) in the intact hippocampus: network and intracellular mechanisms. J Neurosci 15, 30–46.[Abstract]


    Acknowledgements
 
This research was supported by National Institute of Neurological Disorders and Stroke (NINDS) grant 2 R01 NS36981.




This article has been cited by other articles:


Home page
J. Neurophysiol.Home page
K. M. Rodgers, A. M. Benison, and D. S. Barth
Two-Dimensional Coincidence Detection in the Vibrissa/Barrel Field
J Neurophysiol, October 1, 2006; 96(4): 1981 - 1990.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
561/1/205    most recent
jphysiol.2004.075325v1
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Staba, R. J
Right arrow Articles by Barth, D. S
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Staba, R. J
Right arrow Articles by Barth, D. S


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS