Realness of face images can be decoded from non-linear modulation of EEG responses.

Artificially created human faces play an increasingly important role in our digital world. However, the so-called uncanny valley effect may cause people to perceive highly, yet not perfectly human-like faces as eerie, bringing challenges to the interaction with virtual agents. At the same time, the neurocognitive underpinnings of the uncanny valley effect remain elusive. Here, we utilized an electroencephalography (EEG) dataset of steady-state visual evoked potentials (SSVEP) in which participants were presented with human face images of different stylization levels ranging from simplistic cartoons to actual photographs. Assessing neuronal responses both in frequency and time domain, we found a non-linear relationship between SSVEP amplitudes and stylization level, that is, the most stylized cartoon images and the real photographs evoked stronger responses than images with medium stylization. Moreover, realness of even highly similar stylization levels could be decoded from the EEG data with task-related component analysis (TRCA). Importantly, we also account for confounding factors, such as the size of the stimulus face's eyes, which previously have not been adequately addressed. Together, this study provides a basis for future research and neuronal benchmarking of real-time detection of face realness regarding three aspects: SSVEP-based neural markers, efficient classification methods, and low-level stimulus confounders.

Decreased long-range temporal correlations in the resting-state functional magnetic resonance imaging blood-oxygen-level-dependent signal reflect motor sequence learning up to 2 weeks following training.

Decreased long-range temporal correlations (LRTC) in brain signals can be used to measure cognitive effort during task execution. Here, we examined how learning a motor sequence affects long-range temporal memory within resting-state functional magnetic resonance imaging signal. Using the Hurst exponent (HE), we estimated voxel-wise LRTC and assessed changes over 5 consecutive days of training, followed by a retention scan 12 days later. The experimental group learned a complex visuomotor sequence while a complementary control group performed tightly matched movements. An interaction analysis revealed that HE decreases were specific to the complex sequence and occurred in well-known motor sequence learning associated regions including left supplementary motor area, left premotor cortex, left M1, left pars opercularis, bilateral thalamus, and right striatum. Five regions exhibited moderate to strong negative correlations with overall behavioral performance improvements. Following learning, HE values returned to pretraining levels in some regions, whereas in others, they remained decreased even 2 weeks after training. Our study presents new evidence of HE's possible relevance for functional plasticity during the resting-state and suggests that a cortical subset of sequence-specific regions may continue to represent a functional signature of learning reflected in decreased long-range temporal dependence after a period of inactivity.

MEANSP: How Many Channels are Needed to Predict the Performance of a SMR-Based BCI?

Predicting whether a particular individual would reach an adequate control of a Brain-Computer Interface (BCI) has many practical advantages. On the one hand, participants with low predicted performance could be trained with specifically designed sessions and avoid frustrating experiments; on the other hand, planning time and resources would be more efficient; and finally, the variables related to an accurate prediction could be manipulated to improve the prospective BCI performance. To this end, several predictors have been proposed in the literature, most of them based on the power estimation of EEG signals at the specific frequency bands. Many of these studies evaluate their predictors in relatively small datasets and/or using a relatively high number of channels. In this manuscript, we propose a novel predictor called [Formula: see text] to predict the performance of participants using BCIs that are based on the modulation of sensorimotor rhythms. This novel predictor has been positively evaluated using only 2, 3, 4 or 5 channels. [Formula: see text] has shown to perform as well as or better than other state-of-the-art predictors. The best sets of different number of channels are also provided, which have been tested in two different settings to prove their robustness. The proposed predictor has been successfully evaluated using two large-scale datasets containing 150 and 80 participants, respectively. We also discuss predictor thresholds for users to expect good performance in feedback experiments and show the advantages in comparison to a competing algorithm.

Event-related modulation of alpha rhythm explains the auditory P300-evoked response in EEG.

Evoked responses and oscillations represent two major electrophysiological phenomena in the human brain yet the link between them remains rather obscure. Here we show how most frequently studied EEG signals: the P300-evoked response and alpha oscillations (8-12 Hz) can be linked with the baseline-shift mechanism. This mechanism states that oscillations generate evoked responses if oscillations have a non-zero mean and their amplitude is modulated by the stimulus. Therefore, the following predictions should hold: (1) the temporal evolution of P300 and alpha amplitude is similar, (2) spatial localisations of the P300 and alpha amplitude modulation overlap, (3) oscillations are non-zero mean, (4) P300 and alpha amplitude correlate with cognitive scores in a similar fashion. To validate these predictions, we analysed the data set of elderly participants (N=2230, 60-82 years old), using (a) resting-state EEG recordings to quantify the mean of oscillations, (b) the event-related data, to extract parameters of P300 and alpha rhythm amplitude envelope. We showed that P300 is indeed linked to alpha rhythm, according to all four predictions. Our results provide an unifying view on the interdependency of evoked responses and neuronal oscillations and suggest that P300, at least partly, is generated by the modulation of alpha oscillations.

Cardiac activity impacts cortical motor excitability.

Human cognition and action can be influenced by internal bodily processes such as heartbeats. For instance, somatosensory perception is impaired both during the systolic phase of the cardiac cycle and when heartbeats evoke stronger cortical responses. Here, we test whether these cardiac effects originate from overall changes in cortical excitability. Cortical and corticospinal excitability were assessed using electroencephalographic and electromyographic responses to transcranial magnetic stimulation while concurrently monitoring cardiac activity with electrocardiography. Cortical and corticospinal excitability were found to be highest during systole and following stronger neural responses to heartbeats. Furthermore, in a motor task, hand-muscle activity and the associated desynchronization of sensorimotor oscillations were stronger during systole. These results suggest that systolic cardiac signals have a facilitatory effect on motor excitability-in contrast to sensory attenuation that was previously reported for somatosensory perception. Thus, it is possible that distinct time windows exist across the cardiac cycle, optimizing either perception or action.

Visual short-term memory-related EEG components in a virtual reality setup.

Virtual reality (VR) offers a powerful tool for investigating cognitive processes, as it allows researchers to gauge behaviors and mental states in complex, yet highly controlled, scenarios. The use of VR head-mounted displays in combination with physiological measures such as EEG presents new challenges and raises the question whether established findings also generalize to a VR setup. Here, we used a VR headset to assess the spatial constraints underlying two well-established EEG correlates of visual short-term memory: the amplitude of the contralateral delay activity (CDA) and the lateralization of induced alpha power during memory retention. We tested observers' visual memory in a change detection task with bilateral stimulus arrays of either two or four items while varying the horizontal eccentricity of the memory arrays (4, 9, or 14 degrees of visual angle). The CDA amplitude differed between high and low memory load at the two smaller eccentricities, but not at the largest eccentricity. Neither memory load nor eccentricity significantly influenced the observed alpha lateralization. We further fitted time-resolved spatial filters to decode memory load from the event-related potential as well as from its time-frequency decomposition. Classification performance during the retention interval was above-chance level for both approaches and did not vary significantly across eccentricities. We conclude that commercial VR hardware can be utilized to study the CDA and lateralized alpha power, and we provide caveats for future studies targeting these EEG markers of visual memory in a VR setup.

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.

(Non)-experiencing the intention to move: On the comparisons between the Readiness Potential onset and Libet's W-time.

A seminal study by Libet et al. (1983) provided a popular approach to compare the introspective timing of movement execution (the M-time) and the intention to move (the W-time) with respect to the onset of the readiness potential (RP). The difference between the W-time and the RP onsets contributed significantly to the current free-will discussion, insofar as it has been repeatedly shown that the RP onset unequivocally precedes the W-time. However, the interpretations of Libet's paradigm continuously attract criticism, questioning the use of both the W-time and the RP onset as indicators of motor intention. In the current study, we further probe whether the W-time is rather an intention-unrelated product of the participant's inference than an unambiguous temporal marker of the intention to move. Using behavioral reports and concurrent multichannel EEG, we investigated the relationship between the W-time and M-time introspective reports in two groups of participants who started an experiment with a series of different reports. Congruently with previous studies, we have shown that the W-time is affected by the experimental procedures: participants who had prior experience reporting the M-time provided significantly earlier W-time. However, contrary to previous papers, we revealed that even naive participants do introspectively differentiate the W-time and the M-time, which suggests that the W-time might actually reflect the intention to move, at least to some extent. We, therefore, suggest that training-based modulation of the W-time values may explain this finding. Moreover, we further confirm the absence of a direct link between the RP onset and the W-time by showing no covariation between them in both experimental groups. In turn, our findings question the overall interpretation of the comparison between these two time points. Overall, our study further emphasizes the ambiguity of Libet's paradigm, and suggests that the relatedness of both the RP and the W-time to the movement initiation processes should not be assumed a priori.

Prior Movement of One Arm Facilitates Motor Adaptation in the Other.

Many movements in daily life are embedded in motion sequences that involve more than one limb, demanding the motor system to monitor and control different body parts in quick succession. During such movements, systematic changes in the environment or the body might require motor adaptation of specific segments. However, previous motor adaptation research has focused primarily on motion sequences produced by a single limb, or on simultaneous movements of several limbs. For example, adaptation to opposing force fields is possible in unimanual reaching tasks when the direction of a prior or subsequent movement is predictive of force field direction. It is unclear, however, whether multilimb sequences can support motor adaptation processes in a similar way. In the present study (38 females, 38 males), we investigated whether reaches can be adapted to different force fields in a bimanual motor sequence when the information about the perturbation is associated with the prior movement direction of the other arm. In addition, we examined whether prior perceptual (visual or proprioceptive) feedback of the opposite arm contributes to force field-specific motor adaptation. Our key finding is that only active participation in the bimanual sequential task supports pronounced adaptation. This result suggests that active segments in bimanual motion sequences are linked across limbs. If there is a consistent association between movement kinematics of the linked and goal movement, the learning process of the goal movement can be facilitated. More generally, if motion sequences are repeated often, prior segments can evoke specific adjustments of subsequent movements.SIGNIFICANCE STATEMENT Movements in a limb's motion sequence can be adjusted based on linked movements. A prerequisite is that kinematics of the linked movements correctly predict which adjustments are needed. We show that use of kinematic information to improve performance is even possible when a prior linked movement is performed with a different limb. For example, a skilled juggler might have learned how to correctly adjust his catching movement of the left hand when the right hand performed a throwing action in a specific way. Linkage is possibly a key mechanism of the human motor system for learning complex bimanual skills. Our study emphasizes that learning of specific movements should not be studied in isolation but within their motor sequence context.

Anterior cingulate and medial prefrontal cortex oscillations underlie learning alterations in trait anxiety in humans.

Anxiety has been linked to altered belief formation and uncertainty estimation, impacting learning. Identifying the neural processes underlying these changes is important for understanding brain pathology. Here, we show that oscillatory activity in the medial prefrontal, anterior cingulate and orbitofrontal cortex (mPFC, ACC, OFC) explains anxiety-related learning alterations. In a magnetoencephalography experiment, two groups of human participants pre-screened with high and low trait anxiety (HTA, LTA: 39) performed a probabilistic reward-based learning task. HTA undermined learning through an overestimation of volatility, leading to faster belief updating, more stochastic decisions and pronounced lose-shift tendencies. On a neural level, we observed increased gamma activity in the ACC, dmPFC, and OFC during encoding of precision-weighted prediction errors in HTA, accompanied by suppressed ACC alpha/beta activity. Our findings support the association between altered learning and belief updating in anxiety and changes in gamma and alpha/beta activity in the ACC, dmPFC, and OFC.

Alterations in rhythmic and non-rhythmic resting-state EEG activity and their link to cognition in older age.

While many structural and biochemical changes in the brain have previously been associated with older age, findings concerning functional properties of neuronal networks, as reflected in their electrophysiological signatures, remain rather controversial. These discrepancies might arise due to several reasons, including diverse factors determining general spectral slowing in the alpha frequency range as well as amplitude mixing between the rhythmic and non-rhythmic parameters. We used a large dataset (N = 1703, mean age 70) to comprehensively investigate age-related alterations in multiple EEG biomarkers taking into account rhythmic and non-rhythmic activity and their individual contributions to cognitive performance. While we found strong evidence for an individual alpha peak frequency (IAF) decline in older age, we did not observe a significant relationship between theta power and age while controlling for IAF. Not only did IAF decline with age, but it was also positively associated with interference resolution in a working memory task primarily in the right and left temporal lobes suggesting its functional role in information sampling. Critically, we did not detect a significant relationship between alpha power and age when controlling for the 1/f spectral slope, while the latter one showed age-related alterations. These findings thus suggest that the entanglement of IAF slowing and power in the theta frequency range, as well as 1/f slope and alpha power measures, might explain inconsistencies reported previously in the literature. Finally, despite the absence of age-related alterations, alpha power was negatively associated with the speed of processing in the right frontal lobe while 1/f slope showed no consistent relationship to cognitive performance. Our results thus demonstrate that multiple electrophysiological features, as well as their interplay, should be considered for the comprehensive assessment of association between age, neuronal activity, and cognitive performance.

Movement decoding using spatio-spectral features of cortical and subcortical local field potentials.

The first commercially sensing enabled deep brain stimulation (DBS) devices for the treatment of movement disorders have recently become available. In the future, such devices could leverage machine learning based brain signal decoding strategies to individualize and adapt therapy in real-time. As multi-channel recordings become available, spatial information may provide an additional advantage for informing machine learning models. To investigate this concept, we compared decoding performances from single channels vs. spatial filtering techniques using intracerebral multitarget electrophysiology in Parkinson's disease patients undergoing DBS implantation. We investigated the feasibility of spatial filtering in invasive neurophysiology and the putative utility of combined cortical ECoG and subthalamic local field potential signals for decoding grip-force, a well-defined and continuous motor readout. We found that adding spatial information to the model can improve decoding (6% gain in decoding), but the spatial patterns and additional benefit was highly individual. Beyond decoding performance results, spatial filters and patterns can be used to obtain meaningful neurophysiological information about the brain networks involved in target behavior. Our results highlight the importance of individualized approaches for brain signal decoding, for which multielectrode recordings and spatial filtering can improve precision medicine approaches for clinical brain computer interfaces.

Identification of spatial patterns with maximum association between power of resting state neural oscillations and trait anxiety.

Anxiety affects approximately 5-10% of the adult population worldwide, placing a large burden on the health systems. Despite its omnipresence and impact on mental and physical health, most of the individuals affected by anxiety do not receive appropriate treatment. Current research in the field of psychiatry emphasizes the need to identify and validate biological markers relevant to this condition. Neurophysiological preclinical studies are a prominent approach to determine brain rhythms that can be reliable markers of key features of anxiety. However, while neuroimaging research consistently implicated prefrontal cortex and subcortical structures, such as amygdala and hippocampus, in anxiety, there is still a lack of consensus on the underlying neurophysiological processes contributing to this condition. Methods allowing non-invasive recording and assessment of cortical processing may provide an opportunity to help identify anxiety signatures that could be used as intervention targets. In this study, we apply Source-Power Comodulation (SPoC) to electroencephalography (EEG) recordings in a sample of participants with different levels of trait anxiety. SPoC was developed to find spatial filters and patterns whose power comodulates with an external variable in individual participants. The obtained patterns can be interpreted neurophysiologically. Here, we extend the use of SPoC to a multi-subject setting and test its validity using simulated data with a realistic head model. Next, we apply our SPoC framework to resting state EEG of 43 human participants for whom trait anxiety scores were available. SPoC inter-subject analysis of narrow frequency band data reveals neurophysiologically meaningful spatial patterns in the theta band (4-7 Hz) that are negatively correlated with anxiety. The outcome is specific to the theta band and not observed in the alpha (8-12 Hz) or beta (13-30 Hz) frequency range. The theta-band spatial pattern is primarily localised to the superior frontal gyrus. We discuss the relevance of our spatial pattern results for the search of biomarkers for anxiety and their application in neurofeedback studies.

Accuracy and precision of navigated transcranial magnetic stimulation.

Objective.Transcranial magnetic stimulation (TMS) induces an electric field (E-field) in the cortex. To facilitate stimulation targeting, image-guided neuronavigation systems have been introduced. Such systems track the placement of the coil with respect to the head and visualize the estimated cortical stimulation location on an anatomical brain image in real time. The accuracy and precision of the neuronavigation is affected by multiple factors. Our aim was to analyze how different factors in TMS neuronavigation affect the accuracy and precision of the coil-head coregistration and the estimated E-field.Approach.By performing simulations, we estimated navigation errors due to distortions in magnetic resonance images (MRIs), head-to-MRI registration (landmark- and surface-based registrations), localization and movement of the head tracker, and localization of the coil tracker. We analyzed the effect of these errors on coil and head coregistration and on the induced E-field as determined with simplistic and realistic head models.Main results.Average total coregistration accuracies were in the range of 2.2-3.6 mm and 1°; precision values were about half of the accuracy values. The coregistration errors were mainly due to head-to-MRI registration with average accuracies 1.5-1.9 mm/0.2-0.4° and precisions 0.5-0.8 mm/0.1-0.2° better with surface-based registration. The other major source of error was the movement of the head tracker with average accuracy of 1.5 mm and precision of 1.1 mm. When assessed within an E-field method, the average accuracies of the peak E-field location, orientation, and magnitude ranged between 1.5 and 5.0 mm, 0.9 and 4.8°, and 4.4 and 8.5% across the E-field models studied. The largest errors were obtained with the landmark-based registration. When computing another accuracy measure with the most realistic E-field model as a reference, the accuracies tended to improve from about 10 mm/15°/25% to about 2 mm/2°/5% when increasing realism of the E-field model.Significance.The results of this comprehensive analysis help TMS operators to recognize the main sources of error in TMS navigation and that the coregistration errors and their effect in the E-field estimation depend on the methods applied. To ensure reliable TMS navigation, we recommend surface-based head-to-MRI registration and realistic models for E-field computations.

Attenuation of the Heartbeat-Evoked Potential in Patients With Atrial Fibrillation.

BACKGROUND: The heartbeat-evoked potential (HEP) is a brain response to each heartbeat, which is thought to reflect cardiac signaling to central autonomic areas and suggested to be a marker of internal body awareness (eg, interoception). OBJECTIVES: Because cardiac communication with central autonomic circuits has been shown to be impaired in patients with atrial fibrillation (AF), we hypothesized that HEPs are attenuated in these patients. METHODS: By simultaneous electroencephalography and electrocardiography recordings, HEP was investigated in 56 individuals with persistent AF and 56 control subjects matched for age, sex, and body mass index. RESULTS: HEP in control subjects was characterized by right frontotemporal negativity peaking around 300 to 550 ms after the R-peak, consistent with previous studies. In comparison with control subjects, HEP amplitudes were attenuated, and HEP amplitude differences remained significant when matching the samples for heart frequency, stroke volume (assessed by echocardiography), systolic blood pressure, and the amplitude of the T-wave. Effect sizes for the group differences were medium to large (Cohen's d between 0.6 and 0.9). EEG source analysis on HEP amplitude differences pointed to a neural representation within the right insular cortex, an area known as a hub for central autonomic control. CONCLUSIONS: The heartbeat-evoked potential is reduced in AF, particularly in the right insula. We speculate that the attenuated HEP in AF may be a marker of impaired heart-brain interactions. Attenuated interoception might furthermore underlie the frequent occurrence of silent AF.

Non-zero mean alpha oscillations revealed with computational model and empirical data.

Ongoing oscillations and evoked responses are two main types of neuronal activity obtained with diverse electrophysiological recordings (EEG/MEG/iEEG/LFP). Although typically studied separately, they might in fact be closely related. One possibility to unite them is to demonstrate that neuronal oscillations have non-zero mean which predicts that stimulus- or task-triggered amplitude modulation of oscillations can contribute to the generation of evoked responses. We validated this mechanism using computational modelling and analysis of a large EEG data set. With a biophysical model, we indeed demonstrated that intracellular currents in the neuron are asymmetric and, consequently, the mean of alpha oscillations is non-zero. To understand the effect that neuronal currents exert on oscillatory mean, we varied several biophysical and morphological properties of neurons in the network, such as voltage-gated channel densities, length of dendrites, and intensity of incoming stimuli. For a very large range of model parameters, we observed evidence for non-zero mean of oscillations. Complimentary, we analysed empirical rest EEG recordings of 90 participants (50 young, 40 elderly) and, with spatio-spectral decomposition, detected at least one spatially-filtred oscillatory component of non-zero mean alpha oscillations in 93% of participants. In order to explain a complex relationship between the dynamics of amplitude-envelope and corresponding baseline shifts, we performed additional simulations with simple oscillators coupled with different time delays. We demonstrated that the extent of spatial synchronisation may obscure macroscopic estimation of alpha rhythm modulation while leaving baseline shifts unchanged. Overall, our results predict that amplitude modulation of neural oscillations should at least partially explain the generation of evoked responses. Therefore, inference about changes in evoked responses with respect to cognitive conditions, age or neuropathologies should be constructed while taking into account oscillatory neuronal dynamics.

Dopaminergic Modulation of Local Non-oscillatory Activity and Global-Network Properties in Parkinson's Disease: An EEG Study.

Dopaminergic medication for Parkinson's disease (PD) modulates neuronal oscillations and functional connectivity (FC) across the basal ganglia-thalamic-cortical circuit. However, the non-oscillatory component of the neuronal activity, potentially indicating a state of excitation/inhibition balance, has not yet been investigated and previous studies have shown inconsistent changes of cortico-cortical connectivity as a response to dopaminergic medication. To further elucidate changes of regional non-oscillatory component of the neuronal power spectra, FC, and to determine which aspects of network organization obtained with graph theory respond to dopaminergic medication, we analyzed a resting-state electroencephalography (EEG) dataset including 15 PD patients during OFF and ON medication conditions. We found that the spectral slope, typically used to quantify the broadband non-oscillatory component of power spectra, steepened particularly in the left central region in the ON compared to OFF condition. In addition, using lagged coherence as a FC measure, we found that the FC in the beta frequency range between centro-parietal and frontal regions was enhanced in the ON compared to the OFF condition. After applying graph theory analysis, we observed that at the lower level of topology the node degree was increased, particularly in the centro-parietal area. Yet, results showed no significant difference in global topological organization between the two conditions: either in global efficiency or clustering coefficient for measuring global and local integration, respectively. Interestingly, we found a close association between local/global spectral slope and functional network global efficiency in the OFF condition, suggesting a crucial role of local non-oscillatory dynamics in forming the functional global integration which characterizes PD. These results provide further evidence and a more complete picture for the engagement of multiple cortical regions at various levels in response to dopaminergic medication in PD.

Separating Neural Oscillations from Aperiodic 1/f Activity: Challenges and Recommendations.

Electrophysiological power spectra typically consist of two components: An aperiodic part usually following an 1/f power law [Formula: see text] and periodic components appearing as spectral peaks. While the investigation of the periodic parts, commonly referred to as neural oscillations, has received considerable attention, the study of the aperiodic part has only recently gained more interest. The periodic part is usually quantified by center frequencies, powers, and bandwidths, while the aperiodic part is parameterized by the y-intercept and the 1/f exponent [Formula: see text]. For investigation of either part, however, it is essential to separate the two components. In this article, we scrutinize two frequently used methods, FOOOF (Fitting Oscillations & One-Over-F) and IRASA (Irregular Resampling Auto-Spectral Analysis), that are commonly used to separate the periodic from the aperiodic component. We evaluate these methods using diverse spectra obtained with electroencephalography (EEG), magnetoencephalography (MEG), and local field potential (LFP) recordings relating to three independent research datasets. Each method and each dataset poses distinct challenges for the extraction of both spectral parts. The specific spectral features hindering the periodic and aperiodic separation are highlighted by simulations of power spectra emphasizing these features. Through comparison with the simulation parameters defined a priori, the parameterization error of each method is quantified. Based on the real and simulated power spectra, we evaluate the advantages of both methods, discuss common challenges, note which spectral features impede the separation, assess the computational costs, and propose recommendations on how to use them.

Harmoni: A method for eliminating spurious interactions due to the harmonic components in neuronal data.

Cross-frequency synchronization (CFS) has been proposed as a mechanism for integrating spatially and spectrally distributed information in the brain. However, investigating CFS in Magneto- and Electroencephalography (MEG/EEG) is hampered by the presence of spurious neuronal interactions due to the non-sinusoidal waveshape of brain oscillations. Such waveshape gives rise to the presence of oscillatory harmonics mimicking genuine neuronal oscillations. Until recently, however, there has been no methodology for removing these harmonics from neuronal data. In order to address this long-standing challenge, we introduce a novel method (called HARMOnic miNImization - Harmoni) that removes the signal components which can be harmonics of a non-sinusoidal signal. Harmoni's working principle is based on the presence of CFS between harmonic components and the fundamental component of a non-sinusoidal signal. We extensively tested Harmoni in realistic EEG simulations. The simulated couplings between the source signals represented genuine and spurious CFS and within-frequency phase synchronization. Using diverse evaluation criteria, including ROC analyses, we showed that the within- and cross-frequency spurious interactions are suppressed significantly, while the genuine activities are not affected. Additionally, we applied Harmoni to real resting-state EEG data revealing intricate remote connectivity patterns which are usually masked by the spurious connections. Given the ubiquity of non-sinusoidal neuronal oscillations in electrophysiological recordings, Harmoni is expected to facilitate novel insights into genuine neuronal interactions in various research fields, and can also serve as a steppingstone towards the development of further signal processing methods aiming at refining within- and cross-frequency synchronization in electrophysiological recordings.

Differential effects of anodal and dual tDCS on sensorimotor functions in chronic hemiparetic stroke patients.

BACKGROUND AND PURPOSE: Previous tDCS studies in chronic stroke patients reported highly inconsistent effects on sensorimotor functions. Underlying reasons could be the selection of different kinematic parameters across studies and for different tDCS setups. We reasoned that tDCS may not simply induce global changes in a beneficial-adverse dichotomy, but rather that different sensorimotor kinematics are differentially affected. Furthermore, the often-postulated higher efficacy of bilateral-dual (bi-tDCS) over unilateral-anodal (ua-tDCS) could not yet be demonstrated consistently either. We investigated the effects of both setups on a wider range of kinematic parameters from standardized robotic tasks in patients with chronic stroke. METHODS: Twenty-four patients with arm hemiparesis received tDCS (20min, 1 mA) concurrent to kinematic assessments in a sham-controlled, cross-over and double-blind clinical trial. Performance was measured on four sensorimotor tasks (reaching, proprioception, cooperative and independent bimanual coordination) from which 30 parameters were extracted. On the group-level, the patterns of changes relative to sham were assessed using paired-samples t-tests and classified as (1) performance increases, (2) decreases and (3) non-significant differences. Correlations between parametric change scores were calculated for each task to assess effects on the individual-level. RESULTS: Both setups induced complex effect patterns with varying proportions of performance increases and decreases. On the group-level, more increases were induced in the reaching and coordination tasks while proprioception and bimanual cooperation were overall negatively affected. Bi-tDCS induced more performance increases and less decreases compared to ua-tDCS. Changes across parameters occurred more homogeneously under bi-tDCS than ua-tDCS, which induced a larger proportion of performance trade-offs. CONCLUSIONS: Our data demonstrate profound tDCS effects on sensorimotor functions post-stroke, lending support for more pronounced and favorable effects of bi-tDCS compared to ua-tDCS. However, no uniformly beneficial pattern was identified. Instead, the modulations varied depending on the task and electrode setup, with increases in certain parameters occurring at the expense of decreases in others.