Modelling and detecting deep brain activity with MEG and EEG

We introduce an anatomical and electrophysiological model of deep brain structures dedicated to magnetoencephalography (MEG) and electroencephalography (EEG) source imaging. So far, most imaging inverse models considered that MEG/EEG surface signals were predominantly produced by cortical, hence superficial, neural currents. Here we question whether crucial deep brain structures such as the basal ganglia and the hippocampus may also contribute to distant, scalp MEG and EEG measurements. We first design a realistic anatomical and electrophysiological model of these structures and subsequently run Monte-Carlo experiments to evaluate the respective sensitivity of the MEG and EEG to signals from deeper origins. Results indicate that MEG/EEG may indeed localize these deeper generators, which is confirmed here from experimental MEG data reporting on the modulation of alpha (10–12 Hz) brain waves.     

Accounting for Linear Transformations of EEG and MEG Data in Source Analysis

Analyses of electro- and magnetoencephalography (EEG, MEG) data often involve a linear modification of signals at the sensor level. Examples include re-referencing of the EEG, computation of synthetic gradiometer in MEG, or the removal of artifactual independent components to clean EEG and MEG data. A question of practical relevance is, if such modifications must be accounted for by adapting the physical forward model (leadfield) before subsequent source analysis. Here, we show that two scenarios need to be differentiated. In the first scenario, which corresponds to re-referencing the EEG and synthetic gradiometer computation in MEG, the leadfield must be adapted before source analysis. In the second scenario, which corresponds to removing artifactual components to ‘clean’ the data, the leadfield must not be changed. We demonstrate and discuss the consequences of wrongly modifying the leadfield in the latter case for an example. Future EEG and MEG studies employing source analyses should carefully consider whether and, if so, how the leadfield must be modified as explicated here.     

A review of EEG and MEG for brainnetome research

The majority of brain activities are performed by functionally integrating separate regions of the brain. Therefore, the synchronous operation of the brain’s multiple regions or neuronal assemblies can be represented as a network with nodes that are interconnected by links. Because of the complexity of brain interactions and their varying effects at different levels of complexity, one of the corresponding authors of this paper recently proposed the brainnetome as a new –ome to explore and integrate the brain network at different scales. Because electroencephalography (EEG) and magnetoencephalography (MEG) are noninvasive and have outstanding temporal resolution and because they are the primary clinical techniques used to capture the dynamics of neuronal connections, they lend themselves to the analysis of the neural networks comprising the brainnetome. Because of EEG/MEG’s applicability to brainnetome analyses, the aim of this review is to identify the procedures that can be used to form a network using EEG/MEG data in sensor or source space and to promote EEG/MEG network analysis for either neuroscience or clinical applications. To accomplish this aim, we show the relationship of the brainnetome to brain networks at the macroscale and provide a systematic review of network construction using EEG and MEG. Some potential applications of the EEG/MEG brainnetome are to use newly developed methods to associate the properties of a brainnetome with indices of cognition or disease conditions. Associations based on EEG/MEG brainnetome analysis may improve the comprehension of the functioning of the brain in neuroscience research or the recognition of abnormal patterns in neurological disease.     

Combined mapping of human auditory EEG and MEG responses

Auditory electric and magnetic P50(m), N1(m) and MMN(m) responses to standard, deviant and novel sounds were studied by recording brain electrical activity with 25 EEG electrodes simultaneously with the corresponding magnetic signals measured with 122 MEG gradiometer coils. The sources of these responses were located on the basis of the MEG responses; all were found to be in the supratemporal plane. The goal of the present paper was to investigate to what degree the source locations and orientations determined from the magnetic data account for the measured EEG signals. It was found that the electric P50, N1 and MMN responses can to a considerable degree be explained by the sources of the corresponding magnetic responses. In addition, source-current components not detectable by MEG were shown to contribute to the measured EEG signals.     

Functional Mapping with Simultaneous MEG and EEG

We use magnetoencephalography (MEG) and electroencephalography (EEG) to locate and determine the temporal evolution in brain areas involved in the processing of simple sensory stimuli. We will use somatosensory stimuli to locate the hand somatosensory areas, auditory stimuli to locate the auditory cortices, visual stimuli in four quadrants of the visual field to locate the early visual areas. These type of experiments are used for functional mapping in epileptic and brain tumor patients to locate eloquent cortices. In basic neuroscience similar experimental protocols are used to study the orchestration of cortical activity. The acquisition protocol includes quality assurance procedures, subject preparation for the combined MEG/EEG study, and acquisition of evoked-response data with somatosensory, auditory, and visual stimuli. We also demonstrate analysis of the data using the equivalent current dipole model and cortically-constrained minimum-norm estimates. Anatomical MRI data are employed in the analysis for visualization and for deriving boundaries of tissue boundaries for forward modeling and cortical location and orientation constraints for the minimum-norm estimates.Download video file.^{(69M, mp4)}     

Replicability of MEG and EEG measures of the auditory N1/N1m-response

We investigated the replicability of the source location, amplitude and latency measures of the auditory evoked N1 (EEG) and N1m (MEG) responses. Each of the 5 subjects was measured 6 times in two recording sessions. Responses to monaural stimuli were recorded from 122 MEG and 64 EEG channels simultaneously. The EEG data were modeled with a symmetrically-located dipole pair. For the MEG data, one dipole in each hemisphere was located independently using a subset of channels. Standard deviation (SD) was used as a measure for replicability. The average SD of the x, y and z coordinates of the contralateral N1m dipole was about 2 mm, whereas the corresponding figures for the ipsilateral N1m and the contra- and ipsilateral N1 were about twice as large. The SDs of the dipole amplitudes and latencies were almost equal with MEG and EEG. The amplitude and latency measures of the MEG field gradient waveforms were almost as replicable as those of the dipole models. The results suggest that both MEG and EEG can be used for investigating the simultaneous activity of the left and right auditory cortices independently, MEG being superior in certain experimental setups.     

Effects of Spatial Pattern Scale of Brain Activity on the Sensitivity of DOT,fMRI, EEG and MEG

The objective of this work is to quantify how patterns of cortical activity at different spatial scales are measured by noninvasive functional neuroimaging sensors. We simulated cortical activation patterns at nine different spatial scales in a realistic head model and propagated this activity to magnetoencephalography (MEG), electroencephalography (EEG), diffuse optical tomography (DOT), and functional magnetic resonance imaging (fMRI) sensors in arrangements that are typically used in functional neuroimaging studies. We estimated contrast transfer functions (CTF), correlation distances in sensor space, and the minimum resolvable spatial scale of cortical activity for each modality. We found that CTF decreases as the spatial extent of cortical activity decreases, and that correlations between nearby sensors depend on the spatial extent of cortical activity. For cortical activity on the intermediate spatial scale of 6.7 cm², the correlation distances (r>0.5) were 1.0 cm for fMRI, 2.0 cm for DOT, 12.8 for EEG, 9.5 cm for MEG magnetometers and 9.7 cm for MEG gradiometers. The resolvable spatial pattern scale was found to be 1.43 cm² for MEG magnetometers, 0.88 cm² for MEG gradiometers, 376 cm² for EEG, 0.75 cm² for DOT, and 0.072 cm² for fMRI. These findings show that sensitivity to cortical activity varies substantially as a function of spatial scale within and between the different imaging modalities. This information should be taken into account when interpreting neuroimaging data and when choosing the number of nodes for network analyses in sensor space.     

Classification methods for ongoing EEG and MEG signals

Classification algorithms help predict the qualitative properties of a subject's mental state by extracting useful information from the highly multivariate non-invasive recordings of his brain activity. In particular, applying them to Magneto-encephalography (MEG) and electro-encephalography (EEG) is a challenging and promising task with prominent practical applications to e.g. Brain Computer Interface (BCI). In this paper, we first review the principles of the major classification techniques and discuss their application to MEG and EEG data classification. Next, we investigate the behavior of classification methods using real data recorded during a MEG visuomotor experiment. In particular, we study the influence of the classification algorithm, of the quantitative functional variables used in this classifier, and of the validation method. In addition, our findings suggest that by investigating the distribution of classifier coefficients, it is possible to infer knowledge and construct functional interpretations of the underlying neural mechanisms of the performed tasks. Finally, the promising results reported here (up to 97% classification accuracy on 1-second time windows) reflect the considerable potential of MEG for the continuous classification of mental states.     

Divergent Cortical Generators of MEG and EEG during Human Sleep Spindles Suggested by Distributed Source Modeling

Background

Sleep spindles are ∼1-second bursts of 10–15 Hz activity, occurring during normal stage 2 sleep. In animals, sleep spindles can be synchronous across multiple cortical and thalamic locations, suggesting a distributed stable phase-locked generating system. The high synchrony of spindles across scalp EEG sites suggests that this may also be true in humans. However, prior MEG studies suggest multiple and varying generators.

Methodology/Principal Findings

We recorded 306 channels of MEG simultaneously with 60 channels of EEG during naturally occurring spindles of stage 2 sleep in 7 healthy subjects. High-resolution structural MRI was obtained in each subject, to define the shells for a boundary element forward solution and to reconstruct the cortex providing the solution space for a noise-normalized minimum norm source estimation procedure. Integrated across the entire duration of all spindles, sources estimated from EEG and MEG are similar, diffuse and widespread, including all lobes from both hemispheres. However, the locations, phase and amplitude of sources simultaneously estimated from MEG versus EEG are highly distinct during the same spindles. Specifically, the sources estimated from EEG are highly synchronous across the cortex, whereas those from MEG rapidly shift in phase, hemisphere, and the location within the hemisphere.

Conclusions/Significance

The heterogeneity of MEG sources implies that multiple generators are active during human sleep spindles. If the source modeling is correct, then EEG spindles are generated by a different, diffusely synchronous system. Animal studies have identified two thalamo-cortical systems, core and matrix, that produce focal or diffuse activation and thus could underlie MEG and EEG spindles, respectively. Alternatively, EEG spindles could reflect overlap at the sensors of the same sources as are seen from the MEG. Although our results generally match human intracranial recordings, additional improvements are possible and simultaneous intra- and extra-cranial measures are needed to test their accuracy.     

Coupled oscillators for modeling and analysis of EEG/MEG oscillations.

This study presents three EEG/MEG applications in which the modeling of oscillatory signal components offers complementary analysis and an improved explanation of the underlying generator structures. Coupled oscillator networks were used for modeling. Parameters of the corresponding ordinary coupled differential equation (ODE) system are identified using EEG/MEG data and the resulting solution yields the modeled signals. This model-related analysis strategy provides information about the coupling quantity and quality between signal components (example 1, neonatal EEG during quiet sleep), allows identification of the possible contribution of hidden generator structures (example 2, 600-Hz MEG oscillations in somatosensory evoked magnetic fields), and can explain complex signal characteristics such as amplitude-frequency coupling and frequency entrainment (example 3, EEG burst patterns in sedated patients).     

Application of modern tests for stationarity to single-trial MEG data

Stationarity is a crucial yet rarely questioned assumption in the analysis of time series of magneto- (MEG) or electroencephalography (EEG). One key drawback of the commonly used tests for stationarity of encephalographic time series is the fact that conclusions on stationarity are only indirectly inferred either from the Gaussianity (e.g. the Shapiro-Wilk test or Kolmogorov-Smirnov test) or the randomness of the time series and the absence of trend using very simple time-series models (e.g. the sign and trend tests by Bendat and Piersol). We present a novel approach to the analysis of the stationarity of MEG and EEG time series by applying modern statistical methods which were specifically developed in econometrics to verify the hypothesis that a time series is stationary. We report our findings of the application of three different tests of stationarity--the Kwiatkowski-Phillips-Schmidt-Schin (KPSS) test for trend or mean stationarity, the Phillips-Perron (PP) test for the presence of a unit root and the White test for homoscedasticity--on an illustrative set of MEG data. For five stimulation sessions, we found already for short epochs of duration of 250 and 500 ms that, although the majority of the studied epochs of single MEG trials were usually mean-stationary (KPSS test and PP test), they were classified as nonstationary due to their heteroscedasticity (White test). We also observed that the presence of external auditory stimulation did not significantly affect the findings regarding the stationarity of the data. We conclude that the combination of these tests allows a refined analysis of the stationarity of MEG and EEG time series.     

Harmony: EEG/MEG Linear Inverse Source Reconstruction in the Anatomical Basis of Spherical Harmonics

EEG/MEG source localization based on a “distributed solution” is severely underdetermined, because the number of sources is much larger than the number of measurements. In particular, this makes the solution strongly affected by sensor noise. A new way to constrain the problem is presented. By using the anatomical basis of spherical harmonics (or spherical splines) instead of single dipoles the dimensionality of the inverse solution is greatly reduced without sacrificing the quality of the data fit. The smoothness of the resulting solution reduces the surface bias and scatter of the sources (incoherency) compared to the popular minimum-norm algorithms where single-dipole basis is used (MNE, depth-weighted MNE, dSPM, sLORETA, LORETA, IBF) and allows to efficiently reduce the effect of sensor noise. This approach, termed Harmony, performed well when applied to experimental data (two exemplars of early evoked potentials) and showed better localization precision and solution coherence than the other tested algorithms when applied to realistically simulated data.     

A Pharmacokinetics-Neural Mass Model (PK-NMM) for the Simulation of EEG Activity during Propofol Anesthesia

Modeling the effects of anesthetic drugs on brain activity is very helpful in understanding anesthesia mechanisms. The aim of this study was to set up a combined model to relate actual drug levels to EEG dynamics and behavioral states during propofol-induced anesthesia. We proposed a new combined theoretical model based on a pharmacokinetics (PK) model and a neural mass model (NMM), which we termed PK-NMM—with the aim of simulating electroencephalogram (EEG) activity during propofol-induced general anesthesia. The PK model was used to derive propofol effect-site drug concentrations (C _eff) based on the actual drug infusion regimen. The NMM model took C _eff as the control parameter to produce simulated EEG-like (sEEG) data. For comparison, we used real prefrontal EEG (rEEG) data of nine volunteers undergoing propofol anesthesia from a previous experiment. To see how well the sEEG could describe the dynamic changes of neural activity during anesthesia, the rEEG data and the sEEG data were compared with respect to: power-frequency plots; nonlinear exponent (permutation entropy (PE)); and bispectral SynchFastSlow (SFS) parameters. We found that the PK-NMM model was able to reproduce anesthesia EEG-like signals based on the estimated drug concentration and patients’ condition. The frequency spectrum indicated that the frequency power peak of the sEEG moved towards the low frequency band as anesthesia deepened. Different anesthetic states could be differentiated by the PE index. The correlation coefficient of PE was 0.80±0.13 (mean±standard deviation) between rEEG and sEEG for all subjects. Additionally, SFS could track the depth of anesthesia and the SFS of rEEG and sEEG were highly correlated with a correlation coefficient of 0.77±0.13. The PK-NMM model could simulate EEG activity and might be a useful tool for understanding the action of propofol on brain activity.     

A fast and reliable method for simultaneous waveform, amplitude and latency estimation of single-trial EEG/MEG data

The amplitude and latency of single-trial EEG/MEG signals may provide valuable information concerning human brain functioning. In this article we propose a new method to reliably estimate single-trial amplitude and latency of EEG/MEG signals. The advantages of the method are fourfold. First, no a-priori specified template function is required. Second, the method allows for multiple signals that may vary independently in amplitude and/or latency. Third, the method is less sensitive to noise as it models data with a parsimonious set of basis functions. Finally, the method is very fast since it is based on an iterative linear least squares algorithm. A simulation study shows that the method yields reliable estimates under different levels of latency variation and signal-to-noise ratio?s. Furthermore, it shows that the existence of multiple signals can be correctly determined. An application to empirical data from a choice reaction time study indicates that the method describes these data accurately.     

Application of a dissimilarity index of EEG and its sub-bands on prediction of induced epileptic seizures from rat's EEG signals

ObjectiveEpileptic seizures are defined as manifest of excessive and hyper-synchronous activity of neurons in the cerebral cortex that cause frequent malfunction of the human central nervous system. Therefore, finding precursors and predictors of epileptic seizure is of utmost clinical relevance to reduce the epileptic seizure induced nervous system malfunction consequences. Researchers for this purpose may even guide us to a deep understanding of the seizure generating mechanisms. The goal of this paper is to predict epileptic seizures in epileptic rats.MethodsSeizures were induced in rats using pentylenetetrazole (PTZ) model. EEG signals in interictal, preictal, ictal and postictal periods were then recorded and analyzed to predict epileptic seizures. Epileptic seizures were predicted by calculating an index in consecutive windows of EEG signal and comparing the index with a threshold. In this work, a newly proposed dissimilarity index called Bhattacharyya Based Dissimilarity Index (BBDI), dynamical similarity index and fuzzy similarity index were investigated.ResultsBBDI, dynamical similarity index and fuzzy similarity index were examined on case and control groups and compared to each other. The results show that BBDI outperforms dynamical and fuzzy similarity indices. In order to improve the results, EEG sub-bands were also analyzed. The best result achieved when the proposed dissimilarity index was applied on Delta sub-band that predicts epileptic seizures in all rats with a mean of 299.5 s.ConclusionThe dissimilarity of neural network activity between reference window and present window of EEG signal has a significant increase prior to an epileptic seizure and the proposed dissimilarity index (BBDI) can reveal this variation to predict epileptic seizures. In addition, analyzing EEG sub-bands results in more accurate information about constituent neuronal activities underlying the EEG since certain changes in EEG signal may be amplified when each sub-band is analyzed separately.SignificanceThis paper presents application of a dissimilarity index (BBDI) on EEG signals and its sub-bands to predict PTZ-induced epileptic seizures in rats. Based on the results of this work, BBDI will predict epileptic seizures more accurately and more reliably compared to current indices that increases epileptic patient comfort and improves patient outcomes.     

Mapping Cortical Dynamics Using Simultaneous MEG/EEG and Anatomically-constrained Minimum-norm Estimates: an Auditory Attention Example

Magneto- and electroencephalography (MEG/EEG) are neuroimaging techniques that provide a high temporal resolution particularly suitable to investigate the cortical networks involved in dynamical perceptual and cognitive tasks, such as attending to different sounds in a cocktail party. Many past studies have employed data recorded at the sensor level only, i.e., the magnetic fields or the electric potentials recorded outside and on the scalp, and have usually focused on activity that is time-locked to the stimulus presentation. This type of event-related field / potential analysis is particularly useful when there are only a small number of distinct dipolar patterns that can be isolated and identified in space and time. Alternatively, by utilizing anatomical information, these distinct field patterns can be localized as current sources on the cortex. However, for a more sustained response that may not be time-locked to a specific stimulus (e.g., in preparation for listening to one of the two simultaneously presented spoken digits based on the cued auditory feature) or may be distributed across multiple spatial locations unknown a priori, the recruitment of a distributed cortical network may not be adequately captured by using a limited number of focal sources.Here, we describe a procedure that employs individual anatomical MRI data to establish a relationship between the sensor information and the dipole activation on the cortex through the use of minimum-norm estimates (MNE). This inverse imaging approach provides us a tool for distributed source analysis. For illustrative purposes, we will describe all procedures using FreeSurfer and MNE software, both freely available. We will summarize the MRI sequences and analysis steps required to produce a forward model that enables us to relate the expected field pattern caused by the dipoles distributed on the cortex onto the M/EEG sensors. Next, we will step through the necessary processes that facilitate us in denoising the sensor data from environmental and physiological contaminants. We will then outline the procedure for combining and mapping MEG/EEG sensor data onto the cortical space, thereby producing a family of time-series of cortical dipole activation on the brain surface (or "brain movies") related to each experimental condition. Finally, we will highlight a few statistical techniques that enable us to make scientific inference across a subject population (i.e., perform group-level analysis) based on a common cortical coordinate space.     

A Physiology-Based Seizure Detection System for Multichannel EEG

Background

Epilepsy is a common chronic neurological disorder characterized by recurrent unprovoked seizures. Electroencephalogram (EEG) signals play a critical role in the diagnosis of epilepsy. Multichannel EEGs contain more information than do single-channel EEGs. Automatic detection algorithms for spikes or seizures have traditionally been implemented on single-channel EEG, and algorithms for multichannel EEG are unavailable.

Methodology

This study proposes a physiology-based detection system for epileptic seizures that uses multichannel EEG signals. The proposed technique was tested on two EEG data sets acquired from 18 patients. Both unipolar and bipolar EEG signals were analyzed. We employed sample entropy (SampEn), statistical values, and concepts used in clinical neurophysiology (e.g., phase reversals and potential fields of a bipolar EEG) to extract the features. We further tested the performance of a genetic algorithm cascaded with a support vector machine and post-classification spike matching.

Principal Findings

We obtained 86.69% spike detection and 99.77% seizure detection for Data Set I. The detection system was further validated using the model trained by Data Set I on Data Set II. The system again showed high performance, with 91.18% detection of spikes and 99.22% seizure detection.

Conclusion

We report a de novo EEG classification system for seizure and spike detection on multichannel EEG that includes physiology-based knowledge to enhance the performance of this type of system.     

Large-Scale Functional Networks Identified from Resting-State EEG Using Spatial ICA

Several methods have been applied to EEG or MEG signals to detect functional networks. In recent works using MEG/EEG and fMRI data, temporal ICA analysis has been used to extract spatial maps of resting-state networks with or without an atlas-based parcellation of the cortex. Since the links between the fMRI signal and the electromagnetic signals are not fully established, and to avoid any bias, we examined whether EEG alone was able to derive the spatial distribution and temporal characteristics of functional networks. To do so, we propose a two-step original method: 1) An individual multi-frequency data analysis including EEG-based source localisation and spatial independent component analysis, which allowed us to characterize the resting-state networks. 2) A group-level analysis involving a hierarchical clustering procedure to identify reproducible large-scale networks across the population. Compared with large-scale resting-state networks obtained with fMRI, the proposed EEG-based analysis revealed smaller independent networks thanks to the high temporal resolution of EEG, hence hierarchical organization of networks. The comparison showed a substantial overlap between EEG and fMRI networks in motor, premotor, sensory, frontal, and parietal areas. However, there were mismatches between EEG-based and fMRI-based networks in temporal areas, presumably resulting from a poor sensitivity of fMRI in these regions or artefacts in the EEG signals. The proposed method opens the way for studying the high temporal dynamics of networks at the source level thanks to the high temporal resolution of EEG. It would then become possible to study detailed measures of the dynamics of connectivity.     

Song Perception by Professional Singers and Actors: An MEG Study

The cortical correlates of speech and music perception are essentially overlapping, and the specific effects of different types of training on these networks remain unknown. We compared two groups of vocally trained professionals for music and speech, singers and actors, using recited and sung rhyme sequences from German art songs with semantic and/ or prosodic/melodic violations (i.e. violations of pitch) of the last word, in order to measure the evoked activation in a magnetoencephalographic (MEG) experiment. MEG data confirmed the existence of intertwined networks for the sung and spoken modality in an early time window after word violation. In essence for this early response, higher activity was measured after melodic/prosodic than semantic violations in predominantly right temporal areas. For singers as well as for actors, modality-specific effects were evident in predominantly left-temporal lateralized activity after semantic expectancy violations in the spoken modality, and right-dominant temporal activity in response to melodic violations in the sung modality. As an indication of a special group-dependent audiation process, higher neuronal activity for singers appeared in a late time window in right temporal and left parietal areas, both after the recited and the sung sequences.