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 4x5x7 positions.

Long‐range glutamatergic and GABAergic projections participate in temporal coordination of neuronal activity in distributed cortical areas. In the hippocampus, GABAergic neurons project to the medial septum and retrohippocampal areas. Many GABAergic projection cells express somatostatin (SOM+) and, together with locally terminating SOM+ bistratified and O‐LM cells, contribute to dendritic inhibition of pyramidal cells. We tested the hypothesis that diversity in SOM+ cells reflects temporal specialization during behavior using extracellular single cell recording and juxtacellular neurobiotin‐labeling in freely moving rats. We have demonstrated that rare GABAergic projection neurons discharge rhythmically and are remarkably diverse. During sharp wave‐ripples, most projection cells, including a novel SOM+ GABAergic back‐projecting cell, increased their activity similar to bistratified cells, but unlike O‐LM cells. During movement, most projection cells discharged along the descending slope of theta cycles, but some fired at the trough jointly with bistratified and O‐LM cells. The specialization of hippocampal SOM+ projection neurons complements the action of local interneurons in differentially phasing inputs from the CA3 area to CA1 pyramidal cell dendrites during sleep and wakefulness. Our observations suggest that GABAergic projection cells mediate the behavior‐ and network state‐dependent binding of neuronal assemblies amongst functionally‐related brain regions by transmitting local rhythmic entrainment of neurons in CA1 to neuronal populations in other areas. © 2016 The Authors Hippocampus Published by Wiley Periodicals, Inc.

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.

Behavioural reactions to sensory stimuli vary with the level of arousal, but little is known about the underlying reorganization of neuronal networks. In this study, we use chronic recordings from the somatosensory regions of the thalamus and cortex of behaving rats together with a novel analysis of functional connectivity to show that during low arousal tactile signals are transmitted via the ventral posteromedial thalamic nucleus (VPM), a first-order thalamic relay, to the primary somatosensory (barrel) cortex and then from the cortex to the posterior medial thalamic nucleus (PoM), which plays a role of a higher-order thalamic relay. By contrast, during high arousal this network scheme is modified and both VPM and PoM transmit peripheral input to the barrel cortex acting as first-order relays. We also show that in urethane anaesthesia PoM is largely excluded from the thalamo-cortical loop. We thus demonstrate a way in which the thalamo-cortical system, despite its fixed anatomy, is capable of dynamically reconfiguring the transmission route of a sensory signal in concert with the behavioural state of an animal.