Cellular, columnar and modular organization of spatial representations in medial entorhinal cortex

Synaptic interactions between stellate cells and parvalbumin interneurons in layer 2 of the medial entorhinal cortex are organized at the scale of grid cell clusters

Interactions between excitatory and inhibitory neurons are critical to computations in cortical circuits but their organization is difficult to assess with standard electrophysiological approaches. Within the medial entorhinal cortex, representation of location by grid and other spatial cells involves circuits in layer 2 in which excitatory stellate cells interact with each other via inhibitory parvalbumin expressing interneurons. Whether this connectivity is structured to support local circuit computations is unclear. Here, we introduce strategies to address the functional organization of excitatory-inhibitory interactions using crossed Cre- and Flp-driver mouse lines to direct targeted presynaptic optogenetic activation and postsynaptic cell identification. We then use simultaneous patch-clamp recordings from postsynaptic neurons to assess their shared input from optically activated presynaptic populations. We find that extensive axonal projections support spatially organized connectivity between stellate cells and parvalbumin interneurons, such that direct connections are often, but not always, shared by nearby neurons, whereas multisynaptic interactions coordinate inputs to neurons with greater spatial separation. We suggest that direct excitatory-inhibitory synaptic interactions may operate at the scale of grid cell clusters, with local modules defined by excitatory-inhibitory connectivity, while indirect interactions may coordinate activity at the scale of grid cell modules.

Simulation of oscillatory dynamics induced by an approximation of grid cell output

Grid cells, in entorhinal cortex (EC) and related structures, signal animal location relative to hexagonal tilings of 2D space. A number of modeling papers have addressed the question of how grid firing behaviors emerge using (for example) ideas borrowed from dynamical systems (attractors) or from coupled oscillator theory. Here we use a different approach: instead of asking how grid behavior emerges, we take as a given the experimentally observed intracellular potentials of superficial medial EC neurons during grid firing. Employing a detailed neural circuit model modified from a lateral EC model, we then ask how the circuit responds when group of medial EC principal neurons exhibit such potentials, simultaneously with a simulated theta frequency input from the septal nuclei. The model predicts the emergence of robust theta-modulated gamma/beta oscillations, suggestive of oscillations observed in an in vitro medial EC experimental model (Cunningham, M.O., Pervouchine, D.D., Racca, C., Kopell, N.J., Davies, C.H., Jones, R.S.G., Traub, R.D., and Whittington, M.A. (2006). Neuronal metabolism governs cortical network response state. Proc. Natl. Acad. Sci. U S A 103: 5597-5601). Such oscillations result because feedback interneurons tightly synchronize with each other - despite the varying phases of the grid cells - and generate a robust inhibition-based rhythm. The lack of spatial specificity of the model interneurons is consistent with the lack of spatial periodicity in parvalbumin interneurons observed by Buetfering, C., Allen, K., and Monyer, H. (2014). Parvalbumin interneurons provide grid cell-driven recurrent inhibition in the medial entorhinal cortex. Nat. Neurosci. 17: 710-718. If in vivo EC gamma rhythms arise during exploration as our model predicts, there could be implications for interpreting disrupted spatial behavior and gamma oscillations in animal models of Alzheimer's disease and schizophrenia. Noting that experimental intracellular grid cell potentials closely resemble cortical Up states and Down states, during which fast oscillations also occur during Up states, we propose that the co-occurrence of slow principal cell depolarizations and fast network oscillations is a general property of the telencephalon, in both waking and sleep states.

Kinetics and Connectivity Properties of Parvalbumin- and Somatostatin-Positive Inhibition in Layer 2/3 Medial Entorhinal Cortex

Parvalbumin-positive (Pvalb⁺) and somatostatin-positive (Sst⁺) cells are the two largest subgroups of inhibitory interneurons. Studies in visual cortex indicate that synaptic connections between Pvalb⁺ cells are common while connections between Sst⁺ interneurons have not been observed. The inhibitory connectivity and kinetics of these two interneuron subpopulations, however, have not been characterized in medial entorhinal cortex (mEC). Using fluorescence-guided paired recordings in mouse brain slices from interneurons and excitatory cells in layer 2/3 mEC, we found that, unlike neocortical measures, Sst⁺ cells inhibit each other, albeit with a lower probability than Pvalb⁺ cells (18% vs 36% for unidirectional connections). Gap junction connections were also more frequent between Pvalb⁺ cells than between Sst⁺ cells. Pvalb⁺ cells inhibited each other with larger conductances, smaller decay time constants, and shorter delays. Similarly, synaptic connections between Pvalb⁺ and excitatory cells were more likely and expressed faster decay times and shorter delays than those between Sst⁺ and excitatory cells. Inhibitory cells exhibited smaller synaptic decay time constants between interneurons than on their excitatory targets. Inhibition between interneurons also depressed faster, and to a greater extent. Finally, inhibition onto layer 2 pyramidal and stellate cells originating from Pvalb⁺ interneurons were very similar, with no significant differences in connection likelihood, inhibitory amplitude, and decay time. A model of short-term depression fitted to the data indicates that recovery time constants for refilling the available pool are in the range of 50-150 ms and that the fraction of the available pool released on each spike is in the range 0.2-0.5.

Entorhinal Layer II Calbindin-Expressing Neurons Originate Widespread Telencephalic and Intrinsic Projections

In the present study we provide the first systematic and quantitative hodological study of the calbindin-expressing (CB+) principal neurons in layer II of the entorhinal cortex and compared the respective projections of the lateral and medial subdivisions of the entorhinal cortex. Using elaborate quantitative retrograde tracing, complemented by anterograde tracing, we report that the layer II CB+ population comprises neurons with diverse, mainly excitatory projections. At least half of them originate local intrinsic and commissural projections which distribute mainly to layer I and II. We further show that long-range CB+ projections from the two entorhinal subdivisions differ substantially in that MEC projections mainly target field CA1 of the hippocampus, whereas LEC CB+ projections distribute much more widely to a substantial number of known forebrain targets. This connectional difference between the CB+ populations in LEC and MEC is reminiscent of the overall projection pattern of the two entorhinal subdivisions.

Structural modularity and grid activity in the medial entorhinal cortex

Following the groundbreaking discovery of grid cells, the medial entorhinal cortex (MEC) has become the focus of intense anatomical, physiological, and computational investigations. Whether and how grid activity maps onto cell types and cortical architecture is still an open question. Fundamental similarities in microcircuits, function, and connectivity suggest a homology between rodent MEC and human posteromedial entorhinal cortex. Both are specialized for spatial processing and display similar cellular organization, consisting of layer 2 pyramidal/calbindin cell patches superimposed on scattered stellate neurons. Recent data indicate the existence of a further nonoverlapping modular system (zinc patches) within the superficial MEC layers. Zinc and calbindin patches have been shown to receive largely segregated inputs from the presubiculum and parasubiculum. Grid cells are also clustered in the MEC, and we discuss possible structure-function schemes on how grid activity could map onto cortical patch systems. We hypothesize that in the superficial layers of the MEC, anatomical location can be predictive of function; thus relating functional properties and neuronal morphologies to the cortical modules will be necessary for resolving how grid activity maps onto cortical architecture. Imaging or cell identification approaches in freely moving animals will be required for testing this hypothesis.

Complementary Modular Microcircuits of the Rat Medial Entorhinal Cortex

The parahippocampal region is organized into different areas, with the medial entorhinal cortex (MEC), presubiculum and parasubiculum prominent in spatial memory. Here, we also describe a region at the extremity of the MEC and bordering the subicular complex, the medial-most part of the entorhinal cortex. While the subdivisions of hippocampus proper form more or less continuous cell sheets, the superficial layers of the parahippocampal region have a distinct modular architecture. We investigate the spatial distribution, laminar position, and putative connectivity of zinc-positive modules in layer 2 of the MEC of rats and relate them to the calbindin-positive patches previously described in the entorhinal cortex. We found that the zinc-positive modules are complementary to the previously described calbindin-positive patches. We also found that inputs from the presubiculum are directed toward the zinc-positive modules while the calbindin-positive patches received inputs from the parasubiculum. Notably, the dendrites of neurons from layers 3 and 5, positive for Purkinje Cell Protein 4 expression, overlap with the zinc modules. Our data thus indicate that these two complementary modular systems, the calbindin patches and zinc modules, are part of parallel information streams in the hippocampal formation.

Dynamics of Population Activity in Rat Sensory Cortex: Network Correlations Predict Anatomical Arrangement and Information Content

To study the spatiotemporal dynamics of neural activity in a cortical population, we implanted a 10 × 10 microelectrode array in the vibrissal cortex of urethane-anesthetized rats. We recorded spontaneous neuronal activity as well as activity evoked in response to sustained and brief sensory stimulation. To quantify the temporal dynamics of activity, we computed the probability distribution function (PDF) of spiking on one electrode given the observation of a spike on another. The spike-triggered PDFs quantified the strength, temporal delay, and temporal precision of correlated activity across electrodes. Nearby cells showed higher levels of correlation at short delays, whereas distant cells showed lower levels of correlation, which tended to occur at longer delays. We found that functional space built based on the strength of pairwise correlations predicted the anatomical arrangement of electrodes. Moreover, the correlation profile of electrode pairs during spontaneous activity predicted the "signal" and "noise" correlations during sensory stimulation. Finally, mutual information analyses revealed that neurons with stronger correlations to the network during spontaneous activity, conveyed higher information about the sensory stimuli in their evoked response. Given the 400-μm-distance between adjacent electrodes, our functional quantifications unravel the spatiotemporal dynamics of activity among nearby and distant cortical columns.

Conserved size and periodicity of pyramidal patches in layer 2 of medial/caudal entorhinal cortex

To understand the structural basis of grid cell activity, we compare medial entorhinal cortex architecture in layer 2 across five mammalian species (Etruscan shrews, mice, rats, Egyptian fruit bats, and humans), bridging ∼100 million years of evolutionary diversity. Principal neurons in layer 2 are divided into two distinct cell types, pyramidal and stellate, based on morphology, immunoreactivity, and functional properties. We confirm the existence of patches of calbindin-positive pyramidal cells across these species, arranged periodically according to analyses techniques like spatial autocorrelation, grid scores, and modifiable areal unit analysis. In rodents, which show sustained theta oscillations in entorhinal cortex, cholinergic innervation targeted calbindin patches. In bats and humans, which only show intermittent entorhinal theta activity, cholinergic innervation avoided calbindin patches. The organization of calbindin-negative and calbindin-positive cells showed marked differences in entorhinal subregions of the human brain. Layer 2 of the rodent medial and the human caudal entorhinal cortex were structurally similar in that in both species patches of calbindin-positive pyramidal cells were superimposed on scattered stellate cells. The number of calbindin-positive neurons in a patch increased from ∼80 in Etruscan shrews to ∼800 in humans, only an ∼10-fold over a 20,000-fold difference in brain size. The relatively constant size of calbindin patches differs from cortical modules such as barrels, which scale with brain size. Thus, selective pressure appears to conserve the distribution of stellate and pyramidal cells, periodic arrangement of calbindin patches, and relatively constant neuron number in calbindin patches in medial/caudal entorhinal cortex.

Running rewires the neuronal network of adult-born dentate granule cells

Exercise improves cognition in humans and animals. Running increases neurogenesis in the dentate gyrus of the hippocampus, a brain area important for learning and memory. It is unclear how running modifies the circuitry of new dentate gyrus neurons to support their role in memory function. Here we combine retroviral labeling with rabies virus mediated trans-synaptic retrograde tracing to define and quantify new neuron afferent inputs in young adult male C57Bl/6 mice, housed with or without a running wheel for one month. Exercise resulted in a shift in new neuron networks that may promote sparse encoding and pattern separation. Neurogenesis increased in the dorsal, but not the ventral, dentate gyrus by three-fold, whereas afferent traced cell labeling doubled in number. Regional analysis indicated that running differentially affected specific inputs. Within the hippocampus the ratio of innervation from inhibitory interneurons and glutamatergic mossy cells to new neurons was reduced. Distal traced cells were located in sub-cortical and cortical regions, including perirhinal, entorhinal and sensory cortices. Innervation from entorhinal cortex (EC) was augmented, in proportion to the running-induced enhancement of adult neurogenesis. Within EC afferent input and short-term synaptic plasticity from lateral entorhinal cortex, considered to convey contextual information to the hippocampus was increased. Furthermore, running upregulated innervation from regions important for spatial memory and theta rhythm generation, including caudo-medial entorhinal cortex and subcortical medial septum, supra- and medial mammillary nuclei. Altogether, running may facilitate contextual, spatial and temporal information encoding by increasing adult hippocampal neurogenesis and by reorganization of new neuron circuitry.

Cryptic laminar and columnar organization in the dorsolateral pallium of a weakly electric fish

In the weakly electric gymnotiform fish, Apteronotus leptorhynchus, the dorsolateral pallium (DL) receives diencephalic inputs representing electrosensory input utilized for communication and navigation. Cell counts reveal that, similar to thalamocortical projections, many more cells are present in DL than in the diencephalic nucleus that provides it with sensory input. DL is implicated in learning and memory and considered homologous to medial and/or dorsal pallium. The gymnotiform DL has an apparently simple architecture with a random distribution of simple multipolar neurons. We used multiple neurotracer injections in order to study the microcircuitry of DL. Surprisingly, we demonstrated that the intrinsic connectivity of DL is highly organized. It consists of orthogonal laminar and vertical excitatory synaptic connections. The laminar synaptic connections are symmetric sparse, random, and drop off exponentially with distance; they parcellate DL into narrow (60 μm) overlapping cryptic layers. At distances greater than 100 μm, the laminar connections generate a strongly connected directed graph architecture within DL. The vertical connectivity suggests that DL is also organized into cryptic columns; these connections are highly asymmetric, with superficial DL cells preferentially projecting towards deeper cells. Our experimental analyses suggest that the overlapping cryptic columns have a width of 100 μm, in agreement with the minimal distance for strong connectivity. The architecture of DL and the expansive representation of its input, taken together with the strong expression of N-methyl-D-aspartate (NMDA) receptors by its cells, are consistent with theoretical ideas concerning the cortical computations of pattern separation and memory storage via bump attractors.

Theta variation and spatiotemporal scaling along the septotemporal axis of the hippocampus

Hippocampal theta has been related to locomotor speed, attention, anxiety, sensorimotor integration and memory among other emergent phenomena. One difficulty in understanding the function of theta is that the hippocampus (HPC) modulates voluntary behavior at the same time that it processes sensory input. Both functions are correlated with characteristic changes in theta indices. The current review highlights a series of studies examining theta local field potential (LFP) signals across the septotemporal or longitudinal axis of the HPC. While the theta signal is coherent throughout the entirety of the HPC, the amplitude, but not the frequency, of theta varies significantly across its three-dimensional expanse. We suggest that the theta signal offers a rich vein of information about how distributed neuronal ensembles support emergent function. Further, we speculate that emergent function across the long axis varies with respect to spatiotemporal scale. Thus, septal HPC processes details of the proximal spatiotemporal environment while more temporal aspects process larger spaces and wider time-scales. The degree to which emergent functions are supported by the synchronization of theta across the septotemporal axis is an open question. Our working model is that theta synchrony serves to bind ensembles representing varying resolutions of spatiotemporal information at interdependent septotemporal areas of the HPC. Such synchrony and cooperative interactions along the septotemporal axis likely support memory formation and subsequent consolidation and retrieval.

Laminar and dorsoventral molecular organization of the medial entorhinal cortex revealed by large-scale anatomical analysis of gene expression

Neural circuits in the medial entorhinal cortex (MEC) encode an animal's position and orientation in space. Within the MEC spatial representations, including grid and directional firing fields, have a laminar and dorsoventral organization that corresponds to a similar topography of neuronal connectivity and cellular properties. Yet, in part due to the challenges of integrating anatomical data at the resolution of cortical layers and borders, we know little about the molecular components underlying this organization. To address this we develop a new computational pipeline for high-throughput analysis and comparison of in situ hybridization (ISH) images at laminar resolution. We apply this pipeline to ISH data for over 16,000 genes in the Allen Brain Atlas and validate our analysis with RNA sequencing of MEC tissue from adult mice. We find that differential gene expression delineates the borders of the MEC with neighboring brain structures and reveals its laminar and dorsoventral organization. We propose a new molecular basis for distinguishing the deep layers of the MEC and show that their similarity to corresponding layers of neocortex is greater than that of superficial layers. Our analysis identifies ion channel-, cell adhesion- and synapse-related genes as candidates for functional differentiation of MEC layers and for encoding of spatial information at different scales along the dorsoventral axis of the MEC. We also reveal laminar organization of genes related to disease pathology and suggest that a high metabolic demand predisposes layer II to neurodegenerative pathology. In principle, our computational pipeline can be applied to high-throughput analysis of many forms of neuroanatomical data. Our results support the hypothesis that differences in gene expression contribute to functional specialization of superficial layers of the MEC and dorsoventral organization of the scale of spatial representations.

Pyramidal and stellate cell specificity of grid and border representations in layer 2 of medial entorhinal cortex

In medial entorhinal cortex, layer 2 principal cells divide into pyramidal neurons (mostly calbindin positive) and dentate gyrus-projecting stellate cells (mostly calbindin negative). We juxtacellularly labeled layer 2 neurons in freely moving animals, but small sample size prevented establishing unequivocal structure-function relationships. We show, however, that spike locking to theta oscillations allows assigning unidentified extracellular recordings to pyramidal and stellate cells with ∼83% and ∼89% specificity, respectively. In pooled anatomically identified and theta-locking-assigned recordings, nonspatial discharges dominated, and weakly hexagonal spatial discharges and head-direction selectivity were observed in both cell types. Clear grid discharges were rare and mostly classified as pyramids (19%, 19/99 putative pyramids versus 3%, 3/94 putative stellates). Most border cells were classified as stellate (11%, 10/94 putative stellates versus 1%, 1/99 putative pyramids). Our data suggest weakly theta-locked stellate border cells provide spatial input to dentate gyrus, whereas strongly theta-locked grid discharges occur mainly in hexagonally arranged pyramidal cell patches and do not feed into dentate gyrus.

•

Stellates and pyramids can be putatively assigned by spike locking to theta

•

Grid cells in entorhinal layer 2 largely correspond to putative pyramidal cells

•

Border responses in layer 2 largely correspond to putative stellate cells

Juxtacellular recording and morphological identification of single neurons in freely moving rats

It is well established that neural circuits consist of a great diversity of cell types, but very little is known about how neuronal diversity contributes to cognition and behavior. One approach to addressing this problem is to directly link cellular diversity to neuronal activity recorded in vivo in behaving animals. Here we describe the technical procedures for obtaining juxtacellular recordings from single neurons in trained rats engaged in exploratory behavior. The recorded neurons can be labeled to allow subsequent anatomical identification. In its current format, the protocol can be used for resolving the cellular identity of spatially modulated neurons (i.e., head-direction cells and grid cells), which form the basis of the animal's internal representation of space, but this approach can easily be extended to other unrestrained behaviors. The procedures described here, from the beginning of animal training to the histological processing of brain sections, can be completed in ≈ 3-4 weeks.

Nonlinear properties of medial entorhinal cortex neurons reveal frequency selectivity during multi-sinusoidal stimulation

The neurons in layer II of the medial entorhinal cortex are part of the grid cell network involved in the representation of space. Many of these neurons are likely to be stellate cells with specific oscillatory and firing properties important for their function. A fundamental understanding of the nonlinear basis of these oscillatory properties is critical for the development of theories of grid cell firing. In order to evaluate the behavior of stellate neurons, measurements of their quadratic responses were used to estimate a second order Volterra kernel. This paper uses an operator theory, termed quadratic sinusoidal analysis (QSA), which quantitatively determines that the quadratic response accounts for a major part of the nonlinearity observed at membrane potential levels characteristic of normal synaptic events. Practically, neurons were probed with multi-sinusoidal stimulations to determine a Hermitian operator that captures the quadratic function in the frequency domain. We have shown that the frequency content of the stimulation plays an important role in the characteristics of the nonlinear response, which can distort the linear response as well. Stimulations with enhanced low frequency amplitudes evoked a different nonlinear response than broadband profiles. The nonlinear analysis was also applied to spike frequencies and it was shown that the nonlinear response of subthreshold membrane potential at resonance frequencies near the threshold is similar to the nonlinear response of spike trains.