Predicting haemodynamic networks using electrophysiology: The role of non-linear and cross-frequency interactions

Spectral-switching analysis reveals real-time neuronal network representations of concurrent spontaneous naturalistic behaviors in human brain

Despite abundant evidence of functional networks in the human brain, their neuronal underpinnings, and relationships to real-time behavior have been challenging to resolve. Analyzing brain-wide intracranial-EEG recordings with video monitoring, acquired in awake subjects during clinical epilepsy evaluation, we discovered the tendency of each brain region to switch back and forth between 2 distinct power spectral densities (PSDs 2-55Hz). We further recognized that this 'spectral switching' occurs synchronously between distant sites, even between regions with differing baseline PSDs, revealing long-range functional networks that would be obscured in analysis of individual frequency bands. Moreover, the real-time PSD-switching dynamics of specific networks exhibited striking alignment with activities such as conversation and hand movements, revealing a multi-threaded functional network representation of concurrent naturalistic behaviors. Network structures and their relationships to behaviors were stable across days, but were altered during N3 sleep. Our results provide a new framework for understanding real-time, brain-wide neural-network dynamics.

Altered correlation of concurrently recorded EEG-fMRI connectomes in temporal lobe epilepsy

Whole-brain functional connectivity networks (connectomes) have been characterized at different scales in humans using EEG and fMRI. Multimodal epileptic networks have also been investigated, but the relationship between EEG and fMRI defined networks on a whole-brain scale is unclear. A unified multimodal connectome description, mapping healthy and pathological networks would close this knowledge gap. Here, we characterize the spatial correlation between the EEG and fMRI connectomes in right and left temporal lobe epilepsy (rTLE/lTLE). From two centers, we acquired resting-state concurrent EEG-fMRI of 35 healthy controls and 34 TLE patients. EEG-fMRI data was projected into the Desikan brain atlas, and functional connectomes from both modalities were correlated. EEG and fMRI connectomes were moderately correlated. This correlation was increased in rTLE when compared to controls for EEG-delta/theta/alpha/beta. Conversely, multimodal correlation in lTLE was decreased in respect to controls for EEG-beta. While the alteration was global in rTLE, in lTLE it was locally linked to the default mode network. The increased multimodal correlation in rTLE and decreased correlation in lTLE suggests a modality-specific lateralized differential reorganization in TLE, which needs to be considered when comparing results from different modalities. Each modality provides distinct information, highlighting the benefit of multimodal assessment in epilepsy.

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.

Data-driven MEG analysis to extract fMRI resting-state networks

The electrophysiological basis of resting-state networks (RSN) is still under debate. In particular, no principled mechanism has been determined that is capable of explaining all RSN equally well. While magnetoencephalography (MEG) and electroencephalography are the methods of choice to determine the electrophysiological basis of RSN, no standard analysis pipeline of RSN yet exists. In this article, we compare the two main existing data-driven analysis strategies for extracting RSNs from MEG data and introduce a third approach. The first approach uses phase-amplitude coupling to determine the RSN. The second approach extracts RSN through an independent component analysis of the Hilbert envelope in different frequency bands, while the third new approach uses a singular value decomposition instead. To evaluate these approaches, we compare the MEG-RSN to the functional magnetic resonance imaging (fMRI)-RSN from the same subjects. Overall, it was possible to extract RSN with MEG using all three techniques, which matched the group-specific fMRI-RSN. Interestingly the new approach based on SVD yielded significantly higher correspondence to five out of seven fMRI-RSN than the two existing approaches. Importantly, with this approach, all networks-except for the visual network-had the highest correspondence to the fMRI networks within one frequency band. Thereby we provide further insights into the electrophysiological underpinnings of the fMRI-RSNs. This knowledge will be important for the analysis of the electrophysiological connectome.

Enhanced brain structure-function tethering in transmodal cortex revealed by high-frequency eigenmodes

While the link between brain structure and function remains an ongoing challenge, the prevailing hypothesis is that the structure-function relationship may itself be gradually decoupling from unimodal to transmodal cortex. However, this hypothesis is constrained by the underlying models which may neglect requisite information. Here we relate structural and functional connectivity derived from diffusion and functional MRI through orthogonal eigenmodes governing frequency-specific diffusion patterns. We find that low-frequency eigenmodes contribute little to functional interactions in transmodal cortex, resulting in divergent structure-function relationships. Conversely, high-frequency eigenmodes predominantly support neuronal coactivation patterns in these areas, inducing structure-function convergence along a unimodal-transmodal hierarchy. High-frequency information, although weak and scattered, could enhance the structure-function tethering, especially in transmodal association cortices. Our findings suggest that the structure-function decoupling may not be an intrinsic property of brain organization, but can be narrowed through multiplexed and regionally specialized spatiotemporal propagation regimes.

Adult lifespan trajectories of neuromagnetic signals and interrelations with cortical thickness

Oscillatory power and phase synchronization map neuronal dynamics and are commonly studied to differentiate the healthy and diseased brain. Yet, little is known about the course and spatial variability of these features from early adulthood into old age. Leveraging magnetoencephalography (MEG) resting-state data in a cross-sectional adult sample (n = 350), we probed lifespan differences (18-88 years) in connectivity and power and interaction effects with sex. Building upon recent attempts to link brain structure and function, we tested the spatial correspondence between age effects on cortical thickness and those on functional networks. We further probed a direct structure-function relationship at the level of the study sample. We found MEG frequency-specific patterns with age and divergence between sexes in low frequencies. Connectivity and power exhibited distinct linear trajectories or turning points at midlife that might reflect different physiological processes. In the delta and beta bands, these age effects corresponded to those on cortical thickness, pointing to co-variation between the modalities across the lifespan. Structure-function coupling was frequency-dependent and observed in unimodal or multimodal regions. Altogether, we provide a comprehensive overview of the topographic functional profile of adulthood that can form a basis for neurocognitive and clinical investigations. This study further sheds new light on how the brain's structural architecture relates to fast oscillatory activity.

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

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

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Connectomics of human electrophysiology

We present both a scientific overview and conceptual positions concerning the challenges and assets of electrophysiological measurements in the search for the nature and functions of the human connectome. We discuss how the field has been inspired by findings and approaches from functional magnetic resonance imaging (fMRI) and informed by a small number of significant multimodal empirical studies, which show that the canonical networks that are commonplace in fMRI are in fact rooted in electrophysiological processes. This review is also an opportunity to produce a brief, up-to-date critical survey of current data modalities and analytical methods available for deriving both static and dynamic connectomes from electrophysiology. We review hurdles that challenge the significance and impact of current electrophysiology connectome research. We then encourage the field to take a leap of faith and embrace the wealth of electrophysiological signals, despite their apparent, disconcerting complexity. Our position is that electrophysiology connectomics is poised to inform testable mechanistic models of information integration in hierarchical brain networks, constructed from observable oscillatory and aperiodic signal components and their polyrhythmic interactions.

Resting state network connectivity is attenuated by fMRI acoustic noise

INTRODUCTION

During the past decades there has been an increasing interest in tracking brain network fluctuations in health and disease by means of resting state functional magnetic resonance imaging (rs-fMRI). Rs-fMRI however does not provide the ideal environmental setting, as participants are continuously exposed to noise generated by MRI coils during acquisition of Echo Planar Imaging (EPI). We investigated the effect of EPI noise on resting state activity and connectivity using magnetoencephalography (MEG), by reproducing the acoustic characteristics of rs-fMRI environment during the recordings. As compared to fMRI, MEG has little sensitivity to brain activity generated in deep brain structures, but has the advantage to capture both the dynamic of cortical magnetic oscillations with high temporal resolution and the slow magnetic fluctuations highly correlated with BOLD signal.

METHODS

Thirty healthy subjects were enrolled in a counterbalanced design study including three conditions: a) silent resting state (Silence), b) resting state upon EPI noise (fMRI), and c) resting state upon white noise (White). White noise was employed to test the specificity of fMRI noise effect. The amplitude envelope correlation (AEC) in alpha band measured the connectivity of seven Resting State Networks (RSN) of interest (default mode network, dorsal attention network, language, left and right auditory and left and right sensory-motor). Vigilance dynamic was estimated from power spectral activity.

RESULTS

fMRI and White acoustic noise consistently reduced connectivity of cortical networks. The effects were widespread, but noise and network specificities were also present. For fMRI noise, decreased connectivity was found in the right auditory and sensory-motor networks. Progressive increase of slow theta-delta activity related to drowsiness was found in all conditions, but was significantly higher for fMRI . Theta-delta significantly and positively correlated with variations of cortical connectivity.

DISCUSSION

rs-fMRI connectivity is biased by unavoidable environmental factors during scanning, which warrant more careful control and improved experimental designs. MEG is free from acoustic noise and allows a sensitive estimation of resting state connectivity in cortical areas. Although underutilized, MEG could overcome issues related to noise during fMRI, in particular when investigation of motor and auditory networks is needed.

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

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

Deficit of Cross-Frequency Integration in Mild Cognitive Impairment and Alzheimer's Disease: A Multilayer Network Approach

BACKGROUND

Studies at specific frequencies have shown abnormalities in brain functional networks among mild cognitive impairment (MCI) and Alzheimer's disease (AD) patients. Previous studies have failed to take into account the possibility that optimal cognitive integration requires interactions between different frequency bands.

PURPOSE

To study whether there is abnormal cross-frequency integration in patients' brains during disease progression.

STUDY TYPE

Retrospective.

POPULATION

Forty-six normal control (NC), 85 patients with MCI, and 31 patients with AD.

FIELD STRENGTH/SEQUENCE

3T.

ASSESSMENT

Multilayer network models were constructed for NC, MCI, and AD, and multilayer participation coefficient (MPC) was used to study the changes of the interlayer relationship in the course of disease development. In addition, MPC and an overlapping degree were combined to classify nodes in the network, and the role of key nodes in the interlayer interaction was mainly observed. Finally, the correlation between multilayer network measures and cognitive function was investigated.

STATISTICAL TESTS

Pearson chi-squared two-tailed test, one-way analysis of variance (ANOVA), nonparametric Spearman correlation coefficient r, and the false discovery rate.

RESULTS

The MPC of the network decreased significantly in MCI (P < 0.05) and AD (P < 0.05). The number of intralayer nodes increased significantly (MCI [P < 0.05], AD [P < 0.05]) and the number of interlayer nodes decreased significantly. Centrality loss between frequencies of a large number of hub nodes, among which the damaged hub nodes included the left hippocampus, left precuneus, right precuneus, left posterior cingulate gyrus, left precentral gyrus, right precentral gyrus, left medial superior frontal gyrus, and right postcentral gyrus. MPC was significantly associated with memory impairment in patients (AD [Spearman's r = 0.526, P < 0.05], MCI [Spearman's r = 0.229, P < 0.05]), and these related regions included damaged hub nodes in patients.

DATA CONCLUSION

In the multilayer networks of patients, there was an obvious deficit in cross-frequency integration and the hub nodes were preferentially damaged. Moreover, these vulnerable hubs are associated with patients' cognitive scores.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY STAGE: 3.

The relationship between EEG and fMRI connectomes is reproducible across simultaneous EEG-fMRI studies from 1.5T to 7T

Both electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are non-invasive methods that show complementary aspects of human brain activity. Despite measuring different proxies of brain activity, both the measured blood-oxygenation (fMRI) and neurophysiological recordings (EEG) are indirectly coupled. The electrophysiological and BOLD signal can map the underlying functional connectivity structure at the whole brain scale at different timescales. Previous work demonstrated a moderate but significant correlation between resting-state functional connectivity of both modalities, however there is a wide range of technical setups to measure simultaneous EEG-fMRI and the reliability of those measures between different setups remains unknown. This is true notably with respect to different magnetic field strengths (low and high field) and different spatial sampling of EEG (medium to high-density electrode coverage). Here, we investigated the reproducibility of the bimodal EEG-fMRI functional connectome in the most comprehensive resting-state simultaneous EEG-fMRI dataset compiled to date including a total of 72 subjects from four different imaging centers. Data was acquired from 1.5T, 3T and 7T scanners with simultaneously recorded EEG using 64 or 256 electrodes. We demonstrate that the whole-brain monomodal connectivity reproducibly correlates across different datasets and that a moderate crossmodal correlation between EEG and fMRI connectivity of r ≈ 0.3 can be reproducibly extracted in low- and high-field scanners. The crossmodal correlation was strongest in the EEG-β frequency band but exists across all frequency bands. Both homotopic and within intrinsic connectivity network (ICN) connections contributed the most to the crossmodal relationship. This study confirms, using a considerably diverse range of recording setups, that simultaneous EEG-fMRI offers a consistent estimate of multimodal functional connectomes in healthy subjects that are dominantly linked through a functional core of ICNs across spanning across the different timescales measured by EEG and fMRI. This opens new avenues for estimating the dynamics of brain function and provides a better understanding of interactions between EEG and fMRI measures. This observed level of reproducibility also defines a baseline for the study of alterations of this coupling in pathological conditions and their role as potential clinical markers.

Dynamics of task-related electrophysiological networks: a benchmarking study

Motor, sensory and cognitive functions rely on dynamic reshaping of functional brain networks. Tracking these rapid changes is crucial to understand information processing in the brain, but challenging due to the great variety of dimensionality reduction methods used at the network-level and the limited evaluation studies. Using Magnetoencephalography (MEG) combined with Source Separation (SS) methods, we present an integrated framework to track fast dynamics of electrophysiological brain networks. We evaluate nine SS methods applied to three independent MEG databases (N=95) during motor and memory tasks. We report differences between these methods at the group and subject level. We seek to help researchers in choosing objectively the appropriate SS method when tracking fast reconfiguration of functional brain networks, due to its enormous benefits in cognitive and clinical neuroscience.

Measuring functional connectivity with wearable MEG

Optically-pumped magnetometers (OPMs) offer the potential for a step change in magnetoencephalography (MEG) enabling wearable systems that provide improved data quality, accommodate any subject group, allow data capture during movement and potentially reduce cost. However, OPM-MEG is a nascent technology and, to realise its potential, it must be shown to facilitate key neuroscientific measurements, such as the characterisation of brain networks. Networks, and the connectivities that underlie them, have become a core area of neuroscientific investigation, and their importance is underscored by many demonstrations of their disruption in brain disorders. Consequently, a demonstration of network measurements using OPM-MEG would be a significant step forward. Here, we aimed to show that a wearable 50-channel OPM-MEG system enables characterisation of the electrophysiological connectome. To this end, we measured connectivity in the resting state and during a visuo-motor task, using both OPM-MEG and a state-of-the-art 275-channel cryogenic MEG device. Our results show that resting-state connectome matrices from OPM and cryogenic systems exhibit a high degree of similarity, with correlation values >70%. In addition, in task data, similar differences in connectivity between individuals (scanned multiple times) were observed in cryogenic and OPM-MEG data, again demonstrating the fidelity of the OPM-MEG device. This is the first demonstration of network connectivity measured using OPM-MEG, and results add weight to the argument that OPMs will ultimately supersede cryogenic sensors for MEG measurement.

Brain dysconnectivity relates to disability and cognitive impairment in multiple sclerosis

The pathophysiology of cognitive dysfunction in multiple sclerosis (MS) is still unclear. This magnetoencephalography (MEG) study investigates the impact of MS on brain resting-state functional connectivity (rsFC) and its relationship to disability and cognitive impairment. We investigated rsFC based on power envelope correlation within and between different frequency bands, in a large cohort of participants consisting of 99 MS patients and 47 healthy subjects. Correlations were investigated between rsFC and outcomes on disability, disease duration and 7 neuropsychological scores within each group, while stringently correcting for multiple comparisons and possible confounding factors. Specific dysconnections correlating with MS-induced physical disability and disease duration were found within the sensorimotor and language networks, respectively. Global network-level reductions in within- and cross-network rsFC were observed in the default-mode network. Healthy subjects and patients significantly differed in their scores on cognitive fatigue and verbal fluency. Healthy subjects and patients showed different correlation patterns between rsFC and cognitive fatigue or verbal fluency, both of which involved a shift in patients from the posterior default-mode network to the language network. Introducing electrophysiological rsFC in a regression model of verbal fluency and cognitive fatigue in MS patients significantly increased the explained variance compared to a regression limited to structural MRI markers (relative thalamic volume and lesion load). This MEG study demonstrates that MS induces distinct changes in the resting-state functional brain architecture that relate to disability, disease duration and specific cognitive functioning alterations. It highlights the potential value of electrophysiological intrinsic rsFC for monitoring the cognitive impairment in patients with MS.

Oscillation-Based Connectivity Architecture Is Dominated by an Intrinsic Spatial Organization, Not Cognitive State or Frequency

Functional connectivity of neural oscillations (oscillation-based FC) is thought to afford dynamic information exchange across task-relevant neural ensembles. Although oscillation-based FC is classically defined relative to a prestimulus baseline, giving rise to rapid, context-dependent changes in individual connections, studies of distributed spatial patterns show that oscillation-based FC is omnipresent, occurring even in the absence of explicit cognitive demands. Thus, the issue of whether oscillation-based FC is primarily shaped by cognitive state or is intrinsic in nature remains open. Accordingly, we sought to reconcile these observations by interrogating the ECoG recordings of 18 presurgical human patients (8 females) for state dependence of oscillation-based FC in five canonical frequency bands across an array of six task states. FC analysis of phase and amplitude coupling revealed a highly similar, largely state-invariant (i.e., intrinsic) spatial component across cognitive states. This spatial organization was shared across all frequency bands. Crucially, however, each band also exhibited temporally independent FC dynamics capable of supporting frequency-specific information exchange. In conclusion, the spatial organization of oscillation-based FC is largely stable over cognitive states (i.e., primarily intrinsic in nature) and shared across frequency bands. Together, our findings converge with previous observations of spatially invariant patterns of FC derived from extremely slow and aperiodic fluctuations in fMRI signals. Our observations indicate that "background" FC should be accounted for in conceptual frameworks of oscillation-based FC targeting task-related changes.SIGNIFICANCE STATEMENT A fundamental property of neural activity is that it is periodic, enabling functional connectivity (FC) between distant regions through coupling of their oscillations. According to task-based studies, such oscillation-based FC is rapid and malleable to meet cognitive task demands. Studying distributed FC patterns instead of FC in a few individual connections, we found that oscillation-based FC is largely stable across various cognitive states and shares a common layout across oscillation frequencies. This stable spatial organization of FC in fast oscillatory brain signals parallels the known stability of fMRI-based intrinsic FC architecture. Despite the observed spatial state and frequency invariance, FC of individual connections was temporally independent between frequency bands, suggesting a putative mechanism for malleable frequency-specific FC to support cognitive tasks.

Relationships Between Neuronal Oscillatory Amplitude and Dynamic Functional Connectivity

Event-related fluctuations of neural oscillatory amplitude are reported widely in the context of cognitive processing and are typically interpreted as a marker of brain "activity". However, the precise nature of these effects remains unclear; in particular, whether such fluctuations reflect local dynamics, integration between regions, or both, is unknown. Here, using magnetoencephalography, we show that movement induced oscillatory modulation is associated with transient connectivity between sensorimotor regions. Further, in resting-state data, we demonstrate a significant association between oscillatory modulation and dynamic connectivity. A confound with such empirical measurements is that increased amplitude necessarily means increased signal-to-noise ratio (SNR): this means that the question of whether amplitude and connectivity are genuinely coupled, or whether increased connectivity is observed purely due to increased SNR is unanswered. Here, we counter this problem by analogy with computational models which show that, in the presence of global network coupling and local multistability, the link between oscillatory modulation and long-range connectivity is a natural consequence of neural networks. Our results provide evidence for the notion that connectivity is mediated by neural oscillations, and suggest that time-frequency spectrograms are not merely a description of local synchrony but also reflect fluctuations in long-range connectivity.

Vascular physiology drives functional brain networks

We present the first evidence for vascular regulation driving fMRI signals in specific functional brain networks. Using concurrent neuronal and vascular stimuli, we collected 30 BOLD fMRI datasets in 10 healthy individuals: a working memory task, flashing checkerboard stimulus, and CO2 inhalation challenge were delivered in concurrent but orthogonal paradigms. The resulting imaging data were averaged together and decomposed using independent component analysis, and three "neuronal networks" were identified as demonstrating maximum temporal correlation with the neuronal stimulus paradigms: Default Mode Network, Task Positive Network, and Visual Network. For each of these, we observed a second network component with high spatial overlap. Using dual regression in the original 30 datasets, we extracted the time-series associated with these network pairs and calculated the percent of variance explained by the neuronal or vascular stimuli using a normalized R2 parameter. In each pairing, one network was dominated by the appropriate neuronal stimulus, and the other was dominated by the vascular stimulus as represented by the end-tidal CO2 time-series recorded in each scan. We acquired a second dataset in 8 of the original participants, where no CO2 challenge was delivered and CO2 levels fluctuated naturally with breathing variations. Although splitting of functional networks was not robust in these data, performing dual regression with the network maps from the original analysis in this new dataset successfully replicated our observations. Thus, in addition to responding to localized metabolic changes, the brain's vasculature may be regulated in a coordinated manner that mimics (and potentially supports) specific functional brain networks. Multi-modal imaging and advances in fMRI acquisition and analysis could facilitate further study of the dual nature of functional brain networks. It will be critical to understand network-specific vascular function, and the behavior of a coupled vascular-neural network, in future studies of brain pathology.

Multilayer MEG functional connectivity as a potential marker for suicidal thoughts in major depressive disorder

Major depressive disorder (MDD) is highly heterogeneous in its clinical presentation. The present exploratory study used magnetoencephalography (MEG) to investigate electrophysiological intrinsic connectivity differences between healthy volunteers and unmedicated participants with treatment-resistant MDD. The study examined canonical frequency bands from delta through gamma. In addition to group comparisons, correlational studies were conducted to determine whether connectivity was related to five symptom factors: depressed mood, tension, negative cognition, suicidal thoughts, and amotivation. The MDD and healthy volunteer groups did not differ significantly at baseline when corrected across all frequencies and clusters, although evidence of generalized slowing in MDD was observed. Notably, however, electrophysiological connectivity was strongly related to suicidal thoughts, particularly as coupling of low frequency power fluctuations (delta and theta) with alpha and beta power. This analysis revealed hub areas underlying this symptom cluster, including left hippocampus, left anterior insula, and bilateral dorsolateral prefrontal cortex. No other symptom cluster demonstrated a relationship with neurophysiological connectivity, suggesting a specificity to these results as markers of suicidal ideation.

Concurrent EEG- and fMRI-derived functional connectomes exhibit linked dynamics

Long-range connectivity has become the most studied feature of human functional Magnetic Resonance Imaging (fMRI), yet the spatial and temporal relationship between its whole-brain dynamics and electrophysiological connectivity remains largely unknown. FMRI-derived functional connectivity exhibits spatial reconfigurations or time-varying dynamics at infraslow (<0.1Hz) speeds. Conversely, electrophysiological connectivity is based on cross-region coupling of fast oscillations (~1-100Hz). It is unclear whether such fast oscillation-based coupling varies at infraslow speeds, temporally coinciding with infraslow dynamics across the fMRI-based connectome. If so, does the association of fMRI-derived and electrophysiological dynamics spatially vary over the connectome across the functionally distinct electrophysiological oscillation bands? In two concurrent electroencephalography (EEG)-fMRI resting-state datasets, oscillation-based coherence in all canonical bands (delta through gamma) indeed reconfigured at infraslow speeds in tandem with fMRI-derived connectivity changes in corresponding region-pairs. Interestingly, irrespective of EEG frequency-band the cross-modal tie of connectivity dynamics comprised a large proportion of connections distributed across the entire connectome. However, there were frequency-specific differences in the relative strength of the cross-modal association. This association was strongest in visual to somatomotor connections for slower EEG-bands, and in connections involving the Default Mode Network for faster EEG-bands. Methodologically, the findings imply that neural connectivity dynamics can be reliably measured by fMRI despite heavy susceptibility to noise, and by EEG despite shortcomings of source reconstruction. Biologically, the findings provide evidence that contrast with known territories of oscillation power, oscillation coupling in all bands slowly reconfigures in a highly distributed manner across the whole-brain connectome.

Rhythmic Network Modulation to Thalamocortical Couplings in Epilepsy

Thalamus interacts with cortical areas, generating oscillations characterized by their rhythm and levels of synchrony. However, little is known of what function the rhythmic dynamic may serve in thalamocortical couplings. This work introduced a general approach to investigate the modulatory contribution of rhythmic scalp network to the thalamo-frontal couplings in juvenile myoclonic epilepsy (JME) and frontal lobe epilepsy (FLE). Here, time-varying rhythmic network was constructed using the adapted directed transfer function between EEG electrodes, and then was applied as a modulator in fMRI-based thalamocortical functional couplings. Furthermore, the relationship between corticocortical connectivity and rhythm-dependent thalamocortical coupling was examined. The results revealed thalamocortical couplings modulated by EEG scalp network have frequency-dependent characteristics. Increased thalamus- sensorimotor network (SMN) and thalamus-default mode network (DMN) couplings in JME were strongly modulated by alpha band. These thalamus-SMN couplings demonstrated enhanced association with SMN-related corticocortical connectivity. In addition, altered theta-dependent and beta-dependent thalamus-frontoparietal network (FPN) couplings were found in FLE. The reduced theta-dependent thalamus-FPN couplings were associated with the decreased FPN-related corticocortical connectivity. This study proposed interactive links between the rhythmic modulation and thalamocortical coupling. The crucial role of SMN and FPN in subcortical-cortical circuit may have implications for intervention in generalized and focal epilepsy.

Intrinsic connectome organization across temporal scales: New insights from cross-modal approaches

The discovery of a stable, whole-brain functional connectivity organization that is largely independent of external events has drastically extended our view of human brain function. However, this discovery has been primarily based on functional magnetic resonance imaging (fMRI). The role of this whole-brain organization in fast oscillation-based connectivity as measured, for example, by electroencephalography (EEG) and magnetoencephalography (MEG) is only beginning to emerge. Here, we review studies of intrinsic connectivity and its whole-brain organization in EEG, MEG, and intracranial electrophysiology with a particular focus on direct comparisons to connectome studies in fMRI. Synthesizing this literature, we conclude that irrespective of temporal scale over four orders of magnitude, intrinsic neurophysiological connectivity shows spatial similarity to the connectivity organization commonly observed in fMRI. A shared structural connectivity basis and cross-frequency coupling are possible mechanisms contributing to this similarity. Acknowledging that a stable whole-brain organization governs long-range coupling across all timescales of neural processing motivates researchers to take "baseline" intrinsic connectivity into account when investigating brain-behavior associations, and further encourages more widespread exploration of functional connectomics approaches beyond fMRI by using EEG and MEG modalities.

The role of transient spectral 'bursts' in functional connectivity: A magnetoencephalography study

Neural oscillations dominate electrophysiological measures of macroscopic brain activity and fluctuations in these rhythms offer an insightful window on cortical excitation, inhibition, and connectivity. However, in recent years the 'classical' picture of smoothly varying oscillations has been challenged by the idea that many 'oscillations' may actually be formed from the recurrence of punctate high-amplitude bursts in activity, whose spectral composition intersects the traditionally defined frequency ranges (e.g. alpha/beta band). This finding offers a new interpretation of measurable brain activity, however neither the methodological means to detect bursts, nor their link to other findings (e.g. connectivity) have been settled. Here, we use a new approach to detect bursts in magnetoencephalography (MEG) data. We show that a time-delay embedded Hidden Markov Model (HMM) can be used to delineate single-region bursts which are in agreement with existing techniques. However, unlike existing techniques, the HMM looks for specific spectral patterns in timecourse data. We characterise the distribution of burst duration, frequency of occurrence and amplitude across the cortex in resting state MEG data. During a motor task we show how the movement related beta decrease and post movement beta rebound are driven by changes in burst occurrence. Finally, we show that the beta band functional connectome can be derived using a simple measure of burst overlap, and that coincident bursts in separate regions correspond to a period of heightened coherence. In summary, this paper offers a new methodology for burst identification and connectivity analysis which will be important for future investigations of neural oscillations.

Multifractal and entropy analysis of resting-state electroencephalography reveals spatial organization in local dynamic functional connectivity

Functional connectivity of the brain fluctuates even in resting-state condition. It has been reported recently that fluctuations of global functional network topology and those of individual connections between brain regions expressed multifractal scaling. To expand on these findings, in this study we investigated if multifractality was indeed an inherent property of dynamic functional connectivity (DFC) on the regional level as well. Furthermore, we explored if local DFC showed region-specific differences in its multifractal and entropy-related features. DFC analyses were performed on 62-channel, resting-state electroencephalography recordings of twelve young, healthy subjects. Surrogate data testing verified the true multifractal nature of regional DFC that could be attributed to the presumed nonlinear nature of the underlying processes. Moreover, we found a characteristic spatial distribution of local connectivity dynamics, in that frontal and occipital regions showed stronger long-range correlation and higher degree of multifractality, whereas the highest values of entropy were found over the central and temporal regions. The revealed topology reflected well the underlying resting-state network organization of the brain. The presented results and the proposed analysis framework could improve our understanding on how resting-state brain activity is spatio-temporally organized and may provide potential biomarkers for future clinical research.

How Sensitive Are Conventional MEG Functional Connectivity Metrics With Sliding Windows to Detect Genuine Fluctuations in Dynamic Functional Connectivity?

Despite advances in the field of dynamic connectivity, fixed sliding window approaches for the detection of fluctuations in functional connectivity are still widely used. The use of conventional connectivity metrics in conjunction with a fixed sliding window comes with the arbitrariness of the chosen window lengths. In this paper we use multivariate autoregressive and neural mass models with a priori defined ground truths to systematically analyze the sensitivity of conventional metrics in combination with different window lengths to detect genuine fluctuations in connectivity for various underlying state durations. Metrics of interest are the coherence, imaginary coherence, phase lag index, phase locking value and the amplitude envelope correlation. We performed analysis for two nodes and at the network level. We demonstrate that these metrics show indeed higher variability for genuine temporal fluctuations in connectivity compared to a static connectivity state superimposed by noise. Overall, the error of the connectivity estimates themselves decreases for longer state durations (order of seconds), while correlations of the connectivity fluctuations with the ground truth was higher for longer state durations. In general, metrics, in combination with a sliding window, perform poorly for very short state durations. Increasing the SNR of the system only leads to a moderate improvement. In addition, at the network level, only longer window widths were sufficient to detect plausible resting state networks that matched the underlying ground truth, especially for the phase locking value, amplitude envelope correlation and coherence. The length of these longer window widths did not necessarily correspond to the underlying state durations. For short window widths resting state network connectivity patterns could not be retrieved. We conclude that fixed sliding window approaches for connectivity can detect modulations of connectivity, but mostly if the underlying dynamics operate on moderate to slow timescales. In practice, this can be a drawback, as state durations can vary significantly in empirical data.

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.

Spatial and spectral trajectories in typical neurodevelopment from childhood to middle age

Detailed characterization of typical human neurodevelopment is key if we are to understand the nature of mental and neurological pathology. While research on the cellular processes of neurodevelopment has made great advances, in vivo human imaging is crucial to understand our uniquely human capabilities, as well as the pathologies that affect them. Using magnetoencephalography data in the largest normative sample currently available (324 participants aged 6-45 years), we assess the developmental trajectory of resting-state oscillatory power and functional connectivity from childhood to middle age. The maturational course of power, indicative of local processing, was found to both increase and decrease in a spectrally dependent fashion. Using the strength of phase-synchrony between parcellated regions, we found significant linear and nonlinear (quadratic and logarithmic) trajectories to be characterized in a spatially heterogeneous frequency-specific manner, such as a superior frontal region with linear and nonlinear trajectories in theta and gamma band respectively. Assessment of global efficiency revealed similar significant nonlinear trajectories across all frequency bands. Our results link with the development of human cognitive abilities; they also highlight the complexity of neurodevelopment and provide quantitative parameters for replication and a robust footing from which clinical research may map pathological deviations from these typical trajectories.

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.

Altered temporal stability in dynamic neural networks underlies connectivity changes in neurodevelopment

Network connectivity is an integral feature of human brain function, and characterising its maturational trajectory is a critical step towards understanding healthy and atypical neurodevelopment. Here, we used magnetoencephalography (MEG) to investigate both stationary (i.e. time averaged) and rapidly modulating (dynamic) electrophysiological connectivity, in participants aged from mid-childhood to early adulthood (youngest participant 9 years old; oldest participant 25 years old). Stationary functional connectivity (measured via inter-regional coordination of neural oscillations) increased with age in the alpha and beta frequency bands, particularly in bilateral parietal and temporo-parietal connections. Our dynamic analysis (also applied to alpha/beta oscillations) revealed the spatiotemporal signatures of 8 dynamic networks; these modulate on a ∼100 ms time scale, and temporal stability in attentional networks was found to increase with age. Significant overlap was found between age-modulated dynamic networks and inter-regional oscillatory coordination, implying that altered network dynamics underlie age related changes in functional connectivity. Our results provide novel insights into brain network electrophysiology, and lay a foundation for future work in childhood disorders.

Interpreting BOLD: towards a dialogue between cognitive and cellular neuroscience

Cognitive neuroscience depends on the use of blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) to probe brain function. Although commonly used as a surrogate measure of neuronal activity, BOLD signals actually reflect changes in brain blood oxygenation. Understanding the mechanisms linking neuronal activity to vascular perfusion is, therefore, critical in interpreting BOLD. Advances in cellular neuroscience demonstrating differences in this neurovascular relationship in different brain regions, conditions or pathologies are often not accounted for when interpreting BOLD. Meanwhile, within cognitive neuroscience, the increasing use of high magnetic field strengths and the development of model-based tasks and analyses have broadened the capability of BOLD signals to inform us about the underlying neuronal activity, but these methods are less well understood by cellular neuroscientists. In 2016, a Royal Society Theo Murphy Meeting brought scientists from the two communities together to discuss these issues. Here, we consolidate the main conclusions arising from that meeting. We discuss areas of consensus about what BOLD fMRI can tell us about underlying neuronal activity, and how advanced modelling techniques have improved our ability to use and interpret BOLD. We also highlight areas of controversy in understanding BOLD and suggest research directions required to resolve these issues.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Adaptive cortical parcellations for source reconstructed EEG/MEG connectomes

There is growing interest in the rich temporal and spectral properties of the functional connectome of the brain that are provided by Electro- and Magnetoencephalography (EEG/MEG). However, the problem of leakage between brain sources that arises when reconstructing brain activity from EEG/MEG recordings outside the head makes it difficult to distinguish true connections from spurious connections, even when connections are based on measures that ignore zero-lag dependencies. In particular, standard anatomical parcellations for potential cortical sources tend to over- or under-sample the real spatial resolution of EEG/MEG. By using information from cross-talk functions (CTFs) that objectively describe leakage for a given sensor configuration and distributed source reconstruction method, we introduce methods for optimising the number of parcels while simultaneously minimising the leakage between them. More specifically, we compare two image segmentation algorithms: 1) a split-and-merge (SaM) algorithm based on standard anatomical parcellations and 2) a region growing (RG) algorithm based on all the brain vertices with no prior parcellation. Interestingly, when applied to minimum-norm reconstructions for EEG/MEG configurations from real data, both algorithms yielded approximately 70 parcels despite their different starting points, suggesting that this reflects the resolution limit of this particular sensor configuration and reconstruction method. Importantly, when compared against standard anatomical parcellations, resolution matrices of adaptive parcellations showed notably higher sensitivity and distinguishability of parcels. Furthermore, extensive simulations of realistic networks revealed significant improvements in network reconstruction accuracies, particularly in reducing false leakage-induced connections. Adaptive parcellations therefore allow a more accurate reconstruction of functional EEG/MEG connectomes.

•

We introduce adaptive cortical parcellation algorithms for E/MEG source estimation.

•

Optimum number, size and locations of parcels are found based on cross-talk functions

•

Algorithms yielded ∼70 distinguishable parcels regardless of the starting point.

•

Parcel resolution matrices were notably improved compared to anatomical atlases.

•

Network reconstruction accuracies of simulated connectomes improved significantly.

Complementary contributions of concurrent EEG and fMRI connectivity for predicting structural connectivity

While averaged dynamics of brain function are known to estimate the underlying structure, the exact relationship between large-scale function and structure remains an unsolved issue in network neuroscience. These complex functional dynamics, measured by EEG and fMRI, are thought to arise from a shared underlying structural architecture, which can be measured by diffusion MRI (dMRI). While simulation and data transformation (e.g. graph theory measures) have been proposed to refine the understanding of the underlying function-structure relationship, the potential complementary and/or independent contribution of EEG and fMRI to this relationship is still poorly understood. As such, we explored this relationship by analyzing the function-structure correlation in fourteen healthy subjects with simultaneous resting-state EEG-fMRI and dMRI acquisitions. We show that the combination of EEG and fMRI connectivity better explains dMRI connectivity and that this represents a genuine model improvement over fMRI-only models for both group-averaged connectivity matrices and at the individual level. Furthermore, this model improves the prediction within each resting-state network. The best model fit to underlying structure is mediated by fMRI and EEG-δ connectivity in combination with Euclidean distance and interhemispheric connectivity with more local contributions of EEG-γ at the scale of resting-state networks. This highlights that the factors mediating the relationship between functional and structural metrics of connectivity are context and scale dependent, influenced by topological, geometric and architectural features. It also suggests that fMRI studies employing simultaneous EEG measures may characterize additional and essential parts of the underlying neuronal activity of the resting-state, which might be of special interest for both clinical studies and the investigation of resting-state dynamics.

Optimising experimental design for MEG resting state functional connectivity measurement

The study of functional connectivity using magnetoencephalography (MEG) is an expanding area of neuroimaging, and adds an extra dimension to the more common assessments made using fMRI. The importance of such metrics is growing, with recent demonstrations of their utility in clinical research, however previous reports suggest that whilst group level resting state connectivity is robust, single session recordings lack repeatability. Such robustness is critical if MEG measures in individual subjects are to prove clinically valuable. In the present paper, we test how practical aspects of experimental design affect the intra-subject repeatability of MEG findings; specifically we assess the effect of co-registration method and data recording duration. We show that the use of a foam head-cast, which is known to improve co-registration accuracy, increased significantly the between session repeatability of both beamformer reconstruction and connectivity estimation. We also show that recording duration is a critical parameter, with large improvements in repeatability apparent when using ten minute, compared to five minute recordings. Further analyses suggest that the origin of this latter effect is not underpinned by technical aspects of source reconstruction, but rather by a genuine effect of brain state; short recordings are simply inefficient at capturing the canonical MEG network in a single subject. Our results provide important insights on experimental design and will prove valuable for future MEG connectivity studies.

Integrating cross-frequency and within band functional networks in resting-state MEG: A multi-layer network approach

Neuronal oscillations exist across a broad frequency spectrum, and are thought to provide a mechanism of interaction between spatially separated brain regions. Since ongoing mental activity necessitates the simultaneous formation of multiple networks, it seems likely that the brain employs interactions within multiple frequency bands, as well as cross-frequency coupling, to support such networks. Here, we propose a multi-layer network framework that elucidates this pan-spectral picture of network interactions. Our network consists of multiple layers (frequency-band specific networks) that influence each other via inter-layer (cross-frequency) coupling. Applying this model to MEG resting-state data and using envelope correlations as connectivity metric, we demonstrate strong dependency between within layer structure and inter-layer coupling, indicating that networks obtained in different frequency bands do not act as independent entities. More specifically, our results suggest that frequency band specific networks are characterised by a common structure seen across all layers, superimposed by layer specific connectivity, and inter-layer coupling is most strongly associated with this common mode. Finally, using a biophysical model, we demonstrate that there are two regimes of multi-layer network behaviour; one in which different layers are independent and a second in which they operate highly dependent. Results suggest that the healthy human brain operates at the transition point between these regimes, allowing for integration and segregation between layers. Overall, our observations show that a complete picture of global brain network connectivity requires integration of connectivity patterns across the full frequency spectrum.

Relationships between cortical myeloarchitecture and electrophysiological networks

The human brain relies upon the dynamic formation and dissolution of a hierarchy of functional networks to support ongoing cognition. However, how functional connectivities underlying such networks are supported by cortical microstructure remains poorly understood. Recent animal work has demonstrated that electrical activity promotes myelination. Inspired by this, we test a hypothesis that gray-matter myelin is related to electrophysiological connectivity. Using ultra-high field MRI and the principle of structural covariance, we derive a structural network showing how myelin density differs across cortical regions and how separate regions can exhibit similar myeloarchitecture. Building upon recent evidence that neural oscillations mediate connectivity, we use magnetoencephalography to elucidate networks that represent the major electrophysiological pathways of communication in the brain. Finally, we show that a significant relationship exists between our functional and structural networks; this relationship differs as a function of neural oscillatory frequency and becomes stronger when integrating oscillations over frequency bands. Our study sheds light on the way in which cortical microstructure supports functional networks. Further, it paves the way for future investigations of the gray-matter structure/function relationship and its breakdown in pathology.

Measurement of dynamic task related functional networks using MEG

The characterisation of dynamic electrophysiological brain networks, which form and dissolve in order to support ongoing cognitive function, is one of the most important goals in neuroscience. Here, we introduce a method for measuring such networks in the human brain using magnetoencephalography (MEG). Previous network analyses look for brain regions that share a common temporal profile of activity. Here distinctly, we exploit the high spatio-temporal resolution of MEG to measure the temporal evolution of connectivity between pairs of parcellated brain regions. We then use an ICA based procedure to identify networks of connections whose temporal dynamics covary. We validate our method using MEG data recorded during a finger movement task, identifying a transient network of connections linking somatosensory and primary motor regions, which modulates during the task. Next, we use our method to image the networks which support cognition during a Sternberg working memory task. We generate a novel neuroscientific picture of cognitive processing, showing the formation and dissolution of multiple networks which relate to semantic processing, pattern recognition and language as well as vision and movement. Our method tracks the dynamics of functional connectivity in the brain on a timescale commensurate to the task they are undertaking.

•

A method is developed to track dynamic electrophysiological networks using MEG.

•

Method based on ICA applied to timecourses measuring evolution of connectivity.

•

Method allows a unique picture of transient networks that support cognition.

•

Method validated in MEG data recorded during a Sternberg working memory task.

•

Sensory networks observed include visual and sensorimotor.

•

Cognitive networks relate to semantic processing, pattern recognition and language.

A multi-layer network approach to MEG connectivity analysis

Recent years have shown the critical importance of inter-regional neural network connectivity in supporting healthy brain function. Such connectivity is measurable using neuroimaging techniques such as MEG, however the richness of the electrophysiological signal makes gaining a complete picture challenging. Specifically, connectivity can be calculated as statistical interdependencies between neural oscillations within a large range of different frequency bands. Further, connectivity can be computed between frequency bands. This pan-spectral network hierarchy likely helps to mediate simultaneous formation of multiple brain networks, which support ongoing task demand. However, to date it has been largely overlooked, with many electrophysiological functional connectivity studies treating individual frequency bands in isolation. Here, we combine oscillatory envelope based functional connectivity metrics with a multi-layer network framework in order to derive a more complete picture of connectivity within and between frequencies. We test this methodology using MEG data recorded during a visuomotor task, highlighting simultaneous and transient formation of motor networks in the beta band, visual networks in the gamma band and a beta to gamma interaction. Having tested our method, we use it to demonstrate differences in occipital alpha band connectivity in patients with schizophrenia compared to healthy controls. We further show that these connectivity differences are predictive of the severity of persistent symptoms of the disease, highlighting their clinical relevance. Our findings demonstrate the unique potential of MEG to characterise neural network formation and dissolution. Further, we add weight to the argument that dysconnectivity is a core feature of the neuropathology underlying schizophrenia.