Using variance information in magnetoencephalography measures of functional connectivity

Spatiotemporal brain hierarchies of auditory memory recognition and predictive coding

Our brain is constantly extracting, predicting, and recognising key spatiotemporal features of the physical world in order to survive. While neural processing of visuospatial patterns has been extensively studied, the hierarchical brain mechanisms underlying conscious recognition of auditory sequences and the associated prediction errors remain elusive. Using magnetoencephalography (MEG), we describe the brain functioning of 83 participants during recognition of previously memorised musical sequences and systematic variations. The results show feedforward connections originating from auditory cortices, and extending to the hippocampus, anterior cingulate gyrus, and medial cingulate gyrus. Simultaneously, we observe backward connections operating in the opposite direction. Throughout the sequences, the hippocampus and cingulate gyrus maintain the same hierarchical level, except for the final tone, where the cingulate gyrus assumes the top position within the hierarchy. The evoked responses of memorised sequences and variations engage the same hierarchical brain network but systematically differ in terms of temporal dynamics, strength, and polarity. Furthermore, induced-response analysis shows that alpha and beta power is stronger for the variations, while gamma power is enhanced for the memorised sequences. This study expands on the predictive coding theory by providing quantitative evidence of hierarchical brain mechanisms during conscious memory and predictive processing of auditory sequences.

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.

Energy landscape of resting magnetoencephalography reveals fronto-parietal network impairments in epilepsy

Juvenile myoclonic epilepsy (JME) is a form of idiopathic generalized epilepsy. It is yet unclear to what extent JME leads to abnormal network activation patterns. Here, we characterized statistical regularities in magnetoencephalograph (MEG) resting-state networks and their differences between JME patients and controls by combining a pairwise maximum entropy model (pMEM) and novel energy landscape analyses for MEG. First, we fitted the pMEM to the MEG oscillatory power in the front-oparietal network (FPN) and other resting-state networks, which provided a good estimation of the occurrence probability of network states. Then, we used energy values derived from the pMEM to depict an energy landscape, with a higher energy state corresponding to a lower occurrence probability. JME patients showed fewer local energy minima than controls and had elevated energy values for the FPN within the theta, beta, and gamma bands. Furthermore, simulations of the fitted pMEM showed that the proportion of time the FPN was occupied within the basins of energy minima was shortened in JME patients. These network alterations were highlighted by significant classification of individual participants employing energy values as multivariate features. Our findings suggested that JME patients had altered multistability in selective functional networks and frequency bands in the fronto-parietal cortices.

Aberrant Whole-Brain Transitions and Dynamics of Spontaneous Network Microstates in Mild Traumatic Brain Injury

Dynamic Functional Connectivity (DFC) analysis is a promising approach for the characterization of brain electrophysiological activity. In this study, we investigated abnormal alterations due to mild Traumatic Brain Injury (mTBI) using DFC of the source reconstructed magnetoencephalographic (MEG) resting-state recordings. Brain activity in several well-known frequency bands was first reconstructed using beamforming of the MEG data to determine ninety anatomical brain regions of interest. A DFC graph was formulated using the imaginary part of phase-locking values, which were obtained from 30 mTBI patients and 50 healthy controls (HC). Subsequently, we estimated normalized Laplacian transformations of individual, statistically and topologically filtered quasi-static graphs. The corresponding eigenvalues of each node synchronization were then computed and through the neural-gas algorithm, we quantized the evolution of the eigenvalues resulting in distinct network microstates (NMstates). The discrimination level between the two groups was assessed using an iterative cross-validation classification scheme with features either the NMstates in each frequency band, or the combination of the so-called chronnectomics (flexibility index, occupancy time of NMstate, and Dwell time) with the complexity index over the evolution of the NMstates across all frequency bands. Classification performance based on chronnectomics showed 80% accuracy, 99% sensitivity, and 49% specificity. However, performance was much higher (accuracy: 91-97%, sensitivity: 100%, and specificity: 77-93%) when focusing on the microstates. Exploring the mean node degree within and between brain anatomical networks (default mode network, frontoparietal, occipital, cingulo-opercular, and sensorimotor), a reduced pattern occurred from lower to higher frequency bands, with statistically significant stronger degrees for the HC than the mTBI group. A higher entropic profile on the temporal evolution of the modularity index was observed for both NMstates for the mTBI group across frequencies. A significant difference in the flexibility index was observed between the two groups for the β frequency band. The latter finding may support a central role of the thalamus impairment in mTBI. The current study considers a complete set of frequency-dependent connectomic markers of mTBI-caused alterations in brain connectivity that potentially could serve as markers to assess the return of an injured subject back to normality.

Synchrony, metastability, dynamic integration, and competition in the spontaneous functional connectivity of the human brain

The human brain is functionally organized into large-scale neural networks that are dynamically interconnected. Multiple short-lived states of resting-state functional connectivity (rsFC) identified transiently synchronized networks and cross-network integration. However, little is known about the way brain couplings covary as rsFC states wax and wane. In this magnetoencephalography study, we explore the synchronization structure among the spontaneous interactions of well-known resting-state networks (RSNs). To do so, we extracted modes of dynamic coupling that reflect rsFC synchrony and analyzed their spatio-temporal features. These modes identified transient, sporadic rsFC changes characterized by the widespread integration of RSNs across the brain, most prominently in the β band. This is in line with the metastable rsFC state model of resting-state dynamics, wherein our modes fit as state transition processes. Furthermore, the default-mode network (DMN) stood out as being structured into competitive cross-network couplings with widespread DMN-RSN interactions, especially among the β-band modes. These results substantiate the theory that the DMN is a core network enabling dynamic global brain integration in the β band.

Dynamics of large-scale electrophysiological networks: A technical review

For several years it has been argued that neural synchronisation is crucial for cognition. The idea that synchronised temporal patterns between different neural groups carries information above and beyond the isolated activity of these groups has inspired a shift in focus in the field of functional neuroimaging. Specifically, investigation into the activation elicited within certain regions by some stimulus or task has, in part, given way to analysis of patterns of co-activation or functional connectivity between distal regions. Recently, the functional connectivity community has been looking beyond the assumptions of stationarity that earlier work was based on, and has introduced methods to incorporate temporal dynamics into the analysis of connectivity. In particular, non-invasive electrophysiological data (magnetoencephalography/electroencephalography (MEG/EEG)), which provides direct measurement of whole-brain activity and rich temporal information, offers an exceptional window into such (potentially fast) brain dynamics. In this review, we discuss challenges, solutions, and a collection of analysis tools that have been developed in recent years to facilitate the investigation of dynamic functional connectivity using these imaging modalities. Further, we discuss the applications of these approaches in the study of cognition and neuropsychiatric disorders. Finally, we review some existing developments that, by using realistic computational models, pursue a deeper understanding of the underlying causes of non-stationary connectivity.

A Comprehensive Review of Magnetoencephalography (MEG) Studies for Brain Functionality in Healthy Aging and Alzheimer's Disease (AD)

Neural oscillations were established with their association with neurophysiological activities and the altered rhythmic patterns are believed to be linked directly to the progression of cognitive decline. Magnetoencephalography (MEG) is a non-invasive technique to record such neuronal activity due to excellent temporal and fair amount of spatial resolution. Single channel, connectivity as well as brain network analysis using MEG data in resting state and task-based experiments were analyzed from existing literature. Single channel analysis studies reported a less complex, more regular and predictable oscillations in Alzheimer's disease (AD) primarily in the left parietal, temporal and occipital regions. Investigations on both functional connectivity (FC) and effective (EC) connectivity analysis demonstrated a loss of connectivity in AD compared to healthy control (HC) subjects found in higher frequency bands. It has been reported from multiplex network of MEG study in AD in the affected regions of hippocampus, posterior default mode network (DMN) and occipital areas, however, conclusions cannot be drawn due to limited availability of clinical literature. Potential utilization of high spatial resolution in MEG likely to provide information related to in-depth brain functioning and underlying factors responsible for changes in neuronal waves in AD. This review is a comprehensive report to investigate diagnostic biomarkers for AD may be identified by from MEG data. It is also important to note that MEG data can also be utilized for the same pursuit in combination with other imaging modalities.

Reliability of Static and Dynamic Network Metrics in the Resting-State: A MEG-Beamformed Connectivity Analysis

The resting activity of the brain can be described by so-called intrinsic connectivity networks (ICNs), which consist of spatially and temporally distributed, but functionally connected, nodes. The coordinated activity of the resting state can be explored via magnetoencephalography (MEG) by studying frequency-dependent functional brain networks at the source level. Although many algorithms for the analysis of brain connectivity have been proposed, the reliability of network metrics derived from both static and dynamic functional connectivity is still unknown. This is a particular problem for studies of associations between ICN metrics and personality variables or other traits, and for studies of differences between patient and control groups, which both depend critically on the reliability of the metrics used. A detailed investigation of the reliability of metrics derived from resting-state MEG repeat scans is therefore a prerequisite for the development of connectomic biomarkers. Here, we first estimated both static (SFC) and dynamic functional connectivity (DFC) after beamforming source reconstruction using the imaginary part of the phase locking index (iPLV) and the correlation of the amplitude envelope (CorEnv). Using our approach, functional network microstates (FCμstates) were derived from the DFC and chronnectomics were computed from the evolution of FCμstates across experimental time. In both temporal scales, the reliability of network metrics (SFC), the FCμstates and the related chronnectomics were evaluated for every frequency band. Chronnectomic statistics and FCμstates were generally more reliable than node-wise static network metrics. CorEnv-based network metrics were more reproducible at the static approach. The reliability of chronnectomics have been evaluated also in a second dataset. This study encourages the analysis of MEG resting-state via DFC.

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.

Resting-state Functional Connectivity is an Age-dependent Predictor of Motor Learning Abilities

This magnetoencephalography study investigates how ageing modulates the relationship between pre-learning resting-state functional connectivity (rsFC) and subsequent learning. Neuromagnetic resting-state activity was recorded 5 min before motor sequence learning in 14 young (19-30 years) and 14 old (66-70 years) participants. We used a seed-based beta-band power envelope correlation approach to estimate rsFC maps, with the seed located in the right primary sensorimotor cortex. In each age group, the relation between individual rsFC and learning performance was investigated using Pearson's correlation analyses. Our results show that rsFC is predictive of subsequent motor sequence learning but involves different cross-network interactions in the two age groups. In young adults, decreased coupling between the sensorimotor network and the cortico-striato-cerebellar network is associated with better motor learning, whereas a similar relation is found in old adults between the sensorimotor, the dorsal-attentional and the DMNs. Additionally, age-related correlational differences were found in the dorsolateral prefrontal cortex, known to subtend attentional and controlled processes. These findings suggest that motor skill learning depends-in an age-dependent manner-on subtle interactions between resting-state networks subtending motor activity on the one hand, and controlled and attentional processes on the other hand.

The effects of AMPA receptor blockade on resting magnetoencephalography recordings

The ionotropic N-methyl-D-aspartate and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors of the glutamatergic neurotransmitter system are of fundamental importance to healthy brain function. Neuroimaging studies in humans have previously been conducted using various drugs that interact with N-methyl-D-aspartate glutamate receptors, but no such studies have investigated AMPA receptor signalling. The recent approval of perampanel (Fycompa) for use in humans provides a means to specifically study the role of AMPA receptors in the pharmacological basis of neuroimaging signals. Twenty male subjects participated in this placebo-controlled crossover study that consisted of two study days separated by a minimum two-week washout period. On one occasion participants ingested a 6 mg dose of perampanel, and on the other a placebo. Ten minutes of wakeful rest was recorded before and after each dose using magnetoencephalography. Subjective ratings of intoxication were significantly higher following drug than placebo. Cluster-based randomisation testing of sensor-level magnetoencephalography data showed significant drug-induced increases in low frequency power (1-4 Hz, 4-8 Hz, 8-13 Hz, 13-30 Hz), along with a significant decrease in the high gamma range (50-90 Hz). We also observed selective increases in functional connectivity in the alpha and beta bands. The findings are consistent with preclinical work and are similar to the spectral profile of other anti-epileptic drugs.

Alzheimer's disease disrupts alpha and beta-band resting-state oscillatory network connectivity

OBJECTIVE

Neuroimaging studies in Alzheimer's disease (AD) yield conflicting results due to selective investigation. We conducted a comprehensive magnetoencephalography study of connectivity changes in AD and healthy ageing in the resting-state.

METHODS

We performed a whole-brain, source-space assessment of oscillatory neural signalling in multiple frequencies comparing AD patients, elderly and young controls. We compared eyes-open and closed group oscillatory envelope activity in networks obtained through temporal independent component analysis, and calculated whole-brain node-based amplitude and phase connectivity.

RESULTS

In bilateral parietotemporal areas, oscillatory envelope amplitude increased with healthy ageing, whereas both local amplitude and node-to-global connectivity decreased with AD. AD-related decreases were spatially specific and restricted to the alpha and beta bands. A significant proportion of the variance in areas of peak group difference was explained by cognitive integrity, in addition to group. None of the groups differed in phase connectivity. Results were highly similar for eyes-open and closed resting-state.

CONCLUSIONS

These results support the disconnection syndrome hypothesis and suggest that AD shows distinct and unique patterns of disrupted neural functioning, rather than accelerated healthy ageing.

SIGNIFICANCE

Whole-brain assessments show that disrupted regional oscillatory envelope amplitude and connectivity in the alpha and beta bands play a key role in AD.

Assessment and elimination of the effects of head movement on MEG resting-state measures of oscillatory brain activity

Magnetoencephalography (MEG) is increasingly being used to study brain function because of its excellent temporal resolution and its direct association with brain activity at the neuronal level. One possible cause of error in the analysis of MEG data comes from the fact that participants, even MEG-experienced ones, move their head in the MEG system. Head movement can cause source localization errors during the analysis of MEG data, which can result in the appearance of source variability that does not reflect brain activity. The MEG community places great importance in eliminating this source of possible errors as is evident, for example, by recent efforts to develop head casts that limit head movement in the MEG system. In this work we use software tools to identify, assess and eliminate from the analysis of MEG data any possible correlations between head movement in the MEG system and widely-used measures of brain activity derived from MEG resting-state recordings. The measures of brain activity we study are a) the Hilbert-transform derived amplitude envelope of the beamformer time series and b) functional networks; both measures derived by MEG resting-state recordings. Ten-minute MEG resting-state recordings were performed on healthy participants, with head position continuously recorded. The sources of the measured magnetic signals were localized via beamformer spatial filtering. Temporal independent component analysis was subsequently used to derive resting-state networks. Significant correlations were observed between the beamformer envelope time series and head movement. The correlations were substantially reduced, and in some cases eliminated, after a participant-specific temporal high-pass filter was applied to those time series. Regressing the head movement metrics out of the beamformer envelope time series had an even stronger effect in reducing these correlations. Correlation trends were also observed between head movement and the activation time series of the default-mode and frontal networks. Regressing the head movement metrics out of the beamformer envelope time series completely eliminated these correlations. Additionally, applying the head movement correction resulted in changes in the network spatial maps for the visual and sensorimotor networks. Our results a) show that the results of MEG resting-state studies that use the above-mentioned analysis methods are confounded by head movement effects, b) suggest that regressing the head movement metrics out of the beamformer envelope time series is a necessary step to be added to these analyses, in order to eliminate the effect that head movement has on the amplitude envelope of beamformer time series and the network time series and c) highlight changes in the connectivity spatial maps when head movement correction is applied.

Abnormal task driven neural oscillations in multiple sclerosis: A visuomotor MEG study

Multiple sclerosis (MS) is a debilitating disease commonly attributed to degradation of white matter myelin. Symptoms include fatigue, as well as problems associated with vision and movement. Although areas of demyelination in white matter are observed routinely in patients undergoing MRI scans, such measures are often a poor predictor of disease severity. For this reason, it is instructive to measure associated changes in brain function. Widespread white‐matter demyelination may lead to delays of propagation of neuronal activity, and with its excellent temporal resolution, magnetoencephalography can be used to probe such delays in controlled conditions (e.g., during a task). In healthy subjects, responses to visuomotor tasks are well documented: in motor cortex, movement elicits a localised decrease in the power of beta band oscillations (event‐related beta desynchronisation) followed by an increase above baseline on movement cessation (post‐movement beta rebound (PMBR)). In visual cortex, visual stimulation generates increased gamma oscillations. In this study, we use a visuomotor paradigm to measure these responses in MS patients and compare them to age‐ and gender‐matched healthy controls. We show a significant increase in the time‐to‐peak of the PMBR in patients which correlates significantly with the symbol digit modalities test: a measure of information processing speed. A significant decrease in the amplitude of visual gamma oscillations in patients is also seen. These findings highlight the potential value of electrophysiological imaging in generating a new understanding of visual disturbances and abnormal motor control in MS patients. Hum Brain Mapp 38:2441–2453, 2017. © 2017 Wiley Periodicals, Inc.

Multimodal Classification of Schizophrenia Patients with MEG and fMRI Data Using Static and Dynamic Connectivity Measures

Mental disorders like schizophrenia are currently diagnosed by physicians/psychiatrists through clinical assessment and their evaluation of patient's self-reported experiences as the illness emerges. There is great interest in identifying biological markers of prognosis at the onset of illness, rather than relying on the evolution of symptoms across time. Functional network connectivity, which indicates a subject's overall level of "synchronicity" of activity between brain regions, demonstrates promise in providing individual subject predictive power. Many previous studies reported functional connectivity changes during resting-state using only functional magnetic resonance imaging (fMRI). Nevertheless, exclusive reliance on fMRI to generate such networks may limit the inference of the underlying dysfunctional connectivity, which is hypothesized to be a factor in patient symptoms, as fMRI measures connectivity via hemodynamics. Therefore, combination of connectivity assessments using fMRI and magnetoencephalography (MEG), which more directly measures neuronal activity, may provide improved classification of schizophrenia than either modality alone. Moreover, recent evidence indicates that metrics of dynamic connectivity may also be critical for understanding pathology in schizophrenia. In this work, we propose a new framework for extraction of important disease related features and classification of patients with schizophrenia based on using both fMRI and MEG to investigate functional network components in the resting state. Results of this study show that the integration of fMRI and MEG provides important information that captures fundamental characteristics of functional network connectivity in schizophrenia and is helpful for prediction of schizophrenia patient group membership. Combined fMRI/MEG methods, using static functional network connectivity analyses, improved classification accuracy relative to use of fMRI or MEG methods alone (by 15 and 12.45%, respectively), while combined fMRI/MEG methods using dynamic functional network connectivity analyses improved classification up to 5.12% relative to use of fMRI alone and up to 17.21% relative to use of MEG alone.

Age-related differences in practice-dependent resting-state functional connectivity related to motor sequence learning

Decreased neural plasticity is observed with healthy ageing in the primary sensorimotor (SM1) cortex thought to participate in motor learning and memory consolidation processes. In the present magnetoencephalography study, the post-training reorganization of resting-state functional connectivity (rsFC) and its relation with motor learning and early consolidation in 14 young (19-30 years) and 14 old (66-70 years) healthy participants were investigated. At the behavioral level, participants were trained on a motor sequence learning task then retested 20-30 min later for transient offline gains in performance. Using a sensorimotor seed-based approach, rsFC relying on beta band power envelope correlation was estimated immediately before and 10 min after the learning episode. Post-training changes in rsFC (from before to after learning) were correlated with motor learning performance and with the offline improvement in performance within the hour after learning. Young and old participants exhibited differential patterns of sensorimotor-related rsFC, bearing specific relationships with motor learning and consolidation. Our findings suggest that rsFC changes following learning reflect the offline processing of the new motor skill and contribute to the early memory consolidation within the hour after learning. Furthermore, differences in post-training changes in rsFC between young and old participants support the hypothesis that ageing modulates the neural circuits underlying the learning of a new motor skill and the early subsequent consolidation stages. Hum Brain Mapp 38:923-937, 2017. © 2016 Wiley Periodicals, Inc.

Magnetoencephalographic and functional MRI connectomics in schizophrenia via intra- and inter-network connectivity

Examination of intrinsic functional connectivity using functional MRI (fMRI) has provided important findings regarding dysconnectivity in schizophrenia. Extending these results using a complementary neuroimaging modality, magnetoencephalography (MEG), we present the first direct comparison of functional connectivity between schizophrenia patients and controls, using these two modalities combined. We developed a novel MEG approach for estimation of networks using MEG that incorporates spatial independent component analysis (ICA) and pairwise correlations between independent component timecourses, to estimate intra- and intern-network connectivity. This analysis enables group-level inference and testing of between-group differences. Resting state MEG and fMRI data were acquired from a large sample of healthy controls (n=45) and schizophrenia patients (n=46). Group spatial ICA was performed on fMRI and MEG data to extract intrinsic fMRI and MEG networks and to compensate for signal leakage in MEG. Similar, but not identical spatial independent components were detected for MEG and fMRI. Analysis of functional network connectivity (FNC; i.e., pairwise correlations in network (ICA component) timecourses) revealed a differential between-modalities pattern, with greater connectivity among occipital networks in fMRI and among frontal networks in MEG. Most importantly, significant differences between controls and patients were observed in both modalities. MEG FNC results in particular indicated dysfunctional hyperconnectivity within frontal and temporal networks in patients, while in fMRI FNC was always greater for controls than for patients. This is the first study to apply group spatial ICA as an approach to leakage correction, and as such our results may be biased by spatial leakage effects. Results suggest that combining these two neuroimaging modalities reveals additional disease-relevant patterns of connectivity that were not detectable with fMRI or MEG alone.

Quantifying the Test-Retest Reliability of Magnetoencephalography Resting-State Functional Connectivity

The coordinated activity of the resting-state brain can be evaluated with magnetoencephalography (MEG) for distinct brain rhythms by performing source reconstruction to estimate the activities of target brain regions and employing one of the many existent functional connectivity (FC) algorithms. Although this procedure has been applied in a great amount of studies both with healthy and pathological populations, the reliability of such FC estimates is unknown, and this impairs the use of resting-state MEG FC at the individual level. In this study, the test-retest reliability of MEG resting FC was evaluated by exploring both within- and between-subject variability in FC in 16 healthy subjects who underwent three resting-state MEG scans. FC was computed after beamforming source reconstruction with four popular FC metrics: phase-locking value (PLV), phase lag index (PLI), direct envelope correlation (d-ecor), and envelope correlation with leakage correction (lc-ecor). Then, test-restest reliability and within- and between-subject agreement were evaluated with the intraclass correlation coefficient (ICC) and Kendall's W, respectively. Reliability was found to depend on the FC metric, the frequency band, and the specific link. As a general trend, greater test-retest reliability was found for PLV in theta to gamma, and for lc-ecor and d-ecor in beta. Further inspection of the ICC distribution revealed that volume conduction effects could be contributing to high ICC in PLV and d-ecor. In addition, stronger links were found to be more reliable. Overall, this encourages the further use of resting-state MEG FC for individual-level studies, especially with PLV or envelope correlation metrics.

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

Understanding the electrophysiological basis of resting state networks (RSNs) in the human brain is a critical step towards elucidating how inter-areal connectivity supports healthy brain function. In recent years, the relationship between RSNs (typically measured using haemodynamic signals) and electrophysiology has been explored using functional Magnetic Resonance Imaging (fMRI) and magnetoencephalography (MEG). Significant progress has been made, with similar spatial structure observable in both modalities. However, there is a pressing need to understand this relationship beyond simple visual similarity of RSN patterns. Here, we introduce a mathematical model to predict fMRI-based RSNs using MEG. Our unique model, based upon a multivariate Taylor series, incorporates both phase and amplitude based MEG connectivity metrics, as well as linear and non-linear interactions within and between neural oscillations measured in multiple frequency bands. We show that including non-linear interactions, multiple frequency bands and cross-frequency terms significantly improves fMRI network prediction. This shows that fMRI connectivity is not only the result of direct electrophysiological connections, but is also driven by the overlap of connectivity profiles between separate regions. Our results indicate that a complete understanding of the electrophysiological basis of RSNs goes beyond simple frequency-specific analysis, and further exploration of non-linear and cross-frequency interactions will shed new light on distributed network connectivity, and its perturbation in pathology.

•

We introduce a mathematical model to predict fMRI-based RSNs using MEG.

•

Our model is based on a multi-variate Taylor series expansion.

•

The electrophysiological basis of RSNs goes beyond frequency-band specific analysis.

•

RSNs result 1) from multiple frequency bands and cross-frequency coupling.

•

RSNs result 2) from direct and shared electrophysiological connectivity.

Multi-Dimensional Dynamics of Human Electromagnetic Brain Activity

Magnetoencephalography (MEG) and electroencephalography (EEG) are invaluable neuroscientific tools for unveiling human neural dynamics in three dimensions (space, time, and frequency), which are associated with a wide variety of perceptions, cognition, and actions. MEG/EEG also provides different categories of neuronal indices including activity magnitude, connectivity, and network properties along the three dimensions. In the last 20 years, interest has increased in inter-regional connectivity and complex network properties assessed by various sophisticated scientific analyses. We herein review the definition, computation, short history, and pros and cons of connectivity and complex network (graph-theory) analyses applied to MEG/EEG signals. We briefly describe recent developments in source reconstruction algorithms essential for source-space connectivity and network analyses. Furthermore, we discuss a relatively novel approach used in MEG/EEG studies to examine the complex dynamics represented by human brain activity. The correct and effective use of these neuronal metrics provides a new insight into the multi-dimensional dynamics of the neural representations of various functions in the complex human brain.

Evidence that Subanesthetic Doses of Ketamine Cause Sustained Disruptions of NMDA and AMPA-Mediated Frontoparietal Connectivity in Humans

Following the discovery of the antidepressant properties of ketamine, there has been a recent resurgence in the interest in this NMDA receptor antagonist. Although detailed animal models of the molecular mechanisms underlying ketamine's effects have emerged, there are few MEG/EEG studies examining the acute subanesthetic effects of ketamine infusion in man. We recorded 275 channel MEG in two experiments (n = 25 human males) examining the effects of subanesthetic ketamine infusion. MEG power spectra revealed a rich set of significant oscillatory changes compared with placebo sessions, including decreases in occipital, parietal, and anterior cingulate alpha power, increases in medial frontal theta power, and increases in parietal and cingulate cortex high gamma power. Each of these spectral effects demonstrated their own set of temporal dynamics. Dynamic causal modeling of frontoparietal connectivity changes with ketamine indicated a decrease in NMDA and AMPA-mediated frontal-to-parietal connectivity. AMPA-mediated connectivity changes were sustained for up to 50 min after ketamine infusion had ceased, by which time perceptual distortions were absent. The results also indicated a decrease in gain of parietal pyramidal cells, which was correlated with participants' self-reports of blissful state. Based on these results, we suggest that the antidepressant effects of ketamine may depend on its ability to change the balance of frontoparietal connectivity patterns.

SIGNIFICANCE STATEMENT

In this paper, we found that subanesthetic doses of ketamine, similar to those used in antidepressant studies, increase anterior theta and gamma power but decrease posterior theta, delta, and alpha power, as revealed by magnetoencephalographic recordings. Dynamic causal modeling of frontoparietal connectivity changes with ketamine indicated a decrease in NMDA and AMPA-mediated frontal-to-parietal connectivity. AMPA-mediated connectivity changes were sustained for up to 50 min after ketamine infusion had ceased, by which time perceptual distortions were absent. The results also indicated a decrease in gain of parietal pyramidal cells, which was correlated with participants' self-reports of blissful state. The alterations in frontoparietal connectivity patterns we observe here may be important in generating the antidepressant response to ketamine.

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

The human brain can be divided into multiple areas, each responsible for different aspects of behaviour. Healthy brain function relies upon efficient connectivity between these areas and, in recent years, neuroimaging has been revolutionised by an ability to estimate this connectivity. In this paper we discuss measurement of network connectivity using magnetoencephalography (MEG), a technique capable of imaging electrophysiological brain activity with good (~5 mm) spatial resolution and excellent (~1 ms) temporal resolution. The rich information content of MEG facilitates many disparate measures of connectivity between spatially separate regions and in this paper we discuss a single metric known as power envelope correlation. We review in detail the methodology required to measure power envelope correlation including (i) projection of MEG data into source space, (ii) removing confounds introduced by the MEG inverse problem and (iii) estimation of connectivity itself. In this way, we aim to provide researchers with a description of the key steps required to assess envelope based functional networks, which are thought to represent an intrinsic mode of coupling in the human brain. We highlight the principal findings of the techniques discussed, and furthermore, we show evidence that this method can probe how the brain forms and dissolves multiple transient networks on a rapid timescale in order to support current processing demand. Overall, power envelope correlation offers a unique and verifiable means to gain novel insights into network coordination and is proving to be of significant value in elucidating the neural dynamics of the human connectome in health and disease.

A geometric correction scheme for spatial leakage effects in MEG/EEG seed-based functional connectivity mapping

Spatial leakage effects are particularly confounding for seed-based investigations of brain networks using source-level electroencephalography (EEG) or magnetoencephalography (MEG). Various methods designed to avoid this issue have been introduced but are limited to particular assumptions about its temporal characteristics. Here, we investigate the usefulness of a model-based geometric correction scheme (GCS) to suppress spatial leakage emanating from the seed location. We analyze its properties theoretically and then assess potential advantages and limitations with simulated and experimental MEG data (resting state and auditory-motor task). To do so, we apply Minimum Norm Estimation (MNE) for source reconstruction and use variation of error parameters, statistical gauging of spatial leakage correction and comparison with signal orthogonalization. Results show that the GCS has a local (i.e., near the seed) effect only, in line with the geometry of MNE spatial leakage, and is able to map spatially all types of brain interactions, including linear correlations eliminated after signal orthogonalization. Furthermore, it is robust against the introduction of forward model errors. On the other hand, the GCS can be affected by local overcorrection effects and seed mislocation. These issues arise with signal orthogonalization too, although significantly less extensively, so the two approaches complement each other. The GCS thus appears to be a valuable addition to the spatial leakage correction toolkits for seed-based FC analyses in source-projected MEG/EEG data.

Resting-state oscillatory dynamics in sensorimotor cortex in benign epilepsy with centro-temporal spikes and typical brain development

Benign Epilepsy with Centro-Temporal Spikes (BECTS) is a common childhood epilepsy associated with deficits in several neurocognitive domains. Neurophysiological studies in BECTS often focus on centro-temporal spikes, but these correlate poorly with morphology and cognitive impairments. To better understand the neural profile of BECTS, we studied background brain oscillations, thought to be integrally involved in neural network communication, in sensorimotor areas. We used independent component analysis of temporally correlated sources on magnetoencephalography recordings to assess sensorimotor resting-state network activity in BECTS patients and typically developing controls. We also investigated the variability of oscillatory characteristics within focal primary motor cortex (M1), localized with a separate finger abduction task. We hypothesized that background oscillations would differ between patients and controls in the sensorimotor network but not elsewhere, especially in the beta band (13-30 Hz) because of its role in network communication and motor processing. The results support our hypothesis: in the sensorimotor network, patients had a greater variability in oscillatory amplitude compared to controls, whereas there was no difference in the visual network. Network measures did not correlate with age. The coefficient of variation of resting M1 peak frequency correlated negatively with age in the beta band only, and was greater than average for a number of patients. Our results point toward a "disorganized" functional sensorimotor network in BECTS, supporting a neurodevelopmental delay in sensorimotor cortex. Our findings further suggest that investigating the variability of oscillatory peak frequency may be a useful tool to investigate deficits of disorganization in neurodevelopmental disorders.

Group differences in MEG-ICA derived resting state networks: Application to major depressive disorder

Functional magnetic resonance imaging (fMRI) studies have revealed the existence of robust, interconnected brain networks exhibiting correlated low frequency fluctuations during rest, which can be derived by examining inherent spatio-temporal patterns in functional scans independent of any a priori model. In order to explore the electrophysiological underpinnings of these networks, analogous techniques have recently been applied to magnetoencephalography (MEG) data, revealing similar networks that exhibit correlated low frequency fluctuations in the power envelope of beta band (14-30Hz) power. However, studies to date using this technique have concentrated on healthy subjects, and no method has yet been presented for group comparisons. We extended the ICA resting state MEG method to enable group comparisons, and demonstrate the technique in a sample of subjects with major depressive disorder (MDD). We found that the intrinsic resting state networks evident in fMRI appeared to be disrupted in individuals with MDD compared to healthy participants, particularly in the subgenual cingulate, although the electrophysiological correlates of this are unknown. Networks extracted from a combined group of healthy and MDD participants were examined for differences between groups. Individuals with MDD showed reduced correlations between the subgenual anterior cingulate (sgACC) and hippocampus in a network with primary nodes in the precentral and middle frontal gyri. Individuals with MDD also showed increased correlations between insulo-temporal nodes and amygdala compared to healthy controls. To further support our methods and findings, we present test/re-test reliability on independent recordings acquired within the same session. Our results demonstrate that group analyses are possible with the resting state MEG-independent component analysis (ICA) technique, highlighting a new pathway for analysis and discovery. This study also provides the first evidence of altered sgACC connectivity with a motor network. This finding, reliable across multiple sessions, suggests that the sgACC may partially mediate the psychomotor symptoms of MDD via synchronized changes in beta-band power, and expands the idea of the sgACC as a hub region mediating cognitive and emotional symptomatic domains in MDD. Findings of increased connectivity between the amygdala and cortical nodes further support the role of amygdalar networks in mediated depressive symptomatology.

CLINICAL TRIALS IDENTIFIER

NCT00024635 (ZIA-MH002927-04).

The effect of ageing on fMRI: Correction for the confounding effects of vascular reactivity evaluated by joint fMRI and MEG in 335 adults

In functional magnetic resonance imaging (fMRI) research one is typically interested in neural activity. However, the blood‐oxygenation level‐dependent (BOLD) signal is a composite of both neural and vascular activity. As factors such as age or medication may alter vascular function, it is essential to account for changes in neurovascular coupling when investigating neurocognitive functioning with fMRI. The resting‐state fluctuation amplitude (RSFA) in the fMRI signal (rsfMRI) has been proposed as an index of vascular reactivity. The RSFA compares favourably with other techniques such as breath‐hold and hypercapnia, but the latter are more difficult to perform in some populations, such as older adults. The RSFA is therefore a candidate for use in adjusting for age‐related changes in vascular reactivity in fMRI studies. The use of RSFA is predicated on its sensitivity to vascular rather than neural factors; however, the extent to which each of these factors contributes to RSFA remains to be characterized. The present work addressed these issues by comparing RSFA (i.e., rsfMRI variability) to proxy measures of (i) cardiovascular function in terms of heart rate (HR) and heart rate variability (HRV) and (ii) neural activity in terms of resting state magnetoencephalography (rsMEG). We derived summary scores of RSFA, a sensorimotor task BOLD activation, cardiovascular function and rsMEG variability for 335 healthy older adults in the population‐based Cambridge Centre for Ageing and Neuroscience cohort (Cam‐CAN; http://www.cam-can.com). Mediation analysis revealed that the effects of ageing on RSFA were significantly mediated by vascular factors, but importantly not by the variability in neuronal activity. Furthermore, the converse effects of ageing on the rsMEG variability were not mediated by vascular factors. We then examined the effect of RSFA scaling of task‐based BOLD in the sensorimotor task. The scaling analysis revealed that much of the effects of age on task‐based activation studies with fMRI do not survive correction for changes in vascular reactivity, and are likely to have been overestimated in previous fMRI studies of ageing. The results from the mediation analysis demonstrate that RSFA is modulated by measures of vascular function and is not driven solely by changes in the variance of neural activity. Based on these findings we propose that the RSFA scaling method is articularly useful in large scale and longitudinal neuroimaging studies of ageing, or with frail participants, where alternative measures of vascular reactivity are impractical. Hum Brain Mapp 36:2248–2269, 2015. © 2015 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

Multimodal based classification of schizophrenia patients

Schizophrenia is currently diagnosed by physicians through clinical assessment and their evaluation of patient's self-reported experiences over the longitudinal course of the illness. There is great interest in identifying biologically based markers at the onset of illness, rather than relying on the evolution of symptoms across time. Functional network connectivity shows promise in providing individual subject predictive power. The majority of previous studies considered the analysis of functional connectivity during resting-state using only fMRI. However, exclusive reliance on fMRI to generate such networks, may limit inference on dysfunctional connectivity, which is hypothesized to underlie patient symptoms. In this work, we propose a framework for classification of schizophrenia patients and healthy control subjects based on using both fMRI and band limited envelope correlation metrics in MEG to interrogate functional network components in the resting state. Our results show that the combination of these two methods provide valuable information that captures fundamental characteristics of brain network connectivity in schizophrenia. Such information is useful for prediction of schizophrenia patients. Classification accuracy performance was improved significantly (up to ≈ 7%) relative to only the fMRI method and (up to ≈ 21%) relative to only the MEG method.

Power spectral aspects of the default mode network in schizophrenia: an MEG study

Background

Symptoms of schizophrenia are related to deficits in self-monitoring function, which may be a consequence of irregularity in aspects of the default mode network (DMN). Schizophrenia can also be characterized by a functional abnormality of the brain activity that is reflected in the resting state. Oscillatory analysis provides an important understanding of resting brain activity. However, conventional methods using electroencephalography are restricted because of low spatial resolution, despite their excellent temporal resolution.

The aim of this study was to investigate resting brain oscillation and the default mode network based on a source space in various frequency bands such as theta, alpha, beta, and gamma using magnetoencephalography. In addition, we investigated whether these resting and DMN activities could distinguish schizophrenia patients from normal controls. To do this, the power spectral density of each frequency band at rest was imaged and compared on a spatially normalized brain template in 20 patients and 20 controls.

Results

The spatial distribution of DMN activity in the alpha band was similar to that found in previous fMRI studies. The posterior cingulate cortex (PCC) and lateral inferior parietal cortex were activated at rest, while the medial prefrontal cortex (MPFC) was deactivated at rest rather than during the task. Although the MPFC and PCC regions exhibited contrasting activation patterns, these two regions were significantly coherent at rest. The DMN and resting activities of the PCC were increased in schizophrenia patients, predominantly in the theta and alpha bands.

Conclusions

By using MEG to identify the DMN regions, predominantly in the alpha band, we found that both resting and DMN activities were augmented in the posterior cingulate in schizophrenia patients. Furthermore, schizophrenia patients exhibited decreased coherence between the PCC and MPFC in the gamma band at rest.

Guiding functional connectivity estimation by structural connectivity in MEG: an application to discrimination of conditions of mild cognitive impairment

Whole brain resting state connectivity is a promising biomarker that might help to obtain an early diagnosis in many neurological diseases, such as dementia. Inferring resting-state connectivity is often based on correlations, which are sensitive to indirect connections, leading to an inaccurate representation of the real backbone of the network. The precision matrix is a better representation for whole brain connectivity, as it considers only direct connections. The network structure can be estimated using the graphical lasso (GL), which achieves sparsity through l₁-regularization on the precision matrix. In this paper, we propose a structural connectivity adaptive version of the GL, where weaker anatomical connections are represented as stronger penalties on the corresponding functional connections. We applied beamformer source reconstruction to the resting state MEG recordings of 81 subjects, where 29 were healthy controls, 22 were single-domain amnestic Mild Cognitive Impaired (MCI), and 30 were multiple-domain amnestic MCI. An atlas-based anatomical parcellation of 66 regions was obtained for each subject, and time series were assigned to each of the regions. The fiber densities between the regions, obtained with deterministic tractography from diffusion-weighted MRI, were used to define the anatomical connectivity. Precision matrices were obtained with the region specific time series in five different frequency bands. We compared our method with the traditional GL and a functional adaptive version of the GL, in terms of log-likelihood and classification accuracies between the three groups. We conclude that introducing an anatomical prior improves the expressivity of the model and, in most cases, leads to a better classification between groups.

•

We propose an anatomy-driven method for functional connectivity estimation in MEG.

•

Structural prior contributes to a better representation of the functional connectivity.

•

The proposed method is shown to be useful as a biomarker for classification of MCI.

Graph theoretical analysis of resting-state MEG data: Identifying interhemispheric connectivity and the default mode

Interhemispheric connectivity with resting state MEG has been elusive, and demonstration of the default mode network (DMN) yet more challenging. Recent seed-based MEG analyses have shown interhemispheric connectivity using power envelope correlations. The purpose of this study is to compare graph theoretic maps of brain connectivity generated using MEG with and without signal leakage correction to evaluate for the presence of interhemispheric connectivity. Eight minutes of resting state eyes-open MEG data were obtained in 22 normal male subjects enrolled in an IRB-approved study (ages 16-18). Data were processed using an in-house automated MEG processing pipeline and projected into standard (MNI) source space at 7mm resolution using a scalar beamformer. Mean beta-band amplitude was sampled at 2.5second epochs from the source space time series. Leakage correction was performed in the time domain of the source space beam formed signal prior to amplitude transformation. Graph theoretic voxel-wise source space correlation connectivity analysis was performed for leakage corrected and uncorrected data. Degree maps were thresholded across subjects for the top 20% of connected nodes to identify hubs. Additional degree maps for sensory, visual, motor, and temporal regions were generated to identify interhemispheric connectivity using laterality indices. Hubs for the uncorrected MEG networks were predominantly symmetric and midline, bearing some resemblance to fMRI networks. These included the cingulate cortex, bilateral inferior frontal lobes, bilateral hippocampal formations and bilateral cerebellar hemispheres. These uncorrected networks however, demonstrated little to no interhemispheric connectivity using the ROI-based degree maps. Leakage corrected MEG data identified the DMN, with hubs in the posterior cingulate and biparietal areas. These corrected networks demonstrated robust interhemispheric connectivity for the ROI-based degree maps. Graph theoretic analysis of MEG resting state data without signal leakage correction can demonstrate symmetric networks with some resemblance to fMRI networks. These networks however, are an artifact of high local correlation from signal leakage and lack interhemispheric connectivity. Following signal leakage correction, MEG hubs emerge in the DMN, with strong interhemispheric connectivity.

Multi-session statistics on beamformed MEG data

Beamforming has been widely adopted as a source reconstruction technique in the analysis of magnetoencephalography data. Most beamforming implementations incorporate a spatially-varying rescaling (which we term weights normalisation) to correct for the inherent depth bias in raw beamformer estimates. Here, we demonstrate that such rescaling can cause critical problems whenever analyses are performed over multiple sessions of separately beamformed data, for example when comparing effect sizes between different populations. Importantly, we show that the weights-normalised beamformer estimates of neural activity can even lead to a reversal in the inferred sign of the effect being measured. We instead recommend that no weights normalisation be carried out; any depth bias is instead accounted for in the calculation of multi-session (e.g. group) statistics. We demonstrate the severity of the weights normalisation confound with a 2-D simulation, and in real MEG data by performing a group statistical analysis to detect differences in alpha power in eyes-closed rest compared with continuous visual stimulation.

Broadband cortical desynchronization underlies the human psychedelic state

Psychedelic drugs produce profound changes in consciousness, but the underlying neurobiological mechanisms for this remain unclear. Spontaneous and induced oscillatory activity was recorded in healthy human participants with magnetoencephalography after intravenous infusion of psilocybin--prodrug of the nonselective serotonin 2A receptor agonist and classic psychedelic psilocin. Psilocybin reduced spontaneous cortical oscillatory power from 1 to 50 Hz in posterior association cortices, and from 8 to 100 Hz in frontal association cortices. Large decreases in oscillatory power were seen in areas of the default-mode network. Independent component analysis was used to identify a number of resting-state networks, and activity in these was similarly decreased after psilocybin. Psilocybin had no effect on low-level visually induced and motor-induced gamma-band oscillations, suggesting that some basic elements of oscillatory brain activity are relatively preserved during the psychedelic experience. Dynamic causal modeling revealed that posterior cingulate cortex desynchronization can be explained by increased excitability of deep-layer pyramidal neurons, which are known to be rich in 5-HT2A receptors. These findings suggest that the subjective effects of psychedelics result from a desynchronization of ongoing oscillatory rhythms in the cortex, likely triggered by 5-HT2A receptor-mediated excitation of deep pyramidal cells.

The contribution of electrophysiology to functional connectivity mapping

A powerful way to probe brain function is to assess the relationship between simultaneous changes in activity across different parts of the brain. In recent years, the temporal activity correlation between brain areas has frequently been taken as a measure of their functional connections. Evaluating 'functional connectivity' in this way is particularly popular in the fMRI community, but has also drawn interest among electrophysiologists. Like hemodynamic fluctuations observed with fMRI, electrophysiological signals display significant temporal fluctuations, even in the absence of a stimulus. These neural fluctuations exhibit a correlational structure over a wide range of spatial and temporal scales. Initial evidence suggests that certain aspects of this correlational structure bear a high correspondence to so-called functional networks defined using fMRI. The growing family of methods to study activity covariation, combined with the diverse neural mechanisms that contribute to the spontaneous fluctuations, has somewhat blurred the operational concept of functional connectivity. What is clear is that spontaneous activity is a conspicuous, energy-consuming feature of the brain. Given its prominence and its practical applications for the functional connectivity mapping of brain networks, it is of increasing importance that we understand its neural origins as well as its contribution to normal brain function.

MEG source imaging method using fast L1 minimum-norm and its applications to signals with brain noise and human resting-state source amplitude images

The present study developed a fast MEG source imaging technique based on Fast Vector-based Spatio-Temporal Analysis using a L1-minimum-norm (Fast-VESTAL) and then used the method to obtain the source amplitude images of resting-state magnetoencephalography (MEG) signals for different frequency bands. The Fast-VESTAL technique consists of two steps. First, L1-minimum-norm MEG source images were obtained for the dominant spatial modes of sensor-waveform covariance matrix. Next, accurate source time-courses with millisecond temporal resolution were obtained using an inverse operator constructed from the spatial source images of Step 1. Using simulations, Fast-VESTAL's performance was assessed for its 1) ability to localize multiple correlated sources; 2) ability to faithfully recover source time-courses; 3) robustness to different SNR conditions including SNR with negative dB levels; 4) capability to handle correlated brain noise; and 5) statistical maps of MEG source images. An objective pre-whitening method was also developed and integrated with Fast-VESTAL to remove correlated brain noise. Fast-VESTAL's performance was then examined in the analysis of human median-nerve MEG responses. The results demonstrated that this method easily distinguished sources in the entire somatosensory network. Next, Fast-VESTAL was applied to obtain the first whole-head MEG source-amplitude images from resting-state signals in 41 healthy control subjects, for all standard frequency bands. Comparisons between resting-state MEG sources images and known neurophysiology were provided. Additionally, in simulations and cases with MEG human responses, the results obtained from using conventional beamformer technique were compared with those from Fast-VESTAL, which highlighted the beamformer's problems of signal leaking and distorted source time-courses.