Adding dynamics to the Human Connectome Project with MEG

Spectral signature of attentional reorienting in the human brain

As we move in the environment, attention shifts to novel objects of interest based on either their sensory salience or behavioral value (reorienting). This study measures with magnetoencephalography (MEG) different properties (amplitude, onset-to-peak duration) of event-related desynchronization/synchronization (ERD/ERS) of oscillatory activity during a visuospatial attention task designed to separate activity related to reorienting vs. maintaining attention to the same location, controlling for target detection and response processes. The oscillatory activity was measured both in fMRI-defined regions of interest (ROIs) of the dorsal attention (DAN) and visual (VIS) networks, previously defined as task-relevant in the same subjects, or whole-brain in a pre-defined set of cortical ROIs encompassing the main brain networks. Reorienting attention (shift cues) as compared to maintaining attention (stay cues) produced a temporal sequence of ERD/ERS modulations at multiple frequencies in specific anatomical regions/networks. An early (∼330 ms), stronger, transient theta ERS occurred in task-relevant (DAN, VIS) and control networks (VAN, CON, FPN), possibly reflecting an alert/reset signal in response to the cue. A more sustained, behaviorally relevant, low-beta band ERD peaking ∼450 ms following shift cues (∼410 for stay cues) localized in frontal and parietal regions of the DAN. This modulation is consistent with a control signal re-routing information across visual hemifields. Contralateral vs. ipsilateral shift cues produced in occipital visual regions a stronger, sustained alpha ERD (peak ∼470 ms) and a longer, transient high beta/gamma ERS (peak ∼490 ms) related to preparatory visual modulations in advance of target occurrence. This is the first description of a cascade of oscillatory processes during attentional reorienting in specific anatomical regions and networks. Among these processes, a behaviorally relevant beta desynchronization in the FEF is likely associated with the control of attention shifts.

The Human Connectome Project: A retrospective

The Human Connectome Project (HCP) was launched in 2010 as an ambitious effort to accelerate advances in human neuroimaging, particularly for measures of brain connectivity; apply these advances to study a large number of healthy young adults; and freely share the data and tools with the scientific community. NIH awarded grants to two consortia; this retrospective focuses on the "WU-Minn-Ox" HCP consortium centered at Washington University, the University of Minnesota, and University of Oxford. In just over 6 years, the WU-Minn-Ox consortium succeeded in its core objectives by: 1) improving MR scanner hardware, pulse sequence design, and image reconstruction methods, 2) acquiring and analyzing multimodal MRI and MEG data of unprecedented quality together with behavioral measures from more than 1100 HCP participants, and 3) freely sharing the data (via the ConnectomeDB database) and associated analysis and visualization tools. To date, more than 27 Petabytes of data have been shared, and 1538 papers acknowledging HCP data use have been published. The "HCP-style" neuroimaging paradigm has emerged as a set of best-practice strategies for optimizing data acquisition and analysis. This article reviews the history of the HCP, including comments on key events and decisions associated with major project components. We discuss several scientific advances using HCP data, including improved cortical parcellations, analyses of connectivity based on functional and diffusion MRI, and analyses of brain-behavior relationships. We also touch upon our efforts to develop and share a variety of associated data processing and analysis tools along with detailed documentation, tutorials, and an educational course to train the next generation of neuroimagers. We conclude with a look forward at opportunities and challenges facing the human neuroimaging field from the perspective of the HCP consortium.

MEGnet: Automatic ICA-based artifact removal for MEG using spatiotemporal convolutional neural networks

Magnetoencephalography (MEG) is a functional neuroimaging tool that records the magnetic fields induced by neuronal activity; however, signal from non-neuronal sources can corrupt the data. Eye-blinks, saccades, and cardiac activity are three of the most common sources of non-neuronal artifacts. They can be measured by affixing eye proximal electrodes, as in electrooculography (EOG), and chest electrodes, as in electrocardiography (ECG), however this complicates imaging setup, decreases patient comfort, and can induce further artifacts from movement. This work proposes an EOG- and ECG-free approach to identify eye-blinks, saccades, and cardiac activity signals for automated artifact suppression. The contribution of this work is three-fold. First, using a data driven, multivariate decomposition approach based on Independent Component Analysis (ICA), a highly accurate artifact classifier is constructed as an amalgam of deep 1-D and 2-D Convolutional Neural Networks (CNNs) to automate the identification and removal of ubiquitous whole brain artifacts including eye-blink, saccade, and cardiac artifacts. The specific architecture of this network is optimized through an unbiased, computer-based hyperparameter random search. Second, visualization methods are applied to the learned abstraction to reveal what features the model uses and to bolster user confidence in the model's training and potential for generalization. Finally, the model is trained and tested on both resting-state and task MEG data from 217 subjects, and achieves a new state-of-the-art in artifact detection accuracy of 98.95% including 96.74% sensitivity and 99.34% specificity on the held out test-set. This work automates MEG processing for both clinical and research use, adapts to the acquired acquisition time, and can obviate the need for EOG or ECG electrodes for artifact detection.

Delineating between-subject heterogeneity in alpha networks with Spatio-Spectral Eigenmodes

•

A data-driven modal decomposition describes oscillations by their resonant frequency, damping time and network structure.

•

We show that the full multivariate transfer function can be rewritten as a linear superposition of these modes.

•

These modal coordinates factorise oscillatory systems without pre-specification of frequency bands or regions of interest.

•

Using these modes, we find a spatial gradient in alpha peak frequency between Occipital and Parietal cortex .

•

This gradient is highly variable between participants, showing shifts in spatial structure and peak frequency.

Between subject variability in the spatial and spectral structure of oscillatory networks can be highly informative but poses a considerable analytic challenge. Here, we describe a data-driven modal decomposition of a multivariate autoregressive model that simultaneously identifies oscillations by their peak frequency, damping time and network structure. We use this decomposition to define a set of Spatio-Spectral Eigenmodes (SSEs) providing a parsimonious description of oscillatory networks. We show that the multivariate system transfer function can be rewritten in these modal coordinates, and that the full transfer function is a linear superposition of all modes in the decomposition. The modal transfer function is a linear summation and therefore allows for single oscillatory signals to be isolated and analysed in terms of their spectral content, spatial distribution and network structure. We validate the method on simulated data and explore the structure of whole brain oscillatory networks in eyes-open resting state MEG data from the Human Connectome Project. We are able to show a wide between participant variability in peak frequency and network structure of alpha oscillations and show a distinction between occipital ’high-frequency alpha’ and parietal ’low-frequency alpha’. The frequency difference between occipital and parietal alpha components is present within individual participants but is partially masked by larger between subject variability; a 10Hz oscillation may represent the high-frequency occipital component in one participant and the low-frequency parietal component in another. This rich characterisation of individual neural phenotypes has the potential to enhance analyses into the relationship between neural dynamics and a person’s behavioural, cognitive or clinical state.

Exploring MEG brain fingerprints: Evaluation, pitfalls, and interpretations

Individual characterization of subjects based on their functional connectome (FC), termed "FC fingerprinting", has become a highly sought-after goal in contemporary neuroscience research. Recent functional magnetic resonance imaging (fMRI) studies have demonstrated unique characterization and accurate identification of individuals as an accomplished task. However, FC fingerprinting in magnetoencephalography (MEG) data is still widely unexplored. Here, we study resting-state MEG data from the Human Connectome Project to assess the MEG FC fingerprinting and its relationship with several factors including amplitude- and phase-coupling functional connectivity measures, spatial leakage correction, frequency bands, and behavioral significance. To this end, we first employ two identification scoring methods, differential identifiability and success rate, to provide quantitative fingerprint scores for each FC measurement. Secondly, we explore the edgewise and nodal MEG fingerprinting patterns across the different frequency bands (delta, theta, alpha, beta, and gamma). Finally, we investigate the cross-modality fingerprinting patterns obtained from MEG and fMRI recordings from the same subjects. We assess the behavioral significance of FC across connectivity measures and imaging modalities using partial least square correlation analyses. Our results suggest that fingerprinting performance is heavily dependent on the functional connectivity measure, frequency band, identification scoring method, and spatial leakage correction. We report higher MEG fingerprinting performances in phase-coupling methods, central frequency bands (alpha and beta), and in the visual, frontoparietal, dorsal-attention, and default-mode networks. Furthermore, cross-modality comparisons reveal a certain degree of spatial concordance in fingerprinting patterns between the MEG and fMRI data, especially in the visual system. Finally, the multivariate correlation analyses show that MEG connectomes have strong behavioral significance, which however depends on the considered connectivity measure and temporal scale. This comprehensive, albeit preliminary investigation of MEG connectome test-retest identifiability offers a first characterization of MEG fingerprinting in relation to different methodological and electrophysiological factors and contributes to the understanding of fingerprinting cross-modal relationships. We hope that this first investigation will contribute to setting the grounds for MEG connectome identification.

Dynamic relationships between spontaneous and evoked electrophysiological activity

Spontaneous neural activity fluctuations have been shown to influence trial-by-trial variation in perceptual, cognitive, and behavioral outcomes. However, the complex electrophysiological mechanisms by which these fluctuations shape stimulus-evoked neural activity remain largely to be explored. Employing a large-scale magnetoencephalographic dataset and an electroencephalographic replication dataset, we investigate the relationship between spontaneous and evoked neural activity across a range of electrophysiological variables. We observe that for high-frequency activity, high pre-stimulus amplitudes lead to greater evoked desynchronization, while for low frequencies, high pre-stimulus amplitudes induce larger degrees of event-related synchronization. We further decompose electrophysiological power into oscillatory and scale-free components, demonstrating different patterns of spontaneous-evoked correlation for each component. Finally, we find correlations between spontaneous and evoked time-domain electrophysiological signals. Overall, we demonstrate that the dynamics of multiple electrophysiological variables exhibit distinct relationships between their spontaneous and evoked activity, a result which carries implications for experimental design and analysis in non-invasive electrophysiology.

Statistical power: Implications for planning MEG studies

Statistical power is key for robust, replicable science. Here, we systematically explored how numbers of trials and subjects affect statistical power in MEG sensor-level data. More specifically, we simulated "experiments" using the MEG resting-state dataset of the Human Connectome Project (HCP). We divided the data in two conditions, injected a dipolar source at a known anatomical location in the "signal condition", but not in the "noise condition", and detected significant differences at sensor level with classical paired t-tests across subjects, using amplitude, squared amplitude, and global field power (GFP) measures. Group-level detectability of these simulated effects varied drastically with anatomical origin. We thus examined in detail which spatial properties of the sources affected detectability, looking specifically at the distance from closest sensor and orientation of the source, and at the variability of these parameters across subjects. In line with previous single-subject studies, we found that the most detectable effects originate from source locations that are closest to the sensors and oriented tangentially with respect to the head surface. In addition, cross-subject variability in orientation also affected group-level detectability, boosting detection in regions where this variability was small and hindering detection in regions where it was large. Incidentally, we observed a considerable covariation of source position, orientation, and their cross-subject variability in individual brain anatomical space, making it difficult to assess the impact of each of these variables independently of one another. We thus also performed simulations where we controlled spatial properties independently of individual anatomy. These additional simulations confirmed the strong impact of distance and orientation and further showed that orientation variability across subjects affects detectability, whereas position variability does not. Importantly, our study indicates that strict unequivocal recommendations as to the ideal number of trials and subjects for any experiment cannot be realistically provided for neurophysiological studies and should be adapted according to the brain regions under study.

Temporal modes of hub synchronization at rest

The brain is a dynamic system that generates a broad repertoire of perceptual, motor, and cognitive states by the integration and segregation of different functional domains represented in large-scale brain networks. However, the fundamental mechanisms underlying brain network integration remain elusive. Here, for the first time to our knowledge, we found that in the resting state the brain visits few synchronization modes defined as clusters of temporally aligned functional hubs. These modes alternate over time and their probability of switching leads to specific temporal loops among them. Notably, although each mode involves a small set of nodes, the brain integration seems highly vulnerable to a simulated attack on this temporal synchronization mechanism. In line with the hypothesis that the resting state represents a prior sculpted by the task activity, the observed synchronization modes might be interpreted as a temporal brain template needed to respond to task/environmental demands .

Temporal hierarchy of intrinsic neural timescales converges with spatial core-periphery organization

The human cortex exhibits intrinsic neural timescales that shape a temporal hierarchy. Whether this temporal hierarchy follows the spatial hierarchy of its topography, namely the core-periphery organization, remains an open issue. Using magnetoencephalography data, we investigate intrinsic neural timescales during rest and task states; we measure the autocorrelation window in short (ACW-50) and, introducing a novel variant, long (ACW-0) windows. We demonstrate longer ACW-50 and ACW-0 in networks located at the core compared to those at the periphery with rest and task states showing a high ACW correlation. Calculating rest-task differences, i.e., subtracting the shared core-periphery organization, reveals task-specific ACW changes in distinct networks. Finally, employing kernel density estimation, machine learning, and simulation, we demonstrate that ACW-0 exhibits better prediction in classifying a region's time window as core or periphery. Overall, our findings provide fundamental insight into how the human cortex's temporal hierarchy converges with its spatial core-periphery hierarchy.

Decoding working memory task condition using magnetoencephalography source level long-range phase coupling patterns

OBJECTIVE

The objective of the study is to identify phase coupling patterns that are shared across subjects via a machine learning approach that utilises source space magnetoencephalography (MEG) phase coupling data from a working memory (WM) task. Indeed, phase coupling of neural oscillations is putatively a key factor for communication between distant brain areas and is therefore crucial in performing cognitive tasks, including WM. Previous studies investigating phase coupling during cognitive tasks have often focused on a few a priori selected brain areas or a specific frequency band, and the need for data-driven approaches has been recognised. Machine learning techniques have emerged as valuable tools for the analysis of neuroimaging data since they catch fine-grained differences in the multivariate signal distribution. Here, we expect that these techniques applied to MEG phase couplings can reveal WM-related processes that are shared across individuals.

APPROACH

We analysed WM data collected as part of the Human Connectome Project. The MEG data were collected while subjects (n = 83) performed N-back WM tasks in two different conditions, namely 2-back (WM condition) and 0-back (control condition). We estimated phase coupling patterns (multivariate phase slope index) for both conditions and for theta, alpha, beta, and gamma bands. The obtained phase coupling data were then used to train a linear support vector machine in order to classify which task condition the subject was performing with an across-subject cross-validation approach. The classification was performed separately based on the data from individual frequency bands and with all bands combined (multiband). Finally, we evaluated the relative importance of the different features (phase couplings) for classification by the means of feature selection probability.

MAIN RESULTS

The WM condition and control condition were successfully classified based on the phase coupling patterns in the theta (62% accuracy) and alpha bands (60% accuracy) separately. Importantly, the multiband classification showed that phase coupling patterns not only in the theta and alpha but also in the gamma bands are related to WM processing, as testified by improvement in classification performance (71%).

SIGNIFICANCE

Our study successfully decoded WM tasks using MEG source space functional connectivity. Our approach, combining across-subject classification and a multidimensional metric recently developed by our group, is able to detect patterns of connectivity that are shared across individuals. In other words, the results are generalisable to new individuals and allow meaningful interpretation of task-relevant phase coupling patterns.

A Whole-Cortex Probabilistic Diffusion Tractography Connectome

The WU-Minn Human Connectome Project (HCP) is a publicly-available dataset containing state-of-the-art structural magnetic resonance imaging (MRI), functional MRI (fMRI), and diffusion MRI (dMRI) for over a thousand healthy subjects. While the planned scope of the HCP included an anatomic connectome, resting-state fMRI (rs-fMRI) forms the bulk of the HCP's current connectomic output. We address this by presenting a full-cortex connectome derived from probabilistic diffusion tractography and organized into the HCP-MMP1.0 atlas. Probabilistic methods and large sample sizes are preferable for whole-connectome mapping as they increase the fidelity of traced low-probability connections. We find that overall, connection strengths are lognormally distributed and decay exponentially with tract length, that connectivity reasonably matches macaque histologic tracing in homologous areas, that contralateral homologs and left-lateralized language areas are hyperconnected, and that hierarchical similarity influences connectivity. We compare the dMRI connectome to existing rs-fMRI and cortico-cortico-evoked potential connectivity matrices and find that it is more similar to the latter. This work helps fulfill the promise of the HCP and will make possible comparisons between the underlying structural connectome and functional connectomes of various modalities, brain states, and clinical conditions.

HD-EEG for tracking sub-second brain dynamics during cognitive tasks

This work provides the community with high-density Electroencephalography (HD-EEG, 256 channels) datasets collected during task-free and task-related paradigms. It includes forty-three healthy participants performing visual naming and spelling tasks, visual and auditory naming tasks and a visual working memory task in addition to resting state. The HD-EEG data are furnished in the Brain Imaging Data Structure (BIDS) format. These datasets can be used to (i) track brain networks dynamics and their rapid reconfigurations at sub-second time scale in different conditions, (naming/spelling/rest) and modalities, (auditory/visual) and compare them to each other, (ii) validate several parameters involved in the methods used to estimate cortical brain networks through scalp EEG, such as the open question of optimal number of channels and number of regions of interest and (iii) allow the reproducibility of results obtained so far using HD-EEG. We hope that delivering these datasets will lead to the development of new methods that can be used to estimate brain cortical networks and to better understand the general functioning of the brain during rest and task. Data are freely available from https://openneuro.org.

Measurement(s)	brain measurement • cognitive behavior trait
Technology Type(s)	electroencephalography (EEG)
Factor Type(s)	task • age • sex
Sample Characteristic - Organism	Homo sapiens

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13560311

Multi-band MEG signatures of BOLD connectivity reorganization during visuospatial attention

The functional architecture of the resting brain, as measured with the blood oxygenation level-dependent functional connectivity (BOLD-FC), is slightly modified during task performance. In previous work, we reported behaviorally relevant BOLD-FC modulations between visual and dorsal attention regions when subjects performed a visuospatial attention task as compared to central fixation (Spadone et al., 2015). Here we use magnetoencephalography (MEG) in the same group of subjects to identify the electrophysiological correlates of the BOLD-FC modulation found in our previous work. While BOLD-FC topography, separately at rest and during visual attention, corresponded to neuromagnetic Band-Limited Power (BLP) correlation in the alpha and beta bands (8-30 Hz), BOLD-FC modulations evoked by performing the visual attention task (Spadone et al. 2015) did not match any specific oscillatory band BLP modulation. Conversely, following the application of an orthogonal spatial decomposition that identifies common inter-subject co-variations, we found that attention-rest BOLD-FC modulations were recapitulated by multi-spectral BLP-FC components. Notably, individual variability of alpha connectivity between Frontal Eye Fields and visual occipital regions, jointly with decreased interaction in the Visual network, correlated with visual discrimination accuracy. In summary, task-rest BOLD connectivity modulations match multi-spectral MEG BLP connectivity.

A Whole-Cortex Probabilistic Diffusion Tractography Connectome

The WU-Minn Human Connectome Project (HCP) is a publicly-available dataset containing state-of-the-art structural magnetic resonance imaging (MRI), functional MRI (fMRI), and diffusion MRI (dMRI) for over a thousand healthy subjects. While the planned scope of the HCP included an anatomic connectome, resting-state fMRI (rs-fMRI) forms the bulk of the HCP's current connectomic output. We address this by presenting a full-cortex connectome derived from probabilistic diffusion tractography and organized into the HCP-MMP1.0 atlas. Probabilistic methods and large sample sizes are preferable for whole-connectome mapping as they increase the fidelity of traced low-probability connections. We find that overall, connection strengths are lognormally distributed and decay exponentially with tract length, that connectivity reasonably matches macaque histologic tracing in homologous areas, that contralateral homologs and left-lateralized language areas are hyperconnected, and that hierarchical similarity influences connectivity. We compare the dMRI connectome to existing rs-fMRI and cortico-cortico-evoked potential connectivity matrices and find that it is more similar to the latter. This work helps fulfill the promise of the HCP and will make possible comparisons between the underlying structural connectome and functional connectomes of various modalities, brain states, and clinical conditions.

A Connectome-Based, Corticothalamic Model of State- and Stimulation-Dependent Modulation of Rhythmic Neural Activity and Connectivity

Rhythmic activity in the brain fluctuates with behaviour and cognitive state, through a combination of coexisting and interacting frequencies. At large spatial scales such as those studied in human M/EEG, measured oscillatory dynamics are believed to arise primarily from a combination of cortical (intracolumnar) and corticothalamic rhythmogenic mechanisms. Whilst considerable progress has been made in characterizing these two types of neural circuit separately, relatively little work has been done that attempts to unify them into a single consistent picture. This is the aim of the present paper. We present and examine a whole-brain, connectome-based neural mass model with detailed long-range cortico-cortical connectivity and strong, recurrent corticothalamic circuitry. This system reproduces a variety of known features of human M/EEG recordings, including spectral peaks at canonical frequencies, and functional connectivity structure that is shaped by the underlying anatomical connectivity. Importantly, our model is able to capture state- (e.g., idling/active) dependent fluctuations in oscillatory activity and the coexistence of multiple oscillatory phenomena, as well as frequency-specific modulation of functional connectivity. We find that increasing the level of sensory drive to the thalamus triggers a suppression of the dominant low frequency rhythms generated by corticothalamic loops, and subsequent disinhibition of higher frequency endogenous rhythmic behaviour of intracolumnar microcircuits. These combine to yield simultaneous decreases in lower frequency and increases in higher frequency components of the M/EEG power spectrum during states of high sensory or cognitive drive. Building on this, we also explored the effect of pulsatile brain stimulation on ongoing oscillatory activity, and evaluated the impact of coexistent frequencies and state-dependent fluctuations on the response of cortical networks. Our results provide new insight into the role played by cortical and corticothalamic circuits in shaping intrinsic brain rhythms, and suggest new directions for brain stimulation therapies aimed at state-and frequency-specific control of oscillatory brain activity.

Dynamics in cortical activity revealed by resting-state MEG rhythms

The brain is a biophysical system subject to information flows that may be thought of as a many-body architecture with a spatiotemporal dynamics described by its neuronal structures. The oscillatory nature of brain activity allows these structures (nodes) to be described as a set of coupled oscillators forming a network where the node dynamics and that of the network topology can be studied. Quantifying its dynamics at various scales is an issue that claims to be explored for several brain activities, e.g., activity at rest. The resting-state (RS) associates the underlying brain dynamics of healthy subjects that are not actively compromised with sensory or cognitive processes. Studying its dynamics is highly non-trivial but opens the door to understand the general principles of brain functioning, as well as to contrast a passive null condition vs the dynamics of pathologies or non-resting activities. Here, we hypothesize about how the spatiotemporal dynamics of cortical fluctuations could be for healthy subjects at RS. To do that, we retrieve the alphabet that reconstructs the dynamics (entropy-complexity) of magnetoencephalography (MEG) signals. We assemble the cortical connectivity to elicit the dynamics in the network topology. We depict an order relation between entropy and complexity for frequency bands that is ubiquitous for different temporal scales. We unveiled that the posterior cortex conglomerates nodes with both stronger dynamics and high clustering for α band. The existence of an order relation between dynamic properties suggests an emergent phenomenon characteristic of each band. Interestingly, we find the posterior cortex as a domain of dual character that plays a cardinal role in both the dynamics and structure regarding the activity at rest. To the best of our knowledge, this is the first study with MEG involving information theory and network science to better understand the dynamics and structure of brain activity at rest for different bands and scales.

Resolving the Connectome, Spectrally-Specific Functional Connectivity Networks and Their Distinct Contributions to Behavior

The resting human brain exhibits spontaneous patterns of activity that reflect features of the underlying neural substrate. Examination of interareal coupling of resting-state oscillatory activity has revealed that the brain’s resting activity is composed of functional networks, whose topographies differ depending on oscillatory frequency, suggesting a role for carrier frequency as a means of creating multiplexed, or functionally segregated, communication channels between brain areas. Using canonical correlation analysis (CCA), we examined spectrally resolved resting-state connectivity patterns derived from magnetoencephalography (MEG) recordings to determine the relationship between connectivity intrinsic to different frequency channels and a battery of over a hundred behavioral and demographic indicators, in a group of 89 young healthy participants. We demonstrate that each of the classical frequency bands in the range 1–40 Hz (δ, θ, α, β, and γ) delineates a subnetwork that is behaviorally relevant, spatially distinct, and whose expression is either negatively or positively predictive of individual traits, with the strongest link in the α-band being negative and networks oscillating at different frequencies, such as θ, β, and γ carrying positive function.

Non-invasive inference of information flow using diffusion MRI, functional MRI, and MEG

OBJECTIVE

To infer information flow in the white matter of the brain and recover cortical activity using functional MRI, diffusion MRI, and MEG without a manual selection of the white matter connections of interest.

APPROACH

A Bayesian network which encodes the priors knowledge of possible brain states is built from imaging data. Diffusion MRI is used to enumerate all possible connections between cortical regions. Functional MRI is used to prune connections without manual intervention and increase the likelihood of specific regions being active. MEG data is used as evidence into this network to obtain a posterior distribution on cortical regions and connections.

MAIN RESULTS

We show that our proposed method is able to identify connections associated with the a sensory-motor task. This allows us to build the Bayesian network with no manual selection of connections of interest. Using sensory-motor MEG evoked response as evidence into this network, our method identified areas known to be involved in a visuomotor task. In addition, information flow along white matter fiber bundles connecting those regions was also recovered.

SIGNIFICANCE

Current methods to estimate white matter information flow are extremely invasive, therefore limiting our understanding of the interaction between cortical regions. The proposed method makes use of functional MRI, diffusion MRI, and M/EEG to infer communication between cortical regions, therefore opening the door to the non-invasive exploration of information flow in the white matter.

Reduced auditory steady state responses in autism spectrum disorder

BACKGROUND

Auditory steady state responses (ASSRs) are elicited by clicktrains or amplitude-modulated tones, which entrain auditory cortex at their specific modulation rate. Previous research has reported reductions in ASSRs at 40 Hz for autism spectrum disorder (ASD) participants and first-degree relatives of people diagnosed with ASD (Mol Autism. 2011;2:11, Biol Psychiatry. 2007;62:192-197).

METHODS

Using a 1.5 s-long auditory clicktrain stimulus, designed to elicit an ASSR at 40 Hz, this study attempted to replicate and extend these findings. Magnetencephalography (MEG) data were collected from 18 adolescent ASD participants and 18 typically developing controls.

RESULTS

The ASSR localised to bilateral primary auditory regions. Regions of interest were thus defined in left and right primary auditory cortex (A1). While the transient gamma-band response (tGBR) from 0-0.1 s following presentation of the clicktrain stimulus was not different between groups, for either left or right A1, the ASD group had reduced oscillatory power at 40 Hz from 0.5 to 1.5 s post-stimulus onset, for both left and right A1. Additionally, the ASD group had reduced inter-trial coherence (phase consistency over trials) at 40 Hz from 0.64-0.82 s for right A1 and 1.04-1.22 s for left A1.

LIMITATIONS

In this study, we did not conduct a clinical autism assessment (e.g. the ADOS), and therefore, it remains unclear whether ASSR power and/or ITC are associated with the clinical symptoms of ASD.

CONCLUSION

Overall, our results support a specific reduction in ASSR oscillatory power and inter-trial coherence in ASD, rather than a generalised deficit in gamma-band responses. We argue that this could reflect a developmentally relevant reduction in non-linear neural processing.

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.

Combining magnetoencephalography with magnetic resonance imaging enhances learning of surrogate-biomarkers

Electrophysiological methods, that is M/EEG, provide unique views into brain health. Yet, when building predictive models from brain data, it is often unclear how electrophysiology should be combined with other neuroimaging methods. Information can be redundant, useful common representations of multimodal data may not be obvious and multimodal data collection can be medically contraindicated, which reduces applicability. Here, we propose a multimodal model to robustly combine MEG, MRI and fMRI for prediction. We focus on age prediction as a surrogate biomarker in 674 subjects from the Cam-CAN dataset. Strikingly, MEG, fMRI and MRI showed additive effects supporting distinct brain-behavior associations. Moreover, the contribution of MEG was best explained by cortical power spectra between 8 and 30 Hz. Finally, we demonstrate that the model preserves benefits of stacking when some data is missing. The proposed framework, hence, enables multimodal learning for a wide range of biomarkers from diverse types of brain signals.

Predictive regression modeling with MEG/EEG: from source power to signals and cognitive states

Predicting biomedical outcomes from Magnetoencephalography and Electroencephalography (M/EEG) is central to applications like decoding, brain-computer-interfaces (BCI) or biomarker development and is facilitated by supervised machine learning. Yet, most of the literature is concerned with classification of outcomes defined at the event-level. Here, we focus on predicting continuous outcomes from M/EEG signal defined at the subject-level, and analyze about 600 MEG recordings from Cam-CAN dataset and about 1000 EEG recordings from TUH dataset. Considering different generative mechanisms for M/EEG signals and the biomedical outcome, we propose statistically-consistent predictive models that avoid source-reconstruction based on the covariance as representation. Our mathematical analysis and ground-truth simulations demonstrated that consistent function approximation can be obtained with supervised spatial filtering or by embedding with Riemannian geometry. Additional simulations revealed that Riemannian methods were more robust to model violations, in particular geometric distortions induced by individual anatomy. To estimate the relative contribution of brain dynamics and anatomy to prediction performance, we propose a novel model inspection procedure based on biophysical forward modeling. Applied to prediction of outcomes at the subject-level, the analysis revealed that the Riemannian model better exploited anatomical information while sensitivity to brain dynamics was similar across methods. We then probed the robustness of the models across different data cleaning options. Environmental denoising was globally important but Riemannian models were strikingly robust and continued performing well even without preprocessing. Our results suggest each method has its niche: supervised spatial filtering is practical for event-level prediction while the Riemannian model may enable simple end-to-end learning.

Transient spectral events in resting state MEG predict individual task responses

Even in response to simple tasks such as hand movement, human brain activity shows remarkable inter-subject variability. Recently, it has been shown that individual spatial variability in fMRI task responses can be predicted from measurements collected at rest; suggesting that the spatial variability is a stable feature, inherent to the individual's brain. However, it is not clear if this is also true for individual variability in the spatio-spectral content of oscillatory brain activity. Here, we show using MEG (N ​= ​89) that we can predict the spatial and spectral content of an individual's task response using features estimated from the individual's resting MEG data. This works by learning when transient spectral 'bursts' or events in the resting state tend to reoccur in the task responses. We applied our method to motor, working memory and language comprehension tasks. All task conditions were predicted significantly above chance. Finally, we found a systematic relationship between genetic similarity (e.g. unrelated subjects vs. twins) and predictability. Our approach can predict individual differences in brain activity and suggests a link between transient spectral events in task and rest that can be captured at the level of individuals.

Magnetoencephalography in Cognitive Neuroscience: A Primer

Magnetoencephalography (MEG) is an invaluable tool to study the dynamics and connectivity of large-scale brain activity and their interactions with the body and the environment in functional and dysfunctional body and brain states. This primer introduces the basic concepts of MEG, discusses its strengths and limitations in comparison to other brain imaging techniques, showcases interesting applications, and projects exciting current trends into the near future, in a way that might more fully exploit the unique capabilities of MEG.

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.

The Impact of the Geometric Correction Scheme on MEG Functional Topology at Rest

Spontaneous activity is correlated across brain regions in large scale networks (RSN) closely resembling those recruited during several behavioral tasks and characterized by functional specialization and dynamic integration. Specifically, MEG studies revealed a set of central regions (dynamic core) possibly facilitating communication among differently specialized brain systems. However, source projected MEG signals, due to the fundamentally ill-posed inverse problem, are affected by spatial leakage, leading to the estimation of spurious, blurred connections that may affect the topological properties of brain networks and their integration. To reduce leakage effects, several correction schemes have been proposed including the Geometric Correction Scheme (GCS) whose theory, simulations and empirical results on topography of a few RSNs were already presented. However, its impact on the estimation of fundamental graph measures used to describe the architecture of interactions among brain regions has not been investigated yet. Here, we estimated dense, MEG band-limited power connectomes in theta, alpha, beta, and gamma bands from 13 healthy subjects (all young adults). We compared the connectivity and topology of MEG uncorrected and GCS-corrected connectomes. The use of GCS considerably reorganized the topology of connectivity, reducing the local, within-hemisphere interactions mainly in the beta and gamma bands and increasing across-hemisphere interactions mainly in the alpha and beta bands. Moreover, the number of hubs decreased in the alpha and beta bands, but the centrality of some fundamental regions such as the Posterior Cingulate Cortex (PCC), Supplementary Motor Area (SMA) and Middle Prefrontal Cortex (MPFC) remained strong in all bands, associated to an increase of the Global Efficiency and a decrease of Modularity. As a comparison, we applied orthogonalization on connectomes and ran the same topological analyses. The correlation values were considerably reduced, and orthogonalization mainly decreased local within-hemisphere interactions in all bands, similarly to GCS. Notably, the centrality of the PCC, SMA and MPFC was preserved in all bands, as for GCS, together with other hubs in the posterior parietal regions. Overall, leakage correction removes spurious local connections, but confirms the role of dynamic hub regions, specifically the anterior and posterior cingulate, in integrating information in the brain at rest.

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.

Heartbeat Induces a Cortical Theta-Synchronized Network in the Resting State

In the resting state, heartbeats evoke cortical responses called heartbeat-evoked responses (HERs), which reflect cortical cardiac interoceptive processing. While previous studies have reported that the heartbeat evokes cortical responses at a regional level, whether the heartbeat induces synchronization between regions to form a network structure remains unknown. Using resting-state MEG data from 85 human subjects of both genders, we first showed that heartbeat increases the phase synchronization between cortical regions in the theta frequency but not in other frequency bands. This increase in synchronization between cortical regions formed a network structure called the heartbeat-induced network (HIN), which did not reflect artificial phase synchronization. In the HIN, the left inferior temporal gyrus and parahippocampal gyrus played a central role as hubs. Furthermore, the HIN was modularized, containing five subnetworks called modules. In particular, module 1 played a central role in between-module interactions in the HIN. Furthermore, synchronization within module 1 had a positive association with the mood of an individual. In this study, we show the existence of the HIN and its network properties, advancing the current understanding of cortical heartbeat processing and its relationship with mood, which was previously confined to region level.

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.

Brainstorm Pipeline Analysis of Resting-State Data From the Open MEG Archive

We present a simple, reproducible analysis pipeline applied to resting-state magnetoencephalography (MEG) data from the Open MEG Archive (OMEGA). The data workflow was implemented with Brainstorm, which like OMEGA is free and openly accessible. The proposed pipeline produces group maps of ongoing brain activity decomposed in the typical frequency bands of electrophysiology. The procedure is presented as a technical proof of concept for streamlining a broader range and more sophisticated studies of resting-state electrophysiological data. It also features the recently introduced extension of the brain imaging data structure (BIDS) to MEG data, highlighting the scalability and generalizability of Brainstorm analytical pipelines to other, and potentially larger data volumes.

A 204-subject multimodal neuroimaging dataset to study language processing

This dataset, colloquially known as the Mother Of Unification Studies (MOUS) dataset, contains multimodal neuroimaging data that has been acquired from 204 healthy human subjects. The neuroimaging protocol consisted of magnetic resonance imaging (MRI) to derive information at high spatial resolution about brain anatomy and structural connections, and functional data during task, and at rest. In addition, magnetoencephalography (MEG) was used to obtain high temporal resolution electrophysiological measurements during task, and at rest. All subjects performed a language task, during which they processed linguistic utterances that either consisted of normal or scrambled sentences. Half of the subjects were reading the stimuli, the other half listened to the stimuli. The resting state measurements consisted of 5 minutes eyes-open for the MEG and 7 minutes eyes-closed for fMRI. The neuroimaging data, as well as the information about the experimental events are shared according to the Brain Imaging Data Structure (BIDS) format. This unprecedented neuroimaging language data collection allows for the investigation of various aspects of the neurobiological correlates of language.

Hierarchical Heterogeneity across Human Cortex Shapes Large-Scale Neural Dynamics

The large-scale organization of dynamical neural activity across cortex emerges through long-range interactions among local circuits. We hypothesized that large-scale dynamics are also shaped by heterogeneity of intrinsic local properties across cortical areas. One key axis along which microcircuit properties are specialized relates to hierarchical levels of cortical organization. We developed a large-scale dynamical circuit model of human cortex that incorporates heterogeneity of local synaptic strengths, following a hierarchical axis inferred from magnetic resonance imaging (MRI)-derived T1- to T2-weighted (T1w/T2w) mapping and fit the model using multimodal neuroimaging data. We found that incorporating hierarchical heterogeneity substantially improves the model fit to functional MRI (fMRI)-measured resting-state functional connectivity and captures sensory-association organization of multiple fMRI features. The model predicts hierarchically organized higher-frequency spectral power, which we tested with resting-state magnetoencephalography. These findings suggest circuit-level mechanisms linking spatiotemporal levels of analysis and highlight the importance of local properties and their hierarchical specialization on the large-scale organization of human cortical dynamics.

Sad faces increase the heartbeat-associated interoceptive information flow within the salience network: a MEG study

The somatic marker hypothesis proposes that the cortical representation of visceral signals is a crucial component of emotional processing. No previous study has investigated the information flow among brain regions that process visceral information during emotional perception. In this magnetoencephalography study of 32 healthy subjects of either sex, heartbeat-evoked responses (HERs), which reflect the cortical processing of heartbeats, were modulated by the perception of a sad face. The modulation effect was localized to the prefrontal cortices, the globus pallidus, and an interoceptive network including the right anterior insula (RAI) and dorsal anterior cingulate cortex (RdACC). Importantly, our Granger causality analysis provides the first evidence for the increased flow of heartbeat information from the RAI to the RdACC during sad face perception. Moreover, using a surrogate R-peak analysis, we have shown that this HER modulation effect was time-locked to heartbeats. These findings advance the understanding of brain-body interactions during emotional processing.

Deep Convolutional Neural Networks for Feature-Less Automatic Classification of Independent Components in Multi-Channel Electrophysiological Brain Recordings

OBJECTIVE

Interpretation of the electroencephalographic (EEG) and magnetoencephalographic (MEG) signals requires off-line artifacts removal. Since artifacts share frequencies with brain activity, filtering is insufficient. Blind source separation, mainly through independent component analysis (ICA), is the gold-standard procedure for the identification of artifacts in multi-dimensional recordings. However, a classification of brain and artifactual independent components (ICs) is still required. Since ICs exhibit recognizable patterns, classification is usually performed by experts' visual inspection. This procedure is time consuming and prone to errors. Automatic ICs classification has been explored, often through complex ICs features extraction prior to classification. Relying on deep-learning ability of self-extracting the features of interest, we investigated the capabilities of convolutional neural networks (CNNs) for off-line, automatic artifact identification through ICs without feature selection.

METHODS

A CNN was applied to spectrum and topography of a large dataset of few thousand samples of ICs obtained from multi-channel EEG and MEG recordings acquired during heterogeneous experimental settings and on different subjects. CNN performances, when applied to EEG, MEG, and combined EEG and MEG ICs, were explored and compared with state-of-the-art feature-based automatic classification.

RESULTS

Beyond state-of-the-art automatic classification accuracies were demonstrated through cross validation (92.4% EEG, 95.4% MEG, 95.6% EEG+MEG).

CONCLUSION

High CNN classification performances were achieved through heuristical selection of machinery hyperparameters and through the CNN self-selection of the features of interest.

SIGNIFICANCE

Considering the large data availability of multi-channel EEG and MEG recordings, CNNs may be suited for classification of ICs of multi-channel brain electrophysiological recordings.

How do spatially distinct frequency specific MEG networks emerge from one underlying structural connectome? The role of the structural eigenmodes

Functional networks obtained from magnetoencephalography (MEG) from different frequency bands show distinct spatial patterns. It remains to be elucidated how distinct spatial patterns in MEG networks emerge given a single underlying structural network. Recent work has suggested that the eigenmodes of the structural network might serve as a basis set for functional network patterns in the case of functional MRI. Here, we take this notion further in the context of frequency band specific MEG networks. We show that a selected set of eigenmodes of the structural network can predict different frequency band specific networks in the resting state, ranging from delta (1-4 Hz) to the high gamma band (40-70 Hz). These predictions outperform predictions based from surrogate data, suggesting a genuine relationship between eigenmodes of the structural network and frequency specific MEG networks. We then show that the relevant set of eigenmodes can be excited in a network of neural mass models using linear stability analysis only by including delays. Excitation of an eigenmode in this context refers to a dynamic instability of a network steady state to a spatial pattern with a corresponding coherent temporal oscillation. Simulations verify the results from linear stability analysis and suggest that theta, alpha and beta band networks emerge very near to the bifurcation. The delta and gamma bands in the resting state emerges further away from the bifurcation. These results show for the first time how delayed interactions can excite the relevant set of eigenmodes that give rise to frequency specific functional connectivity patterns.

Quantifying the performance of MEG source reconstruction using resting state data

In magnetoencephalography (MEG) research there are a variety of inversion methods to transform sensor data into estimates of brain activity. Each new inversion scheme is generally justified against a specific simulated or task scenario. The choice of this scenario will however have a large impact on how well the scheme performs. We describe a method with minimal selection bias to quantify algorithm performance using human resting state data. These recordings provide a generic, heterogeneous, and plentiful functional substrate against which to test different MEG recording and reconstruction approaches. We used a Hidden Markov model to spatio-temporally partition data into self-similar dynamic states. To test the anatomical precision that could be achieved, we then inverted these data onto libraries of systematically distorted subject-specific cortical meshes and compared the quality of the fit using cross validation and a Free energy metric. This revealed which inversion scheme was able to identify the least distorted (most accurate) anatomical models, and allowed us to quantify an upper bound on the mean anatomical distortion accordingly. We used two resting state datasets, one recorded with head-casts and one without. In the head-cast data, the Empirical Bayesian Beamformer (EBB) algorithm showed the best mean anatomical discrimination (3.7 mm) compared with Minimum Norm/LORETA (6.0 mm) and Multiple Sparse Priors (9.4 mm). This pattern was replicated in the second (conventional dataset) although with a marginally poorer (non-significant) prediction of the missing (cross-validated) data. Our findings suggest that the abundant resting state data now commonly available could be used to refine and validate MEG source reconstruction methods and/or recording paradigms.

Discovering dynamic brain networks from big data in rest and task

Brain activity is a dynamic combination of the responses to sensory inputs and its own spontaneous processing. Consequently, such brain activity is continuously changing whether or not one is focusing on an externally imposed task. Previously, we have introduced an analysis method that allows us, using Hidden Markov Models (HMM), to model task or rest brain activity as a dynamic sequence of distinct brain networks, overcoming many of the limitations posed by sliding window approaches. Here, we present an advance that enables the HMM to handle very large amounts of data, making possible the inference of very reproducible and interpretable dynamic brain networks in a range of different datasets, including task, rest, MEG and fMRI, with potentially thousands of subjects. We anticipate that the generation of large and publicly available datasets from initiatives such as the Human Connectome Project and UK Biobank, in combination with computational methods that can work at this scale, will bring a breakthrough in our understanding of brain function in both health and disease.

Frequency-based brain networks: From a multiplex framework to a full multilayer description

We explore how to study dynamical interactions between brain regions by using functional multilayer networks whose layers represent different frequency bands at which a brain operates. Specifically, we investigate the consequences of considering the brain as (i) a multilayer network, in which all brain regions can interact with each other at different frequency bands; and as (ii) a multiplex network, in which interactions between different frequency bands are allowed only within each brain region and not between them. We study the second-smallest eigenvalue λ 2 of the combinatorial supra-Laplacian matrix of both the multiplex and multilayer networks, as λ 2 has been used previously as an indicator of network synchronizability and as a biomarker for several brain diseases. We show that the heterogeneity of interlayer edge weights and, especially, the fraction of missing edges crucially modify the value of λ 2, and we illustrate our results with both synthetic network models and real data obtained from resting-state magnetoencephalography. Our work highlights the differences between using a multiplex approach and a full multilayer approach when studying frequency-based multilayer brain networks.

MEG-BIDS, the brain imaging data structure extended to magnetoencephalography

We present a significant extension of the Brain Imaging Data Structure (BIDS) to support the specific aspects of magnetoencephalography (MEG) data. MEG measures brain activity with millisecond temporal resolution and unique source imaging capabilities. So far, BIDS was a solution to organise magnetic resonance imaging (MRI) data. The nature and acquisition parameters of MRI and MEG data are strongly dissimilar. Although there is no standard data format for MEG, we propose MEG-BIDS as a principled solution to store, organise, process and share the multidimensional data volumes produced by the modality. The standard also includes well-defined metadata, to facilitate future data harmonisation and sharing efforts. This responds to unmet needs from the multimodal neuroimaging community and paves the way to further integration of other techniques in electrophysiology. MEG-BIDS builds on MRI-BIDS, extending BIDS to a multimodal data structure. We feature several data-analytics software that have adopted MEG-BIDS, and a diverse sample of open MEG-BIDS data resources available to everyone.

Predicting seizure outcome of vagus nerve stimulation using MEG-based network topology

Vagus nerve stimulation (VNS) is a low-risk surgical option for patients with drug resistant epilepsy, although it is impossible to predict which patients may respond to VNS treatment. Resting-state magnetoencephalography (rs-MEG) connectivity analysis has been increasingly utilized to investigate the impact of epilepsy on brain networks and identify alteration of these networks after different treatments; however, there is no study to date utilizing this modality to predict the efficacy of VNS treatment. We investigated whether the rs-MEG network topology before VNS implantation can be used to predict efficacy of VNS treatment. Twenty-three patients with epilepsy who had MEG before VNS implantation were included in this study. We also included 89 healthy control subjects from the Human Connectome Project. Using the phase-locking value in the theta, alpha, and beta frequency bands as a measure of rs-MEG functional connectivity, we calculated three global graph measures: modularity, transitivity, and characteristic path length (CPL). Our results revealed that the rs-MEG graph measures were significantly heritable and had an overall good test-retest reliability, and thus these measures may be used as potential biomarkers of the network topology. We found that the modularity and transitivity in VNS responders were significantly larger and smaller, respectively, than those observed in VNS non-responders. We also observed that the modularity and transitivity in three frequency bands and CPL in delta and beta bands were significantly different in controls than those found in responders or non-responders, although the values of the graph measures in controls were closer to those of responders than non-responders. We used the modularity and transitivity as input features of a naïve Bayes classifier, and achieved an accuracy of 87% in classification of non-responders, responders, and controls. The results of this study revealed that MEG-based graph measures are reliable biomarkers, and that these measures may be used to predict seizure outcome of VNS treatment.

Multi-subject hierarchical inverse covariance modelling improves estimation of functional brain networks

A Bayesian model for sparse, hierarchical, inver-covariance estimation is presented, and applied to multi-subject functional connectivity estimation in the human brain. It enables simultaneous inference of the strength of connectivity between brain regions at both subject and population level, and is applicable to fMRI, MEG and EEG data. Two versions of the model can encourage sparse connectivity, either using continuous priors to suppress irrelevant connections, or using an explicit description of the network structure to estimate the connection probability between each pair of regions. A large evaluation of this model, and thirteen methods that represent the state of the art of inverse covariance modelling, is conducted using both simulated and resting-state functional imaging datasets. Our novel Bayesian approach has similar performance to the best extant alternative, Ng et al.'s Sparse Group Gaussian Graphical Model algorithm, which also is based on a hierarchical structure. Using data from the Human Connectome Project, we show that these hierarchical models are able to reduce the measurement error in MEG beta-band functional networks by 10%, producing concomitant increases in estimates of the genetic influence on functional connectivity.

Efficient communication dynamics on macro-connectome, and the propagation speed

Global communication dynamics in the brain can be captured using fMRI, MEG, or electrocorticography (ECoG), and the global slow dynamics often represent anatomical constraints. Complementary single-/multi-unit recordings have described local fast temporal dynamics. However, global fast temporal dynamics remain incompletely understood with considering of anatomical constraints. Therefore, we compared temporal aspects of cross-area propagations of single-unit recordings and ECoG, and investigated their anatomical bases. First, we demonstrated how both evoked and spontaneous ECoGs can accurately predict latencies of single-unit recordings. Next, we estimated the propagation velocity (1.0–1.5 m/s) from brain-wide data and found that it was fairly stable among different conscious levels. We also found that the shortest paths in anatomical topology strongly predicted the latencies. Finally, we demonstrated that Communicability, a novel graph-theoretic measure, is able to quantify that more than 90% of paths should use shortest paths and the remaining are non-shortest walks. These results revealed that macro-connectome is efficiently wired for detailed communication dynamics in the brain.