Source localization of EEG/MEG data by correlating columns of ICA and lead field matrices

Involvement of White Matter Language Tracts in Glioma: Clinical Implications, Operative Management, and Functional Recovery After Injury

To achieve optimal survival and quality of life outcomes in patients with glioma, the extent of tumor resection must be maximized without causing injury to eloquent structures. Preservation of language function is of particular importance to patients and requires careful mapping to reveal the locations of cortical language hubs and their structural and functional connections. Within this language network, accurate mapping of eloquent white matter tracts is critical, given the high risk of permanent neurological impairment if they are injured during surgery. In this review, we start by describing the clinical implications of gliomas involving white matter language tracts. Next, we highlight the advantages and limitations of methods commonly used to identify these tracts during surgery including structural imaging techniques, functional imaging, non-invasive stimulation, and finally, awake craniotomy. We provide a rationale for combining these complementary techniques as part of a multimodal mapping paradigm to optimize postoperative language outcomes. Next, we review local and long-range adaptations that take place as the language network undergoes remodeling after tumor growth and surgical resection. We discuss the probable cellular mechanisms underlying this plasticity with emphasis on the white matter, which until recently was thought to have a limited role in adults. Finally, we provide an overview of emerging developments in targeting the glioma-neuronal network interface to achieve better disease control and promote recovery after injury.

Embedding decomposition for artifacts removal in EEG signals

Electroencephalogram (EEG) recordings are often contaminated with artifacts. Various methods have been developed to eliminate or weaken the influence of artifacts. However, most of them rely on prior experience for analysis. Here, we propose an deep learning framework to separate neural signal and artifacts in the embedding space and reconstruct the denoised signal, which is called DeepSeparator. DeepSeparator employs an encoder to extract and amplify the features in the raw EEG, a module called decomposer to extract the trend, detect and suppress artifact and a decoder to reconstruct the denoised signal. Besides, DeepSeparator can extract the artifact, which largely increases the model interpretability. The proposed method is tested with a semi-synthetic EEG dataset and a real task-related EEG dataset, suggesting that DeepSeparator outperforms the conventional models in both EOG and EMG artifact removal. DeepSeparator can be extended to multi-channel EEG and data with any arbitrary length. It may motivate future developments and application of deep learning-based EEG denoising. The code for DeepSeparator is available at https://github.com/ncclabsustech/DeepSeparator.

A tutorial on generalized eigendecomposition for denoising, contrast enhancement, and dimension reduction in multichannel electrophysiology

The goal of this paper is to present a theoretical and practical introduction to generalized eigendecomposition (GED), which is a robust and flexible framework used for dimension reduction and source separation in multichannel signal processing. In cognitive electrophysiology, GED is used to create spatial filters that maximize a researcher-specified contrast. For example, one may wish to exploit an assumption that different sources have different frequency content, or that sources vary in magnitude across experimental conditions. GED is fast and easy to compute, performs well in simulated and real data, and is easily adaptable to a variety of specific research goals. This paper introduces GED in a way that ties together myriad individual publications and applications of GED in electrophysiology, and provides sample MATLAB and Python code that can be tested and adapted. Practical considerations and issues that often arise in applications are discussed.

Nonlinear Modeling of Cortical Responses to Mechanical Wrist Perturbations Using the NARMAX Method

OBJECTIVE

Nonlinear modeling of cortical responses (EEG) to wrist perturbations allows for the quantification of cortical sensorimotor function in healthy and neurologically impaired individuals. A common model structure reflecting key characteristics shared across healthy individuals may provide a reference for future clinical studies investigating abnormal cortical responses associated with sensorimotor impairments. Thus, the goal of our study is to identify this common model structure and therefore to build a nonlinear dynamic model of cortical responses, using nonlinear autoregressive-moving-average model with exogenous inputs (NARMAX).

METHODS

EEG was recorded from ten participants when receiving continuous wrist perturbations. A common model structure detection method was developed for identifying a common NARMAX model structure across all participants, with individualized parameter values. The results were compared to conventional subject-specific models.

RESULTS

The proposed method achieved 93.91% variance accounted for (VAF) when implementing a one-step-ahead prediction and around 50% VAF for a k-step ahead prediction (k = 3), without a substantial drop of VAF as compare to subject-specific models. The estimated common structure suggests that the measured cortical response is a mixed outcome of the nonlinear transformation of external inputs and local neuronal interactions or inherent neuronal dynamics at the cortex.

CONCLUSION

The proposed method well determined the common characteristics across subjects in the cortical responses to wrist perturbations.

SIGNIFICANCE

It provides new insights into the human sensorimotor nervous system in response to somatosensory inputs and paves the way for future translational studies on assessments of sensorimotor impairments using our modeling approach.

Rhythmic entrainment source separation: Optimizing analyses of neural responses to rhythmic sensory stimulation

Steady-state evoked potentials (SSEPs) are rhythmic brain responses to rhythmic sensory stimulation, and are often used to study perceptual and attentional processes. We present a data analysis method for maximizing the signal-to-noise ratio of the narrow-band steady-state response in the frequency and time-frequency domains. The method, termed rhythmic entrainment source separation (RESS), is based on denoising source separation approaches that take advantage of the simultaneous but differential projection of neural activity to multiple electrodes or sensors. Our approach is a combination and extension of existing multivariate source separation methods. We demonstrate that RESS performs well on both simulated and empirical data, and outperforms conventional SSEP analysis methods based on selecting electrodes with the strongest SSEP response, as well as several other linear spatial filters. We also discuss the potential confound of overfitting, whereby the filter captures noise in absence of a signal. Matlab scripts are available to replicate and extend our simulations and methods. We conclude with some practical advice for optimizing SSEP data analyses and interpreting the results.

Blind Source Separation of Event-Related EEG/MEG

OBJECTIVE

Blind source separation (BSS) can be used to decompose complex electroencephalography (EEG) or magnetoencephalography data into simpler components based on statistical assumptions without using a physical model. Applications include brain-computer interfaces, artifact removal, and identifying parallel neural processes. We wish to address the issue of applying BSS to event-related responses, which is challenging because of nonstationary data.

METHODS

We introduce a new BSS approach called momentary-uncorrelated component analysis (MUCA), which is tailored for event-related multitrial data. The method is based on approximate joint diagonalization of multiple covariance matrices estimated from the data at separate latencies. We further show how to extend the methodology for autocovariance matrices and how to apply BSS methods suitable for piecewise stationary data to event-related responses. We compared several BSS approaches by using simulated EEG as well as measured somatosensory and transcranial magnetic stimulation (TMS) evoked EEG.

RESULTS

Among the compared methods, MUCA was the most tolerant one to noise, TMS artifacts, and other challenges in the data. With measured somatosensory data, over half of the estimated components were found to be similar by MUCA and independent component analysis. MUCA was also stable when tested with several input datasets.

CONCLUSION

MUCA is based on simple assumptions, and the results suggest that MUCA is robust with nonideal data.

SIGNIFICANCE

Event-related responses and BSS are valuable and popular tools in neuroscience. Correctly designed BSS is an efficient way of identifying artifactual and neural processes from nonstationary event-related data.

Ocular and cardiac artifact rejection for real-time analysis in MEG

BACKGROUND

Recently, magnetoencephalography (MEG) based real-time brain computing interfaces (BCI) have been developed to enable novel and promising methods for neuroscience research. It is well known that artifact rejection prior to source localization largely enhances the localization accuracy. However, many BCI approaches neglect real-time artifact removal due to its time consuming process.

NEW METHOD

The method (referred to as ocular and cardiac artifact rejection for real-time analysis, OCARTA) is based on constrained independent component analysis (cICA), where a priori information of the underlying source signals is used to optimize and accelerate signal decomposition. Thereby, prior information is incorporated by using the subject's individual cardiac and ocular activity. The algorithm automatically uses different separation strategies depending on the underlying source activity.

RESULTS

OCARTA was tested and applied to data from three different but most commonly used MEG systems (4D-Neuroimaging, VSM MedTech Inc. and Elekta Neuromag). Ocular and cardiac artifacts were effectively reduced within one iteration at a time delay of 1ms performed on a standard PC (Intel Core i5-2410M).

COMPARISON WITH EXISTING METHODS

The artifact rejection results achieved with OCARTA are in line with the results reported for offline ICA-based artifact rejection methods.

CONCLUSION

Due to the fast and subject-specific signal decomposition the new approach introduced here is capable of real-time ocular and cardiac artifact rejection.

Independent component analysis of instantaneous power-based fMRI

In functional magnetic resonance imaging (fMRI) studies using spatial independent component analysis (sICA) method, a model of "latent variables" is often employed, which is based on the assumption that fMRI data are linear mixtures of statistically independent signals. However, actual fMRI signals are nonlinear and do not automatically meet with the requirement of sICA. To provide a better solution to this problem, we proposed a novel approach termed instantaneous power based fMRI (ip-fMRI) for regularization of fMRI data. Given that the instantaneous power of fMRI signals is a scalar value, it should be a linear mixture that naturally satisfies the "latent variables" model. Based on our simulated data, the curves of accuracy and resulting receiver-operating characteristic curves indicate that the proposed approach is superior to the traditional fMRI in terms of accuracy and specificity by using sICA. Experimental results from human subjects have shown that spatial components of a hand movement task-induced activation reveal a brain network more specific to motor function by ip-fMRI than that by the traditional fMRI. We conclude that ICA decomposition of ip-fMRI may be used to localize energy signal changes in the brain and may have a potential to be applied to detection of brain activity.

Uncovering neural independent components from highly artifactual TMS-evoked EEG data

Transcranial magnetic stimulation (TMS) combined with electroencephalography (EEG) is a powerful tool for studying cortical excitability and connectivity. To enhance the EEG interpretation, independent component analysis (ICA) has been used to separate the data into independent components (ICs). However, TMS can evoke large artifacts in EEG, which may greatly distort the ICA separation. The removal of such artifactual EEG from the data is a difficult task. In this paper we study how badly the large artifacts distort the ICA separation, and whether the distortions could be avoided without removing the artifacts. We first show that, in the ICA separation, the time courses of the ICs are not affected by the large artifacts, but their topographies could be greatly distorted. Next, we show how this distortion can be circumvented. We introduce a novel technique of suppression, by which the EEG data are modified so that the ICA separation of the suppressed data becomes reliable. The suppression, instead of removing the artifactual EEG, rescales all the data to about the same magnitude as the neural EEG. For the suppressed data, ICA returns the original time courses, but instead of the original topographies, it returns modified ones, which can be used, e.g., for the source localization. We present three suppression methods based on principal component analysis, wavelet analysis, and whitening of the data matrix, respectively. We test the methods with numerical simulations. The results show that the suppression improves the source localization.

Blind Image Separation Method Based on Artificial Bee Colony Algorithm

A novel blind image separation algorithm based on artificial bee colony algorithm was proposed. Negative entropy was chosen as the objective function for separation and the separation matrix was transformed using Givens rotation transform method and the number of variable need to be found reduced. Artificial bee colony algorithm was used for solving the objective function and the images could be separated successfully. Simulation results verify the validity of the method.