Abstract
Gamma-frequency oscillations are prominent during active network states in the hippocampus. An intrahippocampal gamma generator has been identified in the CA3 region. To better understand the synaptic mechanisms involved in gamma oscillogenesis, we recorded action potentials and synaptic currents in distinct types of anatomically identified CA3 neurons during carbachol-induced (20–25 μm) gamma oscillations in rat hippocampal slices. We wanted to compare and contrast the relationship between excitatory and inhibitory postsynaptic currents in pyramidal cells and perisomatic-targeting interneurons, cell types implicated in gamma oscillogenesis, as well as in other interneuron subtypes, and to relate synaptic currents to the firing properties of the cells. We found that phasic synaptic input differed between cell classes. Most strikingly, the dominant phasic input to pyramidal neurons was inhibitory, whereas phase-coupled perisomatic-targeting interneurons often received a strong phasic excitatory input. Differences in synaptic input could account for some of the differences in firing rate, action potential phase precision, and mean action potential phase angle, both between individual cells and between cell types. There was a strong positive correlation between the ratio of phasic synaptic excitation to inhibition and firing rate over all neurons and between the phase precision of excitation and action potentials in interneurons. Moreover, mean action potential phase angle correlated with the phase of the peak of the net-estimated synaptic reversal potential in all phase-coupled neurons. The data support a recurrent mechanism of gamma oscillations, whereby spike timing is controlled primarily by inhibition in pyramidal cells and by excitation in interneurons.
Introduction
Gamma-frequency oscillations (30–100 Hz) are prevalent in active cortical networks. Although such network oscillations have been studied extensively in vivo, neither their underlying mechanisms nor functions are well understood (Paulsen and Sejnowski, 2006). In the hippocampal formation, gamma oscillations are often seen during theta activity (4–7 Hz), which occurs during exploratory behavior and rapid eye movement sleep (Buzsaki et al., 1983; Bragin et al., 1995; Csicsvari et al., 2003). Hippocampal gamma oscillations can also occur independently of theta oscillations (Penttonen et al., 1998; Csicsvari et al., 2003) and may be an inherent property of the local network.
In vitro models of gamma oscillations have provided some insights into the mechanisms of gamma oscillogenesis in the hippocampus. Analysis of carbachol (CCh)-induced gamma oscillations in hippocampal slices has identified an intrinsic gamma generator in area CA3 (Fisahn et al., 1998), which also appears to operate in vivo (Csicsvari et al., 2003). The CA3 network architecture is highly recurrent (Amaral et al., 1990): pyramidal cells (PCs) innervate PCs and interneurons, which in turn synapse onto PCs and other interneurons, making recurrent feedback inhibition a strong candidate mechanism for the generation of gamma oscillations in CA3. A PC–interneuron recurrent model of CCh-induced oscillations is attractive because the same mechanism is suggested by the spike order of PCs and interneurons during theta-related gamma oscillations in vivo (Csicsvari et al., 2003), and theta activity itself is believed to depend on septal cholinergic input (Stewart and Fox, 1990). Other gamma models implicate an independent interneuronal network as the primary rhythm generator (Whittington et al., 1995; Whittington and Traub, 2003; Bacci and Huguenard, 2006; Vida et al., 2006).
To better understand the synaptic mechanisms involved in gamma oscillogenesis and to distinguish between these two models, the present study aimed to record synaptic currents in anatomically identified neurons during cholinergically induced gamma oscillations. If CCh-induced oscillations are generated by a recurrent mechanism, interneurons would be expected to receive a dominant phasic excitatory drive and PCs would be expected to receive a dominant phasic inhibitory input, whereas if oscillations result from an interneuronal network, the dominant phasic input to both PCs and interneurons should be inhibitory. We therefore wanted to compare and contrast the relationship between excitatory and inhibitory synaptic currents in PCs and interneurons, focusing on PCs and perisomatic-targeting interneurons, cell types that have been suggested to underlie gamma oscillogenesis in the recurrent model (Fisahn et al., 1998; Gloveli et al., 2005a, b; Mann et al., 2005a, b). Furthermore, we wanted to relate the synaptic input to the firing properties of individual cells, which have been shown to differ between cell classes (Hájos et al., 2004). We found that synaptic input differed between cell classes and that differences in synaptic input could account for some of the differences in both the firing rate and phase coupling. Our results are consistent with a recurrent model for CCh-induced gamma oscillations in CA3.
Materials and Methods
Slice preparation.
Experiments were conducted in accordance with the United Kingdom Animals (Scientific Procedures) Act (1986), the Hungarian Act of Animal Care and Experimentation (1998, XXVIII, section 243/1998), and the guidelines of the institutional ethical codes. Wistar rats (postnatal day 13–20) were deeply anesthetized with isoflurane and decapitated, and the brain was removed and placed in ice-cold cutting solution containing the following (in mm): 252 sucrose, 2.5 KCl, 26 NaHCO3, 0.5 CaCl2, 5 MgCl2, 1.25 NaH2PO4, and 10 glucose. The cutting solution had been bubbled with 95% O2/5% CO2 (carbogen gas) for at least 30 min before use. Horizontal slices (350–400 μm thick) were cut and trimmed of most extrahippocampal areas. The slices were stored in interface conditions with artificial CSF (ACSF) containing the following (in mm): 126 NaCl, 3 KCl, 1.25 NaH2PO4, 2 MgSO4, 2 CaCl2, 26 NaHCO3; and 10 glucose for at least 60 min before being transferred individually to the recording chamber.
Electrophysiological recordings.
The recordings were conducted in a submerged-type recording chamber maintained at 29–33°C. Carbogen-bubbled ACSF was superfused at a flow rate of 3.5–5.5 ml/min. Oscillations were induced by bath application of 20–25 μm CCh (Sigma-Aldrich, St. Louis, MO). Patch pipettes (≈4–8 MΩ) filled with ACSF were used to record field oscillations and action potentials extracellularly. The field pipette was placed in the stratum pyramidale of the CA3 region, whereas a second pipette was used to detect action potentials from individual neurons identified using differential interference contrast microscopy (Axioskop; Zeiss, Jena, Germany). Action potentials were recorded for ∼60 s. The pipette was then withdrawn from the slice, and the cell was patched with a new pipette filled with intracellular solution containing the following: 138 mm K-gluconate, 3 mm CsCl, 10 mm phosphocreatine, 4 mm ATP, 0.4 mm GTP, 10 mm HEPES, 0.2 mm QX-314 [N-(2,6-dimethylphenylcarbamoylmethyl)-triethylammonium bromide], and 3 mg/ml biocytin. Seal resistance before whole-cell access was ≥1 GΩ. Whole-cell series resistance was in the range of 5–15 MΩ. All data were recorded with a Multiclamp 700A amplifier (Molecular Devices, Foster City, CA). During voltage-clamp recordings, series resistance, compensation was applied as necessary using the built-in series resistance compensation of the amplifier (12–17 kHz, 40–75%). Voltage measurements were not corrected for liquid–liquid junction potential. To record EPSCs and avoid interference from inhibitory events, cells were voltage clamped at a nominal holding potential of the estimated reversal potential for IPSCs (−65 to −60 mV of the uncorrected potential). Similarly, IPSCs were recorded at a nominal holding potential of the estimated EPSC reversal potential (0–20 mV of the uncorrected potential). Both field and unit recordings were low-pass filtered at 2 kHz using the built-in Bessel filter of the amplifier. Data were digitized at 5 kHz with either a PCI-6042E board (National Instruments, Austin, TX) and EVAN 1.3 software (courtesy of Prof. I Mody, University of California, Los Angeles, CA) or an ITC-16 board (Instrutech, Port Washington, NY) and Igor software (Wavemetrics, Lake Oswego, OR). All data were analyzed off-line using Igor Pro 5.01.
Event detection and analysis.
The field recording was further filtered using a digital, bidirectional, phase-conserving, low-pass filter at 1 kHz. Extracellular unit recordings were digitally, bidirectionally high-pass filtered at 40 Hz to isolate the spikes. Whole-cell recordings of postsynaptic currents (PSCs) were digitally, bidirectionally high-pass filtered at 1 Hz to filter out slow fluctuations in holding current. Post hoc digital filtering was done using Igor DSP filters (ftp://ftp.wavemetrics.net/IgorPro/User_Contributions/DSPFilters.ipf.zip).
To measure the power of the field oscillation, power spectral density analysis was performed on recordings of ∼60 s, using standard Igor Pro procedures. Time windows of ∼1.5 s with 50% overlap were multiplied by a Hanning window to minimize end-effects before the fast Fourier transform was performed. Gamma power was taken as the area under the power spectral density curve between 20 and 40 Hz.
To examine phase coupling, event times were converted into phase with respect to the background field oscillation using wavelet analysis (Torrence and Compo, 1998; Le Van Quyen et al., 2001). The Morlet wavelet transform of the field recording was examined between 10 and 45 Hz with scales chosen to reflect the equivalent Fourier frequency (see Fig. 1C–E) (Torrence and Compo, 1998). For each time point, the maximum of the wavelet transform magnitude was found, and the corresponding frequency was identified. The phase of the time point was defined in terms of this dominant frequency (Fig. 1D,E). Phase was defined such that −π was associated with the minimum of the oscillation, and a full cycle ran from −π to π (see Fig. 1F). To validate the wavelet technique, a peak-to-peak linear mapping phase analysis of action potentials was also used. A significant positive correlation was found between action potential phase as determined by the two methods (r = 0.83; p < 0.01; n = 35; angular–angular correlation; data not shown).
Phasic charge transfer was calculated as the mean of the integral from a baseline value over each cycle with the baseline defined as follows: (1) an overall initial baseline level estimate was obtained by computing the cycle-averaged event for the entire recording (see below). The maximum (for EPSCs) or minimum (for IPSCs) of this average event was taken as the initial baseline estimate. (2) For each 0.2 s epoch, the mean of all current values exceeding this initial baseline estimate was calculated. This mean was used as the baseline for the epoch. (3) If no current value was found to exceed the initial baseline estimate during the epoch, the baseline estimate would increment negatively for EPSCs and positively for IPSCs until such current values were found. This was repeated for all epochs of the entire recording (see Fig. 2B,C). Event times for PSCs were defined as time of peak current per cycle of the oscillation, and these were converted to wavelet phases as described.
The cycle-averaged field and cycle-averaged events (see Figs. 1G, 3B) were obtained by summing recordings over cycles between −π and π and dividing by the number of cycles. Cycles were scaled to span 2π radians. This phase is referred to as the cycle-averaged phase to distinguish it from the phase as defined by wavelet analysis.
The net apparent reversal potential (Esynrev) was estimated by first converting the cycle-averaged PSCs to excitatory (ge) and inhibitory (gi) conductances using the following equation: where Ie/i is the phasic excitatory/inhibitory current, Vh is the holding potential, and Ee/irev is the estimated reversal potential for the conductance of interest. The conductances were used to estimate the membrane potential of zero synaptic current (Esynrev) by solving the following equation: Event times for action potentials were defined as the time of crossing a voltage threshold set by visual inspection, and the event phase was taken as the wavelet phase of the dominant frequency at this event time. Analysis was done both by using all spikes per cycle and only the first spike per cycle. A significant positive correlation for phase-coupling variables (see below, Statistics) was found between values obtained with the two methods (action potential coupling strength: r = 0.99, p < 0.01, n = 43, Pearson's correlation; mean action potential phase: r = 0.99, p < 0.01, n = 35, angular–angular correlation; data not shown). The results reported here use all events detected.
To evaluate spike-timing precision, spike times (Δti) were defined as the time difference between the preceding minimum of the smoothed field oscillation and the time of threshold crossing. If ΔtI ≥ Ti/2, where Ti is the cycle period, then Δti was taken as the negative time difference from the subsequent minimum. Ti was calculated as the time between two consecutive minima in the field oscillation (see Fig. 4A). To enable comparison between the SDs for spike times (σΔt) and the corresponding SD for phase (σφ), σΔt was normalized by the cycle period as follows: where T̄ is the mean oscillation period. This quantity is comparable to the normalized SD of phase: where σφ is taken as the circular SD (see below, Statistics).
In some slices, phase-coupling strength of action potentials varied over time with the oscillation power such that phase-coupling strength decreased with decreasing oscillation magnitude (data not shown). Therefore, a reduced data set was obtained for the analysis of the relationships between synaptic events and firing properties. Cells for which the wavelet magnitude of the field oscillation changed by >2 SDs between spike train and PSC recordings were excluded from the reduced data set. Results are stated for both the full and reduced data set where applicable.
Statistics.
Coefficients of variation of the magnitude and frequency of the field oscillation (CVmag and CVfreq, respectively) were calculated from the SD of the maximum magnitude and the SD of the corresponding dominant frequency of the wavelet transform for each time point over the entire recording.
Circular statistics were used to determine the phase coupling of events. The phase of each event was represented by a unit vector pointing in the direction of the phase angle. The individual vectors were summed vectorally over the time epoch in question to give a resultant vector, R⃗. The direction of the resultant vector, φ, was used as the mean event phase for the epoch. The normalized length of R⃗, r = |R⃗|/n, where n is the number of events, served as a measure of the strength of the phase coupling. If the events are perfectly phase coupled, r = 1. If events are random, the vectors will cancel, and r → 0.
To determine the significance of the phase coupling, the Rayleigh probability of R⃗ was calculated using the following: where n is the number of spikes, and Z = n r2 (Fisher, 1993).
Before testing equality of means of linear variables, the Levene statistic was evaluated to test homogeneity of variances. ANOVA was used to compare means between homogenous distributions (p > 0.05, Levene statistic), and the more robust Welch statistic was used for non-homogenous distributions. The Bonferroni post hoc test was performed to find significant differences between distributions with equal variance. The Games-Howell post hoc test was performed to find significant differences between distributions with unequal variances. The multi-sample Watson-Williams test (Zar, 1999) was used to test equality of means of angular variables. The tests used in each case are specified in the text.
Before correlation tests for linear variables, the normality of a distribution was tested by the Kolmogorov–Smirnov test. To correlate variables from normal distributions (p > 0.05, Kolmogorov–Smirnov test), the Pearson's correlation coefficient was used. For variables from non-normal distributions, the more robust Spearman's rank correlation coefficient was used.
To correlate angular–angular variables, the following parametric angular–angular correlation coefficient was used: where the ith pair of data are denoted as ai, bi. The confidence interval was determined by a jackknifing procedure (Zar, 1999). All correlation coefficients are quoted as r, and tests used are specified in the text.
The circular SD was taken as follows: where r is the phase-coupling strength (Zar, 1999). Values are given as mean ± SEM, unless stated otherwise.
Anatomical identification.
After intracellular recording and biocytin filling, the slice was fixed in 4% paraformaldehyde in 0.1 m phosphate buffer (PB), pH 7.4, for at least 30 min, followed by washout with PB several times and incubation in cryoprotecting solution (30% sucrose in 0.01 m PB) for 2 h. Slices were then freeze-thawed three times above liquid nitrogen and treated with 1% H2O2 in PB for 15 min to reduce endogenous peroxidase activity. Filled cells were visualized using avidin-biotinylated horseradish peroxidase complex reaction (Vector Laboratories, Burlingame, CA) with nickel-intensified 3,3′-diaminobenzidine (Sigma-Aldrich, St. Louis, MO) as chromogen giving a dark reaction product. After dehydration and embedding in Durcupan (Sigma-Aldrich), neurons were identified based on their dendritic and axonal arborization, and representative cells were reconstructed with the aid of a drawing tube using a 40× objective (see Fig. 2).
Results
Fast network oscillations were induced in the CA3 region of horizontal hippocampal slices by bath application of CCh (20–25 μm). The network oscillations were monitored by recording the field potential using an ACSF-filled micropipette placed in the stratum pyramidale. Action potentials of individual neurons were recorded extracellularly from visually selected neurons using a second ACSF-filled pipette. This second pipette was then withdrawn, and the EPSCs and IPSCs were recorded in whole-cell voltage-clamp mode with a pipette containing intracellular solution. Biocytin in the intracellular solution was used for post hoc anatomical identification of recorded cells. Of 74 neurons, 43 cells were successfully recovered and anatomically identified based on their dendritic and axonal arborization. Only anatomically identified neurons were included in the following analysis. Data reported excludes outliers (defined as >3 SDs from the mean of the measure of interest for each cell type). Two cells were excluded from additional analysis on this basis.
CCh-induced fast network oscillations
To characterize the field oscillations, extracellular recordings of ∼1 min (Fig. 1A) were acquired and used to obtain a power spectrum (Fig. 1B). Calculation of the power spectral density function revealed a peak in the beta–gamma frequency range with a mean peak frequency over all recordings of 28.1 ± 0.6 Hz (fPSD) and mean gamma power of 1782.9 ± 418.7 μV2 (n = 43). Wavelet analysis using a Morlet wavelet basis (Fig. 1C) was used to extract the magnitude and phase of different frequency components (10–45 Hz) of the field oscillation (Fig. 1D–F). Both the maximum field wavelet magnitude and the corresponding dominant wavelet frequency fluctuated over time (CVmag = 24.4 ± 1.1%; CVfreq = 14.1 ± 1.2%; see Materials and Methods). Inspection of Figure 1, A and D, suggests that the fluctuations in the oscillation magnitude might be rhythmical. Indeed, in 28 of the 43 recordings, a power spectrum analysis of the maximum wavelet magnitude over time revealed a low frequency peak (2.8 ± 1.6 Hz; n = 28; data not shown).
To further validate that wavelet analysis provided an appropriate definition of phase, the cycle-averaged field potential was calculated (Fig. 1G) (see Materials and Methods). This had a regular form, providing qualitative support for the method. Moreover, the frequency corresponding to the period of these cycle-averaged events (29.6 ± 0.7 Hz; n = 43) was similar to and correlated with fPSD (r = 0.95; p < 0.01, Pearson's correlation; n = 43).
Comparison of synaptic input in different cell types
To understand the synaptic mechanisms underlying gamma oscillations, PSCs were recorded from PCs (Fig. 2Ai) (n = 17) and perisomatic-targeting interneurons (basket and axo-axonic cells; n = 9) (Fig. 2Aii). Recordings were also obtained from dendritic-targeting interneurons [n = 7, comprising OLM cells, interneurons with a dendritic tree in the stratum oriens projecting to the stratum lacunosum moleculare (n = 3), and radiatum cells/radiatum-lacunosum-moleculare cells (n = 4)] (Fig. 2Aiii) and interneuron-selective interneurons [as defined by Hájos et al. (2004); n = 7] (Fig. 2Aiv) for comparison. Cells were voltage clamped at the estimated inhibitory and excitatory reversal potentials to record EPSCs and IPSCs, respectively (Fig. 2B,C) (see Materials and Methods). This information could reveal differences in the dominant input to cells and thus suggest a functional network architecture underlying gamma oscillations. Hereafter, any reference to excitation or inhibition refers to phasic excitation or phasic inhibition.
If phasic excitation drives cell firing, differences in the excitatory input to the cells may lead to differences in the recruitment of the various subpopulations of neurons. Phasic excitatory charge transfer (Qe) differed significantly between cell types (Fig. 2D) (p < 0.01, Welch statistic; n = 40). Qe was smallest in PCs (Qe = 0.56 ± 0.063 pC; n = 17) and was, on average, largest in perisomatic-targeting interneurons (Qe = 3.1 ± 0.91 pC, n = 9; dendritic-targeting interneurons: Qe = 0.94 ± 0.267 pC, n = 7; interneuron-selective interneurons: Qe = 2.2 ± 0.42 pC, n = 7).
The balance of excitation and inhibition was investigated by comparison of the ratio of phasic excitatory to inhibitory charge (Qe/Qi) between cell types. Phasic inhibition exceeded phasic excitation in all PCs (Qe/Qi = 0.23 ± 0.049; n = 17), whereas the dominant input to interneurons was often excitatory with perisomatic-targeting interneurons showing the largest excitation to inhibition ratio (Qe/Qi = 1.2 ± 0.25, n = 7; dendritic-targeting interneurons: Qe/Qi = 0.86 ± 0.14, n = 7; interneuron-selective interneurons: Qe/Qi = 1.0 ± 0.05, n = 4) (Fig. 2E). The difference in Qe/Qi between PCs and interneuron classes was significant (Fig. 2E) (p < 0.01, Welch statistic and ANOVA with post hoc Games-Howell test; n = 35).
Although the significant difference in Qe/Qi between PCs and interneurons could result from a larger phasic excitatory input to interneurons, a larger phasic inhibitory input to PCs could also contribute to this difference. Comparison of Qi between cell types revealed a significantly larger Qi in PCs (Qi = 3.2 ± 0.30 pC; n = 17) than in interneurons (Fig. 2F) (n = 35; p < 0.01, ANOVA; perisomatic-targeting interneurons: Qi = 1.5 ± 0.43 pC, n = 7; dendritic-targeting interneurons: Qi = 1.2 ± 0.30 pC, n = 7; interneuron-selective interneurons: Qi = 1.8 ± 0.46 pC, n = 4). Interestingly, Qi was found to correlate with field magnitude (r = 0.58; p < 0.01, Pearson's correlation; n = 35) (see supplemental Fig. 1, available at www.jneurosci.org as supplemental material), and this might account for some of the observed difference in Qi between cell type (p < 0.05, ANOVA with field magnitude as covariate; n = 35; partial η2celltype = 0.25; partial η2fieldmag = 0.16). A correlation between charge and field magnitude was not seen for excitatory charge (r = −0.21; p = 0.80, Spearman's correlation; n = 40) (see supplemental Fig. 1, available at www.jneurosci.org as supplemental material). The data support the conclusion that both excitatory and inhibitory inputs differ between cell types.
The majority of cells received significantly phase-coupled (p < 0.05, Rayleigh probability) excitatory and inhibitory synaptic input as measured by the phase of the peak current recorded in each cycle (40 of 40 and 36 of 37, respectively). The coupling strength of the peak of the EPSC (re) did not differ significantly between cell types (Fig. 3Ai) (p = 0.27, ANOVA; n = 40), whereas the coupling strength of inhibitory inputs (ri) was significantly larger in PCs (ri = 0.83 ± 0.024; n = 17) than in interneurons (Fig. 3Aii) (p < 0.01, Welch statistic; n = 37; perisomatic-targeting interneurons: ri = 0.63 ± 0.074, n = 9; dendritic-targeting interneurons: ri = 0.46 ± 0.122, n = 7; interneuron-selective interneurons: ri = 0.57 ± 0.059, n = 4).
In all cells, the phase of the peak of the cycle-averaged EPSC preceded the peak of the cycle-averaged IPSC (Fig. 3B). To capture the temporal interaction between inhibitory and excitatory synaptic conductances, the effective synaptic conductance was estimated by calculation of the net apparent synaptic reversal potential Esynrev (see Materials and Methods). The shape of the Esynrev curve differed between PCs and interneurons (Fig. 3C), with the half-width (measured in radians) of the peak of the Esynrev curve being significantly narrower in PCs (half-width of 0.79 ± 0.157; n = 17) than in interneurons (Fig. 3D) (p < 0.01, ANOVA with Bonferroni post hoc test; n = 35; perisomatic-targeting interneurons: half-width of 2.6 ± 0.32, n = 7; dendritic-targeting interneurons: half-width of 2.6 ± 0.25, n = 7; interneuron-selective interneurons: half-width of 2.7 ± 0.23, n = 4), primarily as a result of the dominant inhibitory input received by PCs.
Firing properties of anatomically identified neurons during fast network oscillations
To characterize the firing properties of the cells, action potentials were recorded extracellularly for ∼1 min (Fig. 4A) before whole-cell recording with another pipette (see Materials and Methods for details). The firing properties of the cells were quantified using three parameters: mean firing rate, action potential coupling strength (rAP), and mean action potential phase (φAP; see Materials and Methods).
Firing rate
PCs fired at a significantly lower rate (6.0 ± 0.82 Hz; n = 17) than interneurons (Fig. 4B) (p < 0.01, Welch statistic and ANOVA with Games-Howell post hoc test; n = 42), whereas perisomatic-targeting (rate, 19 ± 2.3 Hz; n = 10) and interneuron-selective (rate, 21 ± 3.2 Hz; n = 8) interneurons fired at slightly higher rates than dendritic-targeting interneurons (rate, 16 ± 2.8 Hz; n = 7), in agreement with previous results (Hájos et al., 2004).
Phase versus time
The timing of action potentials after synaptic excitation may be the result of the intrinsic properties of the neuron, recurrent inhibitory synaptic input from the network, or both. If intrinsic properties play a dominant role in governing spike timing, then the absolute spike time after the arrival of the excitatory synaptic input should remain constant independent of the oscillation frequency. Alternatively, spike timing might be more tightly controlled by both the excitatory and inhibitory synaptic input. The field oscillation emerges from the total network activity, and therefore if spike timing is predominantly controlled by both the excitatory and inhibitory synaptic input, variability in spike phase relative to the field oscillation should be smaller than variability in absolute time. This provided the motivation to compare spike-timing precision in terms of absolute time and phase. If the oscillation frequency remains constant, the two measures should give comparable results. However, if the frequency varies over the recording, the results will differ. The SD of spike times relative to the field oscillation (σtnorm) and the SD of the spike phases (σtnorm; see Materials and Methods) were used to evaluate precision (Mainen and Sejnowski, 1995). Overall, σtnorm was significantly smaller than σtnorm (Fig. 4C) (p < 0.05, paired-samples t test; n = 42), indicating that both excitatory and inhibitory synaptic inputs are important in the control of spike timing of CA3 neurons. However, σtnorm > σφnorm in cells that were weakly phase coupled (rAP<0.2; see below, Action potential phase-coupling strength).
Action potential phase-coupling strength
Thirty-eight of the 42 cells were significantly phase coupled (p < 0.05, Rayleigh probability). The degree of phase coupling varied between significantly phase-coupled cells. The distribution of rAP appeared bimodal with a minimum around rAP = 0.2 (Fig. 4D, inset). This was set as the cutoff value to classify cells as strongly or weakly phase coupled. Inspection of the spike-phase histograms confirmed this as an appropriate choice (Fig. 4D) [see also Hájos et al. (2004), their Fig. 7C]. The phase histogram of strongly phase-coupled cells (rAP > 0.2) exhibited a clear peak in the phase histogram, whereas the phase histograms of weakly phase-coupled cells did not. The majority of PCs, perisomatic-targeting, and interneuron-selective interneurons were strongly phase coupled (rAP > 0.2; 17 of 17, 8 of 10, and 6 of 8, respectively), whereas this was the case for only 3 of the 7 dendritic-targeting interneurons (Fig. 4D). It has been reported that the precision of phase coupling is lower in PCs than in the majority of interneurons (Hájos et al., 2004). Consistent with this, rAP of PCs tended to be smaller than that of strongly coupled perisomatic-targeting interneurons (0.60 ± 0.04, n = 17 vs 0.66 ± 0.08, n = 8). Action potential phase-coupling strength was, on average, weakest in dendritic-targeting interneurons, consistent with previous findings (Hájos et al., 2004). The difference in rAP between cell types failed to reach statistical significance (Fig. 4E) (p = 0.1, Welch statistic; n = 42).
Mean action potential phase angle
Strongly phase-coupled PCs fired near the minimum of the field oscillation (φAP = −2.1 ± 0.09; n = 17; phase measured in radians), whereas strongly phase-coupled interneurons tended to fire later in the cycle. This difference was significant (Fig. 4F) (p < 0.01, Watson-Williams test; n = 34). Phase-coupled dendritic-targeting interneurons tended to fire later in the cycle (φAP = −1.4 ± 0.36; n = 3) than perisomatic-targeting (φAP = −1.7 ± 0.13; n = 8) and interneuron-selective (φAP = −1.7 ± 0.16; n = 6) interneurons, but this was not statistically significant (Fig. 4F) (p = 0.20 and p = 0.34, respectively, Watson-Williams test). These findings are consistent with the findings of Hájos et al. (2004).
Thus, different cell types display different firing properties during CCh-induced oscillations with respect to firing rate, action potential coupling strength, and mean action potential phase angle.
Synaptic input and firing properties during fast network oscillations
The relationship between synaptic input from the network and firing properties was investigated next in terms of the three parameters discussed, namely firing rate, action potential phase-coupling strength, and mean action potential phase angle.
Synaptic input and firing rate
If phasic excitatory input drives cells to fire, the firing rate would be expected to vary with the magnitude of the excitatory input. Hence, the relationship between Qe and spike rate was investigated. Over all cells, a significant positive correlation was found between rate and excitatory charge (Fig. 5A) (full set: r = 0.51, p < 0.01, n = 39; reduced set: r = 0.61, p < 0.01, n = 31; Spearman's correlation; see Materials and Methods). We have shown that perisomatic-targeting and interneuron-selective interneurons fire at a higher rate than PCs and dendritic-targeting cells (Fig. 4B). If the firing rate is controlled by excitatory synaptic input, perhaps differences in the amount of excitation received by different cell classes could account for differences in the firing rate. Indeed, cell types with high firing rates received, on average, more excitation than cell types with low firing rates (Figs. 2D, 5B).
If inhibition suppresses firing, a negative correlation would be expected between Qi and the firing rate. This relationship was seen in the reduced data set but failed to reach statistical significance in the full data set (Fig. 5C) (reduced set: r = −0.40, p < 0.05, n = 27; full set: r = −0.32, p = 0.07, n = 34; Pearson's correlation). A stronger correlation was found between the firing rate and Qe/Qi (Fig. 5D) (full set: r = 0.73, p < 0.01, n = 34; reduced set: r = 0.79, p < 0.01, n = 27; Pearson's correlation).
Synaptic input and action potential phase-coupling strength
Whereas phase coupling is a robust phenomenon in a large proportion of CA3 neurons, the strength of the phase coupling differed between cells. Could differences in synaptic input account for this variation? We have shown that cell classes differ with respect to the dominant input received from the network (Fig. 2D). Because cell classes also differ in their intrinsic properties, this led to the question of whether cell classes process synaptic input differently. To investigate this, the relationship between Qe/Qi and rAP was assessed for each cell type studied. All PCs received a dominant inhibition and showed strong phase coupling. In perisomatic-targeting interneurons, Qe/Qi differed between strongly and weakly phase-coupled cells. In this cell group, the dominant input in strongly phase-coupled neurons tended to be excitatory or the inputs were approximately equal such that, on average, Qe/Qi > 1, whereas weakly phase-coupled cells received a dominant inhibitory input (Qe/Qi < 1). This relationship was not observed in other cell types. Dendritic-targeting interneurons showed a broad range of Qe/Qi and variable rAP, whereas the ratio of excitation to inhibition in interneuron-selective interneurons was approximately equal, with Qe/Qi values closer to 1 independent of the value of rAP (Fig. 6A).
If the ratio of excitation to inhibition cannot fully explain the coupling strength of the cell, it may be that the total synaptic input is more important. Indeed, it was found that strongly phase-coupled perisomatic-targeting interneurons received, on average, a larger synaptic drive (both inhibitory and excitatory) than weakly phase-coupled cells. This result extended to the other interneuron classes studied, in which a significant correlation was found between Qe and rAP and between Qi and rAP in the full data set (Fig. 6B) (full set: r = 0.66, p < 0.01, n = 23, and r = 0.53, p < 0.05, n = 18, respectively; reduced set: r = 0.67, p < 0.01, n = 19, and r = 0.3, p = 0.27, n = 15, respectively; Pearson's correlation).
If synaptic input controls spike timing, the phase-coupling strength of the synaptic input would be expected to contribute toward the strength of action potential phase coupling. Indeed, in interneurons, a positive correlation was found between rAP and re (Fig. 6Ci) (full set: r = 0.77, p < 0.01, n = 23; reduced set: r = 0.73, p < 0.01, n = 19; Pearson's correlation), and a weaker correlation was found between rAP and ri (Fig. 6Cii) (full set: r = 0.56, p < 0.01, n = 21; reduced set: r = 0.47, p = 0.057, n = 17; Pearson's correlation), suggesting that it is predominantly the timing precision of excitatory input that controls spike-timing precision in interneurons. In PCs, the correlation between rAP and re was weaker than in interneurons (Fig. 6Di) (full set: r = 0.50, p < 0.05, n = 16; reduced set: r = 0.37, p = 0.23, n = 12; Pearson's correlation), whereas no significant correlation was found between rAP and ri (Fig. 6Dii) (full set: r = 0.25, p = 0.35, n = 16; reduced set: r = 0.14, p = 0.66, n = 12; Pearson's correlation).
Because cells receive excitatory and inhibitory inputs concomitantly, the temporal interaction of excitation and inhibition may also contribute toward the action potential coupling strength. If the net resultant synaptic input alone controls spike time, then the width of the peak of Esynrev should give a time window for action potential firing. A correlation analysis of the half-width of Esynrev and rAP revealed no such relationship in interneurons (Fig. 6Ciii) (full set: r = −0.03, p = 0.90, n = 18; reduced set: r = −0.29, p = 0.30, n = 15; Pearson's correlation) or PCs (Fig. 6Diii) (full set: r = 0.03, p = 0.93, n = 16; reduced set: r = 0.03, p = 0.93, n = 12; Pearson's correlation).
We conclude that both the precision and strength of synaptic input contribute toward controlling rAP in interneurons.
Synaptic input and mean action potential phase angle
We have shown that different classes of cell preferentially fire at different phases of the gamma oscillation. Could differences in the timing of the synaptic input account for the differences in mean action potential phase angle? The relationship between φAP and synaptic input was investigated next. φAP in strongly phase-coupled cells was found to either follow (10 of 26) or precede (16 of 26) the peak of the cycle-averaged EPSC and always preceded the peak of the cycle-averaged IPSC (Fig. 7A). When the effects of excitation and inhibition were combined through the estimation of Esynrev, it was found that the peak of Esynrev (φsyn) preceded φAP in the vast majority of strongly phase-coupled cells (24 of 26). In PCs, φsyn reflected the phase of least inhibition, whereas in phase-coupled perisomatic-targeting interneurons and other interneurons receiving a large excitatory input, it was associated with the peak of the excitation (Fig. 7B).
If the timing of excitatory input controls spike timing in all cells, it might be expected that excitation arrives earlier in PCs than in interneurons, because they fire at earlier phases. However, the phase of peak excitation was significantly later in phase-coupled PCs (φE = −1.1 ± 0.08; n = 16) than in phase-coupled interneurons (Fig. 7C) (full set: p < 0.01; perisomatic-targeting interneurons: φE = −1.8 ± 0.25, n = 7; dendritic-targeting interneurons: φE = −1.7 ± 0.04, n = 3; interneuron-selective interneurons: φE = −1.5 ± 0.15, n = 5; reduced set: p < 0.01, Watson-Williams test; n = 26). This suggests that it is not excitation alone that governs precise action potential timing. Indeed, when the effect of excitation and inhibition were combined, φsyn was found to be significantly earlier in phase-coupled PCs (φsyn = −2.7 ± 0.09; n = 16) than in phase-coupled interneurons (Fig. 7D) (full set: p < 0.01; perisomatic-targeting interneurons: φsyn = −2.2 ± 0.20, n = 5; dendritic-targeting interneurons: φsyn = −1.7 ± 0.05, n = 3; interneuron-selective interneurons: φsyn = −2.1 ± 0.86, n = 2; reduced set: p < 0.01, Watson-Williams test; n = 20), resulting in a significant correlation between φsyn and φAP over all cells (Fig. 7E) (full set: r = 0.46, p < 0.01, n = 26; reduced set: r = 0.61, p < 0.01, n = 20; angular–angular correlation).
In conclusion, we have shown that differences in synaptic input are important in controlling the distinct firing properties of neurons during gamma oscillations.
Discussion
It has been reported previously that different subclasses of CA3 neurons differ in their firing properties during cholinergically induced oscillations (McMahon et al., 1998; Hájos et al., 2004). The present study investigated PSCs in anatomically identified CA3 neurons during CCh-induced network oscillations and related the PSCs to the firing properties of the cells. Our main findings were as follows: (1) synaptic input differed between cell classes during cholinergically induced gamma oscillation; and (2) features of synaptic input could be used to predict firing properties, namely firing rate and spike timing (quantified by phase-coupling strength and mean phase angle).
Variability in CCh-induced field oscillation
Wavelet analysis revealed a rhythmicity in the magnitude of the field oscillation in approximately half of the recordings. This could arise from the interference of oscillators, oscillating at similar, but unequal frequencies, or, alternatively, from slow membrane potential oscillations. Cunningham et al. (2003) suggested that Ih-dependent membrane potential oscillations of stellate cells can produce a theta-frequency modulation of gamma oscillations in the entorhinal cortex. Theta-frequency membrane potential oscillations have been recorded in hippocampal PCs (Leung and Yim, 1991), raising the possibility that a similar mechanism may generate the low frequency modulation seen here. These possibilities were not explored further.
Different synaptic input in distinct cell classes and local network architecture
Cell classes differed in their dominant synaptic input. Inhibitory synaptic currents exceeded excitatory synaptic currents in all PCs, whereas excitation was, on average, stronger and often dominated over inhibition in phase-coupled perisomatic-targeting interneurons. In fact, these phase-coupled neurons received the largest excitatory input of all cell types. In line with these results, anatomical investigations have revealed that perisomatic-targeting interneurons containing the calcium binding protein parvalbumin received the largest number of excitatory synapses among all inhibitory cells (Gulyas et al., 1999).
Based on current source density analysis and spike-time measurement, it has been suggested that hippocampal gamma oscillations arise from the disynaptic interaction of excitatory PCs and inhibitory perisomatic-targeting interneurons (Mann et al., 2005a, b). Our data provide the first direct evidence of the synaptic currents in these neuronal populations (Fig. 2) and is consistent with a recurrent inhibitory model of gamma oscillations (Freeman, 1968). This model is supported by four of our findings: (1) the relative spike times of PCs and interneurons (Fig. 4); (2) the relative timing of excitatory and inhibitory events (Fig. 3B); (3) the differences in excitatory and inhibitory charge in different cell types (Fig. 2); and (4) the differences in coupling strength of the IPSC peak in different cell types. The fact that ri was larger in PCs than in interneurons (Fig. 3A) is in agreement with a recurrent model of gamma oscillogenesis in which a strongly phase-coupled interneuron population provides a perisomatic source of inhibition to PCs, whereas interneurons receive inhibition from numerous subpopulations of interneurons with varied phase coupling. Interestingly, Mann et al. (2005b) observed that the current source was not confined to the soma and axon initial segment but extended to the proximal dendritic region. In our sample, we have recorded from two inhibitory neurons that arborize in the proximal dendritic region (data not included), that showed phase coupling, and received a large excitatory current, comparable to perisomatic-targeting interneurons, suggesting that this population might also be actively involved in oscillogenesis.
Synaptic input controls spike output
Although differences in synaptic input exist between cell classes, it appears that the cell classes also differ in their responses to similar synaptic input (Fig. 6A). This raises the question of the importance of synaptic input versus the intrinsic properties of a cell in the control of firing properties during gamma oscillations.
Firing rate
It was found that PCs fire at a significantly lower rate than interneurons. Although intrinsic conductances, such as the one underlying Ih, can contribute to setting a firing rate preference (Maccaferri and McBain, 1996; Pike et al., 2000, Hu et al., 2002), differences in synaptic input could also account for differences in the firing rate. Indeed, over all cells, a stronger positive correlation was found between the firing rate and Qe/Qi than between the firing rate and Qe, suggesting that both phasic excitation and inhibition contribute toward controlling the firing rate. Furthermore, cell types with high firing rates received larger phasic excitatory synaptic input than cell types with lower firing rates, supporting the conclusion that although intrinsic properties set a firing rate preference, synaptic input also clearly contributes to the differences in firing rates of different cell types.
Spike-timing precision
The ability of neurons to fire action potentials reliably and precisely has been shown for artificially generated input (Mainen and Sejnowski, 1995; Hunter et al., 1998; Fellous et al., 2001; Haas and White, 2002). The CCh model allows for a measure of precision under more physiological conditions. A precise spike output could be the result of intrinsic membrane properties, net synaptic input, or an interaction between the two. The SD of the spike phase distribution was smaller than the equivalent measure for absolute spike time in strongly phase-coupled cells. The opposite relationship was seen in weakly phase-coupled cells. This suggests that synaptic input provides a considerable contribution to the control of spike timing in strongly phase-coupled neurons.
Action potential phase-coupling strength
We found that all interneuron types, including dendritic-targeting interneurons, fire phase coupled to the oscillation when they receive sufficiently strong excitatory synaptic input (Fig. 6B). Moreover, perisomatic-targeting interneurons receiving a dominant excitation were strongly phase coupled, whereas those with Qe/Qi < 1 tended not to be. However, the relationship of Qe/Qi to rAP seen in perisomatic-targeting interneurons was not observed in other cell types. The apparent lack of relationship between phasic charge ratio and rAP in dendritic-targeting and interneuron-selective interneurons could be attributable to the relatively small number of neurons sampled, or to the fact that these cells were treated as two populations when, in fact, they are composed of several distinct subpopulations (radiatum and OLM, and IS-1, IS-2, and IS-3, respectively) (Freund and Buzsaki, 1996). Differences in intrinsic properties might also contribute toward the observed differences between cell types in the relationship between Qe/Qi and rAP.
In interneurons, a positive correlation was found between the phase-coupling strength of synaptic current and action potentials, suggesting that a more precise input leads to a more precise output. Moreover, the correlation was stronger for excitatory inputs, corroborating the conclusion that interneuron spike phase precision is primarily controlled by PCs. In PCs, no strong correlation was found between the precision of the peak phase of synaptic input and rAP, suggesting that other factors control the precise time of action potential generation in this cell class.
If the net synaptic input sets the time window for action potential firing, and because PCs had a narrow Esynrev peak compared with interneurons, it might be expected that the coupling of PC firing to the oscillation be stronger than interneuron phase coupling, because PCs had a narrow Esynrev peak compared with interneurons. However, this was not the case, because phase coupling of PCs was weaker than that of phase-coupled interneurons, in line with previous findings (Hájos et al., 2004). Consistent with these findings, Fricker and Miles (2000) demonstrated that most CA1 interneurons fire more precisely in response to excitatory input than do PCs. They attributed this difference to differences in intrinsic properties between these cell types. Similarly, differences in intrinsic properties, as well as differences in the somato-dendritic distribution of synaptic inputs, might account for some of the differences in spike phase precision observed here. Future experiments will be needed to dissociate the relative contributions of synaptic input and intrinsic properties toward spike output.
Mean action potential phase angle
Action potentials occurred, on average, earlier in PCs than in interneurons, yet the peak of synaptic excitation occurred later in PCs. However, φsyn, which approximates the time when phasic excitation is most dominant relative to phasic inhibition, occurred significantly earlier in PCs than in interneurons, suggesting that both excitation and inhibition are important in the control of spike timing. In PCs, φsyn occurred when inhibition is at a minimum, whereas in perisomatic-targeting interneurons, the main target of PCs, φsyn occurred near the peak of the excitation. Hence, during gamma oscillations, PC spike time is controlled by inhibition, whereas excitation is important in the control of perisomatic-targeting interneuron spike timing. It has been reported, however, that the average firing rate of PCs in vivo is significantly lower than that observed here (Csicsvari et al., 2003). Thus, even with a higher recurrent connectivity in the intact brain, Qe might be overestimated in our study, and PC control of perisomatic-targeting interneuron spike timing may be complemented by other mechanisms in vivo, such as interneuronal gap junctions (Traub et al., 2001).
In conclusion, our findings support a model of cholinergically induced gamma oscillations mediated by an interaction of PCs and perisomatic-targeting interneurons. The results highlight the importance of synaptic input in the control of firing in both PCs and interneurons. These results shed light on the mechanism underlying hippocampal gamma oscillations in vitro. A thorough understanding of the mechanism of in vitro oscillogenesis may provide a first step toward understanding the mechanism and function of these oscillations in the intact brain.
Footnotes
-
This work was supported by The Wellcome Trust, the Biotechnology and Biological Sciences Research Council, the British Council, the Alzheimer's Research Trust, and Országos Tudományos Kutatási Alapprogramok Grant T049517. N.H. was awarded the Bolyai János Scholarship. We are grateful for the generous support of Prof. Tamás F. Freund. We also thank János Pálhalmi for contributions during the early stages of this work and Katalin Lengyel for excellent technical assistance.
- Correspondence should be addressed to Dr. Ole Paulsen, Department of Physiology, Anatomy and Genetics, University of Oxford, Sherrington Building, Parks Road, Oxford OX1 3PT, UK. ole.paulsen{at}physiol.ox.ac.uk