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.

[1]  Gary Lynch,et al.  Interactions between recording technique and AMPA receptor modulators , 2002, Brain Research.

[2]  G. Paxinos,et al.  The Rat Brain in Stereotaxic Coordinates , 1983 .

[3]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[4]  Bruce C. Wheeler,et al.  Two-dimensional current source density analysis of propagation delays for components of epileptiform bursts in rat hippocampal slices , 1989, Brain Research.

[5]  T. J. Teyler,et al.  A method for calculating current source density (CSD) analysis without resorting to recording sites outside the sampling volume , 1988, Journal of Neuroscience Methods.

[6]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[7]  William H. Press,et al.  Numerical recipes in C (2nd ed.): the art of scientific computing , 1992 .

[8]  D. Barth,et al.  Three-dimensional analysis of spontaneous and thalamically evoked gamma oscillations in auditory cortex. , 1998, Journal of neurophysiology.

[9]  C. Nicholson,et al.  Experimental optimization of current source-density technique for anuran cerebellum. , 1975, Journal of neurophysiology.

[10]  Klas H. Pettersen,et al.  Current-source density estimation based on inversion of electrostatic forward solution: Effects of finite extent of neuronal activity and conductivity discontinuities , 2006, Journal of Neuroscience Methods.

[11]  U. Mitzdorf Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena. , 1985, Physiological reviews.

[12]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[13]  G Lynch,et al.  Origins and Distribution of Cholinergically Induced β Rhythms in Hippocampal Slices , 2000, The Journal of Neuroscience.

[14]  J. Fermaglich Electric Fields of the Brain: The Neurophysics of EEG , 1982 .

[15]  C. Nicholson,et al.  Theory of current source-density analysis and determination of conductivity tensor for anuran cerebellum. , 1975, Journal of neurophysiology.

[16]  R. Ilmoniemi,et al.  Magnetoencephalography-theory, instrumentation, and applications to noninvasive studies of the working human brain , 1993 .

[17]  Ewa Kublik,et al.  Cortical contribution to sensory volleys recorded at thalamic nuclei of lemniscal and paralemniscal pathways. , 2003, Acta neurobiologiae experimentalis.