Differences in MEG and EEG power-law scaling explained by a coupling between spatial coherence and frequency: a simulation study

Inhibitory neurons and the asymmetric shape of neuronal avalanches

In the last twenty years neuronal avalanches have been deeply investigated, both experimentally and numerically, also framing the results in the context of the avalanche scaling theory. In particular the avalanche shape has recently received a wide attention, also because the existence of a universal shape is an indication of the brain acting at a critical point. Within this scope, the detection of the shape asymmetry and the understanding of the mechanisms leading to it can provide useful insights into brain activity. Experimental data evidence, either symmetric or leftward asymmetry in the shape, results are not confirmed by numerical studies. Here we analyze the role of inhibition, connectivity range, and short term plasticity in determining the avalanche shape in an integrate and fire model. Results indicate that, not only the physiological fraction of inhibitory neurons is crucial to observe leftward asymmetry, but also the different synaptic recovery rates between excitatory and inhibitory neurons, confirming the importance of a dynamic balance between excitation and inhibition in brain activity.

Novel bias-reduced coherence measure for EEG-based speech tracking in listeners with hearing impairment

In the literature, auditory attention is explored through neural speech tracking, primarily entailing modeling and analyzing electroencephalography (EEG) responses to natural speech via linear filtering. Our study takes a novel approach, introducing an enhanced coherence estimation technique to assess the strength of neural speech tracking. This enables effective discrimination between attended and ignored speech. To mitigate the impact of colored noise in EEG, we address two biases-overall coherence-level bias and spectral peak-shifting bias. In a listening study involving 32 participants with hearing impairment, tasked with attending to competing talkers in background noise, our coherence-based method effectively discerns EEG representations of attended and ignored speech. We comprehensively analyze frequency bands, individual frequencies, and EEG channels. Frequency bands of importance are shown to be delta, theta and alpha, and the important EEG channels are the central. Lastly, we showcase coherence differences across different noise reduction settings implemented in hearing aids (HAs), underscoring our method's potential to objectively assess auditory attention and enhance HA efficacy.

Validating MEG source imaging of resting state oscillatory patterns with an intracranial EEG atlas

BACKGROUND

Magnetoencephalography (MEG) is a widely used non-invasive tool to estimate brain activity with high temporal resolution. However, due to the ill-posed nature of the MEG source imaging (MSI) problem, the ability of MSI to identify accurately underlying brain sources along the cortical surface is still uncertain and requires validation.

METHOD

We validated the ability of MSI to estimate the background resting state activity of 45 healthy participants by comparing it to the intracranial EEG (iEEG) atlas (https://mni-open-ieegatlas.

RESEARCH

mcgill.ca/). First, we applied wavelet-based Maximum Entropy on the Mean (wMEM) as an MSI technique. Next, we converted MEG source maps into intracranial space by applying a forward model to the MEG-reconstructed source maps, and estimated virtual iEEG (ViEEG) potentials on each iEEG channel location; we finally quantitatively compared those with actual iEEG signals from the atlas for 38 regions of interest in the canonical frequency bands.

RESULTS

The MEG spectra were more accurately estimated in the lateral regions compared to the medial regions. The regions with higher amplitude in the ViEEG than in the iEEG were more accurately recovered. In the deep regions, MEG-estimated amplitudes were largely underestimated and the spectra were poorly recovered. Overall, our wMEM results were similar to those obtained with minimum norm or beamformer source localization. Moreover, the MEG largely overestimated oscillatory peaks in the alpha band, especially in the anterior and deep regions. This is possibly due to higher phase synchronization of alpha oscillations over extended regions, exceeding the spatial sensitivity of iEEG but detected by MEG. Importantly, we found that MEG-estimated spectra were more comparable to spectra from the iEEG atlas after the aperiodic components were removed.

CONCLUSION

This study identifies brain regions and frequencies for which MEG source analysis is likely to be reliable, a promising step towards resolving the uncertainty in recovering intracerebral activity from non-invasive MEG studies.

Modulation in alpha band activity reflects syntax composition: an MEG study of minimal syntactic binding

Successful sentence comprehension requires the binding, or composition, of multiple words into larger structures to establish meaning. Using magnetoencephalography, we investigated the neural mechanisms involved in binding at the syntax level, in a task where contributions from semantics were minimized. Participants were auditorily presented with minimal sentences that required binding (pronoun and pseudo-verb with the corresponding morphological inflection; "she grushes") and pseudo-verb wordlists that did not require binding ("cugged grushes"). Relative to no binding, we found that syntactic binding was associated with a modulation in alpha band (8-12 Hz) activity in left-lateralized language regions. First, we observed a significantly smaller increase in alpha power around the presentation of the target word ("grushes") that required binding (-0.05 to 0.1 s), which we suggest reflects an expectation of binding to occur. Second, during binding of the target word (0.15-0.25 s), we observed significantly decreased alpha phase-locking between the left inferior frontal gyrus and the left middle/inferior temporal cortex, which we suggest reflects alpha-driven cortical disinhibition serving to strengthen communication within the syntax composition neural network. Altogether, our findings highlight the critical role of rapid spatial-temporal alpha band activity in controlling the allocation, transfer, and coordination of the brain's resources during syntax composition.

Scale-free dynamics of core-periphery topography

The human brain's cerebral cortex exhibits a topographic division into higher-order transmodal core and lower-order unimodal periphery regions. While timescales between the core and periphery region diverge, features of their power spectra, especially scale-free dynamics during resting-state and their mdulation in task states, remain unclear. To answer this question, we investigated the ~1/f-like pink noise manifestation of scale-free dynamics in the core-periphery topography during rest and task states applying infra-slow inter-trial intervals up to 1 min falling inside the BOLD's infra-slow frequency band. The results demonstrate (1) higher resting-state power-law exponent (PLE) in the core compared to the periphery region; (2) significant PLE increases in task across the core and periphery regions; and (3) task-related PLE increases likely followed the task's atypically low event rates, namely the task's periodicity (inter-trial interval = 52-60 s; 0.016-0.019 Hz). A computational model and a replication dataset that used similar infra-slow inter-trial intervals provide further support for our main findings. Altogether, the results show that scale-free dynamics differentiate core and periphery regions in the resting-state and mediate task-related effects.

Scale-free and oscillatory spectral measures of sleep stages in humans

Power spectra of sleep electroencephalograms (EEG) comprise two main components: a decaying power-law corresponding to the aperiodic neural background activity, and spectral peaks present due to neural oscillations. “Traditional” band-based spectral methods ignore this fundamental structure of the EEG spectra and thus are susceptible to misrepresenting the underlying phenomena. A fitting method that attempts to separate and parameterize the aperiodic and periodic spectral components called “fitting oscillations and one over f” (FOOOF) was applied to a set of annotated whole-night sleep EEG recordings of 251 subjects from a wide age range (4–69 years). Most of the extracted parameters exhibited sleep stage sensitivity; significant main effects and interactions of sleep stage, age, sex, and brain region were found. The spectral slope (describing the steepness of the aperiodic component) showed especially large and consistent variability between sleep stages (and low variability between subjects), making it a candidate indicator of sleep states. The limitations and arisen problems of the FOOOF method are also discussed, possible solutions for some of them are suggested.

Scaling of avalanche shape and activity power spectrum in neuronal networks

Many systems in nature exhibit avalanche dynamics with scale-free features. A general scaling theory has been proposed for critical avalanche profiles in crackling noise, predicting the collapse onto a universal avalanche shape, as well as the scaling behavior of the activity power spectrum as Brown noise. Recently, much attention has been given to the profile of neuronal avalanches, measured in neuronal systems in vitro and in vivo. Although a universal profile was evidenced, confirming the validity of the general scaling theory, the parallel study of the power spectrum scaling under the same conditions was not performed. The puzzling observation is that in the majority of healthy neuronal systems the power spectrum exhibits a behavior close to 1/f, rather than Brown, noise. Here we perform a numerical study of the scaling behavior of the avalanche shape and the power spectrum for a model of integrate and fire neurons with a short-term plasticity parameter able to tune the system to criticality. We confirm that, at criticality, the average avalanche size and the avalanche profile fulfill the general avalanche scaling theory. However, the power spectrum consistently exhibits Brown noise behavior, for both fully excitatory networks and systems with 30% inhibitory networks. Conversely, a behavior closer to 1/f noise is observed in systems slightly off criticality. Results suggest that the power spectrum is a good indicator to determine how close neuronal activity is to criticality.

Simulating epileptic seizures using the bidomain model

Epileptic seizures are due to excessive and synchronous neural activity. Extensive modelling of seizures has been done on the neuronal level, but it remains a challenge to scale these models up to whole brain models. Measurements of the brain's activity over several spatiotemporal scales follow a power-law distribution in terms of frequency. During normal brain activity, the power-law exponent is often found to be around 2 for frequencies between a few Hz and up to 150 Hz, but is higher during seizures and for higher frequencies. The Bidomain model has been used with success in modelling the electrical activity of the heart, but has been explored far less in the context of the brain. This study extends previous models of epileptic seizures on the neuronal level to the whole brain using the Bidomain model. Our approach is evaluated in terms of power-law distributions. The electric potentials were simulated in 7 idealized two-dimensional models and 3 three-dimensional patient-specific models derived from magnetic resonance images (MRI). Computed electric potentials were found to follow power-law distributions with slopes ranging from 2 to 5 for frequencies greater than 10-30 Hz.

Critical behaviour of the stochastic Wilson-Cowan model

Spontaneous brain activity is characterized by bursts and avalanche-like dynamics, with scale-free features typical of critical behaviour. The stochastic version of the celebrated Wilson-Cowan model has been widely studied as a system of spiking neurons reproducing non-trivial features of the neural activity, from avalanche dynamics to oscillatory behaviours. However, to what extent such phenomena are related to the presence of a genuine critical point remains elusive. Here we address this central issue, providing analytical results in the linear approximation and extensive numerical analysis. In particular, we present results supporting the existence of a bona fide critical point, where a second-order-like phase transition occurs, characterized by scale-free avalanche dynamics, scaling with the system size and a diverging relaxation time-scale. Moreover, our study shows that the observed critical behaviour falls within the universality class of the mean-field branching process, where the exponents of the avalanche size and duration distributions are, respectively, 3/2 and 2. We also provide an accurate analysis of the system behaviour as a function of the total number of neurons, focusing on the time correlation functions of the firing rate in a wide range of the parameter space.