Abstract
Magnetoencephalography (MEG), which acquires neuromagnetic fields in the brain, is a useful diagnostic tool in presurgical evaluation of epilepsy. Previous studies have shown that MEG affects the planning intracranial electroencephalography placement and correlates with surgical outcomes by using a single dipole model. Spatiotemporal source analysis using distributed source models is an advanced method for analyzing MEG, and has been recently introduced for analyzing epileptic spikes. It has advantages over the conventional single dipole analysis for obtaining accurate sources and understanding the propagation of epileptic spikes. In this article, we review the source analysis methods, describe the techniques of the distributed source analysis, interpretation of source distribution maps, and discuss the benefits and feasibility of this method in evaluation of epilepsy.
Introduction
Magnetoencephalography (MEG) is an important, non-invasive diagnostic tool, which acquires neuromagnetic fields generated in the brain with high spatial and temporal resolution. Clinical usefulness of MEG, especially in presurgical evaluation of epilepsy, is well documented in recent reviews (Stufflebeam et al., 2009; Stufflebeam, 2011). Currently, clinical applications of MEG are divided into two categories: (1) spontaneous brain activity analysis, including epileptic spike mapping, most often for determining an irritative zone, (2) mapping of eloquent cortex, such as primary motor cortex and language area for avoiding postsurgical functional deficits (Stufflebeam et al., 2009; Stufflebeam, 2011).
Source localization of MEG spikes is frequently performed in clinical practice for identifying an irritative zone (Otsubo and Snead, 2001; Chuang et al., 2006; Stufflebeam, 2011), providing spatial information of spike activities at the sensor level. Source analysis of MEG typically incorporates anatomical information derived from each individual’s magnetic resonance imaging (MRI), and calculates the sources of neural activities by applying a certain mathematical model to the measured magnetic fields. These cortical and subcortical sources are visualized on the MRI or MRI-based anatomical atlas, and provide current dipole distribution maps.
Calculating intracranial sources from MEG obtained outside the brain, an example of the inverse problem, is mathematically non-unique and ill-posed. Certain assumptions are necessary for providing the proper source modeling. Thus, many procedures of source analysis have been proposed, such as single dipole, multi-dipole, and distributed source models, which are also applied for the source analysis of scalp electroencephalography (EEG) (for review, see Michel et al., 2004; Plummer et al., 2008).
In this article, we review the source analysis methods, describe the techniques of the distributed source analysis, and discuss the feasibility of this method in evaluation of epileptic spikes.
Single Equivalent Current Dipole Analysis
Single equivalent current dipole (ECD) analysis has been widely used for source localization of epileptic spikes for decades (Otsubo and Snead, 2001; Stefan et al., 2003; Fischer et al., 2005; Chuang et al., 2006). This model assumes that a single dipole source generates all the neuromagnetic fields recorded on the sensors, and is considered physiologically plausible when a limited area of the cortex is synchronously activated. In the analysis, the measured magnetic fields at a given latency are modeled by the best-fitting single dipole. ECDs are typically calculated by using a standard iterative least-square algorithm (Marquardt, 1963; Iwasaki et al., 2002), and several indicators of their reliability are also calculated, such as goodness of fit (GOF) and correlation coefficient. These indicators reflect the concordance between the magnetic fields calculated from the ECD and the actual measurement MEG data. The current dipole moment is represented by the magnitude of the ECD. These indicators and other metrics are used for selecting adequate sources and discarding inadequate ECDs by setting a threshold. Adequate ECDs are mapped on the patient’s MRI, demonstrating the distribution of ECDs (Figure 1). Previous studies have validated ECD analysis in temporal lobe epilepsy (Baumgartner et al., 2000; Iwasaki et al., 2002; Assaf et al., 2004; Reinsberger et al., 2010) and frontal lobe epilepsy (Shiraishi et al., 2001; Genow et al., 2004; Ossenblok et al., 2007). Several studies have shown spatial concordance of ECD distribution and interictal spiking area on intracranial EEG (IEEG) (Mikuni et al., 1997; Oishi et al., 2002).
Figure 1
Although ECD analysis is well established procedure for localizing MEG spikes, a few issues remain. First, the criteria of selecting ECDs may vary from laboratory to laboratory. Single ECDs generally do not provide a GOF of 100%, i.e., 100% of measured magnetic fields can not be explained by a single ECD, due to the oversimplification of this model. Some studies accept ECDs with a GOF >80–90% or correlation coefficient >0.90 (Iwasaki et al., 2002; Genow et al., 2004; Pataraia et al., 2004; Oishi et al., 2006; Knowlton et al., 2009), but other use lower thresholds (Shiraishi et al., 2005a,b, 2011; Reinsberger et al., 2010). The criteria dramatically affect the number of adequate ECDs (Tanaka et al., 2009a), although these thresholds are determined subjectively. Several researchers selected certain sensors for calculating ECDs to obtain adequate ECDs (Iwasaki et al., 2002; Pataraia et al., 2004). Such a process adds another assumption that spikes only appear in a group of restricted sensors. Second, the ECD analysis is sometimes inaccurate (Kobayashi et al., 2005) such as when the signal-to-noise ratio (SNR) is low and the spike is widespread (Shiraishi et al., 2005a; Hara et al., 2007). This may obscure the precise localization of early ictal discharges, which usually have a low SNR (Tanaka et al., 2009a). Third, the single ECD model does not always provide an accurate representation of the time course of the epileptiform discharge. When ECDs are sequentially calculated from the onset to the peak of spike, the SNR is low in the early latency. ECDs obtained from the spike onset may have low GOF values and thus an unstable localization (Kanamori et al., 2013). Several clinical studies calculate ECDs at the spike peak (Iwasaki et al., 2002; Oishi et al., 2006; Reinsberger et al., 2010; Jin et al., 2013) for obtaining a high SNR. However, spikes may be widespread or may propagate, and those ECDs may not accurately identify the onset of the discharge. Moreover, recent studies have proposed that epilepsy is a network disease caused by abnormal neural networks (for review, see Spencer, 2002). Spike propagation is required for understanding the abnormal networks (Tanaka et al., 2010), and advanced analysis procedure is necessary for understanding the spike propagation.
Spatiotemporal Distributed Source Analysis
Distributed source models assume that a certain amount of cortical patches are simultaneously activated (Dale and Sereno, 1993; Hämäläinen and Ilmoniemi, 1994). This model calculates the source distribution by deploying numerous unit dipoles on the cortical surface. The source space is created by using a realistic head model, such as the boundary elemental model (Hämäläinen and Sarvas, 1989; Oostendorp and van Oosterom, 1989; Crouzeix et al., 1999; Fuchs et al., 2001), obtained by cortical surface reconstruction derived from MRI (Dale et al., 1999; Fischl et al., 1999). Theoretically, the source solution gives concordance between the simulated and measured magnetic fields by adjusting the strength and orientation of unit dipoles (Dale and Sereno, 1993; Hämäläinen and Ilmoniemi, 1994). The source distribution can be calculated at any time points of spike, and from which dynamic activation maps along with the time course can be created. Thus, the source distribution maps represent the cortical activation generated by the spikes both spatially and temporally. The distributed source analysis does not require an assumption on the number of dipoles and thresholds of likelihood parameters as used in the single dipole analysis. However, an infinite number of source distributions can generate a similar magnetic field pattern, and further assumption, such as source distribution with minimum overall intensity (L2-norm), is necessary for determining the optimal solution. Thus, many types of analysis have been proposed based on different assumption and modeling, referred as minimum norm estimate (MNE) (Hämäläinen and Ilmoniemi, 1994), dynamic statistical parametric mapping (Dale et al., 2000), Laplacian weighted minimum norm (LORETA) (Pasqual-Marqui et al., 1994), local autoregressive average (LAURA), and EPIFOCUS (Grave de Peralta et al., 2001). Here, we use the MNE as an example of a distributed source solution.
MNE Analysis in Practice
This section explains spike analysis with MNE software (www.mne.org), which mainly provides an MNE and dSPM solution algorithm.
MEG preprocessing – spike selection and clean up
The purpose of preprocessing is to obtain spike data from raw MEG with minimum artifacts. Both individual and averaged spikes can be analyzed by distributed source models. Note that each of these spikes provides identical series of source distribution maps. Visual inspection of spikes is widely performed while automated spike detection is also introduced (Ossadtchi et al., 2004). However, the MEG spike morphology has not been sufficiently described while EEG spikes were well documented in the literature (International Federation of Societies for Clinical Neurophysiology, 1974). A clinical guideline has recommended identifying the MEG spikes based on the principals established for EEG (Bagic et al., 2011). Clarification of MEG spikes will be useful for clinical application of MEG especially in patients with negative EEG findings. Spike selection is critical for obtaining appropriate averaged spikes and spike morphology must be considered.
Artifact reduction is also an important issue for calculating adequate sources of spikes. Exclusion of sensors (channels), which contain continuous artifacts, is helpful for avoiding inadequate affection on the source analysis. Independent component analysis may be useful for removing artifacts specific patterns (e.g., ECG) as used on EEG (Kobayashi et al., 2001). Taulu et al. (2004) (Taulu and Simola, 2006; Taulu and Hari, 2009) proposed signal source separation (SSS) and its temporal extension (tSSS). SSS decomposes MEG signals into two components, which are generated from sources inside and outside of the sensor space, and remove the latter component. The temporally extended SSS also considers temporal signal correlation and is widely used in clinical practice (Medvedovsky et al., 2009; Song et al., 2009; Tanaka et al., 2009b; Jin et al., 2013; Kakisaka et al., 2013) (Figure 2). Motion artifacts distort geographic information of the head and sensors, resulting in unreliable source localization. These movements can be compensated by tracking the head position during acquisition, and this technique is useful for analyzing ictal MEG data, where head deviation sometimes occur during seizures (Medvedovsky et al., 2007; Kakisaka et al., 2012).
Figure 2
MRI preprocessing – creating a head model
Minimum norm estimate constrains source activities to the cortical surface images. The cortical surface is reconstructed from anatomical T1 MRI data, and the reconstruction is the first step of MRI processing. Various software packages, such as Freesurfer (Dale et al., 1999; Fischl et al., 1999) and Brainstorm (Tadel et al., 2011), provide cortical surface reconstructions. The source space is created by using the cortical surface, deploying grid spacing with numerous cortical patches (Crouzeix et al., 1999; Fuchs et al., 2001). Unit current dipoles are distributed in the source space, and the boundary elemental method (BEM) creates a head model for calculating the activation of these dipoles (Hämäläinen and Sarvas, 1989; Oostendorp and van Oosterom, 1989; Crouzeix et al., 1999; Fuchs et al., 2001). A single-layer BEM model is generally used, since the neuromagnetic signals are not affected by the tissue conductivity (Hämäläinen and Sarvas, 1989).
Coregistration of the patient’s MRI with the MEG sensor space is typically performed through the digitization information of the head and head position indicator that generates artificial electric currents in the sensor array. This process estimates the relationship between the source and the sensor spaces (Figure 3).
Figure 3
Inverse solution, creating source maps
The forward solution, which models the magnetic signals generated by unit current dipoles, is obtained by using the coregistration (Hämäläinen and Sarvas, 1989; Oostendorp and van Oosterom, 1989). Inverse solution is calculated based on the forward solution, mapping the strength of each unit dipole (Dale and Sereno, 1993; Dale et al., 2000). Source activation is projected on the cortical surface by applying a certain threshold, showing the source strength with different colors (Figure 4). The strength and extent of activation change along with the time course (Figure 4).
Figure 4
Interpretation of Distributed Source Maps
Distributed source maps can be examined both spatially and temporally; these maps estimate the spatial source distribution and its time course. A cortical parcellation, provided by Freesurfer (Dale et al., 1999; Fischl et al., 1999), is helpful for understanding detailed anatomical location of activation (Figure 5). More importantly, spatial source distribution is considered representing to the cortical extent of spike involvement; activation in a small cortical area suggests that the spiking area is restricted (Hara et al., 2007; Tanaka et al., 2013b). Widespread cortical activation, extending over two lobes or bilaterally, may reflect a large abnormal neural network associated with epilepsy (Shiraishi et al., 2005a,b, 2011).
Figure 5
Time course of source activation provides a way to evaluate spatial source distribution at each time point of the spike (Tanaka et al., 2010, 2012). The interpreter can trace the activation pattern from the spike onset to the peak. The changing pattern along with the time course may reflect physiologic spiking process, such as growing or propagation. Here, several patterns of spatiotemporal source distribution can be recognized; (1) restricted onset activation + restricted peak activation in the onset area, suggesting a highly limited spike involvement, (2) restricted onset activation + widespread peak activation, suggesting broad spike propagation, (3) widespread onset activation + widespread peak activation, suggesting non-localizing, widespread spiking in the cortex, (4) restricted early activation + restricted late activation in the distant area, suggesting an abnormal propagation pathway (Figure 6).
Figure 6
There are still several issues for assuming that the MNE maps accurately show the spike distribution in the whole brain; (1) MEG has different sensitivity to intracranial electric currents depending on regions (de Jongh et al., 2005; Goldenholz et al., 2009), and the information in the region with low sensitivity may be missing. (2) MNE solution does not consider the anatomical connection through white matter tracts. Tractography using diffusion tensor imaging will be useful in combination with MNE maps (Tanaka et al., 2012). (3) There is no established way to determine the threshold objectively, although such a threshold greatly affects the appearance of source maps. Several studies have introduced thresholds which are determined by a quantitative procedure (Gallagher et al., 2012; Tanaka et al., 2013a,b). Further investigation will be necessary to understand how the threshold should be determined.
Benefits of Distributed Source Maps
Previous studies using distributed source analysis mainly focused on providing the evidence of accurate source localization in comparison with IEEG (Waberski et al., 2000; Zumsteg et al., 2006; Tanaka et al., 2010; Kanamori et al., 2013). There are only a limited number of studies describing the benefit of this method in clinical practice (Tanaka et al., 2009a, 2013b). However, the results of these studies suggest benefits to use this technique in clinical practice.
One of the possible benefits is to obtain accurate spike localization. Distributed source analysis likely provides more reasonable solution than a single dipole model, although there is a localization error still observed in various distributed models (Waberski et al., 2000; Silva et al., 2004; Soufflet and Boeijinga, 2005). Previous studies have shown that MNE (and its derivative, dSPM) provides more accurate source localization than single dipoles, by comparing with single photon emission tomography (Shiraishi et al., 2005a), surgical outcome (Tanaka et al., 2009a), and IEEG (Kanamori et al., 2013) in a small group of patients. By using the distributed source analysis, MEG may contribute to the presurgical evaluation of epilepsy more effectively.
The other benefit is to understand the spike propagation. Some epileptic spikes originate from a restricted onset and propagate to other cortical areas (Alarcon et al., 1994, 1997). In these spikes, source localization of the early phase may be more informative than analyzing the peak for identifying the spike origin. Single dipole method sometimes does not provide adequate sources at the early latency as described above (Kanamori et al., 2013), whereas distributed source analysis provides reliable source distribution which can reconstruct the original, small signals of spikes at the sensor level. This ability is also useful for analyzing ictal MEG, which shows only small discharges in the early phase of seizures (Tanaka et al., 2009a). Distributed source analysis has nicely shown the possible onset of spikes with widespread cortical involvement at the peak in the previous studies (Shiraishi et al., 2005a,b; Kanamori et al., 2013).
Future Directions
Single dipole method has been potentially investigated and its benefits and limitations are well understood, and now several laboratories are doing similar studies with MNE and other distributed source solutions. Recent studies have validated MNE source distribution with IEEG. Tanaka et al. (2010) demonstrated that MNE analysis of frontotemporal MEG spikes accurately represents spike propagation as observed in IEEG. Frontoparietal and temporoparietal propagation patterns are also consistent between MEG and IEEG in a series of cases (Kanamori et al., 2013). Such validation will be highly desirable in other propagation involving various regions and in patients with various types of epilepsy.
An important distinction between ECD and MNE is that distributed source maps generally show the source localization of one single spike whereas single dipole maps project many dipoles obtained from different spikes. Thus, the single dipole method provides a viewpoint of a spike population, such as “clustered” or “scattered.” Distributed source maps do not have such mapping procedures that are widely used. Therefore, combined use of single dipole maps and distributed source maps is necessary in the current settings of spike analysis. Development of new mapping techniques, which overview numerous distributed source maps, will be useful for analyzing many spike populations.
Another area under active investigation is the surgical implications of distributed source maps. Several studies have demonstrated that MEG affects the planning IEEG placement and interpretation (Sutherling et al., 2008; Knowlton et al., 2009), and correlates with surgical outcomes (Iwasaki et al., 2002) by using a single dipole model. On the other hand, a recent study has shown that these propagation patterns are highly correlated with surgical outcomes in patients with temporal lobe epilepsy by using MNE (Tanaka et al., 2013b). Comparison of these techniques regarding with surgical outcomes is now becoming better understood.
Conclusion
Spatiotemporal distributed source analysis is highly useful for understanding epileptic spikes. It provides more accurate source localization than a single dipole model in some situations, and may be informative in presurgical evaluation of epilepsy. However, further observations are necessary for establishing its usefulness in clinical practice.
Statements
Acknowledgments
This work was supported by National Institutes of Health (S10RR014978, 1R01NS069696-01A1) and NIBIB (P41EB015896).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AlarconG.Garcia SeoaneJ. J.BinnieC. D.Martin MiguelM. C.JulerJ.PolkeyC. E.et al (1997). Origin and propagation of interictal discharges in the acute electrocorticogram. Implications for pathophysiology and surgical treatment of temporal lobe epilepsy. Brain120, 2259–2282.10.1093/brain/120.12.2259
2
AlarconG.GuyC. N.BinnieC. D.WalkerS. R.ElwesR. D.PolkeyC. E. (1994). Intracerebral propagation of interictal activity in partial epilepsy: implications for source localisation. J. Neurol. Neurosurg. Psychiatry57, 435–449.10.1136/jnnp.57.4.435
3
AssafB. A.KarkarK. M.LaxerK. D.GarciaP. A.AustinE. J.BarbaroN. M.et al (2004). Magnetoencephalography source localization and surgical outcome in temporal lobe epilepsy. Clin. Neurophysiol.115, 2066–2076.10.1016/j.clinph.2004.04.020
4
BagicA. I.KnowltonR. C.RoseD. F.EbersoleJ. S.ACMEGS Clinical Practice Guideline (CPG) Committee. (2011). American clinical magnetoencephalography society clinical practice guideline 1: recording and analysis of spontaneous cerebral activity. J. Clin. Neurophysiol.28, 348–354.10.1097/WNP.0b013e3182272fed
5
BaumgartnerC.PataraiaE.LindingerG.DeeckeL. (2000). Neuromagnetic recordings in temporal lobe epilepsy. J. Clin. Neurophysiol.17, 177–189.10.1097/00004691-200003000-00007
6
ChuangN. A.OtsuboH.PangE. W.ChuangS. H. (2006). Pediatric magnetoencephalography and magnetic source imaging. Neuroimaging Clin. N. Am.16, 193–210.10.1016/j.nic.2005.11.001
7
CrouzeixA.YvertB.BertrandO.PernierJ. (1999). An evaluation of dipole reconstruction accuracy with spherical head model and realistic head models in MEG. Clin. Neurophysiol.110, 2176–2188.10.1016/S1388-2457(99)00174-1
8
DaleA. M.FischlB.SerenoM. I. (1999). Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage9, 179–194.10.1006/nimg.1998.0395
9
DaleA. M.LiuA. K.FischlB. R.BucknerR. L.BelliveauJ. W.LewineJ. D.et al (2000). Dynamic statistical parametric mapping: combining fMRI and MEG for high-resolution imaging of cortical activity. Neuron26, 55–67.10.1016/S0896-6273(00)81138-1
10
DaleA. M.SerenoM. I. (1993). Improved localization of cortical activity by combining EEG and MEG with MRI cortical surface reconstruction: a linear approach. J. Cogn. Neurosci.5, 162–176.10.1162/jocn.1993.5.2.162
11
de JonghA.de MunckJ. C.GonçalvesS. I.OssenblokP. (2005). Differences in MEG/EEG epileptic spike yields explained by regional differences in signal-to-noise ratios. J. Clin. Neurophysiol.22, 153–158.10.1097/01.WNP.0000158947.68733.51
12
FischerM. J.SchelerG.StefanH. (2005). Utilization of magnetoencephalography results to obtain favorable outcomes in epilepsy surgery. Brain128, 153–157.10.1093/brain/awh333
13
FischlB.SerenoM. I.DaleA. M. (1999). Cortical surface-based analysis. II: inflation, flattening, and a surface-based coordinate system. Neuroimage9, 195–207.10.1006/nimg.1998.0396
14
FuchsM.WagnerM.KastnerJ. (2001). Boundary element method volume conductor models for EEG source reconstruction. Clin. Neurophysiol.112, 1400–1407.10.1016/S1388-2457(01)00589-2
15
GallagherA.TanakaN.SuzukiN.LiuH.ThieleE. A.StufflebeamS. M. (2012). Decreased language laterality in tuberous sclerosis complex: a relationship between language dominance and tuber location as well as history of epilepsy. Epilepsy Behav.25, 36–41.10.1016/j.yebeh.2012.06.013
16
GenowA.HummelC.SchelerG.HopfengärtnerR.KaltenhäuserM.BuchfelderM.et al (2004). Epilepsy surgery, resection volume and MSI localization in lesional frontal lobe epilepsy. Neuroimage21, 444–449.10.1016/j.neuroimage.2003.08.029
17
GoldenholzD. M.AhlforsS. P.HämäläinenM. S.SharonD.IshitobiM.VainaL. M.et al (2009). Mapping the signal-to-noise-ratios of cortical sources in magnetoencephalography and electroencephalography. Hum. Brain Map.30, 1077–1086.10.1002/hbm.20571
18
Grave de PeraltaR.GonzalezS.LantzG.MichelC. M.LandisT. (2001). Noninvasive localization of electromagnetic epileptic activity. I. Method descriptions and simulations. Brain Topogr.14, 131–137.10.1023/A:1012944913650
19
HämäläinenM. S.IlmoniemiR. J. (1994). Interpreting magnetic fields of the brain: minimum norm estimates. Med. Biol. Eng. Comput.32, 35–42.10.1007/BF02512476
20
HämäläinenM. S.SarvasJ. (1989). Realistic conductivity geometry model of the human head for interpretation of neuromagnetic data. IEEE Trans. Biomed. Eng.36, 165–171.10.1109/10.16463
21
HaraK.LinF. H.CamposanoS.FoxeD. M.GrantP. E.BourgeoisB. F.et al (2007). Magnetoencephalographic mapping of interictal spike propagation: a technical and clinical report. AJNR Am. J. Neuroradiol.28, 1486–1488.10.3174/ajnr.A0596
22
International Federation of Societies for Clinical Neurophysiology. (1974). A glossary of terms most commonly used by clinical electroencephalographers. Electroencephalogr. Clin. Neurophysiol.37, 538–548.10.1016/0013-4694(74)90099-6
23
IwasakiM.NakasatoN.ShamotoH.NagamatsuK.KannoA.HatanakaK.et al (2002). Surgical implications of neuromagnetic spike localization in temporal lobe epilepsy. Epilepsia43, 415–424.10.1046/j.1528-1157.2002.30801.x
24
JinK.AlexopoulosA. V.MosherJ. C.BurgessR. C. (2013). Implanted medical devices or other strong sources of interference are not barriers to magnetoencephalographic recordings in epilepsy patients. Clin. Neurophysiol.124, 1283–1289.10.1016/j.clinph.2013.04.004
25
KakisakaY.MosherJ. C.WangZ. I.JinK.DubarryA. S.AlexopoulosA. V.et al (2013). Utility of temporally-extended signal space separation algorithm for magnetic noise from vagal nerve stimulators. Clin. Neurophysiol.124, 1277–1282.10.1016/j.clinph.2012.03.082
26
KakisakaY.WangZ. I.MosherJ. C.DubarryA. S.AlexopoulosA. V.EnatsuR.et al (2012). Clinical evidence for the utility of movement compensation algorithm in magnetoencephalography: successful localization during focal seizure. Epilepsy Res.101, 191–196.10.1016/j.eplepsyres.2012.03.014
27
KanamoriY.ShigetoH.HironagaN.HagiwaraK.UeharaT.ChataniH.et al (2013). Minimum norm estimates in MEG can delineate the onset of interictal epileptic discharges: a comparison with ECoG findings. Neuroimage Clin.2, 663–669.10.1016/j.nicl.2013.04.008
28
KnowltonR. C.RazdanS. N.LimdiN.ElgavishR. A.KillenJ.BlountJ.et al (2009). Effect of epilepsy magnetic source imaging on intracranial electrode placement. Ann. Neurol.65, 716–723.10.1002/ana.21660
29
KobayashiK.MerletI.GotmanJ. (2001). Separation of spikes from background by independent component analysis with dipole modeling and comparison to intracranial recording. Clin. Neurophysiol.112, 405–413.10.1016/S1388-2457(01)00457-6
30
KobayashiK.YoshinagaH.OhtsukaY.GotmanJ. (2005). Dipole modeling of epileptic spikes can be accurate or misleading. Epilepsia46, 397–408.10.1111/j.0013-9580.2005.31404.x
31
MarquardtD. W. (1963). An algorithm for least-squares estimation of nonlinear parameters. SIAM J. Appl. Math.11, 431–441.10.1137/0111030
32
MedvedovskyM.TauluS.BikmullinaR.AhonenA.PaetauR. (2009). Fine tuning the correlation limit of spatio-temporal signal space separation for magnetoencephalography. J. Neurosci. Methods177, 203–211.10.1016/j.jneumeth.2008.09.035
33
MedvedovskyM.TauluS.BikmullinaR.PaetauR. (2007). Artifact and head movement compensation in MEG. Neurol. Neurophysiol. Neurosci.29, 4.
34
MichelC. M.MurrayM. M.LantzG.GonzalezS.SpinelliL.Grave de PeraltaR. (2004). EEG source imaging. Clin. Neurophysiol.115, 2195–2222.10.1016/j.clinph.2004.06.001
35
MikuniN.NagamineT.IkedaA.TeradaK.TakiW.KimuraJ.et al (1997). Simultaneous recording of epileptiform discharges by MEG and subdural electrodes in temporal lobe epilepsy. Neuroimage5, 298–306.10.1006/nimg.1997.0272
36
OishiM.KameyamaS.MasudaH.TohyamaJ.KanazawaO.SasagawaM.et al (2006). Single and multiple clusters of magnetoencephalographic dipoles in neocortical epilepsy: significance in characterizing the epileptogenic zone. Epilepsia47, 355–364.10.1111/j.1528-1167.2006.00428.x
37
OishiM.OtsuboH.KameyamaS.MorotaN.MasudaH.KitayamaM.et al (2002). Epileptic spikes: magnetoencephalography versus simultaneous electrocorticography. Epilepsia43, 1390–1395.10.1046/j.1528-1157.2002.10702.x
38
OostendorpT. F.van OosteromA. (1989). Source parameter estimation in inhomogeneous volume conductors of arbitrary shape. IEEE Trans. Biomed. Eng.36, 382–391.10.1109/10.19859
39
OssadtchiA.BailletS.MosherJ. C.ThyerleiD.SutherlingW.LeahyR. M. (2004). Automated interictal spike detection and source localization in magnetoencephalography using independent components analysis and spatio-temporal clustering. Clin. Neurophysiol.115, 508–522.10.1016/j.clinph.2003.10.036
40
OssenblokP.de MunckJ. C.ColonA.DrolsbachW.BoonP. (2007). Magnetoencephalography is more successful for screening and localizing frontal lobe epilepsy than electroencephalography. Epilepsia48, 2139–2149.10.1111/j.1528-1167.2007.01223.x
41
OtsuboH.SneadO. C.III (2001). Magnetoencephalography and magnetic source imaging in children. J. Child Neurol.16, 227–235.10.2310/7010.2001.17982
42
Pasqual-MarquiR. D.MichelC. M.LehmannD. (1994). Low resolution electromagnetic tomography: a new method to localize electrical activity in the brain. Int. J. Psychophysiol.18, 49–65.10.1016/0167-8760(84)90014-X
43
PataraiaE.SimosP. G.CastilloE. M.BillingsleyR. L.SarkariS.WhelessJ. W.et al (2004). Does magnetoencephalography add to scalp video-EEG as a diagnostic tool in epilepsy surgery?Neurology62, 943–948.10.1212/01.WNL.0000115122.81621.FE
44
PlummerC.HarveyA. S.CookM. (2008). EEG source localization in focal epilepsy: where are we now?Epilepsia49, 201–218.10.1111/j.1528-1167.2007.01381.x
45
ReinsbergerC.TanakaN.ColeA. J.LeeJ. W.DworetzkyB. A.BromfieldE. B.et al (2010). Current dipole orientation and distribution of epileptiform activity correlates with cortical thinning in left mesiotemporal epilepsy. Neuroimage52, 1238–1242.10.1016/j.neuroimage.2010.04.264
46
ShiraishiH.AhlforsS. P.StufflebeamS. M.KnakeS.LarssonP. G.HämäläinenM. S.et al (2011). Comparison of three methods for localizing interictal epileptiform discharges with magnetoencephalography. J. Clin. Neurophysiol.28, 431–440.10.1097/WNP.0b013e318231c86f
47
ShiraishiH.AhlforsS. P.StufflebeamS. M.TakanoK.OkajimaM.KnakeS.et al (2005a). Application of magnetoencephalography in epilepsy patients with widespread spike or slow-wave activity. Epilepsia46, 1264–1272.10.1111/j.1528-1167.2005.65504.x
48
ShiraishiH.StufflebeamS. M.KnakeS.AhlforsS. P.SudoA.AsahinaN.et al (2005b). Dynamic statistical parametric mapping for analyzing the magnetoencephalographic epileptiform activity in patients with epilepsy. J. Child Neurol.20, 363–369.10.1177/08830738050200041601
49
ShiraishiH.WatanabeY.WatanabeM.InoueY.FujiwaraT.YagiK. (2001). Interictal and ictal magnetoencephalographic study in patients with medial frontal lobe epilepsy. Epilepsia42, 875–882.10.1046/j.1528-1157.2001.042007875.x
50
SilvaC.MaltezJ. C.TrindadeE.ArriagaA.Ducla-SoaresE. (2004). Evaluation of L1 and L2 minimum norm performances on EEG localizations. Clin. Neurophysiol.115, 1657–1668.10.1016/j.clinph.2004.02.009
51
SongT.CuiL.GaaK.FefferL.TauluS.LeeR. R.et al (2009). Signal space separation algorithm and its application on suppressing artifacts caused by vagus nerve stimulation for magnetoencephalography recordings. J. Clin. Neurophysiol.26, 392–400.10.1097/WNP.0b013e3181c29896
52
SouffletL.BoeijingaP. H. (2005). Linear inverse solutions: simulations from a realistic head model in MEG. Brain Topogr.18, 87–99.10.1007/s10548-005-0278-6
53
SpencerS. S. (2002). Neural networks in human epilepsy: evidence of and implications for treatment. Epilepsia43, 219–227.10.1046/j.1528-1157.2002.26901.x
54
StefanH.HummelC.SchelerG.GenowA.DruschkyK.TilzC.et al (2003). Magnetic brain source imaging of focal epileptic activity: a synopsis of 455 cases. Brain126, 2396–2405.10.1093/brain/awg239
55
StufflebeamS. M. (2011). Clinical magnetoencephalography for neurosurgery. Neurosurg. Clin. N. Am.22, 153–167.10.1016/j.nec.2010.11.006
56
StufflebeamS. M.TanakaN.AhlforsS. P. (2009). Clinical applications of magnetoencephalography. Hum. Brain Mapp.30, 1813–1823.10.1002/hbm.20792
57
SutherlingW. W.MamelakA. N.ThyerleiD.MaleevaT.MinazadY.PhilpottL.et al (2008). Influence of magnetic source imaging for planning intracranial EEG in epilepsy. Neurology71, 990–996.10.1212/01.wnl.0000326591.29858.1a
58
TadelF.BailletS.MosherJ. C.PantazisD.LeahyR. M. (2011). Brainstorm: a user-friendly application for MEG/EEG analysis. Comput. Intell. Neurosci.2011, 879716.10.1155/2011/879716
59
TanakaN.ColeA. J.von PechmannD.WakemanD. G.HämäläinenM. S.LiuH.et al (2009a). Dynamic statistical parametric mapping for analyzing ictal magnetoencephalographic spikes in patients with intractable frontal lobe epilepsy. Epilepsy Res.85, 279–286.10.1016/j.eplepsyres.2009.03.023
60
TanakaN.ThieleE. A.MadsenJ. R.BourgeoisB. F.StufflebeamS. M. (2009b). Magnetoencephalographic analysis in patients with vagus nerve stimulator. Pediatr. Neurol.41, 383–387.10.1016/j.pediatrneurol.2009.06.007
61
TanakaN.GrantP. E.SuzukiN.MadsenJ. R.BerginA. M.HämäläinenM. S.et al (2012). Multimodal imaging of spike propagation: a technical case report. AJNR Am. J. Neuroradiol.33, E82–E84.10.3174/ajnr.A2701
62
TanakaN.HämäläinenM. S.AhlforsS. P.LiuH.MadsenJ. R.BourgeoisB. F.et al (2010). Propagation of epileptic spikes reconstructed from spatiotemporal magnetoencephalographic and electroencephalographic source analysis. Neuroimage50, 217–222.10.1016/j.neuroimage.2009.12.033
63
TanakaN.LiuH.ReinsbergerC.MadsenJ. R.BourgeoisB. F.DworetzkyB. A.et al (2013a). Language lateralization represented by spatiotemporal mapping of magnetoencephalography. AJNR Am. J. Neuroradiol.34, 558–563.10.3174/ajnr.A3233
64
TanakaN.PetersJ. M.ProhlA. K.TakayaS.MadsenJ. R.BourgeoisB. F.et al (2013b). Clinical value of magnetoencephalographic spike propagation represented by spatiotemporal source analysis: correlation with surgical outcome. Epilepsy Res.10.1016/j.eplepsyres.2013.11.006 Available from: http://dx.doi.org/10.1016/j.eplepsyres.2013.11.006,
65
TauluS.HariR. (2009). Removal of magnetoencephalographic artifacts with temporal signal-space separation: demonstration with single-trial auditory-evoked responses. Hum. Brain Mapp.30, 1524–1534.10.1002/hbm.20627
66
TauluS.KajolaM.SimolaJ. (2004). Suppression of interference and artifacts by the signal space separation method. Brain Topogr.16, 269–275.10.1023/B:BRAT.0000032864.93890.f9
67
TauluS.SimolaJ. (2006). Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurements. Phys. Med. Biol.51, 1759–1768.10.1088/0031-9155/51/7/008
68
WaberskiT. D.GobbeléR.HerrendorfG.SteinhoffB. J.KolleR.FuchsM.et al (2000). Source reconstruction of mesial-temporal epileptiform activity: comparison of inverse techniques. Epilepsia41, 1574–1583.10.1111/j.1499-1654.2000.001574.x
69
ZumstegD.FriedmanA.WieserH. G.WennbergR. A. (2006). Propagation of interictal discharges in temporal lobe epilepsy: correlation of spatiotemporal mapping with intracranial foramen ovale electrode recordings. Clin. Neurophysiol.117, 2615–2626.10.1016/j.clinph.2006.07.319
Summary
Keywords
magnetoencephalography, epilepsy, distributed source analysis, spike propagation, minimum norm estimate, epilepsy surgery
Citation
Tanaka N and Stufflebeam SM (2014) Clinical Application of Spatiotemporal Distributed Source Analysis in Presurgical Evaluation of Epilepsy. Front. Hum. Neurosci. 8:62. doi: 10.3389/fnhum.2014.00062
Received
30 November 2013
Accepted
25 January 2014
Published
10 February 2014
Volume
8 - 2014
Edited by
Christos Papadelis, Harvard Medical School, USA
Reviewed by
Mohamad Zakaria Koubeissi, George Washington University, USA; Yosuke Kakisaka, Tohoku University, Japan
Copyright
© 2014 Tanaka and Stufflebeam.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Naoaki Tanaka, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, 149 Thirteenth Street, Suite 2301, Charlestown, MA 02129, USA e-mail: naoro@nmr.mgh.harvard.edu
This article was submitted to the journal Frontiers in Human Neuroscience.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.