Ghost interactions in MEG/EEG source space: A note of caution on inter-areal coupling measures

Simultaneous invasive and non-invasive recordings in humans: A novel Rosetta stone for deciphering brain activity

Simultaneous noninvasive and invasive electrophysiological recordings provide a unique opportunity to achieve a comprehensive understanding of human brain activity, much like a Rosetta stone for human neuroscience. In this review we focus on the increasingly-used powerful combination of intracranial electroencephalography (iEEG) with scalp electroencephalography (EEG) or magnetoencephalography (MEG). We first provide practical insight on how to achieve these technically challenging recordings. We then provide examples from clinical research on how simultaneous recordings are advancing our understanding of epilepsy. This is followed by the illustration of how human neuroscience and methodological advances could benefit from these simultaneous recordings. We conclude with a call for open data sharing and collaboration, while ensuring neuroethical approaches and argue that only with a true collaborative approach the promises of simultaneous recordings will be fulfilled.

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.

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.

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.

OPM-MEG Measuring Phase Synchronization on Source Time Series: Application in Rhythmic Median Nerve Stimulation

The magnetoencephalogram (MEG) based on array optically pumped magnetometers (OPMs) has the potential of replacing conventional cryogenic superconducting quantum interference device. Phase synchronization is a common method for measuring brain oscillations and functional connectivity. Verifying the feasibility and fidelity of OPM-MEG in measuring phase synchronization will help its widespread application in the study of aforementioned neural mechanisms. The analysis method on source-level time series can weaken the influence of instantaneous field spread effect. In this paper, the OPM-MEG was used for measuring the evoked responses of 20Hz rhythmic and arrhythmic median nerve stimulation, and the inter-trial phase synchronization (ITPS) and inter-reginal phase synchronization (IRPS) of primary somatosensory cortex (SI) and secondary somatosensory cortex (SII) were analysed. The results find that under rhythmic condition, the evoked responses of SI and SII show continuous oscillations and the effect of resetting phase. The values of ITPS and IRPS significantly increase at the stimulation frequency of 20Hz and its harmonic of 40Hz, whereas the arrhythmic stimulation does not exhibit this phenomenon. Moreover, in the initial stage of stimulation, the ITPS and IRPS values are significantly higher at Mu rhythm in the rhythmic condition compared to arrhythmic. In conclusion, the results demonstrate the ability of OPM-MEG in measuring phase pattern and functional connectivity on source-level, and may also prove beneficial for the study on the mechanism of rhythmic stimulation therapy for rehabilitation.

Electroencephalography connectome changes in chronic insomnia disorder are correlated with neurochemical signatures

STUDY OBJECTIVES

This study aimed to investigate the alterations in resting-state electroencephalography (EEG) global brain connectivity (GBC) in patients with chronic insomnia disorder (CID) and to explore the correlation between macroscale connectomic variances and microscale neurotransmitter distributions.

METHODS

We acquired 64-channel EEG from 35 female CID patients and 34 healthy females. EEG signals were source-localized using individual brain anatomy and orthogonalized to mitigate volume conduction. Correlation coefficients between band-limited source-space power envelopes of the DK 68 atlas were computed and averaged across regions to determine specific GBC values. A support vector machine (SVM) classifier utilizing GBC features was employed to differentiate CID patients from controls. We further used Neurosynth and a 3D atlas of neurotransmitter receptors/transporters to assess the cognitive functions and neurotransmitter landscape associated with CID cortical abnormality maps, respectively.

RESULTS

CID patients exhibited elevated GBC within the medial prefrontal cortex and limbic cortex, particularly at the gamma carrier frequency, compared to controls (pFDR<0.05). GBC patterns were found to effectively distinguish CID patients from controls with a precision of 90.8% in the SVM model. The cortical abnormality maps were significantly correlated with meta-analytic terms like "cognitive control" and "emotion regulation." Notably, GBC patterns were associated with neurotransmitter profiles (pspin<0.05), with neurotransmitter systems such as norepinephrine, dopamine, and serotonin making significant contributions.

CONCLUSIONS

This work characterizes the EEG connectomic profile of CID, facilitating the cost-effective clinical translation of EEG-derived markers. Additionally, the linkage between GBC patterns and neurotransmitter distribution offers promising avenues for developing targeted treatment strategies for CID.

Characterization of Second-Order Mixing Effects in Reconstructed Cross-Spectra of Random Neural Fields

Functional connectivity in electroencephalography (EEG) and magnetoencephalography (MEG) data is commonly assessed by using measures that are insensitive to instantaneously interacting sources and as such would not give rise to false positive interactions caused by instantaneous mixing of true source signals (first-order mixing). Recent studies, however, have drawn attention to the fact that such measures are still susceptible to instantaneous mixing from lagged sources (i.e. second-order mixing) and that this can lead to a large number of false positive interactions. In this study we relate first- and second-order mixing effects on the cross-spectra of reconstructed source activity to the properties of the resolution operators that are used for the reconstruction. We derive two identities that relate first- and second-order mixing effects to the transformation properties of measurement and source configurations and exploit them to establish several basic properties of signal mixing. First, we provide a characterization of the configurations that are maximally and minimally sensitive to second-order mixing. It turns out that second-order mixing effects are maximal when the measurement locations are far apart and the sources coincide with the measurement locations. Second, we provide a description of second-order mixing effects in the vicinity of the measurement locations in terms of the local geometry of the point-spread functions of the resolution operator. Third, we derive a version of Lagrange's identity for cross-talk functions that establishes the existence of a trade-off between the magnitude of first- and second-order mixing effects. It also shows that, whereas the magnitude of first-order mixing is determined by the inner product of cross-talk functions, the magnitude of second-order mixing is determined by a generalized cross-product of cross-talk functions (the wedge product) which leads to an intuitive geometric understanding of the trade-off. All results are derived within the general framework of random neural fields on cortical manifolds.

Connectivity of neural signals to the primary motor area during preparatory periods for movement following external and internal cues

PURPOSE

We investigated the connectivity of neural signals from movement-related cortical areas to the primary motor area (M1) in the hemisphere contralateral to the movement side during the period of movement-related magnetic fields before movement.

MATERIALS AND METHODS

Participants were 13 healthy adults, and nerual signals were recorded using magnetoencephalography. Spontaneous extension of the right wrist was performed at the participant's own pace and following a visual cue in internal (IC) and external (EC) cue tasks. The connectivity of neural signals to M1 from each movement-related motor area was assessed by Granger causality analysis (GCA). The GCA was performed on the neural activity elicited in a frequency band between 7.8 and 46.9 Hz during the pre-movement periods, which occurred durng the readiness field (RF) and the negative slope prime (NSp). F-values, as connectivity values obtained by GCA, were compared between the EC and IC cue tasks.

RESULTS

For NSp periods, the connectivity of neural signals from the left superior frontal area (SF-L) to M1 was dominant in the IC task, whereas that from the left superior parietal area (SP-L) to M1 was dominant in the EC task. The F value in the GCA from SP-L to M1 was greater in the EC task during RF than in the IC task during equivalent periods.

CONSLUSIONS

In the present study, there were differences in the connectivity of neural signals to M1 between IC and EC tasks. The present results suggested that the pattern of pre-movement neural activity that resulted in a movement was not uniform but differed between movement tasks just before the movement.

A recurrent Hopfield network for estimating meso-scale effective connectivity in MEG

The architecture of communication within the brain, represented by the human connectome, has gained a paramount role in the neuroscience community. Several features of this communication, e.g., the frequency content, spatial topology, and temporal dynamics are currently well established. However, identifying generative models providing the underlying patterns of inhibition/excitation is very challenging. To address this issue, we present a novel generative model to estimate large-scale effective connectivity from MEG. The dynamic evolution of this model is determined by a recurrent Hopfield neural network with asymmetric connections, and thus denoted Recurrent Hopfield Mass Model (RHoMM). Since RHoMM must be applied to binary neurons, it is suitable for analyzing Band Limited Power (BLP) dynamics following a binarization process. We trained RHoMM to predict the MEG dynamics through a gradient descent minimization and we validated it in two steps. First, we showed a significant agreement between the similarity of the effective connectivity patterns and that of the interregional BLP correlation, demonstrating RHoMM's ability to capture individual variability of BLP dynamics. Second, we showed that the simulated BLP correlation connectomes, obtained from RHoMM evolutions of BLP, preserved some important topological features, e.g, the centrality of the real data, assuring the reliability of RHoMM. Compared to other biophysical models, RHoMM is based on recurrent Hopfield neural networks, thus, it has the advantage of being data-driven, less demanding in terms of hyperparameters and scalable to encompass large-scale system interactions. These features are promising for investigating the dynamics of inhibition/excitation at different spatial scales.

Strengthening of alpha synchronization is a neural correlate of cognitive transfer

Cognitive training can lead to improvements in both task-specific strategies and general capacities, such as visuo-spatial working memory (VSWM). The latter emerge slowly and linearly throughout training, in contrast to strategy where changes typically occur within the first days of training. Changes in strategy and capacity have not been separated in prior neuroimaging studies. Here, we used a within-participants design with dense temporal sampling to capture the time dynamics of neural mechanisms associated with change in capacity. In four participants, neural activity was recorded with magnetoencephalography on seven occasions over two months of visuo-spatial working memory training. During scanning, the participants performed a trained visuo-spatial working memory task, a transfer task, and a control task. First, we extracted an individual visuo-spatial working memory-load-dependent synchronization network for each participant. Next, we identified linear changes over time in the network, congruent with the temporal dynamics of capacity change. Three out of four participants showed a gradual strengthening of alpha synchronization. Strengthening of the same connections was also found in the transfer task but not in the control task. This suggests that cognitive transfer occurs through slow, gradual strengthening of alpha synchronization between cortical regions that are vital for both the trained task and the transfer task.

The time-evolving epileptic brain network: concepts, definitions, accomplishments, perspectives

Epilepsy is now considered a network disease that affects the brain across multiple levels of spatial and temporal scales. The paradigm shift from an epileptic focus-a discrete cortical area from which seizures originate-to a widespread epileptic network-spanning lobes and hemispheres-considerably advanced our understanding of epilepsy and continues to influence both research and clinical treatment of this multi-faceted high-impact neurological disorder. The epileptic network, however, is not static but evolves in time which requires novel approaches for an in-depth characterization. In this review, we discuss conceptual basics of network theory and critically examine state-of-the-art recording techniques and analysis tools used to assess and characterize a time-evolving human epileptic brain network. We give an account on current shortcomings and highlight potential developments towards an improved clinical management of epilepsy.

Networks of cortical activity show graded responses to perinatal asphyxia

BACKGROUND

Perinatal asphyxia often leads to hypoxic-ischemic encephalopathy (HIE) with a high risk of neurodevelopmental consequences. While moderate and severe HIE link to high morbidity, less is known about brain effects of perinatal asphyxia with no or only mild HIE. Here, we test the hypothesis that cortical activity networks in the newborn infants show a dose-response to asphyxia.

METHODS

We performed EEG recordings for infants with perinatal asphyxia/HIE of varying severity (n = 52) and controls (n = 53) and examined well-established computational metrics of cortical network activity.

RESULTS

We found graded alterations in cortical activity networks according to severity of asphyxia/HIE. Furthermore, our findings correlated with early clinical recovery measured by the time to attain full oral feeding.

CONCLUSION

We show that both local and large-scale correlated cortical activity are affected by increasing severity of HIE after perinatal asphyxia, suggesting that HIE and perinatal asphyxia are better represented as a continuum rather than the currently used discreet categories. These findings imply that automated computational measures of cortical function may be useful in characterizing the dose effects of adversity in the neonatal brain; such metrics hold promise for benchmarking clinical trials via patient stratification or as early outcome measures.

IMPACT

Perinatal asphyxia causes every fourth neonatal death worldwide and provides a diagnostic and prognostic challenge for the clinician. We report that infants with perinatal asphyxia show specific graded responses in cortical networks according to severity of asphyxia and ensuing hypoxic-ischaemic encephalopathy. Early EEG recording and automated computational measures of brain function have potential to help in clinical evaluation of infants with perinatal asphyxia.

Virtual lesions in MEG reveal increasing vulnerability of the language network from early childhood through adolescence

In childhood, language outcomes following brain injury are inversely related to age. Neuroimaging findings suggest that extensive representation and/or topological redundancy may confer the pediatric advantage. Here, we assess whole brain and language network resilience using in silico attacks, for 85 children participating in a magnetoencephalography (MEG) study. Nodes are targeted based on eigenvector centrality, betweenness centrality, or at random. The size of each connected component is assessed after iterated node removal; the percolation point, or moment of dis-integration, is defined as the first instance where the second largest component peaks in size. To overcome known effects of fixed thresholding on subsequent graph and resilience analyses, we study percolation across all possible network densities, within a Functional Data Analysis (FDA) framework. We observe age-related increases in vulnerability for random and betweenness centrality-based attacks for whole-brain and stories networks (adjusted-p < 0.05). Here we show that changes in topology underlie increasing language network vulnerability in development.

Brain-to-brain communication during musical improvisation: a performance case study

Understanding and predicting others' actions in ecological settings is an important research goal in social neuroscience. Here, we deployed a mobile brain-body imaging (MoBI) methodology to analyze inter-brain communication between professional musicians during a live jazz performance. Specifically, bispectral analysis was conducted to assess the synchronization of scalp electroencephalographic (EEG) signals from three expert musicians during a three-part 45 minute jazz performance, during which a new musician joined every five minutes. The bispectrum was estimated for all musician dyads, electrode combinations, and five frequency bands. The results showed higher bispectrum in the beta and gamma frequency bands (13-50 Hz) when more musicians performed together, and when they played a musical phrase synchronously. Positive bispectrum amplitude changes were found approximately three seconds prior to the identified synchronized performance events suggesting preparatory cortical activity predictive of concerted behavioral action. Moreover, a higher amount of synchronized EEG activity, across electrode regions, was observed as more musicians performed, with inter-brain synchronization between the temporal, parietal, and occipital regions the most frequent. Increased synchrony between the musicians' brain activity reflects shared multi-sensory processing and movement intention in a musical improvisation task.

Deep learning based source imaging provides strong sublobar localization of epileptogenic zone from MEG interictal spikes

Electromagnetic source imaging (ESI) offers unique capability of imaging brain dynamics for studying brain functions and aiding the clinical management of brain disorders. Challenges exist in ESI due to the ill-posedness of the inverse problem and thus the need of modeling the underlying brain dynamics for regularizations. Advances in generative models provide opportunities for more accurate and realistic source modeling that could offer an alternative approach to ESI for modeling the underlying brain dynamics beyond equivalent physical source models. However, it is not straightforward to explicitly formulate the knowledge arising from these generative models within the conventional ESI framework. Here we investigate a novel source imaging framework based on mesoscale neuronal modeling and deep learning (DL) that can learn the sensor-source mapping relationship directly from MEG data for ESI. Two DL-based ESI models were trained based on data generated by neural mass models and either generic or personalized head models. The robustness of the two DL models was evaluated by systematic computer simulations and clinical validation in a cohort of 29 drug-resistant focal epilepsy patients who underwent intracranial EEG (iEEG) evaluation or surgical resection. Results estimated from pre-operative MEG interictal spikes were quantified using the overlap with resection regions and the distance to the seizure-onset zone (SOZ) defined by iEEG recordings. The DL-based ESI provided robust results when no personalized head geometry is considered, reaching a spatial dispersion of 21.90 ± 19.03 mm, sublobar concordance of 83 %, and sublobar sensitivity and specificity of 66 and 97 % respectively. When using personalized head geometry derived from individual patients' MRI in the training data, personalized DL-based ESI model can further improve the performance and reached a spatial dispersion of 8.19 ± 8.14 mm, sublobar concordance of 93 %, and sublobar sensitivity and specificity of 77 and 99 % respectively. When compared to the SOZ, the localization error of the personalized approach is 15.78 ± 5.54 mm, outperforming the conventional benchmarks. This work demonstrates that combining generative models and deep learning enables an accurate and robust imaging of epileptogenic zone from MEG recordings with strong sublobar precision, suggesting its added value to enhancing MEG source localization and imaging, and to epilepsy source localization and other clinical applications.

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.

Graph neural networks in EEG spike detection

OBJECTIVE

This study develops new machine learning architectures that are more adept at detecting interictal epileptiform discharges (IEDs) in scalp EEG. A comparison of results using the average precision (AP) metric is made with the proposed models on two datasets obtained from Baptist Hospital of Miami and Temple University Hospital.

METHODS

Applying graph neural networks (GNNs) on functional connectivity (FC) maps of different frequency sub-bands to yield a novel architecture we call FC-GNN. Attention mechanism is applied on a complete graph to let the neural network select its important edges, hence bypassing the extraction of features, a model we refer to as CA-GNN.

RESULTS

On the Baptist Hospital dataset, the results were as follows: Vanilla Self-Attention →0.9029±0.0431, Hierarchical Attention →0.8546±0.0587, Vanilla Visual Geometry Group (VGG) →0.92±0.0618, Satelight →0.9219±0.046, FC-GNN →0.9731±0.0187, and CA-GNN →0.9788±0.0125. In the same order, the results on the Temple University Hospital dataset are 0.9692, 0.9113, 0.97, 0.9575, 0.963, and 0.9879.

CONCLUSION

Based on the good results they yield, GNNs prove to have a strong potential in detecting epileptogenic activity.

SIGNIFICANCE

This study opens the door for the discovery of the powerful role played by GNNs in capturing IEDs, which is an essential step for identifying the epileptogenic networks of the affected brain and hence improving the prospects for more accurate 3D source localization.

Enhanced behavioral performance through interareal gamma and beta synchronization

Cognitive functioning requires coordination between brain areas. Between visual areas, feedforward gamma synchronization improves behavioral performance. Here, we investigate whether similar principles hold across brain regions and frequency bands, using simultaneous electrocorticographic recordings from 15 areas of two macaque monkeys during performance of a selective attention task. Short behavioral reaction times (RTs), suggesting efficient interareal communication, occurred when occipital areas V1, V2, V4, and DP showed gamma synchronization, and fronto-central areas S1, 5, F1, F2, and F4 showed beta synchronization. For both area clusters and corresponding frequency bands, deviations from the typically observed phase relations increased RTs. Across clusters and frequency bands, good phase relations occurred in a correlated manner specifically when they processed the behaviorally relevant stimulus. Furthermore, the fronto-central cluster exerted a beta-band influence onto the occipital cluster whose strength predicted short RTs. These results suggest that local gamma and beta synchronization and their inter-regional coordination jointly improve behavioral performance.

Neonatal cortical activity organizes into transient network states that are affected by vigilance states and brain injury

Early neurodevelopment is critically dependent on the structure and dynamics of spontaneous neuronal activity; however, the natural organization of newborn cortical networks is poorly understood. Recent adult studies suggest that spontaneous cortical activity exhibits discrete network states with physiological correlates. Here, we studied newborn cortical activity during sleep using hidden Markov modeling to determine the presence of such discrete neonatal cortical states (NCS) in 107 newborn infants, with 47 of them presenting with a perinatal brain injury. Our results show that neonatal cortical activity organizes into four discrete NCSs that are present in both cardinal sleep states of a newborn infant, active and quiet sleep, respectively. These NCSs exhibit state-specific spectral and functional network characteristics. The sleep states exhibit different NCS dynamics, with quiet sleep presenting higher fronto-temporal activity and a stronger brain-wide neuronal coupling. Brain injury was associated with prolonged lifetimes of the transient NCSs, suggesting lowered dynamics, or flexibility, in the cortical networks. Taken together, the findings suggest that spontaneously occurring transient network states are already present at birth, with significant physiological and pathological correlates; this NCS analysis framework can be fully automatized, and it holds promise for offering an objective, global level measure of early brain function for benchmarking neurodevelopmental or clinical research.

The influence of inter-regional delays in generating large-scale brain networks of phase synchronization

Large-scale networks of phase synchronization are considered to regulate the communication between brain regions fundamental to cognitive function, but the mapping to their structural substrates, i.e., the structure-function relationship, remains poorly understood. Biophysical Network Models (BNMs) have demonstrated the influences of local oscillatory activity and inter-regional anatomical connections in generating alpha-band (8-12 Hz) networks of phase synchronization observed with Electroencephalography (EEG) and Magnetoencephalography (MEG). Yet, the influence of inter-regional conduction delays remains unknown. In this study, we compared a BNM with standard "distance-dependent delays", which assumes constant conduction velocity, to BNMs with delays specified by two alternative methods accounting for spatially varying conduction velocities, "isochronous delays" and "mixed delays". We followed the Approximate Bayesian Computation (ABC) workflow, i) specifying neurophysiologically informed prior distributions of BNM parameters, ii) verifying the suitability of the prior distributions with Prior Predictive Checks, iii) fitting each of the three BNMs to alpha-band MEG resting-state data (N = 75) with Bayesian optimization for Likelihood-Free Inference (BOLFI), and iv) choosing between the fitted BNMs with ABC model comparison on a separate MEG dataset (N = 30). Prior Predictive Checks revealed the range of dynamics generated by each of the BNMs to encompass those seen in the MEG data, suggesting the suitability of the prior distributions. Fitting the models to MEG data yielded reliable posterior distributions of the parameters of each of the BNMs. Finally, model comparison revealed the BNM with "distance-dependent delays", as the most probable to describe the generation of alpha-band networks of phase synchronization seen in MEG. These findings suggest that distance-dependent delays might contribute to the neocortical architecture of human alpha-band networks of phase synchronization. Hence, our study illuminates the role of inter-regional delays in generating the large-scale networks of phase synchronization that might subserve the communication between regions vital to cognition.

Identifying good practices for detecting inter-regional linear functional connectivity from EEG

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

Brain criticality predicts individual levels of inter-areal synchronization in human electrophysiological data

Neuronal oscillations and their synchronization between brain areas are fundamental for healthy brain function. Yet, synchronization levels exhibit large inter-individual variability that is associated with behavioral variability. We test whether individual synchronization levels are predicted by individual brain states along an extended regime of critical-like dynamics - the Griffiths phase (GP). We use computational modelling to assess how synchronization is dependent on brain criticality indexed by long-range temporal correlations (LRTCs). We analyze LRTCs and synchronization of oscillations from resting-state magnetoencephalography and stereo-electroencephalography data. Synchronization and LRTCs are both positively linearly and quadratically correlated among healthy subjects, while in epileptogenic areas they are negatively linearly correlated. These results show that variability in synchronization levels is explained by the individual position along the GP with healthy brain areas operating in its subcritical and epileptogenic areas in its supercritical side. We suggest that the GP is fundamental for brain function allowing individual variability while retaining functional advantages of criticality.

Structural basis of envelope and phase intrinsic coupling modes in the cerebral cortex

Intrinsic coupling modes (ICMs) can be observed in ongoing brain activity at multiple spatial and temporal scales. Two families of ICMs can be distinguished: phase and envelope ICMs. The principles that shape these ICMs remain partly elusive, in particular their relation to the underlying brain structure. Here we explored structure-function relationships in the ferret brain between ICMs quantified from ongoing brain activity recorded with chronically implanted micro-ECoG arrays and structural connectivity (SC) obtained from high-resolution diffusion MRI tractography. Large-scale computational models were used to explore the ability to predict both types of ICMs. Importantly, all investigations were conducted with ICM measures that are sensitive or insensitive to volume conduction effects. The results show that both types of ICMs are significantly related to SC, except for phase ICMs when using measures removing zero-lag coupling. The correlation between SC and ICMs increases with increasing frequency which is accompanied by reduced delays. Computational models produced results that were highly dependent on the specific parameter settings. The most consistent predictions were derived from measures solely based on SC. Overall, the results demonstrate that patterns of cortical functional coupling as reflected in both phase and envelope ICMs are both related, albeit to different degrees, to the underlying structural connectivity in the cerebral cortex.

Identifying oscillatory brain networks with hidden Gaussian graphical spectral models of MEEG

Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflection of functional brain network activity with a temporal resolution that no other in vivo neuroimage may provide. A second step estimating functional connectivity from such activity pseudodata unveil the oscillatory brain networks that strongly correlate with all cognition and behavior. Simulations of such MEG or EEG inverse problem also reveal estimation errors of the functional connectivity determined by any of the state-of-the-art inverse solutions. We disclose a significant cause of estimation errors originating from misspecification of the functional network model incorporated into either inverse solution steps. We introduce the Bayesian identification of a Hidden Gaussian Graphical Spectral (HIGGS) model specifying such oscillatory brain networks model. In human EEG alpha rhythm simulations, the estimation errors measured as ROC performance do not surpass 2% in our HIGGS inverse solution and reach 20% in state-of-the-art methods. Macaque simultaneous EEG/ECoG recordings provide experimental confirmation for our results with 1/3 times larger congruence according to Riemannian distances than state-of-the-art methods.

Sensor array design of optically pumped magnetometers for accurately estimating source currents

An optically pumped magnetometer (OPM) is a new generation of magnetoencephalography (MEG) devices that is small, light, and works at room temperature. Due to these characteristics, OPMs enable flexible and wearable MEG systems. On the other hand, if we have a limited number of OPM sensors, we need to carefully design their sensor arrays depending on our purposes and regions of interests (ROIs). In this study, we propose a method that designs OPM sensor arrays for accurately estimating the cortical currents at the ROIs. Based on the resolution matrix of minimum norm estimate (MNE), our method sequentially determines the position of each sensor to optimize its inverse filter pointing to the ROIs and suppressing the signal leakage from the other areas. We call this method the Sensor array Optimization based on Resolution Matrix (SORM). We conducted simple and realistic simulation tests to evaluate its characteristics and efficacy for real OPM-MEG data. SORM designed the sensor arrays so that their leadfield matrices had high effective ranks as well as high sensitivities to ROIs. Although SORM is based on MNE, the sensor arrays designed by SORM were effective not only when we estimated the cortical currents by MNE but also when we did so by other methods. With real OPM-MEG data we confirmed its validity for real data. These analyses suggest that SORM is especially useful when we want to accurately estimate ROIs' activities with a limited number of OPM sensors, such as brain-machine interfaces and diagnosing brain diseases.

Modules in connectomes of phase-synchronization comprise anatomically contiguous, functionally related regions

Modules in brain functional connectomes are essential to balancing segregation and integration of neuronal activity. Connectomes are the complete set of pairwise connections between brain regions. Non-invasive Electroencephalography (EEG) and Magnetoencephalography (MEG) have been used to identify modules in connectomes of phase-synchronization. However, their resolution is suboptimal because of spurious phase-synchronization due to EEG volume conduction or MEG field spread. Here, we used invasive, intracerebral recordings from stereo-electroencephalography (SEEG, N = 67), to identify modules in connectomes of phase-synchronization. To generate SEEG-based group-level connectomes affected only minimally by volume conduction, we used submillimeter accurate localization of SEEG contacts and referenced electrode contacts in cortical gray matter to their closest contacts in white matter. Combining community detection methods with consensus clustering, we found that the connectomes of phase-synchronization were characterized by distinct and stable modules at multiple spatial scales, across frequencies from 3 to 320 Hz. These modules were highly similar within canonical frequency bands. Unlike the distributed brain systems identified with functional Magnetic Resonance Imaging (fMRI), modules up to the high-gamma frequency band comprised only anatomically contiguous regions. Notably, the identified modules comprised cortical regions involved in shared repertoires of sensorimotor and cognitive functions including memory, language and attention. These results suggest that the identified modules represent functionally specialised brain systems, which only partially overlap with the brain systems reported with fMRI. Hence, these modules might regulate the balance between functional segregation and functional integration through phase-synchronization.

Time-Lagged Multidimensional Pattern Connectivity (TL-MDPC): An EEG/MEG pattern transformation based functional connectivity metric

Functional and effective connectivity methods are essential to study the complex information flow in brain networks underlying human cognition. Only recently have connectivity methods begun to emerge that make use of the full multidimensional information contained in patterns of brain activation, rather than unidimensional summary measures of these patterns. To date, these methods have mostly been applied to fMRI data, and no method allows vertex-to-vertex transformations with the temporal specificity of EEG/MEG data. Here, we introduce time-lagged multidimensional pattern connectivity (TL-MDPC) as a novel bivariate functional connectivity metric for EEG/MEG research. TL-MDPC estimates the vertex-to-vertex transformations among multiple brain regions and across different latency ranges. It determines how well patterns in ROI X at time point tx can linearly predict patterns of ROI Y at time point ty. In the present study, we use simulations to demonstrate TL-MDPC's increased sensitivity to multidimensional effects compared to a unidimensional approach across realistic choices of number of trials and signal-to-noise ratios. We applied TL-MDPC, as well as its unidimensional counterpart, to an existing dataset varying the depth of semantic processing of visually presented words by contrasting a semantic decision and a lexical decision task. TL-MDPC detected significant effects beginning very early on, and showed stronger task modulations than the unidimensional approach, suggesting that it is capable of capturing more information. With TL-MDPC only, we observed rich connectivity between core semantic representation (left and right anterior temporal lobes) and semantic control (inferior frontal gyrus and posterior temporal cortex) areas with greater semantic demands. TL-MDPC is a promising approach to identify multidimensional connectivity patterns, typically missed by unidimensional approaches.

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.

Minimizing the distortions in electrophysiological source imaging of cortical oscillatory activity via Spectral Structured Sparse Bayesian Learning

Oscillatory processes at all spatial scales and on all frequencies underpin brain function. Electrophysiological Source Imaging (ESI) is the data-driven brain imaging modality that provides the inverse solutions to the source processes of the EEG, MEG, or ECoG data. This study aimed to carry out an ESI of the source cross-spectrum while controlling common distortions of the estimates. As with all ESI-related problems under realistic settings, the main obstacle we faced is a severely ill-conditioned and high-dimensional inverse problem. Therefore, we opted for Bayesian inverse solutions that posited a priori probabilities on the source process. Indeed, rigorously specifying both the likelihoods and a priori probabilities of the problem leads to the proper Bayesian inverse problem of cross-spectral matrices. These inverse solutions are our formal definition for cross-spectral ESI (cESI), which requires a priori of the source cross-spectrum to counter the severe ill-condition and high-dimensionality of matrices. However, inverse solutions for this problem were NP-hard to tackle or approximated within iterations with bad-conditioned matrices in the standard ESI setup. We introduce cESI with a joint a priori probability upon the source cross-spectrum to avoid these problems. cESI inverse solutions are low-dimensional ones for the set of random vector instances and not random matrices. We achieved cESI inverse solutions through the variational approximations via our Spectral Structured Sparse Bayesian Learning (ssSBL) algorithm https://github.com/CCC-members/Spectral-Structured-Sparse-Bayesian-Learning. We compared low-density EEG (10-20 system) ssSBL inverse solutions with reference cESIs for two experiments: (a) high-density MEG that were used to simulate EEG and (b) high-density macaque ECoG that were recorded simultaneously with EEG. The ssSBL resulted in two orders of magnitude with less distortion than the state-of-the-art ESI methods. Our cESI toolbox, including the ssSBL method, is available at https://github.com/CCC-members/BC-VARETA_Toolbox.

Enhanced Behavioral Performance through Interareal Gamma and Beta Synchronization

Cognitive functioning requires coordination between brain areas. Between visual areas, feedforward gamma synchronization improves behavioral performance. Here, we investigate whether similar principles hold across brain regions and frequency bands, using simultaneous local field potential recordings from 15 areas during performance of a selective attention task. Short behavioral reaction times (RTs), an index of efficient interareal communication, occurred when occipital areas V1, V2, V4, DP showed gamma synchronization, and fronto-central areas S1, 5, F1, F2, F4 showed beta synchronization. For both area clusters and corresponding frequency bands, deviations from the typically observed phase relations increased RTs. Across clusters and frequency bands, good phase relations occurred in a correlated manner specifically when they processed the behaviorally relevant stimulus. Furthermore, the fronto-central cluster exerted a beta-band influence onto the occipital cluster whose strength predicted short RTs. These results suggest that local gamma and beta synchronization and their inter-regional coordination jointly improve behavioral performance.

Communication dynamics in the human connectome shape the cortex-wide propagation of direct electrical stimulation

Communication between gray matter regions underpins all facets of brain function. We study inter-areal communication in the human brain using intracranial EEG recordings, acquired following 29,055 single-pulse direct electrical stimulations in a total of 550 individuals across 20 medical centers (average of 87 ± 37 electrode contacts per subject). We found that network communication models-computed on structural connectivity inferred from diffusion MRI-can explain the causal propagation of focal stimuli, measured at millisecond timescales. Building on this finding, we show that a parsimonious statistical model comprising structural, functional, and spatial factors can accurately and robustly predict cortex-wide effects of brain stimulation (R²=46% in data from held-out medical centers). Our work contributes toward the biological validation of concepts in network neuroscience and provides insight into how connectome topology shapes polysynaptic inter-areal signaling. We anticipate that our findings will have implications for research on neural communication and the design of brain stimulation paradigms.

Women and men have a similar potential for malevolent creativity - But their underlying brain mechanisms are different

Research interest in gender differences in aggression and creative ideation cumulates in the phenomenon of malevolent creativity. Taking another critical step in understanding malevolent creativity, we investigated gender differences in brain activation and functional coupling of cortical sites in the EEG alpha band while n = 88 women and men purposefully generated malevolent creative ideas for taking revenge on others. Results showed that malevolent creativity performance between the genders was similar; however, their underlying EEG patterns were markedly different. While women exhibited a steep decrease of task-related alpha power from frontal to left central-temporal, men's malevolent creative ideation was characterized by a more diffuse pattern of task-related alpha power changes, along with decreased frontal-central coupling. Per interpretation, women's malevolent creative thinking may more strongly rely on controlled semantic memory retrieval and novel re-combination of social/relationship information, while men may utilize more automatic motor-related imagery that may predominantly facilitate physical revenge ideation. Our findings add novel evidence to the idea that women and men engage different neurocognitive strategies to achieve similar creative performance and may help to further illuminate the darker side of creative ideation.

Multistage classification identifies altered cortical phase- and amplitude-coupling in Multiple Sclerosis

Distinguishing groups of subjects or experimental conditions in a high-dimensional feature space is a common goal in modern neuroimaging studies. Successful classification depends on the selection of relevant features as not every neuronal signal component or parameter is informative about the research question at hand. Here, we developed a novel unsupervised multistage analysis approach that combines dimensionality reduction, bootstrap aggregating and multivariate classification to select relevant neuronal features. We tested the approach by identifying changes of brain-wide electrophysiological coupling in Multiple Sclerosis. Multiple Sclerosis is a demyelinating disease of the central nervous system that can result in cognitive decline and physical disability. However, related changes in large-scale brain interactions remain poorly understood and corresponding non-invasive biomarkers are sparse. We thus compared brain-wide phase- and amplitude-coupling of frequency specific neuronal activity in relapsing-remitting Multiple Sclerosis patients (n = 17) and healthy controls (n = 17) using magnetoencephalography. Changes in this dataset included both, increased and decreased phase- and amplitude-coupling in wide-spread, bilateral neuronal networks across a broad range of frequencies. These changes allowed to successfully classify patients and controls with an accuracy of 84%. Furthermore, classification confidence predicted behavioral scores of disease severity. In sum, our results unravel systematic changes of large-scale phase- and amplitude coupling in Multiple Sclerosis. Furthermore, our results establish a new analysis approach to efficiently contrast high-dimensional neuroimaging data between experimental groups or conditions.

Structural connectome constrained graphical lasso for MEG partial coherence

Structural connectivity provides the backbone for communication between neural populations. Since axonal transmission occurs on a millisecond time scale, measures of M/EEG functional connectivity sensitive to phase synchronization, such as coherence, are expected to reflect structural connectivity. We develop a model of MEG functional connectivity whose edges are constrained by the structural connectome. The edge strengths are defined by partial coherence, a measure of conditional dependence. We build a new method-the adaptive graphical lasso (AGL)-to fit the partial coherence to perform inference on the hypothesis that the structural connectome is reflected in MEG functional connectivity. In simulations, we demonstrate that the structural connectivity's influence on the partial coherence can be inferred using the AGL. Further, we show that fitting the partial coherence is superior to alternative methods at recovering the structural connectome, even after the source localization estimates required to map MEG from sensors to the cortex. Finally, we show how partial coherence can be used to explore how distinct parts of the structural connectome contribute to MEG functional connectivity in different frequency bands. Partial coherence offers better estimates of the strength of direct functional connections and consequently a potentially better estimate of network structure.

Facilitating early parent-infant emotional connection improves cortical networks in preterm infants

Exposure to environmental adversities during early brain development, such as preterm birth, can affect early brain organization. Here, we studied whether development of cortical activity networks in preterm infants may be improved by a multimodal environmental enrichment via bedside facilitation of mother-infant emotional connection. We examined functional cortico-cortical connectivity at term age using high-density electroencephalography recordings in infants participating in a randomized controlled trial of Family Nurture Intervention (FNI). Our results identify several large-scale, frequency-specific network effects of FNI, most extensively in the alpha frequency in fronto-central cortical regions. The connectivity strength in this network was correlated to later neurocognitive performance, and it was comparable to healthy term-born infants rather than the infants receiving standard care. These findings suggest that preterm neurodevelopmental care can be improved by a biologically driven environmental enrichment, such as early facilitation of direct human connection.

Genetic polymorphisms in COMT and BDNF influence synchronization dynamics of human neuronal oscillations

Neuronal oscillations, their inter-areal synchronization, and scale-free dynamics constitute fundamental mechanisms for cognition by regulating communication in neuronal networks. These oscillatory dynamics have large inter-individual variability that is partly heritable. We hypothesized that this variability could be partially explained by genetic polymorphisms in neuromodulatory genes. We recorded resting-state magnetoencephalography (MEG) from 82 healthy participants and investigated whether oscillation dynamics were influenced by genetic polymorphisms in catechol-O-methyltransferase (COMT) Val¹⁵⁸Met and brain-derived neurotrophic factor (BDNF) Val⁶⁶Met. Both COMT and BDNF polymorphisms influenced local oscillation amplitudes and their long-range temporal correlations (LRTCs), while only BDNF polymorphism affected the strength of large-scale synchronization. Our findings demonstrate that COMT and BDNF genetic polymorphisms contribute to inter-individual variability in neuronal oscillation dynamics. Comparison of these results to computational modeling of near-critical synchronization dynamics further suggested that COMT and BDNF polymorphisms influenced local oscillations by modulating the excitation-inhibition balance according to the brain criticality framework.

•

Human local oscillation dynamics is influenced by polymorphisms in COMT and BNDF

•

COMT and BDNF influence oscillation amplitudes and long-range temporal correlations

•

BDNF polymorphism affected the strength of large-scale synchronization

•

Framework of brain criticality links COMT and BDNF with local E/I-balance

Variations of Resting-State EEG-Based Functional Networks in Brain Maturation From Early Childhood to Adolescence

Background and Purpose

Alterations in human brain functional networks with maturation have been explored extensively in numerous electroencephalography (EEG) and functional magnetic resonance imaging studies. It is known that the age-related changes in the functional networks occurring prior to adulthood deviate from ordinary trajectories of network-based brain maturation across the adult lifespan.

Methods

This study investigated the longitudinal evolution of resting-state EEG-based functional networks from early childhood to adolescence among 212 pediatric patients (age 12.2±3.5 years, range 4.4–17.9) in 6 frequency bands using 8 types of functional connectivity measures in the amplitude, frequency, and phase domains.

Results

Electrophysiological aspects of network-based pediatric brain maturation were characterized by increases in both functional segregation and integration up to middle adolescence. EEG oscillations in the upper alpha band reflected the age-related increases in mean node strengths and mean clustering coefficients and a decrease in the characteristic path lengths better than did those in the other frequency bands, especially for the phase-domain functional connectivity. The frequency-band-specific age-related changes in the global network metrics were influenced more by volume-conduction effects than by the domain specificity of the functional connectivity measures.

Conclusions

We believe that this is the first study to reveal EEG-based functional network properties during preadult brain maturation based on various functional connectivity measures. The findings potentially have clinical applications in the diagnosis and treatment of age-related brain disorders.

More creative ideas are associated with increased right posterior power and frontal-parietal/occipital coupling in the upper alpha band: A within-subjects study

The neurophysiological investigation of creative idea generation is a growing research area. EEG studies congruently reported the sensitivity of upper alpha power (10-12 Hz) for the creative ideation process and its outcome. However, the majority of studies were between-subject design studies and research directly comparing the neurophysiological activation pattern when generating more and less creative ideas within a person are rare. Therefore, the present study was specifically focused on investigating brain activation patterns associated with the generation of more vs. less creative ideas. We applied an alternate uses task (AU-task; i.e., finding original uses for everyday objects such as a brick) in a sample of 74 participants and recorded the brain activation during the AU-task and reference period. A portable EEG system with 21 dry electrodes arranged in the international 10-20 system and linked ear as reference was used. We found a higher increase of upper alpha power during creative ideation (relative to reference period, i.e., task-related power, TRP) over right posterior sites when people generated more compared to less creative ideas. This was accompanied by an increase of functional coupling (i.e., task-related coherence increase) between frontal and parietal/occipital sites, which suggests higher internal attention and more control over sensory processes. Taken together, these findings complement the existing creativity research literature and indicate the importance of alpha power for the creative ideation process also within people.

Mixtures of large-scale dynamic functional brain network modes

Accurate temporal modelling of functional brain networks is essential in the quest for understanding how such networks facilitate cognition. Researchers are beginning to adopt time-varying analyses for electrophysiological data that capture highly dynamic processes on the order of milliseconds. Typically, these approaches, such as clustering of functional connectivity profiles and Hidden Markov Modelling (HMM), assume mutual exclusivity of networks over time. Whilst a powerful constraint, this assumption may be compromising the ability of these approaches to describe the data effectively. Here, we propose a new generative model for functional connectivity as a time-varying linear mixture of spatially distributed statistical "modes". The temporal evolution of this mixture is governed by a recurrent neural network, which enables the model to generate data with a rich temporal structure. We use a Bayesian framework known as amortised variational inference to learn model parameters from observed data. We call the approach DyNeMo (for Dynamic Network Modes), and show using simulations it outperforms the HMM when the assumption of mutual exclusivity is violated. In resting-state MEG, DyNeMo reveals a mixture of modes that activate on fast time scales of 100-150 ms, which is similar to state lifetimes found using an HMM. In task MEG data, DyNeMo finds modes with plausible, task-dependent evoked responses without any knowledge of the task timings. Overall, DyNeMo provides decompositions that are an approximate remapping of the HMM's while showing improvements in overall explanatory power. However, the magnitude of the improvements suggests that the HMM's assumption of mutual exclusivity can be reasonable in practice. Nonetheless, DyNeMo provides a flexible framework for implementing and assessing future modelling developments.

Good scientific practice in EEG and MEG research: Progress and perspectives

Good scientific practice (GSP) refers to both explicit and implicit rules, recommendations, and guidelines that help scientists to produce work that is of the highest quality at any given time, and to efficiently share that work with the community for further scrutiny or utilization. For experimental research using magneto- and electroencephalography (MEEG), GSP includes specific standards and guidelines for technical competence, which are periodically updated and adapted to new findings. However, GSP also needs to be regularly revisited in a broader light. At the LiveMEEG 2020 conference, a reflection on GSP was fostered that included explicitly documented guidelines and technical advances, but also emphasized intangible GSP: a general awareness of personal, organizational, and societal realities and how they can influence MEEG research. This article provides an extensive report on most of the LiveMEEG contributions and new literature, with the additional aim to synthesize ongoing cultural changes in GSP. It first covers GSP with respect to cognitive biases and logical fallacies, pre-registration as a tool to avoid those and other early pitfalls, and a number of resources to enable collaborative and reproducible research as a general approach to minimize misconceptions. Second, it covers GSP with respect to data acquisition, analysis, reporting, and sharing, including new tools and frameworks to support collaborative work. Finally, GSP is considered in light of ethical implications of MEEG research and the resulting responsibility that scientists have to engage with societal challenges. Considering among other things the benefits of peer review and open access at all stages, the need to coordinate larger international projects, the complexity of MEEG subject matter, and today's prioritization of fairness, privacy, and the environment, we find that current GSP tends to favor collective and cooperative work, for both scientific and for societal reasons.

Periventricular nodular heterotopia is coupled with the neocortex during resting and task states

Periventricular nodular heterotopia (PVNH) is a well-defined developmental disorder characterized by failed neuronal migration, which forms ectopic neuronal nodules along the ventricular walls. Previous studies mainly focus on clinical symptoms caused by the PVNH tissue, such as seizures. However, little is known about whether and how neurons in the PVNH tissue functionally communicate with neurons in the neocortex. To probe this, we applied magnetoencephalography (MEG) and stereo-electroencephalography (sEEG) recordings to patients with PVNH during resting and task states. By estimating frequency-resolved phase coupling strength of the source-reconstructed neural activities, we found that the PVNH tissue was spontaneously coupled with the neocortex in the α-β frequency range, which was consistent with the synchronization pattern within the neocortical network. Furthermore, the coupling strength between PVNH and sensory areas effectively modulated the local neural activity in sensory areas. In both MEG and sEEG visual experiments, the PVNH tissue exhibited visual-evoked responses, with a similar pattern and latency as the ipsilateral visual cortex. These findings demonstrate that PVNH is functionally integrated into cognition-related cortical circuits, suggesting a co-development perspective of ectopic neurons after their migration failure.

EEG functional connectivity analysis in the source space

There is a growing interest in using electroencephalography (EEG) and source modeling to investigate functional interactions among cortical processes, particularly when dealing with pediatric populations. This paper introduces two pipelines that have been recently used to conduct EEG FC analysis in the cortical source space. The analytic streams of these pipelines can be summarized into the following steps: 1) cortical source reconstruction of high-density EEG data using realistic magnetic resonance imaging (MRI) models created with age-appropriate MRI templates; 2) segmentation of reconstructed source activities into brain regions of interest; and 3) estimation of FC in age-related frequency bands using robust EEG FC measures, such as weighted phase lag index and orthogonalized power envelope correlation. In this paper we demonstrate the two pipelines with resting-state EEG data collected from children at 12 and 36 months of age. We also discuss the advantages and limitations of the methods/techniques integrated into the pipelines. Given there is a need in the research community for open-access analytic toolkits that can be used for pediatric EEG data, programs and codes used for the current analysis are made available to the public.

NLGC: Network Localized Granger Causality with Application to MEG Directional Functional Connectivity Analysis

Identifying the directed connectivity that underlie networked activity between different cortical areas is critical for understanding the neural mechanisms behind sensory processing. Granger causality (GC) is widely used for this purpose in functional magnetic resonance imaging analysis, but there the temporal resolution is low, making it difficult to capture the millisecond-scale interactions underlying sensory processing. Magnetoencephalography (MEG) has millisecond resolution, but only provides low-dimensional sensor-level linear mixtures of neural sources, which makes GC inference challenging. Conventional methods proceed in two stages: First, cortical sources are estimated from MEG using a source localization technique, followed by GC inference among the estimated sources. However, the spatiotemporal biases in estimating sources propagate into the subsequent GC analysis stage, may result in both false alarms and missing true GC links. Here, we introduce the Network Localized Granger Causality (NLGC) inference paradigm, which models the source dynamics as latent sparse multivariate autoregressive processes and estimates their parameters directly from the MEG measurements, integrated with source localization, and employs the resulting parameter estimates to produce a precise statistical characterization of the detected GC links. We offer several theoretical and algorithmic innovations within NLGC and further examine its utility via comprehensive simulations and application to MEG data from an auditory task involving tone processing from both younger and older participants. Our simulation studies reveal that NLGC is markedly robust with respect to model mismatch, network size, and low signal-to-noise ratio, whereas the conventional two-stage methods result in high false alarms and mis-detections. We also demonstrate the advantages of NLGC in revealing the cortical network-level characterization of neural activity during tone processing and resting state by delineating task- and age-related connectivity changes.

Towards an objective evaluation of EEG/MEG source estimation methods - The linear approach

The spatial resolution of EEG/MEG source estimates, often described in terms of source leakage in the context of the inverse problem, poses constraints on the inferences that can be drawn from EEG/MEG source estimation results. Software packages for EEG/MEG data analysis offer a large choice of source estimation methods but few tools to experimental researchers for methods evaluation and comparison. Here, we describe a framework and tools for objective and intuitive resolution analysis of EEG/MEG source estimation based on linear systems analysis, and apply those to the most widely used distributed source estimation methods such as L2-minimum-norm estimation (L2-MNE) and linearly constrained minimum variance (LCMV) beamformers. Within this framework it is possible to define resolution metrics that define meaningful aspects of source estimation results (such as localization accuracy in terms of peak localization error, PLE, and spatial extent in terms of spatial deviation, SD) that are relevant to the task at hand and can easily be visualized. At the core of this framework is the resolution matrix, which describes the potential leakage from and into point sources (point-spread and cross-talk functions, or PSFs and CTFs, respectively). Importantly, for linear methods these functions allow generalizations to multiple sources or complex source distributions. This paper provides a tutorial-style introduction into linear EEG/MEG source estimation and resolution analysis aimed at experimental (rather than methods-oriented) researchers. We used this framework to demonstrate how L2-MNE-type as well as LCMV beamforming methods can be evaluated in practice using software tools that have only recently become available for routine use. Our novel methods comparison includes PLE and SD for a larger number of methods than in similar previous studies, such as unweighted, depth-weighted and normalized L2-MNE methods (including dSPM, sLORETA, eLORETA) and two LCMV beamformers. The results demonstrate that some methods can achieve low and even zero PLE for PSFs. However, their SD as well as both PLE and SD for CTFs are far less optimal for all methods, in particular for deep cortical areas. We hope that our paper will encourage EEG/MEG researchers to apply this approach to their own tasks at hand.

Auditory driven gamma synchrony is associated with cortical thickness in widespread cortical areas

OBJECTIVE

Gamma synchrony is a fundamental functional property of the cerebral cortex, impaired in multiple neuropsychiatric conditions (i.e. schizophrenia, Alzheimer's disease, stroke etc.). Auditory stimulation in the gamma range allows to drive gamma synchrony of the entire cortical mantle and to estimate the efficiency of the mechanisms sustaining it. As gamma synchrony depends strongly on the interplay between parvalbumin-positive interneurons and pyramidal neurons, we hypothesize an association between cortical thickness and gamma synchrony. To test this hypothesis, we employed a combined magnetoencephalography (MEG) - Magnetic Resonance Imaging (MRI) study.

METHODS

Cortical thickness was estimated from anatomical MRI scans. MEG measurements related to exposure of 40 Hz amplitude modulated tones were projected onto the cortical surface. Two measures of cortical synchrony were considered: (a) inter-trial phase consistency at 40 Hz, providing a vertex-wise estimation of gamma synchronization, and (b) phase-locking values between primary auditory cortices and whole cortical mantle, providing a measure of long-range cortical synchrony. A correlation between cortical thickness and synchronization measures was then calculated for 72 MRI-MEG scans.

RESULTS

Both inter-trial phase consistency and phase locking values showed a significant positive correlation with cortical thickness. For inter-trial phase consistency, clusters of strong associations were found in the temporal and frontal lobes, especially in the bilateral auditory and pre-motor cortices. Higher phase-locking values corresponded to higher cortical thickness in the frontal, temporal, occipital and parietal lobes.

DISCUSSION AND CONCLUSIONS

In healthy subjects, a thicker cortex corresponds to higher gamma synchrony and connectivity in the primary auditory cortex and beyond, likely reflecting underlying cell density involved in gamma circuitries. This result hints towards an involvement of gamma synchrony together with underlying brain structure in brain areas for higher order cognitive functions. This study contributes to the understanding of inherent cortical functional and structural brain properties, which might in turn constitute the basis for the definition of useful biomarkers in patients showing aberrant gamma synchronization.

Connectivity dynamics and cognitive variability during aging

Aging is associated with cognitive changes, with strong variations across individuals. One way to characterize this individual variability is to use techniques such as magnetoencephalography (MEG) to measure the dynamics of neural synchronization between brain regions, and the variability of this connectivity over time. Indeed, few studies have focused on fluctuations in the dynamics of brain networks over time and their evolution with age. We therefore characterize aging effects on MEG phase synchrony in healthy young and older adults from the Cam-CAN database. Age-related changes were observed, with an increase in the variability of brain synchronization, as well as a reversal of the direction of information transfer in the default mode network (DMN), in the delta frequency band. These changes in functional connectivity were associated with cognitive decline. Results suggest that advancing age is accompanied by a functional disorganization of dynamic networks, with a loss of communication stability and a decrease in the information transmitted.

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.

State Changes During Resting-State (Magneto)encephalographic Studies: The Effect of Drowsiness on Spectral, Connectivity, and Network Analyses

Background

A common problem in resting-state neuroimaging studies is that subjects become drowsy or fall asleep. Although this could drastically affect neurophysiological measurements, such as magnetoencephalography (MEG), its specific impact remains understudied. We aimed to systematically investigate how often drowsiness is present during resting-state MEG recordings, and how the state changes alter quantitative estimates of oscillatory activity, functional connectivity, and network topology.

Methods

About 8-min MEG recordings of 19 healthy subjects, split into ~13-s epochs, were scored for the presence of eyes-open (EO), alert eyes-closed (A-EC), or drowsy eyes-closed (D-EC) states. After projection to source-space, results of spectral, functional connectivity, and network analyses in 6 canonical frequency bands were compared between these states on a global and regional levels. Functional connectivity was analyzed using the phase lag index (PLI) and corrected amplitude envelope correlation (AECc), and network topology was analyzed using the minimum spanning tree (MST).

Results

Drowsiness was present in >55% of all epochs that did not fulfill the AASM criteria for sleep. There were clear differences in spectral results between the states (A-EC vs. D-EC) and conditions (EO vs. A-EC). The influence of state and condition was far less pronounced for connectivity analyses, with only minimal differences between D-EC and EO in the AECc in the delta band. There were no effects of drowsiness on any of the MST measures.

Conclusions

Drowsiness during eyes-closed resting-state MEG recordings is present in the majority of epochs, despite the instructions to stay awake. This has considerable influence on spectral properties, but much less so on functional connectivity and network topology. These findings are important for interpreting the results of EEG/MEG studies using spectral analyses in neurological disease, where recordings should be evaluated for the presence of drowsiness. For connectivity analyses or studies on network topology, this seems of far less importance.

Predicting time-resolved electrophysiological brain networks from structural eigenmodes

How temporal modulations in functional interactions are shaped by the underlying anatomical connections remains an open question. Here, we analyse the role of structural eigenmodes, in the formation and dissolution of temporally evolving functional brain networks using resting-state magnetoencephalography and diffusion magnetic resonance imaging data at the individual subject level. Our results show that even at short timescales, phase and amplitude connectivity can partly be expressed by structural eigenmodes, but hardly by direct structural connections. Albeit a stronger relationship was found between structural eigenmodes and time-resolved amplitude connectivity. Time-resolved connectivity for both phase and amplitude was mostly characterised by a stationary process, superimposed with very brief periods that showed deviations from this stationary process. For these brief periods, dynamic network states were extracted that showed different expressions of eigenmodes. Furthermore, the eigenmode expression was related to overall cognitive performance and co-occurred with fluctuations in community structure of functional networks. These results implicate that ongoing time-resolved resting-state networks, even at short timescales, can to some extent be understood in terms of activation and deactivation of structural eigenmodes and that these eigenmodes play a role in the dynamic integration and segregation of information across the cortex, subserving cognitive functions.

Ongoing Brain Activity and Its Role in Cognition: Dual versus Baseline Models

What is the role of the brain's ongoing activity for cognition? The predominant perspectives associate ongoing brain activity with resting state, the default-mode network (DMN), and internally oriented mentation. This triad is often contrasted with task states, non-DMN brain networks, and externally oriented mentation, together comprising a "dual model" of brain and cognition. In opposition to this duality, however, we propose that ongoing brain activity serves as a neuronal baseline; this builds upon Raichle's original search for the default mode of brain function that extended beyond the canonical default-mode brain regions. That entails what we refer to as the "baseline model." Akin to an internal biological clock for the rest of the organism, the ongoing brain activity may serve as an internal point of reference or standard by providing a shared neural code for the brain's rest as well as task states, including their associated cognition. Such shared neural code is manifest in the spatiotemporal organization of the brain's ongoing activity, including its global signal topography and dynamics like intrinsic neural timescales. We conclude that recent empirical evidence supports a baseline model over the dual model; the ongoing activity provides a global shared neural code that allows integrating the brain's rest and task states, its DMN and non-DMN, and internally and externally oriented cognition.