Bayesian model averaging in EEG/MEG imaging

Early attentional processing and cortical remapping strategies of tactile stimuli in adults with an early and late-onset visual impairment: A cross-sectional study

Neuroplastic changes appear in people with visual impairment (VI) and they show greater tactile abilities. Improvements in performance could be associated with the development of enhanced early attentional processes based on neuroplasticity. Currently, the various early attentional and cortical remapping strategies that are utilized by people with early (EB) and late-onset blindness (LB) remain unclear. Thus, more research is required to develop effective rehabilitation programs and substitution devices. Our objective was to explore the differences in spatial tactile brain processing in adults with EB, LB and a sighted control group (CG). In this cross-sectional study 27 participants with VI were categorized into EB (n = 14) and LB (n = 13) groups. They were then compared with a CG (n = 15). A vibrotactile device and event-related potentials (ERPs) were utilized while participants performed a spatial tactile line recognition task. The P100 latency and cortical areas of maximal activity were analyzed during the task. The three groups had no statistical differences in P100 latency (p>0.05). All subjects showed significant activation in the right superior frontal areas. Only individuals with VI activated the left superior frontal regions. In EB subjects, a higher activation was found in the mid-frontal and occipital areas. A higher activation of the mid-frontal, anterior cingulate cortex and orbitofrontal zones was observed in LB participants. Compared to the CG, LB individuals showed greater activity in the left orbitofrontal zone, while EB exhibited greater activity in the right superior parietal cortex. The EB had greater activity in the left orbitofrontal region compared to the LB. People with VI may not have faster early attentional processing. EB subjects activate the occipital lobe and right superior parietal cortex during tactile stimulation because of an early lack of visual stimuli and a multimodal information processing. In individuals with LB and EB the orbitofrontal area is activated, suggesting greater emotional processing.

Computational Synaptic Modeling of Pitch and Duration Mismatch Negativity in First-Episode Psychosis Reveals Selective Dysfunction of the N-Methyl-D-Aspartate Receptor

Mismatch negativity (MMN) to pitch (pMMN) and to duration (dMMN) deviant stimuli is significantly more attenuated in long-term psychotic illness compared to first-episode psychosis (FEP). It was recently shown that source-modeling of magnetically recorded MMN increases the detection of left auditory cortex MMN deficits in FEP, and that computational circuit modeling of electrically recorded MMN also reveals left-hemisphere auditory cortex abnormalities. Computational modeling using dynamic causal modeling (DCM) can also be used to infer synaptic activity from EEG-based scalp recordings. We measured pMMN and dMMN with EEG from 26 FEP and 26 matched healthy controls (HCs) and used a DCM conductance-based neural mass model including α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid, N-methyl-D-Aspartate (NMDA), and Gamma-aminobutyric acid receptors to identify any changes in effective connectivity and receptor rate constants in FEP. We modeled MMN sources in bilateral A1, superior temporal gyrus, and inferior frontal gyrus (IFG). No model parameters distinguished groups for pMMN. For dMMN, reduced NMDA receptor activity in right IFG in FEP was detected. This finding is in line with literature of prefrontal NMDA receptor hypofunction in chronic schizophrenia and suggests impaired NMDA-induced synaptic plasticity may be present at psychosis onset where scalp dMMN is only moderately reduced. To the best of our knowledge, this is the first report of impaired NMDA receptor activity in FEP found through computational modeling of dMMN and shows the potential of DCM to non-invasively reveal synaptic-level abnormalities that underly subtle functional auditory processing deficits in early psychosis.

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.

Identifying oscillatory brain networks with hidden Gaussian graphical spectral models of MEEG

Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflection of functional brain network activity with a temporal resolution that no other in vivo neuroimage may provide. A second step estimating functional connectivity from such activity pseudodata unveil the oscillatory brain networks that strongly correlate with all cognition and behavior. Simulations of such MEG or EEG inverse problem also reveal estimation errors of the functional connectivity determined by any of the state-of-the-art inverse solutions. We disclose a significant cause of estimation errors originating from misspecification of the functional network model incorporated into either inverse solution steps. We introduce the Bayesian identification of a Hidden Gaussian Graphical Spectral (HIGGS) model specifying such oscillatory brain networks model. In human EEG alpha rhythm simulations, the estimation errors measured as ROC performance do not surpass 2% in our HIGGS inverse solution and reach 20% in state-of-the-art methods. Macaque simultaneous EEG/ECoG recordings provide experimental confirmation for our results with 1/3 times larger congruence according to Riemannian distances than state-of-the-art methods.

Harmonized multi-metric and multi-centric assessment of EEG source space connectivity for dementia characterization

Introduction

Harmonization protocols that address batch effects and cross-site methodological differences in multi-center studies are critical for strengthening electroencephalography (EEG) signatures of functional connectivity (FC) as potential dementia biomarkers.

Methods

We implemented an automatic processing pipeline incorporating electrode layout integrations, patient-control normalizations, and multi-metric EEG source space connectomics analyses.

Results

Spline interpolations of EEG signals onto a head mesh model with 6067 virtual electrodes resulted in an effective method for integrating electrode layouts. Z-score transformations of EEG time series resulted in source space connectivity matrices with high bilateral symmetry, reinforced long-range connections, and diminished short-range functional interactions. A composite FC metric allowed for accurate multicentric classifications of Alzheimer's disease and behavioral variant frontotemporal dementia.

Discussion

Harmonized multi-metric analysis of EEG source space connectivity can address data heterogeneities in multi-centric studies, representing a powerful tool for accurately characterizing dementia.

Automatic Quality Control of the numerical accuracy of EEG Lead fields

Precise individualized EEG source localization is predicated on having accurate subject-specific Lead Fields (LFs) obtained from their Magnetic Resonance Images (MRI). LF calculation is a complex process involving several error-prone steps that start with obtaining a realistic head model from the MRI and finalizing with computationally expensive solvers such as the Boundary Element Method (BEM) or Finite Element Method (FEM). Current Big-Data applications require the calculation of batches of hundreds or thousands of LFs. LF. Quality Control is conventionally checked subjectively by experts, a procedure not feasible in practice for larger batches. To facilitate this step, we introduce the Lead Field Automatic-Quality Control Index (LF-AQI) that flags LF with potential errors. We base our LF-AQI on the assumption that LFs obtained from simpler head models, i.e., the homogeneous head model LF (HHM-LF) or spherical head model LF (SHM-LF), deviate only moderately from a "good" realistic test LF. Since these simpler LFs are easier to compute and check for errors, they may serve as "reference LF" to detect anomalous realistic test LF. We investigated this assumption by comparing correlation-based channel ρ_min(ref,test)and source τ_min(ref,test) similarity indices (SI) between "gold standards," i.e., very accurate FEM and BEM LFs, and the proposed references (HHM-LF and SHM-LF). Surprisingly we found that the most uncomplicated possible reference, HHM-LF had high SI values with the gold standards-leading us to explore further use of the channel ρ_min(HHM-LF,test)and source τ_min(HHM-LF,test)_{SI as a basis for} our LF-AQI. Indeed, these SI successfully detected five simulated scenarios of LFs artifacts. This result encouraged us to evaluate the SI on a large dataset and thus define our LF-AQI. We thus computed the SI of 1251 LFs obtained from the Child Mind Institute (CMI) MRI dataset. When ρ_min(HHM-LF,test)and source τ_min(HHM-LF,test) were plotted for all test subjects on a 2D space, most were tightly clustered around the median of a high similarity centroid (HSC), except for a smaller proportion of outliers. We define the LF-AQI for a given LF as the log Euclidean distance between its SI and the HSC median. To automatically detect outliers, the threshold is at the 90^th percentile of the CMI LF-AQIs (-0.9755). LF-AQI greater than this threshold flag individual LF to be checked. The robustness of this LF-AQI screening was checked by repeated out-of-sample validation. Strikingly, minor corrections in re-processing the flagged cases eliminated their status as outliers. Furthermore, the "doubtful" labels assigned by LF-AQI were validated by neuroscience students using a Likert scale questionnaire designed to manually check the LF's quality. Item Response Theory (IRT) analysis was applied to the questionnaire results to compute an optimized model and a latent variable θ for that model. A linear mixed model (LMM) between the θ and LF-AQI resulted in an effect with a Cohen's d value of 1.3 and a p-value <0.001, thus validating the correspondence of LF-AQI with the visual quality control. We provide an open-source pipeline to implement both LF calculation and its quality control to allow further evaluation of our index.

Comparing encoding mechanisms in realistic virtual reality and conventional 2D laboratory settings: Event-related potentials in a repetition suppression paradigm

Although the human brain is adapted to function within three-dimensional environments, conventional laboratory research commonly investigates cognitive mechanisms in a reductionist approach using two-dimensional stimuli. However, findings regarding mnemonic processes indicate that realistic experiences in Virtual Reality (VR) are stored in richer and more intertwined engrams than those obtained from the conventional laboratory. Our study aimed to further investigate the generalizability of laboratory findings and to differentiate whether the processes underlying memory formation differ between VR and the conventional laboratory already in early encoding stages. Therefore, we investigated the Repetition Suppression (RS) effect as a correlate of the earliest instance of mnemonic processes under conventional laboratory conditions and in a realistic virtual environment. Analyses of event-related potentials (ERPs) indicate that the ERP deflections at several electrode clusters were lower in VR compared to the PC condition. These results indicate an optimized distribution of cognitive resources in realistic contexts. The typical RS effect was replicated under both conditions at most electrode clusters for a late time window. Additionally, a specific RS effect was found in VR at anterior electrodes for a later time window, indicating more extensive encoding processes in VR compared to the laboratory. Specifically, electrotomographic results (VARETA) indicate multimodal integration involving a broad cortical network and higher cognitive processes during the encoding of realistic objects. Our data suggest that object perception under realistic conditions, in contrast to the conventional laboratory, requires multisensory integration involving an interconnected functional system, facilitating the formation of intertwined memory traces in realistic environments.

The IAS-MEEG Package: A Flexible Inverse Source Reconstruction Platform for Reconstruction and Visualization of Brain Activity from M/EEG Data

We present a standalone Matlab software platform complete with visualization for the reconstruction of the neural activity in the brain from MEG or EEG data. The underlying inversion combines hierarchical Bayesian models and Krylov subspace iterative least squares solvers. The Bayesian framework of the underlying inversion algorithm allows to account for anatomical information and possible a priori belief about the focality of the reconstruction. The computational efficiency makes the software suitable for the reconstruction of lengthy time series on standard computing equipment. The algorithm requires minimal user provided input parameters, although the user can express the desired focality and accuracy of the solution. The code has been designed so as to favor the parallelization performed automatically by Matlab, according to the resources of the host computer. We demonstrate the flexibility of the platform by reconstructing activity patterns with supports of different sizes from MEG and EEG data. Moreover, we show that the software reconstructs well activity patches located either in the subcortical brain structures or on the cortex. The inverse solver and visualization modules can be used either individually or in combination. We also provide a version of the inverse solver that can be used within Brainstorm toolbox. All the software is available online by Github, including the Brainstorm plugin, with accompanying documentation and test data.

Enhanced Mirror Neuron Network Activity and Effective Connectivity during Live Interaction Among Female Subjects

Facial expressions are indispensable in daily human communication. Previous neuroimaging studies investigating facial expression processing have presented pre-recorded stimuli and lacked live face-to-face interaction. Our paradigm alternated between presentations of real-time model performance and pre-recorded videos of dynamic facial expressions to participants. Simultaneous functional magnetic resonance imaging (fMRI) and facial electromyography activity recordings, as well as post-scan valence and arousal ratings were acquired from 44 female participants. Live facial expressions enhanced the subjective valence and arousal ratings as well as facial muscular responses. Live performances showed greater engagement of the right posterior superior temporal sulcus (pSTS), right inferior frontal gyrus (IFG), right amygdala and right fusiform gyrus, and modulated the effective connectivity within the right mirror neuron system (IFG, pSTS, and right inferior parietal lobule). A support vector machine algorithm could classify multivoxel activation patterns in brain regions involved in dynamic facial expression processing in the mentalizing networks (anterior and posterior cingulate cortex). These results indicate that live social interaction modulates the activity and connectivity of the right mirror neuron system and enhances spontaneous mimicry, further facilitating emotional contagion.

Evaluating Effective Connectivity of Trust in Human-Automation Interaction: A Dynamic Causal Modeling (DCM) Study

OBJECTIVE

Using dynamic causal modeling (DCM), we examined how credibility and reliability affected the way brain regions exert causal influence over each other-effective connectivity (EC)-in the context of trust in automation.

BACKGROUND

Multiple brain regions of the central executive network (CEN) and default mode network (DMN) have been implicated in trust judgment. However, the neural correlates of trust judgment are still relatively unexplored in terms of the directed information flow between brain regions.

METHOD

Sixteen participants observed the performance of four computer algorithms, which differed in credibility and reliability, of the system monitoring subtask of the Air Force Multi-Attribute Task Battery (AF-MATB). Using six brain regions of the CEN and DMN commonly identified to be activated in human trust, a total of 30 (forward, backward, and lateral) connection models were developed. Bayesian model averaging (BMA) was used to quantify the connectivity strength among the brain regions.

RESULTS

Relative to the high trust condition, low trust showed unique presence of specific connections, greater connectivity strengths from the prefrontal cortex, and greater network complexity. High trust condition showed no backward connections.

CONCLUSION

Results indicated that trust and distrust can be two distinctive neural processes in human-automation interaction-distrust being a more complex network than trust, possibly due to the increased cognitive load.

APPLICATION

The causal architecture of distributed brain regions inferred using DCM can help not only in the design of a balanced human-automation interface design but also in the proper use of automation in real-life situations.

Allostatic-Interoceptive Overload in Frontotemporal Dementia

BACKGROUND

The predictive coding theory of allostatic-interoceptive load states that brain networks mediating autonomic regulation and interoceptive-exteroceptive balance regulate the internal milieu to anticipate future needs and environmental demands. These functions seem to be distinctly compromised in behavioral variant frontotemporal dementia (bvFTD), including alterations of the allostatic-interoceptive network (AIN). Here, we hypothesize that bvFTD is typified by an allostatic-interoceptive overload.

METHODS

We assessed resting-state heartbeat evoked potential (rsHEP) modulation as well as its behavioral and multimodal neuroimaging correlates in patients with bvFTD relative to healthy control subjects and patients with Alzheimer's disease (N = 94). We measured 1) resting-state electroencephalography (to assess the rsHEP, prompted by visceral inputs and modulated by internal body sensing), 2) associations between rsHEP and its neural generators (source location), 3) cognitive disturbances (cognitive state, executive functions, facial emotion recognition), 4) brain atrophy, and 5) resting-state functional magnetic resonance imaging functional connectivity (AIN vs. control networks).

RESULTS

Relative to healthy control subjects and patients with Alzheimer's disease, patients with bvFTD presented more negative rsHEP amplitudes with sources in critical hubs of the AIN (insula, amygdala, somatosensory cortex, hippocampus, anterior cingulate cortex). This exacerbated rsHEP modulation selectively predicted the patients' cognitive profile (including cognitive decline, executive dysfunction, and emotional impairments). In addition, increased rsHEP modulation in bvFTD was associated with decreased brain volume and connectivity of the AIN. Machine learning results confirmed AIN specificity in predicting the bvFTD group.

CONCLUSIONS

Altogether, these results suggest that bvFTD may be characterized by an allostatic-interoceptive overload manifested in ongoing electrophysiological markers, brain atrophy, functional networks, and cognition.

Early protein energy malnutrition impacts life-long developmental trajectories of the sources of EEG rhythmic activity

Protein Energy Malnutrition (PEM) has lifelong consequences on brain development and cognitive function. We studied the lifelong developmental trajectories of resting-state EEG source activity in 66 individuals with histories of Protein Energy Malnutrition (PEM) limited to the first year of life and in 83 matched classmate controls (CON) who are all participants of the 49 years longitudinal Barbados Nutrition Study (BNS). qEEGt source z-spectra measured deviation from normative values of EEG rhythmic activity sources at 5-11 years of age and 40 years later at 45-51 years of age. The PEM group showed qEEGt abnormalities in childhood, including a developmental delay in alpha rhythm maturation and an insufficient decrease in beta activity. These profiles may be correlated with accelerated cognitive decline.

EEG effective connectivity during the first year of life mirrors brain synaptogenesis, myelination, and early right hemisphere predominance

INTRODUCTION

The maturation of electroencephalogram (EEG) effective connectivity in healthy infants during the first year of life is described.

METHODS

Participants: A cross-sectional sample of 125 healthy at-term infants, from 0 to 12 months of age, underwent EEG in a state of quiet sleep.

PROCEDURES

The EEG primary currents at the source were described with the sLoreta method. An unmixing algorithm was applied to reduce the leakage, and the isolated effective coherence, a direct and directed measurement of information flow, was calculated.

RESULTS AND DISCUSSION

Initially, the highest indices of connectivity are at the subcortical nuclei, continuing to the parietal lobe, predominantly the right hemisphere, then expanding to temporal, occipital, and finally the frontal areas, which is consistent with the myelination process. Age-related connectivity changes were mostly long-range and bilateral. Connections increased with age, mainly in the right hemisphere, while they mainly decreased in the left hemisphere. Increased connectivity from 20 to 30 Hz, mostly at the right hemisphere. These findings were consistent with right hemisphere predominance during the first three years of life. Theta and alpha connections showed the greatest changes with age. Strong connectivity was found between the parietal, temporal, and occipital regions to the frontal lobes, responsible for executive functions and consistent with behavioral development during the first year. The thalamus exchanges information bidirectionally with all cortical regions and frequency bands.

CONCLUSIONS

The maturation of EEG connectivity during the first year in healthy infants is very consistent with synaptogenesis, reductions in synaptogenesis, myelination, and functional and behavioral development.

Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean

Functional near-infrared spectroscopy (fNIRS) measures the hemoglobin concentration changes associated with neuronal activity. Diffuse optical tomography (DOT) consists of reconstructing the optical density changes measured from scalp channels to the oxy-/deoxy-hemoglobin concentration changes within the cortical regions. In the present study, we adapted a nonlinear source localization method developed and validated in the context of Electro- and Magneto-Encephalography (EEG/MEG): the Maximum Entropy on the Mean (MEM), to solve the inverse problem of DOT reconstruction. We first introduced depth weighting strategy within the MEM framework for DOT reconstruction to avoid biasing the reconstruction results of DOT towards superficial regions. We also proposed a new initialization of the MEM model improving the temporal accuracy of the original MEM framework. To evaluate MEM performance and compare with widely used depth weighted Minimum Norm Estimate (MNE) inverse solution, we applied a realistic simulation scheme which contained 4000 simulations generated by 250 different seeds at different locations and 4 spatial extents ranging from 3 to 40[Formula: see text] along the cortical surface. Our results showed that overall MEM provided more accurate DOT reconstructions than MNE. Moreover, we found that MEM was remained particularly robust in low signal-to-noise ratio (SNR) conditions. The proposed method was further illustrated by comparing to functional Magnetic Resonance Imaging (fMRI) activation maps, on real data involving finger tapping tasks with two different montages. The results showed that MEM provided more accurate HbO and HbR reconstructions in spatial agreement with the main fMRI cluster, when compared to MNE.

Optimizing EEG Source Reconstruction with Concurrent fMRI-Derived Spatial Priors

Reconstructing EEG sources involves a complex pipeline, with the inverse problem being the most challenging. Multiple inversion algorithms are being continuously developed, aiming to tackle the non-uniqueness of this problem, which has been shown to be partially circumvented by including prior information in the inverse models. Despite a few efforts, there are still current and persistent controversies regarding the inversion algorithm of choice and the optimal set of spatial priors to be included in the inversion models. The use of simultaneous EEG-fMRI data is one approach to tackle this problem. The spatial resolution of fMRI makes fMRI derived spatial priors very convenient for EEG reconstruction, however, only task activation maps and resting-state networks (RSNs) have been explored so far, overlooking the recent, but already accepted, notion that brain networks exhibit dynamic functional connectivity fluctuations. The lack of a systematic comparison between different source reconstruction algorithms, considering potentially more brain-informative priors such as fMRI, motivates the search for better reconstruction models. Using simultaneous EEG-fMRI data, here we compared four different inversion algorithms (minimum norm, MN; low resolution electromagnetic tomography, LORETA; empirical Bayes beamformer, EBB; and multiple sparse priors, MSP) under a Bayesian framework (as implemented in SPM), each with three different sets of priors consisting of: (1) those specific to the algorithm; (2) those specific to the algorithm plus fMRI task activation maps and RSNs; and (3) those specific to the algorithm plus fMRI task activation maps and RSNs and network modules of task-related dFC states estimated from the dFC fluctuations. The quality of the reconstructed EEG sources was quantified in terms of model-based metrics, namely the expectation of the posterior probability P(model|data) and variance explained of the inversion models, and the overlap/proportion of brain regions known to be involved in the visual perception tasks that the participants were submitted to, and RSN templates, with/within EEG source components. Model-based metrics suggested that model parsimony is preferred, with the combination MSP and priors specific to this algorithm exhibiting the best performance. However, optimal overlap/proportion values were found using EBB and priors specific to this algorithm and fMRI task activation maps and RSNs or MSP and considering all the priors (algorithm priors, fMRI task activation maps and RSNs and dFC state modules), respectively, indicating that fMRI spatial priors, including dFC state modules, might contain useful information to recover EEG source components reflecting neuronal activity of interest. Our main results show that providing fMRI spatial derived priors that reflect the dynamics of the brain might be useful to map neuronal activity more accurately from EEG-fMRI. Furthermore, this work paves the way towards a more informative selection of the optimal EEG source reconstruction approach, which may be critical in future studies.

Dementia ConnEEGtome: Towards multicentric harmonization of EEG connectivity in neurodegeneration

The proposal to use brain connectivity as a biomarker for dementia phenotyping can be potentiated by conducting large-scale multicentric studies using high-density electroencephalography (hd- EEG). Nevertheless, several barriers preclude the development of a systematic "ConnEEGtome" in dementia research. Here we review critical sources of variability in EEG connectivity studies, and provide general guidelines for multicentric protocol harmonization. We describe how results can be impacted by the choice for data acquisition, and signal processing workflows. The implementation of a particular processing pipeline is conditional upon assumptions made by researchers about the nature of EEG. Due to these assumptions, EEG connectivity metrics are typically applicable to restricted scenarios, e.g., to a particular neurocognitive disorder. "Ground truths" for the choice of processing workflow and connectivity analysis are impractical. Consequently, efforts should be directed to harmonizing experimental procedures, data acquisition, and the first steps of the preprocessing pipeline. Conducting multiple analyses of the same data and a proper integration of the results need to be considered in additional processing steps. Furthermore, instead of using a single connectivity measure, using a composite metric combining different connectivity measures brings a powerful strategy to scale up the replicability of multicentric EEG connectivity studies. These composite metrics can boost the predictive strength of diagnostic tools for dementia. Moreover, the implementation of multi-feature machine learning classification systems that include EEG-based connectivity analyses may help to exploit the potential of multicentric studies combining clinical-cognitive, molecular, genetics, and neuroimaging data towards a multi-dimensional characterization of the dementia.

Evaluation of a personalized functional near infra-red optical tomography workflow using maximum entropy on the mean

In the present study, we proposed and evaluated a workflow of personalized near infra-red optical tomography (NIROT) using functional near-infrared spectroscopy (fNIRS) for spatiotemporal imaging of cortical hemodynamic fluctuations. The proposed workflow from fNIRS data acquisition to local 3D reconstruction consists of: (a) the personalized optimal montage maximizing fNIRS channel sensitivity to a predefined targeted brain region; (b) the optimized fNIRS data acquisition involving installation of optodes and digitalization of their positions using a neuronavigation system; and (c) the 3D local reconstruction using maximum entropy on the mean (MEM) to accurately estimate the location and spatial extent of fNIRS hemodynamic fluctuations along the cortical surface. The workflow was evaluated on finger-tapping fNIRS data acquired from 10 healthy subjects for whom we estimated the reconstructed NIROT spatiotemporal images and compared with functional magnetic resonance imaging (fMRI) results from the same individuals. Using the fMRI activation maps as our reference, we quantitatively compared the performance of two NIROT approaches, the MEM framework and the conventional minimum norm estimation (MNE) method. Quantitative comparisons were performed at both single subject and group-level. Overall, our results suggested that MEM provided better spatial accuracy than MNE, while both methods offered similar temporal accuracy when reconstructing oxygenated (HbO) and deoxygenated hemoglobin (HbR) concentration changes evoked by finger-tapping. Our proposed complete workflow was made available in the brainstorm fNIRS processing plugin-NIRSTORM, thus providing the opportunity for other researchers to further apply it to other tasks and on larger populations.

Unification of sparse Bayesian learning algorithms for electromagnetic brain imaging with the majorization minimization framework

Methods for electro- or magnetoencephalography (EEG/MEG) based brain source imaging (BSI) using sparse Bayesian learning (SBL) have been demonstrated to achieve excellent performance in situations with low numbers of distinct active sources, such as event-related designs. This paper extends the theory and practice of SBL in three important ways. First, we reformulate three existing SBL algorithms under the majorization-minimization (MM) framework. This unification perspective not only provides a useful theoretical framework for comparing different algorithms in terms of their convergence behavior, but also provides a principled recipe for constructing novel algorithms with specific properties by designing appropriate bounds of the Bayesian marginal likelihood function. Second, building on the MM principle, we propose a novel method called LowSNR-BSI that achieves favorable source reconstruction performance in low signal-to-noise-ratio (SNR) settings. Third, precise knowledge of the noise level is a crucial requirement for accurate source reconstruction. Here we present a novel principled technique to accurately learn the noise variance from the data either jointly within the source reconstruction procedure or using one of two proposed cross-validation strategies. Empirically, we could show that the monotonous convergence behavior predicted from MM theory is confirmed in numerical experiments. Using simulations, we further demonstrate the advantage of LowSNR-BSI over conventional SBL in low-SNR regimes, and the advantage of learned noise levels over estimates derived from baseline data. To demonstrate the usefulness of our novel approach, we show neurophysiologically plausible source reconstructions on averaged auditory evoked potential data.

Estimating multiple latencies in the auditory system from auditory steady-state responses on a single EEG channel

The latency of the auditory steady-state response (ASSR) may provide valuable information regarding the integrity of the auditory system, as it could potentially reveal the presence of multiple intracerebral sources. To estimate multiple latencies from high-order ASSRs, we propose a novel two-stage procedure that consists of a nonparametric estimation method, called apparent latency from phase coherence (ALPC), followed by a heuristic sequential forward selection algorithm (SFS). Compared with existing methods, ALPC-SFS requires few prior assumptions, and is straightforward to implement for higher-order nonlinear responses to multi-cosine sound complexes with their initial phases set to zero. It systematically evaluates the nonlinear components of the ASSRs by estimating multiple latencies, automatically identifies involved ASSR components, and reports a latency consistency index. To verify the proposed method, we performed simulations for several scenarios: two nonlinear subsystems with different or overlapping outputs. We compared the results from our method with predictions from existing, parametric methods. We also recorded the EEG from ten normal-hearing adults by bilaterally presenting superimposed tones with four frequencies that evoke a unique set of ASSRs. From these ASSRs, two major latencies were found to be stable across subjects on repeated measurement days. The two latencies are dominated by low-frequency (LF) (near 40 Hz, at around 41-52 ms) and high-frequency (HF) (> 80 Hz, at around 21-27 ms) ASSR components. The frontal-central brain region showed longer latencies on LF components, but shorter latencies on HF components, when compared with temporal-lobe regions. In conclusion, the proposed nonparametric ALPC-SFS method, applied to zero-phase, multi-cosine sound complexes is more suitable for evaluating embedded nonlinear systems underlying ASSRs than existing methods. It may therefore be a promising objective measure for hearing performance and auditory cortex (dys)function.

Different and common brain signals of altered neurocognitive mechanisms for unfamiliar face processing in acquired and developmental prosopagnosia

Neuropsychological studies have shown that prosopagnosic individuals perceive face structure in an atypical way. This might preclude the formation of appropriate face representations and, consequently, hamper effective recognition. The present ERP study, in combination with Bayesian source reconstruction, investigates how information related to both external (E) and internal (I) features was processed by E.C. and I.P., suffering from acquired and developmental prosopagnosia, respectively. They carried out a face-feature matching task with new faces. E.C. showed poor performance and remarkable lack of early face-sensitive P1, N170 and P2 responses on right (damaged) posterior cortex. Although she presented the expected mismatch effect to target faces in the E-I sequence, it was of shorter duration than in Controls, and involved left parietal, right frontocentral and dorsofrontal regions, suggestive of reduced neural circuitry to process face configurations. In turn, I.P. performed efficiently but with a remarkable bias to give "match" responses. His face-sensitive potentials P1-N170 were comparable to those from Controls, however, he showed no subsequent P2 response and a mismatch effect only in the I-E sequence, reflecting activation confined to those regions that sustain typically the initial stages of face processing. Relevantly, neither of the prosopagnosics exhibited conspicuous P3 responses to features acting as primes, indicating that diagnostic information for constructing face representations could not be sufficiently attended nor deeply encoded. Our findings suggest a different locus for altered neurocognitive mechanisms in the face network in participants with different types of prosopagnosia, but common indicators of a deficient allocation of attentional resources for further recognition.

Accuracy and spatial properties of distributed magnetic source imaging techniques in the investigation of focal epilepsy patients

Source localization of interictal epileptiform discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. We aimed to compare the performance of four different distributed magnetic source imaging (dMSI) approaches: Minimum norm estimate (MNE), dynamic statistical parametric mapping (dSPM), standardized low-resolution electromagnetic tomography (sLORETA), and coherent maximum entropy on the mean (cMEM). We also evaluated whether a simple average of maps obtained from multiple inverse solutions (Ave) can improve localization accuracy. We analyzed dMSI of 206 IEDs derived from magnetoencephalography recordings in 28 focal epilepsy patients who had a well-defined focus determined through intracranial EEG (iEEG), epileptogenic MRI lesions or surgical resection. dMSI accuracy and spatial properties were quantitatively estimated as: (a) distance from the epilepsy focus, (b) reproducibility, (c) spatial dispersion (SD), (d) map extension, and (e) effect of thresholding on map properties. Clinical performance was excellent for all methods (median distance from the focus MNE = 2.4 mm; sLORETA = 3.5 mm; cMEM = 3.5 mm; dSPM = 6.8 mm, Ave = 0 mm). Ave showed the lowest distance between the map maximum and epilepsy focus (Dmin lower than cMEM, MNE, and dSPM, p = .021, p = .008, p < .001, respectively). cMEM showed the best spatial features, with lowest SD outside the focus (SD lower than all other methods, p < .001 consistently) and high contrast between the generator and surrounding regions. The average map Ave provided the best localization accuracy, whereas cMEM exhibited the lowest amount of spurious distant activity. dMSI techniques have the potential to significantly improve identification of iEEG targets and to guide surgical planning, especially when multiple methods are combined.

A N400 ERP Study in letter recognition after passive tactile stimulation training in blind children and sighted controls

BACKGROUND

We previously demonstrated that using a sensory substitution device (SSD) for one week, tactile stimulation results in faster activation of lateral occipital complex in blind children than in seeing controls.

OBJECTIVE

We used long-term haptic tactile stimulation training with an SSD to test if it results in stable cross-modal reassignment of visual pathways after six months, to provide high level processing of tactile semantic content.

METHODS

We enrolled 12 blind and 12 sighted children. The SSD transforms images to a stimulation matrix in contact with the dominant hand. Subjects underwent twice-daily training sessions, 5 days/week for six months. Children were asked to describe line orientation, name letters, and read words. ERP sessions were performed at baseline and 6 months to analyze the N400 ERP component and reaction times (RT). N400 sources were estimated with Low Resolution Electromagnetic Tomography (LORETA). SPM8 was used to make population-level inferences.

RESULTS

We found no group differences in RTs, accuracy of identifications, N400 latencies or distributions with the line task at 1 week or at 6 months. RTs on the letter recognition task were also similar. After 6 months, behavioral training increased accurate letter identification in both seeing and blind children (Chi 2 = 11906.934, p = 0.000), but the increase was larger in blind children (Chi 2 = 8.272, p = 0.004). Behavioral training shifted peak N400 amplitude to left occipital and bilateral parietal cortices in blind children, but to left precentral and postcentral and bilateral occipital cortices in sighted controls.

CONCLUSIONS

Blind children learn to recognize SSD-delivered letters better than seeing controls and had greater N400 amplitude in the occipital region. To the best of our knowledge, our results provide the first published example of standard letter recognition (not Braille) by children with blindness using a tactile delivery system.

Robust Empirical Bayesian Reconstruction of Distributed Sources for Electromagnetic Brain Imaging

Electromagnetic brain imaging is the reconstruction of brain activity from non-invasive recordings of the magnetic fields and electric potentials. An enduring challenge in this imaging modality is estimating the number, location, and time course of sources, especially for the reconstruction of distributed brain sources with complex spatial extent. Here, we introduce a novel robust empirical Bayesian algorithm that enables better reconstruction of distributed brain source activity with two key ideas: kernel smoothing and hyperparameter tiling. Since the proposed algorithm builds upon many of the performance features of the sparse source reconstruction algorithm - Champagne and we refer to this algorithm as Smooth Champagne. Smooth Champagne is robust to the effects of high levels of noise, interference, and highly correlated brain source activity. Simulations demonstrate excellent performance of Smooth Champagne when compared to benchmark algorithms in accurately determining the spatial extent of distributed source activity. Smooth Champagne also accurately reconstructs real MEG and EEG data.

Improved EEG source localization with Bayesian uncertainty modelling of unknown skull conductivity

Electroencephalography (EEG) source imaging is an ill-posed inverse problem that requires accurate conductivity modelling of the head tissues, especially the skull. Unfortunately, the conductivity values are difficult to determine in vivo. In this paper, we show that the exact knowledge of the skull conductivity is not always necessary when the Bayesian approximation error (BAE) approach is exploited. In BAE, we first postulate a probability distribution for the skull conductivity that describes our (lack of) knowledge on its value, and model the effects of this uncertainty on EEG recordings with the help of an additive error term in the observation model. Before the Bayesian inference, the likelihood is marginalized over this error term. Thus, in the inversion we estimate only our primary unknown, the source distribution. We quantified the improvements in the source localization when the proposed Bayesian modelling was used in the presence of different skull conductivity errors and levels of measurement noise. Based on the results, BAE was able to improve the source localization accuracy, particularly when the unknown (true) skull conductivity was much lower than the expected standard conductivity value. The source locations that gained the highest improvements were shallow and originally exhibited the largest localization errors. In our case study, the benefits of BAE became negligible when the signal-to-noise ratio dropped to 20 dB.

Brain Activity Mapping from MEG Data via a Hierarchical Bayesian Algorithm with Automatic Depth Weighting

A recently proposed iterated alternating sequential (IAS) MEG inverse solver algorithm, based on the coupling of a hierarchical Bayesian model with computationally efficient Krylov subspace linear solver, has been shown to perform well for both superficial and deep brain sources. However, a systematic study of its ability to correctly identify active brain regions is still missing. We propose novel statistical protocols to quantify the performance of MEG inverse solvers, focusing in particular on how their accuracy and precision at identifying active brain regions. We use these protocols for a systematic study of the performance of the IAS MEG inverse solver, comparing it with three standard inversion methods, wMNE, dSPM, and sLORETA. To avoid the bias of anecdotal tests towards a particular algorithm, the proposed protocols are Monte Carlo sampling based, generating an ensemble of activity patches in each brain region identified in a given atlas. The performance in correctly identifying the active areas is measured by how much, on average, the reconstructed activity is concentrated in the brain region of the simulated active patch. The analysis is based on Bayes factors, interpreting the estimated current activity as data for testing the hypothesis that the active brain region is correctly identified, versus the hypothesis of any erroneous attribution. The methodology allows the presence of a single or several simultaneous activity regions, without assuming that the number of active regions is known. The testing protocols suggest that the IAS solver performs well with both with cortical and subcortical activity estimation.

Comparison of source localization techniques in diffuse optical tomography for fNIRS application using a realistic head model

Functional near-infrared spectroscopy (fNIRS) is a non-invasive imaging technique that elicits growing interest for research and clinical applications. In the last decade, efforts have been made to develop a mathematical framework in order to image the effective sources of hemoglobin variations in brain tissues. Different approaches can be used to impose additional information or constraints when reconstructing the cerebral images of an ill-posed problem. The goal of this study is to compare the performance and limitations of several source localization techniques in the context of fNIRS tomography using individual anatomical magnetic resonance imaging (MRI) to model light propagation. The forward problem is solved using a Monte Carlo simulation of light propagation in the tissues. The inverse problem has been linearized using the Rytov approximation. Then, Tikhonov regularization applied to least squares, truncated singular value decomposition, back-projection, L1-norm regularization, minimum norm estimates, low resolution electromagnetic tomography and Bayesian model averaging techniques are compared using a receiver operating characteristic analysis, blurring and localization error measures. Using realistic simulations (n = 450) and data acquired from a human participant, this study depicts how these source localization techniques behave in a human head fNIRS tomography. When compared to other methods, Bayesian model averaging is proposed as a promising method in DOT and shows great potential to improve specificity, accuracy, as well as to reduce blurring and localization error even in presence of noise and deep sources. Classical reconstruction methods, such as regularized least squares, offer better sensitivity but higher blurring; while more novel L1-based method provides sparse solutions with small blurring and high specificity but lower sensitivity. The application of these methods is also demonstrated experimentally using visual fNIRS experiment with adult participant.

Electrophysiological Source Imaging: A Noninvasive Window to Brain Dynamics

Brain activity and connectivity are distributed in the three-dimensional space and evolve in time. It is important to image brain dynamics with high spatial and temporal resolution. Electroencephalography (EEG) and magnetoencephalography (MEG) are noninvasive measurements associated with complex neural activations and interactions that encode brain functions. Electrophysiological source imaging estimates the underlying brain electrical sources from EEG and MEG measurements. It offers increasingly improved spatial resolution and intrinsically high temporal resolution for imaging large-scale brain activity and connectivity on a wide range of timescales. Integration of electrophysiological source imaging and functional magnetic resonance imaging could further enhance spatiotemporal resolution and specificity to an extent that is not attainable with either technique alone. We review methodological developments in electrophysiological source imaging over the past three decades and envision its future advancement into a powerful functional neuroimaging technology for basic and clinical neuroscience applications.

Unified Bayesian Estimator of EEG Reference at Infinity: rREST (Regularized Reference Electrode Standardization Technique)

The choice of reference for the electroencephalogram (EEG) is a long-lasting unsolved issue resulting in inconsistent usages and endless debates. Currently, both the average reference (AR) and the reference electrode standardization technique (REST) are two primary, apparently irreconcilable contenders. We propose a theoretical framework to resolve this reference issue by formulating both (a) estimation of potentials at infinity, and (b) determination of the reference, as a unified Bayesian linear inverse problem, which can be solved by maximum a posterior estimation. We find that AR and REST are very particular cases of this unified framework: AR results from biophysically non-informative prior; while REST utilizes the prior based on the EEG generative model. To allow for simultaneous denoising and reference estimation, we develop the regularized versions of AR and REST, named rAR and rREST, respectively. Both depend on a regularization parameter that is the noise to signal variance ratio. Traditional and new estimators are evaluated with this framework, by both simulations and analysis of real resting EEGs. Toward this end, we leverage the MRI and EEG data from 89 subjects which participated in the Cuban Human Brain Mapping Project. Generated artificial EEGs—with a known ground truth, show that relative error in estimating the EEG potentials at infinity is lowest for rREST. It also reveals that realistic volume conductor models improve the performances of REST and rREST. Importantly, for practical applications, it is shown that an average lead field gives the results comparable to the individual lead field. Finally, it is shown that the selection of the regularization parameter with Generalized Cross-Validation (GCV) is close to the “oracle” choice based on the ground truth. When evaluated with the real 89 resting state EEGs, rREST consistently yields the lowest GCV. This study provides a novel perspective to the EEG reference problem by means of a unified inverse solution framework. It may allow additional principled theoretical formulations and numerical evaluation of performance.

Fast and robust Block-Sparse Bayesian learning for EEG source imaging

We propose a new Sparse Bayesian Learning (SBL) algorithm that can deliver fast, block-sparse, and robust solutions to the EEG source imaging (ESI) problem in the presence of noisy measurements. Current implementations of the SBL framework are computationally expensive and typically handle fluctuations in the measurement noise using different heuristics that are unsuitable for real-time imaging applications. We address these shortcomings by decoupling the estimation of the sensor noise covariance and the sparsity profile of the sources, thereby yielding an efficient two-stage algorithm. In the first stage, we optimize a simplified non-sparse generative model to get an estimate of the sensor noise covariance and a good initialization of the group-sparsity profile of the sources. Sources obtained at this stage are equivalent to those estimated with the popular inverse method LORETA. In the second stage, we apply a fast SBL algorithm with the noise covariance fixed to the value obtained in the first stage to efficiently shrink to zero groups of sources that are irrelevant for explaining the EEG measurements. In addition, we derive an initialization to the first stage of the algorithm that is optimal in the least squares sense, which prevents delays due to suboptimal initial conditions. We validate our method on both simulated and real EEG data. Simulations show that the method is robust to measurement noise and performs well in real-time, with faster performance than two state of the art SBL solvers. On real error-related negativity EEG data, we obtain source images in agreement with the experimental literature. The method shows promise for real-time neuroimaging and brain-machine interface applications.

Source Reconstruction of Brain Potentials Using Bayesian Model Averaging to Analyze Face Intra-Domain vs. Face-Occupation Cross-Domain Processing

We investigated the neural correlates of the access to and retrieval of face structure information in contrast to those concerning the access to and retrieval of person-related verbal information, triggered by faces. We experimentally induced stimulus familiarity via a systematic learning procedure including faces with and without associated verbal information. Then, we recorded event-related potentials (ERPs) in both intra-domain (face-feature) and cross-domain (face-occupation) matching tasks while N400-like responses were elicited by incorrect eyes-eyebrows completions and occupations, respectively. A novel Bayesian source reconstruction approach plus conjunction analysis of group effects revealed that in both cases the generated N170s were of similar amplitude but had different neural origin. Thus, whereas the N170 of faces was associated predominantly to right fusiform and occipital regions (the so-called “Fusiform Face Area”, “FFA” and “Occipital Face Area”, “OFA”, respectively), the N170 of occupations was associated to a bilateral very posterior activity, suggestive of basic perceptual processes. Importantly, the right-sided perceptual P200 and the face-related N250 were evoked exclusively in the intra-domain task, with sources in OFA and extensively in the fusiform region, respectively. Regarding later latencies, the intra-domain N400 seemed to be generated in right posterior brain regions encompassing mainly OFA and, to some extent, the FFA, likely reflecting neural operations triggered by structural incongruities. In turn, the cross-domain N400 was related to more anterior left-sided fusiform and temporal inferior sources, paralleling those described previously for the classic verbal N400. These results support the existence of differentiated neural streams for face structure and person-related verbal processing triggered by faces, which can be activated differentially according to specific task demands.

Dynamic Causal Modeling of Preclinical Autosomal-Dominant Alzheimer's Disease

Dynamic causal modeling (DCM) is a framework for making inferences about changes in brain connectivity using neuroimaging data. We fitted DCMs to high-density EEG data from subjects performing a semantic picture matching task. The subjects are carriers of the PSEN1 mutation, which leads to early onset Alzheimer's disease, but at the time of EEG acquisition in 1999, these subjects were cognitively unimpaired. We asked 1) what is the optimal model architecture for explaining the event-related potentials in this population, 2) which connections are different between this Presymptomatic Carrier (PreC) group and a Non-Carrier (NonC) group performing the same task, and 3) which network connections are predictive of subsequent Mini-Mental State Exam (MMSE) trajectories. We found 1) a model with hierarchical rather than lateral connections between hemispheres to be optimal, 2) that a pathway from right inferotemporal cortex (IT) to left medial temporal lobe (MTL) was preferentially activated by incongruent items for subjects in the PreC group but not the NonC group, and 3) that increased effective connectivity among left MTL, right IT, and right MTL was predictive of subsequent MMSE scores.

EEG Spectral Generators Involved in Motor Imagery: A swLORETA Study

In order to characterize the neural generators of the brain oscillations related to motor imagery (MI), we investigated the cortical, subcortical, and cerebellar localizations of their respective electroencephalogram (EEG) spectral power and phase locking modulations. The MI task consisted in throwing a ball with the dominant upper limb while in a standing posture, within an ecological virtual reality (VR) environment (tennis court). The MI was triggered by the visual cues common to the control condition, during which the participant remained mentally passive. As previously developed, our paradigm considers the confounding problem that the reference condition allows two complementary analyses: one which uses the baseline before the occurrence of the visual cues in the MI and control resting conditions respectively; and the other which compares the analog periods between the MI and the control resting-state conditions. We demonstrate that MI activates specific, complex brain networks for the power and phase modulations of the EEG oscillations. An early (225 ms) delta phase-locking related to MI was generated in the thalamus and cerebellum and was followed (480 ms) by phase-locking in theta and alpha oscillations, generated in specific cortical areas and the cerebellum. Phase-locking preceded the power modulations (mainly alpha-beta ERD), whose cortical generators were situated in the frontal BA45, BA11, BA10, central BA6, lateral BA13, and posterior cortex BA2. Cerebellar-thalamic involvement through phase-locking is discussed as an underlying mechanism for recruiting at later stages the cortical areas involved in a cognitive role during MI.

Sparsity enables estimation of both subcortical and cortical activity from MEG and EEG

Subcortical structures play a critical role in brain function. However, options for assessing electrophysiological activity in these structures are limited. Electromagnetic fields generated by neuronal activity in subcortical structures can be recorded noninvasively, using magnetoencephalography (MEG) and electroencephalography (EEG). However, these subcortical signals are much weaker than those generated by cortical activity. In addition, we show here that it is difficult to resolve subcortical sources because distributed cortical activity can explain the MEG and EEG patterns generated by deep sources. We then demonstrate that if the cortical activity is spatially sparse, both cortical and subcortical sources can be resolved with M/EEG. Building on this insight, we develop a hierarchical sparse inverse solution for M/EEG. We assess the performance of this algorithm on realistic simulations and auditory evoked response data, and show that thalamic and brainstem sources can be correctly estimated in the presence of cortical activity. Our work provides alternative perspectives and tools for characterizing electrophysiological activity in subcortical structures in the human brain.

Spatio Temporal EEG Source Imaging with the Hierarchical Bayesian Elastic Net and Elitist Lasso Models

The estimation of EEG generating sources constitutes an Inverse Problem (IP) in Neuroscience. This is an ill-posed problem due to the non-uniqueness of the solution and regularization or prior information is needed to undertake Electrophysiology Source Imaging. Structured Sparsity priors can be attained through combinations of (L1 norm-based) and (L2 norm-based) constraints such as the Elastic Net (ENET) and Elitist Lasso (ELASSO) models. The former model is used to find solutions with a small number of smooth nonzero patches, while the latter imposes different degrees of sparsity simultaneously along different dimensions of the spatio-temporal matrix solutions. Both models have been addressed within the penalized regression approach, where the regularization parameters are selected heuristically, leading usually to non-optimal and computationally expensive solutions. The existing Bayesian formulation of ENET allows hyperparameter learning, but using the computationally intensive Monte Carlo/Expectation Maximization methods, which makes impractical its application to the EEG IP. While the ELASSO have not been considered before into the Bayesian context. In this work, we attempt to solve the EEG IP using a Bayesian framework for ENET and ELASSO models. We propose a Structured Sparse Bayesian Learning algorithm based on combining the Empirical Bayes and the iterative coordinate descent procedures to estimate both the parameters and hyperparameters. Using realistic simulations and avoiding the inverse crime we illustrate that our methods are able to recover complicated source setups more accurately and with a more robust estimation of the hyperparameters and behavior under different sparsity scenarios than classical LORETA, ENET and LASSO Fusion solutions. We also solve the EEG IP using data from a visual attention experiment, finding more interpretable neurophysiological patterns with our methods. The Matlab codes used in this work, including Simulations, Methods, Quality Measures and Visualization Routines are freely available in a public website.

Functional connectivity analysis in EEG source space: The choice of method

Functional connectivity (FC) is among the most informative features derived from EEG. However, the most straightforward sensor-space analysis of FC is unreliable owing to volume conductance effects. An alternative-source-space analysis of FC-is optimal for high- and mid-density EEG (hdEEG, mdEEG); however, it is questionable for widely used low-density EEG (ldEEG) because of inadequate surface sampling. Here, using simulations, we investigate the performance of the two source FC methods, the inverse-based source FC (ISFC) and the cortical partial coherence (CPC). To examine the effects of localization errors of the inverse method on the FC estimation, we simulated an oscillatory source with varying locations and SNRs. To compare the FC estimations by the two methods, we simulated two synchronized sources with varying between-source distance and SNR. The simulations were implemented for hdEEG, mdEEG, and ldEEG. We showed that the performance of both methods deteriorates for deep sources owing to their inaccurate localization and smoothing. The accuracy of both methods improves with the increasing between-source distance. The best ISFC performance was achieved using hd/mdEEG, while the best CPC performance was observed with ldEEG. In conclusion, with hdEEG, ISFC outperforms CPC and therefore should be the preferred method. In the studies based on ldEEG, the CPC is a method of choice.

EEG-Based Quantification of Cortical Current Density and Dynamic Causal Connectivity Generalized across Subjects Performing BCI-Monitored Cognitive Tasks

Quantification of dynamic causal interactions among brain regions constitutes an important component of conducting research and developing applications in experimental and translational neuroscience. Furthermore, cortical networks with dynamic causal connectivity in brain-computer interface (BCI) applications offer a more comprehensive view of brain states implicated in behavior than do individual brain regions. However, models of cortical network dynamics are difficult to generalize across subjects because current electroencephalography (EEG) signal analysis techniques are limited in their ability to reliably localize sources across subjects. We propose an algorithmic and computational framework for identifying cortical networks across subjects in which dynamic causal connectivity is modeled among user-selected cortical regions of interest (ROIs). We demonstrate the strength of the proposed framework using a "reach/saccade to spatial target" cognitive task performed by 10 right-handed individuals. Modeling of causal cortical interactions was accomplished through measurement of cortical activity using (EEG), application of independent component clustering to identify cortical ROIs as network nodes, estimation of cortical current density using cortically constrained low resolution electromagnetic brain tomography (cLORETA), multivariate autoregressive (MVAR) modeling of representative cortical activity signals from each ROI, and quantification of the dynamic causal interaction among the identified ROIs using the Short-time direct Directed Transfer function (SdDTF). The resulting cortical network and the computed causal dynamics among its nodes exhibited physiologically plausible behavior, consistent with past results reported in the literature. This physiological plausibility of the results strengthens the framework's applicability in reliably capturing complex brain functionality, which is required by applications, such as diagnostics and BCI.

The impact of MEG source reconstruction method on source-space connectivity estimation: A comparison between minimum-norm solution and beamforming

Despite numerous important contributions, the investigation of brain connectivity with magnetoencephalography (MEG) still faces multiple challenges. One critical aspect of source-level connectivity, largely overlooked in the literature, is the putative effect of the choice of the inverse method on the subsequent cortico-cortical coupling analysis. We set out to investigate the impact of three inverse methods on source coherence detection using simulated MEG data. To this end, thousands of randomly located pairs of sources were created. Several parameters were manipulated, including inter- and intra-source correlation strength, source size and spatial configuration. The simulated pairs of sources were then used to generate sensor-level MEG measurements at varying signal-to-noise ratios (SNR). Next, the source level power and coherence maps were calculated using three methods (a) L2-Minimum-Norm Estimate (MNE), (b) Linearly Constrained Minimum Variance (LCMV) beamforming, and (c) Dynamic Imaging of Coherent Sources (DICS) beamforming. The performances of the methods were evaluated using Receiver Operating Characteristic (ROC) curves. The results indicate that beamformers perform better than MNE for coherence reconstructions if the interacting cortical sources consist of point-like sources. On the other hand, MNE provides better connectivity estimation than beamformers, if the interacting sources are simulated as extended cortical patches, where each patch consists of dipoles with identical time series (high intra-patch coherence). However, the performance of the beamformers for interacting patches improves substantially if each patch of active cortex is simulated with only partly coherent time series (partial intra-patch coherence). These results demonstrate that the choice of the inverse method impacts the results of MEG source-space coherence analysis, and that the optimal choice of the inverse solution depends on the spatial and synchronization profile of the interacting cortical sources. The insights revealed here can guide method selection and help improve data interpretation regarding MEG connectivity estimation.

Oscillatory EEG dynamics underlying automatic chunking during sentence processing

Sentences are easier to remember than random word sequences, likely because linguistic regularities facilitate chunking of words into meaningful groups. The present electroencephalography study investigated the neural oscillations modulated by this so-called sentence superiority effect during the encoding and maintenance of sentence fragments versus word lists. We hypothesized a chunking-related modulation of neural processing during the encoding and retention of sentences (i.e., sentence fragments) as compared to word lists. Time-frequency analysis revealed a two-fold oscillatory pattern for the memorization of sentences: Sentence encoding was accompanied by higher delta amplitude (4Hz), originating both from regions processing syntax as well as semantics (bilateral superior/middle temporal regions and fusiform gyrus). Subsequent sentence retention was reflected in decreased theta (6Hz) and beta/gamma (27-32Hz) amplitude instead. Notably, whether participants simply read or properly memorized the sentences did not impact chunking-related activity during encoding. Therefore, we argue that the sentence superiority effect is grounded in highly automatized language processing mechanisms, which generate meaningful memory chunks irrespective of task demands.

Real-time EEG Source-mapping Toolbox (REST): Online ICA and source localization

The Electroencephalogram (EEG) is a noninvasive functional brain activity recording method that shows promise for becoming a 3-D cortical imaging modality with high temporal resolution. Currently, most of the tools developed for EEG analysis focus mainly on offline processing. This study introduces and demonstrates the Real-time EEG Source-mapping Toolbox (REST), an extension to the widely distributed EEGLAB environment. REST allows blind source separation of EEG data in real-time using Online Recursive Independent Component Analysis (ORICA), plus near real-time localization of separated sources. Two source localization methods are available to fit equivalent current dipoles or estimate spatial source distributions of selected sources. Selected measures of raw EEG data or component activations (e.g. time series of the data, spectral changes over time, equivalent current dipoles, etc.) can be visualized in near real-time. Finally, this study demonstrates the accuracy and functionality of REST with data from two experiments and discusses some relevant applications.

Attentional Bias Predicts Increased Reward Salience and Risk Taking in Bipolar Disorder

BACKGROUND

There is amassing evidence that risky decision-making in bipolar disorder is related to reward-based differences in frontostriatal regions. However, the roles of early attentional and later cognitive processes remain unclear, limiting theoretical understanding and development of targeted interventions.

METHODS

Twenty euthymic bipolar disorder and 19 matched control participants played a Roulette task in which they won and lost money. Event-related potentials and source analysis were used to quantify predominantly sensory-attentional (N1), motivational salience (feedback-related negativities [FRN]), and cognitive appraisal (P300) stages of processing. We predicted that the bipolar disorder group would show increased N1, consistent with increased attentional orienting, and reduced FRN, consistent with a bias to perceive outcomes more favorably.

RESULTS

As predicted, the bipolar disorder group showed increased N1 and reduced FRN but no differences in P300. N1 amplitude was additionally associated with real-life risk taking, and N1 source activity was reduced in visual cortex but increased activity in precuneus, frontopolar, and premotor cortex, compared to those of controls.

CONCLUSIONS

These findings demonstrate an early attentional bias to reward that potentially drives risk taking by priming approach behavior and elevating reward salience in the frontostriatal pathway. Although later cognitive appraisals of these inputs may be relatively intact in remission, interventions targeting attention orienting may also be effective in long-term reduction of relapse.

Electrical source imaging of interictal spikes using multiple sparse volumetric priors for presurgical epileptogenic focus localization

Electrical source imaging of interictal spikes observed in EEG recordings of patients with refractory epilepsy provides useful information to localize the epileptogenic focus during the presurgical evaluation. However, the selection of the time points or time epochs of the spikes in order to estimate the origin of the activity remains a challenge. In this study, we consider a Bayesian EEG source imaging technique for distributed sources, i.e. the multiple volumetric sparse priors (MSVP) approach. The approach allows to estimate the time courses of the intensity of the sources corresponding with a specific time epoch of the spike. Based on presurgical averaged interictal spikes in six patients who were successfully treated with surgery, we estimated the time courses of the source intensities for three different time epochs: (i) an epoch starting 50 ms before the spike peak and ending at 50% of the spike peak during the rising phase of the spike, (ii) an epoch starting 50 ms before the spike peak and ending at the spike peak and (iii) an epoch containing the full spike time period starting 50 ms before the spike peak and ending 230 ms after the spike peak. To identify the primary source of the spike activity, the source with the maximum energy from 50 ms before the spike peak till 50% of the spike peak was subsequently selected for each of the time windows. For comparison, the activity at the spike peaks and at 50% of the peaks was localized using the LORETA inversion technique and an ECD approach. Both patient-specific spherical forward models and patient-specific 5-layered finite difference models were considered to evaluate the influence of the forward model. Based on the resected zones in each of the patients, extracted from post-operative MR images, we compared the distances to the resection border of the estimated activity. Using the spherical models, the distances to the resection border for the MSVP approach and each of the different time epochs were in the same range as the LORETA and ECD techniques. We found distances smaller than 23 mm, with robust results for all the patients. For the finite difference models, we found that the distances to the resection border for the MSVP inversions of the full spike time epochs were generally smaller compared to the MSVP inversions of the time epochs before the spike peak. The results also suggest that the inversions using the finite difference models resulted in slightly smaller distances to the resection border compared to the spherical models. The results we obtained are promising because the MSVP approach allows to study the network of the estimated source-intensities and allows to characterize the spatial extent of the underlying sources.

•

A Bayesian ESI technique is evaluated to localize interictal spike activity.

•

Averaged spikes in six patients were used that were seizure free after surgery.

•

We compared the technique with the LORETA an ECD technique.

•

We evaluated both spherical and 5-layered finite difference forward models.

•

Our approach is potentially useful to delineate the irritative zone.

Electrical source localization by LORETA in patients with epilepsy: Confirmation by postoperative MRI

BACKGROUND

Few studies have been conducted that have compared electrical source localization (ESL) results obtained by analyzing ictal patterns in scalp electroencephalogram (EEG) with the brain areas that are found to be responsible for seizures using other brain imaging techniques. Additionally, adequate studies have not been performed to confirm the accuracy of ESL methods.

MATERIALS AND METHODS

In this study, ESL was conducted using LORETA (Low Resolution Brain Electromagnetic Tomography) in 9 patients with lesions apparent on magnetic resonance imaging (MRI) and in 6 patients who did not exhibit lesions on their MRIs. EEGs of patients who underwent surgery for epilepsy and had follow-ups for at least 1 year after operations were analyzed for ictal spike, rhythmic, paroxysmal fast, and obscured EEG activities. Epileptogenic zones identified in postoperative MRIs were then compared with localizations obtained by LORETA model we employed.

RESULTS

We found that brain areas determined via ESL were in concordance with resected brain areas for 13 of the 15 patients evaluated, and those 13 patients were post-operatively determined as being seizure-free.

CONCLUSION

ESL, which is a noninvasive technique, may contribute to the correct delineation of epileptogenic zones in patients who will eventually undergo surgery to treat epilepsy, (regardless of neuroimaging status). Moreover, ESL may aid in deciding on the number and localization of intracranial electrodes to be used in patients who are candidates for invasive recording.

Real-Time Adaptive EEG Source Separation Using Online Recursive Independent Component Analysis

Independent component analysis (ICA) has been widely applied to electroencephalographic (EEG) biosignal processing and brain-computer interfaces. The practical use of ICA, however, is limited by its computational complexity, data requirements for convergence, and assumption of data stationarity, especially for high-density data. Here we study and validate an optimized online recursive ICA algorithm (ORICA) with online recursive least squares (RLS) whitening for blind source separation of high-density EEG data, which offers instantaneous incremental convergence upon presentation of new data. Empirical results of this study demonstrate the algorithm's: 1) suitability for accurate and efficient source identification in high-density (64-channel) realistically-simulated EEG data; 2) capability to detect and adapt to nonstationarity in 64-ch simulated EEG data; and 3) utility for rapidly extracting principal brain and artifact sources in real 61-channel EEG data recorded by a dry and wearable EEG system in a cognitive experiment. ORICA was implemented as functions in BCILAB and EEGLAB and was integrated in an open-source Real-time EEG Source-mapping Toolbox (REST), supporting applications in ICA-based online artifact rejection, feature extraction for real-time biosignal monitoring in clinical environments, and adaptable classifications in brain-computer interfaces.

Bayesian model reduction and empirical Bayes for group (DCM) studies

This technical note describes some Bayesian procedures for the analysis of group studies that use nonlinear models at the first (within-subject) level - e.g., dynamic causal models - and linear models at subsequent (between-subject) levels. Its focus is on using Bayesian model reduction to finesse the inversion of multiple models of a single dataset or a single (hierarchical or empirical Bayes) model of multiple datasets. These applications of Bayesian model reduction allow one to consider parametric random effects and make inferences about group effects very efficiently (in a few seconds). We provide the relatively straightforward theoretical background to these procedures and illustrate their application using a worked example. This example uses a simulated mismatch negativity study of schizophrenia. We illustrate the robustness of Bayesian model reduction to violations of the (commonly used) Laplace assumption in dynamic causal modelling and show how its recursive application can facilitate both classical and Bayesian inference about group differences. Finally, we consider the application of these empirical Bayesian procedures to classification and prediction.

Temporal dynamics of contingency extraction from tonal and verbal auditory sequences

Consecutive sound events are often to some degree predictive of each other. Here we investigated the brain's capacity to detect contingencies between consecutive sounds by means of electroencephalography (EEG) during passive listening. Contingencies were embedded either within tonal or verbal stimuli. Contingency extraction was measured indirectly via the elicitation of the mismatch negativity (MMN) component of the event-related potential (ERP) by contingency violations. MMN results indicate that structurally identical forms of predictability can be extracted from both tonal and verbal stimuli. We also found similar generators to underlie the processing of contingency violations across stimulus types, as well as similar performance in an active-listening follow-up test. However, the process of passive contingency extraction was considerably slower (twice as many rule exemplars were needed) for verbal than for tonal stimuli These results suggest caution in transferring findings on complex predictive regularity processing obtained with tonal stimuli directly to the speech domain.