Replay as wavefronts and theta sequences as bump oscillations in a grid cell attractor network

  1. Louis Kang  Is a corresponding author
  2. Michael R DeWeese
  1. University of California, Berkeley, United States


Grid cells fire in sequences that represent rapid trajectories in space. During locomotion, theta sequences encode sweeps in position starting slightly behind the animal and ending ahead of it. During quiescence and slow wave sleep, bouts of synchronized activity represent long trajectories called replays, which are well-established in place cells and have been recently reported in grid cells. Theta sequences and replay are hypothesized to facilitate many cognitive functions, but their underlying mechanisms are unknown. One mechanism proposed for grid cell formation is the continuous attractor network. We demonstrate that this established architecture naturally produces theta sequences and replay as distinct consequences of modulating external input. Driving inhibitory interneurons at the theta frequency causes attractor bumps to oscillate in speed and size, which gives rise to theta sequences and phase precession, respectively. Decreasing input drive to all neurons produces traveling wavefronts of activity that are decoded as replays.

Data availability

Source code for the simulations have been included as supporting files.

Article and author information

Author details

  1. Louis Kang

    Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, United States
    For correspondence
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5702-2740
  2. Michael R DeWeese

    Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, United States
    Competing interests
    The authors declare that no competing interests exist.


Adolph C. and Mary Sprague Miller Institute for Basic Research in Science, University of California Berkeley (Postdoctoral fellowship)

  • Louis Kang

Army Research Office (W911NF-13-1-0390)

  • Michael R DeWeese

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Sachin Deshmukh, Indian Institute of Science Bangalore, India

Publication history

  1. Received: March 1, 2019
  2. Accepted: November 15, 2019
  3. Accepted Manuscript published: November 18, 2019 (version 1)
  4. Version of Record published: December 9, 2019 (version 2)


© 2019, Kang & DeWeese

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.


  • 1,651
    Page views
  • 228
  • 7

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Louis Kang
  2. Michael R DeWeese
Replay as wavefronts and theta sequences as bump oscillations in a grid cell attractor network
eLife 8:e46351.

Further reading

    1. Neuroscience
    Benjamin D Pedigo, Mike Powell ... Joshua T Vogelstein
    Research Article

    Comparing connectomes can help explain how neural connectivity is related to genetics, disease, development, learning, and behavior. However, making statistical inferences about the significance and nature of differences between two networks is an open problem, and such analysis has not been extensively applied to nanoscale connectomes. Here, we investigate this problem via a case study on the bilateral symmetry of a larval Drosophila brain connectome. We translate notions of'bilateral symmetry' to generative models of the network structure of the left and right hemispheres, allowing us to test and refine our understanding of symmetry. We find significant differences in connection probabilities both across the entire left and right networks and between specific cell types. By rescaling connection probabilities or removing certain edges based on weight, we also present adjusted definitions of bilateral symmetry exhibited by this connectome. This work shows how statistical inferences from networks can inform the study of connectomes, facilitating future comparisons of neural structures.

    1. Neuroscience
    Sinisa Pajevic, Dietmar Plenz ... R Douglas Fields
    Research Article

    Temporal synchrony of signals arriving from different neurons or brain regions is essential for proper neural processing. Nevertheless, it is not well understood how such synchrony is achieved and maintained in a complex network of time-delayed neural interactions. Myelin plasticity, accomplished by oligodendrocytes (OLs), has been suggested as an efficient mechanism for controlling timing in brain communications through adaptive changes of axonal conduction velocity and consequently conduction time delays, or latencies; however, local rules and feedback mechanisms that OLs use to achieve synchronization are not known. We propose a mathematical model of oligodendrocyte-mediated myelin plasticity (OMP) in which OLs play an active role in providing such feedback. This is achieved without using arrival times at the synapse or modulatory signaling from astrocytes; instead, it relies on the presence of global and transient OL responses to local action potentials in the axons they myelinate. While inspired by OL morphology, we provide the theoretical underpinnings that motivated the model and explore its performance for a wide range of its parameters. Our results indicate that when the characteristic time of OL’s transient intracellular responses to neural spikes is between 10 and 40 ms and the firing rates in individual axons are relatively low (⪅ 10 Hz), the OMP model efficiently synchronizes correlated and time-locked signals while latencies in axons carrying independent signals are unaffected. This suggests a novel form of selective synchronization in the CNS in which oligodendrocytes play an active role by modulating the conduction delays of correlated spike trains as they traverse to their targets.