RELICA: a method for estimating the reliability of independent components

How does Independent Component Analysis Preprocessing Affect EEG Microstates?

Over recent years, electroencephalographic (EEG) microstates have been increasingly used to investigate, at a millisecond scale, the temporal dynamics of large-scale brain networks. By studying their topography and chronological sequence, microstates research has contributed to the understanding of the brain's functional organization at rest and its alteration in neurological or mental disorders. Artifact removal strategies, which differ from study to study, may alter microstates topographies and features, possibly reducing the generalizability and comparability of results across research groups. The aim of this work was therefore to test the reliability of the microstate extraction process and the stability of microstate features against different strategies of EEG data preprocessing with Independent Component Analysis (ICA) to remove artifacts embedded in the data. A normative resting state EEG dataset was used where subjects alternate eyes-open (EO) and eyes-closed (EC) periods. Four strategies were tested: (i) avoiding ICA preprocessing altogether, (ii) removing ocular artifacts only, (iii) removing all reliably identified physiological/non physiological artifacts, (iv) retaining only reliably identified brain ICs. Results show that skipping the removal of ocular artifacts affects the stability of microstate evaluation criteria, microstate topographies and greatly reduces the statistical power of EO/EC microstate features comparisons, however differences are not as prominent with more aggressive preprocessing. Provided a good-quality dataset is recorded, and ocular artifacts are removed, microstates topographies and features can capture brain-related physiological data and are robust to artifacts, independently of the level of preprocessing, paving the way to automatized microstate extraction pipelines.

Spatial and temporal patterns of brain neural activity mediating human thermal sensations

This study aimed to elucidate the spatial and temporal patterns of brain neural activity that are associated with cold and hot sensations. Participants (n = 20) sat in a controlled room with their eyes closed and received local thermal stimuli to the right fingers using a Peltier apparatus. The thermal stimuli were repeated 40 times using a paired-thermal stimulus paradigm, comprising a 15 s-reference stimulus (32 °C), followed by 10 s-conditioned stimuli (24 °C and 40 °C, cold and hot conditions, respectively), for which 15-channel electroencephalography (EEG) signals were continuously monitored. To identify the patterns of brain neural activity, an independent component (IC) analysis was applied to the preprocessed EEG data. The equivalent current dipole locations were estimated, followed by clustering of the ICs with a dipole residual variance of <15 %. Subsequently, event-related spectral perturbations were analyzed in each identified cluster to calculate the power changes across specific frequency ranges. The right precentral gyrus, precuneus, medial frontal gyrus, middle frontal gyrus, superior frontal gyrus, cuneus, cingulate gyrus, left precentral gyrus, middle occipital gyrus, and cingulate gyrus were activated in both cold and hot conditions. In most activated regions, EEG power temporal changes were observed across the frequency ranges and were different between the two conditions. These results may suggest that cold and hot sensations are processed through different temporal brain neural activity patterns in overlapping brain regions.

Electrophysiological markers of vestibular-mediated self-motion perception - A pilot study

Peripheral vestibular activation results in multi-level responses, from brainstem-mediated reflexes (e.g. vestibular ocular reflex - VOR) to perception of self-motion. While VOR responses indicate preserved vestibular peripheral and brainstem functioning, there are no automated measures of vestibular perception of self-motion - important since some patients with brain disconnection syndromes manifest a vestibular agnosia (intact VOR but impaired self-motion perception). Electroencephalography ('EEG') - may provide a surrogate marker of vestibular perception of self-motion. A related objective is obtaining an EEG marker of vestibular sensory signal processing, distinct from vestibular-motion perception. We performed a pilot study comparing EEG responses in the dark when healthy participants sat in a vibrationless computer-controlled motorised rotating chair moving at near threshold of self-motion perception, versus a second situation in which subjects sat in the chair at rest in the dark who could be induced (or not) into falsely perceiving self-motion. In both conditions subjects could perceive self-motion perception, but in the second there was no bottom-up reflex-brainstem activation. Time-frequency analyses showed: (i) alpha frequency band activity is linked to vestibular sensory-signal activation; and (ii) theta band activity is a marker of vestibular-mediated self-motion perception. Consistent with emerging animal data, our findings support the role of theta activity in the processing of self-motion perception.

Linking Dementia Pathology and Alteration in Brain Activation to Complex Daily Functional Decline During the Preclinical Dementia Stages: Protocol for a Prospective Observational Cohort Study

BACKGROUND

Progressive difficulty in performing everyday functional activities is a key diagnostic feature of dementia syndromes. However, not much is known about the neural signature of functional decline, particularly during the very early stages of dementia. Early intervention before overt impairment is observed offers the best hope of reducing the burdens of Alzheimer disease (AD) and other dementias. However, to justify early intervention, those at risk need to be detected earlier and more accurately. The decline in complex daily function (CdF) such as managing medications has been reported to precede impairment in basic activities of daily living (eg, eating and dressing).

OBJECTIVE

Our goal is to establish the neural signature of decline in CdF during the preclinical dementia period.

METHODS

Gait is central to many CdF and community-based activities. Hence, to elucidate the neural signature of CdF, we validated a novel electroencephalographic approach to measuring gait-related brain activation while participants perform complex gait-based functional tasks. We hypothesize that dementia-related pathology during the preclinical period activates a unique gait-related electroencephalographic (grEEG) pattern that predicts a subsequent decline in CdF.

RESULTS

We provide preliminary findings showing that older adults reporting CdF limitations can be characterized by a unique gait-related neural signature: weaker sensorimotor and stronger motor control activation. This subsample also had smaller brain volume and white matter hyperintensities in regions affected early by dementia and engaged in less physical exercise. We propose a prospective observational cohort study in cognitively unimpaired older adults with and without subclinical AD (plasma amyloid-β) and vascular (white matter hyperintensities) pathologies. We aim to (1) establish the unique grEEG activation as the neural signature and predictor of decline in CdF during the preclinical dementia period; (2) determine associations between dementia-related pathologies and incidence of the neural signature of CdF; and (3) establish associations between a dementia risk factor, physical inactivity, and the neural signature of CdF.

CONCLUSIONS

By establishing the clinical relevance and biological basis of the neural signature of CdF decline, we aim to improve prediction during the preclinical stages of ADs and other dementias. Our approach has important research and translational implications because grEEG protocols are relatively inexpensive and portable, and predicting CdF decline may have real-world benefits.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

DERR1-10.2196/56726.

Verification of Source Activations in a 3D Brain Model Using 'CLEVER' Algorithm for Mental Arithmetic Conditions

Background: Living conditions are becoming challenging day by day. Mental stress on individuals is increasing due to multiple reasons. As mental stress is a major cause of mental illness, it must be detected at the earliest to prevent serious conditions such as depression and anxiety.

PURPOSE

The focus of this study is to detect the exact location of the source which causes such damage. In this article, we analyse the mental conditions of subjects under a workload of performing mental arithmetic calculations for various frequency bands and plot the topography to understand the areas of active potentials.

METHODS

We propose a Novel Cluster Ensemble Verifier (CLEVER) algorithm, which combines two different techniques: clustering and source localisation. The proposed algorithm is highly efficient in identifying the exact location of the source. It is seen that the topographic plots of the independent component analysis (ICA), which has the maximum percentage of relative variance, correlates to the cluster generated. We are able to give the percentage-wise contribution of every component which is responsible for brain source activation with less time complexity.

RESULTS

Out of 72 subjects, in 67 subjects, 299 out of 433 components originate from the occipital and parietal areas of the brain with a maximum power of 43.5 µv2. As an example, the relative variance of one component is found to be contributing up to 74.03% to source activations. Clusters show similarity across the subjects in the parietal and occipital areas of the brain. The dataset used for experimentation is EEGMAT from Physionet's repository. The computation time for the algorithms is 17.6 ± 3.2 minutes. Conclusion: Findings show that during mental arithmetic calculations, both occipital and parietal areas of the brain are involved. As the data is acquired by orally mentioning the mathematical problem, subjects tend to visualise the numbers while finding the solution, which is reflected in the occipital area of the brain. CLEVER algorithm verifies the origin of the activity in the occipital and parietal areas of the brain.

The effect of age on alpha rhythms in the human brain derived from source localized resting-state electroencephalography

With increasing age, peak alpha frequency (PAF) is slowed, and alpha power is reduced during resting-states with eyes closed. These age-related changes are evident across the whole scalp but remained unclear at the source level. The purpose of this study was to determine whether age impacts the power and frequency of the dominant alpha rhythm equally across source generators or whether the impact of age varies across sources. A total of 28 young adults and 26 elderly adults were recruited. High-density EEG was recorded for 10 mins with eyes closed. Single dipoles for each independent component were localized and clustered based on their anatomical label, resulting in 36 clusters. Meta-analyses were then conducted to assess effect sizes for PAF and power at PAF for all 36 clusters. Subgroup analyses were then implemented for frontal, sensorimotor, parietal, temporal, and occipital regions. The results of the meta-analyses showed that the elderly group exhibited slower PAF and less power at PAF compared to the young group. Subgroup analyses revealed age effects on PAF in parietal (g = 0.38), temporal (g = 0.65), and occipital regions (g = 1.04), with the largest effects observed in occipital regions. For power at PAF, age effects were observed in sensorimotor (g = 0.84) and parietal regions (g = 0.80), with the sensorimotor region showing the largest effect. Our findings show that age-related slowing and attenuation of the alpha rhythm manifests differentially across cortical regions, with sensorimotor and occipital regions most susceptible to age effects.

One out of ten independent components shows flipped polarity with poorer data quality: EEG database study

Independent component analysis (ICA) is widely used today for scalp-recorded EEG analysis. One of the limitations of ICA-based analysis is polarity indeterminacy. It is not easy to find detailed documentations that explains engineering solutions of how the polarity indeterminacy is addressed in a given implementation. We investigated how it is implemented in the case of EEGLAB and also the relation between the outcome of the polarity determination and classification of independent components (ICs) in terms of the estimated nature of the sources (brain, muscle, eye, etc.) using an open database of n = 212 EEG dataset of resting state recordings. We found that (1) about 91% of ICs showed positive-dominant IC scalp topographies; (2) positive-dominant ICs were more associated with brain-originated signals; (3) positive-dominant ICs showed more radial (peaked at 10-30 degrees deviations from the radial axis) dipolar projection pattern with less residual variance from fitting the equivalent current dipole. In conclusion, using the EEGLAB's default ICA algorithm, one out of 10 ICs results in flipping its polarity to negative, which is associated with non-radial dipole orientation with higher residual variance. Thus, we determined EEGLAB biases toward positive polarity in decomposing high-quality brain ICs.

Cortico-muscular connectivity is modulated by passive and active Lokomat-assisted Gait

The effects of robotic-assisted gait (RAG) training, besides conventional therapy, on neuroplasticity mechanisms and cortical integration in locomotion are still uncertain. To advance our knowledge on the matter, we determined the involvement of motor cortical areas in the control of muscle activity in healthy subjects, during RAG with Lokomat, both with maximal guidance force (100 GF-passive RAG) and without guidance force (0 GF-active RAG) as customary in rehabilitation treatments. We applied a novel cortico-muscular connectivity estimation procedure, based on Partial Directed Coherence, to jointly study source localized EEG and EMG activity during rest (standing) and active/passive RAG. We found greater cortico-cortical connectivity, with higher path length and tendency toward segregation during rest than in both RAG conditions, for all frequency bands except for delta. We also found higher cortico-muscular connectivity in distal muscles during swing (0 GF), and stance (100 GF), highlighting the importance of direct supraspinal control to maintain balance, even when gait is supported by a robotic exoskeleton. Source-localized connectivity shows that this control is driven mainly by the parietal and frontal lobes. The involvement of many cortical areas also in passive RAG (100 GF) justifies the use of the 100 GF RAG training for neurorehabilitation, with the aim of enhancing cortical-muscle connections and driving neural plasticity in neurological patients.

Exploring Spatio-Spectral Electroencephalogram Modulations of Imbuing Emotional Intent During Active Piano Playing

Imbuing emotional intent serves as a crucial modulator of music improvisation during active musical instrument playing. However, most improvisation-related neural endeavors have been gained without considering the emotional context. This study attempts to exploit reproducible spatio-spectral electroencephalogram (EEG) oscillations of emotional intent using a data-driven independent component analysis framework in an ecological multiday piano playing experiment. Through the four-day 32-ch EEG dataset of 10 professional players, we showed that EEG patterns were substantially affected by both intra- and inter-individual variability underlying the emotional intent of the dichotomized valence (positive vs. negative) and arousal (high vs. low) categories. Less than half (3-4) of the 10 participants analogously exhibited day-reproducible ( ≥ three days) spectral modulations at the right frontal beta in response to the valence contrast as well as the frontal central gamma and the superior parietal alpha to the arousal counterpart. In particular, the frontal engagement facilitates a better understanding of the frontal cortex (e.g., dorsolateral prefrontal cortex and anterior cingulate cortex) and its role in intervening emotional processes and expressing spectral signatures that are relatively resistant to natural EEG variability. Such ecologically vivid EEG findings may lead to better understanding of the development of a brain-computer music interface infrastructure capable of guiding the training, performance, and appreciation for emotional improvisatory status or actuating music interaction via emotional context.

Feasibility of Using Source-Level Brain Computer Interface for People with Multiple Sclerosis

This work evaluates the feasibility of using a source level Brain-Computer Interface (BCI) for people with Multiple Sclerosis (MS). Data used was previously collected EEG of eight participants (one participant with MS and seven neurotypical participants) who performed imagined movement of the right and left hand. Equivalent current dipole cluster fitting was used to assess related brain activity at the source level and assessed using dipole location and power spectrum analysis. Dipole clusters were resolved within the motor cortices with some notable spatial difference between the MS and control participants. Neural sources that generate motor imagery originated from similar motor areas in the participant with MS compared to the neurotypical participants. Power spectral analysis indicated a reduced level of alpha power in the participant with MS during imagery tasks compared to neurotypical participants. Power in the beta band may be used to distinguish between left and right imagined movement for users with MS in BCI applications.Clinical Relevance- This paper demonstrates the cortical areas activated during imagined BCI-type tasks in a participant with Multiple Sclerosis (MS), and is a proof of concept for translating BCI research to potential users with MS.

Microsynt: exploring the syntax of EEG microstates

Microstates represent electroencephalographic (EEG) activity as a sequence of switching, transient, metastable states. Growing evidence suggests the useful information on brain states is to be found in the higher-order temporal structure of these sequences. Instead of focusing on transition probabilities, here we propose "Microsynt", a method designed to highlight higher-order interactions that form a preliminary step towards understanding the syntax of microstate sequences of any length and complexity. Microsynt extracts an optimal vocabulary of "words" based on the length and complexity of the full sequence of microstates. Words are then sorted into classes of entropy and their representativeness within each class is statistically compared with surrogate and theoretical vocabularies. We applied the method on EEG data previously collected from healthy subjects undergoing propofol anaesthesia, and compared their "fully awake" (BASE) and "fully unconscious" (DEEP) conditions. Results show that microstate sequences, even at rest, are not random but tend to behave in a more predictable way, favoring simpler sub-sequences, or "words". Contrary to high-entropy words, lowest-entropy binary microstate loops are prominent and favored on average 10 times more than what is theoretically expected. Progressing from BASE to DEEP, the representation of low-entropy words increases while that of high-entropy words decreases. During the awake state, sequences of microstates tend to be attracted towards "A - B - C" microstate hubs, and most prominently A - B binary loops. Conversely, with full unconsciousness, sequences of microstates are attracted towards "C - D - E" hubs, and most prominently C - E binary loops, confirming the putative relation of microstates A and B to externally-oriented cognitive processes and microstate C and E to internally-generated mental activity. Microsynt can form a syntactic signature of microstate sequences that can be used to reliably differentiate two or more conditions.

Hybrid Genetic Algorithm for Clustering IC Topographies of EEGs

Clustering of independent component (IC) topographies of Electroencephalograms (EEG) is an effective way to find brain-generated IC processes associated with a population of interest, particularly for those cases where event-related potential features are not available. This paper proposes a novel algorithm for the clustering of these IC topographies and compares its results with the most currently used clustering algorithms. In this study, 32-electrode EEG signals were recorded at a sampling rate of 500 Hz for 48 participants. EEG signals were pre-processed and IC topographies computed using the AMICA algorithm. The algorithm implements a hybrid approach where genetic algorithms are used to compute more accurate versions of the centroids and the final clusters after a pre-clustering phase based on spectral clustering. The algorithm automatically selects the optimum number of clusters by using a fitness function that involves local-density along with compactness and separation criteria. Specific internal validation metrics adapted to the use of the absolute correlation coefficient as the similarity measure are defined for the benchmarking process. Assessed results across different ICA decompositions and groups of subjects show that the proposed clustering algorithm significantly outperforms the (baseline) clustering algorithms provided by the software EEGLAB, including CORRMAP.

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

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

Neural signature of mobility-related everyday function in older adults at-risk of cognitive impairment

Assessment of everyday activities is central to the diagnosis of dementia. Yet, little is known about brain processes associated with everyday functional limitations, particularly during early stages of cognitive decline. Twenty-six older adults (mean = 74.9 y) were stratified by risk using the Montreal Cognitive Assessment battery (MoCA, range: 0- 30) to classify individuals as higher (22-26) and lower risk (27+) of cognitive impairment. We investigated everyday function using a gait task designed to destabilize posture and applied Mobile Brain/Body Imaging. We predicted that participants would increase step width to gain stability, yet the underlying neural signatures would be different for lower versus higher risk individuals. Step width and fronto-parietal activation increased during visually perturbed input. Frontomedial theta increased in higher risk individuals during perturbed and unperturbed inputs. Left sensorimotor beta decreased in lower risk individuals during visually perturbed input. Modulations in theta and beta power were associated with MoCA scores. Our findings suggest that older adults at-risk of cognitive impairment can be characterized by a unique neural signature of everyday function.

Rapid Interactions of Widespread Brain Networks Characterize Semantic Cognition

Language comprehension requires the rapid retrieval and integration of contextually appropriate concepts ("semantic cognition"). Current neurobiological models of semantic cognition are limited by the spatial and temporal restrictions of single-modality neuroimaging and lesion approaches. This is a major impediment given the rapid sequence of processing steps that have to be coordinated to accurately comprehend language. Through the use of fused functional magnetic resonance imaging and electroencephalography analysis in humans (n = 26 adults; 15 females), we elucidate a temporally and spatially specific neurobiological model for real-time semantic cognition. We find that semantic cognition in the context of language comprehension is supported by trade-offs between widespread neural networks over the course of milliseconds. Incorporation of spatial and temporal characteristics, as well as behavioral measures, provide convergent evidence for the following progression: a hippocampal/anterior temporal phonological semantic retrieval network (peaking at ∼300 ms after the sentence final word); a frontotemporal thematic semantic network (∼400 ms); a hippocampal memory update network (∼500 ms); an inferior frontal semantic syntactic reappraisal network (∼600 ms); and nodes of the default mode network associated with conceptual coherence (∼750 ms). Additionally, in typical adults, mediatory relationships among these networks are significantly predictive of language comprehension ability. These findings provide a conceptual and methodological framework for the examination of speech and language disorders, with additional implications for the characterization of cognitive processes and clinical populations in other cognitive domains.SIGNIFICANCE STATEMENT The present study identifies a real-time neurobiological model of the meaning processes required during language comprehension (i.e., "semantic cognition"). Using a novel application of fused magnetic resonance imaging and electroencephalography in humans, we found that semantic cognition during language comprehension is supported by a rapid progression of widespread neural networks related to meaning, meaning integration, memory, reappraisal, and conceptual cohesion. Relationships among these systems were predictive of individuals' language comprehension efficiency. Our findings are the first to use fused neuroimaging analysis to elucidate language processes. In so doing, this study provides a new conceptual and methodological framework in which to characterize language processes and guide the treatment of speech and language deficits/disorders.

A ventral stream-prefrontal cortex processing cascade enables working memory gating dynamics

The representation of incoming information, goals and the flexible processing of these are required for cognitive control. Efficient mechanisms are needed to decide when it is important that novel information enters working memory (WM) and when these WM 'gates' have to be closed. Compared to neural foundations of maintaining information in WM, considerably less is known about what neural mechanisms underlie the representational dynamics during WM gating. Using different EEG analysis methods, we trace the path of mental representations along the human cortex during WM gate opening and closing. We show temporally nested representational dynamics during WM gate opening and closing depending on multiple independent neural activity profiles. These activity profiles are attributable to a ventral stream-prefrontal cortex processing cascade. The representational dynamics start in the ventral stream during WM gate opening and WM gate closing before prefrontal cortical regions are modulated. A regional specific activity profile is shown within the prefrontal cortex depending on whether WM gates are opened or closed, matching overarching concepts of prefrontal cortex functions. The study closes an essential conceptual gap detailing the neural dynamics underlying how mental representations drive the WM gate to open or close to enable WM functions such as updating and maintenance.

EEG microstate dynamics indicate a U-shaped path to propofol-induced loss of consciousness

Evidence suggests that the stream of consciousness is parsed into transient brain states manifesting themselves as discrete spatiotemporal patterns of global neuronal activity. Electroencephalographical (EEG) microstates are proposed as the neurophysiological correlates of these transiently stable brain states that last for fractions of seconds. To further understand the link between EEG microstate dynamics and consciousness, we continuously recorded high-density EEG in 23 surgical patients from their awake state to unconsciousness, induced by step-wise increasing concentrations of the intravenous anesthetic propofol. Besides the conventional parameters of microstate dynamics, we introduce a new implementation of a method to estimate the complexity of microstate sequences. The brain activity under the surgical anesthesia showed a decreased sequence complexity of the stereotypical microstates, which became sparser and longer-lasting. However, we observed an initial increase in microstates' temporal dynamics and complexity with increasing depth of sedation leading to a distinctive "U-shape" that may be linked to the paradoxical excitation induced by moderate levels of propofol. Our results support the idea that the brain is in a metastable state under normal conditions, balancing between order and chaos in order to flexibly switch from one state to another. The temporal dynamics of EEG microstates indicate changes of this critical balance between stability and transition that lead to altered states of consciousness.

Interpersonal Mindfulness Questionnaire: Scale Development and Validation

Objectives

Previously developed mindfulness measures focused on its intrapersonal dimensions and did not measure the interpersonal aspects of mindfulness. Furthermore, recently developed interpersonal mindfulness measures were either specific to a certain context (e.g., parenting, conjugal, teaching) or omitted/minimized the role of the body in the interpersonal dynamic. The proposed Interpersonal Mindfulness Questionnaire (IMQ) aims to operationalize the theoretical notion of embodied and embedded mindfulness by grounding it into four dimensions, each representing a set of skills that can be cultivated through training and practice: (1) Detachment from the Mind, (2) Body-Anchored Presence, (3) Attention to and Awareness of the Other Person, and (4) Mindful Responding.

Methods

The IMQ subscales were developed through consultations with a panel of eight graduate students and ten experts in the field. Three studies were conducted to evaluate the construct, internal consistency, reliability, convergent validity, and utility of the IMQ.

Results

Findings from the three studies supported the proposed four subscales of IMQ and suggested that these four subscales are independent and supported by convergent evidence. In addition, results suggested that IMQ subscales' scores are sensitive to meditation experience and are associated with better intrapersonal and interpersonal outcomes.

Conclusions

IMQ subscales are valid and are consistent with the proposed embodied and embedded conception of interpersonal mindfulness. IMQ subscales are associated with intrapersonal mindfulness, but not strongly enough to be conceived as the same phenomenon. Limitations, as well as theoretical and practical implications of IMQ subscales, are thoroughly discussed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12671-022-01855-1.

Large-scale EEG neural network changes in response to therapeutic TMS

BACKGROUND

Transcranial magnetic stimulation (TMS) is an effective therapy for patients with treatment-resistant depression. TMS likely induces functional connectivity changes in aberrant circuits implicated in depression. Electroencephalography (EEG) "microstates" are topographies hypothesized to represent large-scale resting networks. Canonical microstates have recently been proposed as markers for major depressive disorder (MDD), but it is not known if or how they change following TMS.

METHODS

Resting EEG was obtained from 49 MDD patients at baseline and following six weeks of daily TMS. Polarity-insensitive modified k-means clustering was used to segment EEGs into constituent microstates. Microstates were localized via sLORETA. Repeated-measures mixed models tested for within-subject differences over time and t-tests compared microstate features between TMS responder and non-responder groups.

RESULTS

Six microstates (MS-1 - MS-6) were identified from all available EEG data. Clinical response to TMS was associated with increases in features of MS-2, along with decreased metrics of MS-3. Nonresponders showed no significant changes in any microstate. Change in occurrence and coverage of both MS-2 (increased) and MS-3 (decreased) correlated with symptom change magnitude over the course of TMS treatment.

CONCLUSIONS

We identified EEG microstates associated with clinical improvement following a course of TMS therapy. Results suggest selective modulation of resting networks observable by EEG, which is inexpensive and easily acquired in the clinic setting.

Removal of movement-induced EEG artifacts: current state of the art and guidelines

Electroencephalography (EEG) is a non-invasive technique used to record cortical neurons' electrical activity using electrodes placed on the scalp. It has become a promising avenue for research beyond state-of-the-art EEG research that is conducted under static conditions. EEG signals are always contaminated by artifacts and other physiological signals. Artifact contamination increases with the intensity of movement. In the last decade (since 2010), researchers have started to implement EEG measurements in dynamic setups to increase the overall ecological validity of the studies. Many different methods are used to remove non-brain activity from the EEG signal, and there are no clear guidelines on which method should be used in dynamic setups and for specific movement intensities. Currently, the most common methods for removing artifacts in movement studies are methods based on independent component analysis (ICA). However, the choice of method for artifact removal depends on the type and intensity of movement, which affects the characteristics of the artifacts and the EEG parameters of interest. When dealing with EEG under non-static conditions, special care must be taken already in the designing period of an experiment. Software and hardware solutions must be combined to achieve sufficient removal of unwanted signals from EEG measurements. We have provided recommendations for the use of each method depending on the intensity of the movement and highlighted the advantages and disadvantages of the methods. However, due to the current gap in the literature, further development and evaluation of methods for artifact removal in EEG data during locomotion is needed.

Connectivity modulations induced by reach&grasp movements: a multidimensional approach

Reach&grasp requires highly coordinated activation of different brain areas. We investigated whether reach&grasp kinematics is associated to EEG-based networks changes. We enrolled 10 healthy subjects. We analyzed the reach&grasp kinematics of 15 reach&grasp movements performed with each upper limb. Simultaneously, we obtained a 64-channel EEG, synchronized with the reach&grasp movement time points. We elaborated EEG signals with EEGLAB 12 in order to obtain event related synchronization/desynchronization (ERS/ERD) and lagged linear coherence between Brodmann areas. Finally, we evaluated network topology via sLORETA software, measuring network local and global efficiency (clustering and path length) and the overall balance (small-worldness). We observed a widespread ERD in α and β bands during reach&grasp, especially in the centro-parietal regions of the hemisphere contralateral to the movement. Regarding functional connectivity, we observed an α lagged linear coherence reduction among Brodmann areas contralateral to the arm involved in the reach&grasp movement. Interestingly, left arm movement determined widespread changes of α lagged linear coherence, specifically among right occipital regions, insular cortex and somatosensory cortex, while the right arm movement exerted a restricted contralateral sensory-motor cortex modulation. Finally, no change between rest and movement was found for clustering, path length and small-worldness. Through a synchronized acquisition, we explored the cortical correlates of the reach&grasp movement. Despite EEG perturbations, suggesting that the non-dominant reach&grasp network has a complex architecture probably linked to the necessity of a higher visual control, the pivotal topological measures of network local and global efficiency remained unaffected.

Distinct modulation of mu and beta rhythm desynchronization during observation of embodied fake hand rotation

The rubber hand illusion (RHI) is a phenomenon whereby participants recognize a fake hand as their own. Studies have examined the effects of observing fake hand movements after the RHI on brain sensorimotor activity, although results remain controversial. To address these discrepancies, we investigated the effects of observation of fake hand rotation after the RHI on sensorimotor mu (μ: 8-13 Hz) and beta (β: 15-25 Hz) rhythm event-related desynchronization (ERD) using electroencephalography (EEG). Questionnaire results and proprioceptive drift revealed that the RHI occurred in participants when their invisible hand and fake visible hand were stroked synchronously but not during asynchronous stroking. Independent component (IC) clustering from EEG data during movement observation identified three IC clusters, including the right sensorimotor, left sensorimotor, and left occipital cluster. In the right sensorimotor cluster, we observed distinct modulation of μ and β ERD during fake hand rotation. Illusory ownership over the fake hand enhanced μ ERD but inversely attenuated β ERD. Further, the extent of μ ERD correlated with proprioceptive drift, but not with questionnaire ratings, whereas the converse results were obtained for β ERD. No ownership-dependent ERD modulation was detected in the left sensorimotor cluster. Alpha (α: 8-13 Hz) rhythm ERD of the left occipital cluster was smaller in the synchronous condition than in the asynchronous condition, but α ERD was not correlated with questionnaire rating or drift. These findings suggest that observing embodied fake hand rotation induces distinct cortical processing in sensorimotor brain areas.

Connectivity Measures Differentiate Cortical and Subcortical Sub-Acute Ischemic Stroke Patients

Brain lesions caused by cerebral ischemia lead to network disturbances in both hemispheres, causing a subsequent reorganization of functional connectivity both locally and remotely with respect to the injury. Quantitative electroencephalography (qEEG) methods have long been used for exploring brain electrical activity and functional connectivity modifications after stroke. However, results obtained so far are not univocal. Here, we used basic and advanced EEG methods to characterize how brain activity and functional connectivity change after stroke. Thirty-three unilateral post stroke patients in the sub-acute phase and ten neurologically intact age-matched right-handed subjects were enrolled. Patients were subdivided into two groups based on lesion location: cortico-subcortical (CS, n = 18) and subcortical (S, n = 15), respectively. Stroke patients were evaluated in the period ranging from 45 days since the acute event (T0) up to 3 months after stroke (T1) with both neurophysiological (resting state EEG) and clinical assessment (Barthel Index, BI) measures, while healthy subjects were evaluated once. Brain power at T0 was similar between the two groups of patients in all frequency bands considered (δ, θ, α, and β). However, evolution of θ-band power over time was different, with a normalization only in the CS group. Instead, average connectivity and specific network measures (Integration, Segregation, and Small-worldness) in the β-band at T0 were significantly different between the two groups. The connectivity and network measures at T0 also appear to have a predictive role in functional recovery (BI T1-T0), again group-dependent. The results obtained in this study showed that connectivity measures and correlations between EEG features and recovery depend on lesion location. These data, if confirmed in further studies, on the one hand could explain the heterogeneity of results so far observed in previous studies, on the other hand they could be used by researchers as biomarkers predicting spontaneous recovery, to select homogenous groups of patients for the inclusion in clinical trials.

Source Localization of EEG Brainwaves Activities via Mother Wavelets Families for SWT Decomposition

A Brain-Computer Interface (BCI) is a system used to communicate with an external world through the brain activity. The brain activity is measured by electroencephalography (EEG) signal and then processed by a BCI system. EEG source reconstruction could be a way to improve the accuracy of EEG classification in EEG based brain-computer interface (BCI). The source localization of the human brain activities can be an important resource for the recognition of the cognitive state, medical disorders, and a better understanding of the brain in general. In this study, we have compared 51 mother wavelets taken from 7 different wavelet families, which are applied to a Stationary Wavelet Transform (SWT) decomposition of an EEG signal. This process includes Haar, Symlets, Daubechies, Coiflets, Discrete Meyer, Biorthogonal, and reverse Biorthogonal wavelet families in extracting five different brainwave subbands for source localization. For this process, we used the Independent Component Analysis (ICA) for feature extraction followed by the Boundary Element Model (BEM) and the Equivalent Current Dipole (ECD) for the forward and inverse problem solutions. The evaluation results in investigating the optimal mother wavelet for source localization eventually identified the sym20 mother wavelet as the best choice followed by bior6.8 and coif5.

EEG changes reflecting pain: is alpha suppression better than gamma enhancement?

OBJECTIVE

Suppression of alpha and enhancement of gamma electroencephalographic (EEG) power have both been suggested as objective indicators of cortical pain processing. While gamma activity has been emphasized as the best potential marker, its spectral overlap with pain-related muscular responses is a potential drawback. Since muscle contractions are almost universal concomitants of physical pain, here we investigated alpha and gamma scalp-recorded activities during either tonic pain or voluntary facial grimaces mimicking those triggered by pain.

METHODS

High-density EEG (128 electrodes) was recorded while 14 healthy participants either underwent a cold pressor test (painful hand immersion in 10 °C water) or produced stereotyped facial/nuchal contractions (grimaces) mimicking those evoked by pain. The scalp distribution of spectral EEG changes was quantified via vector-transformation of maps and compared between the pain and grimacing conditions by calculating the cosine of the angle between the two corresponding topographies.

RESULTS

Painful stimuli significantly enhanced gamma power bilaterally in fronto-temporal regions and decreased alpha power in the contralateral central scalp. Sustained cervico-facial contractions (grimaces) gave also rise to significant gamma power increase in fronto-temporal regions but did not decrease central scalp alpha. While changes in alpha topography significantly differed between the pain and grimace situations, the scalp topography of gamma power was statistically indistinguishable from that occurring during grimaces.

CONCLUSION

Gamma power induced by painful stimuli or voluntary facial-cervical muscle contractions had overlapping topography. Pain-related alpha decrease in contralateral central scalp was less disturbed by muscle activity and may therefore prove more discriminant as an ancillary pain biomarker.

Stance leg and surface stability modulate cortical activity during human single leg stance

Mobile Electroencephalography (EEG) provides insights into cortical contributions to postural control. Although changes in theta (4–8 Hz) and alpha frequency power (8–12 Hz) were shown to reflect attentional and sensorimotor processing during balance tasks, information about the effect of stance leg on cortical processing related to postural control is lacking. Therefore, the aim was to examine patterns of cortical activity during single-leg stance with varying surface stability. EEG and force plate data from 21 healthy males (22.43 ± 2.23 years) was recorded during unipedal stance (left/right) on a stable and unstable surface. Using source-space analysis, power spectral density was analyzed in the theta, alpha-1 (8–10 Hz) and alpha-2 (10–12 Hz) frequency bands. Repeated measures ANOVA with the factors leg and surface stability revealed significant interaction effects in the left (p = 0.045, η_p² = 0.13) and right motor clusters (F = 16.156; p = 0.001, η_p² = 0.41). Furthermore, significant main effects for surface stability were observed for the fronto-central cluster (theta), left and right motor (alpha-1), as well as for the right parieto-occipital cluster (alpha-1/alpha-2). Leg dependent changes in alpha-2 power may indicate lateralized patterns of cortical processing in motor areas during single-leg stance. Future studies may therefore consider lateralized patterns of cortical activity for the interpretation of postural deficiencies in unilateral lower limb injuries.

Association between long-range temporal correlations in intrinsic EEG activity and subjective sense of identity

The long-range temporal correlation (LRTC) in resting-state intrinsic brain activity is known to be associated with temporal behavioral patterns, including decision making based on internal criteria such as self-knowledge. However, the association between the neuronal LRTC and the subjective sense of identity remains to be explored; in other words, whether our subjective sense of consistent self across time relates to the temporal consistency of neural activity. The present study examined the relationship between the LRTC of resting-state scalp electroencephalography (EEG) and a subjective sense of identity measured by the Erikson Psychosocial Stage Inventory (EPSI). Consistent with our prediction based on previous studies of neuronal-behavioral relationships, the frontocentral alpha LRTC correlated negatively with identity confusion. Moreover, from the descriptive analyses, centroparietal beta LRTC showed negative correlations with identity confusion, and frontal theta LRTC showed positive relationships with identity synthesis. These results suggest that more temporal consistency (reversely, less random noise) in intrinsic brain activity is associated with less confused and better-synthesized identity. Our data provide further evidence that the LRTC of intrinsic brain activity might serve as a noise suppression mechanism at the psychological level.

The open EEGLAB portal Interface: High-Performance computing with EEGLAB

EEGLAB signal processing environment is currently the leading open-source software for processing electroencephalographic (EEG) data. The Neuroscience Gateway (NSG, nsgportal.org) is a web and API-based portal allowing users to easily run a variety of neuroscience-related software on high-performance computing (HPC) resources in the U.S. XSEDE network. We have reported recently (Delorme et al., 2019) on the Open EEGLAB Portal expansion of the free NSG services to allow the neuroscience community to build and run MATLAB pipelines using the EEGLAB tool environment. We are now releasing an EEGLAB plug-in, nsgportal, that interfaces EEGLAB with NSG directly from within EEGLAB running on MATLAB on any personal lab computer. The plug-in features a flexible MATLAB graphical user interface (GUI) that allows users to easily submit, interact with, and manage NSG jobs, and to retrieve and examine their results. Command line nsgportal tools supporting these GUI functionalities allow EEGLAB users and plug-in tool developers to build largely automated functions and workflows that include optional NSG job submission and processing. Here we present details on nsgportal implementation and documentation, provide user tutorials on example applications, and show sample test results comparing computation times using HPC versus laptop processing.

Consistency of independent component analysis for FMRI

BACKGROUND

Independent component analysis (ICA) has been widely used for blind source separation in the field of medical imaging. However, despite of previous substantial efforts, the stability of ICA components remains a critical issue which has not been adequately addressed, despite numerous previous efforts. Most critical is the inconsistency of some of the extracted components when ICA is run with different model orders (MOs).

NEW METHOD

In this study, a novel method of determining the consistency of component analysis (CoCA) is proposed to evaluate the consistency of extracted components with different model orders. In the method, "consistent components" (CCs) are defined as those which can be extracted repeatably over a range of model orders.

RESULT

The efficacy of the method was evaluated with simulation data and fMRI datasets. With our method, the simulation result showed a clear difference of consistency between ground truths and noise.

COMPARISON WITH EXISTING METHODS

The information criteria were implemented to provide suggestions for the optimal model order, where some of the ICs were revealed inconsistent in our proposed method.

CONCLUSIONS

This method provided an objective protocol for choosing CCs of an ICA decomposition of a data matrix, independent of model order. This is especially useful with high model orders, where noise or other disturbances could possibly lead to an instability of the components.

Hybrid Human-Machine Interface for Gait Decoding Through Bayesian Fusion of EEG and EMG Classifiers

Despite the advances in the field of brain computer interfaces (BCI), the use of the sole electroencephalography (EEG) signal to control walking rehabilitation devices is currently not viable in clinical settings, due to its unreliability. Hybrid interfaces (hHMIs) represent a very recent solution to enhance the performance of single-signal approaches. These are classification approaches that combine multiple human-machine interfaces, normally including at least one BCI with other biosignals, such as the electromyography (EMG). However, their use for the decoding of gait activity is still limited. In this work, we propose and evaluate a hybrid human-machine interface (hHMI) to decode walking phases of both legs from the Bayesian fusion of EEG and EMG signals. The proposed hHMI significantly outperforms its single-signal counterparts, by providing high and stable performance even when the reliability of the muscular activity is compromised temporarily (e.g., fatigue) or permanently (e.g., weakness). Indeed, the hybrid approach shows a smooth degradation of classification performance after temporary EMG alteration, with more than 75% of accuracy at 30% of EMG amplitude, with respect to the EMG classifier whose performance decreases below 60% of accuracy. Moreover, the fusion of EEG and EMG information helps keeping a stable recognition rate of each gait phase of more than 80% independently on the permanent level of EMG degradation. From our study and findings from the literature, we suggest that the use of hybrid interfaces may be the key to enhance the usability of technologies restoring or assisting the locomotion on a wider population of patients in clinical applications and outside the laboratory environment.

Identifying key factors for improving ICA-based decomposition of EEG data in mobile and stationary experiments

Recent developments in EEG hardware and analyses approaches allow for recordings in both stationary and mobile settings. Irrespective of the experimental setting, EEG recordings are contaminated with noise that has to be removed before the data can be functionally interpreted. Independent component analysis (ICA) is a commonly used tool to remove artifacts such as eye movement, muscle activity, and external noise from the data and to analyze activity on the level of EEG effective brain sources. The effectiveness of filtering the data is one key preprocessing step to improve the decomposition that has been investigated previously. However, no study thus far compared the different requirements of mobile and stationary experiments regarding the preprocessing for ICA decomposition. We thus evaluated how movement in EEG experiments, the number of channels, and the high-pass filter cutoff during preprocessing influence the ICA decomposition. We found that for commonly used settings (stationary experiment, 64 channels, 0.5 Hz filter), the ICA results are acceptable. However, high-pass filters of up to 2 Hz cut-off frequency should be used in mobile experiments, and more channels require a higher filter to reach an optimal decomposition. Fewer brain ICs were found in mobile experiments, but cleaning the data with ICA has been proved to be important and functional even with low-density channel setups. Based on the results, we provide guidelines for different experimental settings that improve the ICA decomposition.

The influence of motor tasks and cut-off parameter selection on artifact subspace reconstruction in EEG recordings

Advances in EEG filtering algorithms enable analysis of EEG recorded during motor tasks. Although methods such as artifact subspace reconstruction (ASR) can remove transient artifacts automatically, there is virtually no knowledge about how the vigor of bodily movements affects ASRs performance and optimal cut-off parameter selection process. We compared the ratios of removed and reconstructed EEG recorded during a cognitive task, single-leg stance, and fast walking using ASR with 10 cut-off parameters versus visual inspection. Furthermore, we used the repeatability and dipolarity of independent components to assess their quality and an automatic classification tool to assess the number of brain-related independent components. The cut-off parameter equivalent to the ratio of EEG removed in manual cleaning was strictest for the walking task. The quality index of independent components, calculated using RELICA, reached a maximum plateau for cut-off parameters of 10 and higher across all tasks while dipolarity was largely unaffected. The number of independent components within each task remained constant, regardless of the cut-off parameter used. Surprisingly, ASR performed better in motor tasks compared with non-movement tasks. The quality index seemed to be more sensitive to changes induced by ASR compared to dipolarity. There was no benefit of using cut-off parameters less than 10. Graphical abstract The graphical abstract shows the three tasks performed during EEG recording, the two processing pipelines (manual and artifact subspace reconstruction), and the metrics the conclusion is based on.

Removal of eye-blinking artifacts by ICA in cross-modal long-term EEG recording

Independent Component Analysis (ICA) has became the most popular method to remove eye-blinking artifacts from electroencephalogram (EEG) recording. For long term EEG recording, ICA was commonly considered to costing a lot of computation time. Furthermore, with no ground truth, the discussion about the quality of ICA decomposition in a nonstationary environment was specious. In this study, we investigated the "signal" (P300 waveform) and the "noise" (averaged eye-blinking artifacts) on a cross-modal long-term EEG recording to evaluate the efficiency and effectiveness of different methods on ICA eye-blinking artifacts removal. As a result, it was found that, firstly, down sampling is an effective way to reduce the computation time in ICA. Appropriate down sampling ratio could speed up ICA computation 200 times and keep the decomposition performance stable, in which the computation time of ICA decomposition on a 2800 s EEG recording was less than 5 s. Secondly, dimension reduction by PCA was also a way to improve the efficiency and effectiveness of ICA. Finally, the comparison by cropping the dataset indicated that performing ICA on each run of the experiment separately would achieve a better result for eye-blinking artifacts removal than using all the EEG data input for ICA.

Modulations of Inter-Hemispherical Phase Coupling in Human Single Leg Stance

INTRODUCTION

Recent findings from neuroimaging studies provided initial insights into cortical contributions to postural control. These studies observed enhanced cortical activation and connectivity when task-difficulty and postural instability increased. However, little attention has been paid to the allocation of cortical networks appearing with a decreasing base of support from bipedal to single leg stance. Therefore, the aim of the present study was to investigate modulations of functional connectivity from bipedal to single leg stance.

EXPERIMENTAL PROCEDURES

Cortical activity during bipedal and single leg stance (left) was investigated in 15 male subjects using 128 channel mobile electroencephalography (EEG), while standing on a triaxial force plate. Power spectral density was calculated for theta (4-7 Hz), alpha-1 (8-10 Hz) and alpha-2 (10-12 Hz) frequency bands. Estimations of the phase lag index (PLI) were conducted as a measure of functional connectivity. Moreover, postural control was analyzed by the area of sway and sway velocity.

RESULTS

The results demonstrated a significantly increased area of sway and decreased alpha-2 power in single leg compared to bipedal stance. Furthermore, PLIs within the alpha-2 frequency band showed significantly decreased inter-hemispherical phase coupling in single leg stance, associated with connections involving the left motor region.

DISCUSSION

Altogether, the present findings may indicate modulations of cortical contributions in single leg compared to bipedal stance. The present data suggest that decreased inter-hemispherical functional connectivity, in conjunction with a global increase in cortical excitability, may indicate enhanced alertness and task-specific selective inhibition of motor networks in favor of postural control.

Characterization of multi-channel intraneural stimulation in transradial amputees

Although peripheral nerve stimulation using intraneural electrodes has been shown to be an effective and reliable solution to restore sensory feedback after hand loss, there have been no reports on the characterization of multi-channel stimulation. A deeper understanding of how the simultaneous stimulation of multiple electrode channels affects the evoked sensations should help in improving the definition of encoding strategies for bidirectional prostheses. We characterized the sensations evoked by simultaneous stimulation of median and ulnar nerves (multi-channel configuration) in four transradial amputees who had been implanted with four TIMEs (Transverse Intrafascicular Multichannel Electrodes). The results were compared with the characterization of single-channel stimulation. The sensations were characterized in terms of location, extent, type, and intensity. Combining two or more single-channel configurations caused a linear combination of the sensation locations and types perceived with such single-channel stimulations. Interestingly, this was also true when two active sites from the same nerve were stimulated. When stimulating in multi-channel configuration, the charge needed from each electrode channel to evoke a sensation was significantly lower than the one needed in single-channel configuration (sensory facilitation). This result was also supported by electroencephalography (EEG) recordings during nerve stimulation. Somatosensory potentials evoked by multi-channel stimulation confirmed that sensations in the amputated hand were perceived by the subjects and that a perceptual sensory facilitation occurred. Our results should help the future development of more efficient bidirectional prostheses by providing guidelines for the development of more complex stimulation approaches to effectively restore multiple sensations at the same time.

Pathway-specific modulatory effects of neuromuscular electrical stimulation during pedaling in chronic stroke survivors

Background

Neuromuscular electrical stimulation (NMES) is extensively used in stroke motor rehabilitation. How it promotes motor recovery remains only partially understood. NMES could change muscular properties, produce altered sensory inputs, and modulate fluctuations of cortical activities; but the potential contribution from cortico-muscular couplings during NMES synchronized with dynamic movement has rarely been discussed.

Method

We investigated cortico-muscular interactions during passive, active, and NMES rhythmic pedaling in healthy subjects and chronic stroke survivors. EEG (128 channels), EMG (4 unilateral lower limb muscles) and movement parameters were measured during 3 sessions of constant-speed pedaling. Sensory-level NMES (20 mA) was applied to the muscles, and cyclic stimulation patterns were synchronized with the EMG during pedaling cycles. Adaptive mixture independent component analysis was utilized to determine the movement-related electro-cortical sources and the source dipole clusters. A directed cortico-muscular coupling analysis was conducted between representative source clusters and the EMGs using generalized partial directed coherence (GPDC). The bidirectional GPDC was compared across muscles and pedaling sessions for post-stroke and healthy subjects.

Results

Directed cortico-muscular coupling of NMES cycling was more similar to that of active pedaling than to that of passive pedaling for the tested muscles. For healthy subjects, sensory-level NMES could modulate GPDC of both ascending and descending pathways. Whereas for stroke survivors, NMES could modulate GPDC of only the ascending pathways.

Conclusions

By clarifying how NMES influences neuromuscular control during pedaling in healthy and post-stroke subjects, our results indicate the potential limitation of sensory-level NMES in promoting sensorimotor recovery in chronic stroke survivors.

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.

Sensorimotor and Posterior Brain Activations During the Observation of Illusory Embodied Fake Hand Movement

In the rubber hand illusion (RHI), the subject recognizes a fake hand as his or her own. We recently found that the observation of embodied fake hand movement elicited mu-rhythm (8–13 Hz) desynchronization on electroencephalography (EEG), suggesting brain activation in the sensorimotor regions. However, it is known that mu-rhythm desynchronization during action observation is confounded with occipital alpha-rhythm desynchronization, which is modulated by attention. This study examined the independence of brain activities in the sensorimotor and occipital regions relating to the movement observation under the RHI. The invisible participant’s left and fake right hands were stroked simultaneously, which was interrupted by unexpected fake hand movements. A mirror-reversed image of the fake hand was shown on a monitor in front of the participant with a delay of 80, 280, or 480 ms. Illusion strength decreased as a function of the delay. EEG independent component analysis (ICA) and ICA clustering revealed six clusters with observation-induced desynchronization of 8–13 Hz frequency band. In the right sensorimotor cluster, mu-rhythm desynchronization was the greatest under the 80-ms delay, while alpha-rhythm desynchronization of the occipital clusters did not show delay-dependence. These results suggest that brain activation in the sensorimotor areas (i.e., mu-rhythm desynchronization) induced by embodied fake hand movement is independent of that in the occipital areas (alpha-rhythm desynchronization).

Trial-by-trial source-resolved EEG responses to gait task challenges predict subsequent step adaptation

A growing body of evidence indicates a pivotal role of cognition and in particular executive function in gait control and fall prevention. In a recent gait study using electroencephalographic (EEG) imaging, we provided direct proof for cortical top-down inhibitory control in step adaptation. A crucial part of motor inhibition is recognizing stimuli that signal the need to inhibit or adjust motor actions such as steps during walking. One of the EEG signatures of performance monitoring in response to events signaling the need to adjust motor responses, are error-related potential (error-ERP) features. To examine whether error-ERP features may index executive control during gait adaptation, we analyzed high-density (108-channel) EEG data from an auditory gait pacing study. Participants (N = 18) walking on a steadily moving treadmill were asked to step in time to an auditory cue tone sequence, and then to quickly adapt their step length and rate, to regain step-cue synchrony following occasional unexpected shifts in the pacing cue train to a faster or slower cue tempo. Decomposition of the continuous EEG data by independent component analysis revealed a negative deflection in the source-resolved event-related potential (ERP) time locked to 'late' cue tones marking a shift to a slower cue tempo. This vertex-negative ERP feature, localized primarily to posterior medial frontal cortex (pMFC) and peaking 250 ms after the onset of the tempo-shift cue, we here refer to as the step-cue delay negativity (SDN). SDN source, timing, and polarity resemble other error-related ERP features, e.g., the Error-Related Negativity (ERN) and Feedback-Related Negativity (FRN) in (seated) button press response tasks. In single trials, SDN amplitude varied with the magnitude of the cue latency deviation (the time interval between the expected and actual cue onsets). Regression analysis also identified linear coupling between SDN amplitude and the subsequent speed of gait tempo adaptation (as measured by the increase in length of the ensuing adaptation step). The SDN in this paradigm thus seems both to index the perceived need for and the subsequent magnitude of the immediate gait adjustment, consistent with performance-monitoring models. Future research might investigate relationships of these control processes to the impairment of gait adjustment in motor disorders and cognitive decline, for example to develop a biomarker for fall risk prediction in early-stage Parkinson's.

How are Muscle Synergies Affected by Electromyography Pre-Processing?

Muscle synergies have been used for decades to explain a variety of motor behaviors, both in humans and animals and, more recently, to steer rehabilitation strategies. However, many sources of variability such as factorization algorithms, criteria for dimensionality reduction and data pre-processing constitute a major obstacle to the successful comparison of the results obtained by different research groups. Starting from the canonical EMG processing we determined how variations in filter cut-off frequencies and normalization methods, commonly found in literature, affect synergy weights and inter-subject similarity (ISS) using experimental data related to a 15-muscles upper-limb reaching task. Synergy weights were not significantly altered by either normalization (maximum voluntary contraction - MVC - or maximum amplitude of the signal - SELF) or band-pass filter ([20-500 Hz] or [50-500] Hz). Normalization did, however, alter the amount of variance explained by a set of synergies, which is a criterion often used for model order selection. Comparing different low-pass (LP) filters (0.5 Hz, 4 Hz, 10 Hz, 20 Hz cut-offs) we showed that increasing the low pass filter cut-off had the effect of decreasing the variance accounted for by a set number of synergies and affected individual muscle contributions. Extreme smoothing (i.e., LP cut-off 0.5 Hz) enhanced the contrast between active and inactive muscles but had an unpredictable effect on the ISS. The results presented here constitute a further step towards a thoughtful EMG pre-processing for the extraction of muscle synergies.

Mining Within-Trial Oscillatory Brain Dynamics to Address the Variability of Optimized Spatial Filters

Data-driven spatial filtering algorithms optimize scores, such as the contrast between two conditions to extract oscillatory brain signal components. Most machine learning approaches for the filter estimation, however, disregard within-trial temporal dynamics and are extremely sensitive to changes in training data and involved hyperparameters. This leads to highly variable solutions and impedes the selection of a suitable candidate for, e.g., neurotechnological applications. Fostering component introspection, we propose to embrace this variability by condensing the functional signatures of a large set of oscillatory components into homogeneous clusters, each representing specific within-trial envelope dynamics. The proposed method is exemplified by and evaluated on a complex hand force task with a rich within-trial structure. Based on electroencephalography data of 18 healthy subjects, we found that the components' distinct temporal envelope dynamics are highly subject-specific. On average, we obtained seven clusters per subject, which were strictly confined regarding their underlying frequency bands. As the analysis method is not limited to a specific spatial filtering algorithm, it could be utilized for a wide range of neurotechnological applications, e.g., to select and monitor functionally relevant features for brain-computer interface protocols in stroke rehabilitation.

A visual working memory dataset collection with bootstrap Independent Component Analysis for comparison of electroencephalographic preprocessing pipelines

Here we present an electroencephalographic (EEG) collection of 71-channel datasets recorded from 14 subjects (7 males, 7 females, aged 20-40 years) while performing a visual working memory task with a T set of 150 Independent Component Analysis (ICA) decompositions by Extended Infomax using RELICA, each on a bootstrap resampling of the data. These data are linked to the paper "Applying dimension reduction to EEG data by Principal Component Analysis reduces the quality of its subsequent Independent Component decomposition" [1]. Independent components (ICs) are clustered within subject and thereby associated with a quality index (QIc) measure of their stability to data resampling. Sets of single ICA decompositions obtained after applying Principal Component Analysis (PCA) to the data to perform dimension reduction retaining (85%, 95%, 99%) of data variance are also included, as are the positions of the best fitting equivalent dipoles for ICs whose scalp projections are compatible with a compact brain source. These bootstrap ICs may be used as benchmarks for different data preprocessing pipelines and/or ICA algorithms, allowing investigation of the effects that noise or insufficient data have on the quality of ICA decompositions.

fMRI BOLD Correlates of EEG Independent Components: Spatial Correspondence With the Default Mode Network

Goal: We aimed to identify electroencephalographic (EEG) signal fluctuations within independent components (ICs) that correlate to spontaneous blood oxygenation level dependent (BOLD) activity in regions of the default mode network (DMN) during eyes-closed resting state. Methods: We analyzed simultaneously acquired EEG and functional magnetic resonance imaging (fMRI) eyes-closed resting state data in a convenience sample of 30 participants. IC analysis (ICA) was used to decompose the EEG time-series and common ICs were identified using data-driven IC clustering across subjects. The IC time courses were filtered into seven frequency bands, convolved with a hemeodynamic response function (HRF) and used to model spontaneous fMRI signal fluctuations across the brain. In parallel, group ICA analysis was used to decompose the fMRI signal into ICs from which the DMN was identified. Frequency and IC cluster associated hemeodynamic correlation maps obtained from the regression analysis were spatially correlated with the DMN. To investigate the reliability of our findings, the analyses were repeated with data collected from the same subjects 1 year later. Results: Our results indicate a relationship between power fluctuations in the delta, theta, beta and gamma frequency range and the DMN in different EEG ICs in our sample as shown by small to moderate spatial correlations at the first measurement (0.234 < |r| < 0.346, p < 0.0001). Furthermore, activity within an EEG component commonly identified as eye movements correlates with BOLD activity within regions of the DMN. In addition, we demonstrate that correlations between EEG ICs and the BOLD signal during rest are in part stable across time. Discussion: We show that ICA source separated EEG signals can be used to investigate electrophysiological correlates of the DMN. The relationship between the eye movement component and the DMN points to a behavioral association between DMN activity and the level of eye movement or the presence of neuronal activity in this component. Previous findings of an association between frontal midline theta activity and the DMN were replicated.

Automatic Independent Component Scalp Map Analysis of Electroencephalogram During Motor Preparation

This work presents a method for automatic independent component (IC) scalp map analysis of electroencephalogram during motor preparation in visuomotor tasks. The strength of this approach is the analysis of the IC scalp maps based on the apriori given mask. This uses an image processing approach, comparable to visual classification used by experts, to automate the selection of relevant ICs in visuomotor tasks. Thirty iterations of the Infomax ICA algorithm were used to test the reliability of the ICs. ICs above 95% quality index were used for IC scalp topography image analysis. Here, we used a linkage-clustering algorithm for IC clustering and gap statistic to estimate the number of clusters. After classifying the components with our approach, the labels were compared to those from well-known MARA ("Multiple Artifact Rejection Algorithm") - an open-source EEGLAB plug-in. It was found that 334 of the 568 labels were in-agreement. MARA labeled 81 out of the 177 source-related components, and 238 out of the 319 non-source-related components, as artifacts. Here, the strength of our approach lies in using an image-processing algorithm to identify the task-specific ICs whereas MARA focuses on the automatic classification of the artifactual ICs by combining stereotyped artifact-specific spatial and temporal features that depend on the electrode montage. After "artefactual" ICs are removed, task-specific ICs still remains to be identified from the remaining "good" ICs where our scalp topography image analysis approach can be applied. Our IC scalp topography image analysis is focused on task-specific IC selection based on an apriori mask, which is not limited to specific EEG features and/or electrode configurations for high-density EEG.

Data-driven body-machine interface for the accurate control of drones

The teleoperation of nonhumanoid robots is often a demanding task, as most current control interfaces rely on mappings between the operator’s and the robot’s actions, which are determined by the design and characteristics of the interface, and may therefore be challenging to master. Here, we describe a structured methodology to identify common patterns in spontaneous interaction behaviors, to implement embodied user interfaces, and to select the appropriate sensor type and positioning. Using this method, we developed an intuitive, gesture-based control interface for real and simulated drones, which outperformed a standard joystick in terms of learning time and steering abilities. Implementing this procedure to identify body-machine patterns for specific applications could support the development of more intuitive and effective interfaces.

The accurate teleoperation of robotic devices requires simple, yet intuitive and reliable control interfaces. However, current human–machine interfaces (HMIs) often fail to fulfill these characteristics, leading to systems requiring an intensive practice to reach a sufficient operation expertise. Here, we present a systematic methodology to identify the spontaneous gesture-based interaction strategies of naive individuals with a distant device, and to exploit this information to develop a data-driven body–machine interface (BoMI) to efficiently control this device. We applied this approach to the specific case of drone steering and derived a simple control method relying on upper-body motion. The identified BoMI allowed participants with no prior experience to rapidly master the control of both simulated and real drones, outperforming joystick users, and comparing with the control ability reached by participants using the bird-like flight simulator Birdly.

Bilateral cortical representation of tactile roughness

Roughness is the most important feature for texture discrimination. Here we investigate how the bilateral cortical representation of touch is modulated by tactile roughness by analyzing the neural responses elicited by stimuli with various coarseness levels ranging from fine to medium. A prolonged stimulation was delivered to 10 healthy subjects by passively sliding tactile stimuli under the fingertip while recording the EEG to study the modulation of Somatosensory Evoked Potentials (SEPs) as well as activity in the theta and alpha bands. Elicited long-latency SEPs, namely bilateral P100-N140 and frontal P240 were consistent across stimuli. On the contrary, the temporal lag N140 - P240 was nonlinearly modulated both in contralateral and ipsilateral sides, in agreement with literature. Using a time-frequency analysis approach, we identified a theta band power increase in the [0 0.5]s interval and a partially overlapped power decrease in the alpha band which lasted throughout the stimulation. The estimated time these two phenomena were overlapped was comparable across stimuli, whereas a linear decrease in alpha band amplitude was reported when increasing the stimulus roughness in both contralateral and ipsilateral sides. This study showed that the selected tactile stimuli generated physiological bilateral responses that were modulated in a diversified way according to the stimulus roughness and side. Specifically, we identified sensory processing features (i.e., theta and alpha time overlap) invariant to the stimulus roughness (i.e., associated to a basic cortical mechanism of touch) and roughness-dependent cortical outputs comparable in the contralateral and ipsilateral sides that confirm a bilateral processing of tactile information.

Subanaesthetic ketamine and altered states of consciousness in humans

BACKGROUND

Despite its designation as a 'dissociative anaesthetic,' the dissociative and psychoactive effects of ketamine remain incompletely understood. The goal of this study was to characterise the subjective experiences and accompanying EEG changes with subanaesthetic doses of ketamine.

METHODS

High-density EEG was recorded in 15 human volunteers before, during, and after subanaesthetic ketamine infusion (0.5 mg kg^-1 over 40 min), with self-reported measures of altered states of consciousness obtained after ketamine exposure. Sensor- and source-level EEG changes were analysed with a focus on spectral power and regional changes.

RESULTS

Ketamine-induced altered states were characterised predominantly by dissociative experiences such as disembodiment and ego transcendence; sensory disturbances were also common. Ketamine broadly decreased low-frequency power, with mean reductions largest at alpha (8-12 Hz) in parietal (-0.94 dB, P<0.001) and occipital (-1.8 dB, P<0.001) channel clusters. Significant decreases in alpha were identified in the precuneus and temporal-parietal junction.

CONCLUSIONS

Ketamine induces altered states of consciousness during periods of reduced alpha power in the precuneus and temporal-parietal junction. Modulation of these temporal-parietal loci are candidate mechanisms of the psychoactive effects of ketamine, given that this region is involved in multisensory integration, body representation, and consciousness.