Extracting functional components of neural dynamics with Independent Component Analysis and inverse Current Source Density

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. The contributions from different cortical layers within V1 could however not be robustly separated and identified with LPA. This is likely due to substantial synchrony in population firing rates across layers, which may be reduced with other stimulus protocols in the future. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

Phase locking of hippocampal CA3 neurons to distal CA1 theta oscillations selectively predicts memory performance

How the coordination of neuronal spiking and brain rhythms between hippocampal subregions supports memory function remains elusive. We studied the interregional coordination of CA3 neuronal spiking with CA1 theta oscillations by recording electrophysiological signals along the proximodistal axis of the hippocampus in rats that were performing a high-memory-demand recognition memory task adapted from humans. We found that CA3 population spiking occurs preferentially at the peak of distal CA1 theta oscillations when memory was tested but only when previously encountered stimuli were presented. In addition, decoding analyses revealed that only population cell firing of proximal CA3 together with that of distal CA1 can predict performance at test in the present non-spatial task. Overall, our work demonstrates an important role for the synchronization of CA3 neuronal activity with CA1 theta oscillations during memory testing.

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.

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

Resolving the mesoscopic missing link: Biophysical modeling of EEG from cortical columns in primates

Event-related potentials (ERP) are among the most widely measured indices for studying human cognition. While their timing and magnitude provide valuable insights, their usefulness is limited by our understanding of their neural generators at the circuit level. Inverse source localization offers insights into such generators, but their solutions are not unique. To address this problem, scientists have assumed the source space generating such signals comprises a set of discrete equivalent current dipoles, representing the activity of small cortical regions. Based on this notion, theoretical studies have employed forward modeling of scalp potentials to understand how changes in circuit-level dynamics translate into macroscopic ERPs. However, experimental validation is lacking because it requires in vivo measurements of intracranial brain sources. Laminar local field potentials (LFP) offer a mechanism for estimating intracranial current sources. Yet, a theoretical link between LFPs and intracranial brain sources is missing. Here, we present a forward modeling approach for estimating mesoscopic intracranial brain sources from LFPs and predict their contribution to macroscopic ERPs. We evaluate the accuracy of this LFP-based representation of brain sources utilizing synthetic laminar neurophysiological measurements and then demonstrate the power of the approach in vivo to clarify the source of a representative cognitive ERP component. To that end, LFP was measured across the cortical layers of visual area V4 in macaque monkeys performing an attention demanding task. We show that area V4 generates dipoles through layer-specific transsynaptic currents that biophysically recapitulate the ERP component through the detailed forward modeling. The constraints imposed on EEG production by this method also revealed an important dissociation between computational and biophysical contributors. As such, this approach represents an important bridge between laminar microcircuitry, through the mesoscopic activity of cortical columns to the patterns of EEG we measure at the scalp.

Modular Data Acquisition System for Recording Activity and Electrical Stimulation of Brain Tissue Using Dedicated Electronics

In this paper, we present a modular Data Acquisition (DAQ) system for simultaneous electrical stimulation and recording of brain activity. The DAQ system is designed to work with custom-designed Application Specific Integrated Circuit (ASIC) called Neurostim-3 and a variety of commercially available Multi-Electrode Arrays (MEAs). The system can control simultaneously up to 512 independent bidirectional i.e., input-output channels. We present in-depth insight into both hardware and software architectures and discuss relationships between cooperating parts of that system. The particular focus of this study was the exploration of efficient software design so that it could perform all its tasks in real-time using a standard Personal Computer (PC) without the need for data precomputation even for the most demanding experiment scenarios. Not only do we show bare performance metrics, but we also used this software to characterise signal processing capabilities of Neurostim-3 (e.g., gain linearity, transmission band) so that to obtain information on how well it can handle neural signals in real-world applications. The results indicate that each Neurostim-3 channel exhibits signal gain linearity in a wide range of input signal amplitudes. Moreover, their high-pass cut-off frequency gets close to 0.6Hz making it suitable for recording both Local Field Potential (LFP) and spiking brain activity signals. Additionally, the current stimulation circuitry was checked in terms of the ability to reproduce complex patterns. Finally, we present data acquired using our system from the experiments on a living rat’s brain, which proved we obtained physiological data from non-stimulated and stimulated tissue. The presented results lead us to conclude that our hardware and software can work efficiently and effectively in tandem giving valuable insights into how information is being processed by the brain.

Local and Volume-Conducted Contributions to Cortical Field Potentials

Brain field potentials (FPs) can reach far from their sources, making difficult to know which waves come from where. We show that modern algorithms efficiently segregate the local and remote contributions to cortical FPs by recovering the generator-specific spatial voltage profiles. We investigated experimentally and numerically the local and remote origin of FPs in different cortical areas in anesthetized rats. All cortices examined show significant state, layer, and region dependent contribution of remote activity, while the voltage profiles help identify their subcortical or remote cortical origin. Co-activation of different cortical modules can be discriminated by the distinctive spatial features of the corresponding profiles. All frequency bands contain remote activity, thus influencing the FP time course, in cases drastically. The reach of different FP patterns is boosted by spatial coherence and curved geometry of the sources. For instance, slow cortical oscillations reached the entire brain, while hippocampal theta reached only some portions of the cortex. In anterior cortices, most alpha oscillations have a remote origin, while in the visual cortex the remote theta and gamma even surpass the local contribution. The quantitative approach to local and distant FP contributions helps to refine functional connectivity among cortical regions, and their relation to behavior.

Different theta frameworks coexist in the rat hippocampus and are coordinated during memory-guided and novelty tasks

Hippocampal firing is organized in theta sequences controlled by internal memory processes and by external sensory cues, but how these computations are coordinated is not fully understood. Although theta activity is commonly studied as a unique coherent oscillation, it is the result of complex interactions between different rhythm generators. Here, by separating hippocampal theta activity in three different current generators, we found epochs with variable theta frequency and phase coupling, suggesting flexible interactions between theta generators. We found that epochs of highly synchronized theta rhythmicity preferentially occurred during behavioral tasks requiring coordination between internal memory representations and incoming sensory information. In addition, we found that gamma oscillations were associated with specific theta generators and the strength of theta-gamma coupling predicted the synchronization between theta generators. We propose a mechanism for segregating or integrating hippocampal computations based on the flexible coordination of different theta frameworks to accommodate the cognitive needs.

Closed-loop control of gamma oscillations in the amygdala demonstrates their role in spatial memory consolidation

Gamma is a ubiquitous brain rhythm hypothesized to support cognitive, perceptual, and mnemonic functions by coordinating neuronal interactions. While much correlational evidence supports this hypothesis, direct experimental tests have been lacking. Since gamma occurs as brief bursts of varying frequencies and durations, most existing approaches to manipulate gamma are either too slow, delivered irrespective of the rhythm's presence, not spectrally specific, or unsuitable for bidirectional modulation. Here, we overcome these limitations with an approach that accurately detects and modulates endogenous gamma oscillations, using closed-loop signal processing and optogenetic stimulation. We first show that the rat basolateral amygdala (BLA) exhibits prominent gamma oscillations during the consolidation of contextual memories. We then boost or diminish gamma during consolidation, in turn enhancing or impairing subsequent memory strength. Overall, our study establishes the role of gamma oscillations in memory consolidation and introduces a versatile method for studying fast network rhythms in vivo.

Theta Oscillations Alternate With High Amplitude Neocortical Population Within Synchronized States

Synchronized states are marked by large-amplitude low-frequency oscillations in the cortex. These states can be seen during quiet waking or slow-wave sleep. Within synchronized states, previous studies have noted a plethora of different types of activity, including delta oscillations (0.5-4 Hz) and slow oscillations (<1 Hz) in the neocortex and large- and small- irregular activity in the hippocampus. However, it is not still fully characterized how neural populations contribute to the synchronized state. Here we apply independent component analysis to parse which populations are involved in different kinds of neocortical activity, and find two populations that alternate throughout synchronized states. One population broadly affects neocortical deep layers, and is associated with larger amplitude slower neocortical oscillations. The other population exhibits theta-frequency oscillations that are not easily observed in raw field potential recordings. These theta oscillations apparently come from below the neocortex, suggesting hippocampal origin, and are associated with smaller amplitude faster neocortical oscillations. Relative involvement of these two alternating populations may indicate different modes of operation within synchronized states.

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.

Hybrid Scheme for Modeling Local Field Potentials from Point-Neuron Networks

With rapidly advancing multi-electrode recording technology, the local field potential (LFP) has again become a popular measure of neuronal activity in both research and clinical applications. Proper understanding of the LFP requires detailed mathematical modeling incorporating the anatomical and electrophysiological features of neurons near the recording electrode, as well as synaptic inputs from the entire network. Here we propose a hybrid modeling scheme combining efficient point-neuron network models with biophysical principles underlying LFP generation by real neurons. The LFP predictions rely on populations of network-equivalent multicompartment neuron models with layer-specific synaptic connectivity, can be used with an arbitrary number of point-neuron network populations, and allows for a full separation of simulated network dynamics and LFPs. We apply the scheme to a full-scale cortical network model for a ∼1 mm² patch of primary visual cortex, predict laminar LFPs for different network states, assess the relative LFP contribution from different laminar populations, and investigate effects of input correlations and neuron density on the LFP. The generic nature of the hybrid scheme and its public implementation in hybridLFPy form the basis for LFP predictions from other and larger point-neuron network models, as well as extensions of the current application with additional biological detail.

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.

New uses of LFPs: Pathway-specific threads obtained through spatial discrimination

Local field potentials (LFPs) reflect the coordinated firing of functional neural assemblies during information coding and transfer across neural networks. As such, it was proposed that the extraordinary variety of cytoarchitectonic elements in the brain is responsible for the wide range of amplitudes and for the coverage of field potentials, which in most cases receive contributions from multiple pathways and populations. The influence of spatial factors overrides the bold interpretations of customary measurements, such as the amplitude and polarity, to the point that their cellular interpretation is one of the hardest tasks in Neurophysiology. Temporal patterns and frequency bands are not exclusive to pathways but rather, the spatial configuration of the voltage gradients created by each pathway is highly specific and may be used advantageously. Recent technical and analytical advances now make it possible to separate and then reconstruct activity for specific pathways. In this review, we discuss how spatial features specific to cells and populations define the amplitude and extension of LFPs, why they become virtually indecipherable when several pathways are co-activated, and then we present the recent advances regarding their disentanglement using spatial discrimination techniques. The pathway-specific threads of LFPs have a simple cellular interpretation, and the temporal fluctuations obtained can be applied to a variety of new experimental objectives and improve existing approaches. Among others, they facilitate the parallel readout of activity in several populations over multiple time scales correlating them with behavior. Also, they access information contained in irregular fluctuations, facilitating the testing of ongoing plasticity. In addition, they open the way to unravel the synaptic nature of rhythmic oscillations, as well as the dynamic relationships between multiple oscillatory activities. The challenge of understanding which waves belong to which populations, and the pathways that provoke them, may soon be overcome.

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.

Ocular and cardiac artifact rejection for real-time analysis in MEG

BACKGROUND

Recently, magnetoencephalography (MEG) based real-time brain computing interfaces (BCI) have been developed to enable novel and promising methods for neuroscience research. It is well known that artifact rejection prior to source localization largely enhances the localization accuracy. However, many BCI approaches neglect real-time artifact removal due to its time consuming process.

NEW METHOD

The method (referred to as ocular and cardiac artifact rejection for real-time analysis, OCARTA) is based on constrained independent component analysis (cICA), where a priori information of the underlying source signals is used to optimize and accelerate signal decomposition. Thereby, prior information is incorporated by using the subject's individual cardiac and ocular activity. The algorithm automatically uses different separation strategies depending on the underlying source activity.

RESULTS

OCARTA was tested and applied to data from three different but most commonly used MEG systems (4D-Neuroimaging, VSM MedTech Inc. and Elekta Neuromag). Ocular and cardiac artifacts were effectively reduced within one iteration at a time delay of 1ms performed on a standard PC (Intel Core i5-2410M).

COMPARISON WITH EXISTING METHODS

The artifact rejection results achieved with OCARTA are in line with the results reported for offline ICA-based artifact rejection methods.

CONCLUSION

Due to the fast and subject-specific signal decomposition the new approach introduced here is capable of real-time ocular and cardiac artifact rejection.

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.

Determining the true polarity and amplitude of synaptic currents underlying gamma oscillations of local field potentials

Fluctuations in successive waves of oscillatory local field potentials (LFPs) reflect the ongoing processing of neuron populations. However, their amplitude, polarity and synaptic origin are uncertain due to the blending of electric fields produced by multiple converging inputs, and the lack of a baseline in standard AC-coupled recordings. Consequently, the estimation of underlying currents by laminar analysis yields spurious sequences of inward and outward currents. We devised a combined analytical/experimental approach that is suitable to study laminated structures. The approach was essayed on an experimental oscillatory LFP as the Schaffer-CA1 gamma input in anesthetized rats, and it was verified by parallel processing of model LFPs obtained through a realistic CA1 aggregate of compartmental units. This approach requires laminar LFP recordings and the isolation of the oscillatory input from other converging pathways, which was achieved through an independent component analysis. It also allows the spatial and temporal components of pathway-specific LFPs to be separated. While reconstructed Schaffer-specific LFPs still show spurious inward/outward current sequences, these were clearly stratified into distinct subcellular domains. These spatial bands guided the localized delivery of neurotransmitter blockers in experiments. As expected, only Glutamate but not GABA blockers abolished Schaffer LFPs when applied to the active but not passive subcellular domains of pyramidal cells. The known chemical nature of the oscillatory LFP allowed an empirical offset of the temporal component of Schaffer LFPs, such that following reconstruction they yield only sinks or sources at the appropriate sites. In terms of number and polarity, some waves increased and others decreased proportional to the concomitant inputs in native multisynaptic LFPs. Interestingly, the processing also retrieved the initiation time for each wave, which can be used to discriminate afferent from postsynaptic cells in standard spike-phase correlations. The applicability of this approach to other pathways and structures is discussed.

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.

Common atlas format and 3D brain atlas reconstructor: infrastructure for constructing 3D brain atlases

One of the challenges of modern neuroscience is integrating voluminous data of diferent modalities derived from a variety of specimens. This task requires a common spatial framework that can be provided by brain atlases. The first atlases were limited to two-dimentional presentation of structural data. Recently, attempts at creating 3D atlases have been made to offer navigation within non-standard anatomical planes and improve capability of localization of different types of data within the brain volume. The 3D atlases available so far have been created using frameworks which make it difficult for other researchers to replicate the results. To facilitate reproducible research and data sharing in the field we propose an SVG-based Common Atlas Format (CAF) to store 2D atlas delineations or other compatible data and 3D Brain Atlas Reconstructor (3dBAR), software dedicated to automated reconstruction of three-dimensional brain structures from 2D atlas data. The basic functionality is provided by (1) a set of parsers which translate various atlases from a number of formats into the CAF, and (2) a module generating 3D models from CAF datasets. The whole reconstruction process is reproducible and can easily be configured, tracked and reviewed, which facilitates fixing errors. Manual corrections can be made when automatic reconstruction is not sufficient. The software was designed to simplify interoperability with other neuroinformatics tools by using open file formats. The content can easily be exchanged at any stage of data processing. The framework allows for the addition of new public or proprietary content.

On the estimation of population-specific synaptic currents from laminar multielectrode recordings

Multielectrode array recordings of extracellular electrical field potentials along the depth axis of the cerebral cortex are gaining popularity as an approach for investigating the activity of cortical neuronal circuits. The low-frequency band of extracellular potential, i.e., the local field potential (LFP), is assumed to reflect synaptic activity and can be used to extract the laminar current source density (CSD) profile. However, physiological interpretation of the CSD profile is uncertain because it does not disambiguate synaptic inputs from passive return currents and does not identify population-specific contributions to the signal. These limitations prevent interpretation of the CSD in terms of synaptic functional connectivity in the columnar microcircuit. Here we present a novel anatomically informed model for decomposing the LFP signal into population-specific contributions and for estimating the corresponding activated synaptic projections. This involves a linear forward model, which predicts the population-specific laminar LFP in response to synaptic inputs applied at different positions along each population and a linear inverse model, which reconstructs laminar profiles of synaptic inputs from laminar LFP data based on the forward model. Assuming spatially smooth synaptic inputs within individual populations, the model decomposes the columnar LFP into population-specific contributions and estimates the corresponding laminar profiles of synaptic input as a function of time. It should be noted that constant synaptic currents at all positions along a neuronal population cannot be reconstructed, as this does not result in a change in extracellular potential. However, constraining the solution using a priori knowledge of the spatial distribution of synaptic connectivity provides the further advantage of estimating the strength of active synaptic projections from the columnar LFP profile thus fully specifying synaptic inputs.

Kernel current source density method

Local field potentials (LFP), the low-frequency part of extracellular electrical recordings, are a measure of the neural activity reflecting dendritic processing of synaptic inputs to neuronal populations. To localize synaptic dynamics, it is convenient, whenever possible, to estimate the density of transmembrane current sources (CSD) generating the LFP. In this work, we propose a new framework, the kernel current source density method (kCSD), for nonparametric estimation of CSD from LFP recorded from arbitrarily distributed electrodes using kernel methods. We test specific implementations of this framework on model data measured with one-, two-, and three-dimensional multielectrode setups. We compare these methods with the traditional approach through numerical approximation of the Laplacian and with the recently developed inverse current source density methods (iCSD). We show that iCSD is a special case of kCSD. The proposed method opens up new experimental possibilities for CSD analysis from existing or new recordings on arbitrarily distributed electrodes (not necessarily on a grid), which can be obtained in extracellular recordings of single unit activity with multiple electrodes.

Generalized theory for current-source-density analysis in brain tissue

The current-source density (CSD) analysis is a widely used method in brain electrophysiology, but this method rests on a series of assumptions, namely that the surrounding extracellular medium is resistive and uniform, and in some versions of the theory, that the current sources are exclusively made by dipoles. Because of these assumptions, this standard model does not correctly describe the contributions of monopolar sources or of nonresistive aspects of the extracellular medium. We propose here a general framework to model electric fields and potentials resulting from current source densities, without relying on the above assumptions. We develop a mean-field formalism that is a generalization of the standard model and that can directly incorporate nonresistive (nonohmic) properties of the extracellular medium, such as ionic diffusion effects. This formalism recovers the classic results of the standard model such as the CSD analysis, but in addition, we provide expressions to generalize the CSD approach to situations with nonresistive media and arbitrarily complex multipolar configurations of current sources. We found that the power spectrum of the signal contains the signature of the nature of current sources and extracellular medium, which provides a direct way to estimate those properties from experimental data and, in particular, estimate the possible contribution of electric monopoles.

Cross-trial correlation analysis of evoked potentials reveals arousal-related attenuation of thalamo-cortical coupling

We describe a computational method for assessing functional connectivity in sensory neuronal networks. The method, which we term cross-trial correlation, can be applied to signals representing local field potentials (LFPs) evoked by sensory stimulations and utilizes their trial-to-trial variability. A set of single trial samples of a given post-stimulus latency from consecutive evoked potentials (EPs) recorded at a given site is correlated with such sets for all other latencies and recording sites. The results of this computation reveal how neuronal activities at various sites and latencies correspond to activation of other sites at other latencies. The method was used to investigate the functional connectivity of thalamo-cortical network of somatosensory system in behaving rats at two levels of alertness: habituated and aroused. We analyzed potentials evoked by vibrissal deflections recorded simultaneously from the ventrobasal thalamus and barrel cortex. The cross-trial correlation analysis applied to the early post-stimulus period (<25 ms) showed that the magnitude of the population spike recorded in the thalamus at 5 ms post-stimulus correlated with the cortical activation at 6-13 ms post-stimulus. This correlation value was reduced at 6-9 ms, i.e. at early postsynaptic cortical response, with increased level of the animals' arousal. Similarly, the aroused state diminished positive thalamo-cortical correlation for subsequent early EP waves, whereas the efficacy of an indirect cortico-fugal inhibition (over 15 ms) did not change significantly. Thus we were able to characterize the state related changes of functional connections within the thalamo-cortical network of behaving animals.