Estimation of a large-scale generalized Volterra model for neural ensembles with group lasso and local coordinate descent

Estimation of neural models based on observed spike timing faces challenges as the amount of recorded units increases, especially when identifying detailed model features. Given that neural regions are generally sparsely connected, input selection is a critical step in model estimation but oftentimes computationally and theoretically challenging. In this paper, we detail an efficient methodology for estimating a sparse, nonlinear dynamical multiple-input, single-output model (MISO) applicable to large-scale (n > 50) single-unit recorded activity. The main contribution of this paper is the complete implementation of a principled group-lasso and local coordinate descent (LCD) algorithm into a generalized Volterra model (GVM) framework to achieve efficient sparse model estimation. Input selection is achieved with group-lasso by simultaneously selecting groups of parameters that are associated with each input. LCD yields efficient computation as the amount of inputs and parameters increase. We investigate and validate the performance of this estimation procedure with the application to a 64 input simulated model.

A 3-D admittance-level computational model of a rat hippocampus for improving prosthetic design

Hippocampal prosthetic devices have been developed to bridge the gap between functioning portions of the hippocampus, in order to restore lost memory functionality in those suffering from brain injury or diseases. One approach taken in recent neuroprosthetic design is to use a multi-input, multi-output device that reads data from the CA3 in the hippocampus and electrically stimulates the CA1 in an attempt to mimic the appropriate firing pattern that would occur naturally between the two areas. However, further study needs to be conducted in order to optimize electrode placement, pulse magnitude, and shape for creating the appropriate firing pattern. This paper describes the creation and implementation of an anatomically correct 3D model of the hippocampus to simulate the electric field patterns and axonal activation from electrical stimulation due to an implanted electrode array. The activating function was applied to the voltage results to determine the firing patterns in possible axon locations within the CA1.

Unstable periodic orbits in human epileptic hippocampal slices

Inter-ictal activity is studied in hippocampal slices resected from patients with epilepsy using local field potential recording. Inter-ictal activity in the dentate gyrus (DG) is induced by high-potassium (8 mM), low-magnesium (0.25 mM) aCSF with additional 100 μM 4-aminopyridine(4-AP). The dynamics of the inter-ictal activity is investigated by developing the first return map with inter-pulse intervals. Unstable periodic orbits (UPOs) are detected in the hippocampal slice at the DG area according to both the topological recurrence method and the periodic orbit transform method. Surrogate analysis suggests the presence of UPOs in hippocampal slices from patients with epilepsy. This finding also suggests that inter-ictal activity is a chaotic system and will allow us to apply chaos control techniques to manipulate inter-ictal activity.

Interactions between Inhibitory Interneurons and Excitatory Associational Circuitry in Determining Spatio-Temporal Dynamics of Hippocampal Dentate Granule Cells: A Large-Scale Computational Study

This paper reports on findings from a million-cell granule cell model of the rat dentate gyrus that was used to explore the contributions of local interneuronal and associational circuits to network-level activity. The model contains experimentally derived morphological parameters for granule cells, which each contain approximately 200 compartments, and biophysical parameters for granule cells, basket cells, and mossy cells that were based both on electrophysiological data and previously published models. Synaptic input to cells in the model consisted of glutamatergic AMPA-like EPSPs and GABAergic-like IPSPs from excitatory and inhibitory neurons, respectively. The main source of input to the model was from layer II entorhinal cortical neurons. Network connectivity was constrained by the topography of the system, and was derived from axonal transport studies, which provided details about the spatial spread of axonal terminal fields, as well as how subregions of the medial and lateral entorhinal cortices project to subregions of the dentate gyrus. Results of this study show that strong feedback inhibition from the basket cell population can cause high-frequency rhythmicity in granule cells, while the strength of feedforward inhibition serves to scale the total amount of granule cell activity. Results furthermore show that the topography of local interneuronal circuits can have just as strong an impact on the development of spatio-temporal clusters in the granule cell population as the perforant path topography does, both sharpening existing clusters and introducing new ones with a greater spatial extent. Finally, results show that the interactions between the inhibitory and associational loops can cause high frequency oscillations that are modulated by a low-frequency oscillatory signal. These results serve to further illustrate the importance of topographical constraints on a global signal processing feature of a neural network, while also illustrating how rich spatio-temporal and oscillatory dynamics can evolve from a relatively small number of interacting local circuits.

Synaptic Efficacy as a Function of Ionotropic Receptor Distribution: A Computational Study

Glutamatergic synapses are the most prevalent functional elements of information processing in the brain. Changes in pre-synaptic activity and in the function of various post-synaptic elements contribute to generate a large variety of synaptic responses. Previous studies have explored postsynaptic factors responsible for regulating synaptic strength variations, but have given far less importance to synaptic geometry, and more specifically to the subcellular distribution of ionotropic receptors. We analyzed the functional effects resulting from changing the subsynaptic localization of ionotropic receptors by using a hippocampal synaptic computational framework. The present study was performed using the EONS (Elementary Objects of the Nervous System) synaptic modeling platform, which was specifically developed to explore the roles of subsynaptic elements as well as their interactions, and that of synaptic geometry. More specifically, we determined the effects of changing the localization of ionotropic receptors relative to the presynaptic glutamate release site, on synaptic efficacy and its variations following single pulse and paired-pulse stimulation protocols. The results indicate that changes in synaptic geometry do have consequences on synaptic efficacy and its dynamics.

Generalized Volterra kernel model identification of spike-timing-dependent plasticity from simulated spiking activity

This paper presents a methodology to estimate a learning rule that governs activity-dependent plasticity from behaviorally recorded spiking events. To demonstrate this framework, we simulate a probabilistic spiking neuron with spike-timing-dependent plasticity (STDP) and estimate all model parameters from the simulated spiking data. In the neuron model, output spiking activity is generated by the combination of noise, feedback from the output, and an input-feedforward component whose magnitude is modulated by synaptic weight. The synaptic weight is calculated with STDP with the following features: (1) weight change based on the relative timing of input-output spike pairs, (2) prolonged plasticity induction, and (3) considerations for system stability. Estimation of all model parameters is achieved iteratively by formulating the model as a generalized linear model with Volterra kernels and basis function expansion. Successful estimation of all model parameters in this study demonstrates the feasibility of this approach for in-vivo experimental studies. Furthermore, the consideration of system stability and prolonged plasticity induction enhances the ability to capture how STDP affects a neural population's signal transformation properties over a realistic time course. Plasticity characterization with this estimation method could yield insights into functional implications of STDP and be incorporated into a cortical prosthesis.

Distributed encoding of spatial and object categories in primate hippocampal microcircuits

The primate hippocampus plays critical roles in the encoding, representation, categorization and retrieval of cognitive information. Such cognitive abilities may use the transformational input-output properties of hippocampal laminar microcircuitry to generate spatial representations and to categorize features of objects, images, and their numeric characteristics. Four nonhuman primates were trained in a delayed-match-to-sample (DMS) task while multi-neuron activity was simultaneously recorded from the CA1 and CA3 hippocampal cell fields. The results show differential encoding of spatial location and categorization of images presented as relevant stimuli in the task. Individual hippocampal cells encoded visual stimuli only on specific types of trials in which retention of either, the Sample image, or the spatial position of the Sample image indicated at the beginning of the trial, was required. Consistent with such encoding, it was shown that patterned microstimulation applied during Sample image presentation facilitated selection of either Sample image spatial locations or types of images, during the Match phase of the task. These findings support the existence of specific codes for spatial and numeric object representations in primate hippocampus which can be applied on differentially signaled trials. Moreover, the transformational properties of hippocampal microcircuitry, together with the patterned microstimulation are supporting the practical importance of this approach for cognitive enhancement and rehabilitation, needed for memory neuroprosthetics.

Volterra representation enables modeling of complex synaptic nonlinear dynamics in large-scale simulations

Chemical synapses are comprised of a wide collection of intricate signaling pathways involving complex dynamics. These mechanisms are often reduced to simple spikes or exponential representations in order to enable computer simulations at higher spatial levels of complexity. However, these representations cannot capture important nonlinear dynamics found in synaptic transmission. Here, we propose an input-output (IO) synapse model capable of generating complex nonlinear dynamics while maintaining low computational complexity. This IO synapse model is an extension of a detailed mechanistic glutamatergic synapse model capable of capturing the input-output relationships of the mechanistic model using the Volterra functional power series. We demonstrate that the IO synapse model is able to successfully track the nonlinear dynamics of the synapse up to the third order with high accuracy. We also evaluate the accuracy of the IO synapse model at different input frequencies and compared its performance with that of kinetic models in compartmental neuron models. Our results demonstrate that the IO synapse model is capable of efficiently replicating complex nonlinear dynamics that were represented in the original mechanistic model and provide a method to replicate complex and diverse synaptic transmission within neuron network simulations.

Hippocampal closed-loop modeling and implications for seizure stimulation design

OBJECTIVE

Traditional hippocampal modeling has focused on the series of feedforward synapses known as the trisynaptic pathway. However, feedback connections from CA1 back to the hippocampus through the entorhinal cortex (EC) actually make the hippocampus a closed-loop system. By constructing a functional closed-loop model of the hippocampus, one may learn how both physiological and epileptic oscillations emerge and design efficient neurostimulation patterns to abate such oscillations.

APPROACH

Point process input-output models where estimated from recorded rodent hippocampal data to describe the nonlinear dynamical transformation from CA3 → CA1, via the schaffer-collateral synapse, and CA1 → CA3 via the EC. Each Volterra-like subsystem was composed of linear dynamics (principal dynamic modes) followed by static nonlinearities. The two subsystems were then wired together to produce the full closed-loop model of the hippocampus.

MAIN RESULTS

Closed-loop connectivity was found to be necessary for the emergence of theta resonances as seen in recorded data, thus validating the model. The model was then used to identify frequency parameters for the design of neurostimulation patterns to abate seizures.

SIGNIFICANCE

Deep-brain stimulation (DBS) is a new and promising therapy for intractable seizures. Currently, there is no efficient way to determine optimal frequency parameters for DBS, or even whether periodic or broadband stimuli are optimal. Data-based computational models have the potential to be used as a testbed for designing optimal DBS patterns for individual patients. However, in order for these models to be successful they must incorporate the complex closed-loop structure of the seizure focus. This study serves as a proof-of-concept of using such models to design efficient personalized DBS patterns for epilepsy.

An Application Specific Instruction Set Processor (ASIP) for Adaptive Filters in Neural Prosthetics

Neural coding is an essential process for neuroprosthetic design, in which adaptive filters have been widely utilized. In a practical application, it is needed to switch between different filters, which could be based on continuous observations or point process, when the neuron models, conditions, or system requirements have changed. As candidates of coding chip for neural prostheses, low-power general purpose processors are not computationally efficient especially for large scale neural population coding. Application specific integrated circuits (ASICs) do not have flexibility to switch between different adaptive filters while the cost for design and fabrication is formidable. In this research work, we explore an application specific instruction set processor (ASIP) for adaptive filters in neural decoding activity. The proposed architecture focuses on efficient computation for the most time-consuming matrix/vector operations among commonly used adaptive filters, being able to provide both flexibility and throughput. Evaluation and implementation results are provided to demonstrate that the proposed ASIP design is area-efficient while being competitive to commercial CPUs in computational performance.

Laguerre-Volterra identification of spike-timing-dependent plasticity from spiking activity: a simulation study

This paper presents a Laguerre-Volterra methodology for identifying a plasticity learning rule from spiking neural data with four components: 1) By analyzing input-output spiking data, the effective contribution of an input on the output firing probability can be quantified with weighted Volterra kernels. 2) The weight of these Volterra kernels can be tracked over time using the stochastic state point processing filtering algorithm (SSPPF) 3) Plasticity system Volterra kernels can be estimated by treating the tracked change in weight over time as the plasticity system output and the spike timing data as the input. 4) Laguerre expansion of all Volterra kernels allows for minimization of open parameters during estimation steps. A single input spiking neuron with Spike-timing-dependent plasticity (STDP) and prolonged STDP induction is simulated. Using the spiking data from this simulation, the amplitude of the STDP learning rule and the time course of the induction is accurately estimated. This framework can be applied to identify plasticity for more complicated plasticity paradigms and is applicable to in vivo data.

A Million-Plus Neuron Model of the Hippocampal Dentate Gyrus: Critical Role for Topography in Determining Spatiotemporal Network Dynamics

GOAL

This paper describes a million-plus granule cell compartmental model of the rat hippocampal dentate gyrus, including excitatory, perforant path input from the entorhinal cortex, and feedforward and feedback inhibitory input from dentate interneurons.

METHODS

The model includes experimentally determined morphological and biophysical properties of granule cells, together with glutamatergic AMPA-like EPSP and GABAergic GABAA-like IPSP synaptic excitatory and inhibitory inputs, respectively. Each granule cell was composed of approximately 200 compartments having passive and active conductances distributed throughout the somatic and dendritic regions. Modeling excitatory input from the entorhinal cortex was guided by axonal transport studies documenting the topographical organization of projections from subregions of the medial and lateral entorhinal cortex, plus other important details of the distribution of glutamatergic inputs to the dentate gyrus. Information contained within previously published maps of this major hippocampal afferent were systematically converted to scales that allowed the topographical distribution and relative synaptic densities of perforant path inputs to be quantitatively estimated for inclusion in the current model.

RESULTS

Results showed that when medial and lateral entorhinal cortical neurons maintained Poisson random firing, dentate granule cells expressed, throughout the million-cell network, a robust nonrandom pattern of spiking best described as a spatiotemporal "clustering." To identify the network property or properties responsible for generating such firing "clusters," we progressively eliminated from the model key mechanisms, such as feedforward and feedback inhibition, intrinsic membrane properties underlying rhythmic burst firing, and/or topographical organization of entorhinal afferents.

CONCLUSION

Findings conclusively identified topographical organization of inputs as the key element responsible for generating a spatiotemporal distribution of clustered firing. These results uncover a functional organization of perforant path afferents to the dentate gyrus not previously recognized: topography-dependent clusters of granule cell activity as "functional units" or "channels" that organize the processing of entorhinal signals. This modeling study also reveals for the first time how a global signal processing feature of a neural network can evolve from one of its underlying structural characteristics.

Implementation of the excitatory entorhinal-dentate-CA3 topography in a large-scale computational model of the rat hippocampus

The topography, or the anatomical connectivity, of the excitatory entorhinal-dentate-CA3 circuit of the rat hippocampus has been implemented for a large-scale, biologically realistic, computational model of the rat hippocampus. The implementation thus far covers only the excitatory synapses for the principal neurons in the hippocampal subregions. Starting from layer II of the entorhinal cortex, the projection of their perforant path axons has been mapped across the full extent of the dentate gyrus as well as to the CA3. The mossy fiber axon trajectories from the dentate granule cells to the CA3 pyramidal cells have been derived, incorporating the transverse route the fibers take through the CA3c and CA3b and the septo-temporal turn in the CA3a. The extensive arborization of the CA3 pyramidal axons have been modeled using 2-D, skewed Gaussian distributions which have been parametrized to exhibit the differences that exist among the CA3a, CA3b, and CA3c auto-associational projections. Using the limited samples available from the literature, key parameters for each projection have been interpolated as a function of transverse and/or septo-temporal position in order to create a more complete representation of the topography.

A testbed to explore the optimal electrical stimulation parameters for suppressing inter-ictal spikes in human hippocampal slices

New interventions using neuromodulatory devices such as vagus nerve stimulation, deep brain stimulation and responsive neurostimulation are available or under study for the treatment of refractory epilepsy. Since the actual mechanisms of the onset and termination of the seizure are still unclear, most researchers or clinicians determine the optimal stimulation parameters through trial-and-error procedures. It is necessary to further explore what types of electrical stimulation parameters (these may include stimulation frequency, amplitude, duration, interval pattern, and location) constitute a set of optimal stimulation paradigms to suppress seizures. In a previous study, we developed an in vitro epilepsy model using hippocampal slices from patients suffering from mesial temporal lobe epilepsy. Using a planar multi-electrode array system, inter-ictal activity from human hippocampal slices was consistently recorded. In this study, we have further transferred this in vitro seizure model to a testbed for exploring the possible neurostimulation paradigms to inhibit inter-ictal spikes. The methodology used to collect the electrophysiological data, the approach to apply different electrical stimulation parameters to the slices are provided in this paper. The results show that this experimental testbed will provide a platform for testing the optimal stimulation parameters of seizure cessation. We expect this testbed will expedite the process for identifying the most effective parameters, and may ultimately be used to guide programming of new stimulating paradigms for neuromodulatory devices.

Topography-dependent spatio-temporal correlations in the entorhinal-dentate-CA3 circuit in a large-scale computational model of the Rat Hippocampus

The correlation due to different topographies was characterized in a large-scale, biologically-realistic, computational model of the rat hippocampus using a spatio-temporal correlation analysis. The effect of the topographical projection between the following subregions of the hippocampus was investigated: the entorhinal to dentate projection, the entorhinal to CA3 projection, and the mossy fiber to CA3 projection. Through this work, analysis was performed on the individual and combined effects of these projections on the activity of the principal neurons of the dentate gyrus and CA3. The simulations show that uncorrelated input transmitted through the entorhinal-to-dentate or entorhinal-to-CA3 projection causes spatio-temporally correlated activity in the principal neurons that manifest as spike clusters. However, if the mossy fiber system provides uncorrelated input to the CA3, then the CA3 activity remains uncorrelated. When considering the transfer of correlation through the dentate, this analysis suggests that the mossy fiber system do not imbue any correlation to the activity as it propagates from the granule cells of the dentate to the CA3. With the spatio-temporal correlation analysis, the influence of each topographical projection on the transfer of correlation can be investigated as additional subregions and neuron types are added to the large-scale model.

Sparse generalized volterra model of human hippocampal spike train transformation for memory prostheses

In order to build hippocampal prostheses for restoring memory functions, we build multi-input, multi-output (MIMO) nonlinear dynamical models of the human hippocampus. Spike trains are recorded from the hippocampal CA3 and CA1 regions of epileptic patients performing a memory-dependent delayed match-to-sample task. Using CA3 and CA1 spike trains as inputs and outputs respectively, second-order sparse generalized Laguerre-Volterra models are estimated with group lasso and local coordinate descent methods to capture the nonlinear dynamics underlying the spike train transformations. These models can accurately predict the CA1 spike trains based on the ongoing CA3 spike trains and thus will serve as the computational basis of the hippocampal memory prosthesis.

Neurons and networks organizing and sequencing memories

Hippocampal CA1 and CA3 neurons sampled randomly in large numbers in primate brain show conclusive examples of hierarchical encoding of task specific information. Hierarchical encoding allows multi-task utilization of the same hippocampal neural networks via distributed firing between neurons that respond to subsets, attributes or "categories" of stimulus features which can be applied in events in different contexts. In addition, such networks are uniquely adaptable to neural systems unrestricted by rigid synaptic architecture (i.e. columns, layers or "patches") which physically limits the number of possible task-specific interactions between neurons. Also hierarchical encoding is not random; it requires multiple exposures to the same types of relevant events to elevate synaptic connectivity between neurons for different stimulus features that occur in different task-dependent contexts. The large number of cells within associated hierarchical circuits in structures such as hippocampus provides efficient processing of information relevant to common memory-dependent behavioral decisions within different contextual circumstances. This article is part of a Special Issue entitled SI: Brain and Memory.

System identification of point-process neural systems using probability based Volterra kernels

BACKGROUND

Neural information processing involves a series of nonlinear dynamical input/output transformations between the spike trains of neurons/neuronal ensembles. Understanding and quantifying these transformations is critical both for understanding neural physiology such as short-term potentiation and for developing cognitive neural prosthetics.

NEW METHOD

A novel method for estimating Volterra kernels for systems with point-process inputs and outputs is developed based on elementary probability theory. These Probability Based Volterra (PBV) kernels essentially describe the probability of an output spike given q input spikes at various lags t1, t2, …, tq.

RESULTS

The PBV kernels are used to estimate both synthetic systems where ground truth is available and data from the CA3 and CA1 regions rodent hippocampus. The PBV kernels give excellent predictive results in both cases. Furthermore, they are shown to be quite robust to noise and to have good convergence and overfitting properties. Through a slight modification, the PBV kernels are shown to also deal well with correlated point-process inputs.

COMPARISON WITH EXISTING METHODS

The PBV kernels were compared with kernels estimated through least squares estimation (LSE) and through the Laguerre expansion technique (LET). The LSE kernels were shown to fair very poorly with real data due to the large amount of input noise. Although the LET kernels gave the best predictive results in all cases, they require prior parameter estimation. It was shown how the PBV and LET methods can be combined synergistically to maximize performance.

CONCLUSIONS

The PBV kernels provide a novel and intuitive method of characterizing point-process input-output nonlinear systems.

Identification of functional synaptic plasticity from spiking activities using nonlinear dynamical modeling

This paper presents a systems identification approach for studying the long-term synaptic plasticity using natural spiking activities. This approach consists of three modeling steps. First, a multi-input, single-output (MISO), nonlinear dynamical spiking neuron model is formulated to estimate and represent the synaptic strength in means of functional connectivity between input and output neurons. Second, this MISO model is extended to a nonstationary form to track the time-varying properties of the synaptic strength. Finally, a Volterra modeling method is used to extract the synaptic learning rule, e.g., spike-timing-dependent plasticity, for the explanation of the input-output nonstationarity as the consequence of the past input-output spiking patterns. This framework is developed to study the underlying mechanisms of learning and memory formation in behaving animals, and may serve as the computational basis for building the next-generation adaptive cortical prostheses.

Model-based asessment of an in-vivo predictive relationship from CA1 to CA3 in the rodent hippocampus

Although an anatomical connection from CA1 to CA3 via the Entorhinal Cortex (EC) and through backprojecting interneurons has long been known it exist, it has never been examined quantitatively on the single neuron level, in the in-vivo nonpatholgical, nonperturbed brain. Here, single spike activity was recorded using a multi-electrode array from the CA3 and CA1 areas of the rodent hippocampus (N = 7) during a behavioral task. The predictive power from CA3→CA1 and CA1→CA3 was examined by constructing Multivariate Autoregressive (MVAR) models from recorded neurons in both directions. All nonsignificant inputs and models were identified and removed by means of Monte Carlo simulation methods. It was found that 121/166 (73 %) CA3→CA1 models and 96/145 (66 %) CA1→CA3 models had significant predictive power, thus confirming a predictive 'Granger' causal relationship from CA1 to CA3. This relationship is thought to be caused by a combination of truly causal connections such as the CA1→EC→CA3 pathway and common inputs such as those from the Septum. All MVAR models were then examined in the frequency domain and it was found that CA3 kernels had significantly more power in the theta and beta range than those of CA1, confirming CA3's role as an endogenous hippocampal pacemaker.

An in vitro seizure model from human hippocampal slices using multi-electrode arrays

Temporal lobe epilepsy is a neurological condition marked by seizures, typically accompanied by large amplitude synchronous electrophysiological discharges, affecting a variety of mental and physical functions. The neurobiological mechanisms responsible for the onset and termination of seizures are still unclear. While pharmacological therapies can suppress the symptoms of seizures, typically 30% of patients do not respond well to drug control. Unilateral temporal lobectomy, a procedure in which a substantial part of the hippocampal formation and surrounding tissue is removed, is a common surgical treatment for medically refractory epilepsy. In this study, we have developed an in vitro model of epilepsy using human hippocampal slices resected from patients suffering from intractable mesial temporal lobe epilepsy. We show that using a planar multi-electrode array system, spatio-temporal inter-ictal like activity can be consistently recorded in high-potassium (8 mM), low-magnesium (0.25 mM) artificial cerebral spinal fluid with 4-aminopyridine (100 μM) added. The induced epileptiform discharges can be recorded in different subregions of the hippocampus, including dentate, CA1 and subiculum. This new paradigm will allow the study of seizure generation in different subregions of hippocampus simultaneously, as well as propagation of seizure activity throughout the intrinsic circuitry of hippocampus. This experimental model also should provide insights into seizure control and prevention, while providing a platform to develop novel, anti-seizure therapeutics.

Insights on synaptic paired-pulse response using parametric and non-parametric models

Paired-pulse protocol is a well-established stimulation pattern used to characterize short-term changes in synaptic potency. Due to the experimental difficulty in accessing and measuring responses and interactions between subsynaptic elements, understanding the mechanisms that shape synaptic response is extremely challenging. We already proposed to address this issue and gain insights on the matter using a complex integrated modeling platform called EONS (Elementary Objects of the Nervous System). The use of this parametric platform provided us with insightful information on the subsynaptic components and how their interactions shape synaptic dynamics. We herein propose to add and combine a non-parametric model to (i) simplify the modeling framework, the number of underlying parameters and the overall computational complexity while faithfully maintaining the desirable synaptic behavior and (ii) provide a clear and concise framework to characterize AMPA and NMDA contributions to the observed paired-pulse responses.

Extraction and restoration of hippocampal spatial memories with non-linear dynamical modeling

To build a cognitive prosthesis that can replace the memory function of the hippocampus, it is essential to model the input-output function of the damaged hippocampal region, so the prosthetic device can stimulate the downstream hippocampal region, e.g., CA1, with the output signal, e.g., CA1 spike trains, predicted from the ongoing input signal, e.g., CA3 spike trains, and the identified input-output function, e.g., CA3-CA1 model. In order for the downstream region to form appropriate long-term memories based on the restored output signal, furthermore, the output signal should contain sufficient information about the memories that the animal has formed. In this study, we verify this premise by applying regression and classification modelings of the spatio-temporal patterns of spike trains to the hippocampal CA3 and CA1 data recorded from rats performing a memory-dependent delayed non-match-to-sample (DNMS) task. The regression model is essentially the multiple-input, multiple-output (MIMO) non-linear dynamical model of spike train transformation. It predicts the output spike trains based on the input spike trains and thus restores the output signal. In addition, the classification model interprets the signal by relating the spatio-temporal patterns to the memory events. We have found that: (1) both hippocampal CA3 and CA1 spike trains contain sufficient information for predicting the locations of the sample responses (i.e., left and right memories) during the DNMS task; and more importantly (2) the CA1 spike trains predicted from the CA3 spike trains by the MIMO model also are sufficient for predicting the locations on a single-trial basis. These results show quantitatively that, with a moderate number of unitary recordings from the hippocampus, the MIMO non-linear dynamical model is able to extract and restore spatial memory information for the formation of long-term memories and thus can serve as the computational basis of the hippocampal memory prosthesis.

Prefrontal cortical microcircuits bind perception to executive control

During the perception-to-action cycle, our cerebral cortex mediates the interactions between the environment and the perceptual-executive systems of the brain. At the top of the executive hierarchy, prefrontal cortical microcircuits are assumed to bind perceptual and executive control information to guide goal-driven behavior. Here, we tested this hypothesis by comparing simultaneously recorded neuron firing in prefrontal cortical layers and the caudate-putamen of rhesus monkeys, trained in a spatial-versus-object, rule-based match-to-sample task. We found that during the perception and executive selection phases, cell firing in the localized prefrontal layers and caudate-putamen region exhibited similar location preferences on spatial-trials, but less on object- trials. Then, we facilitated the perceptual-executive circuit by stimulating the prefrontal infra-granular-layers with patterns previously derived from supra-granular-layers, and produced stimulation-induced spatial preference in percent correct performance on spatial trials, similar to neural tuning. These results show that inter-laminar prefrontal microcircuits play causal roles to the perception-to-action cycle.

Recording place cells from multiple sub-regions of the rat hippocampus with a customized micro-electrode array

The hippocampus is a subcortical structure which is involved in memory function. There is a considerable amount of evidence available which indicates that the hippocampal system is necessary for effective spatial learning in rodents and short-term topographical memory in human. Recordings of neural activities from the hippocampus of behaving animals can help us to understand how spatial information is encoded and processed by the hippocampus. In this work, we designed a triple-region microelectrode array (MEA) which took into concern the anatomical structures of the rat hippocampus. The array was composed of 16 stainless steel wires which were arranged into three groups that differed in length. Each group targeted one subregion of the hippocampus. The array was chronically implanted into the rat hippocampus through craniotomy. Neural activities were monitored both during the implantation and after recovery. The triple-region MEA was capable of recording unitary activities from multiple subregions of the rat hippocampus and the spatial distribution of firing rates were analyzed while the animal freely explored in the environment.