A bicoherence approach to analyze multi-dimensional cross-frequency coupling in EEG/MEG data

We introduce a blockwise generalisation of the Antisymmetric Cross-Bicoherence (ACB), a statistical method based on bispectral analysis. The Multi-dimensional ACB (MACB) is an approach that aims at detecting quadratic lagged phase-interactions between vector time series in the frequency domain. Such a coupling can be empirically observed in functional neuroimaging data, e.g., in electro/magnetoencephalographic signals. MACB is invariant under orthogonal trasformations of the data, which makes it independent, e.g., on the choice of the physical coordinate system in the neuro-electromagnetic inverse procedure. In extensive synthetic experiments, we prove that MACB performance is significantly better than that obtained by ACB. Specifically, the shorter the data length, or the higher the dimension of the single data space, the larger the difference between the two methods.

Frontal theta oscillations during emotion regulation in people with borderline personality disorder

BACKGROUND

Borderline personality disorder (BPD) is a severe psychiatric disorder conceptualised as a disorder of emotion regulation. Emotion regulation has been linked to a frontolimbic network comprising the dorsolateral prefrontal cortex and the amygdala, which apparently synchronises its activity via oscillatory coupling in the theta frequency range.

AIMS

To analyse whether there are distinct differences in theta oscillatory coupling in frontal brain regions between individuals with BPD and matched controls during emotion regulation by cognitive reappraisal.

METHOD

Electroencephalogram (EEG) recordings were performed in 25 women diagnosed with BPD and 25 matched controls during a cognitive reappraisal task in which participants were instructed to downregulate negative emotions evoked by aversive visual stimuli. Between- and within-group time-frequency analyses were conducted to analyse regulation-associated theta activity (3.5-8.5 Hz).

RESULTS

Oscillatory theta activity differed between the participants with BPD and matched controls during cognitive reappraisal. Regulation-associated theta increases were lower in frontal regions in the BPD cohort compared with matched controls. Functional connectivity analysis for regulation-associated changes in the theta frequency band revealed a lower multivariate interaction measure (MIM) increase in frontal brain regions in persons with BPD compared with matched controls.

CONCLUSIONS

Our findings support the notion of alterations in a frontal theta network in BPD, which may be underlying core symptoms of the disorder such as deficits in emotion regulation. The results add to the growing body of evidence for altered oscillatory brain dynamics in psychiatric populations, which might be investigated as individualised treatment targets using non-invasive stimulation methods.

Distinct and common mechanisms of cross-model semantic conflict and response conflict in an auditory relevant task

The mechanisms of semantic conflict and response conflict in the Stroop task have mainly been investigated in the visual modality. However, the understanding of these mechanisms in cross-modal modalities remains limited. In this electroencephalography (EEG) study, an audiovisual 2-1 mapping Stroop task was utilized to investigate whether distinct and/or common neural mechanisms underlie cross-modal semantic conflict and response conflict. The response time data showed significant effects on both cross-modal semantic and response conflicts. Interestingly, the magnitude of semantic conflict was found to be smaller in the fast response time bins than in the slow response time bins, whereas no such difference was observed for response conflict. The EEG data demonstrated that cross-modal semantic conflict specifically increased the N450 amplitude. However, cross-modal response conflict specifically enhanced theta band power and theta phase synchronization between the medial frontal cortex (MFC) and lateral prefrontal electrodes as well as between the MFC and motor electrodes. In addition, both cross-modal semantic conflict and response conflict led to a decrease in P3 amplitude. Taken together, these findings provide cross-modal evidence for domain-specific mechanism in conflict detection and suggest both domain-specific and domain-general mechanisms exist in conflict resolution.

Leading and following: Noise differently affects semantic and acoustic processing during naturalistic speech comprehension

Despite the distortion of speech signals caused by unavoidable noise in daily life, our ability to comprehend speech in noisy environments is relatively stable. However, the neural mechanisms underlying reliable speech-in-noise comprehension remain to be elucidated. The present study investigated the neural tracking of acoustic and semantic speech information during noisy naturalistic speech comprehension. Participants listened to narrative audio recordings mixed with spectrally matched stationary noise at three signal-to-ratio (SNR) levels (no noise, 3 dB, -3 dB), and 60-channel electroencephalography (EEG) signals were recorded. A temporal response function (TRF) method was employed to derive event-related-like responses to the continuous speech stream at both the acoustic and the semantic levels. Whereas the amplitude envelope of the naturalistic speech was taken as the acoustic feature, word entropy and word surprisal were extracted via the natural language processing method as two semantic features. Theta-band frontocentral TRF responses to the acoustic feature were observed at around 400 ms following speech fluctuation onset over all three SNR levels, and the response latencies were more delayed with increasing noise. Delta-band frontal TRF responses to the semantic feature of word entropy were observed at around 200 to 600 ms leading to speech fluctuation onset over all three SNR levels. The response latencies became more leading with increasing noise and decreasing speech comprehension and intelligibility. While the following responses to speech acoustics were consistent with previous studies, our study revealed the robustness of leading responses to speech semantics, which suggests a possible predictive mechanism at the semantic level for maintaining reliable speech comprehension in noisy environments.

EEG-based speaker-listener neural coupling reflects speech-selective attentional mechanisms beyond the speech stimulus

When we pay attention to someone, do we focus only on the sound they make, the word they use, or do we form a mental space shared with the speaker we want to pay attention to? Some would argue that the human language is no other than a simple signal, but others claim that human beings understand each other because they form a shared mental ground between the speaker and the listener. Our study aimed to explore the neural mechanisms of speech-selective attention by investigating the electroencephalogram-based neural coupling between the speaker and the listener in a cocktail party paradigm. The temporal response function method was employed to reveal how the listener was coupled to the speaker at the neural level. The results showed that the neural coupling between the listener and the attended speaker peaked 5 s before speech onset at the delta band over the left frontal region, and was correlated with speech comprehension performance. In contrast, the attentional processing of speech acoustics and semantics occurred primarily at a later stage after speech onset and was not significantly correlated with comprehension performance. These findings suggest a predictive mechanism to achieve speaker-listener neural coupling for successful speech comprehension.

Novel multivariate methods to track frequency shifts of neural oscillations in EEG/MEG recordings

Instantaneous and peak frequency changes in neural oscillations have been linked to many perceptual, motor, and cognitive processes. Yet, the majority of such studies have been performed in sensor space and only occasionally in source space. Furthermore, both terms have been used interchangeably in the literature, although they do not reflect the same aspect of neural oscillations. In this paper, we discuss the relation between instantaneous frequency, peak frequency, and local frequency, the latter also known as spectral centroid. Furthermore, we propose and validate three different methods to extract source signals from multichannel data whose (instantaneous, local, or peak) frequency estimate is maximally correlated to an experimental variable of interest. Results show that the local frequency might be a better estimate of frequency variability than instantaneous frequency under conditions with low signal-to-noise ratio. Additionally, the source separation methods based on local and peak frequency estimates, called LFD and PFD respectively, provide more stable estimates than the decomposition based on instantaneous frequency. In particular, LFD and PFD are able to recover the sources of interest in simulations performed with a realistic head model, providing higher correlations with an experimental variable than multiple linear regression. Finally, we also tested all decomposition methods on real EEG data from a steady-state visual evoked potential paradigm and show that the recovered sources are located in areas similar to those previously reported in other studies, thus providing further validation of the proposed methods.

Comparison of transcranial brain stimulation approaches: prefrontal theta alternating current stimulation enhances working memory performance

Introduction

One of the most important cognitive functions in our everyday life is the working memory (WM). In several neuropsychiatric diseases such as ADHD or schizophrenia WM deficits can be observed, making it an attractive target for non-invasive brain stimulation methods like transcranial electrical stimulation (tES). However, the literature shows rather heterogeneous results of tES effects on WM performance. fMRI meta-analyses have identified a WM network including frontoparietal brain areas such as the dorsolateral prefrontal cortex (DLPFC) and the posterior parietal cortex (PPC). Neurophysiological studies revealed oscillatory activity in the theta band frequency range to be of crucial functional relevance for WM processes. Based on this, transcranial alternating current stimulation (tACS) in the theta frequency range targeting DLPFC and PPC in a spatially optimized way might further improve effects of tES on WM performance.

Methods

Sixteen healthy subjects were stimulated with varying stimulation settings on four different days in a counterbalanced within-subject design. These setups included the application of (1) tACS with a frequency of 5 Hz (theta frequency range) over the left DLPFC and (2) the right superior parietal cortex, (3) transcranial direct current stimulation (tDCS) of the DLPFC and (4) a sham stimulation condition during the online performance of a visual delayed-match-to-sample task with varying working memory load. We introduce a procedure to calculate an optimal tES model revealing optimized high-density setups for the present study for 3 cathodes and 1 anode and stimulation currents of 1.5 mA.

Results

A significant interaction effect of stimulation type and load condition on working memory capacity was found. This was reflected by a significant improvement of WM performance in the high load condition during tACS over the left DLPFC compared with sham stimulation, which was not the case for our parietal tACS or tDCS setup.

Discussion

Working memory performance can be improved with optimized high-definition tACS with a frequency of 5 Hz over the left DLPFC. The conception of different mechanisms underlying transcranial electrical stimulation with alternating and direct currents is supported by these results. Patients suffering from working memory impairments due to neuropsychiatric diseases might potentially benefit from this brain stimulation approach.

Topological reorganization of brain network might contribute to the resilience of cognitive functioning in mildly disabled relapsing remitting multiple sclerosis

Multiple sclerosis (MS) is an inflammatory and demyelinating disease which leads to impairment in several functional systems including cognition. Alteration of brain networks is linked to disability and its progression. However, results are mostly cross-sectional and yet contradictory as putative adaptive and maladaptive mechanisms were found. Here, we aimed to explore longitudinal reorganization of brain networks over 2-years by combining diffusion tensor imaging (DTI), resting-state functional MRI (fMRI), magnetoencephalography (MEG), and a comprehensive neuropsychological-battery. In 37 relapsing-remitting MS (RRMS) and 39 healthy-controls, cognition remained stable over-time. We reconstructed network models based on the three modalities and analyzed connectivity in relation to the hierarchical topology and functional subnetworks. Network models were compared across modalities and in their association with cognition using linear-mixed-effect-regression models. Loss of hub connectivity and global reduction was observed on a structural level over-years (p < .010), which was similar for functional MEG-networks but not for fMRI-networks. Structural hub connectivity increased in controls (p = .044), suggesting a physiological mechanism of healthy aging. Despite a general loss in structural connectivity in RRMS, hub connectivity was preserved (p = .002) over-time in default-mode-network (DMN). MEG-networks were similar to DTI and weakly correlated with fMRI in MS (p < .050). Lower structural (β between .23-.33) and both lower (β between .40-.59) and higher functional connectivity (β = -.54) in DMN was associated with poorer performance in attention and memory in RRMS (p < .001). MEG-networks involved no association with cognition. Here, cognitive stability despite ongoing neurodegeneration might indicate a resilience mechanism of DMN hubs mimicking a physiological reorganization observed in healthy aging.

Speaker-listener neural coupling reveals a right-lateralized mechanism for non-native speech-in-noise comprehension

While the increasingly globalized world has brought more and more demands for non-native language communication, the prevalence of background noise in everyday life poses a great challenge to non-native speech comprehension. The present study employed an interbrain approach based on functional near-infrared spectroscopy (fNIRS) to explore how people adapt to comprehend non-native speech information in noise. A group of Korean participants who acquired Chinese as their non-native language was invited to listen to Chinese narratives at 4 noise levels (no noise, 2 dB, -6 dB, and - 9 dB). These narratives were real-life stories spoken by native Chinese speakers. Processing of the non-native speech was associated with significant fNIRS-based listener-speaker neural couplings mainly over the right hemisphere at both the listener's and the speaker's sides. More importantly, the neural couplings from the listener's right superior temporal gyrus, the right middle temporal gyrus, as well as the right postcentral gyrus were found to be positively correlated with their individual comprehension performance at the strongest noise level (-9 dB). These results provide interbrain evidence in support of the right-lateralized mechanism for non-native speech processing and suggest that both an auditory-based and a sensorimotor-based mechanism contributed to the non-native speech-in-noise comprehension.

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.

Preparatory delta phase response is correlated with naturalistic speech comprehension performance

While human speech comprehension is thought to be an active process that involves top-down predictions, it remains unclear how predictive information is used to prepare for the processing of upcoming speech information. We aimed to identify the neural signatures of the preparatory processing of upcoming speech. Participants selectively attended to one of two competing naturalistic, narrative speech streams, and a temporal response function (TRF) method was applied to derive event-related-like neural responses from electroencephalographic data. The phase responses to the attended speech at the delta band (1-4 Hz) were correlated with the comprehension performance of individual participants, with a latency of - 200-0 ms relative to the onset of speech amplitude envelope fluctuations over the fronto-central and left-lateralized parietal electrodes. The phase responses to the attended speech at the alpha band also correlated with comprehension performance but with a latency of 650-980 ms post-onset over the fronto-central electrodes. Distinct neural signatures were found for the attentional modulation, taking the form of TRF-based amplitude responses at a latency of 240-320 ms post-onset over the left-lateralized fronto-central and occipital electrodes. Our findings reveal how the brain gets prepared to process an upcoming speech in a continuous, naturalistic speech context.

Speaker-Listener Neural Coupling Reveals an Adaptive Mechanism for Speech Comprehension in a Noisy Environment

Comprehending speech in noise is an essential cognitive skill for verbal communication. However, it remains unclear how our brain adapts to the noisy environment to achieve comprehension. The present study investigated the neural mechanisms of speech comprehension in noise using an functional near-infrared spectroscopy-based inter-brain approach. A group of speakers was invited to tell real-life stories. The recorded speech audios were added with meaningless white noise at four signal-to-noise levels and then played to listeners. Results showed that speaker-listener neural couplings of listener's left inferior frontal gyri (IFG), that is, sensorimotor system, and right middle temporal gyri (MTG), angular gyri (AG), that is, auditory system, were significantly higher in listening conditions than in the baseline. More importantly, the correlation between neural coupling of listener's left IFG and the comprehension performance gradually became more positive with increasing noise level, indicating an adaptive role of sensorimotor system in noisy speech comprehension; however, the top behavioral correlations for the coupling of listener's right MTG and AG were only obtained in mild noise conditions, indicating a different and less robust mechanism. To sum up, speaker-listener coupling analysis provides added value and new sight to understand the neural mechanism of speech-in-noise comprehension.

Circuit mechanisms for the chemical modulation of cortex-wide network interactions and behavioral variability

Influential theories postulate distinct roles of catecholamines and acetylcholine in cognition and behavior. However, previous physiological work reported similar effects of these neuromodulators on the response properties (specifically, the gain) of individual cortical neurons. Here, we show a double dissociation between the effects of catecholamines and acetylcholine at the level of large-scale interactions between cortical areas in humans. A pharmacological boost of catecholamine levels increased cortex-wide interactions during a visual task, but not rest. An acetylcholine boost decreased interactions during rest, but not task. Cortical circuit modeling explained this dissociation by differential changes in two circuit properties: the local excitation-inhibition balance (more strongly increased by catecholamines) and intracortical transmission (more strongly reduced by acetylcholine). The inferred catecholaminergic mechanism also predicted noisier decision-making, which we confirmed for both perceptual and value-based choice behavior. Our work highlights specific circuit mechanisms for shaping cortical network interactions and behavioral variability by key neuromodulatory systems.

Ketamine Alters Functional Gamma and Theta Resting-State Connectivity in Healthy Humans: Implications for Schizophrenia Treatment Targeting the Glutamate System

Disturbed functional connectivity is assumed to cause neurocognitive deficits in patients suffering from schizophrenia. A Glutamate N-methyl-D-aspartate receptor (NMDAR) dysfunction has been suggested as a possible mechanism underlying altered connectivity in schizophrenia, especially in the gamma- and theta-frequency range. The present study aimed to investigate the effects of the NMDAR-antagonist ketamine on resting-state power, functional connectivity, and schizophrenia-like psychopathological changes in healthy volunteers. In a placebo-controlled crossover design, 25 healthy subjects were recorded using resting-state 64-channel-electroencephalography (EEG) (eyes closed). The imaginary coherence-based Multivariate Interaction Measure (MIM) was used to measure gamma and theta connectivity across 80 cortical regions. The network-based statistic was applied to identify involved networks under ketamine. Psychopathology was assessed with the Positive and Negative Syndrome Scale (PANSS) and the 5-Dimensional Altered States of Consciousness Rating Scale (5D-ASC). Ketamine caused an increase in all PANSS (p < 0.001) as well as 5D-ASC scores (p < 0.01). Significant increases in resting-state gamma and theta power were observed under ketamine compared to placebo (p < 0.05). The source-space analysis revealed two distinct networks with an increased mean functional gamma- or theta-band connectivity during the ketamine session. The gamma-network consisted of midline regions, the cuneus, the precuneus, and the bilateral posterior cingulate cortices, while the theta-band network involved the Heschl gyrus, midline regions, the insula, and the middle cingulate cortex. The current source density (CSD) within the gamma-band correlated negatively with the PANSS negative symptom score, and the activity within the gamma-band network correlated negatively with the subjective changed meaning of percepts subscale of the 5D-ASC. These results are in line with resting-state patterns seen in people who have schizophrenia and argue for a crucial role of the glutamate system in mediating dysfunctional gamma- and theta-band-connectivity in schizophrenia. Resting-state networks could serve as biomarkers for the response to glutamatergic drugs or drug development efforts within the glutamate system.

Defining the filter parameters for phase-amplitude coupling from a bispectral point of view

BACKGROUND

Two measures of cross-frequency coupling are phase-amplitude coupling (PAC) and bicoherence. The estimation of PAC with meaningful bandwidth for the high-frequency amplitude is crucial in order to avoid misinterpretations. While recommendations on the bandwidth of PAC's amplitude component exist, there is no consensus yet. Theoretical relationships between PAC and bicoherence can provide insights on how to set PAC's filters.

NEW METHOD

To illustrate this, PAC estimated from simulated and empirical data are compared to the bispectrum. We used simulations replicated from earlier studies and empirical data from human electro-encephalography and rat local field potentials. PAC's amplitude component was estimated using a filter bandwidth with a ratio of (1) 2:1, (2) 1:1, or (3) 0.5:1 relative to the phase frequency.

RESULTS

For both simulated and empirical data, PAC was smeared over a broad frequency range and not present when the estimates comprised a 2:1- and 0.5:1-ratio, respectively. In contrast, the 1:1-ratio accurately avoids smearing and results in clear signals of cross-frequency coupling. Bicoherence estimates were found to be essentially identical to PAC calculated with the recommended frequency setting.

COMPARISON WITH EXISTING METHOD(S)

Earlier recommendations on filter settings of PAC lead to estimates which are smeared in the frequency domain, which makes it difficult to identify cross-frequency coupling of neural processes operating in narrow frequency bands.

CONCLUSIONS

We conclude that smearing of PAC estimates can be avoided with a different choice of filter settings by theoretically relating PAC to bicoherence.

Modality-specific neural mechanisms of cognitive control in a Stroop-like task

The successful resolution of ever-changing conflicting contexts requires efficient cognitive control. Previous studies have found similar neural patterns in conflict processing for different modalities using an event-related potential (ERP) approach and have concluded that cognitive control is supramodal. However, recent behavioral studies have found that conflict adaptation (a phenomenon with the reduction of congruency effect in the current trial after an incongruent trial as compared with a congruent trial) could not transfer across visual and auditory modalities and suggested that cognitive control is modality-specific, challenging the supramodal view. These discrepancies may have also arisen from methodological differences across studies. The current study examined the electroencephalographic profiles of a Stroop-like task to elucidate the modality-specific neural mechanisms of cognitive control. Participants were instructed to respond to a target always coming from the visual modality while disregarding the distractor coming from either the auditory or the visual modality. The results revealed significant congruency effects on both behavioral indices, i.e., reaction time and error rate, and ERP components, including the P3 and the conflict slow potential. Besides, the congruency effects on the amplitude of the P3 showed a negative correlation with reaction time, indicating an intrinsic link between these neural and behavioral indices. Furthermore, in the modality-repetition condition, conflict adaptation effects were significant on both reaction time and P3 amplitude, and the reaction time could be predicted by the P3 amplitude, while such effects were not observed in the modality-alternation condition. The time-frequency analysis also showed that conflict adaptation occurred in the modality-repetition condition, but not in the modality-alternation condition in low frequency bands, including the theta (4-8 Hz), alpha (8-12 Hz), and beta1 (12-20 Hz) bands. Taken together, our results revealed modality-specific patterns of the conflict adaptation effects on the P3 amplitude and oscillatory power (in theta, alpha, and beta1 bands), providing neural evidence for the modality specificity of cognitive control and expanding the boundaries of cognitive control.

Beyond Pairwise Interactions: The Totally Antisymmetric Part of the Bispectrum as Coupling Measure of at Least Three Interacting Sources

In this paper we make two contributions to the analysis of brain oscillations with CFC techniques. First, we introduce a new bispectral CFC measure which is selective to couplings between three or more brain sources. This measure can be derived from ordinary cross-bispectra by performing a total-antisymmetrization operation on them. Significant coupling values can then be attributed to at least three interacting signals. This selectivity to the number of sources can be helpful to test hypotheses on the number of brain sources involved in the generation of commonly observed brain oscillations, such as the alpha rhythm. In a second step we present the correct empirical distribution for the coupling measure, which is necessary to properly assess the significance of coupling results. More importantly however, this corrected statistic is not limited to our particular measure, but holds for all complex-valued coupling estimators. We illustrate how the very common misassumption of empirical normality of such estimators can lead to a systematic underestimation of p-values, the breakdown of multiple comparison control procedures and in consequence a drastic inflation of the number of false positives.

Mathematical Relations Between Measures of Brain Connectivity Estimated From Electrophysiological Recordings for Gaussian Distributed Data

A large variety of methods exist to estimate brain coupling in the frequency domain from electrophysiological data measured, e.g., by EEG and MEG. Those data are to reasonable approximation, though certainly not perfectly, Gaussian distributed. This work is based on the well-known fact that for Gaussian distributed data, the cross-spectrum completely determines all statistical properties. In particular, for an infinite number of data, all normalized coupling measures at a given frequency are a function of complex coherency. However, it is largely unknown what the functional relations are. We here present those functional relations for six different measures: the weighted phase lag index, the phase lag index, the absolute value and imaginary part of the phase locking value (PLV), power envelope correlation, and power envelope correlation with correction for artifacts of volume conduction. With the exception of PLV, the final results are simple closed form formulas. In an excursion we also discuss differences between short time Fourier transformation and Hilbert transformation for estimations in the frequency domain. We tested in simulations of linear and non-linear dynamical systems and for empirical resting state EEG on sensor level to what extent a model, namely the respective function of coherency, can explain the observed couplings. For empirical data we found that for measures of phase-phase coupling deviations from the model are in general minor, while power envelope correlations systematically deviate from the model for all frequencies. For power envelope correlation with correction for artifacts of volume conduction the model cannot explain the observed couplings at all. We also analyzed power envelope correlation as a function of time and frequency in an event related experiment using a stroop reaction task and found significant event related deviations mostly in the alpha range.

Partial Similarity Reveals Dynamics in Brainstem-Midbrain Networks during Trigeminal Nociception

Imaging studies help us understand the important role of brainstem and midbrain regions in human trigeminal pain processing without solving the question of how these regions actually interact. In the current study, we describe this connectivity and its dynamics during nociception with a novel analytical approach called Partial Similarity (PS). We developed PS specifically to estimate the communication between individual hubs of the network in contrast to the overall communication within that network. Partial Similarity works on trial-to-trial variance of neuronal activity acquired with functional magnetic resonance imaging. It discovers direct communication between two hubs considering the remainder of the network as confounds. A similar method to PS is Representational Similarity, which works with ordinary correlations and does not consider any external influence on the communication between two hubs. Particularly the combination of Representational Similarity and Partial Similarity analysis unravels brainstem dynamics involved in trigeminal pain using the spinal trigeminal nucleus (STN)—the first relay station of peripheral trigeminal input—as a seed region. The combination of both methods can be valuable tools in discovering the network dynamics in fMRI and an important instrument for future insight into the nature of various neurological diseases like primary headaches.

Aerobic Exercise Induces Functional and Structural Reorganization of CNS Networks in Multiple Sclerosis: A Randomized Controlled Trial

Objectives: Evidence from animal studies suggests that aerobic exercise may promote neuroplasticity and could, therefore, provide therapeutic benefits for neurological diseases such as multiple sclerosis (MS). However, the effects of exercise in human CNS disorders on the topology of brain networks, which might serve as an outcome at the interface between biology and clinical performance, remain poorly understood.

Methods: We investigated functional and structural networks in patients with relapsing-remitting MS in a clinical trial of standardized aerobic exercise. Fifty-seven patients were randomly assigned to moderate-intensity exercise for 3 months or a non-exercise control group. We reconstructed functional networks based on resting-state functional magnetic resonance imaging (MRI) and used probabilistic tractography on diffusion-weighted imaging data for structural networks.

Results: At baseline, compared to 30 healthy controls, patients exhibited decreased structural connectivity that was most pronounced in hub regions of the brain. Vice versa, functional connectivity was increased in hubs. After 3 months, we observed hub independent increased functional connectivity in the exercise group while the control group presented a loss of functional hub connectivity. On a structural level, the control group remained unchanged, while the exercise group had also increased connectivity. Increased clustering of hubs indicates a better structural integration and internal connectivity at the top of the network hierarchy.

Conclusion: Increased functional connectivity of hubs contrasts a loss of structural connectivity in relapsing-remitting MS. Under an exercise condition, a further hub independent increase of functional connectivity seems to translate in higher structural connectivity of the whole brain.

Intrinsic 40Hz-phase asymmetries predict tACS effects during conscious auditory perception

Synchronized oscillatory gamma-band activity (30-100Hz) has been suggested to constitute a key mechanism to dynamically orchestrate sensory information integration across multiple spatio-temporal scales. We here tested whether interhemispheric functional connectivity and ensuing auditory perception can selectively be modulated by high-density transcranial alternating current stimulation (HD-tACS). For this purpose, we applied multi-site HD-tACS at 40Hz bilaterally with a phase lag of 180° and recorded a 64-channel EEG to study the oscillatory phase dynamics at the source-space level during a dichotic listening (DL) task in twenty-six healthy participants. In this study, we revealed an oscillatory phase signature at 40Hz which reflects different temporal profiles of the phase asymmetries during left and right ear percept. Here we report that 180°-tACS did not affect the right ear advantage during DL at group level. However, a follow-up analysis revealed that the intrinsic phase asymmetries during sham-tACS determined the directionality of the behavioral modulations: While a shift to left ear percept was associated with augmented interhemispheric asymmetry (closer to 180°), a shift to right ear processing was elicited in subjects with lower asymmetry (closer to 0°). Crucially, the modulation of the interhemispheric network dynamics depended on the deviation of the tACS-induced phase-lag from the intrinsic phase asymmetry. Our characterization of the oscillatory network trends is giving rise to the importance of phase-specific gamma-band coupling during ambiguous auditory perception, and emphasizes the necessity to address the inter-individual variability of phase asymmetries in future studies by tailored stimulation protocols.

Visual input drives increased occipital responsiveness and harmonized oscillations in multiple cortical areas in migraineurs

Migraineurs are hypersensitive for most sensory domains like visual, auditory or somatosensory processing even outside of attacks. This behavioral peculiarity is mirrored by findings of cortical hyper-responsivity already in the interictal state. Using repetitive visual stimulation to elicit steady state visually evoked potentials (SSVEP) in 30 interictal episodic migraineurs and 30 controls we show hyper-responsivity of the visual cortex in the migraineurs. Additionally, the occipital regions were remarkably stronger coupled to the temporal, premotor and the anterior cingulate cortex than in headache free controls. These data suggest harmonized oscillations of different cortical areas as a response to visual input which might be driven by the cuneus. Furthermore, the increased coupling is modulated by the current state of the migraine cycle as the coupling was significantly stronger in patients with longer interictal periods.

•

Migraineurs visual cortex compared to controls is hyper-responsiveness in response to repetitive visual stimulation.

•

The hyper-responsiveness is stronger coupled to temporal, premotor and anterior cingulate cortex than in controls.

•

This increased coupling is modulated by the current state of the migraine cycle.

•

Our data suggest that visual input in migraineurs leads to harmonized oscillations of multiple cortical areas.

Catecholamines alter the intrinsic variability of cortical population activity and perception

The ascending modulatory systems of the brain stem are powerful regulators of global brain state. Disturbances of these systems are implicated in several major neuropsychiatric disorders. Yet, how these systems interact with specific neural computations in the cerebral cortex to shape perception, cognition, and behavior remains poorly understood. Here, we probed into the effect of two such systems, the catecholaminergic (dopaminergic and noradrenergic) and cholinergic systems, on an important aspect of cortical computation: its intrinsic variability. To this end, we combined placebo-controlled pharmacological intervention in humans, recordings of cortical population activity using magnetoencephalography (MEG), and psychophysical measurements of the perception of ambiguous visual input. A low-dose catecholaminergic, but not cholinergic, manipulation altered the rate of spontaneous perceptual fluctuations as well as the temporal structure of "scale-free" population activity of large swaths of the visual and parietal cortices. Computational analyses indicate that both effects were consistent with an increase in excitatory relative to inhibitory activity in the cortical areas underlying visual perceptual inference. We propose that catecholamines regulate the variability of perception and cognition through dynamically changing the cortical excitation-inhibition ratio. The combined readout of fluctuations in perception and cortical activity we established here may prove useful as an efficient and easily accessible marker of altered cortical computation in neuropsychiatric disorders.

Spectral fingerprints of large-scale cortical dynamics during ambiguous motion perception

Ambiguous stimuli have been widely used to study the neuronal correlates of consciousness. Recently, it has been suggested that conscious perception might arise from the dynamic interplay of functionally specialized but widely distributed cortical areas. While previous research mainly focused on phase coupling as a correlate of cortical communication, more recent findings indicated that additional coupling modes might coexist and possibly subserve distinct cortical functions. Here, we studied two coupling modes, namely phase and envelope coupling, which might differ in their origins, putative functions and dynamics. Therefore, we recorded 128-channel EEG while participants performed a bistable motion task and utilized state-of-the-art source-space connectivity analysis techniques to study the functional relevance of different coupling modes for cortical communication. Our results indicate that gamma-band phase coupling in extrastriate visual cortex might mediate the integration of visual tokens into a moving stimulus during ambiguous visual stimulation. Furthermore, our results suggest that long-range fronto-occipital gamma-band envelope coupling sustains the horizontal percept during ambiguous motion perception. Additionally, our results support the idea that local parieto-occipital alpha-band phase coupling controls the inter-hemispheric information transfer. These findings provide correlative evidence for the notion that synchronized oscillatory brain activity reflects the processing of sensory input as well as the information integration across several spatiotemporal scales. The results indicate that distinct coupling modes are involved in different cortical computations and that the rich spatiotemporal correlation structure of the brain might constitute the functional architecture for cortical processing and specific multi-site communication. Hum Brain Mapp 37:4099-4111, 2016. © 2016 Wiley Periodicals, Inc.

Acupuncture analgesia involves modulation of pain-induced gamma oscillations and cortical network connectivity

Recent studies support the view that cortical sensory, limbic and executive networks and the autonomic nervous system might interact in distinct manners under the influence of acupuncture to modulate pain. We performed a double-blind crossover design study to investigate subjective ratings, EEG and ECG following experimental laser pain under the influence of sham and verum acupuncture in 26 healthy volunteers. We analyzed neuronal oscillations and inter-regional coherence in the gamma band of 128-channel-EEG recordings as well as heart rate variability (HRV) on two experimental days. Pain ratings and pain-induced gamma oscillations together with vagally-mediated power in the high-frequency bandwidth (vmHF) of HRV decreased significantly stronger during verum than sham acupuncture. Gamma oscillations were localized in the prefrontal cortex (PFC), mid-cingulate cortex (MCC), primary somatosensory cortex and insula. Reductions of pain ratings and vmHF-power were significantly correlated with increase of connectivity between the insula and MCC. In contrast, connectivity between left and right PFC and between PFC and insula correlated positively with vmHF-power without a relationship to acupuncture analgesia. Overall, these findings highlight the influence of the insula in integrating activity in limbic-saliency networks with vagally mediated homeostatic control to mediate antinociception under the influence of acupuncture.

Dynamic reconfiguration of cortical functional connectivity across brain states

Throughout each day, the brain displays transient changes in state, as evidenced by shifts in behavior and vigilance. While the electrophysiological correlates of brain states have been studied for some time, it remains unclear how large-scale cortico-cortical functional connectivity systematically reconfigures across states. Here, we investigate state-dependent shifts in cortical functional connectivity by recording local field potentials (LFPs) during spontaneous behavioral transitions in the ferret using chronically implanted micro-electrocorticographic (µECoG) arrays positioned over occipital, parietal, and temporal cortical regions. To objectively classify brain state, we describe a data-driven approach that projects time-varying LFP spectral properties into brain state space. Distinct brain states displayed markedly different patterns of cross-frequency phase-amplitude coupling and inter-electrode phase synchronization across several LFP frequency bands. The largest across-state differences in functional connectivity were observed between periods of presumed slow-wave and rapid-eye-movement-sleep/active-state, which were characterized by the contrasting phenomena of cortical network fragmentation and global synchronization, respectively. Collectively, our data provide strong evidence that large-scale functional interactions in the brain dynamically reconfigure across behavioral states.

The Callosal Relay Model of Interhemispheric Communication: New Evidence from Effective Connectivity Analysis

Interhemispheric auditory connectivity via the corpus callosum has been demonstrated to be important for normal speech processing. According to the callosal relay model, directed information flow from the right to the left auditory cortex has been suggested, but this has not yet been proven. For this purpose, 33 healthy participants were investigated with 64-channel EEG while performing the dichotic listening task in which two different consonant–vowel syllables were presented simultaneously to the left (LE) and right ear (RE). eLORETA source estimation was used to investigate the functional (lagged phase synchronization/LPS) and effective (isolated effective coherence/ICoh) connectivity between right and left primary (PAC) and secondary auditory cortices (SAC) in the gamma-band (30–100 Hz) during right and left ear reports. The major finding was a significantly increased effective connectivity in the gamma-band from the right to the left SAC during conscious perception of LE stimuli. In addition, effective and functional connectivity was significantly enhanced during LE as compared to RE reports. These findings give novel insight into transcallosal information transfer during auditory perception by showing that LE performance requires causal interhemispheric inputs from the right to the left auditory cortices, and that this interaction is mediated by synchronized gamma-band oscillations.

Quality assessment of 3D visualizations with vertical disparity: An ERP approach

In an objective approach for the assessment of quality of experience the neural correlates of EEG data are studied when stereoscopic images are presented in three different conditions containing vertical disparity. These conditions are compared to a similar image in 2D both on the channel level by studying the ERP components and on the source level by the localization of the corresponding ERP component. Our findings posit that P1 component in the occipital cortex has significantly increased in amplitude for 3D condition without vertical disparity compared to the 2D condition. According to previous studies, this component increases when depth information are added to the stimulus which is in line with our findings. However the amplitude of this component has significantly decreased for 3D condition with maximum vertical disparity compared to the 3D condition without vertical disparity. We have concluded that the perception of stereoscopic depth by subjects have decreased in this case due to the distortion introduced by vertical disparity. The underlying sources corresponding to P1 component are localized. Except for the power differences, the source locations do not differ for different conditions.

Objective quality assessment of stereoscopic images with vertical disparity using EEG

OBJECTIVE

Neurophysiological correlates of vertical disparity in 3D images are studied in an objective approach using EEG technique. These disparities are known to negatively affect the quality of experience and to cause visual discomfort in stereoscopic visualizations.

APPROACH

We have presented four conditions to subjects: one in 2D and three conditions in 3D, one without vertical disparity and two with different vertical disparity levels. Event related potentials (ERPs) are measured for each condition and the differences between ERP components are studied. Analysis is also performed on the induced potentials in the time frequency domain.

MAIN RESULTS

Results show that there is a significant increase in the amplitude of P1 components in 3D conditions in comparison to 2D. These results are consistent with previous studies which have shown that P1 amplitude increases due to the depth perception in 3D compared to 2D. However the amplitude is significantly smaller for maximum vertical disparity (3D-3) in comparison to 3D with no vertical disparity. Our results therefore suggest that the vertical disparity in 3D-3 condition decreases the perception of depth compared to other 3D conditions and the amplitude of P1 component can be used as a discriminative feature.

SIGNIFICANCE

The results show that the P1 component increases in amplitude due to the depth perception in the 3D stimuli compared to the 2D stimulus. On the other hand the vertical disparity in the stereoscopic images is studied here. We suggest that the amplitude of P1 component is modulated with this parameter and decreases due to the decrease in the perception of depth.

Spectral Variability in the Aged Brain during Fine Motor Control

Physiological aging is paralleled by a decline of fine motor skills accompanied by structural and functional alterations of the underlying brain network. Here, we aim to investigate age-related changes in the spectral distribution of neuronal oscillations during fine skilled motor function. We employ the concept of spectral entropy in order to describe the flatness and peaked-ness of a frequency spectrum to quantify changes in the spectral distribution of the oscillatory motor response in the aged brain. Electroencephalogram was recorded in elderly (n = 32) and young (n = 34) participants who performed either a cued finger movement or a pinch or a whole hand grip task with their dominant right hand. Whereas young participant showed distinct, well-defined movement-related power decreases in the alpha and upper beta band, elderly participants exhibited a flat broadband, frequency-unspecific power desynchronization. This broadband response was reflected by an increase of spectral entropy over sensorimotor and frontal areas in the aged brain. Neuronal activation patterns differed between motor tasks in the young brain, while the aged brain showed a similar activation pattern in all tasks. Moreover, we found a wider recruitment of the cortical motor network in the aged brain. The present study adds to the understanding of age-related changes of neural coding during skilled motor behavior, revealing a less predictable signal with great variability across frequencies in a wide cortical motor network in the aged brain. The increase in entropy in the aged brain could be a reflection of random noise-like activity or could represent a compensatory mechanism that serves a functional role.

Bispectral pairwise interacting source analysis for identifying systems of cross-frequency interacting brain sources from electroencephalographic or magnetoencephalographic signals

Brain cognitive functions arise through the coordinated activity of several brain regions, which actually form complex dynamical systems operating at multiple frequencies. These systems often consist of interacting subsystems, whose characterization is of importance for a complete understanding of the brain interaction processes. To address this issue, we present a technique, namely the bispectral pairwise interacting source analysis (biPISA), for analyzing systems of cross-frequency interacting brain sources when multichannel electroencephalographic (EEG) or magnetoencephalographic (MEG) data are available. Specifically, the biPISA makes it possible to identify one or many subsystems of cross-frequency interacting sources by decomposing the antisymmetric components of the cross-bispectra between EEG or MEG signals, based on the assumption that interactions are pairwise. Thanks to the properties of the antisymmetric components of the cross-bispectra, biPISA is also robust to spurious interactions arising from mixing artifacts, i.e., volume conduction or field spread, which always affect EEG or MEG functional connectivity estimates. This method is an extension of the pairwise interacting source analysis (PISA), which was originally introduced for investigating interactions at the same frequency, to the study of cross-frequency interactions. The effectiveness of this approach is demonstrated in simulations for up to three interacting source pairs and for real MEG recordings of spontaneous brain activity. Simulations show that the performances of biPISA in estimating the phase difference between the interacting sources are affected by the increasing level of noise rather than by the number of the interacting subsystems. The analysis of real MEG data reveals an interaction between two pairs of sources of central mu and beta rhythms, localizing in the proximity of the left and right central sulci.

Using reciprocity for relating the simulation of transcranial current stimulation to the EEG forward problem

To explore the relationship between transcranial current stimulation (tCS) and the electroencephalography (EEG) forward problem, we investigate and compare accuracy and efficiency of a reciprocal and a direct EEG forward approach for dipolar primary current sources both based on the finite element method (FEM), namely the adjoint approach (AA) and the partial integration approach in conjunction with a transfer matrix concept (PI). By analyzing numerical results, comparing to analytically derived EEG forward potentials and estimating computational complexity in spherical shell models, AA turns out to be essentially identical to PI. It is then proven that AA and PI are also algebraically identical even for general head models. This relation offers a direct link between the EEG forward problem and tCS. We then demonstrate how the quasi-analytical EEG forward solutions in sphere models can be used to validate the numerical accuracies of FEM-based tCS simulation approaches. These approaches differ with respect to the ease with which they can be employed for realistic head modeling based on MRI-derived segmentations. We show that while the accuracy of the most easy to realize approach based on regular hexahedral elements is already quite high, it can be significantly improved if a geometry-adaptation of the elements is employed in conjunction with an isoparametric FEM approach. While the latter approach does not involve any additional difficulties for the user, it reaches the high accuracies of surface-segmentation based tetrahedral FEM, which is considerably more difficult to implement and topologically less flexible in practice. Finally, in a highly realistic head volume conductor model and when compared to the regular alternative, the geometry-adapted hexahedral FEM is shown to result in significant changes in tCS current flow orientation and magnitude up to 45° and a factor of 1.66, respectively.

Increased Resting-State Gamma-Band Connectivity in First-Episode Schizophrenia

BACKGROUND

Schizophrenia has long been suggested to represent a disorder with prominent neural dysconnectivity. Gamma-band oscillations are highly relevant in this context, due both to their proposed involvement in neuronal synchronization and to their association with neurotransmitter systems relevant for schizophrenia. Several task-related studies have confirmed reduced power and synchronization of gamma-band oscillations in schizophrenia, but it has been suggested that these findings might not apply to the resting state. The present study aimed to investigate resting-state gamma-band connectivity in patients with schizophrenia.

METHODS

Sixty-four channel resting-state electroencephalography (eyes closed) was recorded in 22 patients with first-episode schizophrenia and 22 healthy controls matched for age and gender. Orthogonalized power envelope correlation was used as a measure of connectivity across 80 cortical regions at 40 Hz. Mean connectivity at each region was compared across groups using the nonparametric randomization approach. Additionally, the network-based statistic was applied to identify affected networks in patients.

RESULTS

Patients displayed increased mean functional gamma-band connectivity compared to controls in the left rolandic operculum. Network-based analyses indicated increased connectivity in patients within a strongly lateralized network consisting mainly of left inferior frontal/orbitofrontal, lateral and medial temporal, and inferior parietal areas. Within this network, gamma-band connectivity was higher in patients with low positive and disorganization symptom levels.

CONCLUSIONS

The present study provides a link between resting-state gamma-band connectivity and the core symptoms of schizophrenia. The observed findings are different than those reported by task-related studies, suggesting that resting-state studies might reveal new aspects in the pathophysiology of schizophrenia.

Resting-state theta-band connectivity and verbal memory in schizophrenia and in the high-risk state

BACKGROUND

Disturbed functional connectivity is assumed to underlie neurocognitive deficits in patients with schizophrenia. As neurocognitive deficits are already present in the high-risk state, identification of the neural networks involved in this core feature of schizophrenia is essential to our understanding of the disorder. Resting-state studies enable such investigations, while at the same time avoiding the known confounder of impaired task performance in patients. The aim of the present study was to investigate EEG resting-state connectivity in high-risk individuals (HR) compared to first episode patients with schizophrenia (SZ) and to healthy controls (HC), and its association with cognitive deficits.

METHODS

64-channel resting-state EEG recordings (eyes closed) were obtained for 28 HR, 19 stable SZ, and 23 HC, matched for age, education, and parental education. The imaginary coherence-based multivariate interaction measure (MIM) was used as a measure of connectivity across 80 cortical regions and six frequency bands. Mean connectivity at each region was compared across groups using the non-parametric randomization approach. Additionally, the network-based statistic was applied to identify affected networks in patients.

RESULTS

SZ displayed increased theta-band resting-state MIM connectivity across midline, sensorimotor, orbitofrontal regions and the left temporoparietal junction. HR displayed intermediate theta-band connectivity patterns that did not differ from either SZ or HC. Mean theta-band connectivity within the above network partially mediated verbal memory deficits in SZ and HR.

CONCLUSIONS

Aberrant theta-band connectivity may represent a trait characteristic of schizophrenia associated with neurocognitive deficits. As such, it might constitute a promising target for novel treatment applications.

Selective modulation of interhemispheric functional connectivity by HD-tACS shapes perception

Oscillatory neuronal synchronization between cortical areas has been suggested to constitute a flexible mechanism to coordinate information flow in the human cerebral cortex. However, it remains unclear whether synchronized neuronal activity merely represents an epiphenomenon or whether it is causally involved in the selective gating of information. Here, we combined bilateral high-density transcranial alternating current stimulation (HD-tACS) at 40 Hz with simultaneous electroencephalographic (EEG) recordings to study immediate electrophysiological effects during the selective entrainment of oscillatory gamma-band signatures. We found that interhemispheric functional connectivity was modulated in a predictable, phase-specific way: In-phase stimulation enhanced synchronization, anti-phase stimulation impaired functional coupling. Perceptual correlates of these connectivity changes were found in an ambiguous motion task, which strongly support the functional relevance of long-range neuronal coupling. Additionally, our results revealed a decrease in oscillatory alpha power in response to the entrainment of gamma band signatures. This finding provides causal evidence for the antagonistic role of alpha and gamma oscillations in the parieto-occipital cortex and confirms that the observed gamma band modulations were physiological in nature. Our results demonstrate that synchronized cortical network activity across several spatiotemporal scales is essential for conscious perception and cognition.

Wedge MUSIC: a novel approach to examine experimental differences of brain source connectivity patterns from EEG/MEG data

We introduce a novel method to estimate bivariate synchronization, i.e. interacting brain sources at a specific frequency or band, from MEG or EEG data robust to artifacts of volume conduction. The data driven calculation is solely based on the imaginary part of the cross-spectrum as opposed to the imaginary part of coherency. In principle, the method quantifies how strong a synchronization between a distinct pair of brain sources is present in the data. As an input of the method all pairs of pre-defined locations inside the brain can be used which is computationally exhaustive. In contrast to that, reference sources can be used that have been identified by any source reconstruction technique in a prior analysis step. We introduce different variants of the method and evaluate the performance in simulations. As a particular advantage of the proposed methodology, we demonstrate that the novel approach is capable of investigating differences in brain source interactions between experimental conditions or with respect to a certain baseline. For measured data, we first show the application on resting state MEG data where we find locally synchronized sources in the motor-cortex based on the sensorimotor idle rhythms. Finally, we show an example on EEG motor imagery data where we contrast hand and foot movements. Here, we also find local interactions in the expected brain areas.

Intrinsic coupling modes: multiscale interactions in ongoing brain activity

Intrinsic coupling constitutes a key feature of ongoing brain activity, which exhibits rich spatiotemporal patterning and contains information that influences cognitive processing. We discuss evidence for two distinct types of intrinsic coupling modes which seem to reflect the operation of different coupling mechanisms. One type arises from phase coupling of band-limited oscillatory signals, whereas the other results from coupled aperiodic fluctuations of signal envelopes. The two coupling modes differ in their dynamics, their origins, and their putative functions and with respect to their alteration in neuropsychiatric disorders. We propose that the concept of intrinsic coupling modes can provide a unifying framework for capturing the dynamics of intrinsically generated neuronal interactions at multiple spatial and temporal scales.