Toroidal topology of population activity in grid cells.

The medial entorhinal cortex is part of a neural system for mapping the position of an individual within a physical environment1. Grid cells, a key component of this system, fire in a characteristic hexagonal pattern of locations2, and are organized in modules3 that collectively form a population code for the animal's allocentric position1. The invariance of the correlation structure of this population code across environments4,5 and behavioural states6,7, independent of specific sensory inputs, has pointed to intrinsic, recurrently connected continuous attractor networks (CANs) as a possible substrate of the grid pattern1,8-11. However, whether grid cell networks show continuous attractor dynamics, and how they interface with inputs from the environment, has remained unclear owing to the small samples of cells obtained so far. Here, using simultaneous recordings from many hundreds of grid cells and subsequent topological data analysis, we show that the joint activity of grid cells from an individual module resides on a toroidal manifold, as expected in a two-dimensional CAN. Positions on the torus correspond to positions of the moving animal in the environment. Individual cells are preferentially active at singular positions on the torus. Their positions are maintained between environments and from wakefulness to sleep, as predicted by CAN models for grid cells but not by alternative feedforward models12. This demonstration of network dynamics on a toroidal manifold provides a population-level visualization of CAN dynamics in grid cells.

Decoding of Neural Data Using Cohomological Feature Extraction.

We introduce a novel data-driven approach to discover and decode features in the neural code coming from large population neural recordings with minimal assumptions, using cohomological feature extraction. We apply our approach to neural recordings of mice moving freely in a box, where we find a circular feature. We then observe that the decoded value corresponds well to the head direction of the mouse. Thus, we capture head direction cells and decode the head direction from the neural population activity without having to process the mouse's behavior. Interestingly, the decoded values convey more information about the neural activity than the tracked head direction does, with differences that have some spatial organization. Finally, we note that the residual population activity, after the head direction has been accounted for, retains some low-dimensional structure that is correlated with the speed of the mouse.

Using persistent homology to reveal hidden covariates in systems governed by the kinetic Ising model.

We propose a method, based on persistent homology, to uncover topological properties of a priori unknown covariates in a system governed by the kinetic Ising model with time-varying external fields. As its starting point the method takes observations of the system under study, a list of suspected or known covariates, and observations of those covariates. We infer away the contributions of the suspected or known covariates, after which persistent homology reveals topological information about unknown remaining covariates. Our motivating example system is the activity of neurons tuned to the covariates physical position and head direction, but the method is far more general.

Synthesising Topological Links.

We discuss the chemical synthesis of topological links, in particular higher order links which have the Brunnian property (namely that removal of any one component unlinks the entire system). Furthermore, we suggest how to obtain both two dimensional and three dimensional objects (surfaces and solids, respectively) which also have this Brunnian property.

Higher order Boolean networks as models of cell state dynamics.

The regulation of the cell state is a complex process involving several components. These complex dynamics can be modeled using Boolean networks, allowing us to explain the existence of different cell states and the transition between them. Boolean models have been introduced both as specific examples and as ensemble or distribution network models. However, current ensemble Boolean network models do not make a systematic distinction between different cell components such as epigenetic factors, gene and transcription factors. Consequently, we still do not understand their relative contributions in controlling the cell fate. In this work we introduce and study higher order Boolean networks, which feature an explicit distinction between the different cell components and the types of interactions between them. We show that the stability of the cell state dynamics can be determined solving the eigenvalue problem of a matrix representing the regulatory interactions and their strengths. The qualitative analysis of this problem indicates that, in addition to the classification into stable and chaotic regimes, the cell state can be simple or complex depending on whether it can be deduced from the independent study of its components or not. Finally, we illustrate how the model can be expanded considering higher levels and higher order dynamics.