Abstract
A hippocampal place cell fires at an increasingly earlier phase in relation to the extracellular theta rhythm as a rodent moves through the place field. The present report presents a compartment model of a CA1 pyramidal cell that explains the increase in amplitude and the phase precession of intracellular theta oscillations, with the assumption that the cell receives an asymmetric ramp depolarization (<10 mV) in the place field and rhythmic inhibitory and/or excitatory synaptic driving. Intracellular subthreshold membrane potential oscillations (MPOs) increase in amplitude and frequency, and show phase precession within the place field. Theta phase precession and MPO power and frequency increase in the place field are caused by a shift in excitatory–inhibitory response, intrinsic theta-frequency resonance, and intrinsic oscillations that depend on voltage-dependent persistent Na+ and slowly inactivating K+ currents, but not on Ih. Phase precession is diminished when theta-frequency resonance is decreased. Simulated spikes fire near the peak of MPOs and precess similarly as the MPOs. The phase of the MPOs/spikes codes for distance in a one-dimensional place field, and phase precession is only weakly dependent on firing rate, running speed, or the duration needed to cross the place field. In addition, phase precession within the place field resumes quickly after disruption by maximal afferent pulse stimulation.
Introduction
A place cell fires at a high rate when an animal runs through a particular location in the environment (O'Keefe and Nadel, 1978; Muller, 1996), with action potentials occurring progressively earlier in the theta-frequency cycle (O'Keefe and Recce, 1993; Skaggs et al., 1996). Theta (4–10 Hz) rhythm in the hippocampus and related structures (Vanderwolf, 1988; Stewart and Fox, 1990; Bland and Colom, 1993; Buzsáki, 2002) serves as a reference for a phase code (Hopfield, 1995) of an animal's location at different running velocities (Huxter et al., 2003).
Harvey et al. (2009) showed three main “subthreshold signatures” of place cells by whole-cell recording of hippocampal neurons in head-fixed mice running through a virtual reality track: (1) an asymmetric ramp-like depolarization, (2) an increase in the amplitude of intracellular theta oscillations, and (3) a phase precession of the peaks of the membrane potential oscillations (MPOs) with the extracellularly recorded theta rhythm. Since spikes fire near the peak of the MPOs (Harvey et al., 2009), spike precession in the place field is determined by MPO precession. Furthermore, the increase in amplitude of the MPOs constrains the participation of the rhythmic EPSPs in intracellular theta generation, since rhythmic EPSPs decrease with depolarization while rhythmic IPSPs increase (Leung and Yim, 1986; Fox, 1989; Ylinen et al., 1995; Bland et al., 2002).
Various models have been proposed for phase precession of place cell firing. Linear interference models propose MPOs of higher frequency at the dendrites compared with the soma (Burgess and O'Keefe, 1996; Lengyel et al., 2003). Other models are not explicit as to the increase in MPO frequency but propose that increasing dendritic excitation over somatic inhibition fires somatic spikes earlier in each theta cycle (Kamondi et al., 1998; Magee, 2001; Harris et al., 2002; Losonczy et al., 2010). In addition, a slow ramp-like depolarization would decrease the latency of spike onset (Harris et al., 2002; Mehta et al., 2002; McLelland and Paulsen, 2009). Some models emphasize that a neural network is critical for phase precession (Jensen and Lisman, 1996; Tsodyks et al., 1996; Wallenstein and Hasselmo, 1997; Bose and Recce, 2001; Hasselmo and Eichenbaum, 2005; Geisler et al., 2007, 2010), but few models address fully the intrinsic properties of single neurons.
Intrinsic theta-frequency MPOs were induced by a small depolarization in entorhinal cortex stellate cells (Alonso and Llinás, 1989) and hippocampal CA1/CA3 pyramidal cells (Leung and Yim, 1991). In hippocampal pyramidal cells, a persistent Na+ current INaP and a slowly inactivating K+ current IKS or IM are required for generating low-threshold MPOs (Leung and Yim, 1991) and somatic resonance at theta frequency (Leung and Yu, 1998; Hu et al., 2002, 2007, 2009).
In the present report, electrophysiological responses of a CA1 pyramidal cell are simulated by a compartment model with INaP and IKS, in addition to other intrinsic currents. When a depolarization ramp and rhythmic synaptic driving are assumed, the model shows phase precession, and frequency and amplitude increase of the subthreshold MPOs.
Materials and Methods
Compartment neuron model.
Electrophysiological responses of a CA1 pyramidal cell were simulated using a compartment model represented by first-order differential equations. The main model consists of 38 compartments (Fig. 1). Each compartment is assumed to be a cylindrical cable, and the ith compartment is described by the following:
where Vi−1, Vi, and Vi+1 are membrane potential of the (i − 1)th, ith, and (i + 1)th compartments, respectively; Gmi and Ci are membrane conductance and capacitance of the ith compartments, respectively; Gli,i+1 is longitudinal conductance between the ith and the (i + 1)th compartments; and El is leak current reversal potential, normally −70 mV but adjusted to keep resting membrane potential near −66 mV.
Model of CA1 pyramidal cell consists of 38 compartments representing different parts of the cell. An axon [compartment (C) 18, 37, 38], two basal dendritic branches (C1–8, C9–16), and an apical dendritic branch (C19–36) representing the combined apical dendritic tree, originate from the soma (C17). Persistent Na+ and slow K+-inactivating currents are assumed only at the proximal compartments (C17–19; colored gray), while thick circles indicate the excitable compartments with Na+ action potential current and IK. The dendrites have Ih, IA conductance that increase distally. For the present simulations, sinusoidal inhibitory currents (assumed to start at phase 0°) drive proximal compartments C17 and C19, and 60° phase-advanced sinusoidal excitatory currents drive the distal dendrites C35 and C36. An asymmetric depolarizing ramp (duration 0.512 to 2.048 s, peaking at 75% of the duration) is applied either at the soma (as an injected current), midapical dendrites (C22–23), or distal dendrites (C35–36) as an excitatory synaptic input. Bottom right inset, Assumed GABAA, GABAB, and excitatory (glutamatergic) synaptic conductances at different compartments; small symbol, standard strength; large symbol, double strength (for details, see Materials and Methods and Table 1).
Specific membrane capacitance (Cm) of 1 μF/cm2, specific membrane resistance (Rm) of 30,000 Ωcm2 (Gm = 0.033 mS/cm2), and longitudinal resistivity Ri of 200 Ωcm are assumed. All compartments are assumed to have the same electronic length of 0.03λ, and Gl = Gm/(0.03)2 = 1111*Gm. Gmi/Ci = 0.033 ms−1 is the same for all compartments, but coupling strengths between compartments, or Gl/Ci (Eq. 1) may differ because of different Ci (Table 1). Division by Ci in Equation 1 converts all currents and conductances as per area since Ci = Cm*area of the compartment, where Cm is specific membrane conductance. Units used are Cm (microfarads), voltage (millivolts), time (milliseconds), conductance (millisiemens per square centimeter), and current density (microamperes per square centimeter).
Compartment parameters
The intrinsic currents INa (fast Na+ transient current) and IDR (delayed rectifier current), assumed for thick-circle compartments (14–23, 37, 38) in Figure 1, follow the description of Hoffman et al. (1997) as follows. INa = GNa * m3 h * (V − ENa), where GNa is maximal Na+ conductance, and ENa is 58 mV. Parameters m and h are originally defined by Hodgkin and Huxley, and referred to specifically as INa_m and INa_h, respectively. Each parameter (m or h) is defined by voltage-dependent rate constants α and β that determine the time constant τ and steady-state (subscript ss) value.
IDR= delayed rectifier K+ current = GK * n4 (V − EK),
For axonal compartments 18, 37, and 38, VNaD is replaced by VNaA (Table 1) and no IA is assumed, but GK is eightfold that of the other compartments. The voltage for half-maximal activation of INa at the axon and dendrites are named VNaA and VNaD, respectively, and will be referred to as the axonal (VNaA) and dendritic (VNaD) spike threshold.
IA (A current) and Ih (H current) follow the properties described by Leung and Peloquin (2006) as follows:
where fKA is a gain factor (channel density) that increases 1.57-fold for every 100 μm increase in depth from the soma toward the apical dendrites, as used by Hoffman et al. (1997), and 2.57-fold for every 100 μm increase from the soma toward the basal dendrite. Depths of the dendritic compartments and fKA are listed in Table 1. IA parameters are defined as follows:
The proximal compartments (Table 1) have higher voltage threshold for IA activation, as defined by the αmA(V) and βmA(V) parameters below (Hoffman et al., 1997):
The hyperpolarization-induced cation current Ih is described by a single parameter κ (Gasparini et al., 2004) such that
fh is a gain factor (channel density) that is 1 at the soma, increases twofold from 0 to 100 μm, and another twofold from 100 to 200 μm (Magee, 1998), as listed for the compartments in Table 1.
With an electrotonic length constant λ = 1375 μm, the soma compartment 17 is assumed to be a cylinder of 8 μm diameter and 51 μm length (0.03λ), and the unbranched apical dendritic shaft (compartments 19 and 20) of 5 μm diameter and 82 μm length is equivalent to 0.06λ. Two primary basal dendritic branches are assumed to be 3.2 μm diameter. The actual depths of the dendritic compartments with respect to the alvear surface (Table 1) are relevant for field potential generation and not for intracellular voltage simulations. These compartment depths assume that (1) bifurcation of the primary and secondary dendrites occurs every 80–100 μm and satisfies the Σd3/2 rule of Rall (1964), and (2) the two basal dendritic branches are angled at 30° from the vertical, while two apical dendritic branches angled at 15° come from the apical dendritic shaft represented by compartment 20. The two daughter apical branches are not shown in Figure 1 but are lumped into a single compartment 21, since the two compartments together behave electrotonically as a cable of uniform diameter with the mother branch (compartment 20). Similar assumptions are made for further bifurcations along the apical and basal dendritic branches.
INaP = GNaP * m' * (V − ENa) = persistent (non-inactivating) Na+ current, where GNaP = maximal persistent Na+ conductance (French et al., 1990, Magistretti and Alonso, 1999; Vervaeke et al., 2006).
An INaP slow inactivation (time constant, ∼6 s), as described by Magistretti and Alonso (1999), is not used since its inclusion does not affect the main results.
IKS = GKS * n′ (V − EK) = M-type K+ current is described according to Yamada et al. (1989), with dn'/dt = (n'ss − ns)/τn',
The time constants for IKS are calculated as 33 and 24 ms at holding voltages of −45 and −55 mV, respectively, assumed for 37°C. If a temperature factor Q10 of 5 is assumed, time constants in the present model are of similar magnitudes as those experimentally measured by Halliwell and Adams (1982) in CA1 pyramidal cells (92 and 47 ms at 30°C) and by Wang et al. (1998) in dissociated rat sympathetic ganglion (132 and 115 ms at room temperature).
Simulation of spiking.
Subthreshold MPOs without spikes are simulated by using axonal (VNaA) and dendritic spike threshold (VNaD) of greater than −32 mV, unless otherwise indicated. This is necessary for showing the properties of MPOs below the spike threshold. (GNa/Ci), (GKA/Ci), and (GDR/Ci) are 60, 24, and 15 ms−1, respectively. For the normal spiking condition, which simulates both MPOs and spike precession, the spike thresholds are set at VNaD = −34 and VNaA = −39, and Vnoise is added to VNaA, with Vnoise assumed to be Poisson-distributed with a mean of 4 mV and to change every millisecond. The Poisson-distributed Vnoise adds variability to spiking, giving it a low rate at rest (0.05 to 0.1 Hz at approximately −65 mV somatic membrane potential), and not necessarily firing at consecutive MPO peaks; (GDR/Ci) is increased to 30 ms−1 when Vnoise is added, to limit burst activity, but all other parameters remain the same.
Excitatory and inhibitory driving and current injection.
Step currents of amplitude SA and >2 s duration are injected into the soma. Injected somatic currents also include a ZAP function (Hutcheon and Yarom, 2000) with time t, of the form SA*sin(2πf(t)*t) over a time period T, where SA = amplitude, and the swept frequency increases from zero to a maximum of fm, f(t) = fm*t/T, for 0 <= t <= T; with fm = 20 or 40 Hz. The input current and soma voltage output are resampled every millisecond (for at least 4096 points), and the amplitude and phase of the complex impedance (voltage/current) are estimated from Fourier transform.
An asymmetric ramp-like excitation (Mehta et al., 2002; Harvey et al., 2009), of duration T and peaking at 0.75T from onset is assumed. This ramp is applied with peak conductance RpG at the distal dendrites (compartments 35 and 36), or at the middle dendrites (compartments 22 and 23), or as an injected current (peak RpI at the soma); the ramp amplitude is adjusted to yield <15 mV depolarization at the soma. Noise inputs are added to most simulations. For a distal ramp, Poisson-distributed impulses (mean interval, 5 ms) are applied to excite the distal dendrites with a peak conductance of ExG and single decay rate constant of 0.35 ms−1. For a proximal current ramp, current noise is added with Poisson distribution (with a mean interval of 5 ms), amplitude ExI and decay rate constant 0.35 ms−1. Typical noise magnitudes are ExG = 0.05 to 0.4, and ExI = 5 (for somatic ramp), or ExI = 0.2*SA (for ZAP function input).
A rhythmic inhibitory conductance Si * sin(2πf*t) is applied to the proximal compartments, and an excitatory conductance driving function Se * sin(2πf*t + φ) is applied to the dendrites (Fig. 1), where f = 5–12 Hz, Si is peak inhibitory conductance, and Se is peak excitatory conductance. The driving functions are sinusoidal, assuming increase and decrease from a tonic level. Proximal-inhibitory and distal-excitatory driving functions are consistent with models of generation of extracellular theta rhythm (Leung, 1984; Kamondi et al., 1998). The excitatory drive is assumed to be 60° (=φ) advanced of the inhibitory drive (Leung, 1984). This results in firing of pyramidal cells at ∼250° phase during baseline (see Fig. 5B), while maximal inhibition occurs at 0 or 360°. For display in the figures below, the rhythmic inhibitory driving is reversed in polarity [i.e., shown as −Si * sin(2πf*t)] to facilitate a direct comparison with the membrane potential response (since positive or high inhibition results in negative membrane potential response).
The MATLAB differential equation solver [ODE23 using Runge–Kutta (2,3) formula] uses variable time steps (<0.09 ms), and the resulting membrane voltages are linearly extrapolated and sampled at 1 ms interval. If spikes are not present, the MPO signal is the soma voltage after low-pass (<16 Hz) digital filtering that does not shift the phase. In the presence of spikes, defined as sharp transients above −30 mV in the wide band simulated soma voltage, the spike peaks are first identified. A MPO signal without spikes is derived from the 1 ms resampled voltage data, after each spike peak and two adjacent points (total of 3 ms duration) are removed and replaced by the average voltage before and after the removed points. The latter procedure removes the spike peaks but spike afterpotentials may remain. The peaks of the inhibitory driving function are used as the reference (zero phase) for measuring phase shift of the MPO or spike peaks. Power of the MPOs at the soma is estimated by Fourier transform during baseline (called base power) and for the whole duration of the depolarizing ramp (called place power); the same duration (typically 1024 ms) is used to construct the power spectra for baseline and place.
To study whether phase precession is disrupted by evoked inhibition and excitation of CA3 afferents, single-pulse stimulation (impulse function) of the associational and commissural afferents is assumed to activate excitatory conductance in the midapical dendrites (Fig. 1, inset; Table 1) with maximal conductance of 1.13 mS/cm2 at compartments 21, 22, 25, and 26, and 2.27 mS/cm2 at compartments 23 and 24. Single-pulse evoked GABAA and GABAB inhibitory conductances are described by the function (t − td) * exp(−t − td/τ); td = 2, τ = 33.3 ms for GABAA conductance function acting on proximal compartments (Fig. 1, inset) that peaks at 2.11 mS/cm2. For GABAB conductance, td = 20 and τ = 192 ms, and it peaks at 0.18 mS/cm2 in dendritic compartments with standard strength and 0.36 mS/cm2 in compartments with double strength; GABAB conductance strengths are listed in Table 1.
Results
Phase code of place location
To clearly visualize the dynamics of the subthreshold MPOs, simulations using a high spike threshold (VNaD = VNaA = −28) are first presented. During baseline, before an animal moves into the place field of the simulated neuron (Fig. 2A), combined excitatory and inhibitory rhythmic inputs result in MPO peaks of ∼290° phase with respect to the inhibitory conductance driving function. In the place field, a distal-dendritric excitatory ramp is assumed. Phase precession of the MPO peaks starts approximately one theta cycle after ramp onset and continues after the peak of the ramp (Fig. 2A). In addition, the amplitude of the MPOs gradually increases during the ramp (Fig. 2A). After the place field is traversed, rapid restoration of the baseline phase occurs. If the duration of running through the field is longer (Fig. 2B, 1024 ms; C, 2048 ms), the decrease of phase with time becomes more gradual (Table 2). When MPO phase is plotted as a fraction of the animal's distance in a one-dimensional place field, assuming constant speed, a similar phase precession function with distance is found for different running speeds (Fig. 2D). The phase decrease with distance is approximately linear with time for the initial 85% of the place field, and then the phase may increase slightly during the last 15% of the field (Fig. 2C,D, open arrow). For a 512 ms run duration, no phase increase is found; this relates to the observation that the phase decrease appears to persist for one more theta cycle, regardless of ramp duration, after the peak of the ramp at 75% of the field. The slope of the linear regression line changes slowly with duration (run speed)—it decreases 27% from 238° to 174° per place field while the duration changes fourfold (Fig. 2E).
Subthreshold model of place cell show phase precession that depends on animal's position in the place field. An animal is assumed to run through a fixed one-dimensional place field at three different velocities, traversing the field in 0.512 s (A), 1.024 s (B), and 2.048 s (C), respectively. The depolarizing ramp conductance (RpG = 1.5) is assumed to be the same function except for time expansion, excitatory synaptic noise (Poisson's distributed with 5 ms mean) is added to the distal dendrites, and rhythmic (9 Hz) distal-dendritic excitatory (Se = 0.1) and proximal inhibitory driving (Si = 0.2) are given. A, In each panel, reversed inhibitory (−inh) and excitatory (exc) conductance drive function, phase (filled circles, axis on right), depolarizing ramp in the place field (thick red line on the time axis), and the somatic membrane potential (unfiltered green trace, and low-pass filtered, red dashed trace) are plotted. The peaks of the “−inh” function are ±180°, and the phase of each peak of the MPOs is measured with reference to this function. B, C, Phase decrease within the place field is smaller for slower runs, and with a slight reversal of phase decline (arrow). D, Phase plotted with the place field normalized to unity shows similar phase precession across all speeds during the initial 85% of the place field and a small, reverse precession (arrow) during the rest of the place field. E, Magnitude of linear regression slope of phase in the initial 85% of place field (units of degrees/field) plotted with the duration of run through the place field.
Phase precession of subthreshold MPOs in preferred model (GNaP = 0.8; GKS = 8) with different input parameters, and depolarization ramp induced by distal dendritic excitation
Evoked inhibition only transiently disrupts phase precession
To test whether phase precession is disrupted by a strong afferent stimulation (Zugaro et al., 2005), the subthreshold model is subjected to an afferent stimulation pulse that generates an EPSP followed by an IPSPA mediated by GABAA receptors, with or without a GABAB receptor-mediated IPSPB. During the strong IPSPA, all synaptic responses are shunted and phase could not be estimated from the MPO peaks for approximately two theta cycles. In experimental situation, spikes were silent for 200–250 ms (Zugaro et al., 2005). As the IPSPA subsides, the simulated MPO amplitudes increase and phase precession continues, as if no IPSPA had occurred (Fig. 3A). If a prolonged IPSPB is evoked, the MPO phase is disrupted for approximately seven theta cycles (∼0.8 s) after the stimulus and then the phase and its precession resume at the normal magnitude and rate (Fig. 3B,D). Disruption of phase coding by an IPSPB has not been shown experimentally. The effect of the evoked inhibition does not depend on the timing of the inhibition in relation to the ramp. At all onset times, an evoked IPSPA disrupted the phase precession for approximately two theta cycles, and an evoked IPSPB disrupted the phase for approximately seven theta cycles.
Perturbation of the subthreshold model with evoked inhibition only transiently disrupts phase precession of the subthreshold membrane potential oscillations. Afferent pulse stimulation given at 600 ms (up arrow), after onset of the ramp (place field) at 512 ms, disrupts phase for ∼200 ms when GABAA inhibitory conductance is evoked (A), and for ∼800 ms when both GABAA and GABAB inhibitory conductances are evoked (B). C, IPSP with GABAA conductance alone (IPSP-A; red trace), or with GABAB conductance (purple dotted trace; IPSP-A+B); a 30 mV EPSP precedes the IPSP. D, Phase precession during the place field without afferent stimulation (“control”; two simulations with open circles), with evoked IPSP-A (two simulations with square symbols), or evoked IPSP-A+B (solid circles). Regression line for control data has slope of −70°/s, R2 = 0.8. GABAA conductance is assumed for the soma and proximal dendrites and GABAB conductance for the middendrites (Fig. 1, inset).
Phase precession of spike follows that of MPOs
When spiking occurs, they show a gradual phase precession in the place field, with each spike generally firing near an MPO peak (Fig. 4B), and thus, spikes and the filtered MPO peaks precess at similar rates (Fig. 4B,C; Table 3). When spikes are suppressed by setting a high spike threshold (Fig. 4A), subthreshold MPOs (sub MPO) still precess at a similar rate as the filtered MPOs during normal spiking (Fig. 4, compare A, B, and group data plotted in C). Despite the asymmetric depolarization ramp, the average spike rate within the field is symmetric (Fig. 4C, inset), as was observed experimentally (O'Keefe and Recce, 1993). As shown for subthreshold MPOs (Fig. 2D), MPOs and spike precession during normal spiking condition show a small dependence on run duration (Fig. 4D). Phase precession of spikes is generally similar to that of MPOs for run duration of 1–2 s, but it is more variable (Fig. 4D, Table 3) as expected from the addition of spike threshold noise (see Materials and Methods).
Spike precession resembles MPO peak precession in spiking model (Se = 0.2; Si = 0.4), during 512 ms baseline and 1024 ms place field (red line on time axis). A, From top to bottom traces, depolarizing ramp with distal apical dendritic noise, excitatory (exc) and reversed inhibitory (−inh) drive, phase at MPO peaks (filled blue circles) or spikes (filled red circles), continuous traces (bottom) of unfiltered (green) and filtered (0–16 Hz) somatic potentials. Spiking suppressed with high axon spike threshold (VNaA = VNaD = −34 mV). B, MPO traces and phase with normal spike thresholds (VNaA = −39 mV; VNaD = −34 mV), and single or double spikes fire near the peak of the MPOs. C, Phase plot, showing the average phase of the subthreshold MPO peaks (sub MPO) with high axon spike threshold (12 trials, X symbols linked by a dash line), and scatter plot of the phase of MPO (open circles) and spike peaks (solid square symbols) with normal axon spike threshold (12 trials). Note high and low spike threshold simulations yield a similar linear regression slope of the phase precession for MPO peaks (dotted line; −146°/s; R2 = 0.55) and spike peaks (thick line; −138°/s; R2 = 0.51) with the latter line displaced downward from the former. The inset shows the firing rate running average (in spikes/second, integrated every 0.2 s). D, Magnitude of linear regression slope of phase of the MPO peaks or spikes during the initial 75% of place field (units of degrees/field) for different run duration. E, Regression slopes for a fixed (1024 ms) duration run plotted versus different average firing rate within the place field simulated by varying spike thresholds. Thresholds from low to high firing rate are, respectively, VNaA = −36.5, −37.5, −39 (above with VNaD = −34), then VNaA = −39 and −40 (both with VNaD = −36); units are in millivolts.
Phase precession parameters simulated in model with Poisson distributed axon spike noise, for different rhythmic inputs, excitatory input with strength Se and inhibitory input with strength Si, and with or without Ih conductance
The relationship between spike rate and phase precession in the place field is investigated by varying spike thresholds at the axon and/or dendrites (VNaA and VNaD do not vary in time in the simulations). Varying spike thresholds changes the mean firing rate, but theta phase precession of the spikes and MPOs in the place field is found with all firing rates, except spike phase precession is minimal at a firing rate of <2/s (Fig. 4E). There is a small increase of precession per field with firing rate change from 3 to 13 Hz, but the simulated data are consistent with the presence of spike phase precession during runs with low and high firing rates (Huxter et al., 2003).
A small phase precession reversal after the peak of the depolarizing ramp is also found for the simulated spikes during long-duration runs (Figs. 4C, 5B,D, open arrow). However, experimentally, spike phase precession within the place field was reported to continue after spike firing decreased (O'Keefe and Recce, 1993; Huxter et al., 2003), and some cells even showed steeper phase precession near the end of the place field [Skaggs et al. (1996), their Fig. 7]. In the present model, continuous spike precession without reversal can be achieved if inhibitory driving is terminated, for example, at a time near the peak of the depolarizing ramp [Fig. 5C, asterisk (*)]. The latter results in earlier cessation of firing (Fig. 5, compare C with B; D, inset), but for the few spikes that fire after the ramp peak (Fig. 5C,D, solid arrow), they fire near the peak of a membrane potential that oscillates at ∼12 Hz. These ∼12 Hz oscillations are near the natural frequency of MPOs induced by the depolarization ramp (Fig. 5A); and the oscillations are there whether excitatory driving is present or not. A fast phase precession at a rate of 0.85°/ms (Fig. 5D, dashed line at 2000–2200 ms) is caused by the difference between the driving frequency (9 Hz) and the natural frequency.
Early termination of inhibitory driving results in intrinsic oscillations and accelerated spike phase precession. A, Membrane potential shows ∼13 Hz oscillations in response to 2048 ms duration depolarizing ramp alone, after a baseline of 512 ms. B, Depolarizing ramp with inhibitory and excitatory driving induce spike precession, which reverses (open arrow) at ramp repolarization. C, Early termination of inhibitory driving (*) at the peak of the depolarization ramp (time, 2048 ms) allows natural oscillations (∼12 Hz) to emerge immediately after the asterisk (*), with spikes arising from the MPO peak (solid arrow). D, Group data for spike phase versus time with “full” duration (filled square symbols) and “short” (early terminated) inhibitory driving (unfilled circle symbols). Spike firing rates for both conditions are shown in inset; 10 simulations for each condition. Base phase indicates average phase of MPO peaks during baseline. Linear phase precession of spikes occur for both conditions during the rising phase of the depolarization ramp (solid red line; −80°/s; R2 = 0.61). After the depolarization peak, reversed phase precession for full duration driving (open arrow) does not happen with short duration driving which results in a faster phase precession (solid arrow) after the ramp peak at >2048 ms (dotted red line; −852°/s; R2 = 0.44).
Rhythmic inhibition and depolarization result in power increase and precession of MPOs
As described above, when combined rhythmic excitation and inhibition are given during the depolarizing ramp, there is a gradual phase precession of the MPOs and spikes (Fig. 6A1). The amplitude of the MPOs increases in the place field, as shown by the power spectrum in the place field (place) compared with that during baseline (base) (Fig. 6A2). In the place field compared with baseline, the peak frequency of the MPOs is increased [Fig. 6A2, asterisk (*)], together with a clear second harmonic of the theta driving frequency (Fig. 6A2, arrow).
Inhibitory or excitatory rhythmic synaptic driving induces weak phase precession of the MPOs in the spiking model. Left column 1, Phase measures at the MPO peaks (filled circles) are plotted with the intracellular potential (lowest trace) and ramp/sinusoidal driving. Right column 2, Power spectrum of the MPOs in simulated place field (“place”; solid red trace), compared with that during baseline (“base”; dashed blue trace). A, MPO and spike precession are observed after combined excitatory and inhibitory driving, but much weaker after only excitatory (exc) driving (B), or only inhibitory (inh) driving (C). D, Same as A after combined excitatory and inhibitory driving but with Ih = 0; MPO power is increased, particularly during baseline. Column 1 shows an increase in power and frequency of MPOs in the place field in A and D (*), and an increase in second harmonic of theta (arrow) in A, C, and D.
When driven by rhythmic excitation only, phase precession of MPOs during the depolarizing ramp is small (Fig. 6B1), compared with combined excitatory and inhibitory driving (Fig. 6A, Table 3). Rhythmic excitation results in MPOs with decreased power/amplitude during the ramp depolarization (Fig. 6B2) compared with baseline. The latter is caused by a decrease in electromotive force of excitation (reversal potential, 0 mV; baseline resting potential, −66 mV). Few spikes fire within the place field, and there is no significant spike phase precession when multiple simulations are grouped, or even when spike thresholds are decreased to allow for a higher firing rate.
When driven only by rhythmic inhibition, phase precession of MPOs and spikes in the place field is modest but consistently observed (Fig. 6C1, Table 3). The MPO amplitude and power increase during the depolarizing ramp (Fig. 6C2) because of an increase in the electromotive force of inhibition (reversal potential assumed for Cl− is −72 mV) and intrinsic resonance (see below).
The phase precession of MPOs decreases slightly when the rhythmic modulation strength increases (Table 2). Driving frequency within 5–11 Hz results in similar theta phase precession in the place field, while driving at 13–15 Hz gives maximal MPO amplitudes but relatively lower phase precession during the second half of the place field (data not shown). The 9 Hz data presented are representative of the 5–11 Hz frequency range.
Nature of the depolarizing ramp and rhythmic excitatory driving
With a distal dendritic depolarizing ramp, which has been assumed for the above simulations, an increase in MPO frequency and power does not occur with rhythmic excitatory driving alone (Figs. 6B, 7A). Since both rhythmic excitatory driving and ramp act on the distal dendrites (compartments 35 and 36; Fig. 7A1), and the excitatory electromotive force decreases with depolarization, a decreased response to rhythmic excitatory driving (Fig. 7A2,A3) during the ramp is observed. If the ramp inputs are shifted to excite the midapical dendrites (compartments 22 and 23), while keeping the average ramp depolarization at the soma the same, decreased response to rhythmic distal-dendritic driving in the place field is still apparent at the midapical compartments, but it is compensated by proximal resonant response such that the somatic MPOs are mainly unchanged in amplitude in the place field (Fig. 7B). More robust increase in somatic MPO amplitude in the place field occurs if the ramp is moved proximally, either becoming a depolarizing current at the soma (Fig. 7C) or a basal-dendritic excitation (data not shown). Voltage-dependent resonance in the place field is evidenced by a large increase in somatic MPO amplitude [Fig. 7C2, asterisk (*)] and an increase in MPO frequency and power (Fig. 7C3). The somatic ramp alone, without rhythmic driving, induced intrinsic oscillations of ∼12 Hz (Fig. 7C2).
Decrease of excitatory driving occurs during a depolarizing ramp. The location of depolarizing ramp changes from distal-dendritic excitation at compartments 35–36 (A), middendritic excitation at compartments 22–23 (B), to depolarizing current at the soma (C). Rhythmic excitatory driving remains at the distal dendrites throughout. Subthreshold MPOs are simulated (VNaA = VNaD = −28) for clarity. Column 1 shows voltage responses in different compartments of the model, 17 (soma), 23 (midapical dendrite), 29 (middistal apical dendrite), and 36 (most distal apical dendrite); red line on time axis (1024 ms duration) indicates active ramp function within the place field. Column 2, Phase and voltage at the soma, in response to rhythmic driving and ramp (solid trace), or ramp only (light dotted trace). Column 3, Power spectra of the soma voltage in the simulated place field (place) or before the place field (base). Occlusion of rhythmic excitatory response is strongest when both rhythmic and ramp excitation are acting on the same compartments (A) than when they are apart (B) and farther apart (C). A robust increase in response amplitude (* in C2 and C3) is found during a somatic ramp, with an increase in response frequency, including a second-harmonic in the place field (C3, arrow), compared with the baseline. Response to the somatic ramp only is reproduced below the time axis in C2, to highlight the intrinsic oscillations (∼12 Hz at *) induced by the ramp only.
High resonance models show MPO power increase and large phase precession
Since INaP and IKS have been implicated in theta-frequency resonance and MPOs (Introduction), different peak conductances of INaP and IKS, GNaP and GKS, respectively, are simulated in subthreshold models. With the parameters presented for the normal GNaP model above (GNaP = 0.8; GKS = 8), step current depolarization of the membrane to approximately −59 mV evoked a damped oscillation (Fig. 8A1), characterized by a relatively sharp MPO power peak at 12 Hz (Fig. 8A2). The half-peak bandwidth (HBW) of the MPO power spectrum measures 7.5 Hz, within the range of HBW determined experimentally by step or ZAP currents in CA1 pyramidal cells in vitro (Leung and Yim, 1991; Leung and Yu, 1998). Applying a swept-frequency sinusoidal current (ZAP function from 0 to 40 Hz) at a mean membrane potential of −59 mV also yields a maximal impedance magnitude at ∼12 Hz (Fig. 8A). Both magnitude and phase shift of the impedance are larger at −59 mV compared with −67 mV resting membrane potential (Fig. 8A2, rest).
Subthreshold membrane potentials are induced by a step depolarizing current, or with current swept from 1 to 40 Hz (Izap, green trace) simulated in models with different conductances for persistent Na+ current (GNaP) and slowly inactivating K+ current (GKS). For all simulations, resting membrane potential ranges from −65 to −66.5 mV. Column 1, The soma (Vstep; red traces) is depolarized ∼8 mV by a step current (onset at red dot; horizontal dotted red line is −60 mV), or induced to spike (black trace) by a 5 ms current (red line underneath spike trace). In another simulation, during the step current depolarization, a zap current (Izap, shown only in A) is applied, giving soma potential Vzap (blue traces). Poisson's noise ExI = 0.2*SA is added to the step and ZAP currents. Column 2, Logarithmic power spectrum of Vstep or impedance (Imp) spectrum estimated by Izap. A, Model with normal GNaP = 0.8, GKS = 8, in units of millisiemens per square centimeter, shows damped oscillation with current step, spike with depolarizing afterpotential (DAP), and a spectral peak at ∼12 Hz in both power and impedance; zap current amplitude SA = 2 μA/cm2. Imp and phase function estimated in column 2 during rest (−65.5 mV) and steady depolarization (−59 mV). B, Model with one-half GNaP (= 0.4), GKS = 8 shows poorly rhythmic Vstep, spike with small DAP, and broad power and impedance peaks; SA = 9. C, Model with zero GNaP (also GKS = 0) shows no detectable rhythm, no spike DAP, power or impedance peaks, SA = 5. D, Model of high GKS (= 24) and GNaP = 0.8 shows no theta-frequency rhythms, spike with large afterhyperpolarization (AHP), and broad power/impedance peaks at 20–25 Hz, SA = 18.
In a model with half normal GNaP (GNaP = 0.4; GKS = 8; Fig. 8B), step depolarization or ZAP current driving results in less rhythmic MPOs and low resonance, with a broad power peak at ∼16 Hz. A model with zero GNaP and GKS shows no rhythmic response and no clear power peaks, with step or ZAP currents (Fig. 8C, zero GNaP model). A high GKS model (GKS = 24; GNaP = 0.8; Fig. 8D) results in a broad response peak at 20–25 Hz, using a step or ZAP current.
Somatic current-induced depolarizing ramp and proximal rhythmic inhibitory driving are used to study phase precession in the above models. Only subthreshold MPOs are illustrated. In a high-resonance model, strong phase precession occurs during the ramp, and a clear increase in MPO power and frequency is observed in the place field compared with baseline (Fig. 9A). Reducing GNaP decreases resonance and the slope of phase precession in the place field (Fig. 9A–C, dashed linear regression line; Table 4, Inhibitory driving alone). However, MPO power in the place field is still higher than baseline in the zero GNaP model (Fig. 9C2). The model with high GKS shows a moderate phase precession and only a small MPO power increase in the place field (Fig. 9D). In the same models, driving by a rhythmic excitatory input generally yields a smaller phase precession than inhibitory drive, but a similar result that phase precession in the place field decreases with resonance of the model is observed (Table 4, Excitatory driving alone). No phase precession is found with purely excitatory diving in the zero GNaP (and zero GKS) model.
Models with strong resonance show stronger phase precession of the subthreshold MPOs. Somatic voltage (column 1) and power spectra in place field and baseline (column 2) are shown. Rhythmic inhibitory drive is given with a depolarizing somatic ramp current; resting membrane potential and average ramp depolarization are kept within 1 mV among the different models (Table 4). The four versions of the model in this figure corresponds to A–D of Figure 8, with normal GNaP, high-resonance model showing an increase in MPO amplitude and frequency, and MPO phase precession (A), and models with little or no theta-frequency resonance showing no increase in MPO amplitude or frequency during the ramp, and a modest phase precession in some cases (B–D).
Phase precession of subthreshold MPOs during a somatic injected depolarizing current ramp of 1024 ms, driven by either rhythmic excitatory or rhythmic inhibitory input at 9 Hz
When normal spike thresholds are assumed in the above models, a brief 5 ms current induces a spike followed by afterpotentials of different shapes in different models (Fig. 8, column 1, spike trace). The normal GNaP model (Fig. 8A1) has the largest depolarizing afterpotential (DAP), while the high GKS model (Fig. 8D1) has the largest afterhyperpolarization (AHP) of medium duration.
In separate simulations, it is shown that reduction or abolition of Ih has no effect on the phase precession of MPOs (Fig. 6D, Table 3), when a similar resting membrane potential is kept. Reduction of Ih increases input resistance and enhances dendritic inputs preferentially, giving larger MPO amplitudes, more so during baseline than in the place field (Fig. 6, compare D with normal Ih responses in A). Similarly, 50% reduction of IA or IK increases input resistance but does not affect phase precession. The small phase precession found with inhibitory driving in the zero GNaP model (Fig. 9C) can be partly attributed to voltage-dependent conductance of the Na+ transient current. A small increase in MPO phase precession rate is shown with increasing spike rate, or decreasing spike threshold (Fig. 4E).
Phase precession explained by passive and active (voltage-dependent) processes
Phase precession in the place field is attributed to both passive and active mechanisms. The passive mechanism occurs as a change in balance between excitatory and inhibitory drives (e.g., simulated in Figs. 2 and 6). During the baseline, the phase in response to a pure inhibitory drive (with respect to the peak of the inhibitory conductance) is ∼180° (hyperpolarization is a negative response; Fig. 6C), and the phase in response to a pure excitatory drive is ∼300° (Fig. 6B). The resultant vector (R1), at resting membrane potential of −67 mV, is approximated by the vector sum of the inhibitory response (vector I1) and the excitatory response (vector E1; Fig. 10A). In the place field, where the depolarizing ramp is assumed to peak at −52.7 mV, the excitatory response (E2) decreases by approximately one-half, because of a decrease in electromotive force and nonlinear shunting (Fig. 7A2), while the inhibitory response (I2) increases approximately sevenfold because of an increase in electromotive force. The resultant vector R2 in the place field thus shifts toward I2, a decrease in phase by 87° from R1.
Explanation of phase precession by passive and active mechanisms. Phase plane diagram of amplitude (in millivolts) and phase change of the inhibitory response (I), excitatory response (E), and their resultant (R) vector, during baseline (resting, −67 mV; vector subscript 1) and peak depolarization in place field (−52.7 mV; vector subscript 2). A, Passive condition, excitatory driving during baseline (E1) is assumed to be 300°, while I1 response is 180°. Change in electromotive force in the place field increases inhibitory response and decreases excitatory response, resulting in 87° phase shift of R2 from R1. B, Active condition with a voltage-dependent impedance and constant current driving (without change in electromotive force); either excitatory or inhibitory response is increased in amplitude and phase shifted by 44°. The bottom inset illustrates the phase plotted with frequency at rest (baseline) and in the place field, measured using ZAP functions (Fig. 8A). C, Combined condition approximated by linear summation of passive and active responses, resultant vector phase change is 132°.
Phase and amplitude change also occurs because of active resonance, contributed by a voltage-dependent intrinsic mechanism. The steady-state phase lag in response to a constant 9 Hz driving current is −6.3° during baseline (−67 mV), and −50.7° at the peak of the depolarization ramp (Fig. 10B, inset), as determined using swept-frequency current (Fig. 8A). This voltage-dependent phase shift occurs with an excitatory drive, resulting in a negative 44° phase shift of E2 with respect to E1 (Fig. 10B). Similarly, I2 is shifted −44° from I1 (Fig. 10B), and the resultant vector (sum of E and I vectors) is also −44° shifted from that during baseline (data not shown). The active resonance contributes an additional 44° to the passive condition (which gives 87°), resulting in a 132° shift of R2 with respect to R1 (Fig. 10C). Phase change because of self-sustaining intrinsic oscillations (Figs. 5A,C, 7C2) has not been included in Figure 10. The latter phase precession is caused by having an intrinsic oscillation of higher frequency than the reference theta signal (9 Hz assumed above).
Single-compartment model of phase precession
To show that the main properties of phase precession of a resonant neuron are not related to an interaction between soma and dendrites and the multiple parameters of the model, simulations are made in a one-compartment model. The latter is the soma compartment of the 38-compartment model, which now also receives a rhythmic excitatory synaptic drive, in addition to the rhythmic inhibitory drive and a somatic current ramp in the following simulation.
The single-compartment model with no GNaP conductance (GNaP = GKS = 0) shows no subthreshold resonance when driven by a ZAP function at a steady-state resting membrane potential of −58.7 mV (Fig. 11A3,B). Adding GNaP = 0.1 (and GKS = 1) conductance bestows a peak response at ∼11 Hz and half-peak bandwidth of 16 Hz (Fig. 11B). When combined rhythmic excitatory and inhibitory inputs are given to the GNaP = 0.1 model with a high spike threshold (VNaD = −32) assumed, gradual phase precession appears soon after the onset of an excitatory ramp (Fig. 11C), and the phase decreases during the depolarizing phase of the ramp. The phase precession of the MPOs remains similar when a normal spike threshold (VNaD = −35) is assumed (Fig. 11C2). The phase decrease within the place field occurs with different run durations, as shown in the plot of phase versus normalized field duration (Fig. 11D). Similar to the 38-compartment model, the magnitude of the average phase precession slope for the MPO peaks is only weakly dependent on duration (Fig. 11E). Spikes ride on top of the MPO peaks, and the linearly fitted phase precession slope of spikes is larger than that for MPO peaks (Fig. 11E; also apparent in C2).
One-compartment model shows run duration-independent phase precession in the place field. A, Voltage response (V) to a swept-frequency sinusoidal current (Izap in row 1), in resonant model with normal GNaP = 0.1 and GKS = 1 in row 2, and model with GNaP = GKS = 0 (row 3); in each case, steady-state depolarizing current maintains depolarization near −60 mV (dotted horizontal line). B, Impedance (voltage/current) as a function of frequency showing resonance with normal GNaP model when depolarized to −59.8 mV (half-peak bandwidth, ∼16 Hz), and no resonance in zero GNaP model. C, Representative time course of phase precession during rhythmic driving by combined rhythmic excitatory (Se = 0.005) and rhythmic inhibitory (Si = 0.04) inputs. An excitatory ramp (RpG = 0.05) depolarizes the membrane to an average of −59.8 mV during the ramp is assumed during the place field, and either subthreshold (VNaD= −32; C1) or normal spike threshold (VNaD = −35 mV, Vnoise = 0, GDR/C1 = 15 ms−1, in C2) is used. Results simulated using a normal spike threshold are summarized in D and E. D, Phase plotted as fraction of the place field shows similar phase precession for runs of three different durations (0.5, 1, and 2.048 s), but slightly lower than that of the 38-compartment model (38-C; symbols with connected dashed line). E, Magnitude of linear regression slope of MPO peaks and spikes during the depolarizing phase of the ramp (°/field) plotted with the duration of run through the place field; all slopes have linear fit of R2 > 0.8. The slope of phase precession is larger for spikes than for MPO peaks, and both slopes show <25% change in magnitude with fourfold change in run duration.
The one-compartment model (Fig. 11C) shows a large increase in MPO amplitude with ramp depolarization in the place field, similar to that during a somatic ramp in the 38-compartment model (Fig. 7C). Impedance spectrum of a single compartment (Fig. 11B) is similar to that of the 38-compartment model (Fig. 8A). However, the single compartment model shows a smaller rate of MPO phase precession, and a smaller total phase decrease in the place field compared with the 38-compartment model (Fig. 11D).
Discussion
Space coding and theta-frequency oscillation
This report illustrates how oscillations can be used to generate neural coding of absolute space (Figs. 2D, 4D, 11D). A spatial signal is coded by a depolarization ramp, and the transduction of this depolarization into a phase code depends on an external oscillator serving as a phase reference as well as a driving input and intrinsic membrane properties. A phase code for MPO is directly related to a phase code for spiking, since spikes tend to fire at the MPO peaks, in support of the notion that precession of MPOs is the cause of spike precession (Harvey et al., 2009).
Like most reliable coding, the present phase code has a considerable safety margin. Phase precession of the MPO peaks is only briefly affected by single-pulse afferent-evoked excitation and inhibition (Fig. 3). It can occur with rhythmic excitatory or/and inhibitory driving of different modulation depths. Phase precession results from different mechanisms induced by the depolarization ramp. The first is a passive phase change caused by a shift in the relative strengths of the excitatory and inhibitory responses (Fig. 10A). The second is an active theta-frequency resonance that mainly depends on the interaction between slow voltage-dependent Na+ and K+ currents (Fig. 10B). The third is a nonlinear interaction of driven and intrinsic oscillations; intrinsic oscillations can be induced by a moderate depolarization of approximately −60 mV alone (Figs. 5A, 7C2). Other than generating phase precession, the active voltage-dependent mechanisms are critical to MPO frequency increase and a second harmonic in the place field (Figs. 6, 9A).
Phase precession of the order of 180°, consistent with the spike phase precession after a single run through the place field (Schmidt et al., 2009), are simulated by passive and resonance mechanisms in the present report. Larger phase precession of 180–360° (O'Keefe and Recce, 1993; Skaggs et al., 1996) may require activation of intrinsic oscillations (Figs. 5C,D, 7C2). A small reversal of phase precession is found at the end of a long-duration simulated run through the place field (Figs. 2D, 5B; but see Figs. 2A, 7A2). Whether this small phase reversal of the MPOs occurs in vivo is not known, but spikes were reported to continue their phase precession despite decreasing firing rate (Huxter et al., 2003). Factors such as spike accommodation will suppress spiking during ramp repolarization (Harris et al., 2002), and cessation of inhibitory driving and release of intrinsic oscillations will maintain continuous spike phase precession (Fig. 5C).
Either excitatory or inhibitory rhythmic driving can result in phase precession in the place field during a depolarizing ramp. With a distal dendritic depolarization ramp, a rhythmic inhibitory drive typically gives stronger phase precession than a rhythmic excitatory drive. This is explained by a depolarization-induced increase in driving current for an inhibitory, but not an excitatory drive. An increase in rhythmic driving increases peak depolarization and gives stronger phase precession via a voltage-dependent mechanism.
Two membrane currents—INaP and IKS—particularly contribute to intrinsic phase precession. Interaction between INaP and IKS has been shown to result in oscillatory response in other studies (Gutfreund et al., 1995; White et al., 1995; Golomb et al., 2006; Prescott and Sejnowski, 2008). In CA1 pyramidal cells in vitro, theta-frequency resonance at a depolarized voltage depended on both INaP and IM (Hu et al., 2002, 2007), and low-threshold subthreshold theta-frequency MPOs were suppressed when INaP was blocked (Leung and Yim, 1991). INaP contributes to the depolarizing afterpotential (Yue et al., 2005), as shown by the simulated action potential (Fig. 8), and the present study suggests that neurons with different INaP to IKS ratio may show different phase precession rate (Fig. 9). The present report further predicts that phase precession of theta-frequency MPOs will be reduced when INaP or IM is blocked. Blockade of Ih in CA1 pyramidal cells did not affect MPOs (Leung and Yim, 1991) or proximal resonance (Hu et al., 2002), and the present model suggests that Ih blockade increases MPO amplitudes but does not affect phase precession (Fig. 6D). Nolan et al. (2004) reported larger theta field potentials in Ih-knock-out mice, but theta phase precession was not studied.
Properties of the external inputs
Two types of external inputs—a depolarization ramp and a rhythmic synaptic input—are assumed. The simulations here do not predict a location for the depolarizing ramp, but they elucidate the consequences of proximal versus distal ramp depolarization. A distal dendritic ramp would suppress rhythmic excitation and increase passive phase precession during combined excitatory–inhibitory driving, since a decrease of excitatory strength will shift the resultant phase more toward the inhibitory vector (Fig. 10A). However, for a fixed current, proximal rather than distal ramp causes a larger depolarization that activates stronger active currents, including INaP and IM, thus giving stronger voltage-dependent phase precession (Fig. 7). As distinct from models in which a depolarizing ramp serves to decrease spike latency (Magee, 2001; Harris et al., 2002; Mehta et al., 2002), the ramp in the present model decreases the period of the MPOs. Since the entorhinal cortex codes for spatial information (Hafting et al., 2005), a distal apical dendritic ramp excitation of CA1 from the entorhinal cortex (Witter, 2007) appears to be more realistic for place coding. Furthermore, other studies suggest that place coding does not require the integrity of CA3 (Brun et al., 2002, 2008). However, during recall, CA3 may be necessary (Wallenstein and Hasselmo, 1997), and CA3 pyramidal cells may provide a ramp excitation of CA1 pyramidal cells at the midapical dendritic synapses (Amaral and Lavenex, 2007; Witter, 2007).
In studies in vivo, injection of current at the soma did not induce MPOs in CA1 neurons of awake mice (Harvey et al., 2009) or affect the phase of MPOs of CA1 neurons during a theta rhythm (of ∼4 Hz) in urethane-anesthetized rats (Bland et al., 2002). It is possible that the conditions during which currents were injected are different from those present within the place field [e.g., lack of IM modulation (below) or nonoptimal theta frequency].
Comparison with previous models
Phase precession, frequency, and power increase of the subthreshold MPOs were intracellular signatures of place cells (Harvey et al., 2009). The present report emphasizes that an intrinsic voltage-dependent mechanism contributes to these signatures. Two recent reports suggest that network models can also account for the intracellular MPO results, as shown in a CA3 network model (Romani et al., 2010) and a nonlinear recurrent network of CA1 neurons (Jayet Bray et al., 2010). Given the different types of neurons that demonstrate theta activity (Bland and Colom, 1993; Buzsáki, 2002) and phase precession (Moser et al., 2008) in the hippocampal formation, including neurons in the dentate gyrus, CA3, and entorhinal cortex, it is likely that more than one mechanism contributes to phase precession. The strength of the present model lies in its close ties with physiological data, and the ability to experimentally test the model by manipulating intrinsic membrane properties of a single CA1 cell.
Interference among two or more independent frequency sources can result in phase precession in hippocampal place cells (Burgess and O'Keefe, 1996; Lengyel et al., 2003) and entorhinal grid cells (Burgess et al., 2007; Hasselmo, 2008). In particular, phase precession in place cells is proposed to result from interference of independent somatic and dendritic oscillations of different frequencies. In the present model, an independent somatic oscillation may be generated by an intrinsic voltage-dependent mechanism involving slowly inactivating, proximal Na+ and K+ currents that are activated by a small depolarization (Leung and Yim, 1991; Hu et al., 2002, 2007). No dendritic oscillations have been included in the present model because they require large depolarization (greater than −50 mV) to activate in vitro (Leung and Yim, 1991) or in vivo (Kamondi et al., 1998). A crucial difference between the present model and previous linear interference models is that intrinsic resonance/oscillations are voltage dependent and not manifested outside of the place field. Thus, there are no off-field oscillations as in a linear interference model (Burgess and O'Keefe, 1996). The contribution of intrinsic resonance/oscillations to MPO phase precession is not included in another somatic-dendritic interaction model (Losonczy et al., 2010).
Hippocampal activation may be needed for optimal coding of space through phase precession. A movement-related theta rhythm indicates hippocampal activation by the medial septum and brainstem (Vanderwolf, 1988; Lee et al., 1994; Leung, 1998; Bland, 2000; Buzsáki, 2002), which involves cholinergic, GABAergic, and glutamatergic afferents (Amaral and Lavenex, 2007). Septal inputs may increase excitability and spiking, as well as modulate subsequent response through intracellular signaling (Krnjević, 1993; Leung and Peloquin, 2010). In particular, muscarinic cholinergic modulation may decrease various K+ conductances in the hippocampus (Krnjević, 1993), including IM, and increase theta-frequency resonance. It is known that place field firing was disrupted by blockade of muscarinic cholinergic receptors (Brazhnik et al., 2004).
In conclusion, it is shown that rhythmic synaptic driving, with the appropriate intrinsic membrane properties, results in robust phase precession of the subthreshold MPOs. An increase in the amplitude (power) and frequency of the subthreshold MPOs in the place field is simulated.
Footnotes
This work was supported by Canadian Institutes of Health Research Grant 15685 and the Natural Sciences and Engineering Research Council of Canada. I thank Gyuri Buzsáki, Brian Bland, and Min-Ching Kuo for comments on this study, and Terrence Sejnowski for providing a stimulating environment for starting the simulations.
- Correspondence should be addressed to Dr. L. Stan Leung, Department of Physiology and Pharmacology, University of Western Ontario, London, ON N6A5C1, Canada. sleung{at}uwo.ca