Ewa Kublik

Inverse Current-Source Density Method in 3D: Reconstruction Fidelity, Boundary Effects, and Influence of Distant Sources

Estimation of the continuous current-source density in bulk tissue from a finite set of electrode measurements is a daunting task. Here we present a methodology which allows such a reconstruction by generalizing the one-dimensional inverse CSD method. The idea is to assume a particular plausible form of CSD within a class described by a number of parameters which can be estimated from available data, for example a set of cubic splines in 3D spanned on a fixed grid of the same size as the set of measurements. To avoid specificity of particular choice of reconstruction grid we add random jitter to the points positions and show that it leads to a correct reconstruction. We propose different ways of improving the quality of reconstruction which take into account the sources located outside the recording region through appropriate boundary treatment. The efficiency of the traditional CSD and variants of inverse CSD methods is compared using several fidelity measures on different test data to investigate when one of the methods is superior to the others. The methods are illustrated with reconstructions of CSD from potentials evoked by stimulation of a bunch of whiskers recorded in a slab of the rat forebrain on a grid of 4×5×7 positions.

Gating of the sensory activity within barrel cortex of the awake rat

Abstract In rat barrel cortex, evoked potentials (EPs) to vibrissa stimulation can be divided into two distinct classes according to the relative contribution of their principal components. Our experiments support the notion that these components can be attributed to activation of two pyramidal cell populations: supra- and infragranular. With well-habituated stimuli EPs are dominated by a component related to the supragranular cells (class 1). However, the first reinforcement of vibrissa stimulation in the classical aversive paradigm favours the appearance of EPs dominated by a component characteristic of infragranular cells which matches with activation in the surround zone of the barrel field (class 2). Similar dynamic changes of the relative occurrence of the two EP classes follow other aversive stimuli, including pressing the animal’s ear and restraining a whisker. We hypothesize that neuromodulatory action elicited by contextual stimulation activates all neurons in the principal barrel column, including those providing an output to the surrounding barrels. In the classical conditioning paradigm this mechanism may lead to experience-dependent changes within the intracortical network.

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

Local field potentials have good temporal resolution but are blurred due to the slow spatial decay of the electric field. For simultaneous recordings on regular grids one can reconstruct efficiently the current sources (CSD) using the inverse Current Source Density method (iCSD). It is possible to decompose the resultant spatiotemporal information about the current dynamics into functional components using Independent Component Analysis (ICA). We show on test data modeling recordings of evoked potentials on a grid of 4×5×7 points that meaningful results are obtained with spatial ICA decomposition of reconstructed CSD. The components obtained through decomposition of CSD are better defined and allow easier physiological interpretation than the results of similar analysis of corresponding evoked potentials in the thalamus. We show that spatiotemporal ICA decompositions can perform better for certain types of sources but it does not seem to be the case for the experimental data studied. Having found the appropriate approach to decomposing neural dynamics into functional components we use the technique to study the somatosensory evoked potentials recorded on a grid spanning a large part of the forebrain. We discuss two example components associated with the first waves of activation of the somatosensory thalamus. We show that the proposed method brings up new, more detailed information on the time and spatial location of specific activity conveyed through various parts of the somatosensory thalamus in the rat.

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.

Identification of principal components in cortical evoked potentials by brief surface cooling.

OBJECTIVES The evoked potential recorded by a single electrode in rats barrel cortex after whisker stimulation was shown to be composed of two main principal components shifted in time by about 3 ms. The purpose of this study was to verify the hypothesis that these components represent activity of supra- and infragranular pyramidal cell classes. RESULTS Our results show that a brief cooling pulse applied to the cortical surface abolishes the shorter latency component, which may therefore be attributed to the response of supragranular pyramidal cells. CONCLUSIONS The longer latency principal component, which disappears only with strong cooling pulses, is proposed to represent postsynaptic activity of infragranular pyramidal neurons.

Transient changes of electrical activity in the rat barrel cortex during conditioning

To reveal the dynamics of neurophysiological changes in the rat barrel cortex induced by conditioned stimulation we recorded the local micro-electroencephalographic (EEG) activity and evoked potentials (EPs) in barrel cortex to stimulation of a single vibrissa before and after pairing it with a mild electric shock applied to the rats tail. Following the introduction of the reinforcing stimulus, the amplitude of the first negative component of evoked potentials in the cortex on the conditioned side grew in relation to the same component of control potentials, evoked by stimulation of the opposite symmetrical vibrissa. This change was accompanied by a latent decrease in spectral power of the EEG within the alpha and beta frequency bands in both hemispheres. The observed changes in both of these electrical manifestations of enhanced neuronal activity reverted after two (EP) or three (EEG) days of conditioning. These results are discussed in relation to the putative activity of neuromodulatory systems.

Spontaneous variability reveals principal components in cortical evoked potentials

USING principal component analysis, we studied trial to trial, spontaneous variability of evoked potentials (EPs) recorded from rat barrel cortex after whisker stimulation. This method allowed for extraction of two distinct components of EP which overlapped in the time domain. Our results are consonant with the previously described depth distribution of current sources and the extracted components can be therefore attributed to activities of two pyramidal cell classes: supra- and infragranular. Qualitatively similar results were found in both anaesthetized and alert animals.

Retinal Origin of Electrically Evoked Potentials in Response to Transcorneal Alternating Current Stimulation in the Rat

PURPOSE Little is known about the physiological mechanisms underlying the reported therapeutic effects of transorbital alternating current stimulation (ACS) in vision restoration, or the origin of the recorded electrically evoked potentials (EEPs) during such stimulation. We examined the issue of EEP origin and electrode configuration for transorbital ACS and characterized the physiological responses to CS in different structures of the visual system. METHODS We recorded visually evoked potentials (VEPs) and EEPs from the rat retina, visual thalamus, tectum, and visual cortex. The VEPs were evoked by light flashes and EEPs were evoked by electric stimuli delivered by two electrodes placed either together on the same eye or on the eyeball and in the neck. Electrically evoked potentials and VEPs were recorded before and after bilateral intraorbital injections of tetrodotoxin that blocked retinal ganglion cell activity. RESULTS Tetrodotoxin abolished VEPs at all levels in the visual pathway, confirming successful blockage of ganglion cell activity. Tetrodotoxin also abolished EEPs and this effect was independent of the stimulating electrode configurations. CONCLUSIONS Transorbital electrically evoked responses in the visual pathway, irrespective of reference electrode placement, are initiated by activation of the retina and not by passive conductance and direct activation of neurons in other visual structures. Thus, placement of stimulating electrodes exclusively around the eyeball may be sufficient to achieve therapeutic effects.

The 5-12 Hz oscillations in the barrel cortex of awake rats--sustained attention during behavioral idling?

OBJECTIVE 5-12 Hz oscillations, observed in cortical LFP of awake rats during quiet immobility, were proposed to be either (i) epileptic events or (ii) physiological alpha-like oscillations, manifesting an idling state of the cortex. We aimed to test this controversy. METHODS We recorded LFP from the barrel cortex of awake Wistar rats, while applying weak tactile (whisker) and stronger arousing (electrical) stimuli. RESULTS We observed a mean effect of desynchronization of the 5-12 Hz rhythm by the weak tactile stimulation. Arousal reduced the incidence of the 5-12 Hz oscillations and increased the desynchronizing power of tactile stimuli. CONCLUSIONS Oscillations that can be disrupted by weak, purely tactile stimulation, and whose incidence is reduced by increased arousal, should be interpreted as a physiological phenomenon typical for behavioral idling while the cerebral cortex maintains sensory sensitivity. SIGNIFICANCE Our results contradict the view that the 5-12 Hz oscillatory activity, often observed in fronto-parietal cortical regions of Wistar rats, represents epileptic discharges. Rather, this activity provides a model for studying the physiology of alpha/mu oscillations.

Sorting functional classes of evoked potentials by wavelets.

Evoked potentials (EPs) recorded within the primary sensory cortex of non-anesthetized rats vary considerably with each peripheral stimulation. We have previously shown that most of this variance reflects the shift of cortical activation between habituated and aroused states. Here we show that a method of matching the potential’s course by wavelet functions can reliably differentiate single EPs and may therefore, be used as a probe for indicating the current activation state of the cortex.