Electrophysiological imaging of brain activity and connectivity-challenges and opportunities

Deactivation in the posterior mid-cingulate cortex reflects perceptual transitions during binocular rivalry: Evidence from simultaneous EEG-fMRI

Binocular rivalry is a phenomenon in which perception spontaneously shifts between two different images that are dichoptically presented to the viewer. By elucidating the cortical networks responsible for these stochastic fluctuations in perception, we can potentially learn much about the neural correlates of visual awareness. We obtained concurrent EEG-fMRI data for a group of 20 healthy human subjects during the continuous presentation of dichoptic visual stimuli. The two eyes’ images were tagged with different temporal frequencies so that eye specific steady-state visual evoked potential (SSVEP) signals could be extracted from the EEG data for direct comparison with changes in fMRI BOLD activity associated with binocular rivalry. We additionally included a smooth replay condition that emulated the perceptual transitions experienced during binocular rivalry as a control stimulus. We evaluated a novel SSVEP-informed fMRI analysis in this study in order to delineate the temporal dynamics of rivalry-related BOLD activity from both an electrophysiological and behavioral perspective. In this manner, we assessed BOLD activity during rivalry that was directly correlated with peaks and crosses of the two rivaling, frequency-tagged SSVEP signals, for comparison with BOLD activity associated with subject reported perceptual transitions. Our findings point to a critical role of a right lateralized fronto-parietal network in the processing of bistable stimuli, given that BOLD activity in the right superior/inferior parietal lobules was significantly elevated throughout binocular rivalry and in particular during perceptual transitions, compared with the replay condition. Based on the SSVEP-informed analysis, rivalry was further associated with significantly enhanced BOLD suppression in the posterior mid-cingulate cortex during perceptual transitions, compared with SSVEP crosses. Overall, this work points to a careful interplay between early visual areas, the right posterior parietal cortex and the mid-cingulate cortex in mediating the spontaneous perceptual changes associated with binocular rivalry and has significant implications for future multimodal imaging studies of perception and awareness.

Reference layer adaptive filtering (RLAF) for EEG artifact reduction in simultaneous EEG-fMRI

OBJECTIVE

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) combines advantages of both methods, namely high temporal resolution of EEG and high spatial resolution of fMRI. However, EEG quality is limited due to severe artifacts caused by fMRI scanners.

APPROACH

To improve EEG data quality substantially, we introduce methods that use a reusable reference layer EEG cap prototype in combination with adaptive filtering. The first method, reference layer adaptive filtering (RLAF), uses adaptive filtering with reference layer artifact data to optimize artifact subtraction from EEG. In the second method, multi band reference layer adaptive filtering (MBRLAF), adaptive filtering is performed on bandwidth limited sub-bands of the EEG and the reference channels.

MAIN RESULTS

The results suggests that RLAF outperforms the baseline method, average artifact subtraction, in all settings and also its direct predecessor, reference layer artifact subtraction (RLAS), in lower (<35 Hz) frequency ranges. MBRLAF is computationally more demanding than RLAF, but highly effective in all EEG frequency ranges. Effectivity is determined by visual inspection, as well as root-mean-square voltage reduction and power reduction of EEG provided that physiological EEG components such as occipital EEG alpha power and visual evoked potentials (VEP) are preserved. We demonstrate that both, RLAF and MBRLAF, improve VEP quality. For that, we calculate the mean-squared-distance of single trial VEP to the mean VEP and estimate single trial VEP classification accuracies. We found that the average mean-squared-distance is lowest and the average classification accuracy is highest after MBLAF. RLAF was second best.

SIGNIFICANCE

In conclusion, the results suggests that RLAF and MBRLAF are potentially very effective in improving EEG quality of simultaneous EEG-fMRI. Highlights We present a new and reusable reference layer cap prototype for simultaneous EEG-fMRI We introduce new algorithms for reducing EEG artifacts due to simultaneous fMRI The algorithms combine a reference layer and adaptive filtering Several evaluation criteria suggest superior effectivity in terms of artifact reduction We demonstrate that physiological EEG components are preserved.

A novel method linking neural connectivity to behavioral fluctuations: Behavior-regressed connectivity

BACKGROUND

During an experimental session, behavioral performance fluctuates, yet most neuroimaging analyses of functional connectivity derive a single connectivity pattern. These conventional connectivity approaches assume that since the underlying behavior of the task remains constant, the connectivity pattern is also constant.

NEW METHOD

We introduce a novel method, behavior-regressed connectivity (BRC), to directly examine behavioral fluctuations within an experimental session and capture their relationship to changes in functional connectivity. This method employs the weighted phase lag index (WPLI) applied to a window of trials with a weighting function. Using two datasets, the BRC results are compared to conventional connectivity results during two time windows: the one second before stimulus onset to identify predictive relationships, and the one second after onset to capture task-dependent relationships.

RESULTS

In both tasks, we replicate the expected results for the conventional connectivity analysis, and extend our understanding of the brain-behavior relationship using the BRC analysis, demonstrating subject-specific BRC maps that correspond to both positive and negative relationships with behavior. Comparison with Existing Method(s): Conventional connectivity analyses assume a consistent relationship between behaviors and functional connectivity, but the BRC method examines performance variability within an experimental session to understand dynamic connectivity and transient behavior.

CONCLUSION

The BRC approach examines connectivity as it covaries with behavior to complement the knowledge of underlying neural activity derived from conventional connectivity analyses. Within this framework, BRC may be implemented for the purpose of understanding performance variability both within and between participants.

Characterization of functional brain activity and connectivity using EEG and fMRI in patients with sickle cell disease

Sickle cell disease (SCD) is a red blood cell disorder that causes many complications including life-long pain. Treatment of pain remains challenging due to a poor understanding of the mechanisms and limitations to characterize and quantify pain. In the present study, we examined simultaneously recording functional MRI (fMRI) and electroencephalogram (EEG) to better understand neural connectivity as a consequence of chronic pain in SCD patients. We performed independent component analysis and seed-based connectivity on fMRI data. Spontaneous power and microstate analysis was performed on EEG-fMRI data. ICA analysis showed that patients lacked activity in the default mode network (DMN) and executive control network compared to controls. EEG-fMRI data revealed that the insula cortex's role in salience increases with age in patients. EEG microstate analysis showed patients had increased activity in pain processing regions. The cerebellum in patients showed a stronger connection to the periaqueductal gray matter (involved in pain inhibition), and negative connections to pain processing areas. These results suggest that patients have reduced activity of DMN and increased activity in pain processing regions during rest. The present findings suggest resting state connectivity differences between patients and controls can be used as novel biomarkers of SCD pain.

•

Simultaneous EEG-fMRI recordings revealed altered connectivity in sickle cell patients.

•

Reduced activity observed in default mode network and executive control network.

•

Patients' salience network strength increases with age; opposite seen in controls.

•

EEG-fMRI parameters reflect disease severity in sickle cell patients.

Apolipoprotein E4 causes early olfactory network abnormalities and short-term olfactory memory impairments

While apolipoprotein (Apo) E4 is linked to increased incidence of Alzheimer's disease (AD), there is growing evidence that it plays a role in functional brain irregularities that are independent of AD pathology. However, ApoE4-driven functional differences within olfactory processing regions have yet to be examined. Utilizing knock-in mice humanized to ApoE4 versus the more common ApoE3, we examined a simple olfactory perceptual memory that relies on the transfer of information from the olfactory bulb (OB) to the piriform cortex (PCX), the primary cortical region involved in higher order olfaction. In addition, we have recorded in vivo resting and odor-evoked local field potentials (LPF) from both brain regions and measured corresponding odor response magnitudes in anesthetized young (6-month-old) and middle-aged (12-month-old) ApoE mice. Young ApoE4 compared to ApoE3 mice exhibited a behavioral olfactory deficit coinciding with hyperactive odor-evoked response magnitudes within the OB that were not observed in older ApoE4 mice. Meanwhile, middle-aged ApoE4 compared to ApoE3 mice exhibited heightened response magnitudes in the PCX without a corresponding olfactory deficit, suggesting a shift with aging in ApoE4-driven effects from OB to PCX. Interestingly, the increased ApoE4-specific response in the PCX at middle-age was primarily due to a dampening of baseline spontaneous activity rather than an increase in evoked response power. Our findings indicate that early ApoE4-driven olfactory memory impairments and OB network abnormalities may be a precursor to later network dysfunction in the PCX, a region that not only is targeted early in AD, but may be selectively vulnerable to ApoE4 genotype.

Combined EMD-sLORETA Analysis of EEG Data Collected during a Contour Integration Task

Lately, Ensemble Empirical Mode Decomposition (EEMD) techniques receive growing interest in biomedical data analysis. Event-Related Modes (ERMs) represent features extracted by an EEMD from electroencephalographic (EEG) recordings. We present a new approach for source localization of EEG data based on combining ERMs with inverse models. As the first step, 64 channel EEG recordings are pooled according to six brain areas and decomposed, by applying an EEMD, into their underlying ERMs. Then, based upon the problem at hand, the most closely related ERM, in terms of frequency and amplitude, is combined with inverse modeling techniques for source localization. More specifically, the standardized low resolution brain electromagnetic tomography (sLORETA) procedure is employed in this work. Accuracy and robustness of the results indicate that this approach deems highly promising in source localization techniques for EEG data.

Noninvasive Electromagnetic Source Imaging and Granger Causality Analysis: An Electrophysiological Connectome (eConnectome) Approach

OBJECTIVE

Combined source-imaging techniques and directional connectivity analysis can provide useful information about the underlying brain networks in a noninvasive fashion. Source-imaging techniques have been used successfully to either determine the source of activity or to extract source time-courses for Granger causality analysis, previously. In this work, we utilize source-imaging algorithms to both find the network nodes [regions of interest (ROI)] and then extract the activation time series for further Granger causality analysis. The aim of this work is to find network nodes objectively from noninvasive electromagnetic signals, extract activation time-courses, and apply Granger analysis on the extracted series to study brain networks under realistic conditions.

METHODS

Source-imaging methods are used to identify network nodes and extract time-courses and then Granger causality analysis is applied to delineate the directional functional connectivity of underlying brain networks. Computer simulations studies where the underlying network (nodes and connectivity pattern) is known were performed; additionally, this approach has been evaluated in partial epilepsy patients to study epilepsy networks from interictal and ictal signals recorded by EEG and/or Magnetoencephalography (MEG).

RESULTS

Localization errors of network nodes are less than 5 mm and normalized connectivity errors of ∼20% in estimating underlying brain networks in simulation studies. Additionally, two focal epilepsy patients were studied and the identified nodes driving the epileptic network were concordant with clinical findings from intracranial recordings or surgical resection.

CONCLUSION

Our study indicates that combined source-imaging algorithms with Granger causality analysis can identify underlying networks precisely (both in terms of network nodes location and internodal connectivity).

SIGNIFICANCE

The combined source imaging and Granger analysis technique is an effective tool for studying normal or pathological brain conditions.

Imaging brain source extent from EEG/MEG by means of an iteratively reweighted edge sparsity minimization (IRES) strategy

Estimating extended brain sources using EEG/MEG source imaging techniques is challenging. EEG and MEG have excellent temporal resolution at millisecond scale but their spatial resolution is limited due to the volume conduction effect. We have exploited sparse signal processing techniques in this study to impose sparsity on the underlying source and its transformation in other domains (mathematical domains, like spatial gradient). Using an iterative reweighting strategy to penalize locations that are less likely to contain any source, it is shown that the proposed iteratively reweighted edge sparsity minimization (IRES) strategy can provide reasonable information regarding the location and extent of the underlying sources. This approach is unique in the sense that it estimates extended sources without the need of subjectively thresholding the solution. The performance of IRES was evaluated in a series of computer simulations. Different parameters such as source location and signal-to-noise ratio were varied and the estimated results were compared to the targets using metrics such as localization error (LE), area under curve (AUC) and overlap between the estimated and simulated sources. It is shown that IRES provides extended solutions which not only localize the source but also provide estimation for the source extent. The performance of IRES was further tested in epileptic patients undergoing intracranial EEG (iEEG) recording for pre-surgical evaluation. IRES was applied to scalp EEGs during interictal spikes, and results were compared with iEEG and surgical resection outcome in the patients. The pilot clinical study results are promising and demonstrate a good concordance between noninvasive IRES source estimation with iEEG and surgical resection outcomes in the same patients. The proposed algorithm, i.e. IRES, estimates extended source solutions from scalp electromagnetic signals which provide relatively accurate information about the location and extent of the underlying source.

The Power of Connecting Dots: Advanced Techniques to Evaluate Brain Functional Connectivity in Humans

Brain functional connectivity estimation allows us to depict patterns of cerebral activity not understandable otherwise with the standard brain imaging techniques such as functional magnetic resonance imaging (fMRI) as well as electro or magnetoencephalography (hr-EEG, MEG). This special issue of the IEEE Transactions on Biomedical Engineering reports a range of methodological innovations toward the estimation of functional connectivity from brain activity data, with emphasis on neuroelectric and hemodynamic imaging modalities. Functional connectivity methodologies enable "connecting of the dots" derived from brain activity observations over multiple distributed sites, as depicted by such fMRI and hr-EEG/MEG devices.

Electrophysiological Source Imaging of Brain Networks Perturbed by Low-Intensity Transcranial Focused Ultrasound

OBJECTIVE

Transcranial focused ultrasound (tFUS) has been introduced as a noninvasive neuromodulation technique with good spatial selectivity. We report an experimental investigation to detect noninvasive electrophysiological response induced by low-intensity tFUS in an in vivo animal model and perform electrophysiological source imaging (ESI) of tFUS-induced brain activity from noninvasive scalp EEG recordings.

METHODS

A single-element ultrasound transducer was used to generate low-intensity tFUS ( ) and induce brain activation at multiple selected sites in an in vivo rat model. Up to 16 scalp electrodes were used to record tFUS-induced EEG. Event-related potentials were analyzed in time, frequency, and spatial domains. Current source distributions were estimated by ESI to reconstruct spatiotemporal distributions of brain activation induced by tFUS.

RESULTS

Neuronal activation was observed following low-intensity tFUS, as correlated to tFUS intensity and sonication duration. ESI revealed initial focal activation in cortical area corresponding to tFUS stimulation site and the activation propagating to surrounding areas over time.

CONCLUSION

The present results demonstrate the feasibility of noninvasively recording brain electrophysiological response in vivo following low-intensity tFUS stimulation, and the feasibility of imaging spatiotemporal distributions of brain activation as induced by tFUS in vivo.

SIGNIFICANCE

The present approach may lead to a new means of imaging brain activity using tFUS perturbation and a closed-loop ESI-guided tFUS neuromodulation modality.

Bayesian electromagnetic spatio-temporal imaging of extended sources with Markov Random Field and temporal basis expansion

Estimating the locations and spatial extents of brain sources poses a long-standing challenge for electroencephalography and magnetoencephalography (E/MEG) source imaging. In the present work, a novel source imaging method, Bayesian Electromagnetic Spatio-Temporal Imaging of Extended Sources (BESTIES), which is built upon a Bayesian framework that determines the spatio-temporal smoothness of source activities in a fully data-driven fashion, is proposed to address this challenge. In particular, a Markov Random Field (MRF), which can precisely capture local cortical interactions, is employed to characterize the spatial smoothness of source activities, the temporal dynamics of which are modeled by a set of temporal basis functions (TBFs). Crucially, all of the unknowns in the MRF and TBF models are learned from the data. To accomplish model inference efficiently on high-resolution source spaces, a scalable algorithm is developed to approximate the posterior distribution of the source activities, which is based on the variational Bayesian inference and convex analysis. The performance of BESTIES is assessed using both simulated and actual human E/MEG data. Compared with L₂-norm constrained methods, BESTIES is superior in reconstructing extended sources with less spatial diffusion and less localization error. By virtue of the MRF, BESTIES also overcomes the drawback of over-focal estimates in sparse constrained methods.

In Vivo Tumour Mapping Using Electrocorticography Alterations During Awake Brain Surgery: A Pilot Study

During awake brain surgery for tumour resection, in situ EEG recording (ECoG) is used to identify eloquent areas surrounding the tumour. We used the ECoG setup to record the electrical activity of cortical and subcortical tumours and then performed frequency and connectivity analyses in order to identify ECoG impairments and map tumours. We selected 16 patients with cortical (8) and subcortical (8) tumours undergoing awake brain surgery. For each patient, we computed the spectral content of tumoural and healthy areas in each frequency band. We computed connectivity of each electrode using connectivity markers (linear and non-linear correlations, phase-locking and coherence). We performed comparisons between healthy and tumour electrodes. The ECoG alterations were used to implement automated classification of the electrodes using clustering or neural network algorithms. ECoG alterations were used to image cortical tumours.Cortical tumours were found to profoundly alter all frequency contents (normalized and absolute power), with an increase in the δ activity and a decreases for the other bands (P < 0.05). Cortical tumour electrodes showed high level of connectivity compared to surrounding electrodes (all markers, P < 0.05). For subcortical tumours, a relative decrease in the γ1 band and in the alpha band in absolute amplitude (P < 0.05) were the only abnormalities. The neural network algorithm classification had a good performance: 93.6 % of the electrodes were classified adequately on a test subject. We found significant spectral and connectivity ECoG changes for cortical tumours, which allowed tumour recognition. Artificial neural algorithm pattern recognition seems promising for electrode classification in awake tumour surgery.

From Structure to Circuits: The Contribution of MEG Connectivity Studies to Functional Neurosurgery

New advances in structural neuroimaging have revealed the intricate and extensive connections within the brain, data which have informed a number of ambitious projects such as the mapping of the human connectome. Elucidation of the structural connections of the brain, at both the macro and micro levels, promises new perspectives on brain structure and function that could translate into improved outcomes in functional neurosurgery. The understanding of neuronal structural connectivity afforded by these data now offers a vista on the brain, in both healthy and diseased states, that could not be seen with traditional neuroimaging. Concurrent with these developments in structural imaging, a complementary modality called magnetoencephalography (MEG) has been garnering great attention because it too holds promise for being able to shed light on the intricacies of functional brain connectivity. MEG is based upon the elemental principle of physics that an electrical current generates a magnetic field. Hence, MEG uses highly sensitive biomagnetometers to measure extracranial magnetic fields produced by intracellular neuronal currents. Put simply then, MEG is a measure of neurophysiological activity, which captures the magnetic fields generated by synchronized intraneuronal electrical activity. As such, MEG recordings offer exquisite resolution in the time and oscillatory domain and, as well, when co-registered with magnetic resonance imaging (MRI), offer excellent resolution in the spatial domain. Recent advances in MEG computational and graph theoretical methods have led to studies of connectivity in the time-frequency domain. As such, MEG can elucidate a neurophysiological-based functional circuitry that may enhance what is seen with MRI connectivity studies. In particular, MEG may offer additional insight not possible by MRI when used to study complex eloquent function, where the precise timing and coordination of brain areas is critical. This article will review the traditional use of MEG for functional neurosurgery, describe recent advances in MEG connectivity analyses, and consider the additional benefits that could be gained with the inclusion of MEG connectivity studies. Since MEG has been most widely applied to the study of epilepsy, we will frame this article within the context of epilepsy surgery and functional neurosurgery for epilepsy.

A Unified Estimation Framework for State-Related Changes in Effective Brain Connectivity

OBJECTIVE

This paper addresses the critical problem of estimating time-evolving effective brain connectivity. Current approaches based on sliding window analysis or time-varying coefficient models do not simultaneously capture both slow and abrupt changes in causal interactions between different brain regions.

METHODS

To overcome these limitations, we develop a unified framework based on a switching vector autoregressive (SVAR) model. Here, the dynamic connectivity regimes are uniquely characterized by distinct vector autoregressive (VAR) processes and allowed to switch between quasi-stationary brain states. The state evolution and the associated directed dependencies are defined by a Markov process and the SVAR parameters. We develop a three-stage estimation algorithm for the SVAR model: 1) feature extraction using time-varying VAR (TV-VAR) coefficients, 2) preliminary regime identification via clustering of the TV-VAR coefficients, 3) refined regime segmentation by Kalman smoothing and parameter estimation via expectation-maximization algorithm under a state-space formulation, using initial estimates from the previous two stages.

RESULTS

The proposed framework is adaptive to state-related changes and gives reliable estimates of effective connectivity. Simulation results show that our method provides accurate regime change-point detection and connectivity estimates. In real applications to brain signals, the approach was able to capture directed connectivity state changes in functional magnetic resonance imaging data linked with changes in stimulus conditions, and in epileptic electroencephalograms, differentiating ictal from nonictal periods.

CONCLUSION

The proposed framework accurately identifies state-dependent changes in brain network and provides estimates of connectivity strength and directionality.

SIGNIFICANCE

The proposed approach is useful in neuroscience studies that investigate the dynamics of underlying brain states.

Spectral and spatial changes of brain rhythmic activity in response to the sustained thermal pain stimulation

The aim of this study was to investigate the neurophysiological correlates of pain caused by sustained thermal stimulation. A group of 21 healthy volunteers was studied. Sixty-four channel continuous electroencephalography (EEG) was recorded while the subject received tonic thermal stimulation. Spectral changes extracted from EEG were quantified and correlated with pain scales reported by subjects, the stimulation intensity, and the time course. Network connectivity was assessed to study the changes in connectivity patterns and strengths among brain regions that have been previously implicated in pain processing. Spectrally, a global reduction in power was observed in the lower spectral range, from delta to alpha, with the most marked changes in the alpha band. Spatially, the contralateral region of the somatosensory cortex, identified using source localization, was most responsive to stimulation status. Maximal desynchrony was observed when stimulation was present. The degree of alpha power reduction was linearly correlated to the pain rating reported by the subjects. Contralateral alpha power changes appeared to be a robust correlate of pain intensity experienced by the subjects. Granger causality analysis showed changes in network level connectivity among pain-related brain regions due to high intensity of pain stimulation versus innocuous warm stimulation. These results imply the possibility of using noninvasive EEG to predict pain intensity and to study the underlying pain processing mechanism in coping with prolonged painful experiences. Once validated in a broader population, the present EEG-based approach may provide an objective measure for better pain management in clinical applications. Hum Brain Mapp 37:2976-2991, 2016. © 2016 Wiley Periodicals, Inc.

Complex discharge-affecting networks in juvenile myoclonic epilepsy: A simultaneous EEG-fMRI study

Juvenile myoclonic epilepsy (JME) is a common subtype of idiopathic generalized epilepsies (IGEs) and is characterized by myoclonic jerks, tonic-clonic seizures and infrequent absence seizures. The network notion has been proposed to better characterize epilepsy. However, many issues remain not fully understood in JME, such as the associations between discharge-affecting networks and the relationships among resting-state networks. In this project, eigenspace maximal information canonical correlation analysis (emiCCA) and functional network connectivity (FNC) analysis were applied to simultaneous EEG-fMRI data from JME patients. The main findings of our study are as follows: discharge-affecting networks comprising the default model (DMN), self-reference (SRN), basal ganglia (BGN) and frontal networks have linear and nonlinear relationships with epileptic discharge information in JME patients; the DMN, SRN and BGN have dense/specific associations with discharge-affecting networks as well as resting-state networks; and compared with controls, significantly increased FNCs between the salience network (SN) and resting-state networks are found in JME patients. These findings suggest that the BGN, DMN and SRN may play intermediary roles in the modulation and propagation of epileptic discharges. These roles further tend to disturb the switching function of the SN in JME patients. We also postulate that emiCCA and FNC analysis may provide a potential analysis platform to provide insights into our understanding of the pathophysiological mechanism of epilepsy subtypes such as JME. Hum Brain Mapp 37:3515-3529, 2016. © 2016 Wiley Periodicals, Inc.

Hybrid EEG-fNIRS Asynchronous Brain-Computer Interface for Multiple Motor Tasks

Non-invasive Brain-Computer Interfaces (BCI) have demonstrated great promise for neuroprosthetics and assistive devices. Here we aim to investigate methods to combine Electroencephalography (EEG) and functional Near-Infrared Spectroscopy (fNIRS) in an asynchronous Sensory Motor rhythm (SMR)-based BCI. We attempted to classify 4 different executed movements, namely, Right-Arm—Left-Arm—Right-Hand—Left-Hand tasks. Previous studies demonstrated the benefit of EEG-fNIRS combination. However, since normally fNIRS hemodynamic response shows a long delay, we investigated new features, involving slope indicators, in order to immediately detect changes in the signals. Moreover, Common Spatial Patterns (CSPs) have been applied to both EEG and fNIRS signals. 15 healthy subjects took part in the experiments and since 25 trials per class were available, CSPs have been regularized with information from the entire population of participants and optimized using genetic algorithms. The different features have been compared in terms of performance and the dynamic accuracy over trials shows that the introduced methods diminish the fNIRS delay in the detection of changes.

EEG Source Imaging Enhances the Decoding of Complex Right-Hand Motor Imagery Tasks

GOAL

Sensorimotor-based brain-computer interfaces (BCIs) have achieved successful control of real and virtual devices in up to three dimensions; however, the traditional sensor-based paradigm limits the intuitive use of these systems. Many control signals for state-of-the-art BCIs involve imagining the movement of body parts that have little to do with the output command, revealing a cognitive disconnection between the user's intent and the action of the end effector. Therefore, there is a need to develop techniques that can identify with high spatial resolution the self-modulated neural activity reflective of the actions of a helpful output device.

METHODS

We extend previous EEG source imaging (ESI) work to decoding natural hand/wrist manipulations by applying a novel technique to classifying four complex motor imaginations of the right hand: flexion, extension, supination, and pronation.

RESULTS

We report an increase of up to 18.6% for individual task classification and 12.7% for overall classification using the proposed ESI approach over the traditional sensor-based method.

CONCLUSION

ESI is able to enhance BCI performance of decoding complex right-hand motor imagery tasks.

SIGNIFICANCE

This study may lead to the development of BCI systems with naturalistic and intuitive motor imaginations, thus facilitating broad use of noninvasive BCIs.

A Space-Time-Frequency Dictionary for Sparse Cortical Source Localization

OBJECTIVE

Cortical source imaging aims at identifying activated cortical areas on the surface of the cortex from the raw electroencephalogram (EEG) data. This problem is ill posed, the number of channels being very low compared to the number of possible source positions.

METHODS

In some realistic physiological situations, the active areas are sparse in space and of short time durations, and the amount of spatio-temporal data to carry the inversion is then limited. In this study, we propose an original data driven space-time-frequency (STF) dictionary which takes into account simultaneously both spatial and time-frequency sparseness while preserving smoothness in the time frequency (i.e., nonstationary smooth time courses in sparse locations). Based on these assumptions, we take benefit of the matching pursuit (MP) framework for selecting the most relevant atoms in this highly redundant dictionary.

RESULTS

We apply two recent MP algorithms, single best replacement (SBR) and source deflated matching pursuit, and we compare the results using a spatial dictionary and the proposed STF dictionary to demonstrate the improvements of our multidimensional approach. We also provide comparison using well-established inversion methods, FOCUSS and RAP-MUSIC, analyzing performances under different degrees of nonstationarity and signal to noise ratio.

CONCLUSION

Our STF dictionary combined with the SBR approach provides robust performances on realistic simulations. From a computational point of view, the algorithm is embedded in the wavelet domain, ensuring high efficiency in term of computation time.

SIGNIFICANCE

The proposed approach ensures fast and accurate sparse cortical localizations on highly nonstationary and noisy data.

Long-Range Reduced Predictive Information Transfers of Autistic Youths in EEG Sensor-Space During Face Processing

The majority of previous functional/effective connectivity studies conducted on the autistic patients converged to the underconnectivity theory of ASD: "long-range underconnectivity and sometimes short-rang overconnectivity". However, to the best of our knowledge the total (linear and nonlinear) predictive information transfers (PITs) of autistic patients have not been investigated yet. Also, EEG data have rarely been used for exploring the information processing deficits in autistic subjects. This study is aimed at comparing the total (linear and nonlinear) PITs of autistic and typically developing healthy youths during human face processing by using EEG data. The ERPs of 12 autistic youths and 19 age-matched healthy control (HC) subjects were recorded while they were watching upright and inverted human face images. The PITs among EEG channels were quantified using two measures separately: transfer entropy with self-prediction optimality (TESPO), and modified transfer entropy with self-prediction optimality (MTESPO). Afterwards, the directed differential connectivity graphs (dDCGs) were constructed to characterize the significant changes in the estimated PITs of autistic subjects compared with HC ones. By using both TESPO and MTESPO, long-range reduction of PITs of ASD group during face processing was revealed (particularly from frontal channels to right temporal channels). Also, it seemed the orientation of face images (upright or upside down) did not modulate the binary pattern of PIT-based dDCGs, significantly. Moreover, compared with TESPO, the results of MTESPO were more compatible with the underconnectivity theory of ASD in the sense that MTESPO showed no long-range increase in PIT. It is also noteworthy that to the best of our knowledge it is the first time that a version of MTE is applied for patients (here ASD) and it is also its first use for EEG data analysis.

Real-Time Neuroimaging and Cognitive Monitoring Using Wearable Dry EEG

GOAL

We present and evaluate a wearable high-density dry-electrode EEG system and an open-source software framework for online neuroimaging and state classification.

METHODS

The system integrates a 64-channel dry EEG form factor with wireless data streaming for online analysis. A real-time software framework is applied, including adaptive artifact rejection, cortical source localization, multivariate effective connectivity inference, data visualization, and cognitive state classification from connectivity features using a constrained logistic regression approach (ProxConn). We evaluate the system identification methods on simulated 64-channel EEG data. Then, we evaluate system performance, using ProxConn and a benchmark ERP method, in classifying response errors in nine subjects using the dry EEG system.

RESULTS

Simulations yielded high accuracy (AUC = 0.97 ± 0.021) for real-time cortical connectivity estimation. Response error classification using cortical effective connectivity [short-time direct-directed transfer function (sdDTF)] was significantly above chance with similar performance (AUC) for cLORETA (0.74 ±0.09) and LCMV (0.72 ±0.08) source localization. Cortical ERP-based classification was equivalent to ProxConn for cLORETA (0.74 ±0.16) but significantly better for LCMV (0.82 ±0.12) .

CONCLUSION

We demonstrated the feasibility for real-time cortical connectivity analysis and cognitive state classification from high-density wearable dry EEG.

SIGNIFICANCE

This paper is the first validated application of these methods to 64-channel dry EEG. This study addresses a need for robust real-time measurement and interpretation of complex brain activity in the dynamic environment of the wearable setting. Such advances can have broad impact in research, medicine, and brain-computer interfaces. The pipelines are made freely available in the open-source SIFT and BCILAB toolboxes.

EEGNET: An Open Source Tool for Analyzing and Visualizing M/EEG Connectome

The brain is a large-scale complex network often referred to as the “connectome”. Exploring the dynamic behavior of the connectome is a challenging issue as both excellent time and space resolution is required. In this context Magneto/Electroencephalography (M/EEG) are effective neuroimaging techniques allowing for analysis of the dynamics of functional brain networks at scalp level and/or at reconstructed sources. However, a tool that can cover all the processing steps of identifying brain networks from M/EEG data is still missing. In this paper, we report a novel software package, called EEGNET, running under MATLAB (Math works, inc), and allowing for analysis and visualization of functional brain networks from M/EEG recordings. EEGNET is developed to analyze networks either at the level of scalp electrodes or at the level of reconstructed cortical sources. It includes i) Basic steps in preprocessing M/EEG signals, ii) the solution of the inverse problem to localize / reconstruct the cortical sources, iii) the computation of functional connectivity among signals collected at surface electrodes or/and time courses of reconstructed sources and iv) the computation of the network measures based on graph theory analysis. EEGNET is the unique tool that combines the M/EEG functional connectivity analysis and the computation of network measures derived from the graph theory. The first version of EEGNET is easy to use, flexible and user friendly. EEGNET is an open source tool and can be freely downloaded from this webpage: https://sites.google.com/site/eegnetworks/.

Systems Neuroengineering: Understanding and Interacting with the Brain

In this paper, we review the current state-of-the-art techniques used for understanding the inner workings of the brain at a systems level. The neural activity that governs our everyday lives involves an intricate coordination of many processes that can be attributed to a variety of brain regions. On the surface, many of these functions can appear to be controlled by specific anatomical structures; however, in reality, numerous dynamic networks within the brain contribute to its function through an interconnected web of neuronal and synaptic pathways. The brain, in its healthy or pathological state, can therefore be best understood by taking a systems-level approach. While numerous neuroengineering technologies exist, we focus here on three major thrusts in the field of systems neuroengineering: neuroimaging, neural interfacing, and neuromodulation. Neuroimaging enables us to delineate the structural and functional organization of the brain, which is key in understanding how the neural system functions in both normal and disease states. Based on such knowledge, devices can be used either to communicate with the neural system, as in neural interface systems, or to modulate brain activity, as in neuromodulation systems. The consideration of these three fields is key to the development and application of neuro-devices. Feedback-based neuro-devices require the ability to sense neural activity (via a neuroimaging modality) through a neural interface (invasive or noninvasive) and ultimately to select a set of stimulation parameters in order to alter neural function via a neuromodulation modality. Systems neuroengineering refers to the use of engineering tools and technologies to image, decode, and modulate the brain in order to comprehend its functions and to repair its dysfunction. Interactions between these fields will help to shape the future of systems neuroengineering-to develop neurotechniques for enhancing the understanding of whole-brain function and dysfunction, and the management of neurological and mental disorders.

Brain Networks Responsible for Sense of Agency: An EEG Study

BACKGROUND

Self-agency (SA) is a person's feeling that his action was generated by himself. The neural substrates of SA have been investigated in many neuroimaging studies, but the functional connectivity of identified regions has rarely been investigated. The goal of this study is to investigate the neural network related to SA.

METHODS

SA of hand movements was modulated with virtual reality. We examined the cortical network relating to SA modulation with electroencephalography (EEG) power spectrum and phase coherence of alpha, beta, and gamma frequency bands in 16 right-handed, healthy volunteers.

RESULTS

In the alpha band, significant relative power changes and phase coherence of alpha band were associated with SA modulation. The relative power decrease over the central, bilateral parietal, and right temporal regions (C4, Pz, P3, P4, T6) became larger as participants more effectively controlled the virtual hand movements. The phase coherence of the alpha band within frontal areas (F7-FP2, F7-Fz) was directly related to changes in SA. The functional connectivity was lower as the participants felt that they could control their virtual hand. In the other frequency bands, significant phase coherences were observed in the frontal (or central) to parietal, temporal, and occipital regions during SA modulation (Fz-O1, F3-O1, Cz-O1, C3-T4L in beta band; FP1-T6, FP1-O2, F7-T4L, F8-Cz in gamma band).

CONCLUSIONS

Our study suggests that alpha band activity may be the main neural oscillation of SA, which suggests that the neural network within the anterior frontal area may be important in the generation of SA.

Estimating underlying neuronal activity from EEG using an iterative sparse technique

In this paper a novel technique for solving the bio-electromagnetic inverse problem is proposed. This method provides information about the location and extent of underlying neuronal activity. This is essential for the presurgical planning for partial epilepsy patients who are resistant to anti-epileptic drugs. The proposed algorithm takes advantage of the fact that neuronal activity transparent to EEG, arises from a spatially extended brain region. This spatial coherence is modeled within the framework of sparse signal processing techniques and makes better use of the limited number of EEG recordings. An iterative data-driven weighting is also introduced to better the extent estimation as well as eliminating the need to threshold estimated solutions.

Noninvasive Brain-Computer Interfaces Based on Sensorimotor Rhythms

Brain-computer interfaces (BCIs) have been explored in the field of neuroengineering to investigate how the brain can use these systems to control external devices. We review the principles and approaches we have taken to develop a sensorimotor rhythm EEG based brain-computer interface (BCI). The methods include developing BCI systems incorporating the control of physical devices to increase user engagement, improving BCI systems by inversely mapping scalp-recorded EEG signals to the cortical source domain, integrating BCI with noninvasive neuromodulation strategies to improve learning, and incorporating mind-body awareness training to enhance BCI learning and performance. The challenges and merits of these strategies are discussed, together with recent findings. Our work indicates that the sensorimotor-rhythm-based noninvasive BCI has the potential to provide communication and control capabilities as an alternative to physiological motor pathways.

STRAPS: A Fully Data-Driven Spatio-Temporally Regularized Algorithm for M/EEG Patch Source Imaging

For M/EEG-based distributed source imaging, it has been established that the L2-norm-based methods are effective in imaging spatially extended sources, whereas the L1-norm-based methods are more suited for estimating focal and sparse sources. However, when the spatial extents of the sources are unknown a priori, the rationale for using either type of methods is not adequately supported. Bayesian inference by exploiting the spatio-temporal information of the patch sources holds great promise as a tool for adaptive source imaging, but both computational and methodological limitations remain to be overcome. In this paper, based on state-space modeling of the M/EEG data, we propose a fully data-driven and scalable algorithm, termed STRAPS, for M/EEG patch source imaging on high-resolution cortices. Unlike the existing algorithms, the recursive penalized least squares (RPLS) procedure is employed to efficiently estimate the source activities as opposed to the computationally demanding Kalman filtering/smoothing. Furthermore, the coefficients of the multivariate autoregressive (MVAR) model characterizing the spatial-temporal dynamics of the patch sources are estimated in a principled manner via empirical Bayes. Extensive numerical experiments demonstrate STRAPS's excellent performance in the estimation of locations, spatial extents and amplitudes of the patch sources with varying spatial extents.

Network community structure detection for directional neural networks inferred from multichannel multisubject EEG data

In many neuroscience applications, one is interested in identifying the functional brain modules from multichannel, multiple subject neuroimaging data. However, most of the existing network community structure detection algorithms are limited to single undirected networks and cannot reveal the common community structure for a collection of directed networks. In this paper, we propose a community detection algorithm for weighted asymmetric (directed) networks representing the effective connectivity in the brain. Moreover, the issue of finding a common community structure across subjects is addressed by maximizing the total modularity of the group. Finally, the proposed community detection algorithm is applied to multichannel multisubject electroencephalogram data.

A nonlinear generalization of spectral Granger causality

Spectral measures of linear Granger causality have been widely applied to study the causal connectivity between time series data in neuroscience, biology, and economics. Traditional Granger causality measures are based on linear autoregressive with exogenous (ARX) inputs models of time series data, which cannot truly reveal nonlinear effects in the data especially in the frequency domain. In this study, it is shown that the classical Geweke's spectral causality measure can be explicitly linked with the output spectra of corresponding restricted and unrestricted time-domain models. The latter representation is then generalized to nonlinear bivariate signals and for the first time nonlinear causality analysis in the frequency domain. This is achieved by using the nonlinear ARX (NARX) modeling of signals, and decomposition of the recently defined output frequency response function which is related to the NARX model.

Effect of EEG electrode number on epileptic source localization in pediatric patients

OBJECTIVE

To investigate the relationship between EEG source localization and the number of scalp EEG recording channels.

METHODS

128 EEG channel recordings of 5 pediatric patients with medically intractable partial epilepsy were used to perform source localization of interictal spikes. The results were compared with surgical resection and intracranial recordings. Various electrode configurations were tested and a series of computer simulations based on a realistic head boundary element model were also performed in order to further validate the clinical findings.

RESULTS

The improvement seen in source localization substantially decreases as the number of electrodes increases. This finding was evaluated using the surgical resection, intracranial recordings and computer simulation. It was also shown in the simulation that increasing the electrode numbers could remedy the localization error of deep sources. A plateauing effect was seen in deep and superficial sources with further increasing the electrode number.

CONCLUSION

The source localization is improved when electrode numbers increase, but the absolute improvement in accuracy decreases with increasing electrode number.

SIGNIFICANCE

Increasing the electrode number helps decrease localization error and thus can more ably assist the physician to better plan for surgical procedures.

Influence of the head model on EEG and MEG source connectivity analyses

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently. In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problems. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only the skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC). Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC.

Asymmetric engagement of amygdala and its gamma connectivity in early emotional face processing

The amygdala has been regarded as a key substrate for emotion processing. However, the engagement of the left and right amygdala during the early perceptual processing of different emotional faces remains unclear. We investigated the temporal profiles of oscillatory gamma activity in the amygdala and effective connectivity of the amygdala with the thalamus and cortical areas during implicit emotion-perceptual tasks using event-related magnetoencephalography (MEG). We found that within 100 ms after stimulus onset the right amygdala habituated to emotional faces rapidly (with duration around 20–30 ms), whereas activity in the left amygdala (with duration around 50–60 ms) sustained longer than that in the right. Our data suggest that the right amygdala could be linked to autonomic arousal generated by facial emotions and the left amygdala might be involved in decoding or evaluating expressive faces in the early perceptual emotion processing. The results of effective connectivity provide evidence that only negative emotional processing engages both cortical and subcortical pathways connected to the right amygdala, representing its evolutional significance (survival). These findings demonstrate the asymmetric engagement of bilateral amygdala in emotional face processing as well as the capability of MEG for assessing thalamo-cortico-limbic circuitry.

Brain-computer interfaces using sensorimotor rhythms: current state and future perspectives

Many studies over the past two decades have shown that people can use brain signals to convey their intent to a computer using brain-computer interfaces (BCIs). BCI systems extract specific features of brain activity and translate them into control signals that drive an output. Recently, a category of BCIs that are built on the rhythmic activity recorded over the sensorimotor cortex, i.e., the sensorimotor rhythm (SMR), has attracted considerable attention among the BCIs that use noninvasive neural recordings, e.g., electroencephalography (EEG), and have demonstrated the capability of multidimensional prosthesis control. This paper reviews the current state and future perspectives of SMR-based BCI and its clinical applications, in particular focusing on the EEG SMR. The characteristic features of SMR from the human brain are described and their underlying neural sources are discussed. The functional components of SMR-based BCI, together with its current clinical applications, are reviewed. Finally, limitations of SMR-BCIs and future outlooks are also discussed.

Characterizing nonlinear relationships in functional imaging data using eigenspace maximal information canonical correlation analysis (emiCCA)

Many important problems in the analysis of neuroimages can be formulated as discovering the relationship between two sets of variables, a task for which linear techniques such as canonical correlation analysis (CCA) have been commonly used. However, to further explore potential nonlinear processes that might co-exist with linear ones in brain function, a more flexible method is required. Here, we propose a new unsupervised and data-driven method, termed the eigenspace maximal information canonical correlation analysis (emiCCA), which is capable of automatically capturing the linear and/or nonlinear relationships between various data sets. A simulation confirmed the superior performance of emiCCA in comparison with linear CCA and kernel CCA (a nonlinear version of CCA). An emiCCA framework for functional magnetic resonance imaging (fMRI) data processing was designed and applied to data from a real motor execution fMRI experiment. This analysis uncovered one linear (in primary motor cortex) and a few nonlinear networks (e.g., in the supplementary motor area, bilateral insula, and cerebellum). This suggests that these various task-related brain areas are part of networks that also contribute to the execution of movements of the hand. These results suggest that emiCCA is a promising technique for exploring various data.

Effective connectivity maps in the swine somatosensory cortex estimated from electrocorticography and validated with intracortical local field potential measurements

Macroscopic techniques are increasingly being used to estimate functional connectivity in the brain, which provides valuable information about brain networks. In any such endeavors it is important to understand capabilities and limitations of each technique through direct validation, which is often lacking. This study evaluated a multiple dipole source analysis technique based on electrocorticography (ECOG) data in estimating effective connectivity maps and validated the technique with intracortical local field potential (LFP) recordings. The study was carried out in an animal model (swine) with a large brain to avoid complications caused by spreading of the volume current. The evaluation was carried out for the cortical projections from the trigeminal nerve and corticocortical connectivity from the first rostrum area (R1) in the primary somatosensory cortex. Stimulation of the snout and layer IV of the R1 did not activate all projection areas in each animal, although whenever an area was activated in a given animal, its location was consistent with the intracortical LFP. The two types of connectivity maps based on ECOG analysis were consistent with each other and also with those estimated from the intracortical LFP, although there were small discrepancies. The discrepancies in mean latency based on ECOG and LFP were all very small and nonsignificant: snout stimulation, -1.1-2.0 msec (contralateral hemisphere) and 3.9-8.5 msec (ipsilateral hemisphere); R1 stimulation, -1.4-2.2 msec for the ipsilateral and 0.6-1.4 msec for the contralateral hemisphere. Dipole source analysis based on ECOG appears to be quite useful for estimating effective connectivity maps in the brain.

Correspondence of electroencephalography and near-infrared spectroscopy sensitivities to the cerebral cortex using a high-density layout

This study investigates the correspondence of the cortical sensitivity of electroencephalography (EEG) and near-infrared spectroscopy (NIRS). EEG forward model sensitivity to the cerebral cortex was calculated for 329 EEG electrodes following the 10-5 EEG positioning system using a segmented structural magnetic resonance imaging scan of a human subject. NIRS forward model sensitivity was calculated for the same subject using 156 NIRS source-detector pairs selected from 32 source and 32 detector optodes positioned on the scalp using a subset of the 10-5 EEG positioning system. Sensitivity correlations between colocalized NIRS source-detector pair groups and EEG channels yielded R = 0.46 ± 0.08. Groups of NIRS source-detector pairs with maximum correlations to EEG electrode sensitivities are tabulated. The mean correlation between the point spread functions for EEG and NIRS regions of interest (ROI) was R = 0.43 ± 0.07. Spherical ROIs with radii of 26 mm yielded the maximum correlation between EEG and NIRS averaged across all cortical mesh nodes. These sensitivity correlations between EEG and NIRS should be taken into account when designing multimodal studies of neurovascular coupling and when using NIRS as a statistical prior for EEG source localization.

Quantification of the effects of volume conduction on the EEG/MEG connectivity estimates: an index of sensitivity to brain interactions

In the context of EEG/MEG, the term 'volume conduction (VC) effects' refers to the recording of an instantaneous linear mixture of multiple brain source activities by each EEG/MEG channel. VC effects may lead to the detection of spurious functional/effective couplings among EEG/MEG channels that are not caused by brain interactions. It is of importance to determine which detected couplings are indicators of brain interactions and which originate from the VC artefacts. In this paper, a quantitative framework is proposed to explore the origin of detected channel couplings by using two types of surrogate datasets. Also, a sensitivity index (called SI) is proposed to compare the power of different connectivity measures to discriminate between the brain interactions and the instantaneous linear mixing effects. We use seven different functional connectivity estimators to evaluate our method on simulation models and resting state EEG data. The error rate of the proposed framework for simulation data by using each of the connectivity estimators is less than 5.2%. Also, SI ranks these connectivity estimators according to their sensitivity to brain interactions in the presence of VC artefacts. As expected, the connectivity measures which are theoretically robust to VC artefacts yield high SI in simulation models and EEG data. In addition, for EEG data in the alpha frequency band the reproducible functional couplings which are indicators of brain interactions are in the back-front directions. This is consistent with the previous studies in this field.

Surface-based mixed effects multilevel analysis of grouped human electrocorticography

Electrocorticography (ECoG) in humans yields data with unmatched spatio-temporal resolution that provides novel insights into cognitive operations. However, the broader application of ECoG has been confounded by difficulties in accurately depicting individual data and performing statistically valid population-level analyses. To overcome these limitations, we developed methods for accurately registering ECoG data to individual cortical topology. We integrated this technique with surface-based co-registration and a mixed-effects multilevel analysis (MEMA) to control for variable cortical surface anatomy and sparse coverage across patients, as well as intra- and inter-subject variability. We applied this surface-based MEMA (SB-MEMA) technique to a face-recognition task dataset (n=22). Compared against existing techniques, SB-MEMA yielded results much more consistent with individual data and with meta-analyses of face-specific activation studies. We anticipate that SB-MEMA will greatly expand the role of ECoG in studies of human cognition, and will enable the generation of population-level brain activity maps and accurate multimodal comparisons.

Functional brain connectivity from EEG in epilepsy: seizure prediction and epileptogenic focus localization

Today, neuroimaging techniques are frequently used to investigate the integration of functionally specialized brain regions in a network. Functional connectivity, which quantifies the statistical dependencies among the dynamics of simultaneously recorded signals, allows to infer the dynamical interactions of segregated brain regions. In this review we discuss how the functional connectivity patterns obtained from intracranial and scalp electroencephalographic (EEG) recordings reveal information about the dynamics of the epileptic brain and can be used to predict upcoming seizures and to localize the seizure onset zone. The added value of extracting information that is not visibly identifiable in the EEG data using functional connectivity analysis is stressed. Despite the fact that many studies have showed promising results, we must conclude that functional connectivity analysis has not made its way into clinical practice yet.

Assessing dynamic spectral causality by lagged adaptive directed transfer function and instantaneous effect factor

It is of significance to assess the dynamic spectral causality among physiological signals. Several practical estimators adapted from spectral Granger causality have been exploited to track dynamic causality based on the framework of time-varying multivariate autoregressive (tvMVAR) models. The nonzero covariance of the model's residuals has been used to describe the instantaneous effect phenomenon in some causality estimators. However, for the situations with Gaussian residuals in some autoregressive models, it is challenging to distinguish the directed instantaneous causality if the sufficient prior information about the "causal ordering" is missing. Here, we propose a new algorithm to assess the time-varying causal ordering of tvMVAR model under the assumption that the signals follow the same acyclic causal ordering for all time lags and to estimate the instantaneous effect factor (IEF) value in order to track the dynamic directed instantaneous connectivity. The time-lagged adaptive directed transfer function (ADTF) is also estimated to assess the lagged causality after removing the instantaneous effect. In this study, we first investigated the performance of the causal-ordering estimation algorithm and the accuracy of IEF value. Then, we presented the results of IEF and time-lagged ADTF method by comparing with the conventional ADTF method through simulations of various propagation models. Statistical analysis results suggest that the new algorithm could accurately estimate the causal ordering and give a good estimation of the IEF values in the Gaussian residual conditions. Meanwhile, the time-lagged ADTF approach is also more accurate in estimating the time-lagged dynamic interactions in a complex nervous system after extracting the instantaneous effect. In addition to the simulation studies, we applied the proposed method to estimate the dynamic spectral causality on real visual evoked potential (VEP) data in a human subject. Its usefulness in time-variant spectral causality assessment was demonstrated through the mutual causality investigation of brain activity during the VEP experiments.

A novel dynamic update framework for epileptic seizure prediction

Epileptic seizure prediction is a difficult problem in clinical applications, and it has the potential to significantly improve the patients' daily lives whose seizures cannot be controlled by either drugs or surgery. However, most current studies of epileptic seizure prediction focus on high sensitivity and low false-positive rate only and lack the flexibility for a variety of epileptic seizures and patients' physical conditions. Therefore, a novel dynamic update framework for epileptic seizure prediction is proposed in this paper. In this framework, two basic sample pools are constructed and updated dynamically. Furthermore, the prediction model can be updated to be the most appropriate one for the prediction of seizures' arrival. Mahalanobis distance is introduced in this part to solve the problem of side information, measuring the distance between two data sets. In addition, a multichannel feature extraction method based on Hilbert-Huang transform and extreme learning machine is utilized to extract the features of a patient's preseizure state against the normal state. At last, a dynamic update epileptic seizure prediction system is built up. Simulations on Freiburg database show that the proposed system has a better performance than the one without update. The research of this paper is significantly helpful for clinical applications, especially for the exploitation of online portable devices.

Noninvasive imaging of the high frequency brain activity in focal epilepsy patients

High-frequency (HF) activity represents a potential biomarker of the epileptogenic zone in epilepsy patients, the removal of which is considered to be crucial for seizure-free surgical outcome. We proposed a high frequency source imaging (HFSI) approach to noninvasively image the brain sources of the scalp-recorded HF EEG activity. Both computer simulation and clinical patient data analysis were performed to investigate the feasibility of using the HFSI approach to image the sources of HF activity from noninvasive scalp EEG recordings. The HF activity was identified from high-density scalp recordings after high-pass filtering the EEG data and the EEG segments with HF activity were concatenated together to form repetitive HF activity. Independent component analysis was utilized to extract the components corresponding to the HF activity. Noninvasive EEG source imaging using realistic geometric boundary element head modeling was then applied to image the sources of the pathological HF brain activity. Five medically intractable focal epilepsy patients were studied and the estimated sources were found to be concordant with the surgical resection or intracranial recordings of the patients. The present study demonstrates, for the first time, that source imaging from the scalp HF activity could help to localize the seizure onset zone and provide a novel noninvasive way of studying the epileptic brain in humans. This study also indicates the potential application of studying HF activity in the presurgical planning of medically intractable epilepsy patients.