Supplementary Materials

Total Page:16

File Type:pdf, Size:1020Kb

Supplementary Materials 1 Supplementary Information 2 Theory of weakly coupled oscillators 3 Detailed reviews and mathematical descriptions of the theory, also its extensions and limitations, 4 can be found in a number of publications1–19. In essence, many oscillatory phenomena in the 5 natural world represent dynamic systems with a limit-cycle attractor2. Even though the 6 underlying system might be complex (e.g. a neuron or neural population), the dynamics of the 7 system can be reduced to a phase-variable if the interaction strength among oscillators is weak. 8 In general, coupled oscillators can interact through adjustments of their amplitudes and phases, 9 yet if interaction strength is weak, amplitude changes are small and play a minor role in the 10 oscillatory dynamics (no strong deviation from the limit-cycle) and only the adjustments of the 11 phase are essential for understanding the behaviour. The transformation of a complex system to a 12 phase variable, if valid, reduces the dimensions of the problem, thereby allowing exact 13 mathematical investigation. The manner by which mutually coupled oscillators adjust their 14 phases, either by phase-delay or phase-advancement, is described by the phase response curve 15 (also called phase resetting curve), the PRC7,9,14–16,19–21. The PRC is important, because if the 16 PRC of a system can be described, the synchronization behaviour of the system can be 17 understood and predicted. Many biological oscillators are inherently noisy or chaotic. Therefore, 18 it is important to take this variability of the dynamics into account for a better understanding of 19 biological data. Cortical neurons in vivo exhibit irregular spiking patterns24 and neural networks 20 oscillations also show significant variability over time25. This type of variation is referred to as 21 phase noise and is distinct from measurement noise with the latter being unrelated to the 22 dynamics of the system. 23 According to the theory of weakly coupled oscillators, the synchronization of two coupled 24 oscillators can be predicted from the forces they exert on each other as a function of their 25 instantaneous phase difference. This function is referred to as the phase-response curve (PRC). 26 Accordingly, the phase evolution of two given cortical V1 locations is reduced to: 27 28 1 112121121 ε H 29 2 2 2ε 12H 12 2 1 2 30 31 where φ1,2 is the phase, 1,2 its temporal derivative, ω1,2 is the preferred frequency, ε12 and ε21 are 32 the interaction strengths, H12 and H21 are the single PRCs and ƞ1,2 is a phase-noise term with ƞ1,2 33 ~ N(0, σ2) N being the normal distribution. The two equations, as given in the main text equation 34 1, can be further simplified to: 35 36 3 εG 37 1 38 where θ = φ1 –φ2 is the phase difference, ∆ω = ω1 – ω2 the detuning, εG(θ)= ε21H21(θ)-ε12H12(-θ) 39 the combined interaction term with ε being the interaction strength and G(θ) the mutual PRC 2 40 and ƞ= ƞ1 –ƞ2 the phase noise with η ~ N(0, 2σ ). 41 42 Analytical derivation of phase-locking and mean phase difference 43 Equation 3 is a stochastic differential equation (a Langevin equation) and was solved as 44 described in 1. Equation 3 can be rewritten in the form of a Fokker-Planck equation that has been 45 developed to give an analytical solution for the evolution of a probability distribution P of a 46 particle influenced by a drag force (first term on the right side of the equation) and a random 47 Gaussian noise process (second term). The drag force is here the combined systematic force of 48 detuning ∆ω and the interaction function εG(θ): 49 ∆ 50 4 51 52 The stationary (time-independent) solution of the Fokker-Planck equation which is: ′ 1 ′ 53 5 exp 54 ∆ ′ 1 ′ 55 6 exp 56 where C is a normalization constant defined by 1. V(θ) represents the influence 57 of systematic force as a function of phase difference. is the phase difference probability 58 distribution and describes how likely a particular phase difference is to occur. A uniform 59 distribution means that every phase difference is equally likely and the oscillator are hence 60 asynchronous. If the distribution approximates a delta distribution (meaning only one phase 61 difference has non-zero probability), then the oscillators are completely synchronized. All other 62 distributions in between signify intermittent (partial) synchronization (also called cycle slipping 63 or phase walk-through, 1,2). To quantifying the narrowness of the distribution, we use the phase- 64 locking value (the mean resultant vector length,3) defined here as: 65 7 | | 66 Further, we were also interested in the mean phase difference, also described as the preferred 67 phase difference, defined here as: 68 8 arg 2 69 A phase difference between oscillators in neural networks implies spike timing differences. It has 70 been shown that spike-timing is an important characteristic in addition to spike synchrony 4–9. 71 72 General behaviour of the model in term of PLV and mean phase difference 73 In Fig.S1 the model’s behaviour is illustrated as a function of detuning ∆ω and interaction 74 strength ε. To understand how the PLV and the mean phase difference change, it is illustrative to 75 consider the noise-free case first. In the noise-free case one needs to solve the equation 3 for 76 zero-points (root or equilibrium points), meaning that the time derivative of the phase difference 77 is zero, ( 0, i.e. zero frequency difference). Let us assume here for illustration that G(θ) is a 78 sine function, as is for example used for the Kuramoto model 10,11. It then becomes (also called 79 the Adler equation, 1): 80 9 0 ∆ 81 The zero-point or the equilibrium is stable if the derivative of the PRC is negative (stable if 82 < 0). The sin function has two zero crossing (at 0 and π with Δω =0), but only one is stable. 83 In the case ε = 0 (no interaction), only the condition of ∆ω = 0 (no detuning) can lead to 84 synchrony. If ε = 0 and ∆ω ≠ 0 (no interaction with detuning), the oscillators will be 85 asynchronous and there will be linear phase precession with the speed determined by the ∆ω 86 value. If ε > 0 (there is interaction), then an equilibrium can only exist if the detuning is not 87 stronger than the interaction strength, denoted by |∆ω|<= ε. This is because for equilibrium the 88 detuning needs to be counter-balanced by the interaction. In the case of |∆ω|> ε (detuning 89 exceeding interaction strength), the interaction cannot counterbalance the detuning and the 90 oscillators start to phase precess. 91 If the oscillators do phase precess but are coupled (ε>0), the phase-precession is non-linear. The 92 precession rate is determined by the detuning (∆ω), the modulation shape (G(θ)), and the 93 modulation amplitude (ε). In this case, the oscillators are in the intermittent synchronization 94 regime. The transition point ∆ω = ε defines the border between high PLV values and low PLV 95 values and defines the borders of the so-called Arnold tongue, the synchronization region in the 96 ∆ω- ε parameter space 1,2. The larger the ε, the larger the detuning can be for oscillators to still 97 synchronize. This leads to a triangle shape of the synchronization region in the ∆ω- ε space as 98 illustrated in Fig.S1A. The mean (preferred) phase difference can also be derived from equation 99 6. For reaching equilibrium the detuning ∆ω and interaction term (εsin(θ)) need to be 100 counterbalanced. For different detuning ∆ω values the counterbalance will depend on different 101 phase difference values of the interaction term εsin(θ). It can be represented graphically 102 (Fig.S1B,2) as the intersection of a horizontal line (∆ω) with the graph of εsin(θ). Notice that the 103 stronger the interaction strength ε is, the smaller the slope of the detuning-to-phase difference 104 translation becomes. Outside of the Arnold tongue (∆ω> ε), there will be phase precession yet 105 with a preference (minimal precession rate) for a phase difference where the PRC is min or max. 106 Phase noise has important effects on the synchronization behavior1. Strictly speaking, the 107 condition of complete synchrony for noisy oscillators does not exist as θ’ cannot remain stable 108 all the time (there will always be small fluctuations). If ∆ω+ƞ < ε, meaning that in spite of 109 frequency variability detuning can be counterbalanced by interaction strength, then oscillators 110 remain close to the equilibrium point and they have high phase-locking. Yet, when noise is 3 111 larger, it is likely that the interaction term cannot counterbalance all the time and the oscillators 112 will phase precess from time to time. However, also in this case, the transition point ∆ω = ε still 113 determines the Arnold tongue border and outlines the transition between the regions of high 114 phase-locking and regions of low phase-locking, but the transition is smoother1,2. In Fig.S1 we 115 show the Arnold tongue in terms of PLV and mean phase-difference for coupled oscillators 116 having moderate levels of white noise (σ=15Hz). 117 118 Biophysical modeling of coupled gamma-generating neural networks 119 To demonstrate that the results from the phase-oscillator equations are generalizable to more 120 biophysically realistic neuronal network oscillations12, we simulated two coupled excitatory- 121 inhibitory spiking neural networks generating pyramidal-interneuron gamma (PING,13) 122 oscillations.
Recommended publications
  • Cognitive Map Formation Through Sequence Encoding by Theta Phase Precession
    LETTER Communicated by A. David Redish Cognitive Map Formation Through Sequence Encoding by Theta Phase Precession Hiroaki Wagatsuma [email protected] Laboratory for Dynamics of Emergent Intelligence, RIKEN BSI, Wako-shi, Saitama 351-0198, Japan Yoko Yamaguchi [email protected] Laboratory for Dynamics of Emergent Intelligence, RIKEN BSI, Wako-shi, Saitama 351- 0198, Japan; College of Science and Engineering, Tokyo Denki University, Hatoyama, Saitama 350-0394, Japan; and CREST, Japan Science and Technology Corporation The rodent hippocampus has been thought to represent the spatial en- vironment as a cognitive map. The associative connections in the hip- pocampus imply that a neural entity represents the map as a geometrical network of hippocampal cells in terms of a chart. According to recent experimental observations, the cells fire successively relative to the theta oscillation of the local field potential, called theta phase precession, when the animal is running. This observation suggests the learning of temporal sequences with asymmetric connections in the hippocampus, but it also gives rather inconsistent implications on the formation of the chart that should consist of symmetric connections for space coding. In this study, we hypothesize that the chart is generated with theta phase coding through the integration of asymmetric connections. Our computer experiments use a hippocampal network model to demonstrate that a geometrical network is formed through running experiences in a few minutes. Asymmetric connections are found to remain and dis- tribute heterogeneously in the network. The obtained network exhibits the spatial localization of activities at each instance as the chart does and their propagation that represents behavioral motions with multidirec- tional properties.
    [Show full text]
  • Position–Theta-Phase Model of Hippocampal Place Cell Activity Applied to Quantification of Running Speed Modulation of Firing Rate
    Position–theta-phase model of hippocampal place cell activity applied to quantification of running speed modulation of firing rate Kathryn McClaina,b,1, David Tingleyb, David J. Heegera,c,1, and György Buzsákia,b,1 aCenter for Neural Science, New York University, New York, NY 10003; bNeuroscience Institute, New York University, New York, NY 10016; and cDepartment of Psychology, New York University, New York, NY 10003 Edited by Terrence J. Sejnowski, Salk Institute for Biological Studies, La Jolla, CA, and approved November 18, 2019 (received for review July 24, 2019) Spiking activity of place cells in the hippocampus encodes the Another challenge in understanding this system is the dynamic animal’s position as it moves through an environment. Within a interaction between the rate code and phase code. As these codes cell’s place field, both the firing rate and the phase of spiking in combine, different formats of information are conveyed simulta- the local theta oscillation contain spatial information. We propose a neously in place cell spiking. The interaction can produce unin- – position theta-phase (PTP) model that captures the simultaneous tuitive, although entirely predictable, results in traditional analyses expression of the firing-rate code and theta-phase code in place cell of place cell activity. These practical challenges have hindered the spiking. This model parametrically characterizes place fields to com- pare across cells, time, and conditions; generates realistic place cell effort to explain variability in hippocampal firing rates. A com- simulation data; and conceptualizes a framework for principled hy- putational tool is needed that accounts for the well-established pothesis testing to identify additional features of place cell activity.
    [Show full text]
  • Neural Mechanism to Simulate a Scale-Invariant Future Timeline
    Neural mechanism to simulate a scale-invariant future Karthik H. Shankar, Inder Singh, and Marc W. Howard1 1Center for Memory and Brain, Initiative for the Physics and Mathematics of Neural Systems, Boston University Predicting future events, and their order, is important for efficient planning. We propose a neural mechanism to non-destructively translate the current state of memory into the future, so as to construct an ordered set of future predictions. This framework applies equally well to translations in time or in one-dimensional position. In a two-layer memory network that encodes the Laplace transform of the external input in real time, translation can be accomplished by modulating the weights between the layers. We propose that within each cycle of hippocampal theta oscillations, the memory state is swept through a range of translations to yield an ordered set of future predictions. We operationalize several neurobiological findings into phenomenological equations constraining translation. Combined with constraints based on physical principles requiring scale-invariance and coherence in translation across memory nodes, the proposition results in Weber-Fechner spacing for the representation of both past (memory) and future (prediction) timelines. The resulting expressions are consistent with findings from phase precession experiments in different regions of the hippocampus and reward systems in the ventral striatum. The model makes several experimental predictions that can be tested with existing technology. I. INTRODUCTION The brain encodes externally observed stimuli in real time and represents information about the current spatial location and temporal history of recent events as activity distributed over neural networks. Although we are physi- cally localized in space and time, it is often useful for us to make decisions based on non-local events, by anticipating events to occur at distant future and remote locations.
    [Show full text]
  • A Temporal Mechansism for Generating the Phase Precession Of
    Journal of Computational Neuroscience 9, 5–30, 2000 °c 2000 Kluwer Academic Publishers. Manufactured in The Netherlands. A Temporal Mechanism for Generating the Phase Precession of Hippocampal Place Cells AMITABHA BOSE, VICTORIA BOOTH AND MICHAEL RECCE Department of Mathematical Sciences, Center for Applied Mathematics and Statistics, New Jersey Institute of Technology, Newark, NJ, 07102-1982 [email protected] Received December 1, 1998; Revised June 21, 1999; Accepted June 25, 1999 Action Editor: John Rinzel Abstract. The phase relationship between the activity of hippocampal place cells and the hippocampal theta rhythm systematically precesses as the animal runs through the region in an environment called the place field of the cell. We present a minimal biophysical model of the phase precession of place cells in region CA3 of the hippocampus. The model describes the dynamics of two coupled point neurons—namely, a pyramidal cell and an interneuron, the latter of which is driven by a pacemaker input. Outside of the place field, the network displays a stable, background firing pattern that is locked to the theta rhythm. The pacemaker input drives the interneuron, which in turn activates the pyramidal cell. A single stimulus to the pyramidal cell from the dentate gyrus, simulating entrance into the place field, reorganizes the functional roles of the cells in the network for a number of cycles of the theta rhythm. In the reorganized network, the pyramidal cell drives the interneuron at a higher frequency than the theta frequency, thus causing a systematic precession relative to the theta input. The frequency of the pyramidal cell can vary to account for changes in the animal’s running speed.
    [Show full text]
  • Theta Phase Precession Beyond the Hippocampus
    Theta phase precession beyond the hippocampus Authors: Sushant Malhotra1,2y, Robert W. A. Cross1y, Matthijs A. A. van der Meer1,3* 1Department of Biology, University of Waterloo, Ontario, Canada 2Systems Design Engineering, University of Waterloo 3Centre for Theoretical Neuroscience, University of Waterloo yThese authors contributed equally. *Correspondence should be addressed to MvdM, Department of Biology, University of Waterloo, 200 Uni- versity Ave W, Waterloo, ON N2L 3G1, Canada. E-mail: [email protected]. Running title: Phase precession beyond the hippocampus 1 Abstract The spike timing of spatially tuned cells throughout the rodent hippocampal formation displays a strikingly robust and precise organization. In individual place cells, spikes precess relative to the theta local field po- tential (6-10 Hz) as an animal traverses a place field. At the population level, theta cycles shape repeated, compressed place cell sequences that correspond to coherent paths. The theta phase precession phenomenon has not only afforded insights into how multiple processing elements in the hippocampal formation inter- act; it is also believed to facilitate hippocampal contributions to rapid learning, navigation, and lookahead. However, theta phase precession is not unique to the hippocampus, suggesting that insights derived from the hippocampal phase precession could elucidate processing in other structures. In this review we consider the implications of extrahippocampal phase precession in terms of mechanisms and functional relevance. We focus on phase precession in the ventral striatum, a prominent output structure of the hippocampus in which phase precession systematically appears in the firing of reward-anticipatory “ramp” neurons. We out- line how ventral striatal phase precession can advance our understanding of behaviors thought to depend on interactions between the hippocampus and the ventral striatum, such as conditioned place preference and context-dependent reinstatement.
    [Show full text]
  • Dendritic Spiking Accounts for Rate and Phase Coding in a Biophysicalmodelof a Hippocampal Place Cell
    ARTICLE IN PRESS Neurocomputing 65–66 (2005) 331–341 www.elsevier.com/locate/neucom Dendritic spiking accounts for rate and phase coding in a biophysicalmodelof a hippocampal place cell Zso´ fia Huhna,Ã,Ma´ te´ Lengyela,b, Gergo+ Orba´ na,Pe´ ter E´ rdia,c aDepartment of Biophysics, KFKI Research Institute for Particle and Nuclear Physics, Hungarian Academy of Sciences, 29-33 Konkoly Thege M. u´t, Budapest H-1121, Hungary bGatsby Computational Neuroscience Unit, University College London, 17 Queen Square, London WC1N 3AR, UK cCenter for Complex Systems Studies, Kalamazoo College, 1200 Academy street, Kalamazoo, MI 49006-3295, USA Available online 15 December 2004 Abstract Hippocampal place cells provide prototypical examples of neurons jointly firing phase and rate-coded spike trains. We propose a biophysicalmechanism accounting for the generation of place cell firing at the single neuron level. An interplay between external theta-modulated excitation impinging the dendrite and intrinsic dendritic spiking as well as between frequency- modulated dendritic spiking and somatic membrane potential oscillations was a key element of the model. Through these interactions robust phase and rate-coded firing emerged in the model place cell, reproducing salient experimentally observed properties of place cell firing. r 2004 Elsevier B.V. All rights reserved. Keywords: Phase precession; Tuning curve; Active dendrite ÃCorresponding author. Fax: +36 1 3922222. E-mail address: zsofi@rmki.kfki.hu (Z. Huhn). 0925-2312/$ - see front matter r 2004 Elsevier B.V. All rights reserved. doi:10.1016/j.neucom.2004.10.026 ARTICLE IN PRESS 332 Z. Huhn et al. / Neurocomputing 65–66 (2005) 331–341 1.
    [Show full text]
  • Phase Precession of Medial Prefrontal Cortical Activity Relative to the Hippocampal Theta Rhythm
    HIPPOCAMPUS 15:867–873 (2005) Phase Precession of Medial Prefrontal Cortical Activity Relative to the Hippocampal Theta Rhythm Matthew W. Jones and Matthew A. Wilson* ABSTRACT: Theta phase-locking and phase precession are two related spike-timing and the phase of the concurrent theta phenomena reflecting coordination of hippocampal place cell firing cycle: place cell spikes are ‘‘phase-locked’’ to theta with the local, ongoing theta rhythm. The mechanisms and functions of both the phenomena remain unclear, though the robust correlation (Buzsaki and Eidelberg, 1983). Thus, the spikes of a between firing phase and location of the animal has lead to the sugges- given neuron are not distributed randomly across the tion that this phase relationship constitutes a temporal code for spatial theta cycle. Rather, they consistently tend to occur information. Recent work has described theta phase-locking in the rat during a restricted phase of the oscillation. Further- medial prefrontal cortex (mPFC), a structure with direct anatomical and more, in CA1, CA3, and the dentate gyrus, phase- functional links to the hippocampus. Here, we describe an initial char- locking to the local theta rhythm is accompanied by acterization of phase precession in the mPFC relative to the CA1 theta rhythm. mPFC phase precession was most robust during behavioral the related phenomenon of phase precession: each epochs known to be associated with enhanced theta-frequency coordi- neuron tends to fire on progressively earlier phases of nation of CA1 and mPFC activities. Precession was coherent across the the concurrent theta cycle as an animal moves through mPFC population, with multiple neurons precessing in parallel as a func- its place field (O’Keefe and Recce, 1993).
    [Show full text]
  • Phase-‐Of-‐Firing Code
    PHASE-OF-FIRING CODE Anna Cattania, Gaute T. Einevollb,c, Stefano Panzeria* a Center for Neuroscience and Cognitive Systems @UniTn, Istituto Italiano di Tecnologia, 38068 Rovereto, Italy b Department of Mathematical Sciences and Technology, Norwegian University of Life Sciences, 1432 Ås, Norway c Department of Physics, University of Oslo, 0316 Oslo, Norway * Corresponding author: [email protected] DEFINITION The phase-of-firing code is a neural coding scheme whereby neurons encode information using the time at which they fire spikes within a cycle of the ongoing oscillatory pattern of network activity. This coding scheme may allow neurons to use their temporal pattern of spikes to encode information that is not encoded in their firing rate. DETAILED DESCRIPTION In this entry, we define phase-of-firing coding, and we outline its properties. We first review early results that suggested the importance of temporal coding, and then we elucidate the insights concerning the central role of phase-of-firing coding in the studies of encoding of spatial variables in the hippocampus and of natural visual stimuli in primary visual and auditory cortices. The phase-of-firing code is a temporal coding scheme that uses the timing of spikes with respect to network fluctuations to encode information about some variables of interest (such as the identity of an external sensory stimulus) that is not encoded by the firing rates of these neurons (Perkel and Bullock, 1968). Figure 1 shows a schematic representation of both a phase-of-firing code and of a firing-rate code (i.e. a code that uses only the total number of spikes fired by the neuron to transmit information).
    [Show full text]
  • The Information Content of Local Field Potentials: Experiments and Models
    The information content of Local Field Potentials: experiments and models Alberto Mazzoni a, Nikos K. Logothetis b, c, Stefano Panzeri a,d * a Center for Neuroscience and Cognitive Systems, Istituto Italiano di Tecnologia, Via Bettini 31, 38068 Rovereto (TN), Italy b Max Planck Institute for Biological Cybernetics, Spemannstrasse 38, 72076 Tübingen, Germany c Division of Imaging Science and Biomedical Engineering, University of Manchester, Manchester M13 9PT, United Kingdom d Institute of Neuroscience and Psychology, University of Glasgow, Glasgow G12 8QB, United Kingdom * Corresponding author: S. Panzeri, email: [email protected] Keywords: Vision, recurrent networks, integrate and fire, LFP modelling, neural code, spike timing, natural movies, mutual information, oscillations, gamma band 1 1. Local Field Potential dynamics offer insights into the function of neural circuits The Local Field Potential (LFP) is a massed neural signal obtained by low pass‐filtering (usually with a cutoff low‐pass frequency in the range of 100‐300 Hz) of the extracellular electrical potential recorded with intracranial electrodes. LFPs have been neglected for a few decades because in‐vivo neurophysiological research focused mostly on isolating action potentials from individual neurons, but the last decade has witnessed a renewed interested in the use of LFPs for studying cortical function, with a large amount of recent empirical and theoretical neurophysiological studies using LFPs to investigate the dynamics and the function of neural circuits under different conditions. There are many reasons why the use of LFP signals has become popular over the last 10 years. Perhaps the most important reason is that LFPs and their different band‐limited components (known e.g.
    [Show full text]
  • Spatial Theta Cells in Competitive Burst Synchronization Networks: Reference Frames from Phase Codes
    bioRxiv preprint doi: https://doi.org/10.1101/211458; this version posted October 30, 2017. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. RUNNING TITLE: Spatial Theta Cells Spatial Theta Cells in Competitive Burst Synchronization Networks: Reference Frames from Phase Codes Joseph D. Monaco1∗, Hugh T. Blair2, Kechen Zhang1;3 1Biomedical Engineering Department, Johns Hopkins University School of Medicine, Baltimore, MD, USA; 2Psychology Department, UCLA, Los Angeles, CA, USA 3Department of Neuroscience, Johns Hopkins University, Baltimore, MD, USA ∗Correspondence: Joseph Monaco Johns Hopkins University School of Medicine Biomedical Engineering Department 720 Rutland Ave Traylor 407 Baltimore, MD 21205, USA [email protected] bioRxiv preprint doi: https://doi.org/10.1101/211458; this version posted October 30, 2017. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. 1 Abstract 2 Spatial cells of the hippocampal formation are embedded in networks of theta cells. The septal 3 theta rhythm (6–10 Hz) organizes the spatial activity of place and grid cells in time, but it remains 4 unclear how spatial reference points organize the temporal activity of theta cells in space. We 5 study spatial theta cells in simulations and single-unit recordings from exploring rats to ask 6 whether temporal phase codes may anchor spatial representations to the outside world.
    [Show full text]
  • Grid Cells in Rat Entorhinal Cortex Encode Physical Space with Independent firing fields and Phase Precession at the Single-Trial Level
    Grid cells in rat entorhinal cortex encode physical space with independent firing fields and phase precession at the single-trial level Eric T. Reifensteina,b, Richard Kempterb, Susanne Schreiberb, Martin B. Stemmlera, and Andreas V. M. Herza,1 aDepartment of Biology II, Ludwig-Maximilians-Universität München, and Bernstein Center for Computational Neuroscience Munich, 82152 Planegg- Martinsried, Germany; and bInstitute for Theoretical Biology, Humboldt-Universität zu Berlin, and Bernstein Center for Computational Neuroscience Berlin, 10115 Berlin, Germany Edited* by John J. Hopfield, Princeton University, Princeton, NJ, and approved March 08, 2012 (received for review June 14, 2011) When a rat moves, grid cells in its entorhinal cortex become active to a particular theta phase of the local field potential (LFP). There is in multiple regions of the external world that form a hexagonal also another aspect of grid-cell activity that needs to be taken into lattice. As the animal traverses one such “firing field,” spikes tend account: From one firing field to the next one visited, the discharge to occur at successively earlier theta phases of the local field po- of a single grid cell may or may not be correlated. Hence it is an open tential. This phenomenon is called phase precession. Here, we question whether the nervous system can actually make use of phase show that spike phases provide 80% more spatial information coding in mEC, notwithstanding the trends seen in pooled data or at than spike counts and that they improve position estimates from the single-run level in hippocampal place cells (13). single neurons down to a few centimeters.
    [Show full text]
  • Position-Theta-Phase Model of Hippocampal Place Cell Activity Applied to Quantification of Running Speed Modulation of firing Rate
    bioRxiv preprint doi: https://doi.org/10.1101/714105; this version posted July 24, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Position-theta-phase model of hippocampal place cell activity applied to quantification of running speed modulation of firing rate Kathryn McClain1,2, David Tingley1, David Heeger*1,3, György Buzsáki*1,2 1. Center for Neural Science, New York University, 4 Washington Pl, New York, NY 10003, USA 2. NYU Neuroscience Institute, 450 East 29th Street, New York, NY 10016, USA. 3. Department of Psychology, New York University, 4 Washington Pl, New York, NY, 10003, United States * Corresponding authors: [email protected], [email protected] bioRxiv preprint doi: https://doi.org/10.1101/714105; this version posted July 24, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Abstract Spiking activity of place cells in the hippocampus encodes the animal’s position as it moves through an environment. Within a cell’s place field, both the firing rate and the phase of spiking in the local theta oscillation contain spatial information. We propose a position-theta-phase (PTP) model that captures the simultaneous expression of the firing-rate code and theta-phase code in place cell spiking.
    [Show full text]