Independent component analysis of high-resolution imaging data identifies distinct functional domains

In the vertebrate brain external stimuli are often represented in distinct functional domains distributed across the cortical surface. Fast imaging techniques used to measure patterns of population activity record movies with many pixels and many frames, i.e., data sets with high dimensionality. Here we demonstrate that principal component analysis (PCA) followed by spatial independent component analysis (sICA), can be exploited to reduce the dimensionality of data sets recorded in the olfactory bulb and the somatosensory cortex of mice as well as the visual cortex of monkeys, without loosing the stimulus-specific responses. Different neuronal populations are separated based on their stimulus-specific spatiotemporal activation. Both, spatial and temporal response characteristics can be objectively obtained, simultaneously. In the olfactory bulb, groups of glomeruli with different response latencies can be identified. This is shown for recordings of olfactory receptor neuron input measured with a calcium-sensitive axon tracer and for network dynamics measured with the voltage-sensitive dye RH 1838. In the somatosensory cortex, barrels responding to the stimulation of single whiskers can be automatically detected. In the visual cortex orientation columns can be extracted. In all cases artifacts due to movement, heartbeat or respiration were separated from the functional signal by sICA and could be removed from the data set. sICA following PCA is therefore a powerful technique for data compression, unbiased analysis and dissection of imaging data of population activity, collected with high spatial and temporal resolution.

[1]  D H HUBEL,et al.  RECEPTIVE FIELDS AND FUNCTIONAL ARCHITECTURE IN TWO NONSTRIATE VISUAL AREAS (18 AND 19) OF THE CAT. , 1965, Journal of neurophysiology.

[2]  Katsuei Shibuki,et al.  Short‐term plasticity visualized with flavoprotein autofluorescence in the somatosensory cortex of anaesthetized rats , 2004, The European journal of neuroscience.

[3]  B. Sakmann,et al.  Dynamic Receptive Fields of Reconstructed Pyramidal Cells in Layers 3 and 2 of Rat Somatosensory Barrel Cortex , 2003, The Journal of physiology.

[4]  D. Chakrabarti,et al.  A fast fixed - point algorithm for independent component analysis , 1997 .

[5]  D. Hubel,et al.  Receptive fields and functional architecture of monkey striate cortex , 1968, The Journal of physiology.

[6]  S Makeig,et al.  Spatially independent activity patterns in functional MRI data during the stroop color-naming task. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Michael P. Stryker,et al.  New Paradigm for Optical Imaging Temporally Encoded Maps of Intrinsic Signal , 2003, Neuron.

[8]  Richard M. Everson,et al.  Independent Component Analysis: Principles and Practice , 2001 .

[9]  K. Obermayer,et al.  Analysis of Calcium Imaging Signals from the Honeybee Brain by Nonlinear Models , 2001, NeuroImage.

[10]  H. Haken,et al.  PHASE TRANSITIONS IN THE HUMAN BRAIN: SPATIAL MODE DYNAMICS , 1992 .

[11]  James S. Schwaber,et al.  Scattered-Light Imaging in Vivo Tracks Fast and Slow Processes of Neurophysiological Activation , 2001, NeuroImage.

[12]  A. Grinvald,et al.  Imaging Cortical Dynamics at High Spatial and Temporal Resolution with Novel Blue Voltage-Sensitive Dyes , 1999, Neuron.

[13]  K. Karhunen Zur Spektraltheorie stochastischer prozesse , 1946 .

[14]  John W. Fisher,et al.  ICA Using Spacings Estimates of Entropy , 2003, J. Mach. Learn. Res..

[15]  Linda B. Buck,et al.  Information coding in the olfactory system: Evidence for a stereotyped and highly organized epitope map in the olfactory bulb , 1994, Cell.

[16]  Wen-Jie Song,et al.  Separation of signal and noise from in vivo optical recording in Guinea pigs using independent component analysis , 2001, Neuroscience Letters.

[17]  Rainer Goebel,et al.  Cortex-based independent component analysis of fMRI time series. , 2004 .

[18]  S Makeig,et al.  Analysis of fMRI data by blind separation into independent spatial components , 1998, Human brain mapping.

[19]  A. Grinvald,et al.  Long-term voltage-sensitive dye imaging reveals cortical dynamics in behaving monkeys. , 2002, Journal of neurophysiology.

[20]  B. Sakmann,et al.  ‐Dynamic representation of whisker deflection by synaptic potentials in spiny stellate and pyramidal cells in the barrels and septa of layer 4 rat somatosensory cortex , 2002, The Journal of physiology.

[21]  K. Obermayer,et al.  Principal Component Analysis and Blind Separation of Sources for Optical Imaging of Intrinsic Signals , 2000, NeuroImage.

[22]  Håkan Johansson,et al.  Modern Techniques in Neuroscience Research , 1999, Springer Berlin Heidelberg.

[23]  Tzyy-Ping Jung,et al.  Independent Component Analysis of Electroencephalographic Data , 1995, NIPS.

[24]  Michael I. Jordan,et al.  Kernel independent component analysis , 2003 .

[25]  Bert Sakmann,et al.  Sub‐ and suprathreshold receptive field properties of pyramidal neurones in layers 5A and 5B of rat somatosensory barrel cortex , 2004, The Journal of physiology.

[26]  H. Kantz,et al.  Nonlinear time series analysis , 1997 .

[27]  L. Cohen,et al.  Representation of Odorants by Receptor Neuron Input to the Mouse Olfactory Bulb , 2001, Neuron.

[28]  A. Grinvald,et al.  Spatio-Temporal Dynamics of Odor Representations in the Mammalian Olfactory Bulb , 2002, Neuron.

[29]  D. Cheyne,et al.  Dynamic cortical activity in the human brain reveals motor equivalence , 1998, Nature.

[30]  Terrence J. Sejnowski,et al.  Blind source separation of more sources than mixtures using overcomplete representations , 1999, IEEE Signal Processing Letters.

[31]  William H. Press,et al.  Numerical recipes in C , 2002 .

[32]  Fabrizio Esposito,et al.  How does spatial extent of fMRI datasets affect independent component analysis decomposition? , 2006, Human brain mapping.

[33]  Rainer W Friedrich,et al.  Temporal Dynamics and Latency Patterns of Receptor Neuron Input to the Olfactory Bulb , 2006, The Journal of Neuroscience.

[34]  Winfried Denk,et al.  Functional organization of sensory input to the olfactory bulb glomerulus analyzed by two-photon calcium imaging , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[35]  Amiram Grinvald,et al.  Iso-orientation domains in cat visual cortex are arranged in pinwheel-like patterns , 1991, Nature.

[36]  A. Grinvald,et al.  Imaging Spatiotemporal Dynamics of Surround Inhibition in the Barrels Somatosensory Cortex , 2003, The Journal of Neuroscience.

[37]  V. Mountcastle Modality and topographic properties of single neurons of cat's somatic sensory cortex. , 1957, Journal of neurophysiology.

[38]  Klaus Obermayer,et al.  Blind signal separation from optical imaging recordings with extended spatial decorrelation , 2000, IEEE Transactions on Biomedical Engineering.

[39]  Jean-Francois Cardoso,et al.  Super-symmetric decomposition of the fourth-order cumulant tensor. Blind identification of more sources than sensors , 1991, [Proceedings] ICASSP 91: 1991 International Conference on Acoustics, Speech, and Signal Processing.

[40]  Richard Axel,et al.  Topographic organization of sensory projections to the olfactory bulb , 1994, Cell.

[41]  T. Sejnowski,et al.  Independent component analysis at the neural cocktail party , 2001, Trends in Neurosciences.

[42]  P. Mitra,et al.  Analysis of dynamic brain imaging data. , 1998, Biophysical journal.

[43]  T. Woolsey,et al.  The structural organization of layer IV in the somatosensory region (S I) of mouse cerebral cortex , 1970 .

[44]  Allan Kardec Barros,et al.  Extraction of event-related signals from multichannel bioelectrical measurements , 2000, IEEE Trans. Biomed. Eng..

[45]  A. Grinvald,et al.  Optical mapping of electrical activity in rat somatosensory and visual cortex , 1985, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[46]  A. Grinvald,et al.  Columnar Resolution of Blood Volume and Oximetry Functional Maps in the Behaving Monkey Implications for fMRI , 2004, Neuron.

[47]  A Grinvald,et al.  In-vivo Optical Imaging of Cortical Architecture and Dynamics , 1999 .

[48]  Jean-Franois Cardoso High-Order Contrasts for Independent Component Analysis , 1999, Neural Computation.

[49]  D. Kleinfeld,et al.  Distributed representation of vibrissa movement in the upper layers of somatosensory cortex revealed with voltage‐sensitive dyes , 1996, The Journal of comparative neurology.

[50]  A. Grinvald,et al.  Dynamics of Ongoing Activity: Explanation of the Large Variability in Evoked Cortical Responses , 1996, Science.

[51]  Rainer Goebel,et al.  Real-time independent component analysis of fMRI time-series , 2003, NeuroImage.

[52]  Takusige Katura,et al.  Isolation of neural activities from respiratory and heartbeat noises for in vivo optical recording in guinea pigs using independent component analysis , 2003, Neuroscience Letters.

[53]  John S George,et al.  A focusing image probe for assessing neural activity in vivo , 1999, Journal of Neuroscience Methods.