Electrophysiological imaging of brain activity and connectivity-challenges and opportunities

Detecting functional hubs of ictogenic networks

Quantitative EEG (qEEG) has modified our understanding of epileptic seizures, shifting our view from the traditionally accepted hyper-synchrony paradigm toward more complex models based on re-organization of functional networks. However, qEEG measurements are so far rarely considered during the clinical decision-making process. To better understand the dynamics of intracranial EEG signals, we examine a functional network derived from the quantification of information flow between intracranial EEG signals. Using transfer entropy, we analyzed 198 seizures from 27 patients undergoing pre-surgical evaluation for pharmaco-resistant epilepsy. During each seizure we considered for each network the in-, out- and total "hubs", defined respectively as the time and the EEG channels with the maximal incoming, outgoing or total (bidirectional) information flow. In the majority of cases we found that the hubs occur around the middle of seizures, and interestingly not at the beginning or end, where the most dramatic EEG signal changes are found by visual inspection. For the patients who then underwent surgery, good postoperative clinical outcome was on average associated with a higher percentage of out- or total-hubs located in the resected area (for out-hubs p = 0.01, for total-hubs p = 0.04). The location of in-hubs showed no clear predictive value. We conclude that the study of functional networks based on qEEG measurements may help to identify brain areas that are critical for seizure generation and are thus potential targets for focused therapeutic interventions.

Simultaneous EEG-fMRI: trial level spatio-temporal fusion for hierarchically reliable information discovery

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) have been pursued in an effort to integrate complementary noninvasive information on brain activity. The primary goal involves better information discovery of the event-related neural activations at a spatial region of the BOLD fluctuation with the temporal resolution of the electrical signal. Many techniques and algorithms have been developed to integrate EEGs and fMRIs; however, the relative reliability of the integrated information is unclear. In this work, we propose a hierarchical framework to ensure the relative reliability of the integrated results and attempt to understand brain activation using this hierarchical ideal. First, spatial Independent Component Analysis (ICA) of fMRI and temporal ICA of EEG were performed to extract features at the trial level. Second, the maximal information coefficient (MIC) was adopted to temporally match them across the modalities for both linear and non-linear associations. Third, fMRI-constrained EEG source imaging was utilized to spatially match components across modalities. The simultaneously occurring events in the above two match steps provided EEG-fMRI spatial-temporal reliable integrated information, resulting in the most reliable components with high spatial and temporal resolution information. The other components discovered in the second or third steps provided second-level complementary information for flexible and cautious explanations. This paper contains two simulations and an example of real data, and the results indicate that the framework is a feasible approach to reveal cognitive processing in the human brain.

Closed-loop brain-machine-body interfaces for noninvasive rehabilitation of movement disorders

Traditional approaches for neurological rehabilitation of patients affected with movement disorders, such as Parkinson's disease (PD), dystonia, and essential tremor (ET) consist mainly of oral medication, physical therapy, and botulinum toxin injections. Recently, the more invasive method of deep brain stimulation (DBS) showed significant improvement of the physical symptoms associated with these disorders. In the past several years, the adoption of feedback control theory helped DBS protocols to take into account the progressive and dynamic nature of these neurological movement disorders that had largely been ignored so far. As a result, a more efficient and effective management of PD cardinal symptoms has emerged. In this paper, we review closed-loop systems for rehabilitation of movement disorders, focusing on PD, for which several invasive and noninvasive methods have been developed during the last decade, reducing the complications and side effects associated with traditional rehabilitation approaches and paving the way for tailored individual therapeutics. We then present a novel, transformative, noninvasive closed-loop framework based on force neurofeedback and discuss several future developments of closed-loop systems that might bring us closer to individualized solutions for neurological rehabilitation of movement disorders.

Whole-scalp EEG mapping of somatosensory evoked potentials in macaque monkeys

High-density scalp EEG recordings are widely used to study whole-brain neuronal networks in humans non-invasively. Here, we validate EEG mapping of somatosensory evoked potentials (SSEPs) in macaque monkeys (Macaca fascicularis) for the long-term investigation of large-scale neuronal networks and their reorganisation after lesions requiring a craniotomy. SSEPs were acquired from 33 scalp electrodes in five adult anaesthetized animals after electrical median or tibial nerve stimulation. SSEP scalp potential maps were identified by cluster analysis and identified in individual recordings. A distributed, linear inverse solution was used to estimate the intracortical sources of the scalp potentials. SSEPs were characterised by a sequence of components with unique scalp topographies. Source analysis confirmed that median nerve SSEP component maps were in accordance with the somatotopic organisation of the sensorimotor cortex. Most importantly, SSEP recordings were stable both intra- and interindividually. We aim to apply this method to the study of recovery and reorganisation of large-scale neuronal networks following a focal cortical lesion requiring a craniotomy. As a prerequisite, the present study demonstrated that a 300-mm² unilateral craniotomy over the sensorimotor cortex necessary to induce a cortical lesion, followed by bone flap repositioning, suture and gap plugging with calcium phosphate cement, did not induce major distortions of the SSEPs. In conclusion, SSEPs can be successfully and reproducibly recorded from high-density EEG caps in macaque monkeys before and after a craniotomy, opening new possibilities for the long-term follow-up of the cortical reorganisation of large-scale networks in macaque monkeys after a cortical lesion.

Lateralization and localization of epilepsy related hemodynamic foci using presurgical fMRI

OBJECTIVE

The aim was to develop a method for the purpose of localizing epilepsy related hemodynamic foci for patients suffering intractable focal epilepsy using task-free fMRI alone.

METHODS

We studied three groups of subjects: patients with intractable focal epilepsy, healthy volunteers performing motor tasks, and healthy volunteers in resting state. We performed spatial independent component analysis (ICA) on the fMRI alone data and developed a set of IC selection criteria to identify epilepsy related ICs. The method was then tested in the two healthy groups.

RESULTS

In seven out of the nine surgery patients, identified ICs were concordant with surgical resection. Our results were also consistent with presurgical evaluation of the remaining one patient without surgery and may explain why she was not suitable for resection treatment. In the motor task study of ten healthy subjects, our method revealed components with concordant spatial and temporal features as expected from the unilateral motor tasks. In the resting state study of seven healthy subjects, the method successfully rejected all components in four out of seven subjects as non-epilepsy related components.

CONCLUSION

These results suggest the lateralization and localization value of fMRI alone in presurgical evaluation for patients with intractable unilateral focal epilepsy.

SIGNIFICANCE

The proposed method is noninvasive in nature and easy to implement. It has the potential to be incorporated in current presurgical workup for treating intractable focal epilepsy patients.

An EEMD-IVA framework for concurrent multidimensional EEG and unidimensional kinematic data analysis

Joint blind source separation (JBSS) is a means to extract common sources simultaneously found across multiple datasets, e.g., electroencephalogram (EEG) and kinematic data jointly recorded during reaching movements. Existing JBSS approaches are designed to handle multidimensional datasets, yet to our knowledge, there is no existing means to examine common components that may be found across a unidimensional dataset and a multidimensional one. In this paper, we propose a simple, yet effective method to achieve the goal of JBSS when concurrent multidimensional EEG and unidimensional kinematic datasets are available, by combining ensemble empirical mode decomposition (EEMD) with independent vector analysis (IVA). We demonstrate the performance of the proposed method through numerical simulations and application to data collected from reaching movements in Parkinson's disease. The proposed method is a promising JBSS tool for real-world biomedical signal processing applications.

Effective connectivity of hippocampal neural network and its alteration in Mg2+-free epilepsy model

Understanding the connectivity of the brain neural network and its evolution in epileptiform discharges is meaningful in the epilepsy researches and treatments. In the present study, epileptiform discharges were induced in rat hippocampal slices perfused with Mg2+-free artificial cerebrospinal fluid. The effective connectivity of the hippocampal neural network was studied by comparing the normal and epileptiform discharges recorded by a microelectrode array. The neural network connectivity was constructed by using partial directed coherence and analyzed by graph theory. The transition of the hippocampal network topology from control to epileptiform discharges was demonstrated. Firstly, differences existed in both the averaged in- and out-degree between nodes in the pyramidal cell layer and the granule cell layer, which indicated an information flow from the pyramidal cell layer to the granule cell layer during epileptiform discharges, whereas no consistent information flow was observed in control. Secondly, the neural network showed different small-worldness in the early, middle and late stages of the epileptiform discharges, whereas the control network did not show the small-world property. Thirdly, the network connectivity began to change earlier than the appearance of epileptiform discharges and lasted several seconds after the epileptiform discharges disappeared. These results revealed the important network bases underlying the transition from normal to epileptiform discharges in hippocampal slices. Additionally, this work indicated that the network analysis might provide a useful tool to evaluate the neural network and help to improve the prediction of seizures.

Applying EEG phase synchronization measures to non-linearly coupled neural mass models

BACKGROUND

Recent neuroimaging analyses aim to understand how information is integrated across brain regions that have traditionally been studied in isolation; however, detecting functional connectivity networks in experimental EEG recordings is a non-trivial task.

NEW METHOD

We use neural mass models to simulate 10-s trials with coupling between 1-3 and 5-8s and compare how well three phase-based connectivity measures recover this connectivity pattern across a set of experimentally relevant conditions: variable oscillation frequency and power spectrum, feed forward connections with or without feedback, and simulated signals with and without volume conduction.

RESULTS

Overall, the results highlight successful detection of the onset and offset of significant synchronizations for a majority of the 28 simulated configurations; however, the tested phase measures sometimes differ in their sensitivity and specificity to the underlying connectivity.

COMPARISON WITH EXISTING METHODS

Prior work has shown that these phase measures perform well on signals generated by a computational model of coupled oscillators. In this work we extend previous studies by exploring the performance of these measures on a different class of computational models, and we compare the methods on 28 variations that capture a set of experimentally relevant conditions.

CONCLUSIONS

Our results underscore that no single phase synchronization measure is substantially better than all others, and experimental investigations will likely benefit from combining a set of measures together that are chosen based on both the experimental question of interest, the signal to noise ratio in the EEG data, and the approach used for statistical significance.

Long-range neural activity evoked by premotor cortex stimulation: a TMS/EEG co-registration study

The premotor cortex is one of the fundamental structures composing the neural networks of the human brain. It is implicated in many behaviors and cognitive tasks, ranging from movement to attention and eye-related activity. Therefore, neural circuits that are related to premotor cortex have been studied to clarify their connectivity and/or role in different tasks. In the present work, we aimed to investigate the propagation of the neural activity evoked in the dorsal premotor cortex using transcranial magnetic stimulation/electroencephalography (TMS/EEG). Toward this end, interest was focused on the neural dynamics elicited in long-ranging temporal and spatial networks. Twelve healthy volunteers underwent a single-pulse TMS protocol in a resting condition with eyes closed, and the evoked activity, measured by EEG, was compared to a sham condition in a time window ranging from 45 ms to about 200 ms after TMS. Spatial and temporal investigations were carried out with sLORETA. TMS was found to induce propagation of neural activity mainly in the contralateral sensorimotor and frontal cortices, at about 130 ms after delivery of the stimulus. Different types of analyses showed propagated activity also in posterior, mainly visual, regions, in a time window between 70 and 130 ms. Finally, a likely “rebounding” activation of the sensorimotor and frontal regions, was observed in various time ranges. Taken together, the present findings further characterize the neural circuits that are driven by dorsal premotor cortex activation in healthy humans.

Grand challenges in mapping the human brain: NSF workshop report

This report summarizes the outcomes of the NSF Workshop on Mapping and Engineering the Brain, held at Arlington, VA, during August 13-14, 2013. Three grand challenges were identified, including high spatiotemporal resolution neuroimaging, perturbation-based neuroimaging, and neuroimaging in naturalistic environments. It was highlighted that each grand challenge requires groundbreaking discoveries, enabling technologies, appropriate knowledge transfer, and multi- and transdisciplinary education and training for success.

Anesthesia differentially modulates spontaneous network dynamics by cortical area and layer

Anesthesia is widely used in medicine and research to achieve altered states of consciousness and cognition. Whereas changes to macroscopic cortical activity patterns by anesthesia measured at the spatial resolution of electroencephalography have been widely studied, modulation of mesoscopic and microscopic network dynamics by anesthesia remain poorly understood. To address this gap in knowledge, we recorded spontaneous mesoscopic (local field potential) and microscopic (multiunit activity) network dynamics in primary visual cortex (V1) and prefrontal cortex (PFC) of awake and isoflurane anesthetized ferrets (Mustela putoris furo). This approach allowed for examination of activity as a function of cortical area, cortical layer, and anesthetic depth with much higher spatial and temporal resolution than in previous studies. We hypothesized that a primary sensory area and an association cortical area would exhibit different patterns of network modulation by anesthesia due to their different functional roles. Indeed, we found effects specific to cortical area and cortical layer. V1 exhibited minimal changes in rhythmic structure with anesthesia but differential modulation of input layer IV. In contrast, anesthesia profoundly altered spectral power in PFC, with more uniform modulation across cortical layers. Our results demonstrate that anesthesia modulates spontaneous cortical activity in an area- and layer-specific manner. These finding provide the basis for 1) refining anesthesia monitoring algorithms, 2) reevaluating the large number of systems neuroscience studies performed in anesthetized animals, and 3) increasing our understanding of differential dynamics across cortical layers and areas.

Functional brain network organisation of children between 2 and 5 years derived from reconstructed activity of cortical sources of high-density EEG recordings

There is increasing interest in applying connectivity analysis to brain measures (Rubinov and Sporns, 2010), but most studies have relied on fMRI, which substantially limits the participant groups and numbers that can be studied. High-density EEG recordings offer a comparatively inexpensive easy-to-use alternative, but require channel-level connectivity analysis which currently lacks a common analytic framework and is very limited in spatial resolution. To address this problem, we have developed a new technique for studies of network development that overcomes the spatial constraint and obtains functional networks of cortical areas by using EEG source reconstruction with age-matched average MRI templates (He et al., 1999). In contrast to previously reported channel-level analysis, this approach provides information about the cortical areas most likely to be involved in the network as well as their functional relationship (Babiloni et al., 2005; De Vico Fallani et al., 2007). In this study, we applied source reconstruction with age-matched templates to task-free high-density EEG recordings in typically-developing children between 2 and 6 years of age (O'Reilly, 2012). Graph theory was then applied to the association strengths of 68 cortical regions of interest based on the Desikan-Killiany atlas. We found linear increases of mean node degree, mean clustering coefficient and maximum betweenness centrality between 2 years and 6 years of age. Characteristic path length was negatively correlated with age. The correlation of the network measures with age indicates network development towards more closely integrated networks similar to reports from other imaging modalities (Fair et al., 2008; Power et al., 2010). We also applied eigenvalue decomposition to obtain functional modules (Clayden et al., 2013). Connection strength within these modules did not change with age, and the modules resembled hub networks previously described for MRI (Hagmann et al., 2010; Power et al., 2010). The high temporal resolution of EEG additionally allowed us to distinguish between frequency bands potentially reflecting dynamic coupling between different neural oscillators. Generally, network parameters were similar for networks based on different frequency bands, but frequency band did emerge as a significant factor for clustering coefficient and characteristic path length. In conclusion, the current analysis shows that source reconstruction of high-density EEG recordings with appropriate head models offers a valuable tool for estimating network parameters in studies of brain development. The findings replicate the pattern of closer functional integration over development described for other imaging modalities (Fair et al., 2008; Power et al., 2010).

Compliant head probe for positioning electroencephalography electrodes and near-infrared spectroscopy optodes

A noninvasive head probe that combines near-infrared spectroscopy (NIRS) and electroencephalography (EEG) for simultaneous measurement of neural dynamics and hemodynamics in the brain is presented. It is composed of a compliant expandable mechanism that accommodates a wide range of head size variation and an elastomeric web that maintains uniform sensor contact pressure on the scalp as the mechanism expands and contracts. The design is intended to help maximize optical and electrical coupling and to maintain stability during head movement. Positioning electrodes at the inion, nasion, central, and preauricular fiducial locations mechanically shapes the probe to place 64 NIRS optodes and 65 EEG electrodes following the 10-5 scalp coordinates. The placement accuracy, precision, and scalp pressure uniformity of the sensors are evaluated. A root-mean-squared (RMS) positional precision of 0.89 ± 0.23 mm, percent arc subdivision RMS accuracy of 0.19 ± 0.15%, and mean normal force on the scalp of 2.28 ± 0.88 N at 5 mm displacement were found. Geometric measurements indicate that the probe will accommodate the full range of adult head sizes. The placement accuracy, precision, and uniformity of sensor contact pressure of the proposed head probe are important determinants of data quality in noninvasive brain monitoring with simultaneous NIRS-EEG.

Dynamic seizure imaging in patients with extratemporal lobe epilepsy

Epilepsy is a common neurological disease that affects about 50 million people worldwide. Extratemporal lobe epilepsy, which represents an important type of epilepsy, may involve seizure activity in various lobes and the surgical treatment in these patients tends to have less favorable surgical outcome. Noninvasive seizure imaging in drug-resistant patients is of vital importance to image the seizure onset zones (SOZs) and understand the mechanisms for an improved treatment plan. In this study, we directly imaged the seizure sources in 8 extratemporal lobe partial epilepsy patients from noninvasive EEG. The surgically resected regions and SOZs identified from intracranial EEG (iEEG) recordings were used to evaluate the source imaging results. All of the eight patients underwent resective surgery and the estimated seizure sources were co-located with the resection zone. Seven of the patients had iEEG recordings available and the source imaging results were concordant with the SOZs marked on the intracranial recording grid. The present results suggest that dynamic seizure imaging could be potentially useful to image the SOZs of extratemporal lobe seizures and help the pre-surgical planning of epilepsy patients.

Identification of feedback loops in neural networks based on multi-step Granger causality

MOTIVATION

Feedback circuits are crucial network motifs, ubiquitously found in many intra- and inter-cellular regulatory networks, and also act as basic building blocks for inducing synchronized bursting behaviors in neural network dynamics. Therefore, the system-level identification of feedback circuits using time-series measurements is critical to understand the underlying regulatory mechanism of synchronized bursting behaviors.

RESULTS

Multi-Step Granger Causality Method (MSGCM) was developed to identify feedback loops embedded in biological networks using time-series experimental measurements. Based on multivariate time-series analysis, MSGCM used a modified Wald test to infer the existence of multi-step Granger causality between a pair of network nodes. A significant bi-directional multi-step Granger causality between two nodes indicated the existence of a feedback loop. This new identification method resolved the drawback of the previous non-causal impulse response component method which was only applicable to networks containing no co-regulatory forward path. MSGCM also significantly improved the ratio of correct identification of feedback loops. In this study, the MSGCM was testified using synthetic pulsed neural network models and also in vitro cultured rat neural networks using multi-electrode array. As a result, we found a large number of feedback loops in the in vitro cultured neural networks with apparent synchronized oscillation, indicating a close relationship between synchronized oscillatory bursting behavior and underlying feedback loops. The MSGCM is an efficient method to investigate feedback loops embedded in in vitro cultured neural networks. The identified feedback loop motifs are considered as an important design principle responsible for the synchronized bursting behavior in neural networks.

Dynamic imaging of seizure activity in pediatric epilepsy patients

OBJECTIVE

To investigate the feasibility of using noninvasive EEG source imaging approach to image continuous seizure activity in pediatric epilepsy patients.

METHODS

Nine pediatric patients with medically intractable epilepsy were included in this study. Eight of the patients had extratemporal lobe epilepsy and one had temporal lobe epilepsy. All of the patients underwent resective surgery and seven of them underwent intracranial EEG (iEEG) monitoring. The ictal EEG was analyzed using a noninvasive dynamic seizure imaging (DSI) approach. The DSI approach separates scalp EEGs into independent components and extracts the spatio-temporal ictal features to achieve dynamic imaging of seizure sources. Surgical resection and intracranial recordings were used to validate the noninvasive imaging results.

RESULTS

The DSI determined seizure onset zones (SOZs) in these patients were localized within or in close vicinity to the surgically resected region. In the seven patients with intracranial monitoring, the estimated seizure onset sources were concordant with the seizure onset zones of iEEG. The DSI also localized the multiple foci involved in the later seizure propagation, which were confirmed by the iEEG recordings.

CONCLUSIONS

Dynamic seizure imaging can noninvasively image the seizure activations in pediatric patients with both temporal and extratemporal lobe epilepsy.

SIGNIFICANCE

EEG seizure imaging can potentially be used to noninvasively image the SOZs and aid the pre-surgical planning in pediatric epilepsy patients.

Inferring functional neural connectivity with phase synchronization analysis: a review of methodology

Functional neural connectivity is drawing increasing attention in neuroscience research. To infer functional connectivity from observed neural signals, various methods have been proposed. Among them, phase synchronization analysis is an important and effective one which examines the relationship of instantaneous phase between neural signals but neglecting the influence of their amplitudes. In this paper, we review the advances in methodologies of phase synchronization analysis. In particular, we discuss the definitions of instantaneous phase, the indexes of phase synchronization and their significance test, the issues that may affect the detection of phase synchronization and the extensions of phase synchronization analysis. In practice, phase synchronization analysis may be affected by observational noise, insufficient samples of the signals, volume conduction, and reference in recording neural signals. We make comments and suggestions on these issues so as to better apply phase synchronization analysis to inferring functional connectivity from neural signals.

Cortical networks of hemianopia stroke patients: a graph theoretical analysis of EEG signals at resting state

Visual cortical stroke patients may have hemianopia symptom, which affects a number of visual functions. Most studies on hemianopia stroke have mainly focused on cortical activation during visual stimulation, leaving the pattern of functional connectivity between different brain regions uncovered yet. In the present study, we investigate the resting neural networks of hemianopia stroke patients by graph theoretical analysis of functional brain networks constructed with phase synchronization indexes of multichannel electroencephalography (EEG) signals. Our results showed that although the global network topological metrics, i.e., weighted clustering coefficient and characteristic path length of patients and healthy controls are comparable, the left primary visual cortex of patients tend to be less active than that of age-matched healthy subjects. However, hemianopia patients showed greater activation in the ipsilesional (left) temporopolar and orbit frontal areas and the contralesional (right) associative visual cortex. These results may offer new insight into neural substrates of the hemianopia stroke, and the further study of neural plasticity and brain reorganization after hemianopia.

Sparse electromagnetic source imaging using EEG and MEG

The present study proposed the combined use of EEG and MEG data in a new sparse electromagnetic source imaging (ESI) technique, i.e., variation-based sparse cortical current density (VB-SCCD) method. Monte Carlo simulations were conducted to investigate the performance of the proposed approach in multiple extended brain activations (up to ten) that were randomly generated. Experimental EEG and MEG data from a face recognition task were further used to evaluate the performance of VB-SCCD. The present results indicate that the proposed approach can accurately reconstruct multiple brain activations and their spatial extents. The source imaging results from real data further demonstrate it is capable to recover networked brain activations involving multiple cortical regions, which are consistent with results from functional magnetic resonance imaging in same task paradigm. The present results further indicate the capability of the proposed approach in reconstructing deep brain sources and temporal dynamics of brain sources at millisecond resolutions. It thus suggests that sparse ESI using combined EEG and MEG is a promising technique probing detailed spatiotemporal brain activations.

Seizure source imaging by means of FINE spatio-temporal dipole localization and directed transfer function in partial epilepsy patients

OBJECTIVE

To investigate the usage of a high-density EEG recording system and source imaging technique for localizing seizure activity in patients with medically intractable partial epilepsy.

METHODS

High-density, 76-channel scalp EEG signals were recorded in 10 patients with partial epilepsy. The patients underwent routine clinical pre-surgical evaluation and all had resective surgery with seizure free outcome. After applying a FINE (first principle vectors) spatio-temporal source localization and DTF (directed transfer function) connectivity analysis approach, ictal sources were imaged. Effects of number of scalp EEG electrodes on the seizure localization were also assessed using 76, 64, 48, 32, and 21 electrodes, respectively.

RESULTS

Surgical resections were used to assess the source imaging results. Results from the 76-channel EEG in the 10 patients showed high correlation with the surgically resected brain regions. The localization of seizure onset zone from 76-channel EEG showed improved source detection accuracy compared to other EEG configurations with fewer electrodes.

CONCLUSIONS

FINE together with DTF was able to localize seizure onset zones of partial epilepsy patients. High-density EEG recording can help achieve improved seizure source imaging.

SIGNIFICANCE

The present results suggest the promise of high-density EEG and electrical source imaging for noninvasively localizing seizure onset zones.

Source connectivity analysis from MEG and its application to epilepsy source localization

We report an approach to perform source connectivity analysis from MEG, and initially evaluate this approach to interictal MEG to localize epileptogenic foci and analyze interictal discharge propagations in patients with medically intractable epilepsy. Cortical activities were reconstructed from MEG using individual realistic geometry boundary element method head models. Directional connectivity among cortical regions of interest was then estimated using directed transfer function. The MEG source connectivity analysis method was implemented in the eConnectome software, which is open-source and freely available at http://econnectome.umn.edu . As an initial evaluation, the method was applied to study MEG interictal spikes from five epilepsy patients. Estimated primary epileptiform sources were consistent with surgically resected regions, suggesting the feasibility of using cortical source connectivity analysis from interictal MEG for potential localization of epileptiform activities.