Kernel current source density method

Community-based reconstruction and simulation of a full-scale model of the rat hippocampus CA1 region

The CA1 region of the hippocampus is one of the most studied regions of the rodent brain, thought to play an important role in cognitive functions such as memory and spatial navigation. Despite a wealth of experimental data on its structure and function, it has been challenging to integrate information obtained from diverse experimental approaches. To address this challenge, we present a community-based, full-scale in silico model of the rat CA1 that integrates a broad range of experimental data, from synapse to network, including the reconstruction of its principal afferents, the Schaffer collaterals, and a model of the effects that acetylcholine has on the system. We tested and validated each model component and the final network model, and made input data, assumptions, and strategies explicit and transparent. The unique flexibility of the model allows scientists to potentially address a range of scientific questions. In this article, we describe the methods used to set up simulations to reproduce in vitro and in vivo experiments. Among several applications in the article, we focus on theta rhythm, a prominent hippocampal oscillation associated with various behavioral correlates and use our computer model to reproduce experimental findings. Finally, we make data, code, and model available through the hippocampushub.eu portal, which also provides an extensive set of analyses of the model and a user-friendly interface to facilitate adoption and usage. This community-based model represents a valuable tool for integrating diverse experimental data and provides a foundation for further research into the complex workings of the hippocampal CA1 region.

Reconciliation of weak pairwise spike-train correlations and highly coherent local field potentials across space

Multi-electrode arrays covering several square millimeters of neural tissue provide simultaneous access to population signals such as extracellular potentials and spiking activity of one hundred or more individual neurons. The interpretation of the recorded data calls for multiscale computational models with corresponding spatial dimensions and signal predictions. Multi-layer spiking neuron network models of local cortical circuits covering about $1\,{\text{mm}^{2}}$ have been developed, integrating experimentally obtained neuron-type-specific connectivity data and reproducing features of observed in-vivo spiking statistics. Local field potentials can be computed from the simulated spiking activity. We here extend a local network and local field potential model to an area of $4\times 4\,{\text{mm}^{2}}$, preserving the neuron density and introducing distance-dependent connection probabilities and conduction delays. We find that the upscaling procedure preserves the overall spiking statistics of the original model and reproduces asynchronous irregular spiking across populations and weak pairwise spike-train correlations in agreement with experimental recordings from sensory cortex. Also compatible with experimental observations, the correlation of local field potential signals is strong and decays over a distance of several hundred micrometers. Enhanced spatial coherence in the low-gamma band around $50\,\text{Hz}$ may explain the recent report of an apparent band-pass filter effect in the spatial reach of the local field potential.

High-Resolution Laminar Identification in Macaque Primary Visual Cortex Using Neuropixels Probes

Laminar electrode arrays allow simultaneous recording of activity of many cortical neurons and assignment to layers using current source density (CSD) analyses. Electrode arrays with 100-micron contact spacing have been used to estimate borders between layer 4 versus superficial or deep layers, but in macaque primary visual cortex (V1) there are far more layers, such as 4A which is only 50-100 microns thick. Neuropixels electrode arrays have 20-micron spacing, and thus could potentially discern thinner layers and more precisely identify laminar borders. Here we show that laminar distributions of CSDs lack consistency and the spatial resolution required for thin layers and accurate layer boundaries. To take full advantage of high density Neuropixels arrays, we have developed approaches based on higher resolution electrical signals and analyses, including spike waveforms and spatial spread, unit density, high-frequency action potential (AP) power spectrum, temporal power change, and coherence spectrum, that afford far higher resolution of laminar distinctions, including the ability to precisely detect the borders of even the thinnest layers of V1.

In vivooptogenetics using a Utah Optrode Array with enhanced light output and spatial selectivity

Objective.Optogenetics allows the manipulation of neural circuitsin vivowith high spatial and temporal precision. However, combining this precision with control over a significant portion of the brain is technologically challenging (especially in larger animal models).Approach.Here, we have developed, optimised, and testedin vivo, the Utah Optrode Array (UOA), an electrically addressable array of optical needles and interstitial sites illuminated by 181μLEDs and used to optogenetically stimulate the brain. The device is specifically designed for non-human primate studies.Main results.Thinning the combinedμLED and needle backplane of the device from 300μm to 230μm improved the efficiency of light delivery to tissue by 80%, allowing lowerμLED drive currents, which improved power management and thermal performance. The spatial selectivity of each site was also improved by integrating an optical interposer to reduce stray light emission. These improvements were achieved using an innovative fabrication method to create an anodically bonded glass/silicon substrate with through-silicon vias etched, forming an optical interposer. Optical modelling was used to demonstrate that the tip structure of the device had a major influence on the illumination pattern. The thermal performance was evaluated through a combination of modelling and experiment, in order to ensure that cortical tissue temperatures did not rise by more than 1 °C. The device was testedin vivoin the visual cortex of macaque expressing ChR2-tdTomato in cortical neurons.Significance.It was shown that the UOA produced the strongest optogenetic response in the region surrounding the needle tips, and that the extent of the optogenetic response matched the predicted illumination profile based on optical modelling-demonstrating the improved spatial selectivity resulting from the optical interposer approach. Furthermore, different needle illumination sites generated different patterns of low-frequency potential activity.

Behavioral control through the direct, focal silencing of neuronal activity

The ability to optically stimulate and inhibit neurons has revolutionized neuroscience research. Here, we present a direct, potent, user-friendly chemical approach for optically silencing neurons. We have rendered saxitoxin (STX), a naturally occurring paralytic agent, transiently inert through chemical protection with a previously undisclosed nitrobenzyl-derived photocleavable group. Exposing the caged toxin, STX-bpc, to a brief (5 ms) pulse of light effects rapid release of a potent STX derivative and transient, spatially precise blockade of voltage-gated sodium channels (NaVs). We demonstrate the efficacy of STX-bpc for parametrically manipulating action potentials in mammalian neurons and brain slice. Additionally, we show the effectiveness of this reagent for silencing neural activity by dissecting sensory-evoked swimming in larval zebrafish. Photo-uncaging of STX-bpc is a straightforward method for non-invasive, reversible, spatiotemporally precise neural silencing without the need for genetic access, thus removing barriers for comparative research.

Local contribution to the somatosensory evoked potentials in rat's thalamus

Local Field Potential (LFP), despite its name, often reflects remote activity. Depending on the orientation and synchrony of their sources, both oscillations and more complex waves may passively spread in brain tissue over long distances and be falsely interpreted as local activity at such distant recording sites. Here we show that the whisker-evoked potentials in the thalamic nuclei are of local origin up to around 6 ms post stimulus, but the later (7-15 ms) wave is overshadowed by a negative component reaching from cortex. This component can be analytically removed and local thalamic LFP can be recovered reliably using Current Source Density analysis. We used model-based kernel CSD (kCSD) method which allowed us to study the contribution of local and distant currents to LFP from rat thalamic nuclei and barrel cortex recorded with multiple, non-linear and non-regular multichannel probes. Importantly, we verified that concurrent recordings from the cortex are not essential for reliable thalamic CSD estimation. The proposed framework can be used to analyze LFP from other brain areas and has consequences for general LFP interpretation and analysis.

In vivo optogenetics using a Utah Optrode Array with enhanced light output and spatial selectivity

Optogenetics allows manipulation of neural circuits in vivo with high spatial and temporal precision. However, combining this precision with control over a significant portion of the brain is technologically challenging (especially in larger animal models). Here, we have developed, optimised, and tested in vivo, the Utah Optrode Array (UOA), an electrically addressable array of optical needles and interstitial sites illuminated by 181 µLEDs and used to optogenetically stimulate the brain. The device is specifically designed for non-human primate studies. Thinning the combined µLED and needle backplane of the device from 300 µm to 230 µm improved the efficiency of light delivery to tissue by 80%, allowing lower µLED drive currents, which improved power management and thermal performance. The spatial selectivity of each site was also improved by integrating an optical interposer to reduce stray light emission. These improvements were achieved using an innovative fabrication method to create an anodically bonded glass/silicon substrate with through-silicon vias etched, forming an optical interposer. Optical modelling was used to demonstrate that the tip structure of the device had a major influence on the illumination pattern. The thermal performance was evaluated through a combination of modelling and experiment, in order to ensure that cortical tissue temperatures did not rise by more than 1°C. The device was tested in vivo in the visual cortex of macaque expressing ChR2-tdTomato in cortical neurons. It was shown that the strongest optogenetic response occurred in the region surrounding the needle tips, and that the extent of the optogenetic response matched the predicted illumination profile based on optical modelling - demonstrating the improved spatial selectivity resulting from the optical interposer approach. Furthermore, different needle illumination sites generated different patterns of low-frequency potential (LFP) activity.

An optrode array for spatiotemporally-precise large-scale optogenetic stimulation of deep cortical layers in non-human primates

Optogenetics has transformed studies of neural circuit function, but remains challenging to apply to non-human primates (NHPs). A major challenge is delivering intense, spatiotemporally-precise, patterned photostimulation across large volumes in deep tissue. Such stimulation is critical, for example, to modulate selectively deep-layer corticocortical feedback circuits. To address this need, we have developed the Utah Optrode Array (UOA), a 10×10 glass needle waveguide array fabricated atop a novel opaque optical interposer, and bonded to an electrically addressable µLED array. In vivo experiments with the UOA demonstrated large-scale, spatiotemporally precise, activation of deep circuits in NHP cortex. Specifically, the UOA permitted both focal (confined to single layers/columns), and widespread (multiple layers/columns) optogenetic activation of deep layer neurons, as assessed with multi-channel laminar electrode arrays, simply by varying the number of activated µLEDs and/or the irradiance. Thus, the UOA represents a powerful optoelectronic device for targeted manipulation of deep-layer circuits in NHP models.

kCSD-python, reliable current source density estimation with quality control

Interpretation of extracellular recordings can be challenging due to the long range of electric field. This challenge can be mitigated by estimating the current source density (CSD). Here we introduce kCSD-python, an open Python package implementing Kernel Current Source Density (kCSD) method and related tools to facilitate CSD analysis of experimental data and the interpretation of results. We show how to counter the limitations imposed by noise and assumptions in the method itself. kCSD-python allows CSD estimation for an arbitrary distribution of electrodes in 1D, 2D, and 3D, assuming distributions of sources in tissue, a slice, or in a single cell, and includes a range of diagnostic aids. We demonstrate its features in a Jupyter Notebook tutorial which illustrates a typical analytical workflow and main functionalities useful in validating analysis results.

Activity in Barrel Cortex Related to Trace Eyeblink Conditioning

In mammals several memory systems are responsible for learning and storage of associative memory. Even apparently simple behavioral tasks, like pavlovian conditioning, have been suggested to engage, for instance, implicit and explicit memory processes. Here, we used single-whisker tactile trace eyeblink conditioning (TTEBC) to investigate learning and its neuronal bases in the mouse barrel column, the primary neocortical tactile representation of one whisker. Behavioral analysis showed that conditioned responses (CRs) are spatially highly restricted; they generalize from the principal whisker only to its direct neighbors. Within the respective neural representation, the principal column and its direct neighbors, spike activity showed a learning-related spike rate suppression starting during the late phase of conditioning stimulus (CS) presentation that was sustained throughout the stimulus-free trace period (Trace). Trial-by-trial analysis showed that learning-related activity was independent from the generation of eyelid movements within a trial, and set in around the steepest part of the learning curve. Optogenetic silencing of responses and their learning-related changes during CS and Trace epochs blocked CR acquisition but not its recall after learning. Silencing during the Trace alone, which carried major parts of the learning-related changes, had no effect. In summary, we demonstrate specific barrel column spike rate plasticity during TTEBC that can be partially decoupled from the CR, the learned eye closure, a hallmark of implicit learning. Our results, thus, point to a possible role of the barrel column in contributing to other kinds of memory as well.

An Optrode Array for Spatiotemporally Precise Large-Scale Optogenetic Stimulation of Deep Cortical Layers in Non-human Primates

Optogenetics has transformed studies of neural circuit function, but remains challenging to apply in non-human primates (NHPs). A major challenge is delivering intense and spatially precise patterned photostimulation across large volumes in deep tissue. Here, we have developed and validated the Utah Optrode Array (UOA) to meet this critical need. The UOA is a 10×10 glass waveguide array bonded to an electrically-addressable μLED array. In vivo electrophysiology and immediate early gene (c-fos) immunohistochemistry demonstrated the UOA allows for large-scale spatiotemporally precise neuromodulation of deep tissue in macaque primary visual cortex. Specifically, the UOA permits both focal (single layers or columns), and large-scale (across multiple layers or columns) photostimulation of deep cortical layers, simply by varying the number of simultaneously activated μLEDs and/or the light irradiance. These results establish the UOA as a powerful tool for studying targeted neural populations within single or across multiple deep layers in complex NHP circuits.

Size tuning of neural response variability in laminar circuits of macaque primary visual cortex

A defining feature of the cortex is its laminar organization, which is likely critical for cortical information processing. For example, visual stimuli of different size evoke distinct patterns of laminar activity. Visual information processing is also influenced by the response variability of individual neurons and the degree to which this variability is correlated among neurons. To elucidate laminar processing, we studied how neural response variability across the layers of macaque primary visual cortex is modulated by visual stimulus size. Our laminar recordings revealed that single neuron response variability and the shared variability among neurons are tuned for stimulus size, and this size-tuning is layer-dependent. In all layers, stimulation of the receptive field (RF) reduced single neuron variability, and the shared variability among neurons, relative to their pre-stimulus values. As the stimulus was enlarged beyond the RF, both single neuron and shared variability increased in supragranular layers, but either did not change or decreased in other layers. Surprisingly, we also found that small visual stimuli could increase variability relative to baseline values. Our results suggest multiple circuits and mechanisms as the source of variability in different layers and call for the development of new models of neural response variability.

Neural correlates in functional brain mapping among breast cancer survivors receiving different chemotherapy regimens: a qEEG/HEG-based investigation

BACKGROUND

Post-chemotherapy cognitive impairment commonly known as 'chemobrain' or 'chemofog' is a well-established clinical disorder affecting various cognitive domains including attention, visuospatial working memory, executive function, etc. Although several studies have confirmed the chemobrain in recent years, scant experiments have evaluated the potential neurotoxicity of different chemotherapy regimens and agents. In this study, we aimed to evaluate the extent of attention deficits, one of the commonly affected cognitive domains, among breast cancer patients treated with different chemotherapy regimens through neuroimaging techniques.

METHODS

Breast cancer patients treated with two commonly prescribed chemotherapy regimens, Adriamycin, Cyclophosphamide and Taxol and Taxotere, Adriamycin and Cyclophosphamide, and healthy volunteers were recruited. Near-infrared hemoencephalography and quantitative electroencephalography assessments were recorded for each participant at rest and during task performance to compare the functional cortical changes associated with each chemotherapy regimen.

RESULTS

Although no differences were observed in hemoencephalography results across groups, the quantitative electroencephalography analysis revealed increased power of high alpha/low beta in left fronto-centro-parietal regions involved in dorsal and ventral attention networks in the Adriamycin, Cyclophosphamide and Taxol-treated group compared with the Taxotere, Adriamycin and Cyclophosphamide and control group. The Adriamycin, Cyclophosphamide and Taxol-treated cases had the highest current source density values in dorsal attention network and ventral attention network and ventral attention network-related centers in 10 and 15 Hz associated with the lowest Z-scored Fast Fourier Transform coherence in the mentioned regions.

CONCLUSIONS

The negatively affected neurocognitive profile in breast cancer patients treated with the Adriamycin, Cyclophosphamide and Taxol regimen proposes presumably neurotoxic sequelae of this chemotherapy regimen as compared with the Taxotere, Adriamycin and Cyclophosphamide regimen.

Implementation of biohybrid olfactory bulb on a high-density CMOS-chip to reveal large-scale spatiotemporal circuit information

Large-scale multi-site biosensors are essential to probe the olfactory bulb (OB) circuitry for understanding the spatiotemporal dynamics of simultaneous discharge patterns. Current ex-vivo biosensing techniques are limited to recording a small set of neurons and cannot provide an adequate resolution, which hinders revealing the fast dynamic underlying the information coding mechanisms in the OB circuit. Here, we demonstrate a novel biohybrid OB-CMOS biosensing platform to decipher the cross-scale dynamics of the OB electrogenesis and quantify the distinct neuronal coding properties. The approach with 4096-microelectrodes offers a non-invasive, label-free, bioelectrical imaging to decode simultaneous firing patterns from thousands of connected neuronal ensembles in acute OB slices. The platform can measure spontaneous and drug-induced extracellular field potential activity with substantially improved spatiotemporal resolution over conventional OB-based biosensors. Also, we employ our OB-CMOS recordings to perform multidimensional analysis to instantiate specific neurophysiological metrics underlying the olfactory spatiotemporal coding that emerged from the OB interconnected layers. Our results delineate the computational implications of large-scale activity patterns in functional olfactory processing. The systematic interplay of the experimental CMOS-base platform architecture and the high-content characterization of the olfactory circuit with various computational analyses endow significant functional interrogations of the OB information processing, high-spatiotemporal connectivity mapping, and global circuit dynamics. Thus, our study can inspire the design of advanced biomimetic olfactory-based biosensors and neuromorphic approaches for diagnostic biomarkers and drug discovery applications.

Cross-population coupling of neural activity based on Gaussian process current source densities

Because local field potentials (LFPs) arise from multiple sources in different spatial locations, they do not easily reveal coordinated activity across neural populations on a trial-to-trial basis. As we show here, however, once disparate source signals are decoupled, their trial-to-trial fluctuations become more accessible, and cross-population correlations become more apparent. To decouple sources we introduce a general framework for estimation of current source densities (CSDs). In this framework, the set of LFPs result from noise being added to the transform of the CSD by a biophysical forward model, while the CSD is considered to be the sum of a zero-mean, stationary, spatiotemporal Gaussian process, having fast and slow components, and a mean function, which is the sum of multiple time-varying functions distributed across space, each varying across trials. We derived biophysical forward models relevant to the data we analyzed. In simulation studies this approach improved identification of source signals compared to existing CSD estimation methods. Using data recorded from primate auditory cortex, we analyzed trial-to-trial fluctuations in both steady-state and task-evoked signals. We found cortical layer-specific phase coupling between two probes and showed that the same analysis applied directly to LFPs did not recover these patterns. We also found task-evoked CSDs to be correlated across probes, at specific cortical depths. Using data from Neuropixels probes in mouse visual areas, we again found evidence for depth-specific phase coupling of primary visual cortex and lateromedial area based on the CSDs.

What we can and what we cannot see with extracellular multielectrodes

Extracellular recording is an accessible technique used in animals and humans to study the brain physiology and pathology. As the number of recording channels and their density grows it is natural to ask how much improvement the additional channels bring in and how we can optimally use the new capabilities for monitoring the brain. Here we show that for any given distribution of electrodes we can establish exactly what information about current sources in the brain can be recovered and what information is strictly unobservable. We demonstrate this in the general setting of previously proposed kernel Current Source Density method and illustrate it with simplified examples as well as using evoked potentials from the barrel cortex obtained with a Neuropixels probe and with compatible model data. We show that with conceptual separation of the estimation space from experimental setup one can recover sources not accessible to standard methods.

Every set of measurements is a window into reality rendering its incomplete or distorted picture. It is often difficult to relate the obtained representation of the world to underlying ground truth. Here we show, for brain electrophysiology, for arbitrary experimental setup (distribution of electrodes), and arbitrary analytical setup (function space of current source densities), that one can identify distributions of current sources which can be recovered precisely, and those which are invisible in the system. This shows what is and what is not observable in the studied system for a given setup, allows to improve the analysis results by modifying analytical setup, and facilitates interpretation of the measured sets of LFP, ECoG and EEG recordings. In particular we show that with conceptual separation of the estimation space from experimental setup one can recover source distributions not accessible to standard methods.

Network and synaptic mechanisms underlying high frequency oscillations in the rat and cat olfactory bulb under ketamine-xylazine anesthesia

Wake-related ketamine-dependent high frequency oscillations (HFO) can be recorded in local field potentials (LFP) from cortical and subcortical regions in rodents. The mechanisms underlying their generation and occurrence in higher mammals are unclear. Unfortunately, anesthetic doses of pure ketamine attenuate HFO, which has precluded their investigation under anesthesia. Here, we show ketamine-xylazine (KX) anesthesia is associated with a prominent 80–130 Hz rhythm in the olfactory bulb (OB) of rats, whereas 30–65 Hz gamma power is diminished. Simultaneous LFP and thermocouple recordings revealed the 80–130 Hz rhythm was dependent on nasal respiration. This rhythm persisted despite surgical excision of the piriform cortex. Silicon probes spanning the dorsoventral aspect of the OB revealed this rhythm was strongest in ventral areas and associated with microcurrent sources about the mitral layer. Pharmacological microinfusion studies revealed dependency on excitatory-inhibitory synaptic activity, but not gap junctions. Finally, a similar rhythm occurred in the OB of KX-anesthetized cats, which shared key features with our rodent studies. We conclude that the activity we report here is driven by nasal airflow, local excitatory-inhibitory interactions, and conserved in higher mammals. Additionally, KX anesthesia is a convenient model to investigate further the mechanisms underlying wake-related ketamine-dependent HFO.

Brain Modeling ToolKit: An open source software suite for multiscale modeling of brain circuits

Experimental studies in neuroscience are producing data at a rapidly increasing rate, providing exciting opportunities and formidable challenges to existing theoretical and modeling approaches. To turn massive datasets into predictive quantitative frameworks, the field needs software solutions for systematic integration of data into realistic, multiscale models. Here we describe the Brain Modeling ToolKit (BMTK), a software suite for building models and performing simulations at multiple levels of resolution, from biophysically detailed multi-compartmental, to point-neuron, to population-statistical approaches. Leveraging the SONATA file format and existing software such as NEURON, NEST, and others, BMTK offers a consistent user experience across multiple levels of resolution. It permits highly sophisticated simulations to be set up with little coding required, thus lowering entry barriers to new users. We illustrate successful applications of BMTK to large-scale simulations of a cortical area. BMTK is an open-source package provided as a resource supporting modeling-based discovery in the community.

Recent Advances in Electrical Neural Interface Engineering: Minimal Invasiveness, Longevity, and Scalability

Electrical neural interfaces serve as direct communication pathways that connect the nervous system with the external world. Technological advances in this domain are providing increasingly more powerful tools to study, restore, and augment neural functions. Yet, the complexities of the nervous system give rise to substantial challenges in the design, fabrication, and system-level integration of these functional devices. In this review, we present snapshots of the latest progresses in electrical neural interfaces, with an emphasis on advances that expand the spatiotemporal resolution and extent of mapping and manipulating brain circuits. We include discussions of large-scale, long-lasting neural recording; wireless, miniaturized implants; signal transmission, amplification, and processing; as well as the integration of interfaces with optical modalities. We outline the background and rationale of these developments and share insights into the future directions and new opportunities they enable.

Application of a flexible polymer microECoG array to map functional coherence in schizophrenia model

Anatomically, connections form the fundamental brain network, functionally the different types of oscillatory electric activities are creating a temporarily connected fraction of the anatomical connectome generating an output to the motor system. Schizophrenia can be considered as a connectome disease, in which the sensory input generates a schizophrenia specific temporary connectome and the signal processing becomes diseased showing hallucinations and adverse behavioral reactions. In this work, flexible, 32-channel polymer microelectrode arrays fabricated by the authors are used to map the functional coherence on large cortical areas during physiological activities in a schizophrenia model in rats.-Fabrication of a flexible microECoG array is shown.-Protocol to use a flexible microECoG is demonstrated to characterize connectome diseases in rats.-Customized method to analyze the functional coherence between different cortical areas during visually evoked potential is detailed.-R-based implementation of the analysis method is presented.

Comparing Open-Source Toolboxes for Processing and Analysis of Spike and Local Field Potentials Data

Analysis of spike and local field potential (LFP) data is an essential part of neuroscientific research. Today there exist many open-source toolboxes for spike and LFP data analysis implementing various functionality. Here we aim to provide a practical guidance for neuroscientists in the choice of an open-source toolbox best satisfying their needs. We overview major open-source toolboxes for spike and LFP data analysis as well as toolboxes with tools for connectivity analysis, dimensionality reduction and generalized linear modeling. We focus on comparing toolboxes functionality, statistical and visualization tools, documentation and support quality. To give a better insight, we compare and illustrate functionality of the toolboxes on open-access dataset or simulated data and make corresponding MATLAB scripts publicly available.

State-aware detection of sensory stimuli in the cortex of the awake mouse

Cortical responses to sensory inputs vary across repeated presentations of identical stimuli, but how this trial-to-trial variability impacts detection of sensory inputs is not fully understood. Using multi-channel local field potential (LFP) recordings in primary somatosensory cortex (S1) of the awake mouse, we optimized a data-driven cortical state classifier to predict single-trial sensory-evoked responses, based on features of the spontaneous, ongoing LFP recorded across cortical layers. Our findings show that, by utilizing an ongoing prediction of the sensory response generated by this state classifier, an ideal observer improves overall detection accuracy and generates robust detection of sensory inputs across various states of ongoing cortical activity in the awake brain, which could have implications for variability in the performance of detection tasks across brain states.

Establishing the link between neural activity and behavior is a central goal of neuroscience. One context in which to examine this link is in a sensory detection task, in which an animal is trained to report the presence of a barely perceptible sensory stimulus. In such tasks, both sensory responses in the brain and behavioral responses are highly variable. A simple hypothesis, originating in signal detection theory, is that perceived inputs generate neural activity that cross some threshold for detection. According to this hypothesis, sensory response variability would predict behavioral variability, but previous studies have not born out this prediction. Further complicating the picture, sensory response variability is partially dependent on the ongoing state of cortical activity, and we wondered whether this could resolve the mismatch between response variability and behavioral variability. Here, we use a computational approach to study an adaptive observer that utilizes an ongoing prediction of sensory responsiveness to detect sensory inputs. This observer has higher overall accuracy than the standard ideal observer. Moreover, because of the adaptation, the observer breaks the direct link between neural and behavioral variability, which could resolve discrepancies arising in past studies. We suggest new experiments to test our theory.

The olfactory bulb is a source of high-frequency oscillations (130-180 Hz) associated with a subanesthetic dose of ketamine in rodents

High-frequency neuronal population oscillations (HFO, 130-180 Hz) are robustly potentiated by subanesthetic doses of ketamine. This frequency band has been recorded in functionally and neuroanatomically diverse cortical and subcortical regions, notably ventral striatal areas. However, the locus of generation remains largely unknown. There is compelling evidence that olfactory regions can drive oscillations in distant areas. Here we tested the hypothesis that the olfactory bulb (OB) is a locus for the generation of HFO following a subanesthetic dose of ketamine. The effect of ketamine on the electrophysiological activity of the OB and ventral striatum of male Wistar rats was examined using field potential and unit recordings, local inhibition, naris blockade, current source density and causality estimates. Ketamine-HFO was of larger magnitude and was phase-advanced in the OB relative to ventral striatum. Granger causality analysis was consistent with the OB as the source of HFO. Unilateral local inhibition of the OB and naris blockade both attenuated HFO recorded locally and in the ventral striatum. Within the OB, current source density analysis revealed HFO current dipoles close to the mitral layer and unit firing of mitral/tufted cells was phase locked to HFO. Our results reveal the OB as a source of ketamine-HFO which can contribute to HFO in the ventral striatum, known to project diffusely to many other brain regions. These findings provide a new conceptual understanding on how changes in olfactory system function may have implications for neurological disorders involving NMDA receptor dysfunction such as schizophrenia and depression.

Distinct Laminar Processing of Local and Global Context in Primate Primary Visual Cortex

Visual perception is affected by spatial context. In visual cortex, neuronal responses to stimuli inside the receptive field (RF) are suppressed by stimuli in the RF surround. To understand the circuits and cortical layers processing spatial context, we simultaneously recorded across all layers of macaque primary visual cortex while presenting stimuli at increasing distances from the recorded cells' RF. We find that near versus far-surround stimuli activate distinct layers, thus revealing unique laminar contributions to the processing of local and global spatial context. Stimuli in the near-surround evoke the earliest subthreshold responses in superficial and upper-deep layers, and earliest suppression of spiking responses in superficial layers. Conversely, far-surround stimuli evoke the earliest subthreshold responses in feedback-recipient layer 1 and lower-deep layers, and earliest suppression of spiking responses almost simultaneously in all layers, except 4C, where suppression emerges last. Our results suggest distinct circuits for local and global signal integration.

Current Source Density Analysis of Electroantennogram Recordings: A Tool for Mapping the Olfactory Response in an Insect Antenna

The set of chemosensory receptors expressed by the olfactory receptor neurons lying in an insect's antennae and maxillary palps define the ability of this insect to perceive the volatile chemicals of its environment. The main two electrophysiological methods of antennal recordings for studying the range of chemicals that activate chemosensory receptors have limitations. Single-sensillum recording (SSR) samples a subset of olfactory receptor neurons and therefore does not reveal the full capacity of an insect to perceive an odor. Electroantennography (EAG), even if less resolutive than SSRs, is sometimes preferred since it samples the activity of a large number of the olfactory receptor neurons. But, at least in flies, the amplitude of the EAG signal is not directly correlated with the degree of sensitivity of the insect to the olfactory compound. Such dual methodology was also used to study mammalian brains, and the current source density (CSD) analysis was developed to bridge the gap between the cellular and the population recordings. This paper details the use of a similar approach adapted to the study of olfactory responses within insects with bulbous antennae. The EAG was recorded at multiple antennal positions and the CSD that generates the EAG potentials were estimated. The method measures the activation of olfactory receptor neurons (ORNs) across the antennae and thus it quantifies the olfactory sensitivity of the insect. It allows a rapid mapping of olfactory responses and thus can be used to guide further SSRs or to determine that two chemicals are detected by independent ORNs. This study further explored biases resulting from a limited number of recording positions or from an approximation of the antennal geometry that should be considered for interpreting the CSD maps. It also shows that the CSD analysis of EAGs is compatible with a gas chromatograph stimulator for analyzing the response to complex odors. Finally, I discuss the origin of the EAG signal in light of the CSD theory.

Deciphering the Contribution of Oriens-Lacunosum/Moleculare (OLM) Cells to Intrinsic θ Rhythms Using Biophysical Local Field Potential (LFP) Models

Oscillations in local field potentials (LFPs) are prevalent and contribute to brain function. An understanding of the cellular correlates and pathways affecting LFPs is needed, but many overlapping pathways in vivo make this difficult to achieve. A prevalent LFP rhythm in the hippocampus associated with memory processing and spatial navigation is the θ (3–12 Hz) oscillation. θ rhythms emerge intrinsically in an in vitro whole hippocampus preparation and this reduced preparation makes it possible to assess the contribution of different cell types to LFP generation. We focus on oriens-lacunosum/moleculare (OLM) cells as a major class of interneurons in the hippocampus. OLM cells can influence pyramidal (PYR) cells through two distinct pathways: by direct inhibition of PYR cell distal dendrites, and by indirect disinhibition of PYR cell proximal dendrites. We use previous inhibitory network models and build biophysical LFP models using volume conductor theory. We examine the effect of OLM cells to ongoing intrinsic LFP θ rhythms by directly comparing our model LFP features with experiment. We find that OLM cell inputs regulate the robustness of LFP responses without affecting their average power and that this robust response depends on coactivation of distal inhibition and basal excitation. We use our models to estimate the spatial extent of the region generating LFP θ rhythms, leading us to predict that about 22,000 PYR cells participate in intrinsic θ generation. Besides obtaining an understanding of OLM cell contributions to intrinsic LFP θ rhythms, our work can help decipher cellular correlates of in vivo LFPs.

waveCSD: A method for estimating transmembrane currents originated from propagating neuronal activity in the neocortex: Application to study cortical spreading depression

BACKGROUND

Recent years have witnessed an upsurge in the development of methods for estimating current source densities (CSDs) in the neocortical tissue from their recorded local field potential (LFP) reflections using microelectrode arrays. Among these, methods utilizing linear arrays work under the assumption that CSDs vary as a function of cortical depth; whereas they are constant in the tangential direction, infinitely or in a confined cylinder. This assumption is violated in the analysis of neuronal activity propagating along the neocortical sheet, e.g. propagation of alpha waves or cortical spreading depression.

NEW METHOD

Here, we developed a novel mathematical method (waveCSD) for CSD analysis of LFPs associated with a planar wave of neocortical neuronal activity propagating at a constant velocity towards a linear probe.

RESULTS

Results show that the algorithm is robust to the presence of noise in LFP data and uncertainties in knowledge of propagation velocity. Also, results show high level of accuracy of the method in a wide range of electrode resolutions. Using in vivo experimental recordings from the rat neocortex, we employed waveCSD to characterize transmembrane currents associated with cortical spreading depressions.

COMPARISON WITH EXISTING METHOD(S)

Simulation results indicate that waveCSD has a significantly higher reconstruction accuracy compared to the widely-used inverse CSD method (iCSD), and the regularized kernel CSD method (kCSD), in the analysis of CSDs originating from propagating neuronal activity.

CONCLUSIONS

The waveCSD method provides a theoretical platform for estimation of transmembrane currents from their LFPs in experimental paradigms involving wave propagation.

Investigating large-scale brain dynamics using field potential recordings: analysis and interpretation

New technologies to record electrical activity from the brain on a massive scale offer tremendous opportunities for discovery. Electrical measurements of large-scale brain dynamics, termed field potentials, are especially important to understanding and treating the human brain. Here, our goal is to provide best practices on how field potential recordings (electroencephalograms, magnetoencephalograms, electrocorticograms and local field potentials) can be analyzed to identify large-scale brain dynamics, and to highlight critical issues and limitations of interpretation in current work. We focus our discussion of analyses around the broad themes of activation, correlation, communication and coding. We provide recommendations for interpreting the data using forward and inverse models. The forward model describes how field potentials are generated by the activity of populations of neurons. The inverse model describes how to infer the activity of populations of neurons from field potential recordings. A recurring theme is the challenge of understanding how field potentials reflect neuronal population activity given the complexity of the underlying brain systems.

Top-down feedback controls spatial summation and response amplitude in primate visual cortex

Sensory information travels along feedforward connections through a hierarchy of cortical areas, which, in turn, send feedback connections to lower-order areas. Feedback has been implicated in attention, expectation, and sensory context, but the mechanisms underlying these diverse feedback functions are unknown. Using specific optogenetic inactivation of feedback connections from the secondary visual area (V2), we show how feedback affects neural responses in the primate primary visual cortex (V1). Reducing feedback activity increases V1 cells' receptive field (RF) size, decreases their responses to stimuli confined to the RF, and increases their responses to stimuli extending into the proximal surround, therefore reducing surround suppression. Moreover, stronger reduction of V2 feedback activity leads to progressive increase in RF size and decrease in response amplitude, an effect predicted by a recurrent network model. Our results indicate that feedback modulates RF size, surround suppression and response amplitude, similar to the modulatory effects of visual spatial attention.

Linear distributed source modeling of local field potentials recorded with intra-cortical electrode arrays

Planar intra-cortical electrode (Utah) arrays provide a unique window into the spatial organization of cortical activity. Reconstruction of the current source density (CSD) underlying such recordings, however, requires “inverting” Poisson’s equation. For inter-laminar recordings, this is commonly done by the CSD method, which consists in taking the second-order spatial derivative of the recorded local field potentials (LFPs). Although the CSD method has been tremendously successful in mapping the current generators underlying inter-laminar LFPs, its application to planar recordings is more challenging. While for inter-laminar recordings the CSD method seems reasonably robust against violations of its assumptions, is it unclear as to what extent this holds for planar recordings. One of the objectives of this study is to characterize the conditions under which the CSD method can be successfully applied to Utah array data. Using forward modeling, we find that for spatially coherent CSDs, the CSD method yields inaccurate reconstructions due to volume-conducted contamination from currents in deeper cortical layers. An alternative approach is to “invert” a constructed forward model. The advantage of this approach is that any a priori knowledge about the geometrical and electrical properties of the tissue can be taken into account. Although several inverse methods have been proposed for LFP data, the applicability of existing electroencephalographic (EEG) and magnetoencephalographic (MEG) inverse methods to LFP data is largely unexplored. Another objective of our study therefore, is to assess the applicability of the most commonly used EEG/MEG inverse methods to Utah array data. Our main conclusion is that these inverse methods provide more accurate CSD reconstructions than the CSD method. We illustrate the inverse methods using event-related potentials recorded from primary visual cortex of a macaque monkey during a motion discrimination task.

Revealing the distribution of transmembrane currents along the dendritic tree of a neuron from extracellular recordings

Revealing the current source distribution along the neuronal membrane is a key step on the way to understanding neural computations; however, the experimental and theoretical tools to achieve sufficient spatiotemporal resolution for the estimation remain to be established. Here, we address this problem using extracellularly recorded potentials with arbitrarily distributed electrodes for a neuron of known morphology. We use simulations of models with varying complexity to validate the proposed method and to give recommendations for experimental applications. The method is applied to in vitro data from rat hippocampus.

Collection of Simulated Data from a Thalamocortical Network Model

A major challenge in experimental data analysis is the validation of analytical methods in a fully controlled scenario where the justification of the interpretation can be made directly and not just by plausibility. In some sciences, this could be a mathematical proof, yet biological systems usually do not satisfy assumptions of mathematical theorems. One solution is to use simulations of realistic models to generate ground truth data. In neuroscience, creating such data requires plausible models of neural activity, access to high performance computers, expertise and time to prepare and run the simulations, and to process the output. To facilitate such validation tests of analytical methods we provide rich data sets including intracellular voltage traces, transmembrane currents, morphologies, and spike times. Moreover, these data can be used to study the effects of different tissue models on the measurement. The data were generated using the largest publicly available multicompartmental model of thalamocortical network (Traub et al., Journal of Neurophysiology, 93(4), 2194–2232 (Traub et al. 2005)), with activity evoked by different thalamic stimuli.

Focal Local Field Potential Signature of the Single-Axon Monosynaptic Thalamocortical Connection

A resurgence has taken place in recent years in the use of the extracellularly recorded local field potential (LFP) to investigate neural network activity. To probe monosynaptic thalamic activation of cortical postsynaptic target cells, so called spike-trigger-averaged LFP (stLFP) signatures have been measured. In these experiments, the cortical LFP is measured by multielectrodes covering several cortical lamina and averaged on spontaneous spikes of thalamocortical (TC) cells. Using a well established forward-modeling scheme, we investigated the biophysical origin of this stLFP signature with simultaneous synaptic activation of cortical layer-4 neurons, mimicking the effect of a single afferent spike from a single TC neuron. Constrained by previously measured intracellular responses of the main postsynaptic target cell types and with biologically plausible assumptions regarding the spatial distribution of thalamic synaptic inputs into layer 4, the model predicted characteristic contributions to monosynaptic stLFP signatures both for the regular-spiking (RS) excitatory neurons and the fast-spiking (FS) inhibitory interneurons. In particular, the FS cells generated stLFP signatures of shorter temporal duration than the RS cells. Added together, a sum of the stLFP signatures of these two principal synaptic targets of TC cells were observed to resemble experimentally measured stLFP signatures. Outside the volume targeted by TC afferents, the resulting postsynaptic LFP signals were found to be sharply attenuated. This implies that such stLFP signatures provide a very local measure of TC synaptic activation, and that newly developed inverse current-source density (CSD)-estimation methods are needed for precise assessment of the underlying spatiotemporal CSD profiles.SIGNIFICANCE STATEMENT Despite its long history and prevalent use, the proper interpretation of the extracellularly recorded local field potential (LFP) is still not fully established. Here we investigate by biophysical modeling the origin of the focal LFP signature of the single-axon monosynaptic thalamocortical connection as measured by spike-trigger-averaging of cortical LFPs on spontaneous spikes of thalamocortical neurons. We find that this LFP signature is well accounted for by a model assuming thalamic projections to two cortical layer-4 cell populations: one excitatory (putatively regular-spiking cells) and one inhibitory (putatively fast-spiking cells). The LFP signature is observed to decay sharply outside the cortical region receiving the thalamocortical projection, implying that it indeed provides a very local measure of thalamocortical synaptic activation.

From Maxwell's equations to the theory of current-source density analysis

Despite the widespread use of current‐source density (CSD) analysis of extracellular potential recordings in the brain, the physical mechanisms responsible for the generation of the signal are still debated. While the extracellular potential is thought to be exclusively generated by the transmembrane currents, recent studies suggest that extracellular diffusive, advective and displacement currents—traditionally neglected—may also contribute considerably toward extracellular potential recordings. Here, we first justify the application of the electro‐quasistatic approximation of Maxwell's equations to describe the electromagnetic field of physiological origin. Subsequently, we perform spatial averaging of currents in neural tissue to arrive at the notion of the CSD and derive an equation relating it to the extracellular potential. We show that, in general, the extracellular potential is determined by the CSD of membrane currents as well as the gradients of the putative extracellular diffusion current. The diffusion current can contribute significantly to the extracellular potential at frequencies less than a few Hertz; in which case it must be subtracted to obtain correct CSD estimates. We also show that the advective and displacement currents in the extracellular space are negligible for physiological frequencies while, within cellular membrane, displacement current contributes toward the CSD as a capacitive current. Taken together, these findings elucidate the relationship between electric currents and the extracellular potential in brain tissue and form the necessary foundation for the analysis of extracellular recordings.

Ion diffusion may introduce spurious current sources in current-source density (CSD) analysis

Standard CSD analysis does not account for ionic diffusion. Using biophysically realistic computer simulations, we show that unaccounted-for diffusive currents can lead to the prediction of spurious current sources. This finding may be of strong interest for in vivo electrophysiologists doing extracellular recordings in general, and CSD analysis in particular.

Current-source density (CSD) analysis is a well-established method for analyzing recorded local field potentials (LFPs), that is, the low-frequency part of extracellular potentials. Standard CSD theory is based on the assumption that all extracellular currents are purely ohmic, and thus neglects the possible impact from ionic diffusion on recorded potentials. However, it has previously been shown that in physiological conditions with large ion-concentration gradients, diffusive currents can evoke slow shifts in extracellular potentials. Using computer simulations, we here show that diffusion-evoked potential shifts can introduce errors in standard CSD analysis, and can lead to prediction of spurious current sources. Further, we here show that the diffusion-evoked prediction errors can be removed by using an improved CSD estimator which accounts for concentration-dependent effects.

NEW & NOTEWORTHY Standard CSD analysis does not account for ionic diffusion. Using biophysically realistic computer simulations, we show that unaccounted-for diffusive currents can lead to the prediction of spurious current sources. This finding may be of strong interest for in vivo electrophysiologists doing extracellular recordings in general, and CSD analysis in particular.

Dipolar extracellular potentials generated by axonal projections

Extracellular field potentials (EFPs) are an important source of information in neuroscience, but their physiological basis is in many cases still a matter of debate. Axonal sources are typically discounted in modeling and data analysis because their contributions are assumed to be negligible. Here, we established experimentally and theoretically that contributions of axons to EFPs can be significant. Modeling action potentials propagating along axons, we showed that EFPs were prominent in the presence of terminal zones where axons branch and terminate in close succession, as found in many brain regions. Our models predicted a dipolar far field and a polarity reversal at the center of the terminal zone. We confirmed these predictions using EFPs from the barn owl auditory brainstem where we recorded in nucleus laminaris using a multielectrode array. These results demonstrate that axonal terminal zones can produce EFPs with considerable amplitude and spatial reach.

Spatial Localization of Sources in the Rat Subthalamic Motor Region Using an Inverse Current Source Density Method

Objective: In this study we introduce the use of the current source density (CSD) method as a way to visualize the spatial organization of evoked responses in the rat subthalamic nucleus (STN) at fixed time stamps resulting from motor cortex stimulation. This method offers opportunities to visualize neuronal input and study the relation between the synaptic input and the neural output of neural populations.

Approach: Motor cortex evoked local field potentials and unit activity were measured in the subthalamic region, with a 3D measurement grid consisting of 320 measurement points and high spatial resolution. This allowed us to visualize the evoked synaptic input by estimating the current source density (CSD) from the measured local field potentials, using the inverse CSD method. At the same time, the neuronal output of the cells within the grid is assessed by calculating post stimulus time histograms.

Main results: The CSD method resulted in clear and distinguishable sources and sinks of the neuronal input activity in the STN after motor cortex stimulation. We showed that the center of the synaptic input of the STN from the motor cortex is located dorsal to the input from globus pallidus.

Significance: For the first time we have performed CSD analysis on motor cortex stimulation evoked LFP responses in the rat STN as a proof of principle. Our results suggest that the CSD method can be used to gain new insights into the spatial extent of synaptic pathways in brain structures.

1D Current Source Density (CSD) Estimation in Inverse Theory: A Unified Framework for Higher-Order Spectral Regularization of Quadrature and Expansion-Type CSD Methods

Estimation of current source density (CSD) from the low-frequency part of extracellular electric potential recordings is an unstable linear inverse problem. To make the estimation possible in an experimental setting where recordings are contaminated with noise, it is necessary to stabilize the inversion. Here we present a unified framework for zero- and higher-order singular-value-decomposition (SVD)-based spectral regularization of 1D (linear) CSD estimation from local field potentials. The framework is based on two general approaches commonly employed for solving inverse problems: quadrature and basis function expansion. We first show that both inverse CSD (iCSD) and kernel CSD (kCSD) fall into the category of basis function expansion methods. We then use these general categories to introduce two new estimation methods, quadrature CSD (qCSD), based on discretizing the CSD integral equation with a chosen quadrature rule, and representer CSD (rCSD), an even-determined basis function expansion method that uses the problem's data kernels (representers) as basis functions. To determine the best candidate methods to use in the analysis of experimental data, we compared the different methods on simulations under three regularization schemes (Tikhonov, tSVD, and dSVD), three regularization parameter selection methods (NCP, L-curve, and GCV), and seven different a priori spatial smoothness constraints on the CSD distribution. This resulted in a comparison of 531 estimation schemes. We evaluated the estimation schemes according to their source reconstruction accuracy by testing them using different simulated noise levels, lateral source diameters, and CSD depth profiles. We found that ranking schemes according to the average error over all tested conditions results in a reproducible ranking, where the top schemes are found to perform well in the majority of tested conditions. However, there is no single best estimation scheme that outperforms all others under all tested conditions. The unified framework we propose expands the set of available estimation methods, provides increased flexibility for 1D CSD estimation in noisy experimental conditions, and allows for a meaningful comparison between estimation schemes.

Generalized Laminar Population Analysis (gLPA) for Interpretation of Multielectrode Data from Cortex

Laminar population analysis (LPA) is a method for analysis of electrical data recorded by linear multielectrodes passing through all lamina of cortex. Like principal components analysis (PCA) and independent components analysis (ICA), LPA offers a way to decompose the data into contributions from separate cortical populations. However, instead of using purely mathematical assumptions in the decomposition, LPA is based on physiological constraints, i.e., that the observed LFP (low-frequency part of signal) is driven by action-potential firing as observed in the MUA (multi-unit activity; high-frequency part of the signal). In the presently developed generalized laminar population analysis (gLPA) the set of basis functions accounting for the LFP data is extended compared to the original LPA, thus allowing for a better fit of the model to experimental data. This enhances the risk for overfitting, however, and we therefore tested various versions of gLPA on virtual LFP data in which we knew the ground truth. These synthetic data were generated by biophysical forward-modeling of electrical signals from network activity in the comprehensive, and well-known, thalamocortical network model developed by Traub and coworkers. The results for the Traub model imply that while the laminar components extracted by the original LPA method overall are in fair agreement with the ground-truth laminar components, the results may be improved by use of gLPA method with two (gLPA-2) or even three (gLPA-3) postsynaptic LFP kernels per laminar population.

Modelling and Analysis of Electrical Potentials Recorded in Microelectrode Arrays (MEAs)

Microelectrode arrays (MEAs), substrate-integrated planar arrays of up to thousands of closely spaced metal electrode contacts, have long been used to record neuronal activity in in vitro brain slices with high spatial and temporal resolution. However, the analysis of the MEA potentials has generally been mainly qualitative. Here we use a biophysical forward-modelling formalism based on the finite element method (FEM) to establish quantitatively accurate links between neural activity in the slice and potentials recorded in the MEA set-up. Then we develop a simpler approach based on the method of images (MoI) from electrostatics, which allows for computation of MEA potentials by simple formulas similar to what is used for homogeneous volume conductors. As we find MoI to give accurate results in most situations of practical interest, including anisotropic slices covered with highly conductive saline and MEA-electrode contacts of sizable physical extensions, a Python software package (ViMEAPy) has been developed to facilitate forward-modelling of MEA potentials generated by biophysically detailed multicompartmental neurons. We apply our scheme to investigate the influence of the MEA set-up on single-neuron spikes as well as on potentials generated by a cortical network comprising more than 3000 model neurons. The generated MEA potentials are substantially affected by both the saline bath covering the brain slice and a (putative) inadvertent saline layer at the interface between the MEA chip and the brain slice. We further explore methods for estimation of current-source density (CSD) from MEA potentials, and find the results to be much less sensitive to the experimental set-up.

GABAergic inhibition shapes interictal dynamics in awake epileptic mice

Epilepsy is characterized by recurrent seizures and brief, synchronous bursts called interictal spikes that are present in-between seizures and observed as transient events in EEG signals. While GABAergic transmission is known to play an important role in shaping healthy brain activity, the role of inhibition in these pathological epileptic dynamics remains unclear. Examining the microcircuits that participate in interictal spikes is thus an important first step towards addressing this issue, as the function of these transient synchronizations in either promoting or prohibiting seizures is currently under debate. To identify the microcircuits recruited in spontaneous interictal spikes in the absence of any proconvulsive drug or anaesthetic agent, we combine a chronic model of epilepsy with in vivo two-photon calcium imaging and multiunit extracellular recordings to map cellular recruitment within large populations of CA1 neurons in mice free to run on a self-paced treadmill. We show that GABAergic neurons, as opposed to their glutamatergic counterparts, are preferentially recruited during spontaneous interictal activity in the CA1 region of the epileptic mouse hippocampus. Although the specific cellular dynamics of interictal spikes are found to be highly variable, they are consistently associated with the activation of GABAergic neurons, resulting in a perisomatic inhibitory restraint that reduces neuronal spiking in the principal cell layer. Given the role of GABAergic neurons in shaping brain activity during normal cognitive function, their aberrant unbalanced recruitment during these transient events could have important downstream effects with clinical implications.

Independent components of neural activity carry information on individual populations

Local field potential (LFP), the low-frequency part of the potential recorded extracellularly in the brain, reflects neural activity at the population level. The interpretation of LFP is complicated because it can mix activity from remote cells, on the order of millimeters from the electrode. To understand better the relation between the recordings and the local activity of cells we used a large-scale network thalamocortical model to compute simultaneous LFP, transmembrane currents, and spiking activity. We used this model to study the information contained in independent components obtained from the reconstructed Current Source Density (CSD), which smooths transmembrane currents, decomposed further with Independent Component Analysis (ICA). We found that the three most robust components matched well the activity of two dominating cell populations: superior pyramidal cells in layer 2/3 (rhythmic spiking) and tufted pyramids from layer 5 (intrinsically bursting). The pyramidal population from layer 2/3 could not be well described as a product of spatial profile and temporal activation, but by a sum of two such products which we recovered in two of the ICA components in our analysis, which correspond to the two first principal components of PCA decomposition of layer 2/3 population activity. At low noise one more cell population could be discerned but it is unlikely that it could be recovered in experiment given typical noise ranges.

Modelling and analysis of local field potentials for studying the function of cortical circuits

The past decade has witnessed a renewed interest in cortical local field potentials (LFPs)--that is, extracellularly recorded potentials with frequencies of up to ~500 Hz. This is due to both the advent of multielectrodes, which has enabled recording of LFPs at tens to hundreds of sites simultaneously, and the insight that LFPs offer a unique window into key integrative synaptic processes in cortical populations. However, owing to its numerous potential neural sources, the LFP is more difficult to interpret than are spikes. Careful mathematical modelling and analysis are needed to take full advantage of the opportunities that this signal offers in understanding signal processing in cortical circuits and, ultimately, the neural basis of perception and cognition.

Frequency dependence of signal power and spatial reach of the local field potential

Despite its century-old use, the interpretation of local field potentials (LFPs), the low-frequency part of electrical signals recorded in the brain, is still debated. In cortex the LFP appears to mainly stem from transmembrane neuronal currents following synaptic input, and obvious questions regarding the 'locality' of the LFP are: What is the size of the signal-generating region, i.e., the spatial reach, around a recording contact? How far does the LFP signal extend outside a synaptically activated neuronal population? And how do the answers depend on the temporal frequency of the LFP signal? Experimental inquiries have given conflicting results, and we here pursue a modeling approach based on a well-established biophysical forward-modeling scheme incorporating detailed reconstructed neuronal morphologies in precise calculations of population LFPs including thousands of neurons. The two key factors determining the frequency dependence of LFP are the spatial decay of the single-neuron LFP contribution and the conversion of synaptic input correlations into correlations between single-neuron LFP contributions. Both factors are seen to give low-pass filtering of the LFP signal power. For uncorrelated input only the first factor is relevant, and here a modest reduction (<50%) in the spatial reach is observed for higher frequencies (>100 Hz) compared to the near-DC ([Formula: see text]) value of about [Formula: see text]. Much larger frequency-dependent effects are seen when populations of pyramidal neurons receive correlated and spatially asymmetric inputs: the low-frequency ([Formula: see text]) LFP power can here be an order of magnitude or more larger than at 60 Hz. Moreover, the low-frequency LFP components have larger spatial reach and extend further outside the active population than high-frequency components. Further, the spatial LFP profiles for such populations typically span the full vertical extent of the dendrites of neurons in the population. Our numerical findings are backed up by an intuitive simplified model for the generation of population LFP.