Evaluation of Electroencephalography Source Localization Algorithms with Multiple Cortical Sources

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Electrical Source Imaging of Somatosensory Evoked Potentials from Intracranial EEG Signals

Stereoelectroencephalography (SEEG) records electrical brain activity with intracerebral electrodes. However, it has an inherently limited spatial coverage. Electrical source imaging (ESI) infers the position of the neural generators from the recorded electric potentials, and thus, could overcome this spatial undersampling problem. Here, we aimed to quantify the accuracy of SEEG ESI under clinical conditions. We measured the somatosensory evoked potential (SEP) in SEEG and in high-density EEG (HD-EEG) in 20 epilepsy surgery patients. To localize the source of the SEP, we employed standardized low resolution brain electromagnetic tomography (sLORETA) and equivalent current dipole (ECD) algorithms. Both sLORETA and ECD converged to similar solutions. Reflecting the large differences in the SEEG implantations, the localization error also varied in a wide range from 0.4 to 10 cm. The SEEG ESI localization error was linearly correlated with the distance from the putative neural source to the most activated contact. We show that it is possible to obtain reliable source reconstructions from SEEG under realistic clinical conditions, provided that the high signal fidelity recording contacts are sufficiently close to the source of the brain activity.

Speaker-induced Suppression in EEG during a Naturalistic Reading and Listening Task

Speaking elicits a suppressed neural response when compared with listening to others' speech, a phenomenon known as speaker-induced suppression (SIS). Previous research has focused on investigating SIS at constrained levels of linguistic representation, such as the individual phoneme and word level. Here, we present scalp EEG data from a dual speech perception and production task where participants read sentences aloud then listened to playback of themselves reading those sentences. Playback was separated into immediate repetition of the previous trial and randomized repetition of a former trial to investigate if forward modeling of responses during passive listening suppresses the neural response. Concurrent EMG was recorded to control for movement artifact during speech production. In line with previous research, ERP analyses at the sentence level demonstrated suppression of early auditory components of the EEG for production compared with perception. To evaluate whether linguistic abstractions (in the form of phonological feature tuning) are suppressed during speech production alongside lower-level acoustic information, we fit linear encoding models that predicted scalp EEG based on phonological features, EMG activity, and task condition. We found that phonological features were encoded similarly between production and perception. However, this similarity was only observed when controlling for movement by using the EMG response as an additional regressor. Our results suggest that SIS operates at a sensory representational level and is dissociated from higher order cognitive and linguistic processing that takes place during speech perception and production. We also detail some important considerations when analyzing EEG during continuous speech production.

Identifying good practices for detecting inter-regional linear functional connectivity from EEG

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

Source space connectomics of neurodegeneration: One-metric approach does not fit all

Brain functional connectivity in dementia has been assessed with dissimilar EEG connectivity metrics and estimation procedures, thereby increasing results' heterogeneity. In this scenario, joint analyses integrating information from different metrics may allow for a more comprehensive characterization of brain functional interactions in different dementia subtypes. To test this hypothesis, resting-state electroencephalogram (rsEEG) was recorded in individuals with Alzheimer's Disease (AD), behavioral variant frontotemporal dementia (bvFTD), and healthy controls (HCs). Whole-brain functional connectivity was estimated in the EEG source space using 101 different types of functional connectivity, capturing linear and nonlinear interactions in both time and frequency-domains. Multivariate machine learning and progressive feature elimination was run to discriminate AD from HCs, and bvFTD from HCs, based on joint analyses of i) EEG frequency bands, ii) complementary frequency-domain metrics (e.g., instantaneous, lagged, and total connectivity), and iii) time-domain metrics with different linearity assumption (e.g., Pearson correlation coefficient and mutual information). <10% of all possible connections were responsible for the differences between patients and controls, and atypical connectivity was never captured by >1/4 of all possible connectivity measures. Joint analyses revealed patterns of hypoconnectivity (patientsHCs) in both groups was mainly identified in frontotemporal regions. These atypicalities were differently captured by frequency- and time-domain connectivity metrics, in a bandwidth-specific fashion. The multi-metric representation of source space whole-brain functional connectivity evidenced the inadequacy of single-metric approaches, and resulted in a valid alternative for the selection problem in EEG connectivity. These joint analyses reveal patterns of brain functional interdependence that are overlooked with single metrics approaches, contributing to a more reliable and interpretable description of atypical functional connectivity in neurodegeneration.

Minimizing the distortions in electrophysiological source imaging of cortical oscillatory activity via Spectral Structured Sparse Bayesian Learning

Oscillatory processes at all spatial scales and on all frequencies underpin brain function. Electrophysiological Source Imaging (ESI) is the data-driven brain imaging modality that provides the inverse solutions to the source processes of the EEG, MEG, or ECoG data. This study aimed to carry out an ESI of the source cross-spectrum while controlling common distortions of the estimates. As with all ESI-related problems under realistic settings, the main obstacle we faced is a severely ill-conditioned and high-dimensional inverse problem. Therefore, we opted for Bayesian inverse solutions that posited a priori probabilities on the source process. Indeed, rigorously specifying both the likelihoods and a priori probabilities of the problem leads to the proper Bayesian inverse problem of cross-spectral matrices. These inverse solutions are our formal definition for cross-spectral ESI (cESI), which requires a priori of the source cross-spectrum to counter the severe ill-condition and high-dimensionality of matrices. However, inverse solutions for this problem were NP-hard to tackle or approximated within iterations with bad-conditioned matrices in the standard ESI setup. We introduce cESI with a joint a priori probability upon the source cross-spectrum to avoid these problems. cESI inverse solutions are low-dimensional ones for the set of random vector instances and not random matrices. We achieved cESI inverse solutions through the variational approximations via our Spectral Structured Sparse Bayesian Learning (ssSBL) algorithm https://github.com/CCC-members/Spectral-Structured-Sparse-Bayesian-Learning. We compared low-density EEG (10-20 system) ssSBL inverse solutions with reference cESIs for two experiments: (a) high-density MEG that were used to simulate EEG and (b) high-density macaque ECoG that were recorded simultaneously with EEG. The ssSBL resulted in two orders of magnitude with less distortion than the state-of-the-art ESI methods. Our cESI toolbox, including the ssSBL method, is available at https://github.com/CCC-members/BC-VARETA_Toolbox.

Effect of channel density, inverse solutions and connectivity measures on EEG resting-state networks reconstruction: A simulation study

Along with the study of brain activity evoked by external stimuli, the past two decades witnessed an increased interest in characterizing the spontaneous brain activity occurring during resting conditions. The identification of connectivity patterns in this so-called "resting-state" has been the subject of a great number of electrophysiology-based studies, using the Electro/Magneto-Encephalography (EEG/MEG) source connectivity method. However, no consensus has been reached yet regarding a unified (if possible) analysis pipeline, and several involved parameters and methods require cautious tuning. This is particularly challenging when different analytical choices induce significant discrepancies in results and drawn conclusions, thereby hindering the reproducibility of neuroimaging research. Hence, our objective in this study was to shed light on the effect of analytical variability on outcome consistency by evaluating the implications of parameters involved in the EEG source connectivity analysis on the accuracy of resting-state networks (RSNs) reconstruction. We simulated, using neural mass models, EEG data corresponding to two RSNs, namely the default mode network (DMN) and dorsal attentional network (DAN). We investigated the impact of five channel densities (19, 32, 64, 128, 256), three inverse solutions (weighted minimum norm estimate (wMNE), exact low-resolution brain electromagnetic tomography (eLORETA), and linearly constrained minimum variance (LCMV) beamforming) and four functional connectivity measures (phase-locking value (PLV), phase-lag index (PLI), and amplitude envelope correlation (AEC) with and without source leakage correction), on the correspondence between reconstructed and reference networks. We showed that, with different analytical choices related to the number of electrodes, source reconstruction algorithm, and functional connectivity measure, high variability is present in the results. More specifically, our results show that a higher number of EEG channels significantly increased the accuracy of the reconstructed networks. Additionally, our results showed significant variability in the performance of the tested inverse solutions and connectivity measures. Such methodological variability and absence of analysis standardization represent a critical issue for neuroimaging studies that should be prioritized. We believe that this work could be useful for the field of electrophysiology connectomics, by increasing awareness regarding the challenge of variability in methodological approaches and its implications on reported results.

Acoustic deprivation modulates central gain in human auditory brainstem and cortex

Beyond reduced audibility, there is convincing evidence that the auditory system adapts according to the principles of homeostatic plasticity in response to a hearing loss. Such compensatory changes include modulation of central auditory gain mechanisms. Earplugging is a common experimental method that has been used to introduce a temporary, reversible hearing loss that induces changes consistent with central gain modulation. In the present study, young, normal-hearing adult participants wore a unilateral earplug for two weeks, during which we measured changes in the acoustic reflex threshold (ART), loudness perception, and cortically-evoked (40 Hz) auditory steady-state response (ASSR) to assess potential modulation in central gain with reduced peripheral input. The ART decreased on average by 8 to 10 dB during the treatment period, with modest increases in loudness perception after one week but not after two weeks of earplug use. Significant changes in both the magnitude and hemispheric laterality of source-localized cortical ASSR measures revealed asymmetrical changes in stimulus-driven cortical activity over time. The ART results following unilateral earplugging are consistent with the literature and suggest that homeostatic plasticity is evident in the brainstem. The novel findings from the cortical ASSR in the present study indicates that reduced peripheral input induces adaptive homeostatic plasticity reflected as both an increase in central gain in the auditory brainstem and reduced cortical activity ipsilateral to the deprived ear. Both the ART and the novel use of the 40-Hz ASSR provide sensitive measures of central gain modulation in the brainstem and cortex of young, normal hearing listeners, and thus may be useful in future studies with other clinical populations.

Multiple functions of the angular gyrus at high temporal resolution

Here, the functions of the angular gyrus (AG) are evaluated in the light of current evidence from transcranial magnetic/electric stimulation (TMS/TES) and EEG/MEG studies. 65 TMS/TES and 52 EEG/MEG studies were examined in this review. TMS/TES literature points to a causal role in semantic processing, word and number processing, attention and visual search, self-guided movement, memory, and self-processing. EEG/MEG studies reported AG effects at latencies varying between 32 and 800 ms in a wide range of domains, with a high probability to detect an effect at 300-350 ms post-stimulus onset. A three-phase unifying model revolving around the process of sensemaking is then suggested: (1) early AG involvement in defining the current context, within the first 200 ms, with a bias toward the right hemisphere; (2) attention re-orientation and retrieval of relevant information within 200-500 ms; and (3) cross-modal integration at late latencies with a bias toward the left hemisphere. This sensemaking process can favour accuracy (e.g. for word and number processing) or plausibility (e.g. for comprehension and social cognition). Such functions of the AG depend on the status of other connected regions. The much-debated semantic role is also discussed as follows: (1) there is a strong TMS/TES evidence for a causal semantic role, (2) current EEG/MEG evidence is however weak, but (3) the existing arguments against a semantic role for the AG are not strong. Some outstanding questions for future research are proposed. This review recognizes that cracking the role(s) of the AG in cognition is possible only when its exact contributions within the default mode network are teased apart.

Dynamic auditory contributions to error detection revealed in the discrimination of Same and Different syllable pairs

During speech production auditory regions operate in concert with the anterior dorsal stream to facilitate online error detection. As the dorsal stream also is known to activate in speech perception, the purpose of the current study was to probe the role of auditory regions in error detection during auditory discrimination tasks as stimuli are encoded and maintained in working memory. A priori assumptions are that sensory mismatch (i.e., error) occurs during the discrimination of Different (mismatched) but not Same (matched) syllable pairs. Independent component analysis was applied to raw EEG data recorded from 42 participants to identify bilateral auditory alpha rhythms, which were decomposed across time and frequency to reveal robust patterns of event related synchronization (ERS; inhibition) and desynchronization (ERD; processing) over the time course of discrimination events. Results were characterized by bilateral peri-stimulus alpha ERD transitioning to alpha ERS in the late trial epoch, with ERD interpreted as evidence of working memory encoding via Analysis by Synthesis and ERS considered evidence of speech-induced-suppression arising during covert articulatory rehearsal to facilitate working memory maintenance. The transition from ERD to ERS occurred later in the left hemisphere in Different trials than in Same trials, with ERD and ERS temporally overlapping during the early post-stimulus window. Results were interpreted to suggest that the sensory mismatch (i.e., error) arising from the comparison of the first and second syllable elicits further processing in the left hemisphere to support working memory encoding and maintenance. Results are consistent with auditory contributions to error detection during both encoding and maintenance stages of working memory, with encoding stage error detection associated with stimulus concordance and maintenance stage error detection associated with task-specific retention demands.

Patient-specific solution of the electrocorticography forward problem in deforming brain

Invasive intracranial electroencephalography (iEEG), or electrocorticography (ECoG), measures electric potential directly on the surface of the brain and can be used to inform treatment planning for epilepsy surgery. Combined with numerical modeling it can further improve accuracy of epilepsy surgery planning. Accurate solution of the iEEG forward problem, which is a crucial prerequisite for solving the iEEG inverse problemin epilepsy seizure onset zone localization, requires accurate representation of the patient's brain geometry and tissue electrical conductivity after implantation of electrodes. However, implantation of subdural grid electrodes causes the brain to deform, which invalidates preoperatively acquired image data. Moreover, postoperative magnetic resonance imaging (MRI) is incompatible with implanted electrodes and computed tomography (CT) has insufficient range of soft tissue contrast, which precludes both MRI and CT from being used to obtain the deformed postoperative geometry. In this paper, we present a biomechanics-based image warping procedure using preoperative MRI for tissue classification and postoperative CT for locating implanted electrodes to perform non-rigid registration of the preoperative image data to the postoperative configuration. We solve the iEEG forward problem on the predicted postoperative geometry using the finite element method (FEM) which accounts for patient-specific inhomogeneity and anisotropy of tissue conductivity. Results for the simulation of a current source in the brain show large differences in electric potential predicted by the models based on the original images and the deformed images corresponding to the brain geometry deformed by placement of invasive electrodes. Computation of the lead field matrix (useful for solution of the iEEG inverse problem) also showed significant differences between the different models. The results suggest that rapid and accurate solution of the forward problem in a deformed brain for a given patient is achievable.

Age-related hearing loss is associated with alterations in temporal envelope processing in different neural generators along the auditory pathway

People with age-related hearing loss suffer from speech understanding difficulties, even after correcting for differences in hearing audibility. These problems are not only attributed to deficits in audibility but are also associated with changes in central temporal processing. The goal of this study is to obtain an understanding of potential alterations in temporal envelope processing for middle-aged and older persons with and without hearing impairment. The time series of activity of subcortical and cortical neural generators was reconstructed using a minimum-norm imaging technique. This novel technique allows for reconstructing a wide range of neural generators with minimal prior assumptions regarding the number and location of the generators. The results indicated that the response strength and phase coherence of middle-aged participants with hearing impairment (HI) were larger than for normal-hearing (NH) ones. In contrast, for the older participants, a significantly smaller response strength and phase coherence were observed in the participants with HI than the NH ones for most modulation frequencies. Hemispheric asymmetry in the response strength was also altered in middle-aged and older participants with hearing impairment and showed asymmetry toward the right hemisphere. Our brain source analyses show that age-related hearing loss is accompanied by changes in the temporal envelope processing, although the nature of these changes varies with age.

Methods Used to Estimate EEG Source-Space Networks: A Comparative Simulation-Based Study

Along with the study of the brain activity evoked by external stimuli, an important advance in current neuroscience involves understanding the spontaneous brain activity that occurs during resting conditions. Interestingly, the identification of the connectivity patterns in "resting-state" has been the subject of a great number of electrophysiology-based studies. In this context, the Electroencephalography (EEG) source connectivity method enables estimating resting-state cortical networks from scalp-EEG recordings. However, there is still no consensus over a unified pipeline adapted in all cases (e.g., type of task, a priori on studied networks) and numerous methodological questions remain unanswered. In order to address this problem, we simulated, using neural mass models, EEG data corresponding to the default mode network (DMN), the most widely studied resting-state network, and tested the effect of different channel densities, two inverse solutions and two functional connectivity measures on the correspondence between the reconstructed networks and the reference networks. Results showed that increasing the number of electrodes enhances the accuracy of the network reconstruction, and that eLORETA/PLV led to better accuracy than other inverse solution/connectivity measure combinations in terms of the correlation between reconstructed and reference connectivity matrices. This work has a wide range of implications in the field of electrophysiology connectomics, and is a step towards a convergence and standardization of approaches in this emerging field.

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.

Functional connectivity using high density EEG shows competitive reliability and agreement across test/retest sessions

BACKGROUND

Electrophysiological resting state functional connectivity using high density electroencephalography (hdEEG) is gaining momentum. The increased resolution offered by hdEEG, usually either 128 or 256 channels, permits source localization of EEG signals on the cortical surface. However, the number of methodological options for the acquisition and analysis of resting state hdEEG is extremely large. These include acquisition duration, eyes open/closed, channel density, source localization methods, and functional connectivity metric.

NEW METHODS

We undertake an extensive examination of the test-retest reliability and methodological agreement of all these options for regional measures of functional connectivity.

RESULTS

Power envelope connectivity shows larger test-retest reliability than imaginary coherence across all bands. While channel density doesn't strongly impact reliability or agreement, source localization methods produce systematically different functional connectivity, highlighting an important obstacle for replicating results in the literature. Most importantly, reliability and agreement often plateaus at or after 6 minutes of acquisition, well beyond the typical duration of 3 minutes. Finally, our study demonstrates that resting EEG can be as or more reliable than resting fMRI acquired in the same individuals.

CONCLUSIONS

The competitive reliability and agreement of power envelope connectivity greatly increases our confidence in measuring resting state connectivity using EEG and its capacity to find individual differences.

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.

Mean-Field Modeling of Brain-Scale Dynamics for the Evaluation of EEG Source-Space Networks

Understanding the dynamics of brain-scale functional networks at rest and during cognitive tasks is the subject of intense research efforts to unveil fundamental principles of brain functions. To estimate these large-scale brain networks, the emergent method called "electroencephalography (EEG) source connectivity" has generated increasing interest in the network neuroscience community, due to its ability to identify cortical brain networks with satisfactory spatio-temporal resolution, while reducing mixing and volume conduction effects. However, no consensus has been reached yet regarding a unified EEG source connectivity pipeline, and several methodological issues have to be carefully accounted to avoid pitfalls. Thus, a validation toolbox that provides flexible "ground truth" models is needed for an objective methods/parameters evaluation and, thereby an optimization of the EEG source connectivity pipeline. In this paper, we show how a recently developed large-scale model of brain-scale activity, named COALIA, can provide to some extent such ground truth by providing realistic simulations of source-level and scalp-level activity. Using a bottom-up approach, the model bridges cortical micro-circuitry and large-scale network dynamics. Here, we provide an example of the potential use of COALIA to analyze, in the context of epileptiform activity, the effect of three key factors involved in the "EEG source connectivity" pipeline: (i) EEG sensors density, (ii) algorithm used to solve the inverse problem, and (iii) functional connectivity measure. Results showed that a high electrode density (at least 64 channels) is required to accurately estimate cortical networks. Regarding the inverse solution/connectivity measure combination, the best performance at high electrode density was obtained using the weighted minimum norm estimate (wMNE) combined with the weighted phase lag index (wPLI). Although those results are specific to the considered aforementioned context (epileptiform activity), we believe that this model-based approach can be successfully applied to other experimental questions/contexts. We aim at presenting a proof-of-concept of the interest of COALIA in the network neuroscience field, and its potential use in optimizing the EEG source-space network estimation pipeline.

Stimulus transformation into motor action: Dynamic graph analysis reveals a posterior-to-anterior shift in brain network communication of older subjects

Cognitive performance slows down with increasing age. This includes cognitive processes that are essential for the performance of a motor act, such as the slowing down in response to an external stimulus. The objective of this study was to identify aging-associated functional changes in the brain networks that are involved in the transformation of external stimuli into motor action. To investigate this topic, we employed dynamic graphs based on phase-locking of Electroencephalography signals recorded from healthy younger and older subjects while performing a simple visually-cued finger-tapping task. The network analysis yielded specific age-related network structures varying in time in the low frequencies (2-7 Hz), which are closely connected to stimulus processing, movement initiation and execution in both age groups. The networks in older subjects, however, contained several additional, particularly interhemispheric, connections and showed an overall increased coupling density. Cluster analyses revealed reduced variability of the subnetworks in older subjects, particularly during movement preparation. In younger subjects, occipital, parietal, sensorimotor and central regions were-temporally arranged in this order-heavily involved in hub nodes. Whereas in older subjects, a hub in frontal regions preceded the noticeably delayed occurrence of sensorimotor hubs, indicating different neural information processing in older subjects. All observed changes in brain network organization, which are based on neural synchronization in the low frequencies, provide a possible neural mechanism underlying previous fMRI data, which report an overactivation, especially in the prefrontal and pre-motor areas, associated with a loss of hemispheric lateralization in older subjects.

Brain mapping of auditory steady-state responses: A broad view of cortical and subcortical sources

Auditory steady-state responses (ASSRs) are evoked brain responses to modulated or repetitive acoustic stimuli. Investigating the underlying neural generators of ASSRs is important to gain in-depth insight into the mechanisms of auditory temporal processing. The aim of this study is to reconstruct an extensive range of neural generators, that is, cortical and subcortical, as well as primary and non-primary ones. This extensive overview of neural generators provides an appropriate basis for studying functional connectivity. To this end, a minimum-norm imaging (MNI) technique is employed. We also present a novel extension to MNI which facilitates source analysis by quantifying the ASSR for each dipole. Results demonstrate that the proposed MNI approach is successful in reconstructing sources located both within (primary) and outside (non-primary) of the auditory cortex (AC). Primary sources are detected in different stimulation conditions (four modulation frequencies and two sides of stimulation), thereby demonstrating the robustness of the approach. This study is one of the first investigations to identify non-primary sources. Moreover, we show that the MNI approach is also capable of reconstructing the subcortical activities of ASSRs. Finally, the results obtained using the MNI approach outperform the group-independent component analysis method on the same data, in terms of detection of sources in the AC, reconstructing the subcortical activities and reducing computational load.

Use of common average reference and large-Laplacian spatial-filters enhances EEG signal-to-noise ratios in intrinsic sensorimotor activity

BACKGROUND

Oscillations in the resting-state scalp electroencephalogram (EEG) represent various intrinsic brain activities. One of the characteristic EEG oscillations is the sensorimotor rhythm (SMR)-with its arch-shaped waveform in alpha- and betabands-that reflect sensorimotor activity. The representation of sensorimotor activity by the SMR depends on the signal-to-noise ratio achieved by EEG spatial filters.

NEW METHOD

We employed simultaneous recording of EEG and functional magnetic resonance imaging, and 10-min resting-state brain activities were recorded in 19 healthy volunteers. To compare the EEG spatial-filtering methods commonly used for extracting sensorimotor cortical activities, we assessed nine different spatial-filters: a default reference of EEG amplifier system, a common average reference (CAR), small-, middle- and large-Laplacian filters, and four types of bipolar manners (C3-Cz, C3-F3, C3-P3, and C3-T7). We identified the brain region that correlated with the EEG-SMR power obtained after each spatial-filtering method was applied. Subsequently, we calculated the proportion of the significant voxels in the sensorimotor cortex as well as the sensorimotor occupancy in all significant regions to examine the sensitivity and specificity of each spatial-filter.

RESULTS

The CAR and large-Laplacian spatial-filters were superior at improving the signal-to-noise ratios for extracting sensorimotor activity from the EEG-SMR signal.

COMPARISON WITH EXISTING METHODS

Our results are consistent with the spatial-filter selection to extract task-dependent activation for better control of EEG-SMR-based interventions. Our approach has the potential to identify the optimal spatial-filter for EEG-SMR.

CONCLUSIONS

Evaluating spatial-filters for extracting spontaneous sensorimotor activity from the EEG is a useful procedure for constructing more effective EEG-SMR-based interventions.

Neural Generators Underlying Temporal Envelope Processing Show Altered Responses and Hemispheric Asymmetry Across Age

Speech understanding problems are highly prevalent in the aging population, even when hearing sensitivity is clinically normal. These difficulties are attributed to changes in central temporal processing with age and can potentially be captured by age-related changes in neural generators. The aim of this study is to investigate age-related changes in a wide range of neural generators during temporal processing in middle-aged and older persons with normal audiometric thresholds. A minimum-norm imaging technique is employed to reconstruct cortical and subcortical neural generators of temporal processing for different acoustic modulations. The results indicate that for relatively slow modulations (<50 Hz), the response strength of neural sources is higher in older adults than in younger ones, while the phase-locking does not change. For faster modulations (80 Hz), both the response strength and the phase-locking of neural sources are reduced in older adults compared to younger ones. These age-related changes in temporal envelope processing of slow and fast acoustic modulations are possibly due to loss of functional inhibition, which is accompanied by aging. Both cortical (primary and non-primary) and subcortical neural generators demonstrate similar age-related changes in response strength and phase-locking. Hemispheric asymmetry is also altered in older adults compared to younger ones. Alterations depend on the modulation frequency and side of stimulation. The current findings at source level could have important implications for the understanding of age-related changes in auditory temporal processing and for developing advanced rehabilitation strategies to address speech understanding difficulties in the aging population.

Hand Focused Upper Extremity Rehabilitation in the Subacute Phase Post-stroke Using Interactive Virtual Environments

Introduction: Innovative motor therapies have attempted to reduce upper extremity impairment after stroke but have not made substantial improvement as over 50% of people post-stroke continue to have sensorimotor deficits affecting their self-care and participation in daily activities. Intervention studies have focused on the role of increased dosing, however recent studies have indicated that timing of rehabilitation interventions may be as important as dosing and importantly, that dosing and timing interact in mediating effectiveness. This study is designed to empirically test dosing and timing.

Methods and Analysis: In this single-blinded, interventional study, subjects will be stratified on two dimensions, impairment level (Fugl-Meyer Upper Extremity Assessment (FM) and presence or absence of Motor Evoked Potentials (MEPs) as follows; (1) Severe, FM score 10–19, MEP+, (2) Severe, FM score 10–19, MEP–, (3) Moderate, FM score 20–49, MEP+, (4) Moderate, FM score 20–49, MEP–. Subjects not eligible for TMS will be assigned to either group 2 (if severe) or group 3 (if moderate). Stratified block randomization will then be used to achieve a balanced assignment. Early Robotic/VR Therapy (EVR) experimental group will receive in-patient usual care therapy plus an extra 10 h of intensive upper extremity therapy focusing on the hand using robotically facilitated rehabilitation interventions presented in virtual environments and initiated 5–30 days post-stroke. Delayed Robotic/VR Therapy (DVR) experimental group will receive the same intervention but initiated 30–60 days post-stroke. Dose-matched usual care group (DMUC) will receive an extra 10 h of usual care initiated 5–30 days post-stroke. Usual Care Group (UC) will receive the usual amount of physical/occupational therapy.

Outcomes: There are clinical, neurophysiological, and kinematic/kinetic measures, plus measures of daily arm use and quality of life. Primary outcome is the Action Research Arm Test (ARAT) measured at 4 months post-stroke.

Discussion: Outcome measures will be assessed to determine whether there is an early time period in which rehabilitation will be most effective, and whether there is a difference in the recapture of premorbid patterns of movement vs. the development of an efficient, but compensatory movement strategy.

Ethical Considerations: The IRBs of New Jersey Institute of Technology, Rutgers University, Northeastern University, and Kessler Foundation reviewed and approved all study protocols. Study was registered in https://ClinicalTrials.gov (NCT03569059) prior to recruitment. Dissemination will include submission to peer-reviewed journals and professional presentations.

sLORETA Source Localisation of Visual Mismatch Negativity in Dyslexic Children During Malay Orthographical Lexicon Stimulations

Background

While there are studies on visual lexical processing in other languages among dyslexics, no studies were done in the Malay language. The origin of visual lexical processing might be different in the Malay language. We aimed to detect the source localisation of visual mismatch negativity (vMMN) during Malay orthographic lexicon stimulations, employing an event-related potential (ERP) study.

Methods

Twelve dyslexic and twelve non-dyslexic children participated in this study. They pushed button '1' when they saw real (meaningful) Malay words and button '2' for pseudowords (meaningless). The source localisation of vMMN was performed in the grand average waveform by applying the standardised low-resolution brain electromagnetic tomography (sLORETA) method using Net Station software.

Results

Left occipital (BA17) and left temporal (BA37) lobes were activated during real words in the non-dyslexic and dyslexic children, respectively. During pseudowords, BA18 and BA17 areas of the left occipital lobe were activated in the non-dyslexic and dyslexic children, separately. vMMN sources were found at the left temporal (BA37) and right frontal (BA11) lobes in non-dyslexic and dyslexic children, respectively.

Conclusion

Right frontal lobe is the decision-making area where vMMN source was found in dyslexic children. We concluded that dyslexic children required the decision-making area to detect Malay real and pseudowords.

Intervention-induced changes in neural connectivity during motor preparation may affect cortical activity at motor execution

Effective interventions have demonstrated the ability to improve motor function by reengaging ipsilesional resources, which appears to be critical and feasible for hand function recovery even in individuals with severe chronic stroke. However, previous studies focus on changes in brain activity related to motor execution. How changes in motor preparation may facilitate these changes at motor execution is still unclear. To address this question, 8 individuals with severe chronic hemiparetic stroke participated in a device-assisted intervention for seven weeks. We then quantified changes in both coupling between regions during motor preparation and changes in topographical cortical activity at motor execution for both hand opening in isolation and together with the shoulder using high-density EEG. We hypothesized that intervention-induced changes in cortico-cortico interactions during motor preparation would lead to changes in activity at motor execution specifically towards an increased reliance on the ipsilesional hemisphere. In agreement with this hypothesis, we found that, following the intervention, individuals displayed a reduction in coupling from ipsilesional M1 to contralesional M1 within gamma frequencies during motor preparation for hand opening. This was followed by a reduction in activity in the contralesional primary sensorimotor cortex during motor execution. Similarly, during lifting and opening, a shift to negative coupling within ipsilesional M1 from gamma to beta frequencies was accompanied by an increase in ipsilesional primary sensorimotor cortex activity following the intervention. Together, these results show that intervention-induced changes in coupling within or between motor regions during motor preparation may affect cortical activity at execution.

Simultaneous human intracerebral stimulation and HD-EEG, ground-truth for source localization methods

Precisely localizing the sources of brain activity as recorded by EEG is a fundamental procedure and a major challenge for both research and clinical practice. Even though many methods and algorithms have been proposed, their relative advantages and limitations are still not well established. Moreover, these methods involve tuning multiple parameters, for which no principled way of selection exists yet. These uncertainties are emphasized due to the lack of ground-truth for their validation and testing. Here we present the Localize-MI dataset, which constitutes the first open dataset that comprises EEG recorded electrical activity originating from precisely known locations inside the brain of living humans. High-density EEG was recorded as single-pulse biphasic currents were delivered at intensities ranging from 0.1 to 5 mA through stereotactically implanted electrodes in diverse brain regions during pre-surgical evaluation of patients with drug-resistant epilepsy. The uses of this dataset range from the estimation of in vivo tissue conductivity to the development, validation and testing of forward and inverse solution methods.

Limited capacity for ipsilateral secondary motor areas to support hand function post-stroke

KEY POINTS

Ipsilateral-projecting corticobulbar pathways, originating primarily from secondary motor areas, innervate the proximal and even distal portions, although they branch more extensively at the spinal cord. It is currently unclear to what extent these ipsilateral secondary motor areas and subsequent cortical projections may contribute to hand function following stroke-induced damage to one hemisphere. In the present study, we provide both structural and functional evidence indicating that individuals increasingly rely on ipsilateral secondary motor areas, although at the detriment of hand function. Increased activity in ipsilateral secondary motor areas was associated with increased involuntary coupling between shoulder abduction and finger flexion, most probably as a result of the low resolution of these pathways, making it increasingly difficult to open the hand. These findings suggest that, although ipsilateral secondary motor areas may support proximal movements, they do not have the capacity to support distal hand function, particularly for hand opening.

ABSTRACT

Recent findings have shown connections of ipsilateral cortico-reticulospinal tract (CRST), predominantly originating from secondary motor areas to not only proximal, but also distal muscles of the arm. Following a unilateral stroke, CRST from the ipsilateral side remains intact and thus has been proposed as a possible backup system for post-stroke rehabilitation even for the hand. We argue that, although CRST from ipsilateral secondary motor areas can provide control for proximal joints, it is insufficient to control either hand or coordinated shoulder and hand movements as a result of its extensive spinal branching compared to contralateral corticospinal tract. To address this issue, we combined magnetic resonance imaging, high-density EEG, and robotics in 17 individuals with severe chronic hemiparetic stroke and 12 age-matched controls. We tested for changes in structural morphometry of the sensorimotor cortex and found that individuals with stroke demonstrated higher grey matter density in secondary motor areas ipsilateral to the paretic arm compared to controls. We then measured cortical activity when participants were attempting to generate hand opening either supported on a table or when lifting against a shoulder abduction load. The addition of shoulder abduction during hand opening increased reliance on ipsilateral secondary motor areas in stroke, but not controls. Crucially, the increased use of ipsilateral secondary motor areas was associated with decreased hand opening ability when lifting the arm as a result of involuntary coupling between the shoulder and wrist/finger flexors. Taken together, this evidence implicates a compensatory role for ipsilateral (i.e. contralesional) secondary motor areas post-stroke, although with no apparent capacity to support hand function.

More Reliable EEG Electrode Digitizing Methods Can Reduce Source Estimation Uncertainty, but Current Methods Already Accurately Identify Brodmann Areas

Electroencephalography (EEG) and source estimation can be used to identify brain areas activated during a task, which could offer greater insight on cortical dynamics. Source estimation requires knowledge of the locations of the EEG electrodes. This could be provided with a template or obtained by digitizing the EEG electrode locations. Operator skill and inherent uncertainties of a digitizing system likely produce a range of digitization reliabilities, which could affect source estimation and the interpretation of the estimated source locations. Here, we compared the reliabilities of five digitizing methods (ultrasound, structured-light 3D scan, infrared 3D scan, motion capture probe, and motion capture) and determined the relationship between digitization reliability and source estimation uncertainty, assuming other contributors to source estimation uncertainty were constant. We digitized a mannequin head using each method five times and quantified the reliability and validity of each method. We created five hundred sets of electrode locations based on our reliability results and applied a dipole fitting algorithm (DIPFIT) to perform source estimation. The motion capture method, which recorded the locations of markers placed directly on the electrodes had the best reliability with an average electrode variability of 0.001 cm. Then, in order of decreasing reliability were the method using a digitizing probe in the motion capture system, an infrared 3D scanner, a structured-light 3D scanner, and an ultrasound digitization system. Unsurprisingly, uncertainty of the estimated source locations increased with greater variability of EEG electrode locations and less reliable digitizing systems. If EEG electrode location variability was ∽1 cm, a single source could shift by as much as 2 cm. To help translate these distances into practical terms, we quantified Brodmann area accuracy for each digitizing method and found that the average Brodmann area accuracy for all digitizing methods was >80%. Using a template of electrode locations reduced the Brodmann area accuracy to ∽50%. Overall, more reliable digitizing methods can reduce source estimation uncertainty, but the significance of the source estimation uncertainty depends on the desired spatial resolution. For accurate Brodmann area identification, any of the digitizing methods tested can be used confidently.

Computationally Efficient Implementation of Aperture Domain Model Image Reconstruction

Aperture domain model image reconstruction (ADMIRE) is a useful tool to mitigate ultrasound imaging artifacts caused by acoustic clutter. However, its lengthy run-time impairs its usefulness. To overcome this drawback, we evaluated the reduced model methods with otherwise similar performance to ADMIRE. We also assessed other approaches to speed up ADMIRE, including the use of different levels of short-time Fourier transform (STFT) window overlap and examining the degrees of freedom of the model fit. In this study, we conducted an analysis of the reduced models, including those using Gram-Schmidt orthonormalization (GSO), singular value decomposition (SVD), and independent component analysis (ICA). We evaluated these reduced models using the data from simulations, experimental phantoms, and in vivo liver scans. We then tested ADMIRE's performance using full, GSO, SVD, and ICA-fourth-order blind identification (ICA-FOBI) models. The results from simulations, experimental phantoms, and in vivo data indicate that a model reduced using the ICA-FOBI method is the most promising for accelerating ADMIRE implementation. In the in vivo liver data, the improvements in contrast relative to delay-and-sum (DAS) using a full model, GSO, SVD, and ICA-FOBI models are 6.29 ± 0.25 dB, 11.88 ± 0.90 dB, 9.01 ± 0.67 dB, and 6.36 ± 0.27 dB, respectively; whereas, the contrast-to-noise ratio (CNR) improvement values in the same order are 0.04 ± 0.06 dB, -1.70 ± 0.17 dB, -1.51 ± 0.19 dB, and 0.12 ± 0.07 dB, respectively. The implementation of ADMIRE using the reduced model methods can decrease ADMIRE's computational complexity over three orders of magnitude. The use of a 50% STFT window overlap reduces ADMIRE's serial run time by more than one order of magnitude without any remarkable loss of image quality, when compared to the use of a 90% window overlap used previously. Based on these findings, a combination of the ICA-FOBI model and the use of a 50% STFT window overlap makes the ADMIRE algorithm more computationally efficient.

Quantitative Evaluation in Estimating Sources Underlying Brain Oscillations Using Current Source Density Methods and Beamformer Approaches

Brain oscillations from EEG and MEG shed light on neurophysiological mechanisms of human behavior. However, to extract information on cortical processing, researchers have to rely on source localization methods that can be very broadly classified into current density estimates such as exact low-resolution brain electromagnetic tomography (eLORETA), minimum norm estimates (MNE), and beamformers such as dynamic imaging of coherent sources (DICS) and linearly constrained minimum variance (LCMV). These algorithms produce a distributed map of brain activity underlying sustained and transient responses during neuroimaging studies of behavior. On the other hand, there are very few comparative analyses that evaluates the “ground truth detection” capabilities of these methods. The current article evaluates the reliability in estimation of sources of spectral event generators in the cortex using a two-pronged approach. First, simulated EEG data with point dipoles and distributed dipoles are used to validate the accuracy and sensitivity of each one of these methods of source localization. The abilities of the techniques were tested by comparing the localization error, focal width, false positive (FP) ratios while detecting already known location of neural activity generators under varying signal-to-noise ratios (SNRs). Second, empirical EEG data during auditory steady state responses (ASSRs) in human participants were used to compare the distributed nature of source localization. All methods were successful in recovery of point sources in favorable signal to noise scenarios and could achieve high hit rates if FPs are ignored. Interestingly, focal activation map is generated by LCMV and DICS when compared to eLORETA while control of FPs is much superior in eLORETA. Subsequently drawbacks and strengths of each method are highlighted with a detailed discussion on how to choose a technique based on empirical requirements.

Joint Encoding of Auditory Timing and Location in Visual Cortex

Co-occurring sounds can facilitate perception of spatially and temporally correspondent visual events. Separate lines of research have identified two putatively distinct neural mechanisms underlying two types of crossmodal facilitations: Whereas crossmodal phase resetting is thought to underlie enhancements based on temporal correspondences, lateralized occipital evoked potentials (ERPs) are thought to reflect enhancements based on spatial correspondences. Here, we sought to clarify the relationship between these two effects to assess whether they reflect two distinct mechanisms or, rather, two facets of the same underlying process. To identify the neural generators of each effect, we examined crossmodal responses to lateralized sounds in visually responsive cortex of 22 patients using electrocorticographic recordings. Auditory-driven phase reset and ERP responses in visual cortex displayed similar topography, revealing significant activity in pericalcarine, inferior occipital-temporal, and posterior parietal cortex, with maximal activity in lateral occipitotemporal cortex (potentially V5/hMT+). Laterality effects showed similar but less widespread topography. To test whether lateralized and nonlateralized components of crossmodal ERPs emerged from common or distinct neural generators, we compared responses throughout visual cortex. Visual electrodes responded to both contralateral and ipsilateral sounds with a contralateral bias, suggesting that previously observed laterality effects do not emerge from a distinct neural generator but rather reflect laterality-biased responses in the same neural populations that produce phase-resetting responses. These results suggest that crossmodal phase reset and ERP responses previously found to reflect spatial and temporal facilitation in visual cortex may reflect the same underlying mechanism. We propose a new unified model to account for these and previous results.

Electrophysiological correlates of top-down attentional modulation in olfaction

The capacity to pay attention is important for the cognitive ability, for example, evaluating an object for its qualities. Attention can selectively prioritize the neural processes that are relevant to a given task. Neuroimaging investigations on human attention are primarily focused on vision to the exclusion of other sensory systems, particularly olfaction. Neural underpinnings of human olfactory attention are still not clearly understood. Here, we combined electroencephalographic measurements of olfactory event related potential with electrical neuroimaging to investigate how the neural responses after inhaling the same odor differ between conditions with varying levels of attention, and, in which brain areas. We examined the neural responses when participants attended to a rose-like odor of phenylethyl alcohol for evaluating its pleasantness versus its passive inhalation. Our results gathered significant evidence for attentional modulation of the olfactory neural response. The most prominent effect was found for the late positive component, P3, of olfactory event related potential within a second from the odor onset. The source reconstruction of this data revealed activations in a distributed network of brain regions predominantly in inferior frontal cortex, insula, and  inferior temporal gyrus. These results suggest that the neuronal modulations from attention to olfactory pleasantness may be subserved by this network.

simBCI - A framework for studying BCI methods by simulated EEG

Brain-Computer Interface (BCI) methods are commonly studied using Electroencephalogram (EEG) data recorded from human experiments. For understanding and developing BCI signal processing techniques real data is costly to obtain and its composition is apriori unknown. The brain mechanisms generating the EEG are not directly observable and their states cannot be uniquely identified from the EEG. Subsequently, we do not have generative ground truth for real data. In this paper we propose a novel convenience framework called simBCI to alleviate testing and studying BCI signal processing methods in simulated, controlled conditions. The framework can be used to generate artificial BCI data and to test classification pipelines with such data. Models and parameters on both data generation and the signal processing side can be iterated over to examine the interplay of different combinations. The framework provides the first time open source implementations of several models and methods. We invite researchers to insert more advanced models. The proposed system does not intend to replace human experiments. Instead, it can be used to discover hypotheses, study algorithms, to educate about BCI, and to debug signal processing pipelines of other BCI systems. The proposed framework is modular, extensible and freely available as open source1. It currently requires Matlab.

Using Event-Related Potentials and Startle to Evaluate Time Course in Anxiety and Depression

The National Institute of Mental Health's Research Domain Criteria initiative is a research framework designed toward understanding psychopathology as abnormalities of dimensional neurobehavioral constructs rather than in terms of DSM-defined categories. Research Domain Criteria constructs within the negative valence domain are particularly relevant for understanding anxiety and depressive disorders, which are pervasive, debilitating, and characterized by negative processing bias. One important direction for Research Domain Criteria research is investigating processes and parameters related to the time course (or chronometry) of negative valenced constructs. Two reliable methods for assessing chronometry are event-related potentials (ERPs) and startle blink. In this qualitative review, we examine ERP and startle studies of individuals with anxiety or depression or individuals vulnerable to affective disorders. The aim of the review is to highlight how these methods can inform the role of chronometry in the spectrum of anxiety and depression. ERP studies examining different chronometry facets of negative valenced responses have shown that transdiagnostic groups of individuals with internalizing psychopathologies exhibit abnormalities at early stages of processing. Startle reactivity studies have robustly differentiated fear-based disorders (e.g., panic disorder, social phobia) from other anxiety disorders (e.g., generalized anxiety disorder) and have also shown that different internalizing phenotypes exhibit different patterns of habituation. Findings lend support to the value of ERP and startle measures in identifying groups that cut across conventional classification systems. We also highlight methodological issues that can aid in the validity and reproducibility of ERP and startle findings and, ultimately, in the goal of developing more precise models of anxiety and depression.

A Time-Varying Connectivity Analysis from Distributed EEG Sources: A Simulation Study

Time-varying connectivity analysis based on sources reconstructed using inverse modeling of electroencephalographic (EEG) data is important to understand the dynamic behaviour of the brain. We simulated cortical data from a visual spatial attention network with a time-varying connectivity structure, and then simulated the propagation to the scalp to obtain EEG data. Distributed EEG source modeling using sLORETA was applied. We compared different dipole (representing a source) selection strategies based on their time series in a region of interest. Next, we estimated multivariate autoregressive (MVAR) parameters using classical Kalman filter and general linear Kalman filter approaches followed by the calculation of partial directed coherence (PDC). MVAR parameters and PDC values for the selected sources were compared with the ground-truth. We found that the best strategy to extract the time series of a region of interest was to select a dipole with time series showing the highest correlation with the average time series in the region of interest. Dipole selection based on power or based on the largest singular value offer comparable alternatives. Among the different Kalman filter approaches, the use of a general linear Kalman filter was preferred to estimate PDC based connectivity except when only a small number of trials are available. In the latter case, the classical Kalman filter can be an alternative.

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.

Neural Plasticity in Moderate to Severe Chronic Stroke Following a Device-Assisted Task-Specific Arm/Hand Intervention

Currently, hand rehabilitation following stroke tends to focus on mildly impaired individuals, partially due to the inability for severely impaired subjects to sufficiently use the paretic hand. Device-assisted interventions offer a means to include this more severe population and show promising behavioral results. However, the ability for this population to demonstrate neural plasticity, a crucial factor in functional recovery following effective post-stroke interventions, remains unclear. This study aimed to investigate neural changes related to hand function induced by a device-assisted task-specific intervention in individuals with moderate to severe chronic stroke (upper extremity Fugl-Meyer < 30). We examined functional cortical reorganization related to paretic hand opening and gray matter (GM) structural changes using a multimodal imaging approach. Individuals demonstrated a shift in cortical activity related to hand opening from the contralesional to the ipsilesional hemisphere following the intervention. This was driven by decreased activity in contralesional primary sensorimotor cortex and increased activity in ipsilesional secondary motor cortex. Additionally, subjects displayed increased GM density in ipsilesional primary sensorimotor cortex and decreased GM density in contralesional primary sensorimotor cortex. These findings suggest that despite moderate to severe chronic impairments, post-stroke participants maintain ability to show cortical reorganization and GM structural changes following a device-assisted task-specific arm/hand intervention. These changes are similar as those reported in post-stroke individuals with mild impairment, suggesting that residual neural plasticity in more severely impaired individuals may have the potential to support improved hand function.

Nonlinear Coupling between Cortical Oscillations and Muscle Activity during Isotonic Wrist Flexion

Coupling between cortical oscillations and muscle activity facilitates neuronal communication during motor control. The linear part of this coupling, known as corticomuscular coherence, has received substantial attention, even though neuronal communication underlying motor control has been demonstrated to be highly nonlinear. A full assessment of corticomuscular coupling, including the nonlinear part, is essential to understand the neuronal communication within the sensorimotor system. In this study, we applied the recently developed n:m coherence method to assess nonlinear corticomuscular coupling during isotonic wrist flexion. The n:m coherence is a generalized metric for quantifying nonlinear cross-frequency coupling as well as linear iso-frequency coupling. By using independent component analysis (ICA) and equivalent current dipole source localization, we identify four sensorimotor related brain areas based on the locations of the dipoles, i.e., the contralateral primary sensorimotor areas, supplementary motor area (SMA), prefrontal area (PFA) and posterior parietal cortex (PPC). For all these areas, linear coupling between electroencephalogram (EEG) and electromyogram (EMG) is present with peaks in the beta band (15–35 Hz), while nonlinear coupling is detected with both integer (1:2, 1:3, 1:4) and non-integer (2:3) harmonics. Significant differences between brain areas is shown in linear coupling with stronger coherence for the primary sensorimotor areas and motor association cortices (SMA, PFA) compared to the sensory association area (PPC); but not for the nonlinear coupling. Moreover, the detected nonlinear coupling is similar to previously reported nonlinear coupling of cortical activity to somatosensory stimuli. We suggest that the descending motor pathways mainly contribute to linear corticomuscular coupling, while nonlinear coupling likely originates from sensory feedback.

Nonlinear Connectivity in the Human Stretch Reflex Assessed by Cross-Frequency Phase Coupling

Communication between neuronal populations is facilitated by synchronization of their oscillatory activity. Although nonlinearity has been observed in the sensorimotor system, its nonlinear connectivity has not been widely investigated yet. This study investigates nonlinear connectivity during the human stretch reflex based on neuronal synchronization. Healthy participants generated isotonic wrist flexion while receiving a periodic mechanical perturbation to the wrist. Using a novel cross-frequency phase coupling metric, we estimate directional nonlinear connectivity, including time delay, from the perturbation to brain and to muscle, as well as from brain to muscle. Nonlinear phase coupling is significantly stronger from the perturbation to the muscle than to the brain, with a shorter time delay. The time delay from the perturbation to the muscle is 33 ms, similar to the reported latency of the spinal stretch reflex at the wrist. Source localization of nonlinear phase coupling from the brain to the muscle suggests activity originating from the motor cortex, although its effect on the stretch reflex is weak. As such nonlinear phase coupling between the perturbation and muscle activity is dominated by the spinal reflex loop. This study provides new evidence of nonlinear neuronal synchronization in the stretch reflex at the wrist joint with respect to spinal and transcortical loops.