Monte Carlo simulation studies of EEG and MEG localization accuracy

Source localization comparison and combination of OPM-MEG and fMRI to detect sensorimotor cortex responses

BACKGROUND AND OBJECTIVES

The exploration of various neuroimaging techniques have become focal points within the field of neuroscience research. Magnetoencephalography based on optically pumped magnetometers (OPM-MEG) has shown significant potential to be the next generation of functional neuroimaging with the advantages of high signal intensity and flexible sensor arrangement. In this study, we constructed a 31-channel OPM-MEG system and performed a preliminary comparison of the temporal and spatial relationship between magnetic responses measured by OPM-MEG and blood-oxygen-level-dependent signals detected by functional magnetic resonance imaging (fMRI) during a grasping task.

METHODS

For OPM-MEG, the β-band (15-30 Hz) oscillatory activities can be reliably detected across multiple subjects and multiple session runs. To effectively localize the inhibitory oscillatory activities, a source power-spectrum ratio-based imaging method was proposed. This approach was compared with conventional source imaging methods, such as minimum norm-type and beamformer methods, and was applied in OPM-MEG source analysis. Subsequently, the spatial and temporal responses at the source-level between OPM-MEG and fMRI were analyzed.

RESULTS

The effectiveness of the proposed method was confirmed through simulations compared to benchmark methods. Our demonstration revealed an average spatial separation of 10.57 ± 4.41 mm between the localization results of OPM-MEG and fMRI across four subjects. Furthermore, the fMRI-constrained OPM-MEG localization results indicated a more focused imaging extent.

CONCLUSIONS

Taken together, the performance exhibited by OPM-MEG positions it as a potential instrument for functional surgery assessment.

Resting state electroencephalographic alpha rhythms are sensitive to Alzheimer's disease mild cognitive impairment progression at a 6-month follow-up

Are posterior resting-state electroencephalographic (rsEEG) alpha rhythms sensitive to the Alzheimer's disease mild cognitive impairment (ADMCI) progression at a 6-month follow-up? Clinical, cerebrospinal, neuroimaging, and rsEEG datasets in 52 ADMCI and 60 Healthy old seniors (equivalent groups for demographic features) were available from an international archive (www.pdwaves.eu). The ADMCI patients were arbitrarily divided into two groups: REACTIVE and UNREACTIVE, based on the reduction (reactivity) in the posterior rsEEG alpha eLORETA source activities from the eyes-closed to eyes-open condition at ≥ -10% and -10%, respectively. 75% of the ADMCI patients were REACTIVE. Compared to the UNREACTIVE group, the REACTIVE group showed (1) less abnormal posterior rsEEG source activity during the eyes-closed condition and (2) a decrease in that activity at the 6-month follow-up. These effects could not be explained by neuroimaging and neuropsychological biomarkers of AD. Such a biomarker might reflect abnormalities in cortical arousal in quiet wakefulness to be used for clinical studies in ADMCI patients using 6-month follow-ups.

Deep EEG source localization via EMD-based fMRI high spatial frequency

Brain imaging with a high-spatiotemporal resolution is crucial for accurate brain-function mapping. Electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are two popular neuroimaging modalities with complementary features that record brain function with high temporal and spatial resolution, respectively. One popular non-invasive way to obtain data with both high spatial and temporal resolutions is to combine the fMRI activation map and EEG data to improve the spatial resolution of the EEG source localization. However, using the whole fMRI map may cause spurious results for the EEG source localization, especially for deep brain regions. Considering the head's conductivity, deep regions' sources with low activity are unlikely to be detected by the EEG electrodes at the scalp. In this study, we use fMRI's high spatial-frequency component to identify the local high-intensity activations that are most likely to be captured by the EEG. The 3D Empirical Mode Decomposition (3D-EMD), a data-driven method, is used to decompose the fMRI map into its spatial-frequency components. Different validation measurements for EEG source localization show improved performance for the EEG inverse-modeling informed by the fMRI's high-frequency spatial component compared to the fMRI-informed EEG source-localization methods. The level of improvement varies depending on the voxels' intensity and their distribution. Our experimental results also support this conclusion.

Oscillatory characteristics of resting-state magnetoencephalography reflect pathological and symptomatic conditions of cognitive impairment

Background

Dementia and mild cognitive impairment are characterised by symptoms of cognitive decline, which are typically assessed using neuropsychological assessments (NPAs), such as the Mini-Mental State Examination (MMSE) and Frontal Assessment Battery (FAB). Magnetoencephalography (MEG) is a novel clinical assessment technique that measures brain activities (summarised as oscillatory parameters), which are associated with symptoms of cognitive impairment. However, the relevance of MEG and regional cerebral blood flow (rCBF) data obtained using single-photon emission computed tomography (SPECT) has not been examined using clinical datasets. Therefore, this study aimed to investigate the relationships among MEG oscillatory parameters, clinically validated biomarkers computed from rCBF, and NPAs using outpatient data retrieved from hospital records.

Methods

Clinical data from 64 individuals with mixed pathological backgrounds were retrieved and analysed. MEG oscillatory parameters, including relative power (RP) from delta to high gamma bands, mean frequency, individual alpha frequency, and Shannon's spectral entropy, were computed for each cortical region. For SPECT data, three pathological parameters-'severity', 'extent', and 'ratio'-were computed using an easy z-score imaging system (eZIS). As for NPAs, the MMSE and FAB scores were retrieved.

Results

MEG oscillatory parameters were correlated with eZIS parameters. The eZIS parameters associated with Alzheimer's disease pathology were reflected in theta power augmentation and slower shift of the alpha peak. Moreover, MEG oscillatory parameters were found to reflect NPAs. Global slowing and loss of diversity in neural oscillatory components correlated with MMSE and FAB scores, whereas the associations between eZIS parameters and NPAs were sparse.

Conclusion

MEG oscillatory parameters correlated with both SPECT (i.e. eZIS) parameters and NPAs, supporting the clinical validity of MEG oscillatory parameters as pathological and symptomatic indicators. The findings indicate that various components of MEG oscillatory characteristics can provide valuable pathological and symptomatic information, making MEG data a rich resource for clinical examinations of patients with cognitive impairments. SPECT (i.e. eZIS) parameters showed no correlations with NPAs. The results contributed to a better understanding of the characteristics of electrophysiological and pathological examinations for patients with cognitive impairments, which will help to facilitate their co-use in clinical application, thereby improving patient care.

Simulation study on high spatio-temporal resolution acousto-electrophysiological neuroimaging

Objective.Acousto-electrophysiological neuroimaging (AENI) is a technique hypothesized to record electrophysiological activity of the brain with millimeter spatial and sub-millisecond temporal resolution. This improvement is obtained by tagging areas with focused ultrasound (fUS). Due to mechanical vibration with respect to the measuring electrodes, the electrical activity of the marked region will be modulated onto the ultrasonic frequency. The region's electrical activity can subsequently be retrieved via demodulation of the measured signal. In this study, the feasibility of this hypothesized technique is tested.Approach.This is done by calculating the forward electroencephalography response under quasi-static assumptions. The head is simplified as a set of concentric spheres. Two sizes are evaluated representing human and mouse brains. Moreover, feasibility is assessed for wet and dry transcranial, and for cortically placed electrodes. The activity sources are modeled by dipoles, with their current intensity profile drawn from a power-law power spectral density.Results.It is shown that mechanical vibration modulates the endogenous activity onto the ultrasonic frequency. The signal strength depends non-linearly on the alignment between dipole orientation, vibration direction and recording point. The strongest signal is measured when these three dependencies are perfectly aligned. The signal strengths are in the pV-range for a dipole moment of 5 nAm and ultrasonic pressures within Food and Drug Administration (FDA)-limits. The endogenous activity can then be accurately reconstructed via demodulation. Two interference types are investigated: vibrational and static. Depending on the vibrational interference, it is shown that millimeter resolution signal detection is possible also for deep brain regions. Subsequently, successful demodulation depends on the static interference, that at MHz-range has to be sub-picovolt.Significance.Our results show that mechanical vibration is a possible underlying mechanism of acousto-electrophyisological neuroimaging. This paper is a first step towards improved understanding of the conditions under which AENI is feasible.

Parietal resting-state EEG alpha source connectivity is associated with subcortical white matter lesions in HIV-positive people

OBJECTIVE

Parietal resting-state electroencephalographic (rsEEG) alpha (8-10 Hz) source connectivity is abnormal in HIV-positive persons. Here we tested whether this abnormality may be associated with subcortical white matter vascular lesions in the cerebral hemispheres.

METHODS

Clinical, rsEEG, and magnetic resonance imaging (MRI) datasets in 38 HIV-positive persons and clinical and rsEEG datasets in 13 healthy controls were analyzed. Radiologists visually evaluated the subcortical white matter hyperintensities from T2-weighted FLAIR MRIs (i.e., Fazekas scale). In parallel, neurophysiologists estimated the eLORETA rsEEG source lagged linear connectivity from parietal cortical regions of interest.

RESULTS

Compared to the HIV participants with no/negligible subcortical white matter hyperintensities, the HIV participants with mild/moderate subcortical white matter hyperintensities showed lower parietal interhemispheric rsEEG alpha lagged linear connectivity. This effect was also observed in HIV-positive persons with unimpaired cognition. This rsEEG marker allowed good discrimination (area under the receiver operating characteristic curve > 0.80) between the HIV-positive individuals with different amounts of subcortical white matter hyperintensities.

CONCLUSIONS

The parietal rsEEG alpha source connectivity is associated with subcortical white matter vascular lesions in HIV-positive persons, even without neurocognitive disorders.

SIGNIFICANCE

Those MRI-rsEEG markers may be used to screen HIV-positive persons at risk of neurocognitive disorders.

Electromagnetic source imaging predicts surgical outcome in children with focal cortical dysplasia

OBJECTIVE

To evaluate the diagnostic accuracy of electromagnetic source imaging (EMSI) in localizing spikes and predict surgical outcome in children with drug resistant epilepsy (DRE) due to focal cortical dysplasia (FCD).

METHODS

We retrospectively analyzed magnetoencephalography (MEG) and high-density (HD-EEG) data from 23 children with FCD-associated DRE who underwent intracranial EEG and surgery. We localized spikes using equivalent current dipole (ECD) fitting, dipole clustering, and dynamical statistical parametric mapping (dSPM) on EMSI, electric source imaging (ESI), and magnetic source imaging (MSI). We calculated the distance from the seizure onset zone (DSOZ) and resection (DRES). We estimated receiver operating characteristic (ROC) curves with Youden's index (J) to predict outcome.

RESULTS

EMSI presented shorter DSOZ (15.18 ± 9.06 mm) and DRES (8.56 ± 6.24 mm) compared to ESI (DSOZ: 25.04 ± 16.20 mm, p < 0.009; DRES: 18.88 ± 17.30 mm, p < 0.03) and MSI (DSOZ: 23.37 ± 8.98 mm, p < 0.03; DRES: 15.51 ± 10.11 mm, p < 0.02) for clustering in patients with good outcome. Clustering showed shorter DSOZ and DRES compared to ECD fitting and dSPM (p < 0.05). EMSI had higher performance as outcome predictor (J = 70.63%) compared to ESI (J = 41.27%) and MSI (J = 33.33%) for clustering.

CONCLUSIONS

EMSI provides superior localization and improved predictive performance than individual modalities.

SIGNIFICANCE

EMSI can help the surgical planning and facilitate the localization of epileptogenic foci.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

MEG language mapping using a novel automatic ECD algorithm in comparison with MNE, dSPM, and DICS beamformer

Introduction

The single equivalent current dipole (sECD) is the standard clinical procedure for presurgical language mapping in epilepsy using magnetoencephalography (MEG). However, the sECD approach has not been widely used in clinical assessments, mainly because it requires subjective judgements in selecting several critical parameters. To address this limitation, we developed an automatic sECD algorithm (AsECDa) for language mapping.

Methods

The localization accuracy of the AsECDa was evaluated using synthetic MEG data. Subsequently, the reliability and efficiency of AsECDa were compared to three other common source localization methods using MEG data recorded during two sessions of a receptive language task in 21 epilepsy patients. These methods include minimum norm estimation (MNE), dynamic statistical parametric mapping (dSPM), and dynamic imaging of coherent sources (DICS) beamformer.

Results

For the synthetic single dipole MEG data with a typical signal-to-noise ratio, the average localization error of AsECDa was less than 2 mm for simulated superficial and deep dipoles. For the patient data, AsECDa showed better test-retest reliability (TRR) of the language laterality index (LI) than MNE, dSPM, and DICS beamformer. Specifically, the LI calculated with AsECDa revealed excellent TRR between the two MEG sessions across all patients (Cor = 0.80), while the LI for MNE, dSPM, DICS-event-related desynchronization (ERD) in the alpha band, and DICS-ERD in the low beta band ranged lower (Cor = 0.71, 0.64, 0.54, and 0.48, respectively). Furthermore, AsECDa identified 38% of patients with atypical language lateralization (i.e., right lateralization or bilateral), compared to 73%, 68%, 55%, and 50% identified by DICS-ERD in the low beta band, DICS-ERD in the alpha band, MNE, and dSPM, respectively. Compared to other methods, AsECDa's results were more consistent with previous studies that reported atypical language lateralization in 20-30% of epilepsy patients.

Discussion

Our study suggests that AsECDa is a promising approach for presurgical language mapping, and its fully automated nature makes it easy to implement and reliable for clinical evaluations.

Time-Lagged Multidimensional Pattern Connectivity (TL-MDPC): An EEG/MEG pattern transformation based functional connectivity metric

Functional and effective connectivity methods are essential to study the complex information flow in brain networks underlying human cognition. Only recently have connectivity methods begun to emerge that make use of the full multidimensional information contained in patterns of brain activation, rather than unidimensional summary measures of these patterns. To date, these methods have mostly been applied to fMRI data, and no method allows vertex-to-vertex transformations with the temporal specificity of EEG/MEG data. Here, we introduce time-lagged multidimensional pattern connectivity (TL-MDPC) as a novel bivariate functional connectivity metric for EEG/MEG research. TL-MDPC estimates the vertex-to-vertex transformations among multiple brain regions and across different latency ranges. It determines how well patterns in ROI X at time point tx can linearly predict patterns of ROI Y at time point ty. In the present study, we use simulations to demonstrate TL-MDPC's increased sensitivity to multidimensional effects compared to a unidimensional approach across realistic choices of number of trials and signal-to-noise ratios. We applied TL-MDPC, as well as its unidimensional counterpart, to an existing dataset varying the depth of semantic processing of visually presented words by contrasting a semantic decision and a lexical decision task. TL-MDPC detected significant effects beginning very early on, and showed stronger task modulations than the unidimensional approach, suggesting that it is capable of capturing more information. With TL-MDPC only, we observed rich connectivity between core semantic representation (left and right anterior temporal lobes) and semantic control (inferior frontal gyrus and posterior temporal cortex) areas with greater semantic demands. TL-MDPC is a promising approach to identify multidimensional connectivity patterns, typically missed by unidimensional approaches.

Multiple Dipole Source Position and Orientation Estimation Using Non-Invasive EEG-like Signals

The problem of precisely estimating the position and orientation of multiple dipoles using synthetic EEG signals is considered in this paper. After determining a proper forward model, a nonlinear constrained optimization problem with regularization is solved, and the results are compared with a widely used research code, namely EEGLAB. A thorough sensitivity analysis of the estimation algorithm to the parameters (such as the number of samples and sensors) in the assumed signal measurement model is conducted. To confirm the efficacy of the proposed source identification algorithm on any category of data sets, three different kinds of data-synthetic model data, visually evoked clinical EEG data, and seizure clinical EEG data are used. Furthermore, the algorithm is tested on both the spherical head model and the realistic head model based on the MNI coordinates. The numerical results and comparisons with the EEGLAB show very good agreement, with little pre-processing required for the acquired data.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Towards an objective evaluation of EEG/MEG source estimation methods - The linear approach

The spatial resolution of EEG/MEG source estimates, often described in terms of source leakage in the context of the inverse problem, poses constraints on the inferences that can be drawn from EEG/MEG source estimation results. Software packages for EEG/MEG data analysis offer a large choice of source estimation methods but few tools to experimental researchers for methods evaluation and comparison. Here, we describe a framework and tools for objective and intuitive resolution analysis of EEG/MEG source estimation based on linear systems analysis, and apply those to the most widely used distributed source estimation methods such as L2-minimum-norm estimation (L2-MNE) and linearly constrained minimum variance (LCMV) beamformers. Within this framework it is possible to define resolution metrics that define meaningful aspects of source estimation results (such as localization accuracy in terms of peak localization error, PLE, and spatial extent in terms of spatial deviation, SD) that are relevant to the task at hand and can easily be visualized. At the core of this framework is the resolution matrix, which describes the potential leakage from and into point sources (point-spread and cross-talk functions, or PSFs and CTFs, respectively). Importantly, for linear methods these functions allow generalizations to multiple sources or complex source distributions. This paper provides a tutorial-style introduction into linear EEG/MEG source estimation and resolution analysis aimed at experimental (rather than methods-oriented) researchers. We used this framework to demonstrate how L2-MNE-type as well as LCMV beamforming methods can be evaluated in practice using software tools that have only recently become available for routine use. Our novel methods comparison includes PLE and SD for a larger number of methods than in similar previous studies, such as unweighted, depth-weighted and normalized L2-MNE methods (including dSPM, sLORETA, eLORETA) and two LCMV beamformers. The results demonstrate that some methods can achieve low and even zero PLE for PSFs. However, their SD as well as both PLE and SD for CTFs are far less optimal for all methods, in particular for deep cortical areas. We hope that our paper will encourage EEG/MEG researchers to apply this approach to their own tasks at hand.

Alzheimer's Disease with Epileptiform EEG Activity: Abnormal Cortical Sources of Resting State Delta Rhythms in Patients with Amnesic Mild Cognitive Impairment

BACKGROUND

Patients with amnesic mild cognitive impairment due to Alzheimer's disease (ADMCI) typically show a "slowing" of cortical resting-state eyes-closed electroencephalographic (rsEEG) rhythms. Some of them also show subclinical, non-convulsive, and epileptiform EEG activity (EEA) with an unclear relationship with that "slowing."

OBJECTIVE

Here we tested the hypothesis that the "slowing" of rsEEG rhythms is related to EEA in ADMCI patients.

METHODS

Clinical and instrumental datasets in 62 ADMCI patients and 38 normal elderly (Nold) subjects were available in a national archive. No participant had received a clinical diagnosis of epilepsy. The eLORETA freeware estimated rsEEG cortical sources. The area under the receiver operating characteristic curve (AUROCC) indexed the accuracy of eLORETA solutions in the classification between ADMCI-EEA and ADMCI-noEEA individuals.

RESULTS

EEA was observed in 15% (N = 8) of the ADMCI patients. The ADMCI-EEA group showed: 1) more abnormal Aβ 42 levels in the cerebrospinal fluid as compared to the ADMCI-noEEA group and 2) higher temporal and occipital delta (<4 Hz) rsEEG source activities as compared to the ADMCI-noEEA and Nold groups. Those source activities showed moderate accuracy (AUROCC = 0.70-0.75) in the discrimination between ADMCI-noEEA versus ADMCI-EEA individuals.

CONCLUSION

It can be speculated that in ADMCI-EEA patients, AD-related amyloid neuropathology may be related to an over-excitation in neurophysiological low-frequency (delta) oscillatory mechanisms underpinning cortical arousal and quiet vigilance.

FLUX: A pipeline for MEG analysis

Magnetoencephalography (MEG) allows for quantifying modulations of human neuronal activity on a millisecond time scale while also making it possible to estimate the location of the underlying neuronal sources. The technique relies heavily on signal processing and source modelling. To this end, there are several open-source toolboxes developed by the community. While these toolboxes are powerful as they provide a wealth of options for analyses, the many options also pose a challenge for reproducible research as well as for researchers new to the field. The FLUX pipeline aims to make the analyses steps and setting explicit for standard analysis done in cognitive neuroscience. It focuses on quantifying and source localization of oscillatory brain activity, but it can also be used for event-related fields and multivariate pattern analysis. The pipeline is derived from the Cogitate consortium addressing a set of concrete cognitive neuroscience questions. Specifically, the pipeline including documented code is defined for MNE Python (a Python toolbox) and FieldTrip (a Matlab toolbox), and a data set on visuospatial attention is used to illustrate the steps. The scripts are provided as notebooks implemented in Jupyter Notebook and MATLAB Live Editor providing explanations, justifications and graphical outputs for the essential steps. Furthermore, we also provide suggestions for text and parameter settings to be used in registrations and publications to improve replicability and facilitate pre-registrations. The FLUX can be used for education either in self-studies or guided workshops. We expect that the FLUX pipeline will strengthen the field of MEG by providing some standardization on the basic analysis steps and by aligning approaches across toolboxes. Furthermore, we also aim to support new researchers entering the field by providing education and training. The FLUX pipeline is not meant to be static; it will evolve with the development of the toolboxes and with new insights. Furthermore, with the anticipated increase in MEG systems based on the Optically Pumped Magnetometers, the pipeline will also evolve to embrace these developments.

Informed MEG/EEG source imaging reveals the locations of interictal spikes missed by SEEG

Determining the accurate locations of interictal spikes has been fundamental in the presurgical evaluation of epilepsy surgery. Stereo-electroencephalography (SEEG) is able to directly record cortical activity and localize interictal spikes. However, the main caveat of SEEG techniques is that they have limited spatial sampling (covering <5% of the whole brain), which may lead to missed spikes originating from brain regions that were not covered by SEEG. To address this problem, we propose a SEEG-informed minimum-norm estimates (SIMNE) method by combining SEEG with magnetoencephalography (MEG) or EEG. Specifically, the spike locations determined by SEEG offer as a priori information to guide MEG source reconstruction. Both computer simulations and experiments using data from five epilepsy patients were conducted to evaluate the performance of SIMNE. Our results demonstrate that SIMNE generates more accurate source estimation than a traditional minimum-norm estimates method and reveals the locations of spikes missed by SEEG, which would improve presurgical evaluation of the epileptogenic zone.

Imaging somatosensory cortex responses measured by OPM-MEG: Variational free energy-based spatial smoothing estimation approach

In recent years, optically pumped magnetometer (OPM)-based magnetoencephalography (MEG) has shown potential for analyzing brain activity. It has a flexible sensor configuration and comparable sensitivity to conventional SQUID-MEG. We constructed a 32-channel OPM-MEG system and used it to measure cortical responses to median and ulnar nerve stimulations. Traditional magnetic source imaging methods tend to blur the spatial extent of sources. Accurate estimation of the spatial size of the source is important for studying the organization of brain somatotopy and for pre-surgical functional mapping. We proposed a new method called variational free energy-based spatial smoothing estimation (FESSE) to enhance the accuracy of mapping somatosensory cortex responses. A series of computer simulations based on the OPM-MEG showed better performance than the three types of competing methods under different levels of signal-to-noise ratios, source patch sizes, and co-registration errors. FESSE was then applied to the source imaging of the OPM-MEG experimental data.

Task modulation of spatiotemporal dynamics in semantic brain networks: An EEG/MEG study

How does brain activity in distributed semantic brain networks evolve over time, and how do these regions interact to retrieve the meaning of words? We compared spatiotemporal brain dynamics between visual lexical and semantic decision tasks (LD and SD), analysing whole-cortex evoked responses and spectral functional connectivity (coherence) in source-estimated electroencephalography and magnetoencephalography (EEG and MEG) recordings. Our evoked analysis revealed generally larger activation for SD compared to LD, starting in primary visual area (PVA) and angular gyrus (AG), followed by left posterior temporal cortex (PTC) and left anterior temporal lobe (ATL). The earliest activation effects in ATL were significantly left-lateralised. Our functional connectivity results showed significant connectivity between left and right ATL, PTC and right ATL in an early time window, as well as between left ATL and IFG in a later time window. The connectivity of AG was comparatively sparse. We quantified the limited spatial resolution of our source estimates via a leakage index for careful interpretation of our results. Our findings suggest that the different demands on semantic information retrieval in lexical and semantic decision tasks first modulate visual and attentional processes, then multimodal semantic information retrieval in the ATLs and finally control regions (PTC and IFG) in order to extract task-relevant semantic features for response selection. Whilst our evoked analysis suggests a dominance of left ATL for semantic processing, our functional connectivity analysis also revealed significant involvement of right ATL in the more demanding semantic task. Our findings demonstrate the complementarity of evoked and functional connectivity analysis, as well as the importance of dynamic information for both types of analyses.

A Long Short-Term Memory Network for Sparse Spatiotemporal EEG Source Imaging

EEG inverse problem is underdetermined, which poses a long standing challenge in Neuroimaging. The combination of source-imaging and analysis of cortical directional networks enables us to noninvasively explore the underlying neural processes. However, existing EEG source imaging approaches mainly focus on performing the direct inverse operation for source estimation, which will be inevitably influenced by noise and the strategy used to find the inverse solution. Here, we develop a new source imaging technique, Deep Brain Neural Network (DeepBraiNNet), for robust sparse spatiotemporal EEG source estimation. In DeepBraiNNet, considering that Recurrent Neural Network (RNN) are usually "deep" in temporal dimension and thus suitable for time sequence modelling, the RNN with Long Short-Term Memory (LSTM) is utilized to approximate the inverse operation for the lead field matrix instead of performing the direct inverse operation, which avoids the possible effect of the direct inverse operation on the underdetermined lead field matrix prone to be influenced by noise. Simulations on various source patterns and noise conditions confirmed that the proposed approach could actually recover the spatiotemporal sources well, outperforming existing state of-the-art methods. DeepBraiNNet also estimated sparse MI related activation patterns when it was applied to a real Motor Imagery dataset, consistent with other findings based on EEG and fMRI. Based on the spatiotemporal sources estimated from DeepBraiNNet, we constructed MI related cortical neural networks, which clearly exhibited strong contralateral network patterns for the two MI tasks. Consequently, DeepBraiNNet may provide an alternative way different from the conventional approaches for spatiotemporal EEG source imaging.

Imaging the extent and location of spatiotemporally distributed epileptiform sources from MEG measurements

Non-invasive MEG/EEG source imaging provides valuable information about the epileptogenic brain areas which can be used to aid presurgical planning in focal epilepsy patients suffering from drug-resistant seizures. However, the source extent estimation for electrophysiological source imaging remains to be a challenge and is usually largely dependent on subjective choice. Our recently developed algorithm, fast spatiotemporal iteratively reweighted edge sparsity minimization (FAST-IRES) strategy, has been shown to objectively estimate extended sources from EEG recording, while it has not been applied to MEG recordings. In this work, through extensive numerical experiments and real data analysis in a group of focal drug-resistant epilepsy patients' interictal spikes, we demonstrated the ability of FAST-IRES algorithm to image the location and extent of underlying epilepsy sources from MEG measurements. Our results indicate the merits of FAST-IRES in imaging the location and extent of epilepsy sources for pre-surgical evaluation from MEG measurements.

Resting State Alpha Electroencephalographic Rhythms Are Affected by Sex in Cognitively Unimpaired Seniors and Patients with Alzheimer's Disease and Amnesic Mild Cognitive Impairment: A Retrospective and Exploratory Study

In the present retrospective and exploratory study, we tested the hypothesis that sex may affect cortical sources of resting state eyes-closed electroencephalographic (rsEEG) rhythms recorded in normal elderly (Nold) seniors and patients with Alzheimer's disease and mild cognitive impairment (ADMCI). Datasets in 69 ADMCI and 57 Nold individuals were taken from an international archive. The rsEEG rhythms were investigated at individual delta, theta, and alpha frequency bands and fixed beta (14-30 Hz) and gamma (30-40 Hz) bands. Each group was stratified into matched females and males. The sex factor affected the magnitude of rsEEG source activities in the Nold seniors. Compared with the males, the females were characterized by greater alpha source activities in all cortical regions. Similarly, the parietal, temporal, and occipital alpha source activities were greater in the ADMCI-females than the males. Notably, the present sex effects did not depend on core genetic (APOE4), neuropathological (Aβ42/phospho-tau ratio in the cerebrospinal fluid), structural neurodegenerative and cerebrovascular (MRI) variables characterizing sporadic AD-related processes in ADMCI seniors. These results suggest the sex factor may significantly affect neurophysiological brain neural oscillatory synchronization mechanisms underpinning the generation of dominant rsEEG alpha rhythms to regulate cortical arousal during quiet vigilance.

A systematic evaluation of source reconstruction of resting MEG of the human brain with a new high-resolution atlas: Performance, precision, and parcellation

Noninvasive functional neuroimaging of the human brain can give crucial insight into the mechanisms that underpin healthy cognition and neurological disorders. Magnetoencephalography (MEG) measures extracranial magnetic fields originating from neuronal activity with high temporal resolution, but requires source reconstruction to make neuroanatomical inferences from these signals. Many source reconstruction algorithms are available, and have been widely evaluated in the context of localizing task-evoked activities. However, no consensus yet exists on the optimum algorithm for resting-state data. Here, we evaluated the performance of six commonly-used source reconstruction algorithms based on minimum-norm and beamforming estimates. Using human resting-state MEG, we compared the algorithms using quantitative metrics, including resolution properties of inverse solutions and explained variance in sensor-level data. Next, we proposed a data-driven approach to reduce the atlas from the Human Connectome Project's multi-modal parcellation of the human cortex based on metrics such as MEG signal-to-noise-ratio and resting-state functional connectivity gradients. This procedure produced a reduced cortical atlas with 230 regions, optimized to match the spatial resolution and the rank of MEG data from the current generation of MEG scanners. Our results show that there is no "one size fits all" algorithm, and make recommendations on the appropriate algorithms depending on the data and aimed analyses. Our comprehensive comparisons and recommendations can serve as a guide for choosing appropriate methodologies in future studies of resting-state MEG.

Inflection across Categories: Tracking Abstract Morphological Processing in Language Production with MEG

Coherent language production requires that speakers adapt words to their grammatical contexts. A fundamental challenge in establishing a functional delineation of this process in the brain is that each linguistic process tends to correlate with numerous others. Our work investigated the neural basis of morphological inflection by measuring magnetoencephalography during the planning of inflected and uninflected utterances that varied across several linguistic dimensions. Results reveal increased activity in the left lateral frontotemporal cortex when inflection is planned, irrespective of phonological specification, syntactic context, or semantic type. Additional findings from univariate and connectivity analyses suggest that the brain distinguishes between different types of inflection. Specifically, planning noun and verb utterances requiring the addition of the suffix -s elicited increased activity in the ventral prefrontal cortex. A broadly distributed effect of syntactic context (verb vs. noun) was also identified. Results from representational similarity analysis indicate that this effect cannot be explained in terms of word meaning. Together, these results 1) offer evidence for a neural representation of abstract inflection that separates from other stimulus properties and 2) challenge theories that emphasize semantic content as a source of verb/noun processing differences.

Testing covariance models for MEG source reconstruction of hippocampal activity

Beamforming is one of the most commonly used source reconstruction methods for magneto- and electroencephalography (M/EEG). One underlying assumption, however, is that distant sources are uncorrelated and here we tested whether this is an appropriate model for the human hippocampal data. We revised the Empirical Bayesian Beamfomer (EBB) to accommodate specific a-priori correlated source models. We showed in simulation that we could use model evidence (as approximated by Free Energy) to distinguish between different correlated and uncorrelated source scenarios. Using group MEG data in which the participants performed a hippocampal-dependent task, we explored the possibility that the hippocampus or the cortex or both were correlated in their activity across hemispheres. We found that incorporating a correlated hippocampal source model significantly improved model evidence. Our findings help to explain why, up until now, the majority of MEG-reported hippocampal activity (typically making use of beamformers) has been estimated as unilateral.

Alpha-band cortico-cortical phase synchronization is associated with effective connectivity in the motor network

OBJECTIVE

Communication-through-coherence proposes that the phase synchronization (PS) of neural oscillations between cortical areas supports neural communication. In this study, we exploited transcranial magnetic stimulation (TMS)-evoked potentials (TEPs) to test this hypothesis at the macroscale level, i.e., whether PS between cortical areas supports interarea communication. TEPs are electroencephalographic (EEG) responses time-locked to TMS pulses reflecting interarea communication, as they are generated by the transmission of neural activity from the stimulated area to connected regions. If interarea PS is important for communication, it should be associated with the TEP amplitude in the connected areas.

METHODS

TMS was delivered over the left primary motor cortex (M1) of fourteen healthy volunteers, and 70-channel EEG was recorded. Early TEP components were source-localized to identify their generators, i.e., distant brain regions activated by M1 through effective connections. Next, linear regressions were used to test the relationship between the TEP amplitude and the pre-stimulus PS between the M1 and the connected regions in four frequency bands (range 4-45 Hz).

RESULTS

Pre-stimulus interarea PS in the alpha-band was positively associated with the amplitude of early TEP components, namely, the N15 (ipsilateral supplementary motor area), P25 (contralateral M1) and P60 (ipsilateral parietal cortex).

CONCLUSIONS

Alpha-band PS predicts the response amplitude of the distant brain regions effectively connected to M1.

SIGNIFICANCE

Our study supports the role of EEG-PS in interarea communication, as theorized by communication-through-coherence.

From letters to composed concepts: A magnetoencephalography study of reading

Language comprehension requires the recognition of individual words and the combination of their meanings to yield complex concepts or interpretations. This combinatory process often requires the insertion of unstated semantic material between words, based on thematic or feature knowledge. For example, the phrase horse barn is not interpreted as a blend of a horse and a barn, but specifically a barn where horses are kept. Previous neuroscientific evidence suggests that left posterior and anterior temporal cortex underpin thematic and feature‐based concept knowledge, respectively, but much remains unclear about how these areas contribute to combinatory language processing. Using magnetoencephalography, we contrasted source‐localized responses to modifier‐noun phrases involving thematic relations versus feature modifications, while also examining how lower‐level orthographic processing fed composition. Participants completed three procedures examining responses to letter‐strings, adjective‐noun phrases, and noun–noun combinations that varied the semantic relations between words. We found that sections of the left anterior temporal lobe, posterior temporal lobe, and cortex surrounding the angular gyrus were all engaged in the minimal composition of adjective‐noun phrases, a more distributed network than in most prior studies of minimal composition. Of these regions, only the left posterior temporal lobe was additionally sensitive to implicit thematic relations between composing words, suggesting that it houses a specialized relational processing component in a wider composition network. We additionally identified a left occipitotemporal progression from orthographic to lexical processing, feeding ventral anterior areas engaged in the combination of word meanings. Finally, by examining source signal leakage, we characterized the degree to which these responses could be distinguished from one another using source estimation.

Resting State Alpha Electroencephalographic Rhythms Are Differently Related to Aging in Cognitively Unimpaired Seniors and Patients with Alzheimer's Disease and Amnesic Mild Cognitive Impairment

BACKGROUND

In relaxed adults, staying in quiet wakefulness at eyes closed is related to the so-called resting state electroencephalographic (rsEEG) rhythms, showing the highest amplitude in posterior areas at alpha frequencies (8-13 Hz).

OBJECTIVE

Here we tested the hypothesis that age may affect rsEEG alpha (8-12 Hz) rhythms recorded in normal elderly (Nold) seniors and patients with mild cognitive impairment due to Alzheimer's disease (ADMCI).

METHODS

Clinical and rsEEG datasets in 63 ADMCI and 60 Nold individuals (matched for demography, education, and gender) were taken from an international archive. The rsEEG rhythms were investigated at individual delta, theta, and alpha frequency bands, as well as fixed beta (14-30 Hz) and gamma (30-40 Hz) bands. Each group was stratified into three subgroups based on age ranges (i.e., tertiles).

RESULTS

As compared to the younger Nold subgroups, the older one showed greater reductions in the rsEEG alpha rhythms with major topographical effects in posterior regions. On the contrary, in relation to the younger ADMCI subgroups, the older one displayed a lesser reduction in those rhythms. Notably, the ADMCI subgroups pointed to similar cerebrospinal fluid AD diagnostic biomarkers, gray and white matter brain lesions revealed by neuroimaging, and clinical and neuropsychological scores.

CONCLUSION

The present results suggest that age may represent a deranging factor for dominant rsEEG alpha rhythms in Nold seniors, while rsEEG alpha rhythms in ADMCI patients may be more affected by the disease variants related to earlier versus later onset of the AD.

A novel beamformer-based imaging of phase-amplitude coupling (BIPAC) unveiling the inter-regional connectivity of emotional prosody processing in women with primary dysmenorrhea

Objective. Neural communication or the interactions of brain regions play a key role in the formation of functional neural networks. A type of neural communication can be measured in the form of phase-amplitude coupling (PAC), which is the coupling between the phase of low-frequency oscillations and the amplitude of high-frequency oscillations. This paper presents a beamformer-based imaging method, beamformer-based imaging of PAC (BIPAC), to quantify the strength of PAC between a seed region and other brain regions.Approach. A dipole is used to model the ensemble of neural activity within a group of nearby neurons and represents a mixture of multiple source components of cortical activity. From ensemble activity at each brain location, the source component with the strongest coupling to the seed activity is extracted, while unrelated components are suppressed to enhance the sensitivity of coupled-source estimation.Main results. In evaluations using simulation data sets, BIPAC proved advantageous with regard to estimation accuracy in source localization, orientation, and coupling strength. BIPAC was also applied to the analysis of magnetoencephalographic signals recorded from women with primary dysmenorrhea in an implicit emotional prosody experiment. In response to negative emotional prosody, auditory areas revealed strong PAC with the ventral auditory stream and occipitoparietal areas in the theta-gamma and alpha-gamma bands, which may respectively indicate the recruitment of auditory sensory memory and attention reorientation. Moreover, patients with more severe pain experience appeared to have stronger coupling between auditory areas and temporoparietal regions.Significance. Our findings indicate that the implicit processing of emotional prosody is altered by menstrual pain experience. The proposed BIPAC is feasible and applicable to imaging inter-regional connectivity based on cross-frequency coupling estimates. The experimental results also demonstrate that BIPAC is capable of revealing autonomous brain processing and neurodynamics, which are more subtle than active and attended task-driven processing.

Abnormalities of Cortical Sources of Resting State Alpha Electroencephalographic Rhythms are Related to Education Attainment in Cognitively Unimpaired Seniors and Patients with Alzheimer's Disease and Amnesic Mild Cognitive Impairment

In normal old (Nold) and Alzheimer's disease (AD) persons, a high cognitive reserve (CR) makes them more resistant and resilient to brain neuropathology and neurodegeneration. Here, we tested whether these effects may affect neurophysiological oscillatory mechanisms generating dominant resting state electroencephalographic (rsEEG) alpha rhythms in Nold and patients with mild cognitive impairment (MCI) due to AD (ADMCI). Data in 60 Nold and 70 ADMCI participants, stratified in higher (Edu+) and lower (Edu-) educational attainment subgroups, were available in an Italian-Turkish archive. The subgroups were matched for age, gender, and education. RsEEG cortical sources were estimated by eLORETA freeware. As compared to the Nold-Edu- subgroup, the Nold-Edu+ subgroup showed greater alpha source activations topographically widespread. On the contrary, in relation to the ADMCI-Edu- subgroup, the ADMCI-Edu+ subgroup displayed lower alpha source activations topographically widespread. Furthermore, the 2 ADMCI subgroups had matched cerebrospinal AD diagnostic biomarkers, brain gray-white matter measures, and neuropsychological scores. The current findings suggest that a high CR may be related to changes in rsEEG alpha rhythms in Nold and ADMCI persons. These changes may underlie neuroprotective effects in Nold seniors and subtend functional compensatory mechanisms unrelated to brain structure alterations in ADMCI patients.

The "when" and "where" of the interplay between attentional capture and response inhibition during a Go/NoGo variant

Inhibitory control relies on attention, inhibition, and other functions that are integrated across neural networks in an interactive manner. Functional MRI studies have provided excellent spatial mapping of the involved regions. However, finer temporal resolution is needed to capture the underlying neural dynamics and the pattern of their functional contributions. Here, we used anatomically-constrained magnetoencephalography (aMEG) which combines MEG with structural MRI to examine how the spatial ("where") and temporal ("when") processing stages and interregional co-oscillations unfold in real time to contribute to inhibitory control. Healthy participants completed a modified Go/NoGo paradigm in which a subset of stimuli was modified to be visually salient (SAL). Compared to the non-modified condition, the SAL manipulation facilitated response withholding on NoGo trials and hindered responding to Go stimuli, reflecting attentional capture effectuated by an orienting response to SAL stimuli. aMEG source estimates indicate SAL stimuli elicited the attentional "circuit breaker" effect through early activity within a right-lateralized network centered around the lateral temporal cortex with additional activity in the pre-supplementary motor area (preSMA) and anterior insula (aINS/FO). Activity of the bilateral inferior frontal cortex responded specifically to inhibitory demands and was generally unaffected by the attentional manipulation. In contrast, early aINS/FO activity was sensitive to stimulus salience while subsequent activity was specific to inhibitory control. Activity estimated to the medial prefrontal cortex including the dorsal anterior cingulate cortex and preSMA reflected an integrative role that was sensitive to both inhibitory and attentional stimulus properties. At the level of neurofunctional networks, neural synchrony in the theta band (4-7 Hz) revealed interactions between principal cortical regions subserving attentional and inhibitory processes. Together, these results underscore the dynamic, integrative processing stages underlying inhibitory control.

Spatial fidelity of MEG/EEG source estimates: A general evaluation approach

Low spatial resolution is often cited as the most critical limitation of magneto- and electroencephalography (MEG and EEG), but a unifying framework for quantifying the spatial fidelity of M/EEG source estimates has yet to be established; previous studies have focused on linear estimation methods under ideal scenarios without noise. Here we present an approach that quantifies the spatial fidelity of M/EEG estimates from simulated patch activations over the entire neocortex superposed on measured resting-state data. This approach grants more generalizability in the evaluation process that allows for, e.g., comparing linear and non-linear estimates in the whole brain for different signal-to-noise ratios (SNR), number of active sources and activation waveforms. Using this framework, we evaluated the MNE, dSPM, sLORETA, eLORETA, and MxNE methods and found that the spatial fidelity varies significantly with SNR, following a largely sigmoidal curve whose shape varies depending on which aspect of spatial fidelity that is being quantified and the source estimation method. We believe that these methods and results will be useful when interpreting M/EEG source estimates as well as in methods development.

Abnormalities of Cortical Sources of Resting State Delta Electroencephalographic Rhythms Are Related to Epileptiform Activity in Patients With Amnesic Mild Cognitive Impairment Not Due to Alzheimer's Disease

In the present exploratory and retrospective study, we hypothesized that cortical sources of resting state eyes-closed electroencephalographic (rsEEG) rhythms might be more abnormal in patients with epileptiform EEG activity (spike-sharp wave discharges, giant spikes) and amnesic mild cognitive impairment not due to Alzheimer's disease (noADMCI-EEA) than matched noADMCI patients without EEA (noADMCI-noEEA). Clinical, neuroimaging, neuropsychological, and rsEEG data in 32 noADMCI and 30 normal elderly (Nold) subjects were available in a national archive. Age, gender, and education were carefully matched among them. No subject had received a clinical diagnosis of epilepsy. Individual alpha frequency peak (IAF) was used to determine the delta, theta, and alpha frequency bands of rsEEG rhythms. Fixed beta and gamma bands were also considered. Regional rsEEG cortical sources were estimated by eLORETA freeware. Area under receiver operating characteristic (AUROC) curves indexed the accuracy of eLORETA solutions in the classification between noADMCI-EEA and noADMCI-noEEA individuals. As novel findings, EEA was observed in 41% of noADMCI patients. Furthermore, these noADMCI-EEA patients showed higher temporal delta source activities as compared to noADMCI-no EEA patients and Nold subjects. Those activities discriminated individuals of the two NoADMCI groups with an accuracy of about 70%. The significant percentage of noADMCI-EEA patients showing EEA and marked abnormalities in temporal rsEEG rhythms at delta frequencies suggest a substantial role of underlying neural hypersynchronization mechanisms in their brain dysfunctions.

Stacked autoencoders as new models for an accurate Alzheimer's disease classification support using resting-state EEG and MRI measurements

OBJECTIVE

This retrospective and exploratory study tested the accuracy of artificial neural networks (ANNs) at detecting Alzheimer's disease patients with dementia (ADD) based on input variables extracted from resting-state electroencephalogram (rsEEG), structural magnetic resonance imaging (sMRI) or both.

METHODS

For the classification exercise, the ANNs had two architectures that included stacked (autoencoding) hidden layers recreating input data in the output. The classification was based on LORETA source estimates from rsEEG activity recorded with 10-20 montage system (19 electrodes) and standard sMRI variables in 89 ADD and 45 healthy control participants taken from a national database.

RESULTS

The ANN with stacked autoencoders and a deep leaning model representing both ADD and control participants showed classification accuracies in discriminating them of 80%, 85%, and 89% using rsEEG, sMRI, and rsEEG + sMRI features, respectively. The two ANNs with stacked autoencoders and a deep leaning model specialized for either ADD or control participants showed classification accuracies of 77%, 83%, and 86% using the same input features.

CONCLUSIONS

The two architectures of ANNs using stacked (autoencoding) hidden layers consistently reached moderate to high accuracy in the discrimination between ADD and healthy control participants as a function of the rsEEG and sMRI features employed.

SIGNIFICANCE

The present results encourage future multi-centric, prospective and longitudinal cross-validation studies using high resolution EEG techniques and harmonized clinical procedures towards clinical applications of the present ANNs.

Electrophysiological evidence of phonemotopic representations of vowels in the primary and secondary auditory cortex

How the brain encodes the speech acoustic signal into phonological representations is a fundamental question for the neurobiology of language. Determining whether this process is characterized by tonotopic maps in primary or secondary auditory areas, with bilateral or leftward activity, remains a long-standing challenge. Magnetoencephalographic studies failed to show hierarchical and asymmetric hints for speech processing. We employed high-density electroencephalography to map the Salento Italian vowel system onto cortical sources using the N1 auditory evoked component. We found evidence that the N1 is characterized by hierarchical and asymmetrical indexes in primary and secondary auditory areas structuring vowel representations. Importantly, the N1 was characterized by early and late phases. The early N1 peaked at 125-135 msec and was localized in the primary auditory cortex; the late N1 peaked at 145-155 msec and was localized in the left superior temporal gyrus. We showed that early in the primary auditory cortex, the cortical spatial arrangements-along the lateral-medial and anterior-posterior gradients-are broadly warped by phonemotopic patterns according to the distinctive feature principle. These phonemotopic patterns are carefully refined in the superior temporal gyrus along the inferior-superior and anterior-posterior gradients. The dynamical and hierarchical interface between primary and secondary auditory areas and the interaction effects between Height and Place features generate the categorical representation of the Salento Italian vowels.

Wearable neuroimaging: Combining and contrasting magnetoencephalography and electroencephalography

One of the most severe limitations of functional neuroimaging techniques, such as magnetoencephalography (MEG), is that participants must maintain a fixed head position during data acquisition. This imposes restrictions on the characteristics of the experimental cohorts that can be scanned and the experimental questions that can be addressed. For these reasons, the use of 'wearable' neuroimaging, in which participants can move freely during scanning, is attractive. The most successful example of wearable neuroimaging is electroencephalography (EEG), which employs lightweight and flexible instrumentation that makes it useable in almost any experimental setting. However, EEG has major technical limitations compared to MEG, and therefore the development of wearable MEG, or hybrid MEG/EEG systems, is a compelling prospect. In this paper, we combine and compare EEG and MEG measurements, the latter made using a new generation of optically-pumped magnetometers (OPMs). We show that these new second generation commercial OPMs, can be mounted on the scalp in an 'EEG-like' cap, enabling the acquisition of high fidelity electrophysiological measurements. We show that these sensors can be used in conjunction with conventional EEG electrodes, offering the potential for the development of hybrid MEG/EEG systems. We compare concurrently measured signals, showing that, whilst both modalities offer high quality data in stationary subjects, OPM-MEG measurements are less sensitive to artefacts produced when subjects move. Finally, we show using simulations that OPM-MEG offers a fundamentally better spatial specificity than EEG. The demonstrated technology holds the potential to revolutionise the utility of functional brain imaging, exploiting the flexibility of wearable systems to facilitate hitherto impractical experimental paradigms.

Simultaneous EEG and fMRI reveals stronger sensitivity to orthographic strings in the left occipito-temporal cortex of typical versus poor beginning readers

The level of reading skills in children and adults is reflected in the strength of preferential neural activation to print. Such preferential activation appears in the N1 event-related potential (ERP) over the occipitotemporal scalp after around 150–250 ms and the corresponding blood oxygen level dependent (BOLD) signal in the ventral occipitotemporal (vOT) cortex. Here, orthography-sensitive (print vs. false font) processing was examined using simultaneous EEG-fMRI in 38 first grade children with poor and typical reading skills, and at varying familial risk for developmental dyslexia. Coarse orthographic sensitivity was observed as an increased activation to print in the N1 ERP and in the BOLD signal of individually varying vOT regions in 57% of beginning readers. Finer differentiation in processing orthographic strings (words vs. nonwords) further occurred in specific vOT clusters. Neither method alone showed robust differences in orthography-sensitive processing between typical and poor reading children. Importantly, using single-trial N1 ERP-informed fMRI analysis, we found differential modulation of the orthography-sensitive BOLD response in the left vOT for typical readers only. This result, thus, confirms subtle functional alterations in a brain structure known to be critical for fluent reading at the very beginning of reading instruction.

The effect of stimulation type, head modeling, and combined EEG and MEG on the source reconstruction of the somatosensory P20/N20 component

Modeling and experimental parameters influence the Electro- (EEG) and Magnetoencephalography (MEG) source analysis of the somatosensory P20/N20 component. In a sensitivity group study, we compare P20/N20 source analysis due to different stimulation type (Electric-Wrist [EW], Braille-Tactile [BT], or Pneumato-Tactile [PT]), measurement modality (combined EEG/MEG - EMEG, EEG, or MEG) and head model (standard or individually skull-conductivity calibrated including brain anisotropic conductivity). Considerable differences between pairs of stimulation types occurred (EW-BT: 8.7 ± 3.3 mm/27.1° ± 16.4°, BT-PT: 9 ± 5 mm/29.9° ± 17.3°, and EW-PT: 9.8 ± 7.4 mm/15.9° ± 16.5° and 75% strength reduction of BT or PT when compared to EW) regardless of the head model used. EMEG has nearly no localization differences to MEG, but large ones to EEG (16.1 ± 4.9 mm), while source orientation differences are non-negligible to both EEG (14° ± 3.7°) and MEG (12.5° ± 10.9°). Our calibration results show a considerable inter-subject variability (3.1-14 mS/m) for skull conductivity. The comparison due to different head model show localization differences smaller for EMEG (EW: 3.4 ± 2.4 mm, BT: 3.7 ± 3.4 mm, and PT: 5.9 ± 6.8 mm) than for EEG (EW: 8.6 ± 8.3 mm, BT: 11.8 ± 6.2 mm, and PT: 10.5 ± 5.3 mm), while source orientation differences for EMEG (EW: 15.4° ± 6.3°, BT: 25.7° ± 15.2° and PT: 14° ± 11.5°) and EEG (EW: 14.6° ± 9.5°, BT: 16.3° ± 11.1° and PT: 12.9° ± 8.9°) are in the same range. Our results show that stimulation type, modality and head modeling all have a non-negligible influence on the source reconstruction of the P20/N20 component. The complementary information of both modalities in EMEG can be exploited on the basis of detailed and individualized head models.

Are Functional Brain Alterations Present in Low Back Pain? A Systematic Review of EEG Studies

This systematic review analyzed available literature on functional brain alterations in low back pain (LBP) measured with electroencephalography (EEG), as until now evidence thereof was unclear. Four electronic databases were systematically searched the 10th of March 2018, resulting in 12 included studies. Studies showed a risk of bias of 37.5 to 75% using the Newcastle-Ottawa Scale for case-control studies. Limited evidence reported higher amplitudes of balance-related potentials and early components of somatosensory evoked potentials (SEP) to noxious stimuli, and altered feedback-related negativity and P300 potentials during decision-making in chronic LBP (CLBP). These findings suggest postural strategies requiring a higher cortical attention-demand, increased sensory-discriminative processing of noxious input, and altered decision-making in CLBP. However, further research is warranted as these inferences were based on single studies. Moderate evidence for unaltered amplitude of late-phase SEPs to noxious stimuli and auditory evoked potentials in LBP implies that the affective-emotional processing of stimuli might be unaffected in LBP. Furthermore, moderate evidence indicated disturbed habituation of somatosensory stimuli in LBP. Most studies examined nonspecific or mixed CLBP populations, hence EEG-quantified brain activity in (sub)acute or recurrent LBP still needs to be explored. PERSPECTIVE: This review presents an overview of the current understanding of the functional LBP brain measured with EEG. The limited evidence in current research suggests altered cortical function regarding balance control, somatosensory processing, and decision making in LBP, and highlights opportunities for future EEG-research.

A Primer on Electroencephalography and Event-Related Potentials for Organizational Neuroscience

Electroencephalography (EEG) was the first of the noninvasive brain measures in neuroscience. Technical advances over the last 100 years or so have rendered EEG a true brain imaging technique. Here, we provide an accessible primer on the biophysics of EEG, on measurement aspects, and on the analysis of EEG data. We use the example of event-related potentials (ERPs), although the issues apply equally to other varieties of EEG signals, and provide an overview of analytic methods at the base of the so-called electrical neuroimaging framework. We detail the interpretational strengths of electrical neuroimaging for organizational researchers and describe some domains of ongoing technical developments. We likewise emphasize practical considerations with the use of EEG in more real-world settings. This primer is intended to provide organizational researchers specifically, and novices more generally, an access point to understanding how EEG may be applied in their research.

Computationally optimized ECoG stimulation with local safety constraints

Direct stimulation of the cortical surface is used clinically for cortical mapping and modulation of local activity. Future applications of cortical modulation and brain-computer interfaces may also use cortical stimulation methods. One common method to deliver current is through electrocorticography (ECoG) stimulation in which a dense array of electrodes are placed subdurally or epidurally to stimulate the cortex. However, proximity to cortical tissue limits the amount of current that can be delivered safely. It may be desirable to deliver higher current to a specific local region of interest (ROI) while limiting current to other local areas more stringently than is guaranteed by global safety limits. Two commonly used global safety constraints bound the total injected current and individual electrode currents. However, these two sets of constraints may not be sufficient to prevent high current density locally (hot-spots). In this work, we propose an efficient approach that prevents current density hot-spots in the entire brain while optimizing ECoG stimulus patterns for targeted stimulation. Specifically, we maximize the current along a particular desired directional field in the ROI while respecting three safety constraints: one on the total injected current, one on individual electrode currents, and the third on the local current density magnitude in the brain. This third set of constraints creates a computational barrier due to the huge number of constraints needed to bound the current density at every point in the entire brain. We overcome this barrier by adopting an efficient two-step approach. In the first step, the proposed method identifies the safe brain region, which cannot contain any hot-spots solely based on the global bounds on total injected current and individual electrode currents. In the second step, the proposed algorithm iteratively adjusts the stimulus pattern to arrive at a solution that exhibits no hot-spots in the remaining brain. We report on simulations on a realistic finite element (FE) head model with five anatomical ROIs and two desired directional fields. We also report on the effect of ROI depth and desired directional field on the focality of the stimulation. Finally, we provide an analysis of optimization runtime as a function of different safety and modeling parameters. Our results suggest that optimized stimulus patterns tend to differ from those used in clinical practice.

Ghost interactions in MEG/EEG source space: A note of caution on inter-areal coupling measures

When combined with source modeling, magneto- (MEG) and electroencephalography (EEG) can be used to study long-range interactions among cortical processes non-invasively. Estimation of such inter-areal connectivity is nevertheless hindered by instantaneous field spread and volume conduction, which artificially introduce linear correlations and impair source separability in cortical current estimates. To overcome the inflating effects of linear source mixing inherent to standard interaction measures, alternative phase- and amplitude-correlation based connectivity measures, such as imaginary coherence and orthogonalized amplitude correlation have been proposed. Being by definition insensitive to zero-lag correlations, these techniques have become increasingly popular in the identification of correlations that cannot be attributed to field spread or volume conduction. We show here, however, that while these measures are immune to the direct effects of linear mixing, they may still reveal large numbers of spurious false positive connections through field spread in the vicinity of true interactions. This fundamental problem affects both region-of-interest-based analyses and all-to-all connectome mappings. Most importantly, beyond defining and illustrating the problem of spurious, or "ghost" interactions, we provide a rigorous quantification of this effect through extensive simulations. Additionally, we further show that signal mixing also significantly limits the separability of neuronal phase and amplitude correlations. We conclude that spurious correlations must be carefully considered in connectivity analyses in MEG/EEG source space even when using measures that are immune to zero-lag correlations.

Functional cortical source connectivity of resting state electroencephalographic alpha rhythms shows similar abnormalities in patients with mild cognitive impairment due to Alzheimer's and Parkinson's diseases

OBJECTIVE

This study tested the hypothesis that markers of functional cortical source connectivity of resting state eyes-closed electroencephalographic (rsEEG) rhythms may be abnormal in subjects with mild cognitive impairment due to Alzheimer's (ADMCI) and Parkinson's (PDMCI) diseases compared to healthy elderly subjects (Nold).

METHODS

rsEEG data had been collected in ADMCI, PDMCI, and Nold subjects (N = 75 for any group). eLORETA freeware estimated functional lagged linear connectivity (LLC) from rsEEG cortical sources. Area under receiver operating characteristic (AUROC) curve indexed the accuracy in the classification of Nold and MCI individuals.

RESULTS

Posterior interhemispheric and widespread intrahemispheric alpha LLC solutions were abnormally lower in both MCI groups compared to the Nold group. At the individual level, AUROC curves of LLC solutions in posterior alpha sources exhibited moderate accuracies (0.70-0.72) in the discrimination of Nold vs. ADMCI-PDMCI individuals. No differences in the LLC solutions were found between the two MCI groups.

CONCLUSIONS

These findings unveil similar abnormalities in functional cortical connectivity estimated in widespread alpha sources in ADMCI and PDMCI. This was true at both group and individual levels.

SIGNIFICANCE

The similar abnormality of alpha source connectivity in ADMCI and PDMCI subjects might reflect common cholinergic impairment.