Mapping the signal-to-noise-ratios of cortical sources in magnetoencephalography and electroencephalography

Voices behind the left shoulder: two patients with right-sided temporal lobe epilepsy

Auditory vocal hallucinations are sometimes observed in temporal-lobe epilepsy, but are a frequent sign of psychosis and may rarely be mistaken for the latter. Here we report two patients who suffered from auditory vocal hallucinations, described as unintelligible human voices perceived at their left side during epileptic seizures. MEG revealed interictal epileptic discharges within the anterior partition of the right superior temporal gyrus; signal-to-noise ratio of these discharges was overall poor in EEG. The findings suggest that auditory vocal hallucinations without verbal content can evolve in the right hemisphere and are probably independent of language lateralization. This is in accordance with evidence from functional imaging, whereas most previous reports of seizures with auditory vocal hallucinations were confined to the left hemisphere.

The temporal evolution of electromagnetic markers sensitive to the capacity limits of visual short-term memory

An electroencephalographic (EEG) marker of the limited contents of human visual short-term memory (VSTM) has previously been described. Termed contralateral delay activity, this consists of a sustained, posterior, negative potential that correlates with memory load and is greatest contralateral to the remembered hemifield. The current investigation replicates this finding and uses magnetoencephalography (MEG) to characterize its magnetic counterparts and their neural generators as they evolve throughout the memory delay. A parametric manipulation of memory load, within and beyond capacity limits, allows separation of signals that asymptote with behavioral VSTM performance from additional responses that contribute to a linear increase with set-size. Both EEG and MEG yielded bilateral signals that track the number of objects held in memory, and contralateral signals that are independent of memory load. In MEG, unlike EEG, the contralateral interaction between hemisphere and item load is much weaker, suggesting that bilateral and contralateral markers of memory load reflect distinct sources to which EEG and MEG are differentially sensitive. Nonetheless, source estimation allowed both the bilateral and the weaker contralateral capacity-limited responses to be localized, along with a load-independent contralateral signal. Sources of global and hemisphere-specific signals all localized to the posterior intraparietal sulcus during the early delay. However the bilateral load response peaked earlier and its generators shifted later in the delay. Therefore the hemifield-specific response may be more closely tied to memory maintenance while the global load response may be involved in initial processing of a limited number of attended objects, such as their individuation or consolidation into memory.

Practical constraints on estimation of source extent with MEG beamformers

We aimed to determine practical constraints on the estimation of the spatial extent of neuronal activation using MEG beamformers. Correct estimation of spatial extent is a pre-requisite for accurate models of electrical activity, allows one to estimate current density, and enables non-invasive monitoring of functional recovery following stroke. The output of an MEG beamformer is maximum when the correct source model is used, so that the spatial extent of a source can in principal be determined through evaluation of different source models with the beamformer. Here, we simulated 275-channel MEG data using sources of varying spatial extents that followed the cortical geometry. These data were subsequently used to estimate the spatial extent of generic disc elements without knowledge of the underlying surface, and we compared these results to estimates based on cortical surface geometry (with and without error in surface location). We found that disc-shaped source models are too simplistic, particularly for areas with high curvature. For areas with low curvature spatial extent was underestimated, although on average there was a linear relationship between the true and estimated extent. In contrast, cortical surface models gave accurate predictions of spatial extent. However, adding small errors (>2 mm) to the estimated location of the cortical surface abolished this relationship between true and estimated extent, implying that accurate co-registration is needed with such models. Our results show that models exploiting surface information are necessary in order to model spatial extent and in turn current density, but in order to render such models applicable in practical situations, the accuracy of the cortical surface model itself needs to improve.

Magnetoencephalography adds to the surgical evaluation process

Summarizing the podium discussion at the AES 2009, strengths and limitations of magnetoencephalography (MEG) are discussed with regard to basic methodological and clinical aspects in routine screening and presurgical evaluation of patients with epilepsies. Current literature and example cases are used to illustrate MEG contribution to clinical decision making, specifically whether a patient with pharmacoresistant epilepsy can move forward to epilepsy surgery. The main conclusion is that the largest role of MEG, as presently performed in the clinical environment, is to increase the number of patients who can go on to surgery, while it should not be used to deny surgery to any patient.

Onset timing of cross-sensory activations and multisensory interactions in auditory and visual sensory cortices

Here we report early cross-sensory activations and audiovisual interactions at the visual and auditory cortices using magnetoencephalography (MEG) to obtain accurate timing information. Data from an identical fMRI experiment were employed to support MEG source localization results. Simple auditory and visual stimuli (300-ms noise bursts and checkerboards) were presented to seven healthy humans. MEG source analysis suggested generators in the auditory and visual sensory cortices for both within-modality and cross-sensory activations. fMRI cross-sensory activations were strong in the visual but almost absent in the auditory cortex; this discrepancy with MEG possibly reflects the influence of acoustical scanner noise in fMRI. In the primary auditory cortices (Heschl's gyrus) the onset of activity to auditory stimuli was observed at 23 ms in both hemispheres, and to visual stimuli at 82 ms in the left and at 75 ms in the right hemisphere. In the primary visual cortex (Calcarine fissure) the activations to visual stimuli started at 43 ms and to auditory stimuli at 53 ms. Cross-sensory activations thus started later than sensory-specific activations, by 55 ms in the auditory cortex and by 10 ms in the visual cortex, suggesting that the origins of the cross-sensory activations may be in the primary sensory cortices of the opposite modality, with conduction delays (from one sensory cortex to another) of 30-35 ms. Audiovisual interactions started at 85 ms in the left auditory, 80 ms in the right auditory and 74 ms in the visual cortex, i.e., 3-21 ms after inputs from the two modalities converged.

Comparison of noise-normalized minimum norm estimates for MEG analysis using multiple resolution metrics

Noise-normalization has been shown to partly compensate for the localization bias towards superficial sources in minimum norm estimation. However, it has been argued that in order to make inferences for the case of multiple sources, localization properties alone are insufficient. Instead, multiple measures of resolution should be applied to both point-spread and cross-talk functions (PSFs and CTFs). Here, we demonstrate that noise-normalization affects the shapes of PSFs, but not of CTFs. We evaluated PSFs and CTFs for the MNE, dSPM and sLORETA inverse operators, on the metrics dipole localization error (DLE), spatial dispersion (SD) and overall amplitude (OA). We used 306-channel MEG configurations obtained from 17 subjects in a real experiment, including individual noise covariance matrices and head geometries. We confirmed that for PSFs DLE improved after noise normalization, and is zero for sLORETA. However, SD was generally lower for the unnormalized MNE. OA distributions were similar for all three methods, indicating that all three methods may greatly underestimate some sources relative to others. The reliability of differences between methods across subjects was demonstrated using distributions of standard deviations and p-values from paired t-tests. As predicted, the shapes of CTFs were the same for all methods, reflecting the general resolution limits of the inverse problem. This means that noise-normalization is of no consequence where linear estimation procedures are used as “spatial filters.” While low DLE is advantageous for the localization of a single source, or possibly a few spatially distinct sources, the benefit for the case of complex source distributions is not obvious. We suggest that software packages for source estimation should include comprehensive tools for evaluating the performance of different methods.

Sensitivity of MEG and EEG to source orientation

An important difference between magnetoencephalography (MEG) and electroencephalography (EEG) is that MEG is insensitive to radially oriented sources. We quantified computationally the dependency of MEG and EEG on the source orientation using a forward model with realistic tissue boundaries. Similar to the simpler case of a spherical head model, in which MEG cannot see radial sources at all, for most cortical locations there was a source orientation to which MEG was insensitive. The median value for the ratio of the signal magnitude for the source orientation of the lowest and the highest sensitivity was 0.06 for MEG and 0.63 for EEG. The difference in the sensitivity to the source orientation is expected to contribute to systematic differences in the signal-to-noise ratio between MEG and EEG.

Patterns of spontaneous magnetoencephalographic activity in patients with schizophrenia

Magnetoencephalography noninvasively measures the magnetic fields produced by the brain. Pertinent research articles from 1993 to 2009 that measured spontaneous, whole-head magnetoencephalography activity in patients with schizophrenia were reviewed. Data on localization of oscillatory activity and correlation of these findings with psychotic symptoms are summarized. Although the variety of measures used by different research groups makes a quantitative meta-analysis difficult, it appears that magnetoencephalography activity in patients may exhibit identifiable patterns, defined by topographic organization and frequency band. Specifically, 11 of the 12 studies showed increased theta (4-8 Hz) and delta (1-4 Hz) band oscillations in the temporal lobes of patients; of the 10 studies that examined the relationship between oscillatory activity and symptomatology, 8 found a positive correlation between temporal lobe theta activity and positive schizophrenic symptoms. Abnormally high frontal delta activity was not seen. These findings are analyzed in comparison with the electroencephalogram literature on schizophrenics, and possible confounds (e.g., medication effects) are discussed. In the future, magnetoencephalography might be used to assist in diagnosis or might be fruitfully used in conjunction with new neuroscience research approaches such as computational modeling, which may be able to link oscillatory activity and cellular-level pathology.

Cancellation of EEG and MEG signals generated by extended and distributed sources

Extracranial patterns of scalp potentials and magnetic fields, as measured with electro- and magnetoencephalography (EEG, MEG), are spatially widespread even when the underlying source in the brain is focal. Therefore, loss in signal magnitude due to cancellation is expected when multiple brain regions are simultaneously active. We characterized these cancellation effects in EEG and MEG using a forward model with sources constrained on an anatomically accurate reconstruction of the cortical surface. Prominent cancellation was found for both EEG and MEG in the case of multiple randomly distributed source dipoles, even when the number of simultaneous dipoles was small. Substantial cancellation occurred also for locally extended patches of simulated activity, when the patches extended to opposite walls of sulci and gyri. For large patches, a difference between EEG and MEG cancellation was seen, presumably due to selective cancellation of tangentially vs. radially oriented sources. Cancellation effects can be of importance when electrophysiological data are related to hemodynamic measures. Furthermore, the selective cancellation may be used to explain some observed differences between EEG and MEG in terms of focal vs. widespread cortical activity.

Regional differences in the sensitivity of MEG for interictal spikes in epilepsy

MEG interictal spikes as recorded in epilepsy patients are a reflection of intracranial interictal activity. This study investigates the relationship between the estimated sources of MEG spikes and the location, distribution and size of interictal spikes in the invasive ECoG of a group of 38 epilepsy patients that are monitored for pre-surgical evaluation. An amplitude/surface area measure is defined to quantify and rank ECoG spikes. It is found that all MEG spikes are associated with an ECoG spike that is among the three highest ranked in a patient. Among the different brain regions considered, the fronto-orbital, inter-hemispheric, tempero-lateral and central regions stand out. In an accompanying simulation study it is shown that for hypothesized extended sources of larger sizes, as suggested by the data, source location, orientation and curvature can partly explain the observed sensitivity of MEG for interictal spikes.

Propagation of epileptic spikes reconstructed from spatiotemporal magnetoencephalographic and electroencephalographic source analysis

The purpose of this study is to assess the accuracy of spatiotemporal source analysis of magnetoencephalography (MEG) and scalp electroencephalography (EEG) for representing the propagation of frontotemporal spikes in patients with partial epilepsy. This study focuses on frontotemporal spikes, which are typically characterized by a preceding anterior temporal peak followed by an ipsilateral inferior frontal peak. Ten patients with frontotemporal spikes on MEG/EEG were studied. We analyzed the propagation of temporal to frontal epileptic spikes on both MEG and EEG independently by using a cortically constrained minimum norm estimate (MNE). Spatiotemporal source distribution of each spike was obtained on the cortical surface derived from the patient's MRI. All patients underwent an extraoperative intracranial EEG (IEEG) recording covering temporal and frontal lobes after presurgical evaluation. We extracted source waveforms of MEG and EEG from the source distribution of interictal spikes at the sites corresponding to the location of intracranial electrodes. The time differences of the ipsilateral temporal and frontal peaks as obtained by MEG, EEG and IEEG were statistically compared in each patient. In all patients, MEG and IEEG showed similar time differences between temporal and frontal peaks. The time differences of EEG spikes were significantly smaller than those of IEEG in nine of ten patients. Spatiotemporal analysis of MEG spikes models the time course of frontotemporal spikes as observed on IEEG more adequately than EEG in our patients. Spatiotemporal source analysis may be useful for planning epilepsy surgery, by predicting the pattern of IEEG spikes.

Combining fMRI with EEG and MEG in order to relate patterns of brain activity to cognition

The common factor that underlies several types of functional brain imaging is the electric current of masses of dendrites. The prodigious demands for the energy that is required to drive the dendritic currents are met by hemodynamic and metabolic responses that are visualized with fMRI and PET techniques. The high current densities in parallel dendritic shafts and the broad distributions of the loop currents outside the dendrites generate both the scalp EEG and the magnetic fields seen in the MEG. The measurements of image intensities and potential fields provide state variables for modeling. The relationships between the intensities of current density and the electric, magnetic, and hemodynamic state variables are complex and far from proportionate. The state variables are complementary, because the information they convey comes from differing albeit overlapping neural populations, so that efforts to cross-validate localization of neural activity relating to specified cognitive behaviors have not always been successful. We propose an alternative way to use the three methods in combination through studies of hemisphere-wide, high-resolution spatiotemporal patterns of neural activity recorded non-invasively and analyzed with multivariate statistics. Success in this proposed endeavor requires specification of what patterns to look for. At the present level of understanding, an appropriate pattern is any significant departure from random noise in the spectral, temporal and spatial domains that can be scaled into the coarse-graining of time by fMRI/BOLD and the coarse-graining of space by EEG and MEG. Here the requisite patterns are predicted to be large-scale spatial amplitude modulation (AM) of synchronized neuronal signals in the beta and gamma ranges that are coordinated but not correlated with fMRI intensities.

Clinical applications of magnetoencephalography

Magnetoencephalography (MEG), in which magnetic fields generated by brain activity are recorded outside of the head, is now in routine clinical practice throughout the world. MEG has become a recognized and vital part of the presurgical evaluation of patients with epilepsy and patients with brain tumors. We review investigations that show an improvement in the postsurgical outcomes of patients with epilepsy by localizing epileptic discharges. We also describe the most common clinical MEG applications that affect the management of patients, and discuss some applications that are close to having a clinical impact on patients.