Kavli Institute for Systems Neuroscience

Episodic memory is defined as the ability to recall events in a spatiotemporal context. Formation of such memories is critically dependent on the hippocampal formation and its inputs from the entorhinal cortex. To be able to support the formation of episodic memories, entorhinal cortex and hippocampal formation should contain a neuronal code that follows several requirements. First, the code should include information about position of the agent ("where"), sequence of events ("when"), and the content of the experience itself ("what"). Second, the code should arise instantly thereby being able to support memory formation of one-shot experiences. For successful encoding and to avoid interference between memories during recall, variations in location, time, or in content of experience should result in unique ensemble activity. Finally, the code should capture several different resolutions of experience so that the necessary details relevant for future memory-based predictions will be stored. We review how neuronal codes in entorhinal cortex and hippocampus follow these requirements and argue that during formation of episodic memories entorhinal cortex provides hippocampus with instant information about ongoing experience. Such information originates from (a) spatially modulated neurons in medial entorhinal cortex, including grid cells, which provide a stable and universal positional metric of the environment; (b) a continuously varying signal in lateral entorhinal cortex providing a code for the temporal progression of events; and (c) entorhinal neurons coding the content of experiences exemplified by object-coding and odor-selective neurons. During formation of episodic memories, information from these systems are thought to be encoded as unique sequential ensemble activity in hippocampus, thereby encoding associations between the content of an event and its spatial and temporal contexts. Upon exposure to parts of the encoded stimuli, activity in these ensembles can be reinstated, leading to reactivation of the encoded activity pattern and memory recollection.

Detecting pathogens and mounting immune responses upon infection is crucial for animal health. However, these responses come at a high metabolic price (McKean and Lazzaro, 2011, Kominsky et al., 2010), and avoiding pathogens before infection may be advantageous. The bacterial endotoxins lipopolysaccharides (LPS) are important immune system infection cues (Abbas et al., 2014), but it remains unknown whether animals possess sensory mechanisms to detect them prior to infection. Here we show that Drosophila melanogaster display strong aversive responses to LPS and that gustatory neurons expressing Gr66a bitter receptors mediate avoidance of LPS in feeding and egg laying assays. We found the expression of the chemosensory cation channel dTRPA1 in these cells to be necessary and sufficient for LPS avoidance. Furthermore, LPS stimulates Drosophila neurons in a TRPA1-dependent manner and activates exogenous dTRPA1 channels in human cells. Our findings demonstrate that flies detect bacterial endotoxins via a gustatory pathway through TRPA1 activation as conserved molecular mechanism.

The enormous potential of modern molecular neuroanatomical tools lies in their ability to determine the precise connectivity of the neuronal cell types comprising the innate circuitry of the brain. We used transgenically targeted viral tracing to identify the monosynaptic inputs to the projection neurons of layer II of medial entorhinal cortex (MEC-LII) in mice. These neurons are not only major inputs to the hippocampus, the structure most clearly implicated in learning and memory, they also are "grid cells." Here we address the question of what kinds of inputs are specifically targeting these MEC-LII cells. Cell-specific infection of MEC-LII with recombinant rabies virus results in unambiguous labeling of monosynaptic inputs. Furthermore, ratios of labeled neurons in different regions are largely consistent between animals, suggesting that label reflects density of innervation. While the results mostly confirm prior anatomical work, they also reveal a novel major direct input to MEC-LII from hippocampal pyramidal neurons. Interestingly, the vast majority of these direct hippocampal inputs arise not from the major hippocampal subfields of CA1 and CA3, but from area CA2, a region that has historically been thought to merely be a transitional zone between CA3 and CA1. We confirmed this unexpected result using conventional tracing techniques in both rats and mice.

High-resolution whole-brain microscopy provides a means for post hoc determination of the location of implanted devices and labelled cell populations that are necessary to interpret in vivo experiments designed to understand brain function. Here we have developed two plugins (brainreg and brainreg-segment) for the Python-based image viewer napari, to accurately map any object in a common coordinate space. We analysed the position of dye-labelled electrode tracks and two-photon imaged cell populations expressing fluorescent proteins. The precise location of probes and cells were physiologically interrogated and revealed accurate segmentation with near-cellular resolution.

The study of complex systems is limited by the fact that only few variables are accessible for modeling and sampling, which are not necessarily the most relevant ones to explain the systems behavior. In addition, empirical data typically under sample the space of possible states. We study a generic framework where a complex system is seen as a system of many interacting degrees of freedom, which are known only in part, that optimize a given function. We show that the underlying distribution with respect to the known variables has the Boltzmann form, with a temperature that depends on the number of unknown variables. In particular, when the unknown part of the objective function decays faster than exponential, the temperature decreases as the number of variables increases. We show in the representative case of the Gaussian distribution, that models are predictable only when the number of relevant variables is less than a critical threshold. As a further consequence, we show that the information that a sample contains on the behavior of the system is quantified by the entropy of the frequency with which different states occur. This allows us to characterize the properties of maximally informative samples: in the under-sampling regime, the most informative frequency size distributions have power law behavior and Zipf's law emerges at the crossover between the under sampled regime and the regime where the sample contains enough statistics to make inference on the behavior of the system. These ideas are illustrated in some applications, showing that they can be used to identify relevant variables or to select most informative representations of data, e.g. in data clustering.

Viewing behavior provides a window into many central aspects of human cognition and health, and it is an important variable of interest or confound in many functional magnetic resonance imaging (fMRI) studies. To make eye tracking freely and widely available for MRI research, we developed DeepMReye, a convolutional neural network (CNN) that decodes gaze position from the magnetic resonance signal of the eyeballs. It performs cameraless eye tracking at subimaging temporal resolution in held-out participants with little training data and across a broad range of scanning protocols. Critically, it works even in existing datasets and when the eyes are closed. Decoded eye movements explain network-wide brain activity also in regions not associated with oculomotor function. This work emphasizes the importance of eye tracking for the interpretation of fMRI results and provides an open source software solution that is widely applicable in research and clinical settings.

Efficient sampling of visual information requires a coordination of eye movements and ongoing brain oscillations. Using intracranial and magnetoencephalography (MEG) recordings, we show that saccades are locked to the phase of visual alpha oscillations and that this coordination is related to successful mnemonic encoding of visual scenes. Furthermore, parahippocampal and retrosplenial cortex involvement in this coordination reflects effective vision-to-memory mapping, highlighting the importance of neural oscillations for the interaction between visual and memory domains.

We study the statistics of spike trains of simultaneously recorded grid cells in freely behaving rats. We evaluate pairwise correlations between these cells and, using a maximum entropy kinetic pairwise model (kinetic Ising model), study their functional connectivity. Even when we account for the covariations in firing rates due to overlapping fields, both the pairwise correlations and functional connections decay as a function of the shortest distance between the vertices of the spatial firing pattern of pairs of grid cells, i.e. their phase difference. They take positive values between cells with nearby phases and approach zero or negative values for larger phase differences. We find similar results also when, in addition to correlations due to overlapping fields, we account for correlations due to theta oscillations and head directional inputs. The inferred connections between neurons in the same module and those from different modules can be both negative and positive, with a mean close to zero, but with the strongest inferred connections found between cells of the same module. Taken together, our results suggest that grid cells in the same module do indeed form a local network of interconnected neurons with a functional connectivity that supports a role for attractor dynamics in the generation of grid pattern.

A key aspect of conceptual knowledge is that it can be flexibly applied at different levels of abstraction, implying a hierarchical organization. It is yet unclear how this hierarchical structure is acquired and represented in the brain. Here we investigate the computations underlying the acquisition and representation of the hierarchical structure of conceptual knowledge in the hippocampal-prefrontal system of 32 human participants (22 females). We assessed the hierarchical nature of learning during a novel tree-like categorization task via computational model comparisons. The winning model allowed to extract and quantify estimates for accumulation and updating of hierarchical compared with single-feature-based concepts from behavior. We find that mPFC tracks accumulation of hierarchical conceptual knowledge over time, and mPFC and hippocampus both support trial-to-trial updating. As a function of those learning parameters, mPFC and hippocampus further show connectivity changes to rostro-lateral PFC, which ultimately represented the hierarchical structure of the concept in the final stages of learning. Our results suggest that mPFC and hippocampus support the integration of accumulated evidence and instantaneous updates into hierarchical concept representations in rostro-lateral PFC. <b>SIGNIFICANCE STATEMENT</b> A hallmark of human cognition is the flexible use of conceptual knowledge at different levels of abstraction, ranging from a coarse category level to a fine-grained subcategory level. While previous work probed the representational geometry of long-term category knowledge, it is unclear how this hierarchical structure inherent to conceptual knowledge is acquired and represented. By combining a novel hierarchical concept learning task with computational modeling of categorization behavior and concurrent fMRI, we differentiate the roles of key concept learning regions in hippocampus and PFC in learning computations and the representation of a hierarchical category structure.

Although a variety of remarkable molecular tools for studying neural circuits have recently been developed, the ability to deploy them in particular neuronal subtypes is limited by the fact that native promoters are almost never specific enough. We recently showed that one can generate transgenic mice with anatomical specificity surpassing that of native promoters by combining enhancers uniquely active in particular brain regions with a heterologous minimal promoter, an approach we call EDGE (Enhancer-Driven Gene Expression). Here we extend this strategy to the generation of viral (rAAV) vectors, showing that some EDGE rAAVs can recapitulate the specificity of the corresponding transgenic lines in wild-type animals, even of another species. This approach thus holds the promise of enabling circuit-specific manipulations in wild-type animals, not only enhancing our understanding of brain function, but perhaps one day even providing novel therapeutic avenues to approach disorders of the brain.

Tract-tracing studies in primates indicate that different subregions of the medial temporal lobe (MTL) are connected with multiple brain regions. However, no clear framework defining the distributed anatomy associated with the human MTL exists. This gap in knowledge originates in notoriously low MRI data quality in the anterior human MTL and in group-level blurring of idiosyncratic anatomy between adjacent brain regions, such as entorhinal and perirhinal cortices, and parahippocampal areas TH/TF. Using MRI, we intensively scanned four human individuals and collected whole-brain data with unprecedented MTL signal quality. Following detailed exploration of cortical networks associated with MTL subregions within each individual, we discovered three biologically meaningful networks associated with the entorhinal cortex, perirhinal cortex, and parahippocampal area TH, respectively. Our findings define the anatomical constraints within which human mnemonic functions must operate and are insightful for examining the evolutionary trajectory of the MTL connectivity across species.

Neurons subject to a common non-stationary input may exhibit a correlated firing behavior. Correlations in the statistics of neural spike trains also arise as the effect of interaction between neurons. Here we show that these two situations can be distinguished, with machine learning techniques, provided the data are rich enough. In order to do this, we study the problem of inferring a kinetic Ising model, stationary or nonstationary, from the available data. We apply the inference procedure to two data sets: one from salamander retinal ganglion cells and the other from a realistic computational cortical network model. We show that many aspects of the concerted activity of the salamander retinal neurons can be traced simply to the external input. A model of non-interacting neurons subject to a non-stationary external field outperforms a model with stationary input with couplings between neurons, even accounting for the differences in the number of model parameters. When couplings are added to the non-stationary model, for the retinal data, little is gained: the inferred couplings are generally not significant. Likewise, the distribution of the sizes of sets of neurons that spike simultaneously and the frequency of spike patterns as function of their rank (Zipf plots) are well-explained by an independent-neuron model with time-dependent external input, and adding connections to such a model does not offer significant improvement. For the cortical model data, robust couplings, well correlated with the real connections, can be inferred using the non-stationary model. Adding connections to this model slightly improves the agreement with the data for the probability of synchronous spikes but hardly affects the Zipf plot.

The amyloid-beta peptide (Aβ) is thought to have prion-like properties promoting its spread throughout the brain in Alzheimer's disease (AD). However, the cellular mechanism(s) of this spread remains unclear. Here, we show an important role of intracellular Aβ in its prion-like spread. We demonstrate that an intracellular source of Aβ can induce amyloid plaques in vivo via hippocampal injection. We show that hippocampal injection of mouse AD brain homogenate not only induces plaques, but also damages interneurons and affects intracellular Aβ levels in synaptically connected brain areas, paralleling cellular changes seen in AD. Furthermore, in a primary neuron AD model, exposure of picomolar amounts of brain-derived Aβ leads to an apparent redistribution of Aβ from soma to processes and dystrophic neurites. We also observe that such neuritic dystrophies associate with plaque formation in AD-transgenic mice. Finally, using cellular models, we propose a mechanism for how intracellular accumulation of Aβ disturbs homeostatic control of Aβ levels and can contribute to the up to 10,000-fold increase of Aβ in the AD brain. Our data indicate an essential role for intracellular prion-like Aβ and its synaptic spread in the pathogenesis of AD.

How do we decide what we do? This is the essence of action control, the process of selecting the most appropriate response among multiple possible choices. Suboptimal action control can involve a failure to initiate or adapt actions, or conversely it can involve making actions impulsively. There has been an increasing focus on the specific role of the subthalamic nucleus (STN) in action control. This has been fueled by the clinical relevance of this basal ganglia nucleus as a target for deep brain stimulation (DBS), primarily in Parkinson's disease but also in obsessive-compulsive disorder. The context of DBS has opened windows to study STN function in ways that link neuroscientific and clinical fields closely together, contributing to an exceptionally high level of two-way translation. In this review, we first outline the role of the STN in both motor and nonmotor action control, and then discuss how these functions might be implemented by neuronal activity in the STN. Gaining a better understanding of these topics will not only provide important insights into the neurophysiology of action control but also the pathophysiological mechanisms relevant for several brain disorders and their therapies.

The brain derives cognitive maps from sensory experience that guide memory formation and behavior. Despite extensive efforts, it still remains unclear how the underlying population activity unfolds during spatial navigation and how it relates to memory performance. To examine these processes, we combined 7T-fMRI with a kernel-based encoding model of virtual navigation to map world-centered directional tuning across the human cortex. First, we present an in-depth analysis of directional tuning in visual, retrosplenial, parahippocampal and medial temporal cortices. Second, we show that tuning strength, width and topology of this directional code during memory-guided navigation depend on successful encoding of the environment. Finally, we show that participants' locomotory state influences this tuning in sensory and mnemonic regions such as the hippocampus. We demonstrate a direct link between neural population tuning and human cognition, where high-level memory processing interacts with network-wide visuospatial coding in the service of behavior.

In the mammalian hippocampus, the dentate gyrus (DG) is characterized by sparse and powerful unidirectional projections to CA3 pyramidal cells, the so-called mossy fibers. Mossy fiber synapses appear to duplicate, in terms of the information they convey, what CA3 cells already receive from entorhinal cortex layer II cells, which project both to the dentate gyrus and to CA3. Computational models of episodic memory have hypothesized that the function of the mossy fibers is to enforce a new, well-separated pattern of activity onto CA3 cells, to represent a new memory, prevailing over the interference produced by the traces of older memories already stored on CA3 recurrent collateral connections. Can this hypothesis apply also to spatial representations, as described by recent neurophysiological recordings in rats? To address this issue quantitatively, we estimate the amount of information DG can impart on a new CA3 pattern of spatial activity, using both mathematical analysis and computer simulations of a simplified model. We confirm that, also in the spatial case, the observed sparse connectivity and level of activity are most appropriate for driving memory storage-and not to initiate retrieval. Surprisingly, the model also indicates that even when DG codes just for space, much of the information it passes on to CA3 acquires a non-spatial and episodic character, akin to that of a random number generator. It is suggested that further hippocampal processing is required to make full spatial use of DG inputs.

Theta oscillations are thought to play a critical role in neuronal information processing, especially in the hippocampal region, where their presence is particularly salient. A detailed description of theta dynamics in this region has revealed not only a consortium of layer-specific theta dipoles, but also within-layer differences in the expression of theta. This complex and articulated arrangement of current flows is reflected in the way neuronal firing is modulated in time. Several models have proposed that these different theta modulators flexibly coordinate hippocampal regions, to support associative memory formation and retrieval. Here, we summarily review different approaches related to this issue and we describe a mechanism, based on experimental and simulation results, for memory retrieval in CA3 involving theta modulation.

We derive and study dynamical TAP equations for Ising spin glasses obeying both synchronous and asynchronous dynamics using a generating functional approach. The system can have an asymmetric coupling matrix, and the external fields can be time-dependent. In the synchronously updated model, the TAP equations take the form of self consistent equations for magnetizations at time $t+1$, given the magnetizations at time $t$. In the asynchronously updated model, the TAP equations determine the time derivatives of the magnetizations at each time, again via self consistent equations, given the current values of the magnetizations. Numerical simulations suggest that the TAP equations become exact for large systems.

We study inference and reconstruction of couplings in a partially observed kinetic Ising model. With hidden spins, calculating the likelihood of a sequence of observed spin configurations requires performing a trace over the configurations of the hidden ones. This, as we show, can be represented as a path integral. Using this representation, we demonstrate that systematic approximate inference and learning rules can be derived using dynamical mean-field theory. Although naive mean-field theory leads to an unstable learning rule, taking into account Gaussian corrections allows learning the couplings involving hidden nodes. It also improves learning of the couplings between the observed nodes compared to when hidden nodes are ignored.

We describe how the couplings in an asynchronous kinetic Ising model can be inferred. We consider two cases: one in which we know both the spin history and the update times and one in which we know only the spin history. For the first case, we show that one can average over all possible choices of update times to obtain a learning rule that depends only on spin correlations and can also be derived from the equations of motion for the correlations. For the second case, the same rule can be derived within a further decoupling approximation. We study all methods numerically for fully asymmetric Sherrington-Kirkpatrick models, varying the data length, system size, temperature, and external field. Good convergence is observed in accordance with the theoretical expectations.