Quantifying the performance of MEG source reconstruction using resting state data.

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

Loss of brain inter-frequency hubs in Alzheimer's disease

Alzheimer's disease (AD) causes alterations of brain network structure and function. The latter consists of connectivity changes between oscillatory processes at different frequency channels. We proposed a multi-layer network approach to analyze multiple-frequency brain networks inferred from magnetoencephalographic recordings during resting-states in AD subjects and age-matched controls. Main results showed that brain networks tend to facilitate information propagation across different frequencies, as measured by the multi-participation coefficient (MPC). However, regional connectivity in AD subjects was abnormally distributed across frequency bands as compared to controls, causing significant decreases of MPC. This effect was mainly localized in association areas and in the cingulate cortex, which acted, in the healthy group, as a true inter-frequency hub. MPC values significantly correlated with memory impairment of AD subjects, as measured by the total recall score. Most predictive regions belonged to components of the default-mode network that are typically affected by atrophy, metabolism disruption and amyloid-beta deposition. We evaluated the diagnostic power of the MPC and we showed that it led to increased classification accuracy (78.39%) and sensitivity (91.11%). These findings shed new light on the brain functional alterations underlying AD and provide analytical tools for identifying multi-frequency neural mechanisms of brain diseases.Comment: 27 pages, 6 figures, 3 tables, 3 supplementary figure

Dynamic Decomposition of Spatiotemporal Neural Signals

Neural signals are characterized by rich temporal and spatiotemporal dynamics that reflect the organization of cortical networks. Theoretical research has shown how neural networks can operate at different dynamic ranges that correspond to specific types of information processing. Here we present a data analysis framework that uses a linearized model of these dynamic states in order to decompose the measured neural signal into a series of components that capture both rhythmic and non-rhythmic neural activity. The method is based on stochastic differential equations and Gaussian process regression. Through computer simulations and analysis of magnetoencephalographic data, we demonstrate the efficacy of the method in identifying meaningful modulations of oscillatory signals corrupted by structured temporal and spatiotemporal noise. These results suggest that the method is particularly suitable for the analysis and interpretation of complex temporal and spatiotemporal neural signals

Whole brain functional connectivity using phase locking measures of resting state magnetoencephalography

The analysis of spontaneous functional connectivity (sFC) reveals the statistical connections between regions of the brain consistent with underlying functional communication networks within the brain. In this work, we describe the implementation of a complete all-to-all network analysis of resting state neuronal activity from magnetoencephalography (MEG). Using graph theory to define networks at the dipole level, we established functionally defined regions by k-means clustering cortical surface locations using Eigenvector centrality (EVC) scores from the all-to-all adjacency model. Permutation testing was used to estimate regions with statistically significant connections compared to empty room data, which adjusts for spatial dependencies introduced by the MEG inverse problem. In order to test this model, we performed a series of numerical simulations investigating the effects of the MEG reconstruction on connectivity estimates. We subsequently applied the approach to subject data to investigate the effectiveness of our method in obtaining whole brain networks. Our findings indicated that our model provides statistically robust estimates of functional region networks. Application of our phase locking network methodology to real data produced networks with similar connectivity to previously published findings, specifically, we found connections between contralateral areas of the arcuate fasciculus that have been previously investigated. The use of data-driven methods for neuroscientific investigations provides a new tool for researchers in identifying and characterizing whole brain functional connectivity networks. © 2014 Schmidt, Ghuman and Huppert

MEG cortical microstates: Spatiotemporal characteristics, dynamic functional connectivity and stimulus-evoked responses

EEG microstate analysis is an approach to study brain states and their fast transitions in healthy cognition and disease. A key limitation of conventional microstate analysis is that it must be performed at the sensor level, and therefore gives limited anatomical insight. Here, we generalise the microstate methodology to be applicable to source-reconstructed electrophysiological data. Using simulations of a neural-mass network model, we first established the validity and robustness of the proposed method. Using MEG resting-state data, we uncovered ten microstates with distinct spatial distributions of cortical activation. Multivariate pattern analysis demonstrated that source-level microstates were associated with distinct functional connectivity patterns. We further demonstrated that the occurrence probability of MEG microstates were altered by auditory stimuli, exhibiting a hyperactivity of the microstate including the auditory cortex. Our results support the use of source-level microstates as a method for investigating brain dynamic activity and connectivity at the millisecond scale

Efficient transfer entropy analysis of non-stationary neural time series

Information theory allows us to investigate information processing in neural systems in terms of information transfer, storage and modification. Especially the measure of information transfer, transfer entropy, has seen a dramatic surge of interest in neuroscience. Estimating transfer entropy from two processes requires the observation of multiple realizations of these processes to estimate associated probability density functions. To obtain these observations, available estimators assume stationarity of processes to allow pooling of observations over time. This assumption however, is a major obstacle to the application of these estimators in neuroscience as observed processes are often non-stationary. As a solution, Gomez-Herrero and colleagues theoretically showed that the stationarity assumption may be avoided by estimating transfer entropy from an ensemble of realizations. Such an ensemble is often readily available in neuroscience experiments in the form of experimental trials. Thus, in this work we combine the ensemble method with a recently proposed transfer entropy estimator to make transfer entropy estimation applicable to non-stationary time series. We present an efficient implementation of the approach that deals with the increased computational demand of the ensemble method's practical application. In particular, we use a massively parallel implementation for a graphics processing unit to handle the computationally most heavy aspects of the ensemble method. We test the performance and robustness of our implementation on data from simulated stochastic processes and demonstrate the method's applicability to magnetoencephalographic data. While we mainly evaluate the proposed method for neuroscientific data, we expect it to be applicable in a variety of fields that are concerned with the analysis of information transfer in complex biological, social, and artificial systems.Comment: 27 pages, 7 figures, submitted to PLOS ON

Long-Range Temporal Correlations in Resting State Beta Oscillations are Reduced in Schizophrenia

Symptoms of schizophrenia (SCZ) are likely to be generated by genetically mediated synaptic dysfunction, which contribute to large-scale functional neural dysconnectivity. Recent electrophysiological studies suggest that this dysconnectivity is present not only at a spatial level but also at a temporal level, operationalized as long-range temporal correlations (LRTCs). Previous research suggests that alpha and beta frequency bands have weaker temporal stability in people with SCZ. This study sought to replicate these findings with high-density electroencephalography (EEG), enabling a spatially more accurate analysis of LRTC differences, and to test associations with characteristic SCZ symptoms and cognitive deficits. A 128-channel EEG was used to record eyes-open resting state brain activity of 23 people with SCZ and 24 matched healthy controls (HCs). LRTCs were derived for alpha (8–12 Hz) and beta (13–25 Hz) frequency bands. As an exploratory analysis, LRTC was source projected using sLoreta. People with SCZ showed an area of significantly reduced beta-band LRTC compared with HCs over bilateral posterior regions. There were no between-group differences in alpha-band activity. Individual symptoms of SCZ were not related to LRTC values nor were cognitive deficits. The study confirms that people with SCZ have reduced temporal stability in the beta frequency band. The absence of group differences in the alpha band may be attributed to the fact that people had, in contrast to previous studies, their eyes open in the current study. Taken together, our study confirms the utility of LRTC as a marker of network instability in people with SCZ and provides a novel empirical perspective for future examinations of network dysfunction salience in SCZ research