MEG source imaging method using fast L1 minimum-norm and its applications to signals with brain noise and human resting-state source amplitude images

The present study developed a fast MEG source imaging technique based on Fast Vector-based Spatio-Temporal Analysis using a L1-minimum-norm (Fast-VESTAL) and then used the method to obtain the source amplitude images of resting-state magnetoencephalography (MEG) signals for different frequency bands. The Fast-VESTAL technique consists of two steps. First, L1-minimum-norm MEG source images were obtained for the dominant spatial modes of sensor-waveform covariance matrix. Next, accurate source time-courses with millisecond temporal resolution were obtained using an inverse operator constructed from the spatial source images of Step 1. Using simulations, Fast-VESTAL's performance was assessed for its 1) ability to localize multiple correlated sources; 2) ability to faithfully recover source time-courses; 3) robustness to different SNR conditions including SNR with negative dB levels; 4) capability to handle correlated brain noise; and 5) statistical maps of MEG source images. An objective pre-whitening method was also developed and integrated with Fast-VESTAL to remove correlated brain noise. Fast-VESTAL's performance was then examined in the analysis of human median-nerve MEG responses. The results demonstrated that this method easily distinguished sources in the entire somatosensory network. Next, Fast-VESTAL was applied to obtain the first whole-head MEG source-amplitude images from resting-state signals in 41 healthy control subjects, for all standard frequency bands. Comparisons between resting-state MEG sources images and known neurophysiology were provided. Additionally, in simulations and cases with MEG human responses, the results obtained from using conventional beamformer technique were compared with those from Fast-VESTAL, which highlighted the beamformer's problems of signal leaking and distorted source time-courses.

[1]  R. Hari,et al.  Functional Organization of the Human First and Second Somatosensory Cortices: a Neuromagnetic Study , 1993, The European journal of neuroscience.

[2]  Gareth R. Barnes,et al.  The use of anatomical constraints with MEG beamformers , 2003, NeuroImage.

[3]  M. Tanaka,et al.  Periodic appearance of theta rhythm in the frontal midline area during performance of a mental task. , 1980, Electroencephalography and clinical neurophysiology.

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

[5]  R D Pascual-Marqui,et al.  Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details. , 2002, Methods and findings in experimental and clinical pharmacology.

[6]  R M Leahy,et al.  EEG source localization and imaging using multiple signal classification approaches. , 1999, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[7]  Alec Marantz,et al.  MEG spatio-temporal analysis using a covariance matrix calculated from nonaveraged multiple-epoch data , 1999 .

[8]  G. Barnes,et al.  Statistical flattening of MEG beamformer images , 2003, Human brain mapping.

[9]  M. Weisend,et al.  Magnetoencephalographic characterization of sleep spindles in humans. , 1997, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[10]  Horst A. Eiselt,et al.  Continuous Optimization Models , 1987 .

[11]  K. Matsuura,et al.  A robust reconstruction of sparse biomagnetic sources , 1997, IEEE Transactions on Biomedical Engineering.

[12]  R. Hari,et al.  Phase locking between human primary and secondary somatosensory cortices , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[13]  J. Mäkelä,et al.  Magnetoencephalographic cortical rhythms. , 1997, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[14]  K. Sekihara,et al.  Noise covariance incorporated MEG-MUSIC algorithm: a method for multiple-dipole estimation tolerant of the influence of background brain activity , 1997, IEEE Transactions on Biomedical Engineering.

[15]  Y. Mizuki,et al.  Appearance of Frontal Midline Theta Rhythm and Personality Traits , 1984, Folia psychiatrica et neurologica japonica.

[16]  Ernst Fernando Lopes Da Silva Niedermeyer,et al.  Electroencephalography, basic principles, clinical applications, and related fields , 1982 .

[17]  D. Cohen Detection of magnetic fields outside the human head produced by alpha rhythm currents. , 1970, Electroencephalography and Clinical Neurophysiology.

[18]  F. D. Silva,et al.  Source localization of MEG sleep spindles and the relation to sources of alpha band rhythms , 2002, Clinical Neurophysiology.

[19]  K Sekihara,et al.  Dynamic activation of distinct cytoarchitectonic areas of the human SI cortex after median nerve stimulation , 2001, Neuroreport.

[20]  Mark W. Woolrich,et al.  Using variance information in magnetoencephalography measures of functional connectivity , 2013, NeuroImage.

[21]  Veikko Jousmäki,et al.  Sensorimotor integration in human primary and secondary somatosensory cortices , 1998, Brain Research.

[22]  A. C. Papanicolaou,et al.  Comparison of single current dipole and Magnetic Field Tomography analyses of the cortical response to auditory stimuli , 2005, Brain Topography.

[23]  W. Drongelen,et al.  Localization of brain electrical activity via linearly constrained minimum variance spatial filtering , 1997, IEEE Transactions on Biomedical Engineering.

[24]  Mingxiong Huang,et al.  A non-invasive method for observing hippocampal function , 2003, Neuroreport.

[25]  B. W. van Dijk,et al.  Feasibility of clinical Magnetoencephalography (MEG) functional mapping in the presence of dental artefacts , 2013, Clinical Neurophysiology.

[26]  J.C. Mosher,et al.  Recursive MUSIC: A framework for EEG and MEG source localization , 1998, IEEE Transactions on Biomedical Engineering.

[27]  E. Somersalo,et al.  Visualization of Magnetoencephalographic Data Using Minimum Current Estimates , 1999, NeuroImage.

[28]  A. Ioannides,et al.  Linear transformations of data space in MEG. , 1999, Physics in medicine and biology.

[29]  Hiroshi Shibasaki,et al.  Second somatosensory area (SII) plays a significant role in selective somatosensory attention. , 2002, Brain research. Cognitive brain research.

[30]  Julia P. Owen,et al.  Performance evaluation of the Champagne source reconstruction algorithm on simulated and real M/EEG data , 2012, NeuroImage.

[31]  Toshiki Yoshimine,et al.  Frequency-dependent spatial distribution of human somatosensory evoked neuromagnetic fields , 2002, Neuroscience Letters.

[32]  Seppo P. Ahlfors,et al.  New Six-Layer Magnetically-Shielded Room for MEG , 2002 .

[33]  Thomas T. Liu,et al.  Accurate reconstruction of temporal correlation for neuronal sources using the enhanced dual-core MEG beamformer , 2011, NeuroImage.

[34]  L. Krubitzer,et al.  Evidence for interhemispheric processing of inputs from the hands in human S2 and PV. , 2001, Journal of neurophysiology.

[35]  Thomas T. Liu,et al.  An automatic MEG low-frequency source imaging approach for detecting injuries in mild and moderate TBI patients with blast and non-blast causes , 2012, NeuroImage.

[36]  Rebecca J. Theilmann,et al.  Dual-Core Beamformer for obtaining highly correlated neuronal networks in MEG , 2011, NeuroImage.

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

[38]  William Gaetz,et al.  PRESURGICAL LOCALIZATION OF PRIMARY MOTOR CORTEX IN PEDIATRIC PATIENTS WITH BRAIN LESIONS BY THE USE OF SPATIALLY FILTERED MAGNETOENCEPHALOGRAPHY , 2009, Neurosurgery.

[39]  D. Cheyne,et al.  Localization of human somatosensory cortex using spatially filtered magnetoencephalography , 2003, Neuroscience Letters.

[40]  Seana Coulson,et al.  Magnetoencephalography reveals early activation of V4 in grapheme-color synesthesia , 2010, NeuroImage.

[41]  Erkki Oja,et al.  Independence: a new criterion for the analysis of the electromagnetic fields in the global brain? , 2000, Neural Networks.

[42]  Jeffrey L. Elman,et al.  A novel integrated MEG and EEG analysis method for dipolar sources , 2007, NeuroImage.

[43]  Gregory A. Miller,et al.  A parietal–frontal network studied by somatosensory oddball MEG responses, and its cross-modal consistency , 2005, NeuroImage.

[44]  A. A. Fife,et al.  Advanced electronics for the CTF MEG system. , 2004, Neurology & clinical neurophysiology : NCN.

[45]  Julia M Stephen,et al.  Investigation of the normal proximal somatomotor system using magnetoencephalography , 2003, Clinical Neurophysiology.

[46]  Y. Mizuki,et al.  Differential responses to mental stress in high and low anxious normal humans assessed by frontal midline theta activity. , 1992, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[47]  D. Cohen Magnetoencephalography: Evidence of Magnetic Fields Produced by Alpha-Rhythm Currents , 1968, Science.

[48]  S Taulu,et al.  MEG recordings of DC fields using the signal space separation method (SSS). , 2004, Neurology & clinical neurophysiology : NCN.

[49]  Darren Price,et al.  Investigating the electrophysiological basis of resting state networks using magnetoencephalography , 2011, Proceedings of the National Academy of Sciences.

[50]  J.C. Mosher,et al.  Multiple dipole modeling and localization from spatio-temporal MEG data , 1992, IEEE Transactions on Biomedical Engineering.

[51]  H. Berger Über das Elektrenkephalogramm des Menschen , 1938, Archiv für Psychiatrie und Nervenkrankheiten.

[52]  S Makeig,et al.  Blind separation of auditory event-related brain responses into independent components. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[53]  Matthew J. Brookes,et al.  Beamformer reconstruction of correlated sources using a modified source model , 2007, NeuroImage.

[54]  N. Takahashi,et al.  Frontal Midline Theta Rhythm in Young Healthy Adults , 1997, Clinical EEG.

[55]  S. Taulu,et al.  Suppression of Interference and Artifacts by the Signal Space Separation Method , 2003, Brain Topography.

[56]  Matthew J. Brookes,et al.  Measuring functional connectivity using MEG: Methodology and comparison with fcMRI , 2011, NeuroImage.

[57]  J. Stephen,et al.  Sources on the anterior and posterior banks of the central sulcus identified from magnetic somatosensory evoked responses using Multi‐Start Spatio‐Temporal localization , 2000, Human brain mapping.

[58]  M. Alexander,et al.  Principles of Neural Science , 1981 .

[59]  S. Schultz Principles of Neural Science, 4th ed. , 2001 .

[60]  Li Cui,et al.  Evaluation of signal space separation via simulation , 2007, Medical & Biological Engineering & Computing.

[61]  Christian Wienbruch,et al.  Abnormal slow wave mapping (ASWAM)—A tool for the investigation of abnormal slow wave activity in the human brain , 2007, Journal of Neuroscience Methods.

[62]  Gregory A. Miller,et al.  Somatosensory System Deficits in Schizophrenia Revealed by MEG during a Median-Nerve Oddball Task , 2009, Brain Topography.

[63]  M. Corbetta,et al.  Large-scale cortical correlation structure of spontaneous oscillatory activity , 2012, Nature Neuroscience.

[64]  Raul Coimbra,et al.  Integrated imaging approach with MEG and DTI to detect mild traumatic brain injury in military and civilian patients. , 2009, Journal of neurotrauma.

[65]  M Huang,et al.  Multi-start downhill simplex method for spatio-temporal source localization in magnetoencephalography. , 1998, Electroencephalography and clinical neurophysiology.

[66]  Gary G. R. Green,et al.  Source stability index: A novel beamforming based localisation metric , 2010, NeuroImage.

[67]  R. Leahy,et al.  EEG and MEG: forward solutions for inverse methods , 1999, IEEE Transactions on Biomedical Engineering.

[68]  Roland R. Lee,et al.  Temporal dynamics of ipsilateral and contralateral motor activity during voluntary finger movement , 2004, Human brain mapping.

[69]  J. Stephen,et al.  Central versus peripheral visual field stimulation results in timing differences in dorsal stream sources as measured with MEG , 2002, Vision Research.

[70]  C. Aine,et al.  Multistart Algorithms for MEG Empirical Data Analysis Reliably Characterize Locations and Time Courses of Multiple Sources , 2000, NeuroImage.

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

[72]  Juan L. P. Soto,et al.  Investigation of cross-frequency phase-amplitude coupling in visuomotor networks using magnetoencephalography , 2012, 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[73]  R. Hari,et al.  Magnetoencephalography in the study of human somatosensory cortical processing. , 1999, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[74]  Erkki Oja,et al.  Independent component approach to the analysis of EEG and MEG recordings , 2000, IEEE Transactions on Biomedical Engineering.

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

[76]  C. Rosenberg,et al.  Electroencephalography: Basic Principles, Clinical Applications, and Related Fields, 3rd Ed. , 1994 .

[77]  Kensuke Sekihara,et al.  Modified beamformers for coherent source region suppression , 2006, IEEE Transactions on Biomedical Engineering.

[78]  Karl J. Friston,et al.  Multiple sparse priors for the M/EEG inverse problem , 2008, NeuroImage.

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

[80]  Kensuke Sekihara,et al.  Performance of Prewhitening Beamforming in MEG Dual Experimental Conditions , 2008, IEEE Transactions on Biomedical Engineering.

[81]  J Gross,et al.  REPRINTS , 1962, The Lancet.