Choosing Strategies to Deal with Artifactual EEG Data in Children with Cognitive Impairment.

Rett syndrome is a disease that involves acute cognitive impairment and, consequently, a complex and varied symptomatology. This study evaluates the EEG signals of twenty-nine patients and classify them according to the level of movement artifact. The main goal is to achieve an artifact rejection strategy that performs well in all signals, regardless of the artifact level. Two different methods have been studied: one based on the data distribution and the other based on the energy function, with entropy as its main component. The method based on the data distribution shows poor performance with signals containing high amplitude outliers. On the contrary, the method based on the energy function is more robust to outliers. As it does not depend on the data distribution, it is not affected by artifactual events. A double rejection strategy has been chosen, first on a motion signal (accelerometer or EEG low-pass filtered between 1 and 10 Hz) and then on the EEG signal. The results showed a higher performance when working combining both artifact rejection methods. The energy-based method, to isolate motion artifacts, and the data-distribution-based method, to eliminate the remaining lower amplitude artifacts were used. In conclusion, a new method that proves to be robust for all types of signals is designed.

Acute Sleep Deprivation Induces a Local Brain Transfer Information Increase in the Frontal Cortex in a Widespread Decrease Context.

Sleep deprivation (SD) has adverse effects on mental and physical health, affecting the cognitive abilities and emotional states. Specifically, cognitive functions and alertness are known to decrease after SD. The aim of this work was to identify the directional information transfer after SD on scalp EEG signals using transfer entropy (TE). Using a robust methodology based on EEG recordings of 18 volunteers deprived from sleep for 36 h, TE and spectral analysis were performed to characterize EEG data acquired every 2 h. Correlation between connectivity measures and subjective somnolence was assessed. In general, TE showed medium- and long-range significant decreases originated at the occipital areas and directed towards different regions, which could be interpreted as the transfer of predictive information from parieto-occipital activity to the rest of the head. Simultaneously, short-range increases were obtained for the frontal areas, following a consistent and robust time course with significant maps after 20 h of sleep deprivation. Changes during sleep deprivation in brain network were measured effectively by TE, which showed increased local connectivity and diminished global integration. TE is an objective measure that could be used as a potential measure of sleep pressure and somnolence with the additional property of directed relationships.

A Dynamic Fitting Strategy for Physiological Models: A Case Study of a Cardiorespiratory Model for the Simulation of Incremental Aerobic Exercise.

Using mathematical models of physiological systems in medicine has allowed for the development of diagnostic, treatment, and medical educational tools. However, their complexity restricts, in most cases, their application for predictive, preventive, and personalized purposes. Although there are strategies that reduce the complexity of applying models based on fitting techniques, most of them are focused on a single instant of time, neglecting the effect of the system's temporal evolution. The objective of this research was to introduce a dynamic fitting strategy for physiological models with an extensive array of parameters and a constrained amount of experimental data. The proposed strategy focused on obtaining better predictions based on the temporal trends in the system's parameters and being capable of predicting future states. The study utilized a cardiorespiratory model as a case study. Experimental data from a longitudinal study of healthy adult subjects undergoing aerobic exercise were used for fitting and validation. The model predictions obtained in a steady state using the proposed strategy and the traditional single-fit approach were compared. The most successful outcomes were primarily linked to the proposed strategy, exhibiting better overall results regarding accuracy and behavior than the traditional population fitting approach at a single instant in time. The results evidenced the usefulness of the dynamic fitting strategy, highlighting its use for predictive, preventive, and personalized applications.

A Novel Strategy to Fit and Validate Physiological Models: A Case Study of a Cardiorespiratory Model for Simulation of Incremental Aerobic Exercise.

Applying complex mathematical models of physiological systems is challenging due to the large number of parameters. Identifying these parameters through experimentation is difficult, and although procedures for fitting and validating models are reported, no integrated strategy exists. Additionally, the complexity of optimization is generally neglected when the number of experimental observations is restricted, obtaining multiple solutions or results without physiological justification. This work proposes a fitting and validation strategy for physiological models with many parameters under various populations, stimuli, and experimental conditions. A cardiorespiratory system model is used as a case study, and the strategy, model, computational implementation, and data analysis are described. Using optimized parameter values, model simulations are compared to those obtained using nominal values, with experimental data as a reference. Overall, a reduction in prediction error is achieved compared to that reported for model building. Furthermore, the behavior and accuracy of all the predictions in the steady state were improved. The results validate the fitted model and provide evidence of the proposed strategy's usefulness.

An Accurate and Real-Time Method for Resolving Superimposed Action Potentials in MultiUnit Recordings.

OBJECTIVE: Spike sorting of muscular and neural recordings requires separating action potentials that overlap in time (superimposed action potentials (APs)). We propose a new algorithm for resolving superimposed action potentials, and we test it on intramuscular EMG (iEMG) and intracortical recordings. METHODS: Discrete-time shifts of the involved APs are first selected based on a heuristic extension of the peel-off algorithm. Then, the time shifts that provide the minimal residual Euclidean norm are identified (Discrete Brute force Correlation (DBC)). The optimal continuous-time shifts are then estimated (High-Resolution BC (HRBC)). In Fusion HRBC (FHRBC), two other cost functions are used. A parallel implementation of the DBC and HRBC algorithms was developed. The performance of the algorithms was assessed on 11,000 simulated iEMG and 14,000 neural recording superpositions, including two to eight APs, and eight experimental iEMG signals containing four to eleven active motor units. The performance of the proposed algorithms was compared with that of the Branch-and-Bound (BB) algorithm using the Rank-Product (RP) method in terms of accuracy and efficiency. RESULTS: The average accuracy of the DBC, HRBC and FHRBC methods on the entire simulated datasets was 92.16±17.70, 93.65±16.89, and 94.90±15.15 (%). The DBC algorithm outperformed the other algorithms based on the RP method. The average accuracy and running time of the DBC algorithm on 10.5 ms superimposed spikes of the experimental signals were 92.1±21.7 (%) and 2.3±15.3 (ms). CONCLUSION AND SIGNIFICANCE: The proposed algorithm is promising for real-time neural decoding, a central problem in neural and muscular decoding and interfacing.

Electroencephalographic assessment in patients with Rett syndrome during cognitive stimulation by means of eye tracking technology and alternative and augmentative communication systems.

Rett syndrome (RTT) is considered a rare disease despite being the leading genetic disorder to cause severe intellectual disability in women. There is no cure for RTT, so the treatment is symptomatic and supporting, requiring a multidisciplinary approach. Occupational therapy can help girls and their families to improve communication, being one of the main concerns when verbal language and intentional hand movement are impaired or lost. This paper presents a pilot study of cognitive training through the combined use of eye-tracking technology (ETT) and augmentative and alternative communication (AAC) using the Peabody Picture Vocabulary Test (PPVT-IV). The objective was to evaluate brain activation by means of electroencephalography (EEG) during the stimulation of non-verbal communication. EEG data were recorded during an eyes-open resting state (EO-RS) period and during cognitive stimulation via AAC activity. To assess their effect, both signals were compared at the spectral level, focusing on frequency, brain symmetry and connectivity. During the task, a redistribution of power towards fast frequency bands was observed, as well as an improvement in the brain symmetry index (BSI) and functional synchronicity through increased coherence. Therefore, the results of the spectral analysis showed a possible deviation from the pathological pattern, manifesting a positive effect in the use of non-verbal cognitive stimulation activities. In conclusion, it was observed that it is possible to establish a cognitive training system that produces brain activation and favors communication and learning despite intentional language loss.Clinical Relevance- This manifests a method of cognitive training that would induce brain activation in RTT patients with absence of intentional communication. The evaluation system through spectral analysis could complement the standardized protocols to asses communication that are based on verbal and motor production.

Spinal Cord Injury Patients Exhibit Changes in Motor-Related Activity and Topographic Distribution.

Spinal Cord Injury (SCI) is a common disease that usually limits the patient's independence by affecting their motor function. SCI patients usually present neuroplasticity, which allows brain signals transmission through spread pathways. Some innovative rehabilitation therapies, such as functional electrical stimulation (FES) or Brain-computer interfaces (BCIs) jointly with motor neuroprostheses, provide hope for functional restoration. BCIs require the analysis of event-related EEG potentials (ERPs). Movement-related cortical potentials (MRCPs) and event-related desynchroni-zation and synchronization (ERD/ERS) are the most commonly studied ERPs during motor activity. ERPs of healthy subjects may vary from SCI patients. Thus, this study aimed to compare ERPs between healthy subjects and SCI patients during upper-limb movements (forearm supination and pronation, and hand open). Differences between controls and SCI patients were shown in terms of ERPs' amplitude as well as in topographic maps. Changes in amplitude were more substantial in ERD potentials than in MRCPs, while topographic maps showed better localization of all features in healthy patients. The level of SCI injury determines the patients' mobility. A comparison between complete, partial and no motor function subjects showed lower values of feature's amplitudes in the latter group.Clinical Relevance- This demonstrates the existence of significant statistical differences between healthy and SCI subjects, and might be helpful when performing SCI rehabilitation techniques such as designing BCI and neuroprostheses, or analyzing and understanding the brain plasticity process.

High-density surface EMG maps from upper-arm and forearm muscles.

BACKGROUND: sEMG signal has been widely used in different applications in kinesiology and rehabilitation as well as in the control of human-machine interfaces. In general, the signals are recorded with bipolar electrodes located in different muscles. However, such configuration may disregard some aspects of the spatial distribution of the potentials like location of innervation zones and the manifestation of inhomogineties in the control of the muscular fibers. On the other hand, the spatial distribution of motor unit action potentials has recently been assessed with activation maps obtained from High Density EMG signals (HD-EMG), these lasts recorded with arrays of closely spaced electrodes. The main objective of this work is to analyze patterns in the activation maps, associating them with four movement directions at the elbow joint and with different strengths of those tasks. Although the activation pattern can be assessed with bipolar electrodes, HD-EMG maps could enable the extraction of features that depend on the spatial distribution of the potentials and on the load-sharing between muscles, in order to have a better differentiation between tasks and effort levels. METHODS: An experimental protocol consisting of isometric contractions at three levels of effort during flexion, extension, supination and pronation at the elbow joint was designed and HD-EMG signals were recorded with 2D electrode arrays on different upper-limb muscles. Techniques for the identification and interpolation of artifacts are explained, as well as a method for the segmentation of the activation areas. In addition, variables related to the intensity and spatial distribution of the maps were obtained, as well as variables associated to signal power of traditional single bipolar recordings. Finally, statistical tests were applied in order to assess differences between information extracted from single bipolar signals or from HD-EMG maps and to analyze differences due to type of task and effort level. RESULTS: Significant differences were observed between EMG signal power obtained from single bipolar configuration and HD-EMG and better results regarding the identification of tasks and effort levels were obtained with the latter. Additionally, average maps for a population of 12 subjects were obtained and differences in the co-activation pattern of muscles were found not only from variables related to the intensity of the maps but also to their spatial distribution. CONCLUSIONS: Intensity and spatial distribution of HD-EMG maps could be useful in applications where the identification of movement intention and its strength is needed, for example in robotic-aided therapies or for devices like powered- prostheses or orthoses. Finally, additional data transformations or other features are necessary in order to improve the performance of tasks identification.

Influence of the number of trials on evoked motor cortical activity in EEG recordings.

Objective. Improvements in electroencephalography enable the study of the localization of active brain regions during motor tasks. Movement-related cortical potentials (MRCPs), and event-related desynchronization (ERD) and synchronization are the main motor-related cortical phenomena/neural correlates observed when a movement is elicited. When assessing neurological diseases, averaging techniques are commonly applied to characterize motor related processes better. In this case, a large number of trials is required to obtain a motor potential that is representative enough of the subject's condition. This study aimed to assess the effect of a limited number of trials on motor-related activity corresponding to different upper limb movements (elbow flexion/extension, pronation/supination and hand open/close).Approach. An open dataset consisting on 15 healthy subjects was used for the analysis. A Monte Carlo simulation approach was applied to analyse, in a robust way, different typical time- and frequency-domain features, topography, and low-resolution electromagnetic tomography.Main results. Grand average potentials, and topographic and tomographic maps showed few differences when using fewer trials, but shifts in the localization of motor-related activity were found for several individuals. MRCP and beta ERD features were more robust to a limited number of trials, yielding differences lower than 20% for cases with 50 trials or more. Strong correlations between features were obtained for subsets above 50 trials. However, the inter-subject variability increased as the number of trials decreased. The elbow flexion/extension movement showed a more robust performance for a limited number of trials, both in population and in individual-based analysis.Significance. Our findings suggested that 50 trials can be an appropriate number to obtain stable motor-related features in terms of differences in the averaged motor features, correlation, and changes in topography and tomography.

Improving the ripple classification in focal pediatric epilepsy: identifying pathological high-frequency oscillations by Gaussian mixture model clustering.

Objective. High-frequency oscillations (HFOs) have emerged as a promising clinical biomarker for presurgical evaluation in childhood epilepsy. HFOs are commonly classified in stereo-encephalography as ripples (80-200 Hz) and fast ripples (200-500 Hz). Ripples are less specific and not so directly associated with epileptogenic activity because of their physiological and pathological origin. The aim of this paper is to distinguish HFOs in the ripple band and to improve the evaluation of the epileptogenic zone (EZ).Approach. This study constitutes a novel modeling approach evaluated in ten patients from Sant Joan de Deu Pediatric Hospital (Barcelona, Spain), with clearly-defined seizure onset zones (SOZ) during presurgical evaluation. A subject-by-subject basis analysis is proposed: a probabilistic Gaussian mixture model (GMM) based on the combination of specific ripple features is applied for estimating physiological and pathological ripple subpopulations.Main Results. Clear pathological and physiological ripples are identified. Features differ considerably among patients showing within-subject variability, suggesting that individual models are more appropriate than a traditional whole-population approach. The difference in rates inside and outside the SOZ for pathological ripples is significantly higher than when considering all the ripples. These significant differences also appear in signal segments without epileptiform activity. Pathological ripple rates show a sharp decline from SOZ to non-SOZ contacts and a gradual decrease with distance.Significance. This novel individual GMM approach improves ripple classification and helps to refine the delineation of the EZ, as well as being appropriate to investigate the interaction of epileptogenic and propagation networks.

Drug effect on EEG connectivity assessed by linear and nonlinear couplings.

Quantitative analysis of human electroencephalogram (EEG) is a valuable method for evaluating psychopharmacological agents. Although the effects of different drug classes on EEG spectra are already known, interactions between brain locations remain unclear. In this work, cross mutual information function and appropriate surrogate data were applied to assess linear and nonlinear couplings between EEG signals. The main goal was to evaluate the pharmacological effects of alprazolam on brain connectivity during wakefulness in healthy volunteers using a cross-over, placebo-controlled design. Eighty-five pairs of EEG leads were selected for the analysis, and connectivity was evaluated inside anterior, central, and posterior zones of the scalp. Connectivity between these zones and interhemispheric connectivity were also measured. Results showed that alprazolam induced significant changes in EEG connectivity in terms of information transfer in comparison with placebo. Trends were opposite depending on the statistical characteristics: decreases in linear connectivity and increases in nonlinear couplings. These effects were generally spread over the entire scalp. Linear changes were negatively correlated, and nonlinear changes were positively correlated with drug plasma concentrations; the latter showed higher correlation coefficients. The use of both linear and nonlinear approaches revealed the importance of assessing changes in EEG connectivity as this can provide interesting information about psychopharmacological effects.

SGM: a novel time-frequency algorithm based on unsupervised learning improves high-frequency oscillation detection in epilepsy.

OBJECTIVE: We propose a novel automated method called the S-Transform Gaussian Mixture detection algorithm (SGM) to detect high-frequency oscillations (HFO) combining the strengths of different families of previously published detectors. APPROACH: This algorithm does not depend on parameter tuning on a subject (or database) basis, uses time-frequency characteristics, and relies on non-supervised classification to determine if the events standing out from the baseline activity are HFO or not. SGM consists of three steps: the first stage computes the signal baseline using the entropy of the autocorrelation; the second uses the S-Transform to obtain several time-frequency features (area, entropy, and time and frequency widths); and in the third stage Gaussian mixture models cluster time-frequency features to decide if events correspond to HFO-like activity. To validate the SGM algorithm we tested its performance in simulated and real environments. MAIN RESULTS: We assessed the algorithm on a publicly available simulated stereoelectroencephalographic (SEEG) database with varying signal-to-noise ratios (SNR), obtaining very good results for medium and high SNR signals. We further tested the SGM algorithm on real signals from patients with focal epilepsy, in which HFO detection was performed visually by experts, yielding a high agreement between experts and SGM. SIGNIFICANCE: The SGM algorithm displayed proper performance in simulated and real environments and therefore can be used for non-supervised detection of HFO. This non-supervised algorithm does not require previous labelling by experts or parameter adjustment depending on the subject or database considered. SGM is not a computationally intensive algorithm, making it suitable to detect and characterize HFO in long-term SEEG recordings.

Influence of ocular filtering in EEG data on the assessment of drug-induced effects on the brain.

Ocular artifacts in EEG signals affect the interpretation of clinical study results. The aim of this study was to assess the influence of automatic ocular filtering procedures in the conclusions drawn from a pharmaco-EEG trial. Regression analysis, gold standard, and blind source separation (BSS), Second Order Blind Identification algorithm, ocular filtering procedures were compared using time, frequency, topographic and tomographic brain mapping approaches and pharmacokinetic-pharmacodynamic (PK-PD) relationships. Data consisted of EEGs obtained from 20 volunteers who received single oral doses of haloperidol 3 mg, risperidone 1 mg, olanzapine 5 mg and placebo in a randomized cross-over double-blind design. Although the BSS-based technique preserved brain activity more than regression analysis in anterior leads, in general, topographic significance probability maps globally showed similar results with both methods for most spectral variables. However, different results were obtained when using whole multi-channel information for studying drug effects in the brain: (i) higher correlations between PK and PD time courses showing that BSS allowed estimation of spectral variables more accurately related to drug effects and (ii) larger and more symmetric drug related tomographic LORETA maps showing that BSS led to results that were more neurophysiopharmacologically sound. Definitely, the BSS-based procedure is an effective and efficient preprocessing method to remove ocular artifacts from EEG data. The selection of the ocular filtering procedure could determine different results whose impact depends on the evaluating tool applied to analyze the pharmaco-EEG data.

Alterations in EEG connectivity in healthy young adults provide an indicator of sleep depth.

Current sleep analyses have used electroencephalography (EEG) to establish sleep intensity through linear and nonlinear measures. Slow wave activity (SWA) and entropy are the most commonly used markers of sleep depth. The purpose of this study is to evaluate changes in brain EEG connectivity during sleep in healthy subjects and compare them with SWA and entropy. Four different connectivity metrics: coherence (MSC), synchronization likelihood (SL), cross mutual information function (CMIF), and phase locking value (PLV), were computed focusing on their correlation with sleep depth. These measures provide different information and perspectives about functional connectivity. All connectivity measures revealed to have functional changes between the different sleep stages. The averaged CMIF seemed to be a more robust connectivity metric to measure sleep depth (correlations of 0.78 and 0.84 with SWA and entropy, respectively), translating greater linear and nonlinear interdependences between brain regions especially during slow wave sleep. Potential changes of brain connectivity were also assessed throughout the night. Connectivity measures indicated a reduction of functional connectivity in N2 as sleep progresses. The validation of connectivity indexes is necessary because they can reveal the interaction between different brain regions in physiological and pathological conditions and help understand the different functions of deep sleep in humans.

A comparative study of automatic techniques for ocular artifact reduction in spontaneous EEG signals based on clinical target variables: a simulation case.

Eye movement artifacts represent a critical issue for quantitative electroencephalography (EEG) analysis and a number of mathematical approaches have been proposed to reduce their contribution in EEG recordings. The aim of this paper was to objectively and quantitatively evaluate the performance of ocular filtering methods with respect to spectral target variables widely used in clinical and functional EEG studies. In particular the following methods were applied: regression analysis and some blind source separation (BSS) techniques based on second-order statistics (PCA, AMUSE and SOBI) and on higher-order statistics (JADE, INFOMAX and FASTICA). Considering blind source decomposition methods, a completely automatic procedure of BSS based on logical rules related to spectral and topographical information was proposed in order to identify the components related to ocular interference. The automatic procedure was applied in different montages of simulated EEG and electrooculography (EOG) recordings: a full montage with 19 EEG and 2 EOG channels, a reduced one with only 6 EEG leads and a third one where EOG channels were not available. Time and frequency results in all of them indicated that AMUSE and SOBI algorithms preserved and recovered more brain activity than the other methods mainly at anterior regions. In the case of full montage: (i) errors were lower than 5% for all spectral variables at anterior sites; and (ii) the highest improvement in the signal-to-artifact (SAR) ratio was obtained up to 40dB at these anterior sites. Finally, we concluded that second-order BSS-based algorithms (AMUSE and SOBI) provided an effective technique for eye movement removal even when EOG recordings were not available or when data length was short.

An Improved Dynamic Model for the Respiratory Response to Exercise.

Respiratory system modeling has been extensively studied in steady-state conditions to simulate sleep disorders, to predict its behavior under ventilatory diseases or stimuli and to simulate its interaction with mechanical ventilation. Nevertheless, the studies focused on the instantaneous response are limited, which restricts its application in clinical practice. The aim of this study is double: firstly, to analyze both dynamic and static responses of two known respiratory models under exercise stimuli by using an incremental exercise stimulus sequence (to analyze the model responses when step inputs are applied) and experimental data (to assess prediction capability of each model). Secondly, to propose changes in the models' structures to improve their transient and stationary responses. The versatility of the resulting model vs. the other two is shown according to the ability to simulate ventilatory stimuli, like exercise, with a proper regulation of the arterial blood gases, suitable constant times and a better adjustment to experimental data. The proposed model adjusts the breathing pattern every respiratory cycle using an optimization criterion based on minimization of work of breathing through regulation of respiratory frequency.

Evaluation of respiratory muscles activity by means of cross mutual information function at different levels of ventilatory effort.

Analysis of respiratory muscles activity is an effective technique for the study of pulmonary diseases such as obstructive sleep apnea syndrome (OSAS). Respiratory diseases, especially those associated with changes in the mechanical properties of the respiratory apparatus, are often associated with disruptions of the normally highly coordinated contractions of respiratory muscles. Due to the complexity of the respiratory control, the assessment of OSAS related dysfunctions by linear methods are not sufficient. Therefore, the objective of this study was the detection of diagnostically relevant nonlinear complex respiratory mechanisms. Two aims of this work were: (1) to assess coordination of respiratory muscles contractions through evaluation of interactions between respiratory signals and myographic signals through nonlinear analysis by means of cross mutual information function (CMIF); (2) to differentiate between functioning of respiratory muscles in patients with OSAS and in normal subjects. Electromyographic (EMG) and mechanomyographic (MMG) signals were recorded from three respiratory muscles: genioglossus, sternomastoid and diaphragm. Inspiratory pressure and flow were also acquired. All signals were measured in eight patients with OSAS and eight healthy subjects during an increased respiratory effort while awake. Several variables were defined and calculated from CMIF in order to describe correlation between signals. The results indicate different nonlinear couplings of respiratory muscles in both populations. This effect is progressively more evident at higher levels of respiratory effort.

Automated detection of epileptic ripples in MEG using beamformer-based virtual sensors.

OBJECTIVE: In epilepsy, high-frequency oscillations (HFOs) are expressively linked to the seizure onset zone (SOZ). The detection of HFOs in the noninvasive signals from scalp electroencephalography (EEG) and magnetoencephalography (MEG) is still a challenging task. The aim of this study was to automate the detection of ripples in MEG signals by reducing the high-frequency noise using beamformer-based virtual sensors (VSs) and applying an automatic procedure for exploring the time-frequency content of the detected events. APPROACH: Two-hundred seconds of MEG signal and simultaneous iEEG were selected from nine patients with refractory epilepsy. A two-stage algorithm was implemented. Firstly, beamforming was applied to the whole head to delimitate the region of interest (ROI) within a coarse grid of MEG-VS. Secondly, a beamformer using a finer grid in the ROI was computed. The automatic detection of ripples was performed using the time-frequency response provided by the Stockwell transform. Performance was evaluated through comparisons with simultaneous iEEG signals. MAIN RESULTS: ROIs were located within the seizure-generating lobes in the nine subjects. Precision and sensitivity values were 79.18% and 68.88%, respectively, by considering iEEG-detected events as benchmarks. A higher number of ripples were detected inside the ROI compared to the same region in the contralateral lobe. SIGNIFICANCE: The evaluation of interictal ripples using non-invasive techniques can help in the delimitation of the epileptogenic zone and guide placement of intracranial electrodes. This is the first study that automatically detects ripples in MEG in the time domain located within the clinically expected epileptic area taking into account the time-frequency characteristics of the events through the whole signal spectrum. The algorithm was tested against intracranial recordings, the current gold standard. Further studies should explore this approach to enable the localization of noninvasively recorded HFOs to help during pre-surgical planning and to reduce the need for invasive diagnostics.

Influence of metallic artifact filtering on MEG signals for source localization during interictal epileptiform activity.

OBJECTIVE: Medical intractable epilepsy is a common condition that affects 40% of epileptic patients that generally have to undergo resective surgery. Magnetoencephalography (MEG) has been increasingly used to identify the epileptogenic foci through equivalent current dipole (ECD) modeling, one of the most accepted methods to obtain an accurate localization of interictal epileptiform discharges (IEDs). Modeling requires that MEG signals are adequately preprocessed to reduce interferences, a task that has been greatly improved by the use of blind source separation (BSS) methods. MEG recordings are highly sensitive to metallic interferences originated inside the head by implanted intracranial electrodes, dental prosthesis, etc and also coming from external sources such as pacemakers or vagal stimulators. To reduce these artifacts, a BSS-based fully automatic procedure was recently developed and validated, showing an effective reduction of metallic artifacts in simulated and real signals (Migliorelli et al 2015 J. Neural Eng. 12 046001). The main objective of this study was to evaluate its effects in the detection of IEDs and ECD modeling of patients with focal epilepsy and metallic interference. APPROACH: A comparison between the resulting positions of ECDs was performed: without removing metallic interference; rejecting only channels with large metallic artifacts; and after BSS-based reduction. Measures of dispersion and distance of ECDs were defined to analyze the results. MAIN RESULTS: The relationship between the artifact-to-signal ratio and ECD fitting showed that higher values of metallic interference produced highly scattered dipoles. Results revealed a significant reduction on dispersion using the BSS-based reduction procedure, yielding feasible locations of ECDs in contrast to the other two approaches. SIGNIFICANCE: The automatic BSS-based method can be applied to MEG datasets affected by metallic artifacts as a processing step to improve the localization of epileptic foci.

Automatic BSS-based filtering of metallic interference in MEG recordings: definition and validation using simulated signals.

OBJECTIVE: One of the principal drawbacks of magnetoencephalography (MEG) is its high sensitivity to metallic artifacts, which come from implanted intracranial electrodes and dental ferromagnetic prosthesis and produce a high distortion that masks cerebral activity. The aim of this study was to develop an automatic algorithm based on blind source separation (BSS) techniques to remove metallic artifacts from MEG signals. APPROACH: Three methods were evaluated: AMUSE, a second-order technique; and INFOMAX and FastICA, both based on high-order statistics. Simulated signals consisting of real artifact-free data mixed with real metallic artifacts were generated to objectively evaluate the effectiveness of BSS and the subsequent interference reduction. A completely automatic detection of metallic-related components was proposed, exploiting the known characteristics of the metallic interference: regularity and low frequency content. MAIN RESULTS: The automatic procedure was applied to the simulated datasets and the three methods exhibited different performances. Results indicated that AMUSE preserved and consequently recovered more brain activity than INFOMAX and FastICA. Normalized mean squared error for AMUSE decomposition remained below 2%, allowing an effective removal of artifactual components. SIGNIFICANCE: To date, the performance of automatic artifact reduction has not been evaluated in MEG recordings. The proposed methodology is based on an automatic algorithm that provides an effective interference removal. This approach can be applied to any MEG dataset affected by metallic artifacts as a processing step, allowing further analysis of unusable or poor quality data.