|Home | About | Journals | Submit | Contact Us | Français|
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.
One of the main differences between magnetoencephalography (MEG) and electroencephalography (EEG) is thought to be that MEG is mainly sensitive to tangentially oriented sources, whereas EEG can detect sources of all orientations (Baillet et al. 2001; Cohen and Cuffin 1991; Hämäläinen et al. 1993). Theoretically, a radially oriented current dipole produces no magnetic field outside a spherically symmetric volume conductor (Baule and McFee 1965; Grynszpan and Geselowitz 1973). However, since the human head is not exactly spherically symmetric, the radial orientation is not well defined, and an approximately radial source in the brain is not necessarily silent in MEG (Cuffin 1990). Therefore, it is possible that all source orientations generate a magnetic field that is detectable by MEG.
The relative magnitude of the MEG signal generated by a radial and a tangential source at the same location was examined experimentally by Melcher and Cohen (1988), who found that the signals differed by a factor of about 6 for several locations in the rabbit brain. Similar values were found in simulations taking into account realistic tissue properties of the human head (Haueisen et al. 1995). Only a small number of source locations, however, was examined in the previous studies. In the present simulation study we mapped the dependency of MEG and EEG signal magnitude on the orientation of a current dipole in a dense grid of source locations in the human cerebral cortex. We examined the ratio of the signals due to a dipole having an orientation with the lowest vs. the highest sensitivity for an array of MEG or EEG sensors. A forward solution for each source location was constructed using a boundary-element method. The source orientations with the lowest and the highest sensitivity were determined with the help of the singular value decomposition (SVD). Our goal was to confirm the existence of an insensitive source orientation for MEG at any given location in the human brain.
The MEG and EEG forward models were based on anatomical data from a 42-year old male volunteer subject. The procedures were approved by the Massachusetts General Hospital Institutional Review Board. High-resolution structural T1-weighted magnetic resonance images (MRIs) were acquired on a 1.5T Siemens scanner (TR = 2530 ms, TE = 3.25 ms, flip angle = 7°, 128 sagittal slices, slice thickness = 1.3 mm, voxel size = 1.3 × 1.0 × 1.3 mm3). A surface mesh of the gray-white matter border of the cerebral cortex with about 300 000 vertices was constructed from the MRI using the Free Surfer software (Dale et al. 1999; Fischl et al. 1999). Each vertex represented a location of an MEG and EEG source.
The conductivity geometry of the head was modeled using a boundary element method (BEM) with three compartments: the brain, the skull, and the scalp (Hämäläinen and Sarvas 1989). The boundaries were determined from the MRI by identifying the inner and outer surfaces of the skull and the outer surface of the skin (Liu et al. 2002). Each surface consisted of 5120 vertices. The linear collocation method was used in the BEM calculations (Mosher et al. 1999). The conductivity ratios of 1:0.0125:1 were assumed for brain:skull:scalp (Geddes and Baker 1967).
The N × 3M-dimensional forward matrix A describes the relationship between the MEG and EEG signals in an array of N sensors and the three components of a dipole moment vector at M locations in the brain. Here M was about 300,000. The MEG sensor array consisted of N = 102 MEG magnetometers, measuring the magnetic field component approximately radial to the surface of the head; the EEG sensor array comprised N = 60 scalp electrodes. The column vectors aj of A correspond to the signal patterns generated by a single source element of unit amplitude. For EEG, the average reference was used. The forward matrix was computed using the MNE software (http://www.nmr.mgh.harvard.edu/martinos/userInfo/data/sofMNE.php).
To examine the sensitivity of the MEG and EEG sensor arrays to sources of different orientations, we computed the SVD of the N × 3-dimensional dipolar gain matrix Ak = [a3(k-1)+1,a3(k-1)+2,a3(k-1)+3] at each location indexed by k = 1, …, M (Huang et al. 2007):
where uk,j and vk,j are the left and right singular vectors, respectively, and λk,1, λk,2, and λk,3 are the singular values. The N-dimensional vector gk,j = λk,juk,j describes the signal pattern in the MEG or EEG sensor array generated by a current dipole of unit amplitude in the orientation vk,j. We denote the largest and smallest singular values, corresponding to the dipole orientation to which the MEG sensor array is the most and the least sensitive, as λk,1 = λk,max and λk,3 = λk,min, respectively. In the spherically symmetric volume conductor model λk,min = 0 for MEG for all source locations, because the radial source component does not generate any MEG signals.
was computed for all the locations on the cortical surface. This measure of relative sensitivity to the “suppressed” source orientation is analogous to the suppression ratio for radial vs. tangential dipoles studied previously (Haueisen et al. 1995; Melcher and Cohen 1988). Melcher and Cohen used the maximum signal power among pairs of planar gradiometers, whereas our definition is similar to the Rp measure of Haueisen et al., which is based on the signal power in an array of sensors. Different definitions of the suppression ratio, however, have been reported to provide converging results (Haueisen et al. 1995). We note that instead of defining the “radial” and “tangential” orientations on the basis of a sphere fit on the local curvature of the inner surface of skull near the source location (Haueisen et al. 1995), the SVD automatically provides the source orientation of lowest and highest signal power in the sensor array.
Distributions of the ratio Rλ across the cerebral cortex are shown in Fig..1.1. For MEG, the median value of Rλ was 0.06 and 95% of the values were below 0.23. The prominence of small values for Rλ for MEG suggests that for most locations on the cortex there was a source orientation with which little if any MEG signals was generated. In contrast, for EEG the median value of Rλ was 0.63 and 95% of the values were above 0.42.
For superficial regions of the lateral surface of the temporal, parietal, and occipital cortices, the ratio Rλ for MEG was typically small, indicating the presence of a source orientation to which MEG is insensitive to. The largest values of the ratio Rλ for MEG were found in the orbitofrontal regions and sulcal regions of the lateral frontal cortex, as well as in the Sylvian fissure including the insula, and the medial surface. Occasional isolated vertices with large values of Rλ were seen at the crowns of gyri; these are likely to be caused by numerical inaccuracy of the forward model for sources that are located close to the inner skull boundary. For EEG, the smallest values for Rλ (~0.2) were found in the orbitofrontal and temporal pole regions, likely reflecting a limited coverage of the EEG scalp electrode array in the face and temple areas. Also for MEG, the results for the orbitofrontal region are likely to be affected by the positioning of the sensor array (Marinkovic et al. 2004).
The sensitivity of MEG and EEG for three orthogonal source orientations at each location on the cortex is shown in Fig. 2. The sensitivity was quantified by the three singular values of the dipole gain matrix Ak. The singular value corresponding to the highest sensitivity (λk,max) showed a strong depth dependency for MEG. The spatial distribution of the smallest singular value (λk,min) for MEG was similar to that of the ratio Rλ shown in Fig. 1, except that the large values of Rλ for MEG occurred in the subcortical regions near the center of the head where λk,max was very small.
The histograms in Fig. 2a illustrate how the sensitivity of MEG for one source orientation (corresponding to λk,min) is typically much smaller than the other orientations. For EEG, the histograms indicate that there is no orientation for which the sensitivity is very small (Fig. 2b). Furthermore, for MEG the distributions of the two largest singular values (λk,1 and λk,2) were similar, whereas for EEG the distribution of the two smallest singular values (λk,2 and λk,3) were similar. The two similar distributions in each case are likely to be mainly due to tangentially oriented sources, which generated similar signal patters rotated by 90° with respect to each other, whereas a radial source generates a signal pattern that is qualitatively different from that from the tangential sources.
To examine the dependency of MEG and EEG signals on the source orientation, we compared the forward solutions for current dipoles with the orientation of the lowest and the highest sensitivity across the cerebral cortex, i.e., dipole orientations generating the smallest and the largest signals. For most locations, especially the superficial ones, there was a “quiet” source orientation for MEG, analogously to the radial orientation in the spherical model. For EEG, the sensitivity depended on the source orientation much less than for MEG.
The mean value for the sensitivity ratio Rλ for MEG (0.06 ≈ 1/17) was somewhat smaller than the experimental value (1/6) observed in the rabbit brain (Melcher and Cohen 1988). This difference is not surprising, given the difficulty of aligning an artificial current dipole exactly along the orientation of the lowest sensitivity. The SVD analysis provided the orientation of the lowest sensitivity based on the signal power in the array of sensors, rather than relying on the “radial” orientation defined in terms of the local curvature of the skull (Haueisen et al. 1995; Melcher and Cohen 1988).
In some regions, sources with the orientation of lowest MEG sensitivity generated substantial MEG signals, up to about 20% of that generated by a source with the orientation of highest MEG sensitivity at the same location. A region with a non-vanishing Rλ for MEG implies that the spherically symmetric volume conductor model is not accurate for that region, warranting the use a more realistic head model in inverse modeling also for MEG. This is particularly important for deep cortical and subcortical regions (Parkkonen et al. 2009; Tesche et al. 1996). Further anatomical information about conductivity properties of the head could be taken into account by using finite element modeling (FEM) (Haueisen et al. 1995; Haueisen et al. 2002; Ramon et al. 2004; Wolters et al. 2006). For example, incorporating the cerebrospinal fluid or anisotropic conductivity in the white matter could affect the sensitivity measures examined in the present study. The regions where the deviations from the spherical model were largest, i.e., in frontal, insular, and medial regions, are in accordance with those found in previous studies comparing dipole localization accuracy in spherically symmetric volume conductor model vs. anatomically more realistic head models (Hämäläinen and Sarvas 1989; Haueisen et al. 2002; Schimpf et al. 2002; Tarkiainen et al. 2003; Wolters et al. 2006). We note that, even though Rλ being close to zero does not guarantee that the spherical model will provide an accurate MEG forward solution for those locations, previous studies have suggested that often the spherical model does indeed provide an accurate forward model for MEG, especially for most of the superficial cortical regions (Huang et al. 1999).
Some caution is necessary in interpreting our results, as numerical accuracy of the forward model for sources close to the skull may suffer from the sparseness of the BEM vertices (Fuchs et al. 2001). Previous studies have suggested that the resolution used here provides a reasonable accuracy for most cortical source locations (Yvert et al. 1996). Furthermore, we selected the inner surface of the cerebral cortex, rather than, e.g., the middle of the gray matter, for the location of the sources, to avoid them being very close to the brain–skull interface when computing the forward model. Despite these precautions, we observed large values of λk,min at some isolated superficial vertices, which presumably were caused by numerical instabilities.
Silent source elements are a special case of silent source distributions defined by the null space of the forward matrix (Mosher et al. 1999; Riera et al. 2006). Extensive coverage of the scalp by the sensor arrays, as well as combined use of EEG and MEG sensors help to minimize the null space. In general, the orientation of the source has a strong influence on the spatial pattern of the EEG and MEG signals. The importance of the source orientation for the interpretation of the data is sometimes overlooked in the common practice of visualizing only the estimated source magnitudes of inverse solutions.
Overall, the present results agree with the commonly held view that an important difference between MEG and EEG is that MEG is mainly sensitive to the tangential source components only, whereas EEG is sensitive to all components. Few cortical sources, however, are expected to have exactly the orientation of the lowest sensitivity in MEG (Hillebrand and Barnes 2002). Therefore, the depth, rather than the orientation, is likely to be the critical factor determining the detectability of an individual current dipole source in MEG. In other words, the sensitivity of MEG to source having the orientation corresponding to λk,min is low at all cortical locations, whereas the sensitivity to the other two orientations depends strongly on the depth (Fig. 2). This is consistent with MEG being most sensitive to superficial sources, and the sensitivity being much reduced for deep sources (Cohen and Cuffin 1983; de Jongh et al. 2005; Goldenholz et al. 2009; Hillebrand and Barnes 2002).
The insensitivity of MEG to one source orientation has a specific consequence for extended patches of cortical activity. For extended sources, selective cancellation of MEG and EEG signals generated by tangential sources may occur (Ahlfors et al. 2010; Eulitz et al. 1997; Freeman et al. 2009). The selective cancellation can affect the relative signal-to-noise ratio (SNR) of MEG and EEG (de Jongh et al. 2005; Goldenholz et al. 2009; Huiskamp et al. 2010). Several studies have reported that epileptic spikes are sometimes detected in either EEG or MEG but not necessarily in both, even when EEG and MEG are recorded simultaneously (e.g., Knake et al. 2006; Lin et al. 2003). Spikes detected in EEG only could be due to sources corresponding to the orientation of the lowest sensitivity in MEG. Spikes that are detected only in MEG, however, cannot be easily explained by the orientation of the source alone. One reason could be spatial under-sampling due to a low number of EEG sensors. More likely, however, the SNR for these spikes is higher in MEG than in EEG because of selective cancellation of background brain noise from tangentially oriented sources, such that in EEG the spikes are obscured by radially oriented background sources (Ahlfors et al. 2010). Therefore, the presence of a source orientation to which MEG is insensitive contributes to the complementary properties of MEG and EEG beyond the case of a single focal source going undetected in MEG because it is exactly radially oriented.
This work was supported by NIH grants NS057500, NS037462, and HD040712. This work was supported in part by The National Center for Research Resources (P41RR14075) and the MIND Research Network. We thank Drs. David Cohen and Maria Mody for valuable discussions and comments on the manuscript.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.