Through synapses to spatial memory maps: a topological model
Abstract
Abstract. Various neurophysiological and cognitive functions are based on transferring information between spiking neurons via a complex system of synaptic connections. In particular, the capacity of presynaptic inputs to influence the postsynaptic outputs—the efficacy of the synapses—plays a principal role in all aspects of hippocampal neurophysiology. However, a direct link between the information processed at the level of individual synapses and the animal’s ability to form memories at the organismal level has not yet been fully understood. Here, we investigate the effect of synaptic transmission probabilities on the ability of the hippocampal place cell ensembles to produce a cognitive map of the environment. Using methods from algebraic topology, we find that weakening synaptic connections increase spatial learning times, produce topological defects in the largescale representation of the ambient space and restrict the range of parameters for which place cell ensembles are capable of producing a map with correct topological structure. On the other hand, the results indicate a possibility of compensatory phenomena, namely that spatial learning deficiencies may be mitigated through enhancement of neuronal activity.
I Introduction
The locationspecific spiking activity of the hippocampal neurons, known as place cells OKeefe (), gives rise to an internalized representation of space—a cognitive map. Each place cell fires a series of action potentials in specific spatial region—its place field, so that the ensemble of such cells produces a “map” of the environment in which they are active (Fig. 1A). By construction, such a map defines the temporal order in which place cells fire as the animal explores the environment, and therefore it can be viewed as a geometric representation of the spatial memory framework encoded by the hippocampus Moser (); Schmidt ().
The exact nature of this framework is currently actively studied both computationally and experimentally. For example, it was demonstrated that if the shape of the environment gradually changes, then the place field map deforms in a way that preserves mutual overlaps, adjacencies, containments, etc., between the place fields Gothard (); Leutgeb (); Wills (); Touretzky (); eLife (). This observation implies that the sequence in which the place cells fire during animal’s navigation remains invariant throughout the reshaping of the arena and suggests that the place cells do not represent precise geometric information, but a set of qualitative connections between portions of the environment—a topological map Poucet (); Alvernhe1 (); eLife ().
From the computational perspective, the topological nature of the cognitive map suggests that the information transmitted via place cell spiking should be amenable to topological analyses. In our previous work PLoS (); Arai (); Basso (); Hoffman (); CAs (), we developed a topological model that allows tracing how the information provided by the individual place cells may combine into a largescale topological map of the navigated space and quantifying the contributions of different neurophysiological parameters. However, previous studies did not include a key physiological aspect—the contribution of the synaptic connections into the processes of assembling the map. Below we will use the topological approach to model how synaptic imperfections can affect the topological structure of the cognitive map, its dynamics and its stability.
The paper is organized as follows. First, we outline the basic ideas and the key concepts used in the topological model—simplexes, simplicial complexes, topological loops, Betti numbers, etc., and explain how these concepts can be applied for describing hippocampal physiology. Second, we outline the parameters of synaptic connectivity and the constructions used to incorporate these parameters into the model. The analyses of the outcomes is given in the Results section and their implications are outlined in the Discussion.
Ii The Model
Topological description of the place cell spiking patterns. It is generally believed that the information encoded by the place cell network is represented by the connectivity between the place fields. A specific link is suggested by the classical AlexandrovČech’s theorem of Algebraic Topology asserts that the pattern of overlaps between regions that cover a space does, in fact, capture its topological structure Alexandroff (); Cech (). The implementation of this theorem is based on constructing the socalled “nerve simplicial complex” , whose vertexes correspond to the individual domains of the cover: onedimensional () links—to their pairwise overlaps, twodimensional () facets—to the triple overlaps and so forth (Fig. 1B). In other words, each order overlap between the place fields is schematically represented by an dimensional simplex , so that the full set of such simplexes incorporates the connectivity structure of the entire place field map Ghrist (); Curto (); PLoS (). According to the AlexandrovČech’s theorem, this complex has the same “topological shape” as , i.e., the same number of pieces, gaps and holes Hatcher (); Alexandrov (), which provides a link between the place cells’ spiking pattern and the topology of ambient space PLoS (); Arai (); Basso (); Hoffman (); CAs (), exploited below.
In general, simplicial complexes provide a convenient framework for describing a wide scope of physiological phenomena. For example, the combinations of the place fields traversed during the rat’s moves correspond to a chain of simplexes that qualitatively represents the shape of the physical trajectory: a closed chain represents a closed physical route, a pair of topologically equivalent chains represent two similar physical paths and so forth Guger (); Brown1 (). The pool of such chains can be used to describe the topological shape of the entire complex—and hence of the corresponding environment. For example, the number of chains that can be deformed into the same vertex defines how many disconnected pieces has. The number of topologically inequivalent chains that contract to a closed sequence of links defines the number of distinct holes that prevent these chains from contracting to vertexes and so forth Hatcher (); Alexandrov (). In the following, we will refer to these two types of chains, counted up to topological equivalence, as to zerodimensional () and onedimensional () “topological loops” (a standard mathematical terminology), evaluate their numbers—in mathematical terms, zeroth and first Betti numbers, and , and use them to describe shapes of the simplicial complexes.
Learning dynamics. To describe how the animal “learns” the environment, one can follow how the nerve complex and its Betti numbers develop in time. In the beginning of exploration, the nerve complex represents connections between the place fields that the animal had time to visit. Such a complex is small and may contain gaps that do not necessarily correspond to physical holes or inaccessible spatial domains of the environment. As the animal continues to navigate, the nerve complex grows and acquires more details; as a result, its the spurious gaps and holes (topological noise) disappear, leaving behind a few persistent ones that represent stable topological information (Fig. 2). The minimal time, , required to recover the correct number of topological loops,
(1) 
can be used as a theoretical estimate of the time needed to learn path connectivity PLoS (). In the case of the environment illustrated on Fig. 1A, with the Betti numbers , , the nerve complex is expected to have the same “topological barcode”: , .
Temporal coactivity complex. From the physiological perspective, the arguments based on the analyses of place fields and trajectories provide only an indirect description of information processing in the brain. In reality, the hippocampus and the downstream brain regions do not have access to the shapes and the locations of the place fields or to other artificial geometric constructs used by experimentalists to visualize their data. Physiologically, the information is represented via neuronal spiking activity: if the animal enters a location where several place fields overlap, then there is a probability, modulated by the rat’s location, that the corresponding place cells will produce spike trains that overlap temporally. This pattern of coactivity signals to the downstream brain areas that the regions encoded by these place cells overlap. Thus, in order to describe the learning process in proper terms, one needs to construct a temporal analogue of the nerve complex based only on the spiking signals, which is, in fact, straightforward. Indeed, one can represent an active place cell, , by a vertex ; a pair of coactive place cells, and —by a bond between the vertices and ; a coactive triple of place cells, , and —by a three vertex simplex and so on Ghrist (); Curto (); PLoS (). This construction produces a timedependent “coactivity complex” —a temporal analogue of the nerve complex constructed above, whose dynamics can also be used to model topological learning, e.g., to compute the learning time from the spiking data, , and so forth PLoS ().
Cell assembly complex. The construction of a temporal complex can be refined to reflect more subtle physiological details, e.g., the functional organization of the hippocampal network. Studies of place cells’ spiking times point out that these neurons tend to fire in “assemblies”—functionally interconnected groups that are believed to synaptically drive a population of “readout” neurons in the downstream networks Harris1 (); Harris2 (); Jackson (); ONeill (); Syntax (). The latter are wired to integrate spiking inputs from their respective cell assemblies and actualize the connectivity relationships between the regions encoded by the corresponding place cells Syntax (); SchemaS ().
This structure can be represented by the cell assembly complex, —a temporal coactivity complex whose maximal simplexes represent cell assemblies, rather than arbitrary combinations of coactive place cells. A convenient implementation of this construction is based on the classical “cognitive graph” model, in which place cells are represented as vertexes of a graph , while the connections (functional or physiological) between pairs of coactive cells are represented by the links, of SchemaS (); Burgess (); Muller (). The place cell assemblies then correspond to fully interconnected subgraphs of , i.e., to its maximal cliques Hoffman (); CAs (). Since a clique , as a combinatorial object, can be viewed as a simplex span by the same sets of vertexes, the collection of cliques of the coactivity graph produces a socalled clique simplicial complex Jonsson (), which represents the population of place cell assemblies and may hence be viewed as a cell assembly complex (Fig. 3).
Phenomenological description of the synaptic parameters. In the previous studies, we demonstrated that such complexes can acquire a correct topological shape in a biologically plausible period of time, in both planar and in voluminous environments, provided that the simulated spiking parameters values fall into the biological range PLoS (); Arai (); Basso (); Hoffman (); CAs (). However, the organization and the dynamics of these complexes did not reflect the parameters of synaptic connectivity, e.g., the mechanisms of transferring, detecting and interpreting neuronal (co)activity in the hippocampus and in the downstream networks. To account for these components, the topological model requires a basic modification: a particular coactivity pattern should be incorporated into an effective coactivity complex not by the virtue of being merely produced, but by the virtue of being produced, transmitted and ultimately detected by a readout neuron. In other words, only detected activity of a place cell should be represented by a vertex ; a detected coactivity of two place cells, and —by a bond , a detected coactivity of three place cells, , and —by a simplex and so on (Fig. 1B). The resulting complex then constitutes a basic phenomenological model of a cognitive map assembled from the spiking inputs transmitted through imperfect synaptic connections.
Statistical approach. The mechanisms of spike generation, transmission and detection are probabilistic in nature. Transmitting action potentials requires producing a sufficient number of synaptic contacts in suitable locations of the postsynaptic neuron’s membrane, releasing proper amount of neurotransmitter at each synapse at suitable times, inducing the excitatory postsynaptic potential (EPSP) of required magnitudes, etc., all of which involve probabilistic mechanisms Soltani (); London (). In addition, there may appear flaws and glitches in the axons, synaptic clefts and in the structure of the postsynaptic membrane’s polarization. Thus, there exists a probability that a connection in a cell assembly will induce sufficient EPSP in the readout neuron’s membrane and a probability that the latter will spike upon receiving the inputs (Fig. 3C).
In principle, these values could be estimated from the synaptic configuration of each individual assembly, which, however, would present a tremendous computational challenge Branco (); Arleo (); Garrido (). In order to avoid such complications, we will assume a basic statistical approach. First, we will regard the probabilities , and as the prime parameters that describe the synaptic connections with the readout neuron. Second, we will view and as random variables, distributed according to a unimodal distribution, and were and are the modes (the characteristic values) and and define the corresponding variances. Third, we will disregard synaptic plasticity processes and assume that the distributions are stationary, i.e., that the modes and the variances are fixed. Fourth, we will assume that both variables are distributed lognormally, as suggested by experimental observations BuzMi1 (); Barbour (); Brunel (). We will also define the variances as functions of the modes, and , which will allow us to exclude nonbiological statistics and to study the topological properties of the emerging cognitive maps as functions of just two parameters, and .
Implementation. In order to isolate the effects of varying transition probabilities while keeping the temporal structure of the presynaptic spike trains “clamped,” we use the spiking data that was precomputed for the “ideal” synaptic connections (), and then screen out some of the spikes, to match each individual transmission probabilities and to simulate the readout neurons’ responses to the igniting cell assemblies with probabilities .
To evaluate the latter, we reasoned as follows. Since in our approach the cell assemblies are modeled as the cliques of the coactivity graph , i.e., as composite objects assembled from pairs of place cells, the probabilities of igniting the higher order place cell combinations can be computed from the pairwise coactivities. Indeed, if the spikes produced by the place cells and are transmitted to the readout neuron with the probabilities and respectively, then the corresponding pairwise coactivity occurs with the probability . The probability of a third order coactivity, e.g., the ignition of a clique is then defined by the probability of transmitting the coactive pairs , , and and detecting the result with the probability ; the probability of igniting the fourth order cliques is defined by the corresponding six coactive pairs and so forth.
With these assumptions, one can test how the spike transmission and detection probabilities affect the emergence of a spatial map, e.g., how synaptic depletion affects spatial learning, how the learning times and the topological structure of the cognitive map depend upon the strengths of synaptic connections between the place cells and the readout neurons, at what point spatial learning may fail, and so on.
Iii Results
Learning times. Lowering the characteristic probability of spike transmissions and the characteristic probability of the readout neurons’ responses produces an uneven delay in spatial learning times (Fig. 4A). If the spike transmission probability is high (typically ), then the small variations of do not inflict a strong impact on , i.e., the time required to learn the spatial map in a network with strong synaptic connections is nearly unaffected by occasional omissions of spikes. On the other hand, as lowers to a certain critical value , the learning times become high and, as drops below , the coactivity complex fails to produce the correct topological shape of the environment in finite time. For the intermediate values, the learning time increases at a power rate,
(2) 
where ranges between and for different values of . The effects produced by the diminishing probability of the postsynaptic neurons’ responses, , are qualitatively similar but weaker than the effects of lowering the spike transmission probability : the learning time shows a weak or no dependence for large (typically ), followed by the power divergence near the critical value,
(3) 
with a small power exponent (Fig. 4B). Lowering both and simultaneously leads to a combined, accelerated increase of the learning time (Fig. 4B).
An implication of this phenomenon is that, and , being independent characteristics of synaptic efficacy, can also compensate for each other’s alterations: the effect of decreasing can be counterbalanced by increasing and vice versa. Indeed, the dependencies (2) and (3) also define the changes of the learning time induced by small variations in the transmission probability,
(4) 
and by the variations of the postsynaptic neuron’s response probability,
(5) 
These relationships imply that the compensation of the changes of the learning time, , is achieved if
(6) 
Notice, that this dependence is independent and nonlinear: given a particular value of , the required compensatory change of depends on the initial values of both and .
Dynamics of the effective coactivity complex. The failures of the learning and memory capacity caused by deterioration of synapses are broadly discussed in the literature Selkoe (); Neves (); Mayford (). However, empirical observations provide only correlative links between these two scopes of phenomena. Indeed, the direct effects of the synaptic changes, e.g. the alterations of EPSP magnitudes, the spike transmission probabilities, the parameters of synaptic plasticity, etc., occur at cellular scale. It therefore remains unclear how such changes may accumulate at the network scale to control the net structure and the dynamics of the largescale memory framework at the organismal level. The topological model allows addressing these questions at a phenomenological level, in terms of the structure of the coactivity complex —its topological shape, its size, the dynamics of its topological loops and so forth, in response to the changes of synaptic parameters.
For example, one can evaluate the statistics of birth () and death () times of the topological loops in the coactivity complex. As shown on (Fig. 4C), the time when spurious loops begin to emerge depend only marginally on spike transmission probability. However, the spurious loops’ disappearance times are impacted much stronger: although shows only weak dependence at high , further suppression of the spike transmissions may double or triple the loops’ disappearance time. The contribution of the decreasing response probability is similar, but at a smaller scale: over the range , the learning time changes only by a few percent (Fig. 4D). Similar effects are indicated by the  and dependencies of the spurious loops’ lengths, which may grow significantly as a result of the diminishing spike transmission probability, but increase only by due to the lowering probability of the readout neuron’s responses (Fig. 4E).
Taken together, these results explicate the power growth of the learning times indicated by (2) and (3) and provide a simple intuitive explanation for the decelerated spatial learning and its eventual failure caused by the synaptic depletion: according to the model, lowering synaptic efficacy stabilizes spurious topological loops in the coactivity complex, making it harder to extract physical information from the transient noise.
Additional perspective on the mechanisms of the cognitive map’s deterioration is produced by analyzing the size of the coactivity complex and the number of the topological loops in it. As shown on Fig. 5A, the decay of causes rapid decay of the coactivity complex’s size: the number of its twodimensional simplexes (i.e., links in the coactivity graph, see below) drops as , where . Diminishing also shrinks the coactivity complex, but at a slower rate, with . However, despite the shrinking size of the coactivity complex, the number of spurious loops in it grows exponentially, , from a few dozen to a few hundred, accompanied by a weak increase (Fig. 5B). Similar effects are produced by the lowering detection probability, but again, at a much smaller scale: the number of loops, , increases by about while the does not change (Fig. 5C).
These outcomes indicate that, as a result of weakening synaptic connections, the spurious topological loops do not only stabilize but also proliferate, thus preventing the effective coactivity complex from capturing the correct topology of the ambient space. In physiological terms, the model predicts that weakening synapses produce large numbers of longerlasting topological defects in the cognitive map, which results in a rapid increase of the time required to learn the topology of the physical environment from poorly communicated spiking inputs.
Critical probabilities. As indicated above, if the synaptic efficacies are too weak, i.e., if either the spike transmission or the postsynaptic response probability drops below their respective critical values, then the effective coactivity complex may disintegrate into a few disconnected pieces and lose its physical shape—a single large piece with a hole in the middle ( and , Fig. 1C), may be replaced by a “spongy” configuration containing several smaller pieces with many holes Ambjorn (); Hamber (). Thus, the cognitive map may appear in two distinct states: for and the spurious topological defects can be separated from the topological signatures of the physical environment, whereas below the critical values, topological noise overwhelms physical information. The transition between these two states is accompanied an increased variability of the learning times (Fig. 4A) and by their power divergence caused by the exponential proliferation of the topological fluctuations in the coactivity complex. These effects suggest that, near and , the coactivity complex may experience a phaselike transition Franzosi1 (); Franzosi2 (); Vaccarino () from a regular state, in which spatial learning is effective to an irregular state, in which spatial learning fails.
Since in most of the studied cases, the critical synaptic transmission probability, , is easier to achieve than the critical probability of the readout neuron’s responses, , we studied the dependence of the former on the ensemble parameters, i.e., on the number of place cells in the ensemble, their mean firing rate and the mean place field size,
(7) 
The results shown on Fig. 6 reveal powerlaw dependencies: , , and a more complex dependence. Since the domain of these dependences covers the experimentally observed range of parameters, the results can be interpreted physiologically. First, if the ensemble firing rates are too low, or if the place fields are too meager, or the number of the active neurons is too small (the left ends of the dependencies shown on Fig. 6), then the corresponding place cell ensemble fails to learn the spatial map of the environment, even if the synaptic connections are nearly perfect (), which corresponds to the results discussed in PLoS (); Arai (); Basso (). As the mean firing rate and the number of active neurons increase, the critical probability steadily decreases, which implies that the synaptic depletion may be compensated by enhancing neuronal activity, as observed in experimental studies Palop (); Hartley (); Busche (). In contrast, the dependence saturates and even reverses its direction for overly large place fields. This, however, is a natural result since poor spatial specificity of the place cells’ spiking should prevent successful spatial leaning even for large PLoS (); Arai (); Basso ().
Electrophysiological studies show that only up to of spikes are transmitted between the neurons in CA1 slices, which is lower than the critical values discussed above BuzMi1 (); Csicsvari (); BuzMi2 (). However, the results shown on Fig. 6C imply that the experimentally observed values of can be achieved for larger values of , i.e., in larger place cell ensembles. Interpolated dependence indicates that the physiological values of , can be achieved for the ensembles of cells, which corresponds to the experimentally observed values Wilson1 (); Leutgeb1 ().
Learning region. One of the key characteristics of the place cell spiking activity produced by the topological model is the range of the spiking parameters, for which the coactivity complex can assume a correct topological shape in a biologically feasible period. Geometrically, this set of parameters forms a domain in the parameter space that we refer to as the learning region, PLoS (). The shape and the size of the learning region varies with the geometric complexity of the environment and the difficulty of the task: the simpler is the environment and easier the task, the larger is , i.e., the wider the range of physiological values that permits learning a map of that space Nithianantharajah (); Hernan (). On the other hand, a larger implies a greater range within which the brain can compensate for physiological variation: if one parameter begins to drive the system outside the learning region, then successful spatial learning can still occur, provided that compensatory changes of other parameters can keep the neuronal ensemble inside . For example, a reduction of the number of active neurons can sometimes be compensated by adjusting the firing rate or the place field size in such a way as to bring their behavior back within the perimeter of the learning region.
Interpreting the parameters of a given place cell ensemble in the context of its placement within or relative to the learning region sheds light on the mechanism of memory failure caused by certain neurophysiological conditions, e.g., by the Alzheimer Disease Cacucci (); LaFerla (), or by aging Robitsek (); Wilson2 () or certain chemicals, e.g., ethanol White (); Matthews (), cannabinoids Robbe1 (); Robbe2 () or methamphetamines Kalechstein (); Silvers (), which appear to disrupt spatial learning by gradually shifting the parameters of spiking activity beyond the learning region. On the other hand, the performance of a place cell ensemble can improve by enhancing place cells spiking activity pharmacologically or by Deep Brain Stimulation Laxton (), or by modulating the hippocampal neural oscillations Shirvalkar (), as the model predicts PLoS (); Arai (); Basso ().
In contrast, diminishing spike transmission probability produces a qualitatively different effect: as shown on Fig. 7, it reduces the learning region from its original (largest) size at to its compete disappearance at the critical value . During this process, the time required to form the cognitive map of the environment progressively increases from a few minutes to over an hour (Fig. 7).
Physiologically, these results suggest that if the synaptic connections are too weak, then the system may fail to form a map not only because the parameters of neuronal firing are pushed beyond a certain “working range,” but also because that range itself may diminish or cease to exist. In particular, the fact that the learning region disappears if the transmission probability drops below the critical value, implies that the deterioration of memory capacity caused by synaptic failure may not be compensated by increasing the place field’s firing rates or by recruiting a larger population of active neurons, i.e., some neuropathological conditions may indeed be primarily “synaptic” in nature Selkoe ().
Deteriorating cognitive graph A simple alternative explanation of these results can be provided in terms of the place cell coactivity statistics. As pointed out in Section II, the collection of the unique pairs of the coactive place cells in a network with ideal synaptic connections () is represented by the coactivity graph . The imperfect synapses diminish the pool of the transmitted and the detected coactive pairs, which then corresponds to a smaller, effective coactivity graph . The corresponding set of higher order coactivities—the effective coactivity complex induced from is a subcomplex of the original coactivity complex, with potentially altered topological properties. The net results discussed above imply that, for high transmission probabilities, the effective coactivity complex retains the original topological shape of , but as diminishes, the effective complex shrinks, acquires multiple topological defects and eventually loses its correct shape, indicating a failure of spatial learning.
An illuminating perspective on the changing structure of the coactivity graph described above is provided by its Forman curvature—a combinatorial analogue of the standard differentialgeometric notion of curvature Forman (); Lewiner (). The Forman curvature is adopted for discrete, combinatorial structures, such as datasets, networks and graphs Weber1 (); Weber2 (); Weber3 (); Sreejith (), and can be flexibly defined in terms of an individual network’s characteristics—the “weights” of its vertexes and edges. Specifically, for an undirected edge with a weight connecting the vertexes and with the weights and it is defined as
(8) 
where the summation goes over the other edges and connecting to and . The curvature associated with a vertex, , equals to the mean curvature of the edges that meet at .
As discussed in Weber1 (); Weber2 (); Weber3 (); Sreejith (), the values and provide a measure of the divergence of information flow across the network, highlighting the most “important” edges and vertexes. Applying these ideas to the case of the coactivity graph, weighing its vertexes with the number of spikes produced by the corresponding place cells and its edges with correlation coefficients between the corresponding pairs of cells, reveals that the distribution of the resulting Forman curvatures follows the structure of the occupancy map (Fig. 8A,B). In other words, the most visited vertexes and edges appear as the most “curved” ones, controlling the flow of information in .
This quantification also allows a natural interpretation of the effective coactivity graph’s dynamics: as the spike transmission probability decreases, sheds the “least important” vertexes and links with low curvatures (Fig. 8C,D,E). Thus, as the synaptic efficacies weaken, the emerging effective coactivity graph reflects only the most persistently firing place cells and the highly correlated pairs of such cells, which can sustain the full topological connectivity information, but only for so long. As the synapses deteriorate below critical value, , the corresponding effective coactivity complex acquires an irreparable amount of topological defects and fails to encode the correct topological map of the environment.
Iv Discussion
Countless observations point out that deteriorations of synapses often accompany memory deficiencies. For example, the recurrent connectivity of CA3 area of the hippocampus and the manytoone projections from the CA3 to the CA1 area Shepherd (); Syntax () suggest that the CA1 cells may provide readouts for the activity of the CA3 place cell assemblies CAs (). Behavioral and cognitive experiments demonstrate that weakening of the synapses between these two areas, a reduction in the number of active neurons in either domain, diminishing neuronal activity and so forth, correlate with learning and memory deficiencies observed, e.g., in Alzheimer’s disease Cacucci (); LaFerla () or in aging subjects Robitsek (); Wilson2 (). However, without a theoretical framework that can link the “synaptic” and the “organismal” scales, the detailed connections between these two scopes of phenomena are hard to trace. For example, if the spike transmission rate in an ensemble of place cells decreases, e.g., by , will the time required to learn the environment increase by , or by ? Does the outcome depend on the “base” level of the transmission probability? Can an increase in learning time caused by synaptic depression always be compensated by increasing the population of active cells, or by elevating their spiking rates? The topological model permits addressing these questions computationally, at a phenomenological level, thus allowing us to move beyond mere correlative descriptions to a deeper understanding of the spatial memory deterioration mechanisms.
The hypothesis about topological nature of the hippocampal map eLife () is broader than the proposed Algebraic Topology (AT) model or the scope of questions that this model allows addressing. For example, the description based on the AT algorithms does not capture biologically relevant metrical differences between topologically equivalent paths or qualitative differences between topologically equivalent environments, e.g., between the widely used W, U or Tmazes, even though such differences are reflected in the place cell spiking patterns and are known to affect animals behavior Frank (); Pastalkova (); Jadhav (). Addressing these differences requires using alternative mathematical apparatuses, e.g., Qualitative Space Representation (QSR) techniques, such as Region Connection Calculi (RCC) Cui (); Hazarika (); Cohn (), which would complement the scope of topological methods used in neuroscience SchemaS (). In the current approach, we use AT instruments to assess a particular scope of questions, namely to estimate the conditions that guarantee structural integrity of the cognitive map and to describe its overall topological shape.
Fundamentally, producing a cognitive map requires two key components: a proper temporal structure of the spike trains and a physiological mechanism for detecting and interpreting neuronal coactivity—a suitable network architecture, a proper distribution of the connectivity strengths, of the parameters of synaptic plasticity, etc. All these components influence spike transmission and detection probabilities, which, in our model, affect the shape and topological structure of the coactivity complex, the statistics of learning times, the structure of the learning region, etc. This produces a quantitative connection between the information processed at the microscopic level (neurons and synapses) and the properties of the largescale representations of space emerging at the organismal level, described here by means of the Persistence Homology theory Zomorodian (); GhristBar (); Edelsbrunner ().
Lastly, it should be pointed out that the lowdimensional components of the coactivity complexes were used above to represent cognitive maps, i.e., frameworks spatial memories. However, the combinations of the coactive place cells, modeled as simplexes of , may represent generic memory elements Wood (); Ginther (). In other words, it can be argued that the net structure of represents not only spatial, but also nonspatial memories—a larger memory framework that can be viewed as a “memory space” Eichenbaum1 (); SchemaM (). Thus, a disintegration of the cell assembly complex caused by deteriorating synapses discussed above may also be viewed as a model of the full memory space decay. From such perspective, it may be noticed that the results of the model parallel the experience of patients acquiring a slowly progressing dementia. For example, the model provides an explanation for the reason cognitive declines often do not manifest until quite a lot of damage has occurred. It also predicts that when the weakening synapses deteriorate beyond the range of parameters within which learning is effective, the damages push the neuronal ensemble beyond the bounds of the learning region. As a result, the failure becomes more frequent, and finally, the brain cannot perform that particular learning task, certain memories or abilities begin to flicker and then are lost mostly for good (Fig. 9).
V Methods
The simulated environment shown on Fig. 1A is designed similarly to the arenas used in typical electrophysiological experiments. Combining such small arenas allows simulating learning in larger, more complex environments Arai (). The simulated trajectory represents nonpreferential, exploratory spatial behavior, with no artificial patterns of moves or favoring of one segment of the environment over another.
Place cell spiking probability was modeled as a Poisson process with the rate
where is the maximal rate of place cell and defines the size of its place field centered at Barbieri (). In an ensemble of place cells, the parameters and , are lognormally distributed with the means and and the variances and . To avoid overly broad or overly narrow distributions, we used additional conditions and , with and PLoS (). In addition, spiking probability was modulated by the wave Arai (); Mizuseki (); Huxter (). The wave also defines the temporal window ms (about two periods) for detecting the place cell spiking coactivity, as suggested by experimental studies Mizuseki (); Huxter (); BuzsakiTh () and by our model Arai (). This value also defines the timestep used in the computations. The place field centers for each computed place field map were randomly and uniformly scattered over the environment.
Place cell ensembles are specified by a triple of parameters and hence the learning region represents a domain of this parameter space. The ensembles studied above contain between and place cells. The ensemble mean peak firing rate ranges from to Hz, and the average place field size ranges between cm and cm ( cm).
Persistent Homology Theory is used to describe the evolving topological shape of the coactivity complexes in terms of their homological invariants Zomorodian (); GhristBar (); Edelsbrunner (). In particular, it allows computing the timedependence of the Betti numbers and deducing the dynamics of its topological loops—their mean lifetimes, their mean lengths, their numbers, etc. Computations were performed using javaplex computational software developed at Stanford University JPlex ().
Acknowledgements
The work was supported by the NSF 1422438 grant.
Vi References
References
 (1) O’Keefe, J., & Nadel, L. The hippocampus as a cognitive map. New York: Clarendon Press; Oxford University Press. (1978).
 (2) Moser, E.I., Kropff, E. & Moser, MB. Place Cells, Grid Cells, and the Brain’s Spatial Representation System. Annu. Rev. Neurosci. 31: 6989 (2008).
 (3) Schmidt, B. & Redish, A.D. Neuroscience: Navigation with a cognitive map. Nature 497: 4243 (2013).
 (4) Gothard, K., Skaggs, W. & McNaughton, B. Dynamics of mismatch correction in the hippocampal ensemble code for space: interaction between path integration and environmental cues. J. Neurosci. 16: 80278040 (1996).
 (5) Leutgeb, J., Leutgeb, S., Treves, A., Meyer, R., Barnes, C., McNaughton, B., Moser, M.B. & Moser E. Progressive transformation of hippocampal neuronal representations in “morphed” environments. Neuron 48: 345358. (2005)
 (6) Wills, T., Lever, C., Cacucci, F., Burgess, N. & O’Keefe J. Attractor dynamics in the hippocampal representation of the local environment. Science 308: 873876. (2005).
 (7) Touretzky, D., Weisman, W., Fuhs, M., Skaggs, W., Fenton, A. & Muller, R. Deforming the hippocampal map. Hippocampus 15: 4155. (2005).
 (8) Dabaghian, Y., Brandt, V., & Frank, L. Reconceiving the hippocampal map as a topological template. eLife. 10.7554/eLife.03476 (2014).
 (9) Poucet, B. & Herrmann, T. Exploratory patterns of rats on a complex maze provide evidence for topological coding. Behav Processes, 53: 155162 (2001).
 (10) Alvernhe, A., Sargolini, F., & Poucet, B. Rats build and update topological representations through exploration. Anim. Cogn. 15: 359368 (2012).
 (11) Wu, X. & Foster, D. Hippocampal replay captures the unique topological structure of a novel environment. J Neurosci., 34: 64596469 (2014).

(12)
De Silva, V. & Ghrist, R.Coverage in sensor networks via persistent homology. Algebraic & Geometric Topology 7: 339â358 (2007).
 (13) Curto, C. & Itskov, V. Cell groups reveal structure of stimulus space, PLoS Comput. Biol., 4: e1000205 (2008).
 (14) Dabaghian, Y., Mémoli, F., Frank, L. & Carlsson, G. A Topological Paradigm for Hippocampal Spatial Map Formation Using Persistent Homology. PLoS Comput. Biol. 8: e1002581 (2012).
 (15) Arai, M., Brandt, V. & Dabaghian, Y. The Effects of Theta Precession on Spatial Learning and Simplicial Complex Dynamics in a Topological Model of the Hippocampal Spatial Map. PLoS Comput. Biol. 10: e1003651 (2014).
 (16) Basso, E., Arai, M. & Dabaghian, Y. Gamma Synchronization Influences Map Formation Time in a Topological Model of Spatial Learning. PLoS Comput. Biol. 12: e1005114 (2016).
 (17) Hoffman, K., Babichev, A. & Dabaghian, Y. A model of topological mapping of space in bat hippocampus, Hippocampus: 26(10): 13451353 (2016).
 (18) Babichev, A., Ji, D., Mémoli, F. & Dabaghian, Y. A Topological Model of the Hippocampal Cell Assembly Network. Front. Comput. Neurosci. 10 (2016).
 (19) Alexandroff, P. Untersuchungen über Gestalt und Lage abgeschlossener Mengen beliebiger Dimension. Annals of Mathematics, 30: 101187 (1928).
 (20) Čech, E. Théorie générale de l’homologie dans un espace quelconque. Fund. Mathematicae, 19: 149183 (1932).
 (21) Hatcher A, Algebraic topology, Cambridge; New York: Cambridge University Press (2002).
 (22) Alexandrov, P.S. Elementary concepts of topology. New York: F. Ungar Pub. Co. (1965).
 (23) Guger, C., Gener, T., Pennartz, C., BrotonsMas, J., Edlinger, G., Bermudez, I., Badia, S., Verschure, P., Schaffelhofer, S. & SanchezVives, M. Realtime Position Reconstruction with Hippocampal Place Cells. Front. Neurosci. 5 (2011).
 (24) Brown, E., Frank, L., Tang, D., Quirk M. & Wilson M. A statistical paradigm for neural spike train decoding applied to position prediction from ensemble firing patterns of rat hippocampal place cells. J. Neurosci., 18: 74117425 (1998).
 (25) Harris, K., Csicsvari, J., Hirase, H., Dragoi, G. & Buzsaki, G. Organization of cell assemblies in the hippocampus. Nature 424: 552556 (2003).
 (26) Harris, K. Neural signatures of cell assembly organization, Nat. Rev. Neurosci., 6: 399407 (2005).
 (27) Jackson, J. & Redish, A. Network dynamics of hippocampal cellassemblies resemble multiple spatial maps within single tasks. Hippocampus 17: 12091229 (2007).
 (28) O’Neill, J., Senior, T., Allen, K., Huxter, J. & Csicsvari, J., Reactivation of experiencedependent cell assembly patterns in the hippocampus. Nat. Neurosci. 11: 209215 (2008).
 (29) Buzsaki, G. Neural syntax: cell assemblies, synapsembles, and readers. Neuron 68: 362385 (2010).
 (30) Babichev, A., Cheng, S. & Dabaghian, Y. Topological schemas of cognitive maps and spatial learning. Front. Comput. Neurosci. 10 (2016).
 (31) Burgess, N. & O’Keefe, J., Cognitive graphs, resistive grids, and the hippocampal representation of space, J. Gen. Physiol., 107: 659662 (1996).
 (32) Muller, R., Stead, M. & Pach, J. The hippocampus as a cognitive graph, J. Gen. Physiol., 107: 663694 (1996).
 (33) Jonsson, J. Simplicial complexes of graphs. Berlin; New York: Springer (2008).
 (34) Soltani, A. & Wang, X. Synaptic computation underlying probabilistic inference. Nat. Neurosci. 13: 112 (2009).
 (35) London, M., Schreibman, A., Hausser, M., Larkum, M. & Segev, I. The information efficacy of a synapse. Nat. Neurosci. 5: 332340 (2002).
 (36) Branco, T., Staras, K., Darcy, K. & Goda, Y. Local Dendritic Activity Sets Release Probability at Hippocampal Synapses. Neuron 59: 475485 (2008).
 (37) Arleo, A., Nieus, T., Bezzi, M., D’Errico, A., D’Angelo, E. & Coenen, O. How Synaptic Release Probability Shapes Neuronal Transmission: InformationTheoretic Analysis in a Cerebellar Granule Cell. Neural Comput. 22: 20312058 (2010).
 (38) Garrido, J., Ros, E. & D’Angelo, E. Spike Timing Regulation on the Millisecond Scale by Distributed Synaptic Plasticity at the Cerebellum Input Stage: A Simulation Study. Front. Comput. Neurosci. 7 (2013).
 (39) Buzsaki, G. & Mizuseki, K. The logdynamic brain: how skewed distributions affect network operations. Nat. Rev. Neurosci. 15(4): 264â278 (2014).
 (40) Barbour, B., Brunel, N., Hakim, V. & Nadal J.P. What can we learn from synaptic weight distributions? Trends Neurosci. 30: 622629 (2007).
 (41) Brunel, N., Hakim, V., Isope, P., Nadal, J.P. & Barbour, B. Optimal Information Storage and the Distribution of Synaptic Weights: Perceptron versus Purkinje Cell. Neuron 43: 745757 (2004).
 (42) Selkoe, D.J. Alzheimer’s Disease Is a Synaptic Failure. Science 298: 789791 (2002).
 (43) Neves, G., Cooke, S., & Bliss, T. Synaptic plasticity, memory and the hippocampus: a neural network approach to causality. Nat. Rev. Neurosci. 9: 6575 (2008).
 (44) Mayford, M., Siegelbaum, S. & Kandel, E. Synapses and Memory Storage. Cold Spring Harbor Perspectives in Biology 4 (2012).
 (45) Ambjørn J, Carfora M, Marzuoli A (1997) The geometry of dynamical triangulations, Berlin, Springer.
 (46) Hamber HW (2009) Quantum gravitation: the Feynman path integral approach, Berlin: Springer. Nuclear Physics B  Proceedings Supplements 94: 689692.
 (47) R. Franzosi, M. Pettini, and L. Spinelli. Topology and phase transitions i. preliminary results. Nuclear Physics B, 782(3): 189 â 218 (2007).
 (48) R. Franzosi and M. Pettini. Topology and phase transitions ii. theorem on a necessary relation. Nuclear Physics B, 782(3): 219 â 240 (2007).
 (49) Donato, I., Gori M., Pettini M., Petri G., De Nigris S., Franzosi R. & Vaccarino F. Persistent homology analysis of phase transitions. Phys. Rev. E 93: 052138 (2016).
 (50) Palop, J., Chin, J., Roberson, E., Wang, J., Thwin, M., BienLy, N., Yoo, J., Ho, K., Yu, G., Kreitzer, A., Finkbeiner, S., Noebels, J. & Mucke, L. Aberrant Excitatory Neuronal Activity and Compensatory Remodeling of Inhibitory Hippocampal Circuits in Mouse Models of Alzheimer’s Disease. Neuron 55: 697711 (2007).
 (51) Hartley, T. & Burgess, N. Complementary memory systems: competition, cooperation and compensation. Trends Neurosci. 28: 169170 (2005).
 (52) Busche, M. & Konnerth, A. Neuronal hyperactivity â A key defect in Alzheimer’s disease? Bioessays 37: 624632 (2015).
 (53) Csicsvari, J., Hirase, H., Czurko, A. & Buzsáki, G. Reliability and State Dependence of Pyramidal CellâInterneuron Synapses in the Hippocampus: an Ensemble Approach in the Behaving Rat. Neuron 21: 179189 (1998).
 (54) Mizuseki, K. & Buzsáki, G. Preconfigured, Skewed Distribution of Firing Rates in the Hippocampus and Entorhinal Cortex. Cell Rep. 4: 10101021 (2013).
 (55) Wilson, M. & McNaughton, B. Dynamics of the hippocampal ensemble code for space. Science 261: 10551058 (1993).
 (56) Leutgeb, S., Leutgeb, J., Treves, A., Moser, M.B. & Moser, E. Distinct Ensemble Codes in Hippocampal Areas CA3 and CA1. Science 305: 12951298 (2004).
 (57) Nithianantharajah, J. & Hannan, A. Enriched environments, experiencedependent plasticity and disorders of the nervous system. Nat. Rev. Neurosci. 7: 697709 (2006).
 (58) Hernan, A., Mahoney, J., Curry, W., Richard, G., Lucas, M., Massey, A., Holmes, G. & Scott, R. Environmental enrichment normalizes hippocampal timing coding in a malformed hippocampus. PLoS One 13: e0191488 (2018).
 (59) Cacucci, F., Yi, M., Wills, T., Chapman, P. & O’Keefe, J. Place cell firing correlates with memory deficits and amyloid plaque burden in Tg2576 Alzheimer mouse model. Proc. Natl. Acad. Sci. 105: 78637868 (2008).
 (60) LaFerla, F. & Oddo, S. Alzheimer’s disease: A, tau and synaptic dysfunction. Trends in Molecular Medicine 11: 170176 (2005).
 (61) Robitsek, R., Fortin, N., Koh, M., Gallagher, M. & Eichenbaum, H. Cognitive aging: a common decline of episodic recollection and spatial memory in rats. J. Neurosci. 28: 89458954 (2008).
 (62) Wilson, I., Ikonen, S., Gureviciene, I., McMahan, R., Gallagher, M., Eichenbaum, H. & Tanila, H. Cognitive aging and the hippocampus: how old rats represent new environments. J. Neurosci. 24: 38703878 (2004).
 (63) White, A. & Best, P. Effects of ethanol on hippocampal placecell and interneuron activity. Brain Res. 876: 154165 (2000).
 (64) Matthews, D., Simson, P. & Best, P. Ethanol alters spatial processing of hippocampal place cells: a mechanism for impaired navigation when intoxicated. Alcohol Clin. Exp. Res. 20: 404407 (1996).
 (65) Robbe, D., Montgomery, S., Thome, A., RuedaOrozco, P., McNaughton & G. Buzsaki. Cannabinoids reveal importance of spike timing coordination in hippocampal function. Nat. Neurosci. 9: 15261533 (2006).
 (66) Robbe, D. & Buzsaki, G. Alteration of theta timescale dynamics of hippocampal place cells by a cannabinoid is associated with memory impairment. J. Neurosci. 29: 1259712605 (2009).
 (67) Kalechstein, A., De la Garza, R., Newton, T., Green, M., Cook, I. & Leuchter, AF. Quantitative EEG abnormalities are associated with memory impairment in recently abstinent methamphetaminedependent individuals. J. Neuropsychiatry Clin. Neurosci. 21: 254258 (2009).
 (68) Silvers, J., Tokunaga, S., Berry, R., White, A. & Matthews, D. Impairments in spatial learning and memory: ethanol, allopregnanolone, and the hippocampus. Brain Res. Rev. 43: 275284 (2003).
 (69) Laxton, A., TangWai, D., McAndrews, M., Zumsteg, D., Wennberg, R., Keren, R., Wherrett, J., Naglie, G., Hamani, C., Smith, G. & Lozano, A. A phase I trial of deep brain stimulation of memory circuits in Alzheimer’s disease. Ann. Neurol. 68: 521534 (2010).
 (70) Shirvalkar, P., Rapp, P. & Shapiro, M. Bidirectional changes to hippocampal thetaâgamma comodulation predict memory for recent spatial episodes. Proc. Nat. Acad. Sci. 107: 70547059 (2010).
 (71) Forman, D. Bochner’s Method for Cell Complexes and Combinatorial Ricci Curvature. Discrete and Computational Geometry 29: 323374 (2003).
 (72) Lewiner, T., Lopes, H. & Tavares, G. Visualizing Forman’s Discrete Vector Field. Visualization and Mathematics III. Berlin, Heidelberg. Springer Berlin Heidelberg. pp. 95112 (2003).
 (73) Weber, M., Jost, J. & Saucan, E. FormanRicci Flow for Change Detection in Large Dynamic Data Sets. Axioms 5: 26 (2016).
 (74) Weber, M., Saucan, E. & Jost, J. Characterizing complex networks with FormanRicci curvature and associated geometric flows. J. Complex Networks 5: 527550 (2017).
 (75) Weber, M., Stelzer, J., Saucan, E., Naitsat, A., Lohmann, G. & Jost, J. Curvaturebased Methods for Brain Network Analysis. arXiv:1707.00180 (2017)
 (76) Sreejith, R., Jost, J., Saucan, E. & Samal, A. Systematic evaluation of a new combinatorial curvature for complex networks. Chaos, Solitons and Fractals 101: 5067 (2017).
 (77) Shepherd, G., The synaptic organization of the brain, Oxford; New York: Oxford University Press (2004).
 (78) Frank, L., Brown, E. and Wilson, M., A comparison of the firing properties of putative excitatory and inhibitory neurons from CA1 and the entorhinal cortex, J. Neurophys., 86(4): 2029–2040 (2001)
 (79) Pastalkova, E., Itskov, V. Amarasingham, A. and Buzsáki, G. Internally generated cell assembly sequences in the rat hippocampus. Science, 321: 1322–1327 (2008).
 (80) Jadhav, S., Kemere, C. German, P. and Frank, L. Awake hippocampal sharpwave ripples support spatial memory, Science, 336: 1454–1458 (2012).
 (81) Cui, Z., Cohn, A., Randell, D. Qualitative and Topological Relationships in Spatial Databases. Proceedings of the Third International Symposium on Advances in Spatial Databases: SpringerVerlag. pp. 296315.(1993)
 (82) Hazarika, S., Cohn, A. Qualitative SpatioTemporal Continuity. Proceedings of the Int. Conference on Spatial Information Theory: Foundations of Geographic Information Science: SpringerVerlag. pp. 92107.(2001)
 (83) Dabaghian Y, Cohn AG, Frank L (2007) Topological maps from signals. Proceedings of the 15th annual ACM international symposium on Advances in geographic information systems. Seattle, Washington: ACM. pp. 14.
 (84) Wood, R., Dudchenko, P., Robitsek, R. & Eichenbaum, H. Hippocampal neurons encode information about different types of memory episodes occurring in the same location. Neuron 27: 623633 (2000).
 (85) Ginther, M., Walsh, D. & Ramus, S. Hippocampal Neurons Encode Different Episodes in an Overlapping Sequence of Odors Task. J. Neurosci. 31: 27062711 (2011).
 (86) Eichenbaum, H., Dudchenko, P., Wood, E., Shapiro, M. & Tanila, H. The hippocampus, memory, and place cells: is it spatial memory or a memory space? Neuron 23: 209226 (1999).
 (87) Babichev, A. & Y. Dabaghian, Y. Topological Schemas of Memory Spaces. Front. Comput. Neurosci. 12 (2018).
 (88) Barbieri, R., Frank, L., Nguyen, D., Quirk, M., Solo, V., Wilson, M. & Brown, E. Dynamic analyses of information encoding in neural ensembles, Neural Comput., 16, pp. 277307 (2004).
 (89) Mizuseki, K., Sirota, A., Pastalkova, E. & Buzsaki, G. Theta oscillations provide temporal windows for local circuit computation in the entorhinalhippocampal loop, Neuron 64: 267280 (2009).
 (90) Huxter, J., Senior, T., Allen, K. & Csicsvari, J. Theta phasespecific codes for twodimensional position, trajectory and heading in the hippocampus. Nat. Neurosci. 11: 587594 (2008).
 (91) Buzsaki, G. Theta rhythm of navigation: link between path integration and landmark navigation, episodic and semantic memory. Hippocampus 15: 827840 (2005).
 (92) Zomorodian, A. & Carlsson, G. Computing persistent homology. Discrete and Computational Geometry 33: 249–274 (2005).
 (93) Ghrist, R. Barcodes: The persistent topology of data, Bull. Amer. Math. Soc., 45: 6175 (2008).
 (94) Edelsbrunner H, Letscher, D., and Zomorodian, A. Topological Persistence and Simplification, Discrete & Computational Geometry 28: 511–533 (2002).
 (95) javaplex freeware, Stanford University, Palo Alto, USA.