Effects of visual inputs on neural dynamics for coding of location and running speed in medial entorhinal cortex
Abstract
Neuronal representations of spatial location and movement speed in the medial entorhinal cortex during the ‘active’ theta state of the brain are important for memory-guided navigation and rely on visual inputs. However, little is known about how visual inputs change neural dynamics as a function of running speed and time. By manipulating visual inputs in mice, we demonstrate that changes in spatial stability of grid cell firing correlate with changes in a proposed speed signal by local field potential theta frequency. In contrast, visual inputs do not alter the running speed-dependent gain in neuronal firing rates. Moreover, we provide evidence that sensory inputs other than visual inputs can support grid cell firing, though less accurately, in complete darkness. Finally, changes in spatial accuracy of grid cell firing on a 10 s time scale suggest that grid cell firing is a function of velocity signals integrated over past time.
Introduction
The spatially periodic firing of grid cells in the medial entorhinal cortex (MEC) (Fyhn et al., 2008; Fyhn et al., 2004; Hafting et al., 2005) is widely recognized as a neuronal correlate of memory-guided and/or self-motion-based spatial navigation (McNaughton et al., 2006). Based on that idea, computational models of grid cells including attractor network models (Burak and Fiete, 2009; Fuhs and Touretzky, 2006), oscillatory interference models (Blair et al., 2008; Burgess, 2008; Burgess et al., 2007; Hasselmo, 2008; Zilli and Hasselmo, 2010), hybrid models of attractor networks and oscillatory interference (Bush and Burgess, 2014), and probabilistic learning models (Cheung, 2016) suggest that grid cells integrate running direction and traveled distance to maintain a spatially periodic firing pattern during path integration (Gil et al., 2018; Jacob et al., 2019). Path integration depends on a velocity signal which depends on vestibular and proprioceptive signals (Chen et al., 2019; Tennant et al., 2018) and also moving sensory features such as optic flow (Campbell et al., 2018; Raudies et al., 2012; Raudies and Hasselmo, 2015). However, any path integration mechanism ultimately suffers from the accumulation of error over distance traveled or elapsed time (Kraus et al., 2015; Sreenivasan and Fiete, 2011) necessitating error correction by static environmental sensory inputs (Evans et al., 2016; Giocomo, 2016) such as visual landmark cues (Campbell et al., 2018; Hafting et al., 2005) or environmental boundaries (Barry et al., 2007; Hardcastle et al., 2015) to maintain an accurate representation of relative self-location within a spatial environment. Experimental data from studies on rats support the hypothesized role of grid cells in performing self-motion-based path integration and error correction by external sensory cues: Self-motion-based signals such as head direction are necessary for grid cell firing (Campbell et al., 2018; Winter et al., 2015b), and grid cell firing persists during darkness with reduced spatial selectivity (Allen et al., 2014; Hafting et al., 2005) reflecting the loss of error correction by visual cues. In contrast, two more recent studies on mice reported a rapid and nearly complete loss of spatially periodic grid cell firing during complete darkness (Chen et al., 2016; Pérez-Escobar et al., 2016) suggesting that spatially periodic grid cell firing and path integration depend on error correction by visual cues. However, motion-related interoceptive representations, such as proprioceptive, vestibular, and motor efference copy signals, can combine with visual representations to form a single multimodal representation supporting path integration (Tcheang et al., 2011) suggesting that external sensory cues other than visual cues may be sufficient to stabilize grid cells and support path integration. In this study, we therefore investigated whether and how the selective removal of visual inputs affects grid cell firing in mice under conditions which leave other landmark features such as walls, auditory, and olfactory cues still available to the animal.
Errors in path integration originate primarily in the velocity signal (Stangl et al., 2020), which combines information about movement direction and movement speed. However, the source of both the directional signal and the speed signal are still under debate. The direction signal has been attributed to the head direction signal terminating in the MEC (Sargolini et al., 2006; Taube et al., 1990a; Taube et al., 1990b; Taube, 2007), but head direction does not always match movement direction (Raudies et al., 2015). Likewise, the dominant source of the speed signal is still under debate. One potential candidate for the speed signal is the linear modulation of firing rates in a subpopulation of ‘speed cells’ in the MEC (Kropff et al., 2015). However, such a speed signal by firing rate has been shown to be more heterogeneous than previously assumed (Hardcastle et al., 2017) and not accurate on sub-second-long time scales (Dannenberg et al., 2019) as would be required for updating grid cell firing in real time. An alternative or complementary oscillatory speed signal is provided by the positive correlation between the theta rhythmic firing frequency of theta modulated MEC neurons and the running speed of an animal and/or the positive correlation between the local field potential (LFP) theta frequency and the running speed of an animal (Hinman et al., 2016; Jeewajee et al., 2008a; McFarland et al., 1975; Rivas et al., 1996; Shin and Talnov, 2001; Sławińska and Kasicki, 1998; Whishaw and Vanderwolf, 1973). The slope of the LFP theta frequency vs. running speed relationship in particular has been linked to spatial cognition and locomotion in both hippocampus and MEC (Korotkova et al., 2018; Monaghan et al., 2017; Wells et al., 2013). Interestingly, environmental novelty selectively reduces the slope of the LFP theta vs. running speed relationship (Jeewajee et al., 2008b; Wells et al., 2013) and also results in an increase in grid scale (Barry et al., 2012) as predicted by oscillatory interference models of grid cell firing (Burgess, 2008). In contrast, the y-intercept of the LFP theta vs. running speed relationship has been linked to emotion-sensitivity, in particular anxiolysis (Korotkova et al., 2018; Monaghan et al., 2017; Wells et al., 2013). However, little is known about the time courses of changes in coding of location by grid cells and coding of running speed by theta oscillations or firing rate (for instance during changes from light to darkness) although changes in grid cell firing, theta rhythmic dynamics, and speed modulation of firing rates may underlie many cognitive functions on a behavioral time scale ranging from seconds to minutes. By examining the temporal correlations between grid cell spatial stability and the two proposed speed signals—by firing rate or theta oscillations—we test if observed changes in grid cell firing can be explained by changes in the speed signal by firing rate or changes in the speed signal by LFP theta frequency. We further test if changes in the representation of a spatial metric by grid cells in the MEC is accompanied by similar changes in the tuning of head direction cells, border cells, or both. In addition, data presented in this study characterize the visual input-evoked changes in neural dynamics in the MEC. These data are relevant for models of path integration and for our understanding of how visual inputs modulate the coding of location and running speed in the entorhinal cortex.
Results
The dynamic range of speed coding by local field potential theta frequency decreases in darkness
To study changes in neural dynamics after transitioning from foraging in an illuminated environment (light) to foraging in complete darkness (dark) and vice versa, we recorded single unit and LFP activity in mice freely exploring an open-field environment during alternating light and dark epochs (Figure 1A). We first studied if the absence of visual inputs during darkness changed the relationship between LFP theta frequency and running speed. We found a large change in the slope and a smaller change in the y-intercept of the LFP theta frequency vs. running speed relationship during the dark condition (t-statistic of interaction effect (slope): 18.045, df = 806, p=2×10−61; t-statistic of change in y-intercept: −3.676, df = 806, p=0.0003; n = 15 mice) (Figure 1B). In contrast, the theta amplitude vs. running speed relationship did not change significantly (t-statistic of interaction effect (slope): 1.7343, df = 806, p=0.0832; t-statistic of change in y-intercept: −1.764, p=0.0781; n = 15 mice) (Figure 1C). The average running speed during light, 11.6 ± 0.5 cm/s, slightly decreased during darkness by 0.89 ± 0.31 cm/s (7.7%) to 10.72 ± 1.4 cm/s (mean ± s.e.m., p=0.0117, n = 15 mice, Wilcoxon’s signed-rank test) (Figure 1D and Figure 1—figure supplement 1A). Similarly, running acceleration also slightly decreased from 20.8 ± 1.9 cm/s2 in the light condition by 1.7 ± 0.7 cm/s2 (8.0%) to 19.1 ± 2.3 cm/s2 in the dark condition (p=0.0070, n = 15 mice, Wilcoxon’s signed-rank test) (Figure 1E). We therefore examined if the observed decrease in the slope of the LFP theta frequency vs. running speed relationship is independent from the decrease in running acceleration during darkness. If truly independent, the y-intercept of the correlation between LFP theta frequency and running acceleration is expected to decrease during the dark condition. In fact, the y-intercept of the LFP theta frequency vs. running acceleration relationship was significantly lower during darkness (t-statistic: 8.115, p=1.16×10−6, n = 15 mice, Student’s paired t-test) (Figure 1F). Furthermore, restricting our analysis of the LFP theta frequency to running speed analysis to time points when running acceleration was near zero (0 ± 5 cm/s2) did not significantly change speed tuning curves of LFP theta frequency (Figure 1G). These data demonstrate that the LFP theta frequency is directly modulated by running speed and that the observed change in the slope of the LFP theta frequency vs. running speed relationship in darkness reflects a change in visual inputs.
The LFP theta frequency vs. running speed relationship has been hypothesized to serve as a speed signal important for coding of location during self-motion-based navigation including path integration. Having established that the removal of visual inputs changes the speed signal by LFP theta frequency, we further examined the time course of that change by analyzing sequential 10 s blocks after the start of light and dark epochs. Averaging data from 78 sessions from 14 mice, we found that the LFP theta frequency vs. running speed relationships in light and dark became more different with elapsed time after start of light and dark epochs at a time scale of tens of seconds (Figure 1H and Figure 1—figure supplement 1B). This divergence appeared to be largely driven by a change in the slope of the theta frequency vs. running speed relationship as opposed to changes in the y-intercept. We therefore quantified the time courses of those changes in slope and y-intercept during light and dark conditions. This analysis revealed a slow increase in slope over time during the light condition (Pearson’s R = 0.79, p=0.0001) and a slow decrease in slope over time during the dark condition (Pearson’s R = −0.79, p=0.0001) (Figure 1I). The increase in slope during the light condition was accompanied by a small, yet significant, decrease in the y-intercept (Pearson’s R = −0.77, p=0.0002) (Figure 1J). In contrast, the y-intercept did not change significantly over time during the dark condition (Pearson’s R = −0.16, p=0.5251) (Figure 1J). In addition to those slow components of change, analyzing the transition points between light and dark conditions revealed fast (<10 s) light-to-dark decreases and dark-to-light increases in the slope (Figure 1K) and y-intercept (Figure 1L) of the LFP theta frequency to running speed relationship.
Changes in theta rhythmic firing as a function of running speed depend on visual inputs
Local field potential activity in the MEC primarily reflects the summed activity of synaptic inputs to the MEC, which is likely to be correlated with firing patterns of MEC neurons. Many neurons in the MEC show a theta rhythmic firing pattern and—similarly to the LFP theta frequency—the frequency of theta rhythmic firing increases with running speed (Hinman et al., 2016). We asked if and how theta rhythmic firing in theta modulated MEC neurons (n = 342 neurons from 14 mice) is affected by removing and reinstating all visual inputs. Towards that aim we used an MLE approach (Climer et al., 2015) to fit a model of theta spiking rhythmicity to the observed spike train autocorrelations and identified the frequency and magnitude of theta rhythmic firing in theta modulated MEC neurons in light and dark conditions (see Figure 2A for one example neuron). 339 out of 342 neurons could be analyzed using this model. We observed a decrease in theta rhythmic firing frequency during the dark condition (Figure 2B and Figure 2—figure supplement 7). Similar to the observed changes in the LFP theta frequency, the decrease in theta rhythmic firing frequency was largely attributed to a decrease in the slope of its relationship to running speed (Figure 2C) as previously shown in rats (Hinman et al., 2016).
However, averaging data across cells can hide differences between cells and important characteristics of the population response. We therefore analyzed changes in theta rhythmic firing on the single cell and population level as a function of binned running speed (immobility, defined as running speed <3 cm/s; 5–10 cm/s; 10–15 cm/s; 15–20 cm/s; 20–25 cm/s) (n = 339 neurons, Figure 2D). To our surprise, we found that neurons continued to fire rhythmically during immobility periods in both light and dark conditions. However, theta rhythmic frequencies varied largely across neurons within a wide range of theta frequencies (4–12 Hz) during immobility. Intriguingly, theta rhythmic firing synchronized across neurons around a narrower theta-frequency range of 6–10 Hz when the mouse started running. Importantly, running speed-dependent theta dynamics differed between light and dark conditions (Figure 2D–G). In order to quantify the changes in theta rhythmic firing dynamics as a function of running speed, we analyzed theta rhythmic firing separately for each running speed bin. In the light condition, the number of significantly theta rhythmic neurons increased sharply from 181 during immobility to 287 at the onset of running. Interestingly, this number then steadily decreased to 242, 215, and 186 for higher running speeds. In contrast, the number of significantly theta rhythmic neurons sharply increased at the onset of running and then saturated close to the maximal number during darkness (157 during immobility, 322, 326, and 328 out of 339 neurons at higher running speeds) (Figure 2E). A linear regression model for running speeds in the range from 5 cm/s to 25 cm/s revealed a significant difference in the number of theta rhythmic neurons as function of running speed between light and dark conditions (p=0.0003). These data demonstrate that the ensemble of theta rhythmic neurons decreased with running speed in the light condition but increased to ceiling during the dark condition (Figure 2E) contributing to the stronger synchrony of theta rhythmic firing observed on the population level (Figure 2D). These results were consistent across different significance levels between alpha = 0.001 and alpha = 0.1 for classification of theta rhythmic neurons (Figure 2—figure supplement 1). We next identified the cells which remained theta rhythmic across all running speed bins and quantified changes in the frequency and magnitude of theta rhythmic firing as a function of running speed during light and dark conditions. The frequency of theta rhythmic firing increased with running speed during the light condition (p=4.2×10–8, linear mixed effects model), but not during the dark condition (p=0.8085, linear mixed effects model; significant difference in slopes between light and dark, p=6.8×10–6) (Figure 2F). Similarly, the magnitude of theta rhythmic firing increased with running speed during the light condition (p=0.0002, linear mixed effects model), but not during the dark condition (p=0.0845, linear mixed effects model; difference in slopes between light and dark, p=9.1×10–6) (Figure 2G). Similar results were obtained when different significance levels were used to identify theta rhythmic neurons (Figure 2—figure supplement 1). Notably, the amplitudes of the average autocorrelogram models across all cells were very similar between light and dark conditions (Figure 2D, third row) despite the observed differences in the magnitudes of theta rhythmic firing on the single cell level (Figure 2G). These data suggest that—during darkness— higher synchrony across MEC neurons on the population level counterbalances the lower magnitude of theta rhythmicity on the single cell level which would explain that we did not observe a change in LFP theta amplitude in the dark condition (Figure 1C).
Population characteristics of running speed-dependent theta rhythmic firing and their changes with removal of visual inputs may differ between principal neurons and interneurons. We therefore classified all theta rhythmic neurons in our data set (n = 342) into putative principal neurons (PNs, n = 273) and putative interneurons (INs, n = 69) based on the spike width and firing rate of these neurons (Figure 2—figure supplement 2E). In general, both the frequency and magnitude of theta rhythmic firing was higher in PNs compared to INs (Figure 2—figure supplements 2C and 3C). In the light condition, the frequency of theta rhythmic firing showed a similar increase as a function of running speed for both PNs and INs but was consistently higher in PNs compared to INs. In contrast, the magnitude of theta rhythmic firing increased as a function of running speed in INs but not in PNs, which showed high theta rhythmic magnitudes at all running speeds. In the dark condition, frequencies of theta rhythmic firing differed less between PNs and INs compared to the light condition but showed an opposite change as a function of running speed with values slightly increasing in INs but slightly decreasing in PNs. In contrast, the magnitude of theta rhythmic firing was significantly higher in PNs compared to INs at all running speeds. Statistical results on changes in the frequency and magnitude of theta rhythmic firing within and between light and dark conditions for PNs and INs are reported in Tables 1 and 2. Statistical results on differences between PNs and INs with respect to changes in theta rhythmic firing as a function of running speed during light and dark conditions are reported in Tables 3 and 4. We next addressed the question whether the observed changes in theta rhythmic spiking activity in PNs was different for different functional cell types such as grid cells, or alternatively was similar across PNs. We found that changes in theta rhythmic spiking activity as a function of running speed during and between light and dark conditions were similar across grid cells, head direction cells, and border cells (Figure 2—figure supplements 4, 5 and 6) and indistinguishable from all PNs, indicating that changes in theta rhythmic activity is ubiquitous across different functional cell types in the medial entorhinal cortex.
To study the time course of changes in theta rhythmic firing, we next analyzed theta rhythmic firing in sequential 10 s blocks after the start of light and dark epochs (Figure 2H). The frequency of theta rhythmic firing was significantly reduced in the dark condition (effect size: –0.381 Hz, t-statistic: –8.598, df: 4316, p=1×10–17, CI: [–0.468 –0.294]). Similar to the fast initial changes observed in the slope of the LFP theta frequency vs. running speed relationship, we observed a fast decrease in the frequency of theta rhythmic firing when transitioning from light to dark and a fast increase when transitioning from dark to light (Figure 2I). In contrast to the observed changes in the slope of the LFP theta frequency vs. running speed relationship, the frequency of theta rhythmic firing did not further change on a time scale of tens of seconds during light or dark epochs. However, differentiating between PNs and INs revealed a decrease in frequency of theta rhythmic firing as a function of time in both light and dark conditions in INs but not PNs (Figure 2—figure supplements 2F and 3E).
Theta-frequency phase locking of neuronal spiking increases during darkness
Theta-frequency phase locking of single neuron spiking is a prominent feature of the hippocampal and entorhinal network (Klausberger et al., 2003; Mizuseki et al., 2009). Increased phase locking has been associated with dynamical switching between encoding and retrieval dynamics (Hasselmo et al., 2002; Siegle and Wilson, 2014), memory performance (Liebe et al., 2012; Rutishauser et al., 2010; Schomburg et al., 2014), and spatial cognitive demand (Barry et al., 2016). We therefore asked if the larger synchronization across neurons observed in the dark condition was associated with an increase of theta-frequency phase locking of theta rhythmic spiking activity. To compare phase locking of firing between light and dark conditions, we computed the mean resultant length (MRL) of theta-phase firing rate maps for the light and dark condition (Figure 3A). Interestingly, phase locking of theta rhythmic neurons increased during the dark condition (median MRL (Light) = 0.29, median MRL (Dark) = 0.37, p=1×10–32, Wilcoxon signed-rank test; Figure 3B). Such an increase is consistent with the increased numbers of theta rhythmic neurons during the dark condition as shown in Figure 2E. No significant differences were observed between PNs and INs (Figure 3—figure supplement 1). Notably, each neuron’s preferred phase of firing remained unchanged between light and dark conditions (Figure 3C & D). We next investigated the time course of those changes in phase locking (Figure 3E). The time course of changes in theta-frequency phase locking was similar to the changes in the frequency of theta rhythmic spiking, showing a fast initial change but no further slow component of change. No differences were observed between PNs and INs.
Visual inputs alter mean firing rates of neurons in the medial entorhinal cortex
Motivated by the differences in LFP activity during light and dark conditions, we asked if these changes were accompanied by changes in firing rates of MEC neurons. We therefore first asked whether and how the removal and reinstating of visual inputs affected firing rates in MEC neurons. Consistent with a previous study (Pérez-Escobar et al., 2016), we found striking differences in mean firing rates between light and dark epochs among a wide variety of MEC neurons including speed cells, head direction cells, grid cells, border cells and interneurons (Figure 4A & B; Figure 4—figure supplement 1). In general, changes in the mean firing rates of MEC neurons could be positively or negatively correlated to visual inputs. Positive modulation by visual inputs was observed ~4 times more frequently than negative modulation (139 out of 842 or 16.5% positively modulated neurons; 32 out of 842 or 3.8% negatively modulated neurons at significance level alpha = 0.05; data from 15 mice). In both cases, changes in firing rate after transitions between light and dark conditions occurred in less than 5 s (Figure 4C & D). Since we observed slight but significant changes in the animals’ mean running speeds during darkness and many neurons’ firing rates are modulated by running speed, we first tested whether the observed light vs. dark changes in mean firing rates of neurons are merely a consequence of the interaction between speed modulation of firing rates and changes in the animal’s running speed. Towards that goal, we computed the dark-to-light ratios of changes in mean firing rates for negatively (0.62 ± 0.16) and positively (1.53 ± 0.28) modulated neurons, respectively, and compared these observed ratios with running speed-adjusted dark-to-light ratios computed from the areas under the speed tuning curves in light and dark conditions (Figure 4—figure supplement 2). After running speed adjustment, dark-to-light ratios were 0.61 ± 0.17 (mean adjustment = –0.014, CI = [–0.002 –0.026], t = –2.31, df = 138, p=0.023, t-test) and 1.36 ± 0.23 (mean adjustment = –0.174, CI = [–0.112 –0.235], t = –5.78, df = 31, p=2.3×10–6) for negatively and positively modulated neurons, respectively. Notably, even after adjustment for running speed, the dark-to-light effect sizes were very similar for negatively and positively modulated neurons, namely a 39% decrease and a 36% increase in firing rates, respectively. These data demonstrate that the observed dark versus light changes in mean firing rates are not an artifact of the interaction between speed modulation of firing rates and changes in the animals’ running speed.
One possible explanation for the changes in mean firing rates during darkness can be a binary change in the available spatial information gained from visual cues. Alternatively, changes in firing rates could be driven by changes in net excitatory drive from areas in the visual cortex due to changes in luminance and retinal stimulation independent of any gain or loss of spatial information. To distinguish between these two alternative hypotheses, we analyzed data from optogenetic control experiments in a subset of two mice, in which mice carried a light fiber attached to their headstage with the fiber tip approximately 1 cm above the animal’s head. Green laser light (561 nm, 1 mW) was emitted from the tip of that light fiber and illuminated the head and immediate surroundings of the animal’s head in alternating epochs. Importantly, this experiment was performed with the ceiling lights turned on during the whole recording session with the head-mounted green laser light being turned on and off in alternating epochs during that session. While being visible to the animal as an egocentrically stationary spot of light, the laser light was too weak to illuminate the open-field environment and thus did not change the visibility of distal visual cues useful for navigational purposes (e.g. walls or the cue card of the arena). Nevertheless, we found strong effects of that kind of visual point stimulation on firing rates of five MEC neurons from two mice during foraging under the light condition (Figure 4E & F). Four out of those five neurons showed positive modulation by the egocentrically stationary spot of laser light, and one showed negative modulation. The effects of presenting an egocentrically stationary visual stimulus during the light condition on baseline firing rates of MEC neurons resembled the effects of transitioning from the dark to the light condition with respect to effect sizes in individual neurons and with respect to the ~4:1 ratio in the number of neurons that are positively and negatively modulated by retinal stimulation. In summary, these data demonstrate that the changes in retinal stimulation due to changes in the amount of spatially non-informative light hitting the retina after transitioning from dark to light and vice versa are sufficient to explain the observed changes in mean firing rates of MEC neurons.
Changes in mean firing rates of speed cells result in changes in the slopes of speed tuning curves without changing the running speed-dependent gain
We next turned to the analysis of the hypothesized speed signal by firing rate, in particular speed cells in the MEC (Kropff et al., 2015), and asked if the speed modulation of firing rates in MEC neurons is changed in complete darkness. The slopes of speed tuning curves have previously been reported to change after the removal of visual inputs and this change has been interpreted as a change in the hypothesized speed signal by speed cells, which has further been indicated to underlie changes in grid cell firing and path integration during darkness (Chen et al., 2016; Pérez-Escobar et al., 2016). An alternative approach, however, analyzes the running speed-dependent gain in firing rates. When speed modulation is analyzed after normalization across conditions, any change in external or internal parameters affecting a neuron’s mean firing rate will inevitably also affect the y-intercept and the slope of the neuron’s speed tuning curve. In contrast, analyzing speed modulation after normalization within conditions allows analyzing the running speed-dependent gain in firing rates independently from changes in a neuron’s mean firing rate. In this study, we do not examine which of the two alternative measures of speed modulation–the slope of the speed tuning curve or the running speed-dependent gain in firing rate—is the more useful measure for speed modulation of firing rates. However, we argue in favor of the second approach, analyzing the running speed-dependent gain in firing rate after normalization of speed tuning curves within conditions, on the grounds that it is standard in the field to compare tuning of other parameters such as head directional tuning by head direction cells or spatial tuning by place cells, grid cells, or border cells after normalization within conditions. As shown above, the mean firing rates of many neurons are responsive to changes in visual inputs. In addition, many speed cells show conjunctive coding properties which may be changed in complete darkness (Hardcastle et al., 2017; Hinman et al., 2016). To our knowledge, previous analyses of changes in firing rate modulation by running speed between two experimental conditions, such as light and dark environments, have not been corrected for possible confounding factors, such as changes in mean firing rate or conjunctive coding properties. In particular, data presented in this study show that changes in visual inputs cause changes in the firing rates of neurons and those changes may be independent from possible additional changes in the firing rate modulation by running speed. The same logic applies for changes in conjunctive coding properties, which may be misinterpreted as changes in the running speed modulation of firing rate if not properly controlled for. To account for possible changes in conjunctive coding properties, such as conjunctive coding for position, head direction, and LFP theta phase, we applied a linear-nonlinear (LN) model (Hardcastle et al., 2017) which provides a method for identification of significantly speed-modulated neurons without restricting the analysis to linear modulation of firing rates and computes speed response curves instead of speed tuning curves. Speed response curves are speed tuning curves corrected for the influence of conjunctive coding properties (Hardcastle et al., 2017). The firing rates of n = 164 neurons out of 495 cells from 15 mice showed significant modulation by running speed in either the light or the dark condition. Consistent with previous reports, speed modulation of firing rates was heterogeneous (Figure 5A & B; Dannenberg et al., 2019; Hardcastle et al., 2017; Hinman et al., 2016). Furthermore, we noticed a distribution of peaks in the speed response curves covering the full range of running speeds (Figure 5A & B). Hierarchical clustering of speed cells resulted in two main clusters of 138 positively and 26 negatively modulated neurons (Figure 5—figure supplement 1). Consistent with previous results (Chen et al., 2016; Pérez-Escobar et al., 2016), we observed a significant decrease in the slope of speed tuning curves of positively speed-modulated neurons during the dark condition if speed response curves were normalized to the mean firing rate across conditions (Figure 5A & C and Table 5). 153 out of the 164 speed-modulated neurons could also be examined for potential co-modulation of firing rates by visual inputs. 34 neurons out of those 153 (22.2%) speed cells showed significant modulation of mean firing rate by visual inputs. This proportion is very similar to the proportion of 171 out of 842 (20.3%) neurons with significant modulation of mean firing rates by visual inputs observed in the whole data set. Notably, this is a conservative estimate because the firing rate of many more neurons may be modulated by visual inputs without that modulation being detected as significant with the given statistical power and an alpha level of 0.05. To address the question whether the observed changes in the slope of speed response curves are caused by changes in the running speed-dependent gain in firing rates or are instead artifacts of changes in the mean firing rate of neurons, we compared speed response curves between light and dark conditions after normalization within conditions as opposed to normalization across conditions. Note that normalization within conditions corrects for changes in mean firing rate between conditions. Interestingly, we found that speed response curves appeared very similar between light and dark conditions after correcting for changes in mean firing rate and conjunctive coding properties (Figure 5B &D). In particular, no significant difference was observed between speed responses of positively speed-modulated neurons in the dark condition compared to the light condition (Table 6). With respect to negatively speed-modulated neurons, we even observed a small increase in the running speed-dependent gain in firing rate during the dark condition (Figure 5D and Table 6). Likewise, speed scores (Pearson’s correlation coefficient between firing rate and running speed) of all speed-modulated neurons were not significantly changed between light and dark conditions on the population level (p=0.912, n = 164, Kolmogorov-Smirnov test) (Figure 5E) and were highly correlated between light and dark conditions on the single cell level (R = 0.86) (Figure 5F). Moreover, the numbers of speed-modulated cells with conjunctive coding of position, head direction and theta phase were unchanged in the dark condition (Figure 5G). In summary, these data demonstrate that the running speed-dependent gains in firing rate of MEC speed cells remain unaltered during the absence of visual inputs, even if mean firing rates and the slopes of speed tuning curves change substantially.
Grid cell spatial firing is maintained but less stable during darkness
While in earlier studies grid cells have been reported to show stable grid cell firing during complete darkness in rats and mice (Allen et al., 2014; Hafting et al., 2005), two more recent studies in mice showed a strong reduction in the hexagonal symmetry of grid cell firing during complete darkness (Chen et al., 2016; Pérez-Escobar et al., 2016). To test if this discrepancy is due to differences in the experimental settings, the focus in this study was on selectively eliminating visual inputs only, while leaving other environmental and/or sensory cues such as walls or odors and auditory cues potentially present in the testing arena intact. We first tested if grid cell firing can emerge in complete darkness. We therefore introduced mice to the familiar open-field arena during complete darkness (the experimenter was wearing night-vision goggles) and recorded grid cells during the initial dark session and a following light and dark session (10–20 min duration each). We found that spatial periodicity of grid cell firing was low—but not completely eliminated—in the initial dark session and then strongly increased during the first light session (Figure 6A & B). Interestingly, after this initial orientation, spatial periodicity of grid cell firing decreased again after removal of visual inputs, but remained significantly higher than in the initial dark session (Chi-square = 10.89, p=0.0043, n = 9 cells from two mice, Friedman test; post-hoc tests revealed p<0.05 comparing initial dark and light, and initial dark and second dark sessions). We next aimed to study the changes in grid cell spatial firing patterns after transitioning from light to dark in more detail. Towards that aim, we recorded grid cells during alternating light/dark epochs, always beginning with a light epoch followed by a dark epoch. 28 cells from six mice were identified as grid cells based on their grid scores (see Materials and methods). To compare spatial firing properties of grid cells between light and dark conditions, we concatenated all light and dark sessions to compute firing rate maps for light and dark conditions. We found that grid cell field centers remained stable, but grid fields appeared wider during complete darkness (Figure 6C). Quantification of spatial periodicity and spatial precision of grid cell firing by computing grid scores and spatial information (see Materials and methods) revealed a significant decrease in both spatial periodicity and spatial precision of spike locations (Figure 6D & E). However, firing rate maps are the results of integrating firing rates of grid cells over long time periods, thereby measuring only spatial variability, but not spatiotemporal variability of grid cell firing. This distinction becomes important when considering two alternative hypotheses with regard to the underlying cause for the observed broadening of grid fields in complete darkness. Grid fields could appear wider because grid cells fire over a longer distance throughout a single grid field traversal. Alternatively, grid fields could appear wider because grid cell firing becomes less stable at the level of single grid field traversals. To distinguish between these two alternative hypotheses, we used two complementary approaches, First, we computed the spatiotemporal correlation between the observed firing rate over time and the firing rate predicted by the animal’s location at any given point in time and the associated value in the firing rate map at that location. The results show a significant reduction in spatiotemporal correlation during complete darkness (Figure 6F). Second, we computed spike-triggered firing rate maps of grid cells for the light and dark condition. Spike-triggered firing rate maps are not affected by slow translational shifts of the grid pattern. Nevertheless, we observed a significant decrease in the grid scores computed from spike-triggered firing rate maps (p=0.0354, n = 15 cells, Wilcoxon signed-rank test) (Figure 6—figure supplements 1 and 2A) suggesting that the phase of each grid map remained stable in darkness. Furthermore, we did not observe a ring-like pattern in the spike-triggered firing rate maps suggesting that the orientation of each grid map remained stable in darkness too. Lastly, we measured the grid spacing in each of the spike-triggered firing rate maps. Grid spacing remained unchanged in darkness (median grid spacing = 42 cm in light and dark, p=1, Wilcoxon signed-rank test, Bayes factor = 0.2005) (Figure 6—figure supplement 2B). Taken together, these data support the hypothesis that grid fields in the firing rate map appear broader because the spatiotemporal accuracy of grid cell firing is reduced without changes in the phase, orientation, or spacing of grid maps.
Visual inputs sharpen the tuning curves of head direction cells and border cells in the medial entorhinal cortex
Grid cell firing may rely on an intact head direction signal, in particular on the firing of head direction cells (Winter et al., 2015a). We therefore investigated whether the reduction in spatial accuracy of grid cell firing was accompanied by changes in the head direction tuning of head direction cells in the MEC. We identified 67 cells with significant unidirectional head direction tuning which were classified as head direction cells (Sargolini et al., 2006; Taube et al., 1990b; Taube et al., 1990a). Interestingly, head direction cells maintained their head direction preference in the dark condition (Figure 7A & B and Figure 7—figure supplement 1) suggesting that animals remained oriented within the environment during complete darkness. However, the mean resultant lengths of head directional firing rate maps were decreased in the dark condition (Cohen’s d = 0.238, p=0.0337, Wilcoxon signed-rank test) suggesting that the accuracy of head directional tuning was reduced in darkness (Figure 7C). In addition to grid cells, border or boundary vector cells (Lever et al., 2010; Solstad et al., 2008) are hypothesized to contribute to a representation of space in the MEC. We identified 27 border cells (Figure 7D and Figure 7—figure supplement 2) and compared their border scores (Solstad et al., 2008) between light and dark conditions. Border scores were significantly decreased during darkness (Cohen’s d = –0.4975, p=0.0325, Wilcoxon signed-rank test) (Figure 7E) suggesting that representations of boundaries in the open-field environment were less accurate during darkness.
Changes in grid cell spatial stability as a function of time correlate with changes in a proposed speed signal by local field potential theta frequency
Our data set included data on n = 7 grid cells from three mice recorded over an exceptionally long time period across multiple days. These data provided sufficient statistical power and sufficient sampling of the environment to analyze the time courses of the observed changes in spatial periodicity, spatial selectivity and spatiotemporal stability of grid cell firing during light and dark conditions. Towards this goal, we first computed firing rate maps for each of the 10 s bins of activity to analyze the time course of spatial periodicity quantified by the grid score (Figure 8A & B). Second, we computed the spatial correlations between the firing rate maps of the 10 s binned data and the overall firing rate map for the light condition (Figure 8C). Third, we computed the spatiotemporal correlation for each of the 10 s bins (Figure 8D). We found a saturating exponential increase in spatial periodicity of grid cell firing (Figure 8B), spatial correlation (Figure 8C), and spatiotemporal correlation (Figure 8D & E) during the light condition. During darkness, we observed a similar, though weaker, trend in the opposite direction. Intriguingly, grid cell stability measured as spatiotemporal correlation of grid cell firing was strongly correlated with the slope of the LFP theta frequency vs. running speed relationship (R = 0.966, p=1.83×10–21, n = 36 (18 time points x two conditions); Figure 8F and Figure 8—figure supplement 1). The time courses of changes in grid cell stability and the slope of the LFP theta frequency vs. running speed relationship were also correlated within light and dark conditions (Light: R = 0.825, p=2.55×10–5, n = 18 time points; Dark: R = 0.490, p=0.039, n = 18 time points) further confirming the link between grid cell stability and the LFP theta frequency to running speed relationship. In contrast, changes in the y-intercept of the LFP theta frequency vs. running speed relationship as a function of time were negatively correlated with changes in grid cell stability in the light condition (R = –0.633, p=0.0048, n = 18 time points) and uncorrelated to changes in grid cell stability in the dark condition (R = 0.056, p=0.83, n = 18 time points) (Figure 8—figure supplement 2). Taken together, these data demonstrate that the representation of space by grid cell firing is correlated to the speed signal by LFP theta frequency and that both of these signals undergo plastic changes at a behavioral time scale of tens of seconds.
Discussion
Data presented in this study show how changes in visual inputs affect firing rates and theta rhythmic firing dynamics in the MEC. Furthermore, these changes correlate with changes in the representation of location by grid cells, representation of environmental boundaries by border cells, and representation of head direction by head direction cells in the MEC. Intriguingly, a proposed speed signal by firing rate—the modulation of firing rate by running speed—remains unaltered if the analysis is adjusted for changes in mean firing rate. However, an alternative speed signal by LFP theta frequency—the linear correlation between LFP theta frequency and running speed—is a function of visual inputs and reduced in darkness. Intriguingly, changes in the slope of the LFP theta frequency vs. running speed relationship are strongly correlated to the changes in spatial stability of grid cell firing. While changes in firing rates of MEC neurons and changes in theta rhythmic firing properties of MEC neurons occur almost instantaneously during the transitions between light and dark, changes in the proposed speed signal by LFP theta frequency and changes in grid cell stability show a bimodal change with a fast (<10 s) component and a further slow (10–180 s) saturating exponential component. Such changes in response to changes in sensory input are longer than previously addressed by theoretical and computational models of grid cell firing. Lastly, we report changes in the running speed-dependent population dynamics of theta rhythmic MEC neurons upon changes in visual inputs. Synchronization of theta rhythmic firing across principal and interneurons increased with the onset of running, and even more so during darkness. The stronger synchronization across neurons during darkness was accompanied by an increase in theta-frequency phase locking of spiking activity (Figure 3). While the frequency of theta rhythmic firing increased with running speed during the light condition in both principal and interneurons, no such changes in theta rhythmic firing frequency as a function of running speed were observed during darkness. Likewise, the magnitude of theta rhythmic firing of interneurons increased as a function of running speed during the light, but not during the dark condition.
It has previously been proposed that disruption of speed signals may cause disruption of spatially periodic grid cell firing in the MEC. Current computational models of grid cell firing and path integration require a linear speed signal (Burak and Fiete, 2009; Burgess, 2008; Burgess et al., 2007; Fuhs and Touretzky, 2006; McNaughton et al., 2006). Two candidates for such a speed signal terminating in the MEC have been proposed: a speed signal by firing rate, most clearly present in a subpopulation of ‘speed cells’ showing positive and linear speed tuning curves (Kropff et al., 2015); and an oscillatory speed signal by changes in LFP theta frequency or by changes in the theta rhythmic firing frequency of individual MEC neurons (Burgess, 2008; Hinman et al., 2016). Both speed signals can change in response to changes of the size and shape of the environment (Munn et al., 2020), as do grid cells (Barry et al., 2007; Krupic et al., 2015). In this study, we examined how changes in spatial stability of grid cells relate to changes in the potential speed signals by firing rate and LFP theta frequency as a function of time.
In order to analyze the time course of grid cell stability over time after transitions between light and dark, we took advantage of long-term recordings of the same grid cells over multiple days providing us with sufficient statistical power. Intriguingly, changes in spatial stability of grid cell firing showed an instantaneous component and a slow saturating exponential component at a behavioral time scale of tens of seconds after transitions between light and dark conditions. Notably, the observed saturating exponential component in changes of grid cell spatially periodic firing supports an earlier report by Pérez-Escobar et al., 2016. We next tested if the observed change in spatial stability of grid cell firing could be explained by a change in the potential speed signal by firing rate. We distinguished the effect of visual inputs on the mean firing rate of neurons from the effect of visual inputs on the running speed-dependent gain in firing rates, which is equivalent to the slope in speed response curves after normalization within conditions. We further corrected for multiplexed or conjunctive coding properties by computing speed response curves using a linear-nonlinear model (Hardcastle et al., 2017). After adjusting for changes in mean firing rates and conjunctive coding properties, we found that visual inputs had no impact on the running speed-dependent gain in firing rates of positive speed cells in the MEC, while the running speed-dependent gain was even slightly increased in negative speed cells. We therefore conclude that changes in the potential speed signal by firing rate cannot explain the change in spatially periodic grid cell firing. A previous study (Pérez-Escobar et al., 2016) has found that speed scores were changed between light and dark conditions. However, that study used a speed score threshold to identify speed cells which is biased for identification of linearly modulated neurons because the speed score is a measure of linear correlation. In contrast, the current study used the LN-model to identify speed cells. As a consequence, our data set includes many speed cells with non-linear speed response curves, where the speed score is a sub-optimal measure. Like the slope of the speed tuning curve, the speed score is also a function of firing rate and previously reported changes in the speed score of speed cells during darkness may have been driven by changes in mean firing rates as opposed to changes in the running speed-dependent gain in firing rate.
We next turned to the analysis of the time courses of changes in the two oscillatory speed signals, namely the speed signal by LFP theta frequency and the speed signal by theta rhythmic firing frequency. Notably, the neuronal distribution of theta rhythmic firing frequencies covers a broad range of theta (4–12 Hz) indicating that rhythmic firing of MEC neurons is only loosely coupled to the LFP theta frequency. Moreover—in contrast to LFP theta rhythmicity—theta rhythmic firing persisted during immobility, though significantly less synchronous across neurons. The observed changes in theta rhythmic firing frequency at the single cell and population levels provide a possible explanation for why the LFP theta amplitude vs. running speed relationship is not changed during complete darkness. During darkness, movement-related theta-rhythmicity magnitude is decreased on the single cell level, while theta synchronicity across neurons is increased on the population level. These two effects would influence the LFP theta amplitude in opposite directions and thus may cancel each other out. Intriguingly, the time course of changes in the slope of the LFP theta frequency vs. running speed signal matched the time course of changes in grid cell spatial stability (Figure 8F). In contrast, changes in theta rhythmic firing frequency showed only the fast (< 10 s) component of change (Figure 2H & I). These data on the parallel time course of changes in the oscillatory speed signal by LFP theta frequency and grid cell stability suggest that grid cell firing may be more strongly associated with changes in the mean attractor state of rhythmic firing in the MEC reflected by changes in the LFP theta signal than with changes in the rhythmic firing of individual neurons. Alternately, the exponential change in LFP theta may reflect a slow change in an oscillatory speed signal arising from other cortical structures such as the medial septum.
Other data support the role of LFP oscillatory dynamics in the generation of the grid cell spatial code. In contrast to the speed signal by firing rate (Dannenberg et al., 2019; Hinman et al., 2016), the theta oscillatory speed signal depends on the MSDB (Chrobak et al., 1989; Dannenberg et al., 2015; Givens and Olton, 1994; Hinman et al., 2016; Lawson and Bland, 1993), and pharmacological inactivation of MSDB neurons strongly reduces both theta oscillations and spatial periodicity of grid cell firing in the MEC (Brandon et al., 2011; Koenig et al., 2011). Pharmacological inactivation of the medial septum also impairs self-motion-based linear distance estimation (Jacob et al., 2017). Furthermore, passive transportation of rats reduces theta oscillatory power and abolishes both spatially periodic grid cell firing and the speed signal by LFP theta frequency (Winter et al., 2015b). Likewise, pharmacological inactivation of the vestibular system or preventing linear/translational head movements in a virtual reality setup abolish the modulation of LFP theta frequency by running speed and are associated with a decrease in the spatial information score and an increase in grid scale (Chen et al., 2018; Jacob et al., 2014). Interestingly, the increase in grid scale observed during path integration on a circular track is associated with a change in the frequency of LFP theta oscillations, but not with a change in the firing patterns of speed and head direction cells (Jacob et al., 2019). Moreover, a recent study demonstrated that artificial low-theta stimulation can affect an animal’s estimation of linear distance during approaches to a goal zone (Mouchati et al., 2020) further strengthening a potential link between LFP theta frequency and path integration. The slope of the LFP theta frequency vs. running speed relationship has thus been linked to spatial cognition and grid cell firing (Korotkova et al., 2018). In addition, speed modulation of hippocampal theta frequency has been shown to correlate with spatial memory performance (Richard et al., 2013). We demonstrate in this study that the slope of the LFP theta frequency to running speed relationship and grid cell stability change as a function of visual inputs and as a function of time in a familiar environment. One possible explanation for the correlation between changes in spatial accuracy of grid cell firing and the slope of the LFP theta frequency is that the latter serves as a speed signal in the path integration process generating the spatially periodic grid cell firing pattern as suggested by oscillatory interference models of grid cell firing (Burgess, 2008; Burgess et al., 2007; Hasselmo, 2008). The reduced slope may then result in a lower resolution for controlling a potential path integration mechanism due to the reduced dynamic range of the oscillatory speed signal resulting in less accurate readout and increased variability. Importantly, this would require a rescaling of the oscillatory speed signal which may be facilitated by knowledge of the distances between the borders of a familiar environment. Previous studies have shown that environmental novelty reduces the slope of the LFP theta frequency to running speed relationship to a similar degree as darkness (Wells et al., 2013) but resulted in an increase in grid scale (Barry et al., 2012) as opposed to an increase in spatial stability. Taken together, these data may indicate that in a novel environment, a reduced slope in the LFP theta frequency to running speed relationship may be the consequence of underestimating the distances between borders resulting in a larger grid scale but that in a familiar environment, the lack of error correction in a path integration mechanism may result in a reduction of the dynamic range in an oscillatory speed signal by LFP theta frequency resulting in less accurate speed coding.
Data presented in this study also show that the reduction in spatial accuracy of grid cell firing is accompanied by a reduction of accuracy in head directional tuning and a reduction of spatial accuracy of border cell firing (Figure 7). An alternative explanation for the observed reduction of spatial accuracy in grid cell firing during darkness is therefore a reduction of accuracy in head directional signals (Winter et al., 2015a), boundary signals, or both due to the loss of visual information. Another possibility is that the frequency of theta rhythmic firing of MEC neurons serves as a speed signal that is integrated over longer time scales than previously accounted for by computational models of grid cells (Dannenberg et al., 2019). In that scenario, both the slope of the LFP theta frequency and spatial periodicity of grid cell firing would depend on the history of a velocity signal, which would be consistent with the observed slow saturating exponential component of change. Such an integration over seconds has previously been shown to optimize the proposed speed code by firing rate (Dannenberg et al., 2019). Notably, the slow component of change in spatial accuracy of grid cell firing cannot be explained by error accumulation alone. While error accumulation could account for the observed slow decrease in spatial stability after removal of visual inputs, error accumulation cannot be the reason for the observed slow increase in spatial accuracy after reinstating visual inputs. If signals required for grid cell firing are integrated over time, however, such a slow increase is expected. Notably, a similar increase in spatial accuracy of grid cell firing was also reported recently by Weber and Sprekeler, 2019. Future experimental and computational modeling studies are needed to address the questions whether and how longer integration time windows could support grid cell firing.
The observed changes in neural correlates of running speed and location during complete darkness may be attributed to either a loss of visual cues for navigational purposes or alternatively to a net decrease in excitatory drive from areas in the visual cortex due to the lack of retinal stimulation (Campbell and Giocomo, 2018). In order to distinguish between these two possibilities, we took advantage of control data from an optogenetic study (Dannenberg et al., 2019), where a visible bright spot of light was moving geostationary to the animal without further illuminating the environment. During this manipulation, the amount of light hitting the retina increased without changing the visibility of static visual landmark cues. This perturbation resulted in very similar changes of firing rates in MEC neurons as observed during the light vs. dark paradigm (Figure 4). Importantly, changes of firing rates, theta rhythmic firing frequency, and theta-phase coupling of rhythmic firing occurred almost instantaneously upon changes in visual inputs (Figure 4C & D and Figure 2I). Taken together, these data suggest that changes in firing rates and theta rhythmic firing are primarily driven by changes in net excitatory input from visual cortex and do not reflect changes in the perception of the spatial aspects of the environment. Instead, our data suggest that changes in firing rates and theta rhythmic firing are driven by changes in retinal stimulation. In contrast, the modulation of firing rates by running speed is unaltered by changes in illumination or presence or absence of visual cues.
Two recent studies on freely moving mice reported a major disruption of spatially periodic grid cell firing in the absence of visual cues (Chen et al., 2016; Pérez-Escobar et al., 2016). However, grid cell firing is intrinsically linked to path integration-based navigation (Gil et al., 2018; Jacob et al., 2019) and results from studies using virtual reality demonstrate that physical motion has a greater influence on grid cell firing in mice than visual inputs (Campbell et al., 2018; Chen et al., 2019; Tennant et al., 2018). While the recording environment in the Chen et al. study did not differ significantly from the recording environment used in this study, the majority of light and dark trials performed in the Chen et al. study was discontinuous and, in most trials, mice were introduced to the recording environment in complete darkness. Likewise, we also found a significantly stronger disruption in spatially periodic grid cell firing when mice were introduced to the recording environment in complete darkness (Figure 6A and B). Conversely, Chen et al. found a less severe disruption in grid cell firing during darkness in the few trials where mice were introduced to a familiar environment with the room lights turned on. Interestingly, Chen et al. reported changes in the preferred head direction angles of head direction cells suggesting that mice were frequently disoriented. The study by Pérez-Escobar et al. used a circular maze environment surrounded by water with no walls to eliminate tactile and olfactory cues. In addition, a white noise generator was used to mask auditory cues. Finally, mice underwent an explicit disorientation procedure before the start of recordings. These recording settings resulted in an almost complete loss of head directional tuning and very likely left animals disoriented. In contrast to those two previous studies, this study investigates the specific contribution of visual inputs on spatially periodic grid cell firing and theta oscillatory network dynamics in the MEC by selectively manipulating visual inputs on alternating epochs but leaving other static environmental and sensory cues such as walls, auditory cues, and potential olfactory cues unaltered. During this more selective manipulation, grid cells in the MEC persisted during complete darkness as shown by unaltered orientation, phase, and scale of grid cell firing fields. However, accuracy of grid cell firing was reduced during darkness. This reduction of accuracy in grid cell firing was accompanied by a reduction in the accuracy of the representation of head direction and environmental boundaries in the MEC. While the preferred angle of head directional firing in head direction cells was not altered during darkness, the variability in head directional firing increased during darkness (Figure 7A–C). Likewise, border cells continued to fire in their respective place fields yet with less spatial accuracy (Figure 7D & E). Taken together, the data presented in this study suggest that animals need an initial orientation in the environment supported by static visual cues but can then use sensory cues other than visual cues, including wall encounters, to roughly maintain that orientation over long time periods. For maintaining a spatial representation in a completely dark environment, vestibular signals, motor efference copies, proprioceptive feedback, and external geometrical cues may be sufficient to maintain the orientation, scale, and phase of grid cell firing fields. During the light condition, static visual cues may offer additional cues to self-location since the distance and direction to these cues systematically change with the animal’s movement, and optic flow may be used for refining an estimate of self-motion.
Theta signals within the hippocampal formation have been suggested to not only reflect alterations in sensorimotor integration and the flow of sensory input, but also distinct cognitive operations (Hinman et al., 2011). Furthermore, much of theta physiology and function seems to be shared across rodents and primates including the correlation of theta frequency with movement speed (Goyal et al., 2020). However, there is currently little data on dissociations between the slope and y-intercept of the theta frequency to running speed relationship. Data presented in this study show that spatial accuracy of grid cell firing is intrinsically linked to the slope but not the y-intercept of the LFP theta frequency vs. running speed relationship. Slow changes in both the slope of the LFP theta frequency vs. running speed relationship and the spatial periodicity of grid cells after transitions between light and dark conditions suggest that the representation of space in a cognitive map may integrate velocity signals over past time and may therefore be history-dependent and updated at a slow behavioral time scale. The result of such an integration would be a representation of space which is malleable by internal memory processes.
Materials and methods
Subjects
Data were collected from a total of 17 adult 3–6 month-old male mice (Mus musculus) of wild type and transgenic backgrounds (Wild type, C57Bl/6J; ChAT-IRES-Cre, B6;129S6-Chattm2(cre)Lowl/J; PV-IRES-Cre, B6;129P2-Pvalbtm1(cre)Arbr/J). Experimental PV-IRES-Cre mice were homozygous and experimental ChAT-IRES-Cre mice were heterozygous. All experimental procedures were approved by the Institutional Animal Care and Use Committee for the Charles River Campus at Boston University. Mice underwent surgery for implantation of 2 movable tetrodes targeting the superficial layers of the MEC in the left hemisphere. After surgery, mice were separated for individual housing and maintained on a reversed 12 hr light/12 hr dark cycle.
Mice were habituated to the experimenter and testing room before the start of all experimental procedures. For the purpose of additional experiments not reported in this study, mice also underwent rAAV injections of Archaerhodopsin T or GCaMP6f into the medial septum and light fiber implantation targeting the medial septum (Dannenberg et al., 2019). Light delivery of ~1–15 mW into the medial septum during one or two test recording sessions did not show any effects and we did not continue light delivery for the rest of the study. We did not notice any differences between data recorded before those test stimulations and after those test stimulations and also did not observe any differences between those mice which had and had not been stimulated. No differences were observed between wild type or transgenic mice. Prior to surgery, mice were housed in Plexiglas cages together with their siblings.
Surgery for microdrive implantation in mice
Request a detailed protocolMice were injected with buprenorphine (0.1 mg/kg, s.c.) and atropine (0.1 mg/kg, i.p.), and chronically implanted under isoflurane anesthesia with 8-channel Axona microdrives (Axona, St. Albans, UK) carrying two movable tetrodes targeting the superficial layers of the dorsal MEC in the left hemisphere. One ground screw was implanted above the cerebellum, five anchoring screws were positioned across the skull, and a craniotomy was performed for implantation of the microdrive, 0.5 mm anterior to the edge of the transverse sinus and 3.4 mm lateral to the midline. Tetrodes were then implanted with a 6° polar angle and lowered 1 mm into the brain. The microdrive was secured to the skull with dental cement. Animals were given buprenorphine (0.1 mg/kg, i.p.), enrofloxacin (7.5 mg/kg, i.p.), and ketoprofen (3 mg/kg, i.p.) during a 5 day postsurgical care period and allowed one week in total to fully recover after surgery before the start of recordings.
In vivo electrophysiological recordings
Request a detailed protocolBefore implantation, impedance of tetrodes (17 µm platinum wires, California Fine Wires) was adjusted by electroplating to <150 kOhm. Tetrodes were lowered daily by 50 µm until the abundance of spikes and the presence of theta rhythmic cells in test recordings indicated that tetrode tips were located near or within the superficial cell layers of the MEC. Neural signals were pre-amplified by unity-gain operational amplifiers located on each animal’s head (head stage) and AC-wire-coupled to an Axona recording system (Axona Limited, St. Albans, UK). Spikes were threshold-detected and recorded at 48 kHz. The LFP signal was recorded single-ended with the ground screw as reference at 4800 Hz. Position and head direction of mice were tracked with headstage-mounted invisible infrared LEDs at 50 Hz. Spikes were clustered into single units using spike clustering software (DacqUSB/Tint, Axona, St. Albans, UK). Further offline analyses were performed using Matlab (MathWorks, Natick, MA, USA), the circular statistics toolbox (Berens, 2009) for Matlab, and R (R Foundation for Statistical Computing, Vienna, Austria).
Behavior during light and dark conditions
Request a detailed protocolAll recordings were performed while animals foraged for randomly scattered small pieces of Froot Loops (Kellogg Company, Battle Creek, MI, USA) in a 1 m square open-field arena with 60 cm high black walls and a white letter size cue card at one wall. Recordings were performed during alternating epochs of open-field illumination (light) and darkness (dark). The dark condition was completely dark from the perspective of the human experimenter after adaptation to the dark condition for at least five minutes. A ceiling light provided full illumination of the open field during the light condition. Turning that ceiling light off resulted in complete darkness, which was achieved by the following measures. The open-field environment was located at the end of the recording room surrounded by three walls with no doors or windows and a thick laser-proof black curtain covering the front side. All LED lights of all electrical equipment were completely shielded with black duck-tape. Tracking of position and head direction of the mouse was achieved by using invisible infrared LEDs. These infrared LEDs were additionally covered with tape to minimize their brightness and preventing a potential spectral bleed into the visible spectrum. The slits around the door to the recording room were covered with black shower curtain material to minimize light leaking into the space within the recording room outside the laser-proof curtain. The computer monitors outside of the curtain were also covered with black shower curtains. Importantly, the experimenter frequently tested for any potential light sources within or light leaks into the space around the open-field environment. Recordings only began when it was completely dark from the perspective of the human eye after having adopted to the dark condition for at least 5 min. Alternating epochs of recordings in light and darkness could vary in duration between 30 s and 30 min; most experiments used 3 min alternating epochs allowing high statistical power due to multiple repetitions from light to darkness in each session. A typical recording session lasted 40–60 min.
Experiments with alternating epochs of light and dark began with the light condition. Experiments with a dark-light-dark paradigm started in complete darkness. For those experiments, mice were brought into the room with their cage covered by a black curtain, and mice were plugged in and introduced into the arena during complete darkness by an experimenter wearing night-vision goggles. The experimenter then left the arena side of the curtain before recordings began.
Histology
Request a detailed protocolAfter the end of electrophysiological data collection, animals were deeply anesthetized by isoflurane or i.p. injection of Euthasol (390 mg/kg), followed by electro-lesions (20-µA, 16 s) through one channel of each tetrode and subsequent transcardial perfusion with saline followed by 10% buffered formalin (SF100-4, ThermoFisher Scientifc). Brains were extracted and stored in fixative for one day. Cresyl violet staining was performed on 50 µm sagittal slices of MEC to visualize tetrode tracks and confirm positions of recording sites in the MEC.
Data analysis
Computation of running speed and firing rate signals
Request a detailed protocolThe animal’s running speed was estimated by applying a Kalman filter to the positional data obtained from the video tracking sampled at 50 Hz (Fyhn et al., 2004; Hinman et al., 2016; Kropff et al., 2015). To obtain an instantaneous firing rate signal matching the sampling rate of the behavioral tracking, we calculated the firing rate as the number of spikes per video frame period (20 ms) and smoothed this signal with a 125 ms Gaussian filter.
Modulation of mean firing rates by visual inputs
Request a detailed protocolIn order to identify neurons, whose firing rates are modulated by the illumination of the environment (consistently either higher or lower firing rates during light or darkness), we computed the average firing rate for each light and dark epoch and tested if the differences between light epochs and the following dark epochs and the differences between dark epochs and the following light epochs reached statistical significance at alpha = 0.05 using Student’s paired t-test. To analyze the time course of changes in mean firing rates after transition from light to dark or from dark to light, we first computed the peristimulus time histogram with 5 s binning around the transition point and then computed the average of z-scored firing rates across neurons. Only cells which showed spiking activity across three light-to-dark and three dark-to-light epoch pairs were included into this analysis.
LFP theta rhythm
Request a detailed protocolLFP data were down-sampled to 600 Hz and power spectral density analyses were performed using custom written MATLAB scripts and the FieldTrip toolbox (Oostenveld et al., 2011). A theta-frequency range was determined for each session based on the range of the theta peak in the power spectral density plot (Dannenberg et al., 2019). The LFP signal was then band-pass-filtered at this theta-frequency range with a third order Butterworth filter, and the instantaneous theta amplitude and phase were determined via a Hilbert transform. The instantaneous theta frequency could then be inferred from the instantaneous phase. To analyze changes in the LFP theta frequency vs. running speed relationship over time, we first down-sampled the LFP theta-frequency signal to 50 Hz in order to match the temporal binning of the behavioral tracking. We next sectioned data into 10 s blocks after transitions between light and dark and pooled data from all sessions and all mice (Figure 1H–L). For each of those 10 s blocks of data, the slope and y-intercept of the correlation between LFP theta frequency and running speed were computed using a linear model. The same analysis as described above was applied to 10 s blocks around the transition points from light to dark or from dark to light (Figure 1K & L).
Theta rhythmic firing
Request a detailed protocolWe used maximum likelihood estimation (MLE) of a parametric model of the interspike intervals underlying the spike time autocorrelogram to analyze theta rhythmic firing, which has been proven a more sensitive and statistically rigorous method than traditional measures of spike time autocorrelogram features (Climer et al., 2015). In particular, this model allowed us to statistically test significance of theta modulation in MEC neurons and compare theta rhythmic firing frequency and theta-rhythmicity magnitude on a single cell level between light and dark conditions. The model tested for rhythmic firing in the frequency range between 1–13 Hz. For the analysis of the relationship between theta rhythmic firing and running speed (Figure 2), we first identified continuous temporal epochs larger than 400 ms, in which the running speed of the animal varied less than ±4 cm/s. We next sorted these epochs according to the average running speed during the epochs into 1 cm/s bins from 1 cm/s to 29 cm/s and used the MLE approach to fit an autocorrelogram model to each of those epochs. The minimum number of interspike intervals which need to be present in the autocorrelogram to compute an accurate model is 100 (Climer et al., 2015), and no autocorrelogram model was computed if the number of interspike intervals was smaller than 100 for the epoch of interest. For each neuron and running speed bin, an average autocorrelogram was computed across the autocorrelograms for each epoch. The maximal run time for the theta rhythmic firing vs. running speed analysis on the research computing cluster was set to two weeks resulting in the successful analysis of 340 out of 342 neurons. For the analysis of the time course of changes in theta rhythmic firing after transition from light to dark or from dark to light, the first two minutes after the transition were sectioned into 10 s blocks (n = 342 neurons, Figure 2H). Alternately, 10 s blocks around the transition points from light to dark or from dark to light were analyzed (n = 342 neurons, Figure 2I). For each neuron, an autocorrelogram model was fitted to each of those 10 s blocks of data using the MLE approach. For each of those 10 s blocks of data, the average theta rhythmic firing frequency across neurons could then be computed from the peaks of the average autocorrelogram model across neurons by taking the inverse of the first peak in the autocorrelogram model. To analyze running speed-dependent theta rhythmic firing on the population level, we first identified continuous epochs of running speed within the ranges of 0–3 cm/s (immobility), 5–10 cm/s, 10–15 cm/s, 15–20 cm/s, and 20–25 cm/s. We next used the MLE approach to model an autocorrelogram for each running speed bin of each neuron. Neurons with sufficient number of spikes resulting in at least 100 interspike intervals to compute an MLE-modeled autocorrelogram in at least one of the running speed bins were included into the analysis (n = 339 neurons).
Analysis of theta-frequency phase locking
Request a detailed protocolThe analysis of theta-frequency phase locking is only valid when LFP theta rhythm is present and a sufficiently large number of spikes is present to estimate the phase distribution. We therefore restricted our analysis to time points where the running speed of the mouse was larger than 5 cm/s and where the sample size was larger than 30 spikes. LFP theta phase and the mean resultant length were computed as described in detail in Dannenberg et al., 2015. In short, the instantaneous LFP theta phase was computed from the Hilbert transform of the LFP signal filtered within the theta-range (6–10 Hz). For higher accuracy, the boundaries of this theta range were adjusted separately for each animal as described in detail in Dannenberg et al., 2019. We next interpolated the LFP theta-phase signal and assigned each spike of a theta rhythmic neuron to the instantaneous LFP theta phase at the time of spiking.
Classification of putative principal neurons and putative interneurons
Request a detailed protocolFor classification of single units into putative principal neurons or interneurons, we computed the spike width and firing rate for each unit. The spike width was defined as the time between the peak and the trough of the spike waveform. The distribution of spike widths was bimodal with a minimum between 0.29 ms and 0.34 ms. Units with a spike width <0.29 ms or a firing rate >10 Hz were classified as putative interneurons, the remaining units were classified as putative principal neurons.
Computation of spatial firing rate maps
Request a detailed protocolOccupancy-normalized spatial firing rate maps were generated by dividing the open-field environment into 3 cm spatial bins. For each spatial bin, an occupancy-normalized firing rate was computed by dividing the total number of spikes falling into that bin by the total time the animal spent in that bin. Spatial rate maps were then smoothed by a 3 cm wide two-dimensional Gaussian kernel.
Grid score and grid cells
Request a detailed protocolThe grid score was used as a measure of hexagonal symmetry of grid fields. The grid score was calculated as follows. We first computed the two-dimensional autocorrelogram of the spatial rate map, extracted a ring that encased the six peaks closest to the center peak excluding the central peak and computed correlations between this modified spatial autocorrelogram and rotated versions of itself after correction of elliptical eccentricity as described in detail in the supplementary methods for computation of gridness measure three in Brandon et al., 2011. If hexagonal symmetry is present, those correlations will peak at multiples of 60 degrees and repeat themselves infinitely at a frequency of 6 peaks per rotational cycle. We therefore can compute the power spectrum of those correlations as a function of rotational angle. The power at 60 degree was then defined as the grid score. Compared to previously used grid scores, this grid score largely reduced the number of false-positive grid cells while reliably assigning high grid scores to clear grid cell examples. This grid score is defined as a measure between 0 and 1, with 0 meaning no hexagonal symmetry of firing fields in the spatial rate map or simply ‘gridness’ and one meaning perfect gridness. Grid cells were defined as cells with a grid score >= 0.19 and at least 300 spikes in either light or dark condition, or across the whole recording session for experiments that started in complete darkness.
Spatial information
Request a detailed protocolWe defined the spatial information as the delta entropy . is the maximal amount of randomness given by a completely random distribution of spike locations in the open-field environment defined by , with , where is the number of spatial bins. is the observed amount of randomness in the distribution of spike locations given by , where is the animal’s location, is the spatial bin, is the number of spatial bins and is the probability that when a spike is observed. can then be understood as an increase in the amount of order in the distribution of spike locations relative to a completely random (uniform) distribution of spike locations.
Spatial correlation
Request a detailed protocolFiring rate map similarities were quantified by the spatial correlation between two firing rate maps. To compute the spatial correlation, the spatial bins of the firing rate maps were transformed into a single vector. The spatial correlation was then computed as the Pearson’s correlation coefficient between two of those spatially binned firing rate vectors.
Spatiotemporal correlation
Request a detailed protocolTo analyze spatial stability of grid cell firing in the time domain, we used the firing rate map of the light condition to compute the expected firing rate for each of the 20 ms time bins of the animal’s trajectory. We next computed the Pearson’s correlation coefficient between the observed 20 ms binned firing rate and the expected 20 ms binned firing rate.
Spike-triggered firing rate maps
Request a detailed protocolSpike-triggered firing rate maps (Bonnevie et al., 2013; Chen et al., 2016; Pérez-Escobar et al., 2016) were constructed as follows. For each spike, occupancy-normalized firing rates were computed from the next 10 s of data and centered to the spatial bin at the time of the spike. All those firing rate maps were then averaged and smoothed with a Gaussian kernel (SD = 3 cm) to create the spike-triggered firing rate map. To compute the occupancy-normalized firing rate as a function of radial distance, the circumferences of circles with increasing radii ranging from zero to 95 cm were overlaid with the spike-triggered firing rate map. For each of those circles, the spatial firing rates of all bins falling on the circle’s circumference were then averaged to compute the spatial firing as a function of radial distance. The location of the first peak in the plot of the firing rate as function of radial distance defined the spacing of the grid cell.
Border score and border cells
Request a detailed protocolA border score was computed as described in detail by Solstad et al., 2008. Border scores ranged from –1 for cells with central firing fields to +1 for cells with firing fields that perfectly aligned with at least one entire wall of the open-field environment. Cells with a border score larger than 0.68 in either the light or dark condition were defined as border cells. The threshold of 0.68 was the 90th percentile of the distribution of border scores from all neurons computed from whole recording sessions.
Computation of head direction firing rate maps
Request a detailed protocolHead directional firing rate maps were generated by dividing the head directional parameter space into six-degree bins and computing the occupancy-normalized firing rate as the total number of spikes falling into that bin divided by the total time the animal spent in that bin.
Head directional tuning and head direction cells
Request a detailed protocolHead direction cells were classified based on the criteria reported in Peyrache et al., 2015 and Viejo and Peyrache, 2020. Only cells with at least 10 spikes in each of the light and dark conditions were included in the analysis. Cells were classified as head direction cells if each of the following was true. First, the inverse of the circular standard deviation was greater than one in either the light or the dark condition; second, the p-value of the Rayleigh test was smaller than 0.001 in either the light or the dark condition; third, the peak firing rate in the head directional field was larger than 1 Hz in either the light or the dark condition.
Identification and classification of speed-modulated neurons
Request a detailed protocolTo identify speed-modulated neurons, we used a linear non-linear Poisson (LN) model (Hardcastle et al., 2017) which tested for significant modulation of firing rates by position, head direction, LFP theta phase, and/or running speed (Figure 5). For a robust analysis, we first screened for cells which showed stable firing across the whole recording session. Cells, whose mean firing rates were below 1 Hz or whose mean firing rates differed more than by a factor of two between the first and second halves of the recording session, were excluded from the analysis. To compare speed modulation of firing rates on the population level, we fitted the data with a quadratic polynomial curve (poly2 fit, MATLAB) and clustered neurons into positively and negatively speed-modulated neurons using hierarchical clustering (clustergram, bioinformatics toolbox, MATLAB).
Speed score
Request a detailed protocolThe speed score was calculated as the Pearson’s product moment correlation between instantaneous firing rate and running speed (Kropff et al., 2015).
Experimental design and statistical analysis
Request a detailed protocolDetails on statistical analyses on each experiment are reported in the Results section. Statistical tests used in this study include Student’s paired t-test and nonparametric statistical tests for repeated measures, namely Wilcoxon’s signed-rank test and Friedman’s test. post-hoc testing was performed if the Friedman test yielded a significant result at alpha = 0.05. Data used for analysis of differences in neural dynamics between light and dark conditions were obtained from recordings in a total of 15 mice. Data included LFP signals from 14 out of those 15 mice and 1267 single units from all 15 mice. Out of 1267 recorded single units, 495 single units with stable firing across the recording session and mean firing rates larger than 1 Hz were included in the LN-model analysis of firing rate modulation by running speed; 842 single units which showed spiking activity over at least three light-to-dark and three dark-to-light epoch pairs were used for the analysis of firing rate modulation by visual inputs. Data with green laser light stimulation (related to Figure 4) were recorded in two additional mice.
Code accessibility
Request a detailed protocolCustom MATLAB code used for data analysis can be found at https://github.com/hasselmonians or can be made available upon request. In particular, MATLAB code for analyzing single unit activity can be found at (https://github.com/hasselmonians/CMBHOME; Chapman and Bogaard, 2020; copy archived at swh:1:rev:fd5829a812414144df7ff76a177eb6b7263dbc5c); MATLAB code for analyzing firing rate modulation by light and dark conditions can be found at (https://github.com/hasselmonians/light-modulation; Dannenberg et al., 2020; copy archived at swh:1:rev:b949888144712b048e20812f4c25c7c3394ee4a6); and MATLAB code for analyzing single unit theta-rhythmicity can be found at (https://github.com/hasselmonians/mle_rhythmicity; Climer, 2020; copy archived at swh:1:rev:7c50ba73fa34dbe2586fb0f9e0204bd1e28f3209). Violin plots were created using R software (R Development Core Team, 2019).
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files. Matlab code for Figure 4 is provided on the laboratory's GitHub page (https://github.com/hasselmonians/light-modulation; copy archived at https://archive.softwareheritage.org/swh:1:rev:b949888144712b048e20812f4c25c7c3394ee4a6/ and https://github.com/hasselmonians/mle_rhythmicity copy; archived at https://archive.softwareheritage.org/swh:1:rev:7c50ba73fa34dbe2586fb0f9e0204bd1e28f3209/). Source data files have been provided for Figures 2, 3, 4, 5, and 6.
References
-
Experience-dependent rescaling of entorhinal gridsNature Neuroscience 10:682–684.https://doi.org/10.1038/nn1905
-
CircStat : a MATLAB toolbox for circular statisticsJournal of Statistical Software 31:1–21.https://doi.org/10.18637/jss.v031.i10
-
Grid cells require excitatory drive from the HippocampusNature Neuroscience 16:309–317.https://doi.org/10.1038/nn.3311
-
Accurate path integration in continuous attractor network models of grid cellsPLOS Computational Biology 5:e1000291.https://doi.org/10.1371/journal.pcbi.1000291
-
An oscillatory interference model of grid cell firingHippocampus 17:801–812.https://doi.org/10.1002/hipo.20327
-
A hybrid oscillatory interference/continuous attractor network model of grid cell firingThe Journal of Neuroscience 34:5065–5079.https://doi.org/10.1523/JNEUROSCI.4017-13.2014
-
Self-motion processing in visual and entorhinal cortices: inputs, integration, and implications for position codingJournal of Neurophysiology 120:2091–2106.https://doi.org/10.1152/jn.00686.2017
-
SoftwareCMBHOME: A Custom MATLAB Class for Neural Data, version swh:1:rev:fd5829a812414144df7ff76a177eb6b7263dbc5cSoftware Heritage.
-
Probabilistic learning by rodent grid cellsPLOS Computational Biology 12:e1005165.https://doi.org/10.1371/journal.pcbi.1005165
-
Intraseptal administration of muscimol produces dose-dependent memory impairments in the ratBehavioral and Neural Biology 52:357–369.https://doi.org/10.1016/S0163-1047(89)90472-X
-
Softwaremle_rhythmicity, version swh:1:rev:7c50ba73fa34dbe2586fb0f9e0204bd1e28f3209Software Heritage.
-
Synergy of direct and indirect cholinergic septo-hippocampal pathways coordinates firing in hippocampal networksThe Journal of Neuroscience 35:8394–8410.https://doi.org/10.1523/JNEUROSCI.4460-14.2015
-
SoftwareLight/Dark Modulation, version swh:1:rev:b949888144712b048e20812f4c25c7c3394ee4a6Software Heritage.
-
How environment and self-motion combine in neural representations of spaceThe Journal of Physiology 594:6535–6546.https://doi.org/10.1113/JP270666
-
A spin glass model of path integration in rat medial entorhinal cortexJournal of Neuroscience 26:4266–4276.https://doi.org/10.1523/JNEUROSCI.4353-05.2006
-
Impaired path integration in mice with disrupted grid cell firingNature Neuroscience 21:81–91.https://doi.org/10.1038/s41593-017-0039-3
-
Environmental boundaries as a mechanism for correcting and anchoring spatial mapsThe Journal of Physiology 594:6501–6511.https://doi.org/10.1113/JP270624
-
Local modulation of basal forebrain: effects on working and reference memoryThe Journal of Neuroscience 14:3578–3587.https://doi.org/10.1523/JNEUROSCI.14-06-03578.1994
-
Septotemporal variation in dynamics of theta: speed and habituationJournal of Neurophysiology 105:2675–2686.https://doi.org/10.1152/jn.00837.2010
-
Vestibular control of entorhinal cortex activity in spatial navigationFrontiers in Integrative Neuroscience 8:38.https://doi.org/10.3389/fnint.2014.00038
-
Medial entorhinal cortex and medial septum contribute to self-motion-based linear distance estimationBrain Structure and Function 222:2727–2742.https://doi.org/10.1007/s00429-017-1368-4
-
Reconciling the different faces of hippocampal theta: the role of theta oscillations in cognitive, emotional and innate behaviorsNeuroscience & Biobehavioral Reviews 85:65–80.https://doi.org/10.1016/j.neubiorev.2017.09.004
-
Relationship between hippocampal theta activity and running speed in the ratJournal of Comparative and Physiological Psychology 88:324–328.https://doi.org/10.1037/h0076177
-
Path integration and the neural basis of the 'cognitive map'Nature Reviews Neuroscience 7:663–678.https://doi.org/10.1038/nrn1932
-
Entorhinal velocity signals reflect environmental geometryNature Neuroscience 23:239–251.https://doi.org/10.1038/s41593-019-0562-5
-
FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological dataComputational Intelligence and Neuroscience 2011:1–9.https://doi.org/10.1155/2011/156869
-
Internally organized mechanisms of the head direction senseNature Neuroscience 18:569–575.https://doi.org/10.1038/nn.3968
-
SoftwareR: A Language and Environment for Statistical ComputingR Foundation for Statistical Computing, Vienna, Austria.
-
Modeling the influence of optic flow on grid cell firing in the absence of other cues1Journal of Computational Neuroscience 33:475–493.https://doi.org/10.1007/s10827-012-0396-6
-
Grid cells generate an analog error-correcting code for singularly precise neural computationNature Neuroscience 14:1330–1337.https://doi.org/10.1038/nn.2901
-
Sources of path integration error in young and aging humansNature Communications 11:2626.https://doi.org/10.1038/s41467-020-15805-9
-
The head direction signal: origins and sensory-motor integrationAnnual Review of Neuroscience 30:181–207.https://doi.org/10.1146/annurev.neuro.29.051605.112854
-
A local measure of symmetry and orientation for individual spikes of grid cellsPLOS Computational Biology 15:e1006804.https://doi.org/10.1371/journal.pcbi.1006804
-
Novelty and anxiolytic drugs dissociate two components of hippocampal theta in behaving ratsJournal of Neuroscience 33:8650–8667.https://doi.org/10.1523/JNEUROSCI.5040-12.2013
-
Passive transport disrupts grid signals in the parahippocampal cortexCurrent Biology 25:2493–2502.https://doi.org/10.1016/j.cub.2015.08.034
-
Coupled noisy spiking neurons as velocity-controlled oscillators in a model of grid cell spatial firingJournal of Neuroscience 30:13850–13860.https://doi.org/10.1523/JNEUROSCI.0547-10.2010
Article and author information
Author details
Funding
National Institutes of Health (R01MH60013)
- Holger Dannenberg
- Hallie Lazaro
- Pranav Nambiar
- Alec Hoyland
- Michael E Hasselmo
National Institutes of Health (R01MH120073)
- Holger Dannenberg
- Hallie Lazaro
- Pranav Nambiar
- Alec Hoyland
- Michael E Hasselmo
National Institutes of Health (R01MH052090)
- Holger Dannenberg
- Hallie Lazaro
- Pranav Nambiar
- Alec Hoyland
- Michael E Hasselmo
National Institutes of Health (K99NS116129)
- Holger Dannenberg
Office of Naval Research (MURI N00014-16-1-2832)
- Holger Dannenberg
- Hallie Lazaro
- Pranav Nambiar
- Alec Hoyland
- Michael E Hasselmo
Office of Naval Research (MURI N00014-19-1-2571)
- Holger Dannenberg
- Hallie Lazaro
- Pranav Nambiar
- Alec Hoyland
- Michael E Hasselmo
Office of Naval Research (DURIP N00014-17-1-2304)
- Holger Dannenberg
- Hallie Lazaro
- Pranav Nambiar
- Alec Hoyland
- Michael E Hasselmo
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the National Institute of Mental Health, grant numbers R01 MH60013, R01 MH120073 and R01 MH052090, by the National Institute of Neurological Disorders and Stroke, grant number K99NS116129, and by the Office of Naval Research MURI N00014-16-1-2832 and MURI N00014-19-1-2571, and Office of Naval Research DURIP N00014-17-1-2304.
Ethics
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All handling of animals and experimental procedures were approved by the Institutional Animal Care and Use Committee (IACUC) for the Charles River Campus at Boston University under protocol #16-008. All surgery was performed under isoflurane anesthesia and buprenorphine analgesia, and every effort was made to minimize suffering.
Copyright
© 2020, Dannenberg et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,389
- views
-
- 293
- downloads
-
- 29
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Specialized chemosensory signals elicit innate social behaviors in individuals of several vertebrate species, a process that is mediated via the accessory olfactory system (AOS). The AOS comprising the peripheral sensory vomeronasal organ has evolved elaborate molecular and cellular mechanisms to detect chemo signals. To gain insight into the cell types, developmental gene expression patterns, and functional differences amongst neurons, we performed single-cell transcriptomics of the mouse vomeronasal sensory epithelium. Our analysis reveals diverse cell types with gene expression patterns specific to each, which we made available as a searchable web resource accessed from https://www.scvnoexplorer.com. Pseudo-time developmental analysis indicates that neurons originating from common progenitors diverge in their gene expression during maturation with transient and persistent transcription factor expression at critical branch points. Comparative analysis across two of the major neuronal subtypes that express divergent GPCR families and the G-protein subunits Gnai2 or Gnao1, reveals significantly higher expression of endoplasmic reticulum (ER) associated genes within Gnao1 neurons. In addition, differences in ER content and prevalence of cubic membrane ER ultrastructure revealed by electron microscopy, indicate fundamental differences in ER function.
-
- Medicine
- Neuroscience
Pain after surgery causes significant suffering. Opioid analgesics cause severe side effects and accidental death. Therefore, there is an urgent need to develop non-opioid therapies for managing post-surgical pain. Local application of Clarix Flo (FLO), a human amniotic membrane (AM) product, attenuated established post-surgical pain hypersensitivity without exhibiting known side effects of opioid use in mice. This effect was achieved through direct inhibition of nociceptive dorsal root ganglion (DRG) neurons via CD44-dependent pathways. We further purified the major matrix component, the heavy chain-hyaluronic acid/pentraxin 3 (HC-HA/PTX3) from human AM that has greater purity and water solubility than FLO. HC-HA/PTX3 replicated FLO-induced neuronal and pain inhibition. Mechanistically, HC-HA/PTX3-induced cytoskeleton rearrangements to inhibit sodium current and high-voltage activated calcium current on nociceptive DRG neurons, suggesting it is a key bioactive component mediating pain relief. Collectively, our findings highlight the potential of naturally derived biologics from human birth tissues as an effective non-opioid treatment for post-surgical pain. Moreover, we unravel the underlying neuronal mechanisms of pain inhibition induced by FLO and HC-HA/PTX3.