Temporal Frequency of Subthreshold Oscillations Scales with Entorhinal Grid Cell Field Spacing

See allHide authors and affiliations

Science  23 Mar 2007:
Vol. 315, Issue 5819, pp. 1719-1722
DOI: 10.1126/science.1139207


Grid cells in layer II of rat entorhinal cortex fire to spatial locations in a repeating hexagonal grid, with smaller spacing between grid fields for neurons inmore dorsal anatomical locations. Data from in vitro whole-cell patch recordings showed differences in frequency of subthreshold membrane potential oscillations in entorhinal neurons that correspond to different positions along the dorsal-to-ventral axis, supporting a model of physiological mechanisms for grid cell responses.

The entorhinal cortex plays an important role in encoding of spatial information (13) and episodic memory (4). Many layer II neurons of rat entorhinal cortex are grid cells, firing when the rat is in an array of spatial locations forming a hexagonal grid within the environment (57). The spacing of firing fields in the grid varies with anatomical position of the cell along the dorsal-ventral axis of entorhinal cortex, as measured by distance from the postrhinal border (5). Neurons closer to the dorsal border of entorhinal cortex have shorter distances between firing fields. Computational models explicitly predict that differences in grid field spacing should correspond to differences in intrinsic frequencies of neurons along the dorsal-to-ventral axis (3, 8). This could provide systematic variation in the gain of a movement-speed signal for path integration (2, 3, 9).

Subthreshold membrane potential oscillations in entorhinal cortical stellate cells (10) arise from a single-cell mechanism involving voltage-sensitive currents (1113) and could contribute to network dynamics (14). We recorded subthreshold oscillations from 57 stellate cells in layer II of medial entorhinal cortex (fig. S1) in slices from different anatomical positions along the dorsal-ventral axis by using whole-cell patch clamp techniques (15). The position of each horizontal slice was measured individually relative to the dorsal surface of the brain (Fig. 1A).

Fig. 1.

Higher frequency of subthreshold oscillations in dorsal versus ventral entorhinal cortex. (A) Top, dorsal view of brain after slicing. Bottom, sagittal view showing anatomical location of horizontal slices A to H. (B) Mean frequency of subthreshold oscillations for neurons from dorsal slices (3.8 to 4.9 mm bregma) and ventral slices (4.9 to 7.1 mm bregma) at –50 mV and –45 mV. Error bars indicate SEM. (C) Subthreshold oscillation frequency plotted versus anatomical distance from dorsal surface. (D) Examples of subthreshold oscillations at –50 mV in dorsal regions (left) and ventral regions (right). (E) Corresponding autocorrelations used to measure frequency in (D).

Stellate cells in dorsal entorhinal cortex showed higher temporal frequencies (f) of subthreshold membrane potential oscillations compared with lower frequencies in cells from more ventral entorhinal slices (Fig. 1B). Dorsal cells (n = 30) were defined as cells recorded in slices taken between 3.8 mm [the border with postrhinal cortex (16)] and 4.9 mm from the dorsal surface of the brain. Ventral cells (n = 27) were defined as cells recorded in slices between 4.9 and 7.1 mm from the dorsal surface. Because the group means of the frequency of subthreshold oscillations recorded from these populations (Fig. 1B) can depend upon the mean membrane potential voltage, we performed this analysis separately for data gathered at different approximate holding membrane potentials of –50 mV and –45 mV. The mean frequency in dorsal cells was significantly higher than the mean frequency in ventral cells for measurements at about –50 mV (dorsal mean f ± SEM = 6.42 ± 0.40 Hz and n = 28; ventral f = 4.23 ± 0.32 Hz and n = 25; P < 0.001) and –45 mV (dorsal mean f ±SEM = 7.18 ± 0.50 Hz and n = 14; ventral f = 4.88± 0.59 Hz and n = 7; P < 0.01). Frequencies were determined by an automated autocorrelation analysis algorithm (15) and were used throughout unless otherwise noted. The significant difference between mean frequencies was replicated with analyses using the peak in the power spectra for data at –50 mV (dorsal f = 4.86 ± 0.37 Hz and n = 28; ventral f = 3.44 ± 0.25 Hz and n = 25; P < 0.01) and –45 mV (dorsal f = 5.54 ± 0.49 Hz and n = 14; ventral f = 3.76 ± 0.31 Hz and n = 8; P < 0.01) (fig. S2). The resting membrane potential, firing threshold, resistance, and age [dorsal mean age ± SEM = 19.8 ± 0.4 days and ventral age = 19.1 ± 0.4 days, P value not significant (NS)] did not contribute to the difference in the frequency of oscillations observed between neurons from dorsal versus ventral slices (15).

Entorhinal neurons showed a systematic difference in subthreshold oscillation frequency when plotted for different locations along the dorsal-ventral axis (Fig. 1C) for data at –50 mV (r = 0.48, slope = –1.26) and at –45 mV (r = 0.60, slope = –2.48). This pattern resembles the difference in spatial periodicity of grid cells previously recorded from layer II at different positions along the dorsal-ventral-axis in awake, behaving animals (5, 7). Examples of subthreshold oscillations recorded from individual entorhinal neurons illustrate the difference in dorsal versus ventral frequency (Fig 1D). 500-ms segments of the autocorrelations computed for the same cells in Fig. 1D demonstrate the difference in peak-to-peak wavelength (Fig. 1E).

The frequency of subthreshold oscillations has been shown to correlate with the peak frequency of membrane potential resonance at rest (–60 to –64 mV) in entorhinal neurons (17). We evaluated the resonant frequency of neurons by delivering a 20-s-long impedance amplitude profile (ZAP) stimulus and measuring the input frequency that caused the largest amplitude depolarization of membrane potential. Individual examples of the response to the ZAP stimulus are shown for four neurons from dorsal entorhinal slices (Fig. 2A) and four neurons from more-ventral slices. The resonant frequency of stellate cells was significantly higher in dorsal cells compared with that in ventral cells (dorsal frequency = 6.13 ± 0.41 Hz and n = 18; ventral frequency = 4.45 ± 0.40 Hz and n = 14; P < 0.01) (Fig. 2B). There was no significant difference in resting potential between the populations used for this analysis (dorsal potential = –60.12 ± 0.26 mV and n = 18; ventral potential = –60.11 ± 0.26 mV and n = 14; P NS). Resonant frequency was systematically higher in dorsal regions when plotted in relation to anatomical location (r = 0.51, slope = –1.18) (Fig. 2C).

Fig. 2.

Differences in resonance properties of neurons in dorsal versus ventral entorhinal cortex. (A) Examples of resonance responses of neurons in dorsal (blue) and ventral entorhinal cortex (red) in response to ZAP stimulus (black). (B) Histogram of mean peak resonance frequency for neurons in dorsal versus ventral entorhinal cortex. (C) Resonance frequency plotted versus anatomical distance from the dorsal surface.

To analyze the relationship of subthreshold oscillations to other intrinsic properties, we measured the change in membrane potential during a hyperpolarizing current injection that usually caused a slow depolarizing shift (a “sag” of the membrane potential). The sag was fit with a dual exponential equation starting just after the trough of the sag and ending near the end of the current injection at steady state potential (Fig. 3, A and B). The faster of the two time constants (τ1) of the sag was measured for hyperpolarizing current steps that ended at steady-state membrane voltages of –69.9 to –65 mV (V1) and at steady-state membrane voltages of –64.9 to –60 (V2). The time constant (τ1) of the sag increased with depth on the dorsal-ventral axis (at V1, r = 0.55 and slope = 9.67; at V2, r = 0.61 and slope = 8.23) (Fig. 3C), with significantly faster values of τ1 in more-dorsal portions compared with more-ventral portions at V1 (dorsal τ1 = 23.37 ± 1.24 ms and n = 17; ventral τ1 = 35.52 ± 1.81 ms, n = 14; P < 0.01) and at V2 (dorsal τ1 = 22.37 ± 1.60 ms and n = 11; ventral τ1 = 34.71 ± 2.68 ms and n = 17; P <0.01). The τ1 of the sag correlated with the frequency of subthreshold oscillations (at V1, r = 0.63 and slope = –4.32; at V2, r = 0.65 and slope = –4.10) (Fig. 3D). Voltage-clamp studies are necessary to confirm underlying currents. A potential candidate is the h current, which has been shown to underlie sag in stellate cells (11, 12).

Fig. 3.

Relationship of subthreshold membrane potential oscillations to sag potential. For single neurons in dorsal (A) and ventral entorhinal cortex (B), traces on left show membrane potential responses to multiple current injection levels. Single traces (right) show sag potentials at higher resolution. (C) τ1 of sag potential (ms) plotted versus anatomical location. (D) τ1 of sag potential plotted versus subthreshold oscillation frequency.

The differences in temporal frequency shown in this study correspond to differences in spatial periodicity of unit firing observed with extracellular recording in awake, behaving animals (57). Neurons at more-dorsal locations show higher intrinsic subthreshold oscillation frequency in vitro and smaller spacing between grid fields in vivo. Plots of the reciprocal of temporal frequency versus anatomical position (Fig. 4A) reveal slopes similar to the slope of grid field spacing (G) relative to anatomical position (Fig. 4B).

Fig. 4.

Subthreshold oscillations may underlie differences in grid field spacing. (A) Reciprocal of oscillation frequency plotted versus anatomical depth for anatomical range matching a previous publication (7) as shown in (B). (C and D) Simulations of grid cell model (3, 8). Top, gray indicates trajectory of simulated rat. Red indicates firing locations. Bottom, firing rate map (red, maximum; blue, no spikes). Grid field spacing is small for mean dorsal oscillation frequency of 6.42 Hz (C) and larger for mean ventral frequency of 4.23 Hz (D).

These data support the prediction of a model (SOM text) (3, 8) related to other models of grid cells and theta phase precession (1820). In this model, grid cell periodicity arises from an interference pattern generated by intrinsic temporal oscillations in the soma and dendrites of a single cell. During simulated rat movement, cells modulated by head direction and speed (7, 21, 22) shift the frequency of dendritic oscillations (consistent with voltage effects on frequency). The grid pattern is the product of interference by three dendritic oscillations, each receiving a different head direction input, shifting in and out of phase with soma oscillations in proportion to distance moved in the preferred direction of each head direction cell. Spiking occurs when all three dendrites are in phase with the soma, causing oscillations to cross threshold (figs. S4 and S5). Spiking does not alter soma or dendritic phase but occurs with theta rhythmicity, consistent with in vivo recordings in entorhinal cortex (57) and hippocampus (2, 9, 18) and potentially causing precession relative to field potential oscillations (3, 8, 19, 20).

The model (3, 8) was modified to include shifts in dendritic frequency proportional to soma frequency and to use a scaling factor, H=fG (fig. S6), determined from the experimental data. Simulations (SOM text) (Fig. 4, C and D) demonstrated that differences in temporal frequency of somatic oscillations result in different grid field spacing. Insertion of the experimentally determined value for subthreshold oscillation frequency (f) at a particular anatomical location resulted in simulated grid cell spacing that matches data (7) on grid cell spacing (G) at the same anatomical location (figs. S4 and S7).

The model demonstrates one possible mechanism for path integration (2) and supports the prediction that systematic variation in gain of a movement-speed signal could underlie differences in grid field spacing (2, 9). The model is compatible with maintenance of grid cell representations by persistent firing (23) or attractor dynamics arising from patterned excitatory connectivity (2, 24, 25). Long-term stability of grid fields could require place cell input dependent on external landmarks (3, 8, 26). Differences in intrinsic frequency along the dorsal-ventral axis of entorhinal cortex could contribute to differences in place field size along the septal-to-temporal axis of the hippocampus (9, 27).

The systematic differences in intrinsic temporal frequency may provide multiple scales for coding of both space and time. Periodic representation of the environment at multiple spatial scales could prove essential to mechanisms of path integration (2, 3, 9, 28), consistent with impairments after entorhinal lesions (1). Coding of continuous dimensions by interacting frequencies could also allow the coding of continuous relative time necessary for episodic memory (29). These results suggest that, beyond simple summation of input, neural processing involves interactions of synaptic input and interference between intrinsic frequencies (8, 9, 18, 19).

Supporting Online Material

Materials and Methods

SOM Text

Figs. S1 to S7

References and Notes

View Abstract

Navigate This Article