Dynamic linear model analysis of optical imaging data acquired from the human neocortex

The amount of light absorbed and scattered by neocortical tissue is altered by neuronal activity. Imaging of intrinsic optical signals (ImIOS), a technique for mapping these activity-evoked optical changes with an imaging detector, has the potential to be useful for both clinical and experimental investigations of the human neocortex. However, its usefulness for human studies is currently limited because intraoperatively acquired ImIOS data is noisy. To improve the reliability and usefulness of ImIOS for human studies, it is desirable to find appropriate methods for the removal of noise artifacts and its statistical analysis. Here we develop a Bayesian, dynamic linear modeling approach that appears to address these problems. A dynamic linear model (DLM) was constructed that included cyclic components to model the heartbeat and respiration artifacts, and a local linear component to model the activity-evoked response. The robustness of the model was tested on a set of ImIOS data acquired from the exposed cortices of six human subjects illuminated with either 535nm or 660nm light. The DLM adequately reduced noise artifacts in these data while reliably preserving their activity-evoked optical responses. To demonstrate how these methods might be used for intraoperative neurosurgical mapping, optical data acquired from a single human subject during direct electrical stimulation of the cortex were quantitatively analyzed. This example showed that the DLM can be used to provide quantitative information about human ImIOS data that is not available through qualitative analysis alone.

[1]  G. Striedter Principles of brain evolution. , 2005 .

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

[3]  G. Ojemann,et al.  Optical imaging of epileptiform and functional activity in human cerebral cortex , 1992, Nature.

[4]  Giovanni Petris,et al.  Dynamic Linear Models with R , 2009 .

[5]  Thierry Blu,et al.  Wavelet-based multi-resolution statistics for optical imaging signals: Application to automated detection of odour activated glomeruli in the mouse olfactory bulb , 2007, NeuroImage.

[6]  G.J. Suaning,et al.  Optical Imaging of Electrically Evoked Visual Signals in Cats: II. ICA "Harmonic Filtering" Noise Reduction , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[7]  D. Hochman,et al.  Intrinsic optical changes in neuronal tissue. Basic mechanisms. , 1997, Neurosurgery clinics of North America.

[8]  G. Mchedlishvili,et al.  The Modular Organization of the Pial Arterial System in Phylogeny , 1984, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[9]  Tadashi Nariai,et al.  Intraoperative intrinsic optical imaging of neuronal activity from subdivisions of the human primary somatosensory cortex. , 2002, Cerebral cortex.

[10]  O. Aguilar,et al.  Bayesian Inference on Latent Structure in Time Series , 1998 .

[11]  Kary L. Myers,et al.  State‐space models for optical imaging , 2007, Statistics in medicine.

[12]  A. Grinvald,et al.  Interactions Between Electrical Activity and Cortical Microcirculation Revealed by Imaging Spectroscopy: Implications for Functional Brain Mapping , 1996, Science.

[13]  Ingo Schießl,et al.  Independent components of the haemodynamic response in intrinsic optical imaging , 2008, NeuroImage.

[14]  L. Cohen Changes in neuron structure during action potential propagation and synaptic transmission. , 1973, Physiological reviews.

[15]  Katsushige Sato,et al.  Functional representation of the finger and face in the human somatosensory cortex: intraoperative intrinsic optical imaging , 2005, NeuroImage.

[16]  Peter A. Bandettini,et al.  Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI , 2006, NeuroImage.

[17]  Dimitrios I. Fotiadis,et al.  Bayesian Methods for fMRI Time-Series Analysis Using a Nonstationary Model for the Noise , 2010, IEEE Transactions on Information Technology in Biomedicine.

[18]  R. Frostig,et al.  Optical imaging of neuronal activity. , 1988, Physiological reviews.

[19]  Jörn Diedrichsen,et al.  Detecting and adjusting for artifacts in fMRI time series data , 2005, NeuroImage.

[20]  Hongtao Ma,et al.  Focal Increases in Perfusion and Decreases in Hemoglobin Oxygenation Precede Seizure Onset in Spontaneous Human Epilepsy , 2007, Epilepsia.

[21]  M. Haglund,et al.  Optical Imaging of Epileptiform Activity in Human Neocortex , 2004, Epilepsia.

[22]  Lawrence L. Wald,et al.  Comparison of physiological noise at 1.5 T, 3 T and 7 T and optimization of fMRI acquisition parameters , 2005, NeuroImage.

[23]  A. Toga,et al.  Temporal spatial differences observed by functional MRI and human intraoperative optical imaging , 2001, NeuroImage.

[24]  J. Besag,et al.  Bayesian Computation and Stochastic Systems , 1995 .

[25]  Michael M Haglund,et al.  Imaging of Intrinsic Optical Signals in Primate Cortex during Epileptiform Activity , 2007, Epilepsia.

[26]  Theodore H Schwartz,et al.  Intraoperative optical imaging of human face cortical topography: a case study , 2004, Neuroreport.

[27]  J. Mayhew,et al.  Cerebral Vasomotion: A 0.1-Hz Oscillation in Reflected Light Imaging of Neural Activity , 1996, NeuroImage.

[28]  Jun S. Liu,et al.  Monte Carlo strategies in scientific computing , 2001 .

[29]  M M Haglund,et al.  Optical imaging of bipolar cortical stimulation. , 1993, Journal of neurosurgery.

[30]  E. Kaplan,et al.  A Principal Components-Based Method for the Detection of Neuronal Activity Maps: Application to Optical Imaging , 2000, NeuroImage.

[31]  Michael A. West,et al.  Evaluation and Comparison of EEG Traces: Latent Structure in Nonstationary Time Series , 1999 .

[32]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.