Measuring electrophysiological connectivity by power envelope correlation: a technical review on MEG methods

Genetic fingerprinting with heritable phenotypes of the resting-state brain network topology

Cognitive, behavioral, and disease traits are influenced by both genetic and environmental factors. Individual differences in these traits have been associated with graph theoretical properties of resting-state networks, indicating that variations in connectome topology may be driven by genetics. In this study, we establish the heritability of global and local graph properties of resting-state networks derived from functional MRI (fMRI) and magnetoencephalography (MEG) using a large sample of twins and non-twin siblings from the Human Connectome Project. We examine the heritability of MEG in the source space, providing a more accurate estimate of genetic influences on electrophysiological networks. Our findings show that most graph measures are more heritable for MEG compared to fMRI and the heritability for MEG is greater for amplitude compared to phase synchrony in the delta, high beta, and gamma frequency bands. This suggests that the fast neuronal dynamics in MEG offer unique insights into the genetic basis of brain network organization. Furthermore, we demonstrate that brain network features can serve as genetic fingerprints to accurately identify pairs of identical twins within a cohort. These results highlight novel opportunities to relate individual connectome signatures to genetic mechanisms underlying brain function.

Altered long-range functional connectivity in PTSD: Role of the infraslow oscillations of cortical activity amplitude envelopes

OBJECTIVE

Coupling between the amplitude envelopes (AEs) of regional cortical activity reflects mechanisms that coordinate the excitability of large-scale cortical networks. We used resting-state MEG recordings to investigate the association between alterations in the coupling of cortical AEs and symptoms of post-traumatic stress disorder (PTSD).

METHODS

Participants (n = 96) were service members with combat exposure and various levels of post-traumatic stress severity (PTSS). We assessed the correlation between PTSS and (1) coupling of broadband cortical AEs of beta band activity, (2) coupling of the low- (<0.5 Hz) and high-frequency (>0.5 Hz) components of the AEs, and (3) their time-varying patterns.

RESULTS

PTSS was associated with widespread hypoconnectivity assessed from the broadband AE fluctuations, which correlated with subscores for the negative thoughts and feelings/emotional numbing (NTF/EN) and hyperarousal clusters of symptoms. Higher NTF/EN scores were also associated with smaller increases in resting-state functional connectivity (rsFC) with time during the recordings. The distinct patterns of rsFC in PTSD were primarily due to differences in the coupling of low-frequency (infraslow) fluctuations of the AEs of beta band activity.

CONCLUSIONS

Our findings implicate the mechanisms underlying the regulation/coupling of infraslow oscillations in the alterations of rsFC assessed from broadband AEs and in PTSD symptomatology.

SIGNIFICANCE

Altered coordination of infraslow amplitude fluctuations across large-scale cortical networks can contribute to network dysfunction and may provide a target for treatment in PTSD.

Automatic detection of ALS from single-trial MEG signals during speech tasks: a pilot study

Amyotrophic lateral sclerosis (ALS) is an idiopathic, fatal, and fast-progressive neurodegenerative disease characterized by the degeneration of motor neurons. ALS patients often experience an initial misdiagnosis or a diagnostic delay due to the current unavailability of an efficient biomarker. Since impaired speech is typical in ALS, we hypothesized that functional differences between healthy and ALS participants during speech tasks can be explained by cortical pattern changes, thereby leading to the identification of a neural biomarker for ALS. In this pilot study, we collected magnetoencephalography (MEG) recordings from three early-diagnosed patients with ALS and three healthy controls during imagined (covert) and overt speech tasks. First, we computed sensor correlations, which showed greater correlations for speakers with ALS than healthy controls. Second, we compared the power of the MEG signals in canonical bands between the two groups, which showed greater dissimilarity in the beta band for ALS participants. Third, we assessed differences in functional connectivity, which showed greater beta band connectivity for ALS than healthy controls. Finally, we performed single-trial classification, which resulted in highest performance with beta band features (∼ 98%). These findings were consistent across trials, phrases, and participants for both imagined and overt speech tasks. Our preliminary results indicate that speech-evoked beta oscillations could be a potential neural biomarker for diagnosing ALS. To our knowledge, this is the first demonstration of the detection of ALS from single-trial neural signals.

Tracking the neurodevelopmental trajectory of beta band oscillations with optically pumped magnetometer-based magnetoencephalography

Neural oscillations mediate the coordination of activity within and between brain networks, supporting cognition and behaviour. How these processes develop throughout childhood is not only an important neuroscientific question but could also shed light on the mechanisms underlying neurological and psychiatric disorders. However, measuring the neurodevelopmental trajectory of oscillations has been hampered by confounds from instrumentation. In this paper, we investigate the suitability of a disruptive new imaging platform - optically pumped magnetometer-based magnetoencephalography (OPM-MEG) - to study oscillations during brain development. We show how a unique 192-channel OPM-MEG device, which is adaptable to head size and robust to participant movement, can be used to collect high-fidelity electrophysiological data in individuals aged between 2 and 34 years. Data were collected during a somatosensory task, and we measured both stimulus-induced modulation of beta oscillations in sensory cortex, and whole-brain connectivity, showing that both modulate significantly with age. Moreover, we show that pan-spectral bursts of electrophysiological activity drive task-induced beta modulation, and that their probability of occurrence and spectral content change with age. Our results offer new insights into the developmental trajectory of beta oscillations and provide clear evidence that OPM-MEG is an ideal platform for studying electrophysiology in neurodevelopment.

The effects of spatial leakage correction on the reliability of EEG-based functional connectivity networks

Electroencephalography (EEG) functional connectivity (FC) estimates are confounded by the volume conduction problem. This effect can be greatly reduced by applying FC measures insensitive to instantaneous, zero-lag dependencies (corrected measures). However, numerous studies showed that FC measures sensitive to volume conduction (uncorrected measures) exhibit higher reliability and higher subject-level identifiability. We tested how source reconstruction contributed to the reliability difference of EEG FC measures on a large (n = 201) resting-state data set testing eight FC measures (including corrected and uncorrected measures). We showed that the high reliability of uncorrected FC measures in resting state partly stems from source reconstruction: idiosyncratic noise patterns define a baseline resting-state functional network that explains a significant portion of the reliability of uncorrected FC measures. This effect remained valid for template head model-based, as well as individual head model-based source reconstruction. Based on our findings we made suggestions how to best use spatial leakage corrected and uncorrected FC measures depending on the main goals of the study.

Dynamic functional connectivity is modulated by the amount of p-Tau231 in blood in cognitively intact participants

INTRODUCTION

Electrophysiology and plasma biomarkers are early and non-invasive candidates for Alzheimer's disease detection. The purpose of this paper is to evaluate changes in dynamic functional connectivity measured with magnetoencephalography, associated with the plasma pathology marker p-tau231 in unimpaired adults.

METHODS

73 individuals were included. Static and dynamic functional connectivity were calculated using leakage corrected amplitude envelope correlation. Each source's strength entropy across trials was calculated. A data-driven statistical analysis was performed to find the association between functional connectivity and plasma p-tau231 levels. Regression models were used to assess the influence of other variables over the clusters' connectivity.

RESULTS

Frontotemporal dynamic connectivity positively associated with p-tau231 levels. Linear regression models identified pathological, functional and structural factors that influence dynamic functional connectivity.

DISCUSSION

These results expand previous literature on dynamic functional connectivity in healthy individuals at risk of AD, highlighting its usefulness as an early, non-invasive and more sensitive biomarker.

Bayesian reduced rank regression models generalizable neural fingerprints that differentiate between individuals in magnetoencephalography data

Recent magnetoencephalography (MEG) studies have reported that functional connectivity (FC) and power spectra can be used as neural fingerprints in differentiating individuals. Such studies have mainly used correlations between measurement sessions to distinguish individuals from each other. However, it has remained unclear whether such correlations might reflect a more generalizable principle of individually distinctive brain patterns. Here, we evaluated a machine-learning based approach, termed latent-noise Bayesian reduced rank regression (BRRR) as a means of modelling individual differences in the resting-state MEG data of the Human Connectome Project (HCP), using FC and power spectra as neural features. First, we verified that BRRR could model and reproduce the differences between metrics that correlation-based fingerprinting yields. We trained BRRR models to distinguish individuals based on data from one measurement and used the models to identify subsequent measurement sessions of those same individuals. The best performing BRRR models, using only 20 spatiospectral components, were able to identify subjects across measurement sessions with over 90% accuracy, approaching the highest correlation-based accuracies. Using cross-validation, we then determined whether that BRRR model could generalize to unseen subjects, successfully classifying the measurement sessions of novel individuals with over 80% accuracy. The results demonstrate that individual neurofunctional differences can be reliably extracted from MEG data with a low-dimensional predictive model and that the model is able to classify novel subjects.

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

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

The neurodevelopmental trajectory of beta band oscillations: an OPM-MEG study

Neural oscillations mediate the coordination of activity within and between brain networks, supporting cognition and behaviour. How these processes develop throughout childhood is not only an important neuroscientific question but could also shed light on the mechanisms underlying neurological and psychiatric disorders. However, measuring the neurodevelopmental trajectory of oscillations has been hampered by confounds from instrumentation. In this paper, we investigate the suitability of a disruptive new imaging platform - Optically Pumped Magnetometer-based magnetoencephalography (OPM-MEG) - to study oscillations during brain development. We show how a unique 192-channel OPM-MEG device, which is adaptable to head size and robust to participant movement, can be used to collect high-fidelity electrophysiological data in individuals aged between 2 and 34 years. Data were collected during a somatosensory task, and we measured both stimulus-induced modulation of beta oscillations in sensory cortex, and whole-brain connectivity, showing that both modulate significantly with age. Moreover, we show that pan-spectral bursts of electrophysiological activity drive task-induced beta modulation, and that their probability of occurrence and spectral content change with age. Our results offer new insights into the developmental trajectory of beta oscillations and provide clear evidence that OPM-MEG is an ideal platform for studying electrophysiology in neurodevelopment.

Network Hyperexcitability in Early-Stage Alzheimer's Disease: Evaluation of Functional Connectivity Biomarkers in a Computational Disease Model

Background

There is increasing evidence from animal and clinical studies that network hyperexcitability (NH) may be an important pathophysiological process and potential target for treatment in early Alzheimer's disease (AD). Measures of functional connectivity (FC) have been proposed as promising biomarkers for NH, but it is unknown which measure has the highest sensitivity for early-stage changes in the excitation/inhibition balance.

Objective

We aim to test the performance of different FC measures in detecting NH at the earliest stage using a computational approach.

Methods

We use a whole brain computational model of activity dependent degeneration to simulate progressive AD pathology and NH. We investigate if and at what stage four measures of FC (amplitude envelope correlation corrected [AECc], phase lag index [PLI], joint permutation entropy [JPE] and a new measure: phase lag time [PLT]) can detect early-stage AD pathophysiology.

Results

The activity dependent degeneration model replicates spectral changes in line with clinical data and demonstrates increasing NH. Compared to relative theta power as a gold standard the AECc and PLI are shown to be less sensitive in detecting early-stage NH and AD-related neurophysiological abnormalities, while the JPE and the PLT show more sensitivity with excellent test characteristics.

Conclusions

Novel FC measures, which are better in detecting rapid fluctuations in neural activity and connectivity, may be superior to well-known measures such as the AECc and PLI in detecting early phase neurophysiological abnormalities and in particular NH in AD. These markers could improve early diagnosis and treatment target identification.

Tuning Minimum-Norm regularization parameters for optimal MEG connectivity estimation

The accurate characterization of cortical functional connectivity from Magnetoencephalography (MEG) data remains a challenging problem due to the subjective nature of the analysis, which requires several decisions at each step of the analysis pipeline, such as the choice of a source estimation algorithm, a connectivity metric and a cortical parcellation, to name but a few. Recent studies have emphasized the importance of selecting the regularization parameter in minimum norm estimates with caution, as variations in its value can result in significant differences in connectivity estimates. In particular, the amount of regularization that is optimal for MEG source estimation can actually be suboptimal for coherence-based MEG connectivity analysis. In this study, we expand upon previous work by examining a broader range of commonly used connectivity metrics, including the imaginary part of coherence, corrected imaginary part of Phase Locking Value, and weighted Phase Lag Index, within a larger and more realistic simulation scenario. Our results show that the best estimate of connectivity is achieved using a regularization parameter that is 1 or 2 orders of magnitude smaller than the one that yields the best source estimation. This remarkable difference may imply that previous work assessing source-space connectivity using minimum-norm may have benefited from using less regularization, as this may have helped reduce false positives. Importantly, we provide the code for MEG data simulation and analysis, offering the research community a valuable open source tool for informed selections of the regularization parameter when using minimum-norm for source space connectivity analyses.

Texture-based brain networks for characterization of healthy subjects from MRI

Brain networks have been widely used to study the relationships between brain regions based on their dynamics using, e.g. fMRI or EEG, and to characterize their real physical connections using DTI. However, few studies have investigated brain networks derived from structural properties; and those have been based on cortical thickness or gray matter volume. The main objective of this work was to investigate the feasibility of obtaining useful information from brain networks derived from structural MRI, using texture features. We also wanted to verify if texture brain networks had any relation with established functional networks. T1-MR images were segmented using AAL and texture parameters from the gray-level co-occurrence matrix were computed for each region, for 760 subjects. Individual texture networks were used to evaluate the structural connections between regions of well-established functional networks; assess possible gender differences; investigate the dependence of texture network measures with age; and single out brain regions with different texture-network characteristics. Although around 70% of texture connections between regions belonging to the default mode, attention, and visual network were greater than the mean connection value, this effect was small (only between 7 and 15% of these connections were larger than one standard deviation), implying that texture-based morphology does not seem to subside function. This differs from cortical thickness-based morphology, which has been shown to relate to functional networks. Seventy-five out of 86 evaluated regions showed significant (ANCOVA, p < 0.05) differences between genders. Forty-four out of 86 regions showed significant (ANCOVA, p < 0.05) dependence with age; however, the R2 indicates that this is not a linear relation. Thalamus and putamen showed a very unique texture-wise structure compared to other analyzed regions. Texture networks were able to provide useful information regarding gender and age-related differences, as well as for singling out specific brain regions. We did not find a morphological texture-based subsidy for the evaluated functional brain networks. In the future, this approach will be extended to neurological patients to investigate the possibility of extracting biomarkers to help monitor disease evolution or treatment effectiveness.

Topological data analysis for revealing dynamic brain reconfiguration in MEG data

In recent years, the focus of the functional connectivity community has shifted from stationary approaches to the ones that include temporal dynamics. Especially, non-invasive electrophysiological data (magnetoencephalography/electroencephalography (MEG/EEG)) with high temporal resolution and good spatial coverage have made it possible to measure the fast alterations in the neural activity in the brain during ongoing cognition. In this article, we analyze dynamic brain reconfiguration using MEG images collected from subjects during the rest and the cognitive tasks. Our proposed topological data analysis method, called Mapper, produces biomarkers that differentiate cognitive tasks without prior spatial and temporal collapse of the data. The suggested method provides an interactive visualization of the rapid fluctuations in electrophysiological data during motor and cognitive tasks; hence, it has the potential to extract clinically relevant information at an individual level without temporal and spatial collapse.

Network Hyperexcitability in Early Alzheimer's Disease: Is Functional Connectivity a Potential Biomarker?

Network hyperexcitability (NH) is an important feature of the pathophysiology of Alzheimer's disease. Functional connectivity (FC) of brain networks has been proposed as a potential biomarker for NH. Here we use a whole brain computational model and resting-state MEG recordings to investigate the relation between hyperexcitability and FC. Oscillatory brain activity was simulated with a Stuart Landau model on a network of 78 interconnected brain regions. FC was quantified with amplitude envelope correlation (AEC) and phase coherence (PC). MEG was recorded in 18 subjects with subjective cognitive decline (SCD) and 18 subjects with mild cognitive impairment (MCI). Functional connectivity was determined with the corrected AECc and phase lag index (PLI), in the 4-8 Hz and the 8-13 Hz bands. The excitation/inhibition balance in the model had a strong effect on both AEC and PC. This effect was different for AEC and PC, and was influenced by structural coupling strength and frequency band. Empirical FC matrices of SCD and MCI showed a good correlation with model FC for AEC, but less so for PC. For AEC the fit was best in the hyperexcitable range. We conclude that FC is sensitive to changes in E/I balance. The AEC was more sensitive than the PLI, and results were better for the thetaband than the alpha band. This conclusion was supported by fitting the model to empirical data. Our study justifies the use of functional connectivity measures as surrogate markers for E/I balance.

Brain-wide network analysis of resting-state neuromagnetic data

The present study performed a brain-wide network analysis of resting-state magnetoencephalograms recorded from 53 healthy participants to visualize elaborate brain maps of phase- and amplitude-derived graph-theory metrics at different frequencies. To achieve this, we conducted a vertex-wise computation of threshold-independent graph metrics by combining proportional thresholding and a conjunction analysis and applied them to a correlation analysis of age and brain networks. Source power showed a frequency-dependent cortical distribution. Threshold-independent graph metrics derived from phase- and amplitude-based connectivity showed similar or different distributions depending on frequency. Vertex-wise age-brain correlation maps revealed that source power at the beta band and the amplitude-based degree at the alpha band changed with age in local regions. The present results indicate that a brain-wide analysis of neuromagnetic data has the potential to reveal neurophysiological network features in the human brain in a resting state.

Dissociation between phase and power correlation networks in the human brain is driven by co-occurrent bursts

Well-known haemodynamic resting-state networks are better mirrored in power correlation networks than phase coupling networks in electrophysiological data. However, what do these power correlation networks reflect? We address this long-outstanding question in neuroscience using rigorous mathematical analysis, biophysical simulations with ground truth and application of these mathematical concepts to empirical magnetoencephalography (MEG) data. Our mathematical derivations show that for two non-Gaussian electrophysiological signals, their power correlation depends on their coherence, cokurtosis and conjugate-coherence. Only coherence and cokurtosis contribute to power correlation networks in MEG data, but cokurtosis is less affected by artefactual signal leakage and better mirrors haemodynamic resting-state networks. Simulations and MEG data show that cokurtosis may reflect co-occurrent bursting events. Our findings shed light on the origin of the complementary nature of power correlation networks to phase coupling networks and suggests that the origin of resting-state networks is partly reflected in co-occurent bursts in neuronal activity.

MEG Oscillatory Slowing in Cognitive Impairment is Associated with the Presence of Subjective Cognitive Decline

The mechanisms behind Alzheimer's disease are not yet fully described, and changes in the electrophysiology of patients across the continuum of the disease could help to understand them. In this work, we study the power spectral distribution of a set of 129 individuals from the Connectomics of Brian Aging and Dementia project.From this sample, we acquired task-free data, with eyes closed, and estimated the power spectral distribution in source space. We compared the spectral profiles of three groups of individuals: 70 healthy controls, 27 patients with amnestic MCI, and 32 individuals showing cognitive impairment without subjective complaints (IWOC).The results showed a slowing of the brain activity in the aMCI patients, when compared to both the healthy controls and the IWOC individuals. These differences appeared both as a decrease in power for high frequency oscillations and an increase in power in alpha oscillations. The slowing of the spectrum was significant mainly in parietal and medial frontal areas.We were able to validate the slowing of the brain activity in individuals with aMCI, appearing in our sample in areas related to the default mode network. However, this pattern did not appear in the IWOC individuals, suggesting that their condition is not part of the AD continuum. This work raises interesting questions about this group of individuals, and the underlying brain mechanisms behind their cognitive impairment.

Modulations of local synchrony over time lead to resting-state functional connectivity in a parsimonious large-scale brain model

Biophysical models of large-scale brain activity are a fundamental tool for understanding the mechanisms underlying the patterns observed with neuroimaging. These models combine a macroscopic description of the within- and between-ensemble dynamics of neurons within a single architecture. A challenge for these models is accounting for modulations of within-ensemble synchrony over time. Such modulations in local synchrony are fundamental for modeling behavioral tasks and resting-state activity. Another challenge comes from the difficulty in parametrizing large scale brain models which hinders researching principles related with between-ensembles differences. Here we derive a parsimonious large scale brain model that can describe fluctuations of local synchrony. Crucially, we do not reduce within-ensemble dynamics to macroscopic variables first, instead we consider within and between-ensemble interactions similarly while preserving their physiological differences. The dynamics of within-ensemble synchrony can be tuned with a parameter which manipulates local connectivity strength. We simulated resting-state static and time-resolved functional connectivity of alpha band envelopes in models with identical and dissimilar local connectivities. We show that functional connectivity emerges when there are high fluctuations of local and global synchrony simultaneously (i.e. metastable dynamics). We also show that for most ensembles, leaning towards local asynchrony or synchrony correlates with the functional connectivity with other ensembles, with the exception of some regions belonging to the default-mode network.

Disrupted local beta band networks in schizophrenia revealed through graph analysis: A magnetoencephalography study

AIMS

Schizophrenia (SZ) is characterized by psychotic symptoms and cognitive impairment, and is hypothesized to be a 'dysconnection' syndrome due to abnormal neural network formation. Although numerous studies have helped elucidate the pathophysiology of SZ, many aspects of the mechanism underlying psychotic symptoms remain unknown. This study used graph theory analysis to evaluate the characteristics of the resting-state network (RSN) in terms of microscale and macroscale indices, and to identify candidates as potential biomarkers of SZ. Specifically, we discriminated topological characteristics in the frequency domain and investigated them in the context of psychotic symptoms in patients with SZ.

METHODS

We performed graph theory analysis of electrophysiological RSN data using magnetoencephalography to compare topological characteristics represented by microscale (degree centrality and clustering coefficient) and macroscale (global efficiency, local efficiency, and small-worldness) indices in 29 patients with SZ and 38 healthy controls. In addition, we investigated the aberrant topological characteristics of the RSN in patients with SZ and their relationship with SZ symptoms.

RESULTS

SZ was associated with a decreased clustering coefficient, local efficiency, and small-worldness, especially in the high beta band. In addition, macroscale changes in the low beta band are closely associated with negative symptoms.

CONCLUSIONS

The local networks of patients with SZ may disintegrate at both the microscale and macroscale levels, mainly in the beta band. Adopting an electrophysiological perspective of SZ as a failure to form local networks in the beta band will provide deeper insights into the pathophysiology of SZ as a 'dysconnection' syndrome.

Brain connectivity changes of propofol-induced altered states of consciousness using High-Density EEG Source Estimation

Through source estimation, high-density electroencephalogram (EEG) signals at scalp level can be converted into signals at cerebral cortex level, which helps to measure cortical activity during anesthesia induced changes in consciousness level to explore the mechanism. In this research, the high-density EEG of propofol-induced consciousness states alterations in 20 healthy adults were converted into cortical signals of 68 regions of interest (ROI), after alpha bandpass filtering, the pairwise orthogonal power envelope connectivity (PEC) was calculated. Then, due to the number of PECs was huge, the least absolute shrinkage and selection operator (LASSO) was used to select as few PECs as possible as the indicators to distinguish baseline (BS) and moderate sedation (MD) states. The results show that most PECs that can be used as indicators are related to ROI related to default mode network (DMN). At the same time, changes of thalamocortical connectivity and frontal-parietal connectivity could be observed, similar to the neuroimaging method of directly measuring cerebral cortical activity. By extracting the PEC as a classifier to classify the BS and MD States, the accuracy could reach more than 70%. Therefore, this method can not only reflect the mechanism of cortical activity alterations induced by anesthetics, but also provide a new idea for monitoring the depth of anesthesia in the future. Clinical Relevance - This shows that the high-density EEG of scalp level can be converted into cortical signals by source estimation, which is similar to the neuroimaging method of directly measuring cortical activity.

Electrophysiological connectivity markers of preserved language functions in post-stroke aphasia

Post-stroke aphasia is a consequence of localized stroke-related damage as well as global disturbances in a highly interactive and bilaterally-distributed language network. Aphasia is increasingly accepted as a network disorder and it should be treated as such when examining the reorganization and recovery mechanisms after stroke. In the current study, we sought to investigate reorganized patterns of electrophysiological connectivity, derived from resting-state magnetoencephalography (rsMEG), in post-stroke chronic (>6 months after onset) aphasia. We implemented amplitude envelope correlations (AEC), a metric of connectivity commonly used to describe slower aspects of interregional communication in resting-state electrophysiological data. The main focus was on identifying the oscillatory frequency bands and frequency-specific spatial topology of connections associated with preserved language abilities after stroke.

RsMEG was recorded for 5 min in 21 chronic stroke survivors with aphasia and in 20 matched healthy controls. Source-level MEG activity was reconstructed and summarized within 72 atlas-defined brain regions (or nodes). A 72 × 72 leakage-corrected connectivity (of AEC) matrix was obtained for frequencies from theta to low-gamma (4–50 Hz). Connectivity was compared between groups, and, the correlations between connectivity and subscale scores from the Western Aphasia Battery (WAB) were evaluated in the stroke group, using partial least squares analyses. Posthoc multiple regression analyses were also conducted on a graph theory measure of node strengths, derived from significant connectivity results, to control for node-wise properties (local spectral power and lesion sizes) and demographic and stroke-related variables.

Connectivity among the left hemisphere regions, i.e. those ipsilateral to the stroke lesion, was greatly reduced in stroke survivors with aphasia compared to matched healthy controls in the alpha (8–13 Hz; p = 0.011) and beta (15–30 Hz; p = 0.001) bands. The spatial topology of hypoconnectivity in the alpha vs. beta bands was distinct, revealing a greater involvement of ventral frontal, temporal and parietal areas in alpha, and dorsal frontal and parietal areas in beta. The node strengths from alpha and beta group differences remained significant after controlling for nodal spectral power. AEC correlations with WAB subscales of object naming and fluency were significant. Greater alpha connectivity was associated with better naming performance (p = 0.045), and greater connectivity in both the alpha (p = 0.033) and beta (p = 0.007) bands was associated with better speech fluency performance. The spatial topology was distinct between these frequency bands. The node strengths remained significant after controlling for age, time post stroke onset, nodal spectral power and nodal lesion sizes.

Our findings provide important insights into the electrophysiological connectivity profiles (frequency and spatial topology) potentially underpinning preserved language abilities in stroke survivors with aphasia.

Interference suppression techniques for OPM-based MEG: Opportunities and challenges

One of the primary technical challenges facing magnetoencephalography (MEG) is that the magnitude of neuromagnetic fields is several orders of magnitude lower than interfering signals. Recently, a new type of sensor has been developed - the optically pumped magnetometer (OPM). These sensors can be placed directly on the scalp and move with the head during participant movement, making them wearable. This opens up a range of exciting experimental and clinical opportunities for OPM-based MEG experiments, including paediatric studies, and the incorporation of naturalistic movements into neuroimaging paradigms. However, OPMs face some unique challenges in terms of interference suppression, especially in situations involving mobile participants, and when OPMs are integrated with electrical equipment required for naturalistic paradigms, such as motion capture systems. Here we briefly review various hardware solutions for OPM interference suppression. We then outline several signal processing strategies aimed at increasing the signal from neuromagnetic sources. These include regression-based strategies, temporal filtering and spatial filtering approaches. The focus is on the practical application of these signal processing algorithms to OPM data. In a similar vein, we include two worked-through experiments using OPM data collected from a whole-head sensor array. These tutorial-style examples illustrate how the steps for suppressing external interference can be implemented, including the associated data and code so that researchers can try the pipelines for themselves. With the popularity of OPM-based MEG rising, there will be an increasing need to deal with interference suppression. We hope this practical paper provides a resource for OPM-based MEG researchers to build upon.

Looking through the windows: a study about the dependency of phase-coupling estimates on the data length

OBJECTIVE

Being able to characterize functional connectivity (FC) state dynamics in a real-time setting, such as in brain-computer interface, neurofeedback or closed-loop neurostimulation frameworks, requires the rapid detection of the statistical dependencies that quantify FC in short windows of data. The aim of this study is to characterize, through extensive realistic simulations, the reliability of FC estimation as a function of the data length. In particular, we focused on FC as measured by phase-coupling (PC) of neuronal oscillations, one of the most functionally relevant neural coupling modes.

APPROACH

We generated synthetic data corresponding to different scenarios by varying the data length, the signal-to-noise ratio, the phase difference value, the spectral analysis approach (Hilbert or Fourier) and the fractional bandwidth. We compared seven PC metrics, i.e. imaginary part of phase locking value (PLV), PLV of orthogonalized signals, phase lag index (PLI), debiased weighted PLI, imaginary part of coherency, coherence of orthogonalized signals and lagged coherence.

MAIN RESULTS

Our findings show that, for a signal-to-noise-ratio of at least 10 dB, a data window that contains 5 to 8 cycles of the oscillation of interest (e.g. a 500-800ms window at 10Hz) is generally required to achieve reliable PC estimates. In general, Hilbert-based approaches were associated with higher performance than Fourier-based approaches. Furthermore, the results suggest that, when the analysis is performed in a narrow frequency range, a larger window is required.

SIGNIFICANCE

The achieved results pave the way to the introduction of best-practice guidelines to be followed when a real-time frequency-specific PC assessment is at target.

MVPAlab: A machine learning decoding toolbox for multidimensional electroencephalography data

BACKGROUND AND OBJECTIVE

The study of brain function has recently expanded from classical univariate to multivariate analyses. These multivariate, machine learning-based algorithms afford neuroscientists extracting more detailed and richer information from the data. However, the implementation of these procedures is usually challenging, especially for researchers with no coding experience. To address this problem, we have developed MVPAlab, a MATLAB-based, flexible decoding toolbox for multidimensional electroencephalography and magnetoencephalography data.

METHODS

The MVPAlab Toolbox implements several machine learning algorithms to compute multivariate pattern analyses, cross-classification, temporal generalization matrices and feature and frequency contribution analyses. It also provides access to an extensive set of preprocessing routines for, among others, data normalization, data smoothing, dimensionality reduction and supertrial generation. To draw statistical inferences at the group level, MVPAlab includes a non-parametric cluster-based permutation approach.

RESULTS

A sample electroencephalography dataset was compiled to test all the MVPAlab main functionalities. Significant clusters (p<0.01) were found for the proposed decoding analyses and different configurations, proving the software capability for discriminating between different experimental conditions.

CONCLUSIONS

This toolbox has been designed to include an easy-to-use and intuitive graphic user interface and data representation software, which makes MVPAlab a very convenient tool for users with few or no previous coding experience. In addition, MVPAlab is not for beginners only, as it implements several high and low-level routines allowing more experienced users to design their own projects in a highly flexible manner.

Reproducibility of graph measures derived from resting-state MEG functional connectivity metrics in sensor and source spaces

Prior studies have used graph analysis of resting‐state magnetoencephalography (MEG) to characterize abnormal brain networks in neurological disorders. However, a present challenge for researchers is the lack of guidance on which network construction strategies to employ. The reproducibility of graph measures is important for their use as clinical biomarkers. Furthermore, global graph measures should ideally not depend on whether the analysis was performed in the sensor or source space. Therefore, MEG data of the 89 healthy subjects of the Human Connectome Project were used to investigate test–retest reliability and sensor versus source association of global graph measures. Atlas‐based beamforming was used for source reconstruction, and functional connectivity (FC) was estimated for both sensor and source signals in six frequency bands using the debiased weighted phase lag index (dwPLI), amplitude envelope correlation (AEC), and leakage‐corrected AEC. Reliability was examined over multiple network density levels achieved with proportional weight and orthogonal minimum spanning tree thresholding. At a 100% density, graph measures for most FC metrics and frequency bands had fair to excellent reliability and significant sensor versus source association. The greatest reliability and sensor versus source association was obtained when using amplitude metrics. Reliability was similar between sensor and source spaces when using amplitude metrics but greater for the source than the sensor space in higher frequency bands when using the dwPLI. These results suggest that graph measures are useful biomarkers, particularly for investigating functional networks based on amplitude synchrony.

Dynamic large-scale connectivity of intrinsic cortical oscillations supports adaptive listening in challenging conditions

In multi-talker situations, individuals adapt behaviorally to this listening challenge mostly with ease, but how do brain neural networks shape this adaptation? We here establish a long-sought link between large-scale neural communications in electrophysiology and behavioral success in the control of attention in difficult listening situations. In an age-varying sample of N = 154 individuals, we find that connectivity between intrinsic neural oscillations extracted from source-reconstructed electroencephalography is regulated according to the listener's goal during a challenging dual-talker task. These dynamics occur as spatially organized modulations in power-envelope correlations of alpha and low-beta neural oscillations during approximately 2-s intervals most critical for listening behavior relative to resting-state baseline. First, left frontoparietal low-beta connectivity (16 to 24 Hz) increased during anticipation and processing of a spatial-attention cue before speech presentation. Second, posterior alpha connectivity (7 to 11 Hz) decreased during comprehension of competing speech, particularly around target-word presentation. Connectivity dynamics of these networks were predictive of individual differences in the speed and accuracy of target-word identification, respectively, but proved unconfounded by changes in neural oscillatory activity strength. Successful adaptation to a listening challenge thus latches onto two distinct yet complementary neural systems: a beta-tuned frontoparietal network enabling the flexible adaptation to attentive listening state and an alpha-tuned posterior network supporting attention to speech.

Mild traumatic brain injury impairs the coordination of intrinsic and motor-related neural dynamics

•

MTBI is poorly understood and lacks objective diagnostic and prognostic tools.

•

Abnormal neural oscillations are found in subjects with a history of mTBI.

•

We identify transient bursts in MEG data using a Hidden Markov Model.

•

We explain a deficit in beta connectivity and power in terms of transient bursts.

•

Data-driven feature selection identifies symptom-relevant functional connections.

Mild traumatic brain injury (mTBI) poses a considerable burden on healthcare systems. Whilst most patients recover quickly, a significant number suffer from sequelae that are not accompanied by measurable structural damage. Understanding the neural underpinnings of these debilitating effects and developing a means to detect injury, would address an important unmet clinical need. It could inform interventions and help predict prognosis. Magnetoencephalography (MEG) affords excellent sensitivity in probing neural function and presents significant promise for assessing mTBI, with abnormal neural oscillations being a potential specific biomarker. However, growing evidence suggests that neural dynamics are (at least in part) driven by transient, pan-spectral bursting and in this paper, we employ this model to investigate mTBI. We applied a Hidden Markov Model to MEG data recorded during resting state and a motor task and show that previous findings of diminished intrinsic beta amplitude in individuals with mTBI are largely due to the reduced beta band spectral content of bursts, and that diminished beta connectivity results from a loss in the temporal coincidence of burst states. In a motor task, mTBI results in diminished burst amplitude, altered modulation of burst probability during movement, and a loss in connectivity in motor networks. These results suggest that, mechanistically, mTBI disrupts the structural framework underlying neural synchrony, which impairs network function. Whilst the damage may be too subtle for structural imaging to see, the functional consequences are detectable and persist after injury. Our work shows that mTBI impairs the dynamic coordination of neural network activity and proposes a potent new method for understanding mTBI.

Cross-Frequency Coupling in Childhood Absence Epilepsy

Objective: Absence seizures are the prototypic primarily generalized seizures, but there is incomplete understanding regarding their generation and maintenance. A core network for absence seizures has been defined, including focal cortical and thalamic regions that have frequency-dependent interactions. The purpose of this study was to investigate within-frequency coupling and cross-frequency coupling (CFC) during human absence seizures, to identify key regions (hubs) within the absence network that contribute to propagation and maintenance. Methods: Thirteen children with new-onset and untreated childhood absence epilepsy had over 60 typical absence seizures during both electroencephalography-functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) recordings. The spatial map of the ictal network was defined using fMRI and used as prior information for MEG connectivity. A multilayer network approach was used to investigate within-frequency coupling and CFC for canonical frequency bands. A rigorous null-modeling approach was used to determine connections outside the noise floor. Results: Strong coupling between beta and gamma frequencies, within the left frontal cortex, and between the left frontal and right parietal regions was observed. There was also strong connectivity between left frontal and right parietal nodes within the gamma band. Multilayer versatility analysis identified a cluster of network hubs in the left frontal region. Interpretation: Cortical regions commonly identified as being critical for absence seizure generation (frontal cortex, precuneus) have strong CFC and within-frequency coupling between beta and gamma bands. As nonpharmacologic treatments, such as neuromodulation, become available for generalized epilepsies, detailed mechanistic understanding of how "diffuse" seizures are generated and maintained will be necessary to provide optimal outcomes.

Testing covariance models for MEG source reconstruction of hippocampal activity

Beamforming is one of the most commonly used source reconstruction methods for magneto- and electroencephalography (M/EEG). One underlying assumption, however, is that distant sources are uncorrelated and here we tested whether this is an appropriate model for the human hippocampal data. We revised the Empirical Bayesian Beamfomer (EBB) to accommodate specific a-priori correlated source models. We showed in simulation that we could use model evidence (as approximated by Free Energy) to distinguish between different correlated and uncorrelated source scenarios. Using group MEG data in which the participants performed a hippocampal-dependent task, we explored the possibility that the hippocampus or the cortex or both were correlated in their activity across hemispheres. We found that incorporating a correlated hippocampal source model significantly improved model evidence. Our findings help to explain why, up until now, the majority of MEG-reported hippocampal activity (typically making use of beamformers) has been estimated as unilateral.

Identification of psychiatric disorder subtypes from functional connectivity patterns in resting-state electroencephalography

The understanding and treatment of psychiatric disorders, which are known to be neurobiologically and clinically heterogeneous, could benefit from the data-driven identification of disease subtypes. Here, we report the identification of two clinically relevant subtypes of post-traumatic stress disorder (PTSD) and major depressive disorder (MDD) on the basis of robust and distinct functional connectivity patterns, prominently within the frontoparietal control network and the default mode network. We identified the disease subtypes by analysing, via unsupervised and supervised machine learning, the power-envelope-based connectivity of signals reconstructed from high-density resting-state electroencephalography in four datasets of patients with PTSD and MDD, and show that the subtypes are transferable across independent datasets recorded under different conditions. The subtype whose functional connectivity differed most from those of healthy controls was less responsive to psychotherapy treatment for PTSD and failed to respond to an antidepressant medication for MDD. By contrast, both subtypes responded equally well to two different forms of repetitive transcranial magnetic stimulation therapy for MDD. Our data-driven approach may constitute a generalizable solution for connectome-based diagnosis.

Mild traumatic brain injury is associated with dysregulated neural network functioning in children and adolescents

Mild traumatic brain injury is highly prevalent in paediatric populations, and can result in chronic physical, cognitive and emotional impairment, known as persistent post-concussive symptoms. Magnetoencephalography has been used to investigate neurophysiological dysregulation in mild traumatic brain injury in adults; however, whether neural dysrhythmia persists in chronic mild traumatic brain injury in children and adolescents is largely unknown. We predicted that children and adolescents would show similar dysfunction as adults, including pathological slow-wave oscillations and maladaptive, frequency-specific, alterations to neural connectivity. Using magnetoencephalography, we investigated regional oscillatory power and distributed brain-wide networks in a cross-sectional sample of children and adolescents in the chronic stages of mild traumatic brain injury. Additionally, we used a machine learning pipeline to identify the most relevant magnetoencephalography features for classifying mild traumatic brain injury and to test the relative classification performance of regional power versus functional coupling. Results revealed that the majority of participants with chronic mild traumatic brain injury reported persistent post-concussive symptoms. For neurophysiological imaging, we found increased regional power in the delta band in chronic mild traumatic brain injury, predominantly in bilateral occipital cortices and in the right inferior temporal gyrus. Those with chronic mild traumatic brain injury also showed dysregulated neuronal coupling, including decreased connectivity in the delta range, as well as hyper-connectivity in the theta, low gamma and high gamma bands, primarily involving frontal, temporal and occipital brain areas. Furthermore, our multivariate classification approach combined with functional connectivity data outperformed regional power in terms of between-group classification accuracy. For the first time, we establish that local and large-scale neural activity are altered in youth in the chronic phase of mild traumatic brain injury, with the majority presenting persistent post-concussive symptoms, and that dysregulated interregional neural communication is a reliable marker of lingering paediatric ‘mild’ traumatic brain injury.

Abnormal meta-state activation of dynamic brain networks across the Alzheimer spectrum

The characterization of the distinct dynamic functional connectivity (dFC) patterns that activate in the brain during rest can help to understand the underlying time-varying network organization. The presence and behavior of these patterns (known as meta-states) have been widely studied by means of functional magnetic resonance imaging (fMRI). However, modalities with high-temporal resolution, such as electroencephalography (EEG), enable the characterization of fast temporally evolving meta-state sequences. Mild cognitive impairment (MCI) and dementia due to Alzheimer's disease (AD) have been shown to disrupt spatially localized activation and dFC between different brain regions, but not much is known about how they affect meta-state network topologies and their network dynamics. The main hypothesis of the study was that MCI and dementia due to AD alter normal meta-state sequences by inducing a loss of structure in their patterns and a reduction of their dynamics. Moreover, we expected that patients with MCI would display more flexible behavior compared to patients with dementia due to AD. Thus, the aim of the current study was twofold: (i) to find repeating, distinctly organized network patterns (meta-states) in neural activity; and (ii) to extract information about meta-state fluctuations and how they are influenced by MCI and dementia due to AD. To accomplish these goals, we present a novel methodology to characterize dynamic meta-states and their temporal fluctuations by capturing aspects based on both their discrete activation and the continuous evolution of their individual strength. These properties were extracted from 60-s resting-state EEG recordings from 67 patients with MCI due to AD, 50 patients with dementia due to AD, and 43 cognitively healthy controls. First, the instantaneous amplitude correlation (IAC) was used to estimate instantaneous functional connectivity with a high temporal resolution. We then extracted meta-states by means of graph community detection based on recurrence plots (RPs), both at the individual- and group-level. Subsequently, a diverse set of properties of the continuous and discrete fluctuation patterns of the meta-states was extracted and analyzed. The main novelty of the methodology lies in the usage of Louvain GJA community detection to extract meta-states from IAC-derived RPs and the extended analysis of their discrete and continuous activation. Our findings showed that distinct dynamic functional connectivity meta-states can be found on the EEG time-scale, and that these were not affected by the oscillatory slowing induced by MCI or dementia due to AD. However, both conditions displayed a loss of meta-state modularity, coupled with shorter dwell times and higher complexity of the meta-state sequences. Furthermore, we found evidence that meta-state sequencing is not entirely random; it shows an underlying structure that is partially lost in MCI and dementia due to AD. These results show evidence that AD progression is associated with alterations in meta-state switching, and a degradation of dynamic brain flexibility.

Multilayer Connector Hub Mapping Reveals Key Brain Regions Supporting Expressive Language

Introduction: How components of the distributed brain networks that support cognition participate in typical functioning remains a largely unanswered question. An important subgroup of regions in the larger network are connector hubs, which are areas that are highly connected to several other functionally specialized sets of regions, and are likely important for sensorimotor integration. The present study attempts to characterize connector hubs involved in typical expressive language functioning using a data-driven, multimodal, full multilayer magnetoencephalography (MEG) connectivity-based pipeline. Methods: Twelve adolescents, 16-18 years of age (five males), participated in this study. Participants underwent MEG scanning during a verb generation task. MEG and structural connectivity were calculated at the whole-brain level. Amplitude/amplitude coupling (AAC) was used to compute functional connections both within and between discrete frequency bins. AAC values were then multiplied by a binary structural connectivity matrix, and then entered into full multilayer network analysis. Initially, hubs were defined based on multilayer versatility and subsequently reranked by a novel measure called delta centrality on interconnectedness (DCI). DCI is defined as the percent change in interfrequency interconnectedness after removal of a hub. Results: We resolved regions that are important for between-frequency communication among other areas during expressive language, with several potential theoretical and clinical applications that can be generalized to other cognitive domains. Conclusion: Our multilayer, data-driven framework captures nonlinear connections that span across scales that are often missed in conventional analyses. The present study suggests that crucial hubs may be conduits for interfrequency communication between action and perception systems that are crucial for typical functioning.

Lag-invariant detection of interactions in spatially-extended systems using linear inverse modeling

Measurements on physical systems result from the systems' activity being converted into sensor measurements by a forward model. In a number of cases, inversion of the forward model is extremely sensitive to perturbations such as sensor noise or numerical errors in the forward model. Regularization is then required, which introduces bias in the reconstruction of the systems' activity. One domain in which this is particularly problematic is the reconstruction of interactions in spatially-extended complex systems such as the human brain. Brain interactions can be reconstructed from non-invasive measurements such as electroencephalography (EEG) or magnetoencephalography (MEG), whose forward models are linear and instantaneous, but have large null-spaces and high condition numbers. This leads to incomplete unmixing of the forward models and hence to spurious interactions. This motivated the development of interaction measures that are exclusively sensitive to lagged, i.e. delayed interactions. The drawback of such measures is that they only detect interactions that have sufficiently large lags and this introduces bias in reconstructed brain networks. We introduce three estimators for linear interactions in spatially-extended systems that are uniformly sensitive to all lags. We derive some basic properties of and relationships between the estimators and evaluate their performance using numerical simulations from a simple benchmark model.

Distinct Patterns of Functional Connectivity During the Comprehension of Natural, Narrative Speech

Recent continuous task studies, such as narrative speech comprehension, show that fluctuations in brain functional connectivity (FC) are altered and enhanced compared to the resting state. Here, we characterized the fluctuations in FC during comprehension of speech and time-reversed speech conditions. The correlations of Hilbert envelope of source-level EEG data were used to quantify FC between spatially separate brain regions. A symmetric multivariate leakage correction was applied to address the signal leakage issue before calculating FC. The dynamic FC was estimated based on a sliding time window. Then, principal component analysis (PCA) was performed on individually concatenated and temporally concatenated FC matrices to identify FC patterns. We observed that the mode of FC induced by speech comprehension can be characterized with a single principal component. The condition-specific FC demonstrated decreased correlations between frontal and parietal brain regions and increased correlations between frontal and temporal brain regions. The fluctuations of the condition-specific FC characterized by a shorter time demonstrated that dynamic FC also exhibited condition specificity over time. The FC is dynamically reorganized and FC dynamic pattern varies along a single mode of variation during speech comprehension. The proposed analysis framework seems valuable for studying the reorganization of brain networks during continuous task experiments.

Frequency-dependent functional connectivity in resting state networks

Functional magnetic resonance imaging studies have documented the resting human brain to be functionally organized in multiple large‐scale networks, called resting‐state networks (RSNs). Other brain imaging techniques, such as electroencephalography (EEG) and magnetoencephalography (MEG), have been used for investigating the electrophysiological basis of RSNs. To date, it is largely unclear how neural oscillations measured with EEG and MEG are related to functional connectivity in the resting state. In addition, it remains to be elucidated whether and how the observed neural oscillations are related to the spatial distribution of the network nodes over the cortex. To address these questions, we examined frequency‐dependent functional connectivity between the main nodes of several RSNs, spanning large part of the cortex. We estimated connectivity using band‐limited power correlations from high‐density EEG data collected in healthy participants. We observed that functional interactions within RSNs are characterized by a specific combination of neuronal oscillations in the alpha (8–13 Hz), beta (13–30 Hz), and gamma (30–80 Hz) bands, which highly depend on the position of the network nodes. This finding may contribute to a better understanding of the mechanisms through which neural oscillations support functional connectivity in the brain.

Using structural connectivity to augment community structure in EEG functional connectivity

Recently, EEG recording techniques and source analysis have improved, making it feasible to tap into fast network dynamics. Yet, analyzing whole-cortex EEG signals in source space is not standard, partly because EEG suffers from volume conduction: Functional connectivity (FC) reflecting genuine functional relationships is impossible to disentangle from spurious FC introduced by volume conduction. Here, we investigate the relationship between white matter structural connectivity (SC) and large-scale network structure encoded in EEG-FC. We start by confirming that FC (power envelope correlations) is predicted by SC beyond the impact of Euclidean distance, in line with the assumption that SC mediates genuine FC. We then use information from white matter structural connectivity in order to smooth the EEG signal in the space spanned by graphs derived from SC. Thereby, FC between nearby, structurally connected brain regions increases while FC between nonconnected regions remains unchanged, resulting in an increase in genuine, SC-mediated FC. We analyze the induced changes in FC, assessing the resemblance between EEG-FC and volume-conduction- free fMRI-FC, and find that smoothing increases resemblance in terms of overall correlation and community structure. This result suggests that our method boosts genuine FC, an outcome that is of interest for many EEG network neuroscience questions.

Magnetoencephalography resting state connectivity patterns as indicatives of surgical outcome in epilepsy patients

OBJECTIVE

Focal epilepsy is a disorder affecting several brain networks; however, epilepsy surgery usually targets a restricted region, the so-called epileptic focus. There is a growing interest in embedding resting state (RS) connectivity analysis into pre-surgical workup.

APPROACH

In this retrospective study, we analyzed Magnetoencephalography (MEG) long-range RS functional connectivity patterns in patients with drug-resistant focal epilepsy. MEG recorded prior to surgery from seven seizure-free (Engel Ia) and five non seizure-free (Engel III or IV) patients were analyzed (minimum 2-years post-surgical follow-up). MEG segments without any detectable epileptic activity were source localized using wavelet-based Maximum Entropy on the Mean method. Amplitude envelope correlation in the theta (4-8 Hz), alpha (8-13 Hz), and beta (13-26 Hz) bands were used for assessing connectivity.

MAIN RESULTS

For seizure-free patients, we found an isolated epileptic network characterized by weaker connections between the brain region where interictal epileptic discharges (IED) are generated and the rest of the cortex, when compared to connectivity between the corresponding contralateral homologous region and the rest of the cortex. Contrarily, non seizure-free patients exhibited a widespread RS epileptic network characterized by stronger connectivity between the IED generator and the rest of the cortex, in comparison to the contralateral region and the cortex. Differences between the two seizure outcome groups concerned mainly distant long-range connections and were found in the alpha-band.

SIGNIFICANCE

Importantly, these connectivity patterns suggest specific mechanisms describing the underlying organization of the epileptic network and were detectable at the individual patient level, supporting the prospect use of MEG connectivity patterns in epilepsy to predict post-surgical seizure outcome.

Discovering dynamic task-modulated functional networks with specific spectral modes using MEG

Efficient neuronal communication between brain regions through oscillatory synchronization at certain frequencies is necessary for cognition. Such synchronized networks are transient and dynamic, established on the timescale of milliseconds in order to support ongoing cognitive operations. However, few studies characterizing dynamic electrophysiological brain networks have simultaneously accounted for temporal non-stationarity, spectral structure, and spatial properties. Here, we propose an analysis framework for characterizing the large-scale phase-coupling network dynamics during task performance using magnetoencephalography (MEG). We exploit the high spatiotemporal resolution of MEG to measure time-frequency dynamics of connectivity between parcellated brain regions, yielding data in tensor format. We then use a tensor component analysis (TCA)-based procedure to identify the spatio-temporal-spectral modes of covariation among separate regions in the human brain. We validate our pipeline using MEG data recorded during a hand movement task, extracting a transient motor network with beta-dominant spectral mode, which is significantly modulated by the movement task. Next, we apply the proposed pipeline to explore brain networks that support cognitive operations during a working memory task. The derived results demonstrate the temporal formation and dissolution of multiple phase-coupled networks with specific spectral modes, which are associated with face recognition, vision, and movement. The proposed pipeline can characterize the spectro-temporal dynamics of functional connectivity in the brain on the subsecond timescale, commensurate with that of cognitive performance.

Whole-brain electrophysiological functional connectivity dynamics in resting-state EEG

OBJECTIVE

Functional connectivity (FC) dynamics have been studied in functional magnetic resonance imaging (fMRI) data, while it is largely unknown in electrophysiological data, e.g. EEG.

APPROACH

The present study proposed a novel analytic framework to study spatiotemporal dynamics of FC (dFC) in resting-state human EEG data, including independent component analysis, cortical source imaging, sliding-window correlation analysis, and k-means clustering.

MAIN RESULTS

Our results confirm that major fMRI intrinsic connectivity networks (ICNs) can be successfully reconstructed from EEG using our analytic framework. Prominent spatial and temporal variability were revealed in these ICNs. The mean dFC spatial patterns of individual ICNs resemble their corresponding static FC (sFC) patterns but show fewer cross-talks among distinct ICNs. Our investigation unveils evidences of time-domain variations in individual ICNs comparable to their mean FC level in terms of magnitude. The major contributors to these variations are from the frequency below 0.0156 Hz, in the similar range of FC dynamics from fMRI data. Among different ICNs, larger temporal variabilities are observed in the frontal attention and auditory/visual ICNs, while sensorimotor, salience, and default model networks showed less. Our analytic framework for the first time revealed quasi-stable states within individual EEG ICNs, with various strengths or spatial patterns that were reliably detected at both group and individual levels. These states all together reveal a more complete picture of EEG ICNs: (1) quasi-stable state spatial patterns as a whole for each EEG ICN are more consistent with the corresponding fMRI ICN in terms of the bilateral distribution and multi-nodes structure; (2) EEG ICNs reveal more transient patterns about within-ICN between-node communications than fMRI ICNs.

SIGNIFICANCE

The present findings highlight the fact that rich temporal and spatial dynamics exist in ICN that can be detected from EEG data. Future studies might extend investigations towards spectral dynamics of EEG ICNs.

Dissociated neuronal phase- and amplitude-coupling patterns in the human brain

Coupling of neuronal oscillations may reflect and facilitate the communication between neuronal populations. Two primary neuronal coupling modes have been described: phase-coupling and amplitude-coupling. Theoretically, both coupling modes are independent, but so far, their neuronal relationship remains unclear. Here, we combined MEG, source-reconstruction and simulations to systematically compare cortical amplitude-coupling and phase-coupling patterns in the human brain. Importantly, we took into account a critical bias of amplitude-coupling measures due to phase-coupling. We found differences between both coupling modes across a broad frequency range and most of the cortex. Furthermore, by combining empirical measurements and simulations we ruled out that these results were caused by methodological biases, but instead reflected genuine neuronal amplitude coupling. Our results show that cortical phase- and amplitude-coupling patterns are non-redundant, which may reflect at least partly distinct neuronal mechanisms. Furthermore, our findings highlight and clarify the compound nature of amplitude coupling measures.

Frequency-Specific Resting Connectome in Bipolar Disorder: An MEG Study

BACKGROUND

Bipolar disorder (BD) is a serious psychiatric disorder that is associated with a high suicide rate, and for which no clinical biomarker has yet been identified. To address this issue, we investigated the use of magnetoencephalography (MEG) as a new prospective tool. MEG has been used to evaluate frequency-specific connectivity between brain regions; however, no previous study has investigated the frequency-specific resting-state connectome in patients with BD. This resting-state MEG study explored the oscillatory representations of clinical symptoms of BD via graph analysis.

METHODS

In this prospective case-control study, 17 patients with BD and 22 healthy controls (HCs) underwent resting-state MEG and evaluations for depressive and manic symptoms. After estimating the source current distribution, orthogonalized envelope correlations between multiple brain regions were evaluated for each frequency band. We separated regions-of-interest into seven left and right network modules, including the frontoparietal network (FPN), limbic network (LM), salience network (SAL), and default mode network (DMN), to compare the intra- and inter-community edges between the two groups.

RESULTS

In the BD group, we found significantly increased inter-community edges of the right LM-right DMN at the gamma band, and decreased inter-community edges of the right SAL-right FPN at the delta band and the left SAL-right SAL at the theta band. Intra-community edges in the left LM at the high beta band were significantly higher in the BD group than in the HC group. The number of connections in the left LM at the high beta band showed positive correlations with the subjective and objective depressive symptoms in the BD group.

CONCLUSION

We introduced graph theory into resting-state MEG studies to investigate the functional connectivity in patients with BD. To the best of our knowledge, this is a novel approach that may be beneficial in the diagnosis of BD. This study describes the spontaneous oscillatory brain networks that compensate for the time-domain issues associated with functional magnetic resonance imaging. These findings suggest that the connectivity of the LM at the beta band may be a good objective biological biomarker of the depressive symptoms associated with BD.

Test-Retest Reliability of Magnetoencephalography Resting-State Functional Connectivity in Schizophrenia

The reliability of magnetoencephalography (MEG) resting-state functional connectivity in schizophrenia (SZ) is unknown as previous research has focused on healthy controls (HC). Here, we examined reliability in 26 participants (13-SZ, 13-HC). Eyes opened and eyes closed resting-state data were collected on 4 separate occasions during 2 visits, 1 week apart. For source modeling, we used minimum norm software to apply dynamic statistical parametric mapping. Source analyses compared the following functional connectivity metrics from each data run: coherence (coh), imaginary coherence (imcoh), pairwise phase consistency (ppc), phase-locking value (plv), phase lag index (pli), weighted phase lag index (wpli), and weighted phase lag index debiased (wpli2). Intraclass correlation coefficients (ICCs) were calculated for whole brain, network, and network pair averages. For reliability, ICCs above 0.75 = excellent, above 0.60 = good, above 0.40 = fair, and below 0.40 = poor reliability. We found the reliability of these metrics varied greatly depending on frequency band, network, network pair, and participant group examined. Broadband (1-58 Hz) whole brain averages in both HC and SZ showed excellent reliability for wpli2, and good to fair reliability for ppc, plv, and coh. Broadband network averages showed excellent to good reliability across 1 hour and 1 week for coh, imcoh, ppc, plv, wpli within default mode, cognitive control, and visual networks in HC, while the same metrics had excellent to fair reliability in SZ. Regional network pair averages showed good to fair reliability for coh, ppc, plv within default mode, cognitive control and visual network pairs in HC and SZ. In general, HC had higher reliability compared to SZ, and the default mode, cognitive control, and visual networks had higher reliability compared to somatosensory and auditory networks. Similar reliability levels occurred for both eyes opened and eyes closed resting-states for most metrics. The functional connectivity metrics of coh, ppc, and plv performed best across 1 hour and 1 week in HC and SZ. We also found that SZ had reduced coh, plv, and ppc in the dmn average and pair values indicating dysconnectivity in SZ. These findings encourage collecting both eyes opened and eyes closed resting-state MEG, while demonstrating that clinical populations may differ in reliability.

Towards OPM-MEG in a virtual reality environment

Virtual reality (VR) provides an immersive environment in which a participant can experience a feeling of presence in a virtual world. Such environments generate strong emotional and physical responses and have been used for wide-ranging applications. The ability to collect functional neuroimaging data whilst a participant is immersed in VR would represent a step change for experimental paradigms; unfortunately, traditional brain imaging requires participants to remain still, limiting the scope of naturalistic interaction within VR. Recently however, a new type of magnetoencephalography (MEG) device has been developed, that employs scalp-mounted optically-pumped magnetometers (OPMs) to measure brain electrophysiology. Lightweight OPMs, coupled with precise control of the background magnetic field, enables participant movement during data acquisition. Here, we exploit this technology to acquire MEG data whilst a participant uses a virtual reality head-mounted display (VRHMD). We show that, despite increased magnetic interference from the VRHMD, we were able to measure modulation of alpha-band oscillations, and the visual evoked field. Moreover, in a VR experiment in which a participant had to move their head to look around a virtual wall and view a visual stimulus, we showed that the measured MEG signals map spatially in accordance with the known organisation of primary visual cortex. This technique could transform the type of neuroscientific experiment that can be undertaken using functional neuroimaging.

Brain Functional Connectivity Through Phase Coupling of Neuronal Oscillations: A Perspective From Magnetoencephalography

Magnetoencephalography has gained an increasing importance in systems neuroscience thanks to the possibility it offers of unraveling brain networks at time-scales relevant to behavior, i.e., frequencies in the 1-100 Hz range, with sufficient spatial resolution. In the first part of this review, we describe, in a unified mathematical framework, a large set of metrics used to estimate MEG functional connectivity at the same or at different frequencies. The different metrics are presented according to their characteristics: same-frequency or cross-frequency, univariate or multivariate, directed or undirected. We focus on phase coupling metrics given that phase coupling of neuronal oscillations is a putative mechanism for inter-areal communication, and that MEG is an ideal tool to non-invasively detect such coupling. In the second part of this review, we present examples of the use of specific phase methods on real MEG data in the context of resting state, visuospatial attention and working memory. Overall, the results of the studies provide evidence for frequency specific and/or cross-frequency brain circuits which partially overlap with brain networks as identified by hemodynamic-based imaging techniques, such as functional Magnetic Resonance (fMRI). Additionally, the relation of these functional brain circuits to anatomy and to behavior highlights the usefulness of MEG phase coupling in systems neuroscience studies. In conclusion, we believe that the field of MEG functional connectivity has made substantial steps forward in the recent years and is now ready for bringing the study of brain networks to a more mechanistic understanding.

Uncovering Dynamic Brain Reconfiguration in MEG Working Memory n-Back Task Using Topological Data Analysis

The increasing availability of high temporal resolution neuroimaging data has increased the efforts to understand the dynamics of neural functions. Until recently, there are few studies on generative models supporting classification and prediction of neural systems compared to the description of the architecture. However, the requirement of collapsing data spatially and temporally in the state-of-the art methods to analyze functional magnetic resonance imaging (fMRI), electroencephalogram (EEG) and magnetoencephalography (MEG) data cause loss of important information. In this study, we addressed this issue using a topological data analysis (TDA) method, called Mapper, which visualizes evolving patterns of brain activity as a mathematical graph. Accordingly, we analyzed preprocessed MEG data of 83 subjects from Human Connectome Project (HCP) collected during working memory n-back task. We examined variation in the dynamics of the brain states with the Mapper graphs, and to determine how this variation relates to measures such as response time and performance. The application of the Mapper method to MEG data detected a novel neuroimaging marker that explained the performance of the participants along with the ground truth of response time. In addition, TDA enabled us to distinguish two task-positive brain activations during 0-back and 2-back tasks, which is hard to detect with the other pipelines that require collapsing the data in the spatial and temporal domain. Further, the Mapper graphs of the individuals also revealed one large group in the middle of the stimulus detecting the high engagement in the brain with fine temporal resolution, which could contribute to increase spatiotemporal resolution by merging different imaging modalities. Hence, our work provides another evidence to the effectiveness of the TDA methods for extracting subtle dynamic properties of high temporal resolution MEG data without the temporal and spatial collapse.