Improving human brain mapping via joint inversion of brain electrodynamics and the BOLD signal

We present several methods to improve the resolution of human brain mapping by combining information obtained from surface electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) of the same participants performing the same task in separate imaging sessions. As an initial step in our methods we used independent component analysis (ICA) to obtain task-related sources for both EEG and fMRI. We then used that information in an integrated cost function that attempts to match both data sources and trades goodness of fit in one regime for another. We compared the performance and drawbacks of each method in localizing sources for a dual visual evoked response experiment, and we contrasted the results of adding fMRI information to simple EEG-only inversion methods. We found that adding fMRI information in a variety of ways gives superior results to classical minimum norm source estimation. Our findings lead us to favor a method which attempts to match EEG scalp dynamics along with voxel power obtained from ICA-processed blood oxygenation level dependent (BOLD) data; this method of joint inversion enables us to treat the two data sources as symmetrically as possible.

[1]  T. Sejnowski,et al.  Removing electroencephalographic artifacts by blind source separation. , 2000, Psychophysiology.

[2]  John Skilling,et al.  Data analysis : a Bayesian tutorial , 1996 .

[3]  R. Quiroga,et al.  Unmixing concurrent EEG-fMRI with parallel independent component analysis. , 2008, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[4]  R. Ilmoniemi,et al.  Interpreting magnetic fields of the brain: minimum norm estimates , 2006, Medical and Biological Engineering and Computing.

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

[6]  Christian Jutten,et al.  Blind separation of sources, part I: An adaptive algorithm based on neuromimetic architecture , 1991, Signal Process..

[7]  Terrence J. Sejnowski,et al.  Complex Independent Component Analysis of Frequency-Domain Electroencephalographic Data , 2003, Neural Networks.

[8]  R. C. Oldfield The assessment and analysis of handedness: the Edinburgh inventory. , 1971, Neuropsychologia.

[9]  Polina Golland,et al.  A distributed spatio-temporal EEG/MEG inverse solver , 2009, NeuroImage.

[10]  P. Berg,et al.  A fast method for forward computation of multiple-shell spherical head models. , 1994, Electroencephalography and clinical neurophysiology.

[11]  T. Sejnowski,et al.  Electroencephalographic Brain Dynamics Following Manually Responded Visual Targets , 2004, PLoS biology.

[12]  T. Mexia,et al.  Author ' s personal copy , 2009 .

[13]  David Poeppel,et al.  Reconstructing spatio-temporal activities of neural sources using an MEG vector beamformer technique , 2001, IEEE Transactions on Biomedical Engineering.

[14]  Bart Vanrumste,et al.  Journal of Neuroengineering and Rehabilitation Open Access Review on Solving the Inverse Problem in Eeg Source Analysis , 2022 .

[15]  O. Sporns,et al.  Complex brain networks: graph theoretical analysis of structural and functional systems , 2009, Nature Reviews Neuroscience.

[16]  S.C. Strother,et al.  Evaluating fMRI preprocessing pipelines , 2006, IEEE Engineering in Medicine and Biology Magazine.

[17]  Vince D. Calhoun,et al.  Neuronal chronometry of target detection: Fusion of hemodynamic and event-related potential data , 2005, NeuroImage.

[18]  Anders M. Dale,et al.  Vector-based spatial–temporal minimum L1-norm solution for MEG , 2006, NeuroImage.

[19]  JuttenChristian,et al.  Blind separation of sources, Part 1 , 1991 .

[20]  T. Sejnowski,et al.  Analysis and visualization of single‐trial event‐related potentials , 2001, Human brain mapping.

[21]  K H Ting,et al.  Automatic correction of artifact from single-trial event-related potentials by blind source separation using second order statistics only. , 2006, Medical engineering & physics.

[22]  A K Liu,et al.  Spatiotemporal imaging of human brain activity using functional MRI constrained magnetoencephalography data: Monte Carlo simulations. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Stephanie Ortigue,et al.  Using ICA and realistic BOLD models to obtain joint EEG/fMRI solutions to the problem of source localization , 2009, NeuroImage.

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

[25]  C. Michel,et al.  Electromagnetic Inverse Solutions in Anatomically Constrained Spherical Head Models , 2004, Brain Topography.

[26]  Stefano Alliney,et al.  An algorithm for the minimization of mixed l1 and l2 norms with application to Bayesian estimation , 1994, IEEE Trans. Signal Process..

[27]  Karl J. Friston,et al.  An empirical Bayesian solution to the source reconstruction problem in EEG , 2005, NeuroImage.

[28]  Peter A. Robinson,et al.  BOLD responses to stimuli: Dependence on frequency, stimulus form, amplitude, and repetition rate , 2006, NeuroImage.

[29]  C. Wunsch The Ocean Circulation Inverse Problem , 1996 .

[30]  J. Gotman,et al.  A system for automatic artifact removal in ictal scalp EEG based on independent component analysis and Bayesian classification , 2006, Clinical Neurophysiology.

[31]  Karl J. Friston,et al.  Anatomically Informed Basis Functions for EEG Source Localization: Combining Functional and Anatomical Constraints , 2002, NeuroImage.

[32]  Lars Kai Hansen,et al.  Optimizing the fMRI data-processing pipeline using prediction and reproducibility performance metrics: I. A preliminary group analysis , 2004, NeuroImage.

[33]  M. Fuchs,et al.  Depth normalization in MEG/EEG current density imaging , 1996, Proceedings of 18th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[34]  D. Lehmann,et al.  Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. , 1994, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[35]  S. Debener,et al.  Mining EEG-fMRI using independent component analysis. , 2009, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[36]  Donald B. Percival,et al.  Spectral Analysis for Physical Applications: Harmonic Analysis , 1993 .

[37]  R Grave de Peralta Menendez,et al.  Imaging the electrical activity of the brain: ELECTRA , 2000, Human brain mapping.

[38]  Olaf Hauk,et al.  Keep it simple: a case for using classical minimum norm estimation in the analysis of EEG and MEG data , 2004, NeuroImage.

[39]  Erkki Oja,et al.  Independent component analysis: algorithms and applications , 2000, Neural Networks.

[40]  J. J. Ermer,et al.  Rapidly recomputable EEG forward models for realistic head shapes. , 2001, Physics in medicine and biology.

[41]  V. Calhoun,et al.  Joint independent component analysis for simultaneous EEG-fMRI: principle and simulation. , 2008, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[42]  Bart Vanrumste,et al.  Review on solving the forward problem in EEG source analysis , 2007, Journal of NeuroEngineering and Rehabilitation.

[43]  Bin He,et al.  fMRI–EEG integrated cortical source imaging by use of time-variant spatial constraints , 2008, NeuroImage.

[44]  A. Walden,et al.  Spectral analysis for physical applications : multitaper and conventional univariate techniques , 1996 .

[45]  Scott T. Grafton,et al.  Spatio-Temporal Dynamics of Human Intention Understanding in Temporo-Parietal Cortex: A Combined EEG/fMRI Repetition Suppression Paradigm , 2009, PloS one.

[46]  John S. Allen,et al.  CODE-2 : moored array and large-scale data report , 1985 .

[47]  James V. Stone Independent Component Analysis: A Tutorial Introduction , 2007 .

[48]  Martin Fodslette Meiller A Scaled Conjugate Gradient Algorithm for Fast Supervised Learning , 1993 .

[49]  Riitta Hari,et al.  Towards natural stimulation in fMRI—Issues of data analysis , 2007, NeuroImage.

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

[51]  R. Buxton,et al.  Dynamics of blood flow and oxygenation changes during brain activation: The balloon model , 1998, Magnetic resonance in medicine.

[52]  E. Halgren,et al.  Dynamic Statistical Parametric Mapping Combining fMRI and MEG for High-Resolution Imaging of Cortical Activity , 2000, Neuron.

[53]  Heather Ting Ma,et al.  Effects of the physiological parameters on the signal-to-noise ratio of single myoelectric channel , 2007, Journal of NeuroEngineering and Rehabilitation.

[54]  W. Freeman,et al.  Combining fMRI with EEG and MEG in order to relate patterns of brain activity to cognition. , 2009, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.