MEG-SIM: a web portal for testing MEG analysis methods using realistic simulated and empirical data

MEG and EEG measure electrophysiological activity in the brain with exquisite temporal resolution. Because of this unique strength relative to noninvasive hemodynamic-based measures (fMRI, PET), the complementary nature of hemodynamic and electrophysiological techniques is becoming more widely recognized (e.g., Human Connectome Project). However, the available analysis methods for solving the inverse problem for MEG and EEG have not been compared and standardized to the extent that they have for fMRI/PET. A number of factors, including the non-uniqueness of the solution to the inverse problem for MEG/EEG, have led to multiple analysis techniques which have not been tested on consistent datasets, making direct comparisons of techniques challenging (or impossible). Since each of the methods is known to have their own set of strengths and weaknesses, it would be beneficial to quantify them. Toward this end, we are announcing the establishment of a website containing an extensive series of realistic simulated data for testing purposes ( http://cobre.mrn.org/megsim/ ). Here, we present: 1) a brief overview of the basic types of inverse procedures; 2) the rationale and description of the testbed created; and 3) cases emphasizing functional connectivity (e.g., oscillatory activity) suitable for a wide assortment of analyses including independent component analysis (ICA), Granger Causality/Directed transfer function, and single-trial analysis.

Different strategies for auditory word recognition in healthy versus normal aging

To explore the effects of commonly encountered pathology on auditory recognition strategies in elderly participants, magnetoencephalographic (MEG) brain activation patterns and performance were examined in 30 elderly [18 controls and 12 elderly with mild cognitive impairment (MCI) or probable Alzheimer's disease (AD)]. It was predicted that participants with known pathology would reveal different networks of brain activation, compared to healthy elderly, which should correlate with poorer performance. Participants heard a list of words representing common objects, twice. After 20 minutes a list of new and old words was presented and participants judged whether each word was heard earlier. MEG responses were analyzed using a semiautomated source modeling procedure. A cluster analysis using all subjects' MEG sources revealed three dominant patterns of activity which correlated with IQ and task performance. The highest performing group revealed activity in premotor, anterior temporal, and superior parietal lobes with little contribution from prefrontal cortex. Performance and brain activation patterns were also compared for individuals with or without abnormalities such as white matter hyperintensities and/or volume reduction evidenced on their MRIs. Memory performance and activation patterns for individuals with white matter hyperintensities resembled the group of MCI/AD patients. These results emphasize the following: (1) general pathology correlates with cognitive decline and (2) full characterization of the health of elderly participants is important in studies of normal aging since random samples from the elderly population are apt to include individuals with subclinical pathology that can affect cognitive performance.

Modelling the magnetic signature of neuronal tissue

Neuronal communication in the brain involves electrochemical currents, which produce magnetic fields. Stimulus-evoked brain responses lead to changes in these fields and can be studied using magneto- and electro-encephalography (MEG/EEG). In this paper we model the spatiotemporal distribution of the magnetic field of a physiologically idealized but anatomically realistic neuron to assess the possibility of using magnetic resonance imaging (MRI) for directly mapping the neuronal currents in the human brain. Our results show that the magnetic field several centimeters from the centre of the neuron is well approximated by a dipole source, but the field close to the neuron is not, a finding particularly important for understanding the possible contrast mechanism underlying the use of MRI to detect and locate these currents. We discuss the importance of the spatiotemporal characteristics of the magnetic field in cortical tissue for evaluating and optimizing an experiment based on this mechanism and establish an upper bound for the expected MRI signal change due to stimulus-induced cortical response. Our simulations show that the expected change of the signal magnitude is 1.6% and its phase shift is 1 degrees . An unexpected finding of this work is that the cortical orientation with respect to the external magnetic field has little effect on the predicted MRI contrast. This encouraging result shows that magnetic resonance contrast directly based on the neuronal currents present in the cortex is theoretically a feasible imaging technique. MRI contrast generation based on neuronal currents depends on the dendritic architecture and we obtained high-resolution optical images of cortical tissue to discuss the spatial structure of the magnetic field in grey matter.

Multidipole analysis of simulated epileptic spikes with real background activity

This simulated magnetoencephalographic study was designed to determine the variability in source parameters with real subject background activity when applying multidipole spatial-temporal dipole analyses, for which the correct model was compared with undermodeled and overmodeled cases. The simulated sources were created from patches of the cortical surface of each subject's MRI. One- and two-source frontal lobe spikes were generated in two cortical regions seen commonly in frontal lobe epilepsy patients tested at our site (orbital frontal and premotor cortex). In general, the modeling results were adequate for the correct model order and the correct model order plus one. In addition, if the localization error was less than 10 mm from the simulated source, the peak latency of the spike and orientation were very reliable, but the peak amplitude was not. The additional source in the overmodeled condition, on the other hand, was not localized reliably across the different epochs within subjects. The results suggest that consistency of the spike localization and inconsistency of other sources will allow one to determine reliably the appropriate model order in real data, and therefore determine single and multifocal spike generators.

Single vs. paired visual stimulation: superposition of early neuromagnetic responses and retinotopy in extrastriate cortex in humans

Neuromagnetic techniques were used in conjunction with magnetic resonance imaging (MRI) techniques to: (1) localize and characterize cortical sources evoked by visual stimuli presented at different locations in the lower right visual field; (2) examine the superposition of cortical responses by comparing the summation of responses to the presentation of single stimuli with responses to paired stimuli; and (3) examine the spatial resolution of magnetoencephalographic (MEG) techniques by comparing the identified source locations evoked by the presentation of single vs. paired stimuli. Using multi-dipole, non-linear minimization analyses, three sources were localized for each stimulus condition during the initial 80-170 ms poststimulus interval for all subjects. In addition to an occipital source, two extrastriate sources were identified: occipital-parietal and occipital-temporal. Each source evidenced a systematic shift in location associated with changes in stimulus placement parallel to the vertical meridian. To our knowledge, this is the first demonstration of retinotopic organization of extrastriate areas, using non-invasive neuromagnetic techniques. The paired presentation of stimuli reflected superposition of the responses evoked by single stimuli but only for early activity up to 150 ms poststimulus. Undersummation was evident after 150 ms. All sources identified for single stimuli were also identified in the paired-stimulus responses; but at the expense of larger errors for some of the estimated parameters.

Multi-start downhill simplex method for spatio-temporal source localization in magnetoencephalography

A multi-start downhill simplex method is examined as a global minimization technique for fitting multidipole, spatio-temporal magnetoencephalography (MEG) data. This procedure has been performed on both simulated and empirical human visual data, known to exhibit complex field patterns due to multiple sources. Unlike some other non-linear fitting techniques the multi-start downhill simplex method does not require users to provide initial guesses for the dipole parameters, hence the fitting procedure is less time-consuming, more objective, and user-friendly. In addition, this method offers more than one adequate solution thus providing a range of uncertainty for the estimated parameters. The Multi-start downhill simplex method is used to fit the non-linear dipole spatial parameters, while the linear temporal parameters are fit using a separate linear fitting procedure. Singular value decomposition (SVD) is also used in order to improve the procedure for determining the adequate number of modeled dipoles.

Nonlinear analysis of biological systems using short M-sequences and sparse-stimulation techniques

The m-sequence pseudorandom signal has been shown to be a more effective probing signal than traditional Gaussian white noise for studying nonlinear biological systems using cross-correlation techniques. The effectiveness is evidenced by the high signal-to-noise (S/N) ratio and speed of data acquisition. However, the "anomalies" that occur in the estimations of the cross-correlations represent an obstacle that prevents m-sequences from being more widely used for studying nonlinear systems. The sparse-stimulation method for measuring system kernels can help alleviate estimation errors caused by anomalies. In this paper, a "padded sparse-stimulation" method is evaluated, a modification of the "inserted sparse-stimulation" technique introduced by Sutter, along with a short m-sequence as a probing signal. Computer simulations show that both the "padded" and "inserted" methods can effectively eliminate the anomalies in the calculation of the second-order kernel, even when short m-sequences were used (length of 1023 for a binary m-sequence, and 728 for a ternary m-sequence). Preliminary experimental data from neuromagnetic studies of the human visual system are also presented, demonstrating that the system kernels can be measured with high signal-to-noise (S/N) ratios using short m-sequences.

Retinotopic organization of human visual cortex: departures from the classical model

Retinotopic mapping strategies similar to those used for invasive electrophysiological studies to identify multiple visual areas in monkeys have been adapted for noninvasive studies in humans, using magnetic recordings of brain activity in conjunction with anatomical magnetic resonance imaging. The retinotopic organization of the primary visual area (V1) in the left hemisphere of human subjects was examined by presenting a small patterned stimuli near the vertical and horizontal meridians in the lower right visual field. In contrast with the classical model of V1 retinotopy, our results suggest that the representation of the horizontal meridian does not necessarily correspond in a one-to-one manner with the base of the calcarine fissure and that some lower field stimuli can activate regions in the lower bank of the fissure. The results also indicate significant individual variability in the details of how V1 maps around the calcarine fissure.