Functional brain alterations in symptomatic dermographism patients-An exploratory magnetoencephalography study.

Symptomatic dermographism (SD) is a common form of urticaria, which is triggered by stroking the skin. Brain involvement in its aetiology was investigated by means of magnetoencephalography (MEG) after provocation with histamine and dermography. Wheals were induced by histamine skin prick test and dermography in twelve SD patients and fourteen controls. Itch severity was scored on a Visual Analogue Scale (VAS). Relative power and functional connectivity (FC) were measured using a 306-channel whole-head MEG system at baseline and 10 min after histamine and dermography, and contrasted between groups and conditions. Furthermore, wheal diameter and itch scores after these procedures were correlated with the MEG values. SD patients had higher itch scores after histamine and dermography. No significant group-differences were observed in relative power or FC for any condition. In both groups, power decreases were mostly observed in the beta band, and power increases in the alpha bands, after provocation, with more regions involved in patients compared to controls. Increased FC was seen after histamine in patients, and after dermography in controls. In patients only, dermography and histamine wheal size correlated with the alpha2 power in the regions of interest that showed significant condition effects after these procedures. Our findings may be cautiously interpreted as aberrant itch processing, and suggest involvement of the central nervous system in the aetiology of SD.

The relationship between pathological brain activity and functional network connectivity in glioma patients.

PURPOSE: Glioma is associated with pathologically high (peri)tumoral brain activity, which relates to faster progression. Functional connectivity is disturbed locally and throughout the entire brain, associating with symptomatology. We, therefore, investigated how local activity and network measures relate to better understand how the intricate relationship between the tumor and the rest of the brain may impact disease and symptom progression. METHODS: We obtained magnetoencephalography in 84 de novo glioma patients and 61 matched healthy controls. The offset of the power spectrum, a proxy of neuronal activity, was calculated for 210 cortical regions. We calculated patients' regional deviations in delta, theta and lower alpha network connectivity as compared to controls, using two network measures: clustering coefficient (local connectivity) and eigenvector centrality (integrative connectivity). We then tested group differences in activity and connectivity between (peri)tumoral, contralateral homologue regions, and the rest of the brain. We also correlated regional offset to connectivity. RESULTS: As expected, patients' (peri)tumoral activity was pathologically high, and patients showed higher clustering and lower centrality than controls. At the group-level, regionally high activity related to high clustering in controls and patients alike. However, within-patient analyses revealed negative associations between regional deviations in brain activity and clustering, such that pathologically high activity coincided with low network clustering, while regions with 'normal' activity levels showed high network clustering. CONCLUSION: Our results indicate that pathological activity and connectivity co-localize in a complex manner in glioma. This insight is relevant to our understanding of disease progression and cognitive symptomatology.

Tau protein spreads through functionally connected neurons in Alzheimer's disease: a combined MEG/PET study.

Recent studies on Alzheimer's disease (AD) suggest that tau proteins spread through the brain following neuronal connections. Several mechanisms could be involved in this process: spreading between brain regions that interact strongly (functional connectivity); through the pattern of anatomical connections (structural connectivity); or simple diffusion. Using magnetoencephalography (MEG), we investigated which spreading pathways influence tau protein spreading by modelling the tau propagation process using an epidemic spreading model. We compared the modelled tau depositions with 18F-flortaucipir PET binding potentials at several stages of the AD continuum. In this cross-sectional study, we analysed source-reconstructed MEG data and dynamic 100-min 18F-flortaucipir PET from 57 subjects positive for amyloid-ß pathology [preclinical AD (n = 16), mild cognitive impairment (MCI) due to AD (n = 16) and AD dementia (n = 25)]. Cognitively healthy subjects without amyloid-ß pathology were included as controls (n = 25). Tau propagation was modelled as an epidemic process (susceptible-infected model) on MEG-based functional networks [in alpha (8-13 Hz) and beta (13-30 Hz) bands], a structural or diffusion network, starting from the middle and inferior temporal lobe. The group-level network of the control group was used as input for the model to predict tau deposition in three stages of the AD continuum. To assess performance, model output was compared to the group-specific tau deposition patterns as measured with 18F-flortaucipir PET. We repeated the analysis by using networks of the preceding disease stage and/or using regions with most observed tau deposition during the preceding stage as seeds. In the preclinical AD stage, the functional networks predicted most of the modelled tau-PET binding potential, with best correlations between model and tau-PET [corrected amplitude envelope correlation (AEC-c) alpha C = 0.584; AEC-c beta C = 0.569], followed by the structural network (C = 0.451) and simple diffusion (C = 0.451). Prediction accuracy declined for the MCI and AD dementia stages, although the correlation between modelled tau and tau-PET binding remained highest for the functional networks (C = 0.384; C = 0.376). Replacing the control-network with the network from the preceding disease stage and/or alternative seeds improved prediction accuracy in MCI but not in the dementia stage. These results suggest that in addition to structural connections, functional connections play an important role in tau spread, and highlight that neuronal dynamics play a key role in promoting this pathological process. Aberrant neuronal communication patterns should be taken into account when identifying targets for future therapy. Our results also suggest that this process is more important in earlier disease stages (preclinical AD/MCI); possibly, in later stages, other processes may be influential.

Neurophysiological MEG markers of cognitive impairment and performance validity in multiple sclerosis.

BACKGROUND: Suboptimal performance during neuropsychological testing frequently occurs in multiple sclerosis (MS), leading to unreliable cognitive outcomes. Neurophysiological alterations correlate with MS-related cognitive impairment, but studies have not yet considered performance validity. OBJECTIVES: To investigate neurophysiological markers of cognitive impairment in MS, while explicitly addressing performance validity. METHODS: Magnetoencephalography recordings, neuropsychological assessments, and performance validity testing were obtained from 90 MS outpatients with cognitive complaints. Spectral and resting-state functional connectivity (rsFC) properties were compared between cognitively impaired (CI), cognitively preserved (CP), and suboptimally performing (SUB) patients using regression models and permutation testing. RESULTS: CI had higher power in low-frequency bands and lower power in high bands compared to CP, indicating neuronal slowing. CI also showed lower beta power compared to SUB. Overall power spectra visually differed between CI and CP, and SUB showed overlap with both groups. CI had lower rsFC than CP and SUB patients. CP and SUB patients showed no differences. CONCLUSION: Neuronal slowing and altered rsFC can be considered cognitive markers in MS. Patients who performed suboptimally showed resemblance with patients with and without cognitive impairments, and although their overall neurophysiological profile was more similar to patients without impairments, it suggests heterogeneity regarding their pathophysiology.

Regional healthy brain activity, glioma occurrence and symptomatology.

It is unclear why exactly gliomas show preferential occurrence in certain brain areas. Increased spiking activity around gliomas leads to faster tumour growth in animal models, while higher non-invasively measured brain activity is related to shorter survival in patients. However, it is unknown how regional intrinsic brain activity, as measured in healthy controls, relates to glioma occurrence. We first investigated whether gliomas occur more frequently in regions with intrinsically higher brain activity. Second, we explored whether intrinsic cortical activity at individual patients' tumour locations relates to tumour and patient characteristics. Across three cross-sectional cohorts, 413 patients were included. Individual tumour masks were created. Intrinsic regional brain activity was assessed through resting-state magnetoencephalography acquired in healthy controls and source-localized to 210 cortical brain regions. Brain activity was operationalized as: (i) broadband power; and (ii) offset of the aperiodic component of the power spectrum, which both reflect neuronal spiking of the underlying neuronal population. We additionally assessed (iii) the slope of the aperiodic component of the power spectrum, which is thought to reflect the neuronal excitation/inhibition ratio. First, correlation coefficients were calculated between group-level regional glioma occurrence, as obtained by concatenating tumour masks across patients, and group-averaged regional intrinsic brain activity. Second, intrinsic brain activity at specific tumour locations was calculated by overlaying patients' individual tumour masks with regional intrinsic brain activity of the controls and was associated with tumour and patient characteristics. As proposed, glioma preferentially occurred in brain regions characterized by higher intrinsic brain activity in controls as reflected by higher offset. Second, intrinsic brain activity at patients' individual tumour locations differed according to glioma subtype and performance status: the most malignant isocitrate dehydrogenase-wild-type glioblastoma patients had the lowest excitation/inhibition ratio at their individual tumour locations as compared to isocitrate dehydrogenase-mutant, 1p/19q-codeleted glioma patients, while a lower excitation/inhibition ratio related to poorer Karnofsky Performance Status, particularly in codeleted glioma patients. In conclusion, gliomas more frequently occur in cortical brain regions with intrinsically higher activity levels, suggesting that more active regions are more vulnerable to glioma development. Moreover, indices of healthy, intrinsic excitation/inhibition ratio at patients' individual tumour locations may capture both tumour biology and patients' performance status. These findings contribute to our understanding of the complex and bidirectional relationship between normal brain functioning and glioma growth, which is at the core of the relatively new field of 'cancer neuroscience'.

Cortical and Subcortical Changes in MEG Activity Reflect Parkinson's Progression over a Period of 7 Years.

In this study of early functional changes in Parkinson's disease (PD), we aimed to provide a comprehensive assessment of the development of changes in both cortical and subcortical neurophysiological brain activity, including their association with clinical measures of disease severity. Repeated resting-state MEG recordings and clinical assessments were obtained in the context of a unique longitudinal cohort study over a seven-year period using a multiple longitudinal design. We used linear mixed-models to analyze the relationship between neurophysiological (spectral power and functional connectivity) and clinical data. At baseline, early-stage (drug-naïve) PD patients demonstrated spectral slowing compared to healthy controls in both subcortical and cortical brain regions, most outspoken in the latter. Over time, spectral slowing progressed in strong association with clinical measures of disease progression (cognitive and motor). Global functional connectivity was not different between groups at baseline and hardly changed over time. Therefore, investigation of associations with clinical measures of disease progression were not deemed useful. An analysis of individual connections demonstrated differences between groups at baseline (higher frontal theta, lower parieto-occipital alpha2 band functional connectivity) and over time in PD patients (increase in frontal delta and theta band functional connectivity). Our results suggest that spectral measures are promising candidates in the search for non-invasive markers of both early-stage PD and of the ongoing disease process.

Cellular Substrates of Functional Network Integration and Memory in Temporal Lobe Epilepsy.

Temporal lobe epilepsy (TLE) patients are at risk of memory deficits, which have been linked to functional network disturbances, particularly of integration of the default mode network (DMN). However, the cellular substrates of functional network integration are unknown. We leverage a unique cross-scale dataset of drug-resistant TLE patients (n = 31), who underwent pseudo resting-state functional magnetic resonance imaging (fMRI), resting-state magnetoencephalography (MEG) and/or neuropsychological testing before neurosurgery. fMRI and MEG underwent atlas-based connectivity analyses. Functional network centrality of the lateral middle temporal gyrus, part of the DMN, was used as a measure of local network integration. Subsequently, non-pathological cortical tissue from this region was used for single cell morphological and electrophysiological patch-clamp analysis, assessing integration in terms of total dendritic length and action potential rise speed. As could be hypothesized, greater network centrality related to better memory performance. Moreover, greater network centrality correlated with more integrative properties at the cellular level across patients. We conclude that individual differences in cognitively relevant functional network integration of a DMN region are mirrored by differences in cellular integrative properties of this region in TLE patients. These findings connect previously separate scales of investigation, increasing translational insight into focal pathology and large-scale network disturbances in TLE.

Predicting time-resolved electrophysiological brain networks from structural eigenmodes.

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

Abnormal meta-state activation of dynamic brain networks across the Alzheimer spectrum.

The characterization of the distinct dynamic functional connectivity (dFC) patterns that activate in the brain during rest can help to understand the underlying time-varying network organization. The presence and behavior of these patterns (known as meta-states) have been widely studied by means of functional magnetic resonance imaging (fMRI). However, modalities with high-temporal resolution, such as electroencephalography (EEG), enable the characterization of fast temporally evolving meta-state sequences. Mild cognitive impairment (MCI) and dementia due to Alzheimer's disease (AD) have been shown to disrupt spatially localized activation and dFC between different brain regions, but not much is known about how they affect meta-state network topologies and their network dynamics. The main hypothesis of the study was that MCI and dementia due to AD alter normal meta-state sequences by inducing a loss of structure in their patterns and a reduction of their dynamics. Moreover, we expected that patients with MCI would display more flexible behavior compared to patients with dementia due to AD. Thus, the aim of the current study was twofold: (i) to find repeating, distinctly organized network patterns (meta-states) in neural activity; and (ii) to extract information about meta-state fluctuations and how they are influenced by MCI and dementia due to AD. To accomplish these goals, we present a novel methodology to characterize dynamic meta-states and their temporal fluctuations by capturing aspects based on both their discrete activation and the continuous evolution of their individual strength. These properties were extracted from 60-s resting-state EEG recordings from 67 patients with MCI due to AD, 50 patients with dementia due to AD, and 43 cognitively healthy controls. First, the instantaneous amplitude correlation (IAC) was used to estimate instantaneous functional connectivity with a high temporal resolution. We then extracted meta-states by means of graph community detection based on recurrence plots (RPs), both at the individual- and group-level. Subsequently, a diverse set of properties of the continuous and discrete fluctuation patterns of the meta-states was extracted and analyzed. The main novelty of the methodology lies in the usage of Louvain GJA community detection to extract meta-states from IAC-derived RPs and the extended analysis of their discrete and continuous activation. Our findings showed that distinct dynamic functional connectivity meta-states can be found on the EEG time-scale, and that these were not affected by the oscillatory slowing induced by MCI or dementia due to AD. However, both conditions displayed a loss of meta-state modularity, coupled with shorter dwell times and higher complexity of the meta-state sequences. Furthermore, we found evidence that meta-state sequencing is not entirely random; it shows an underlying structure that is partially lost in MCI and dementia due to AD. These results show evidence that AD progression is associated with alterations in meta-state switching, and a degradation of dynamic brain flexibility.

Functional brain network organization measured with magnetoencephalography predicts cognitive decline in multiple sclerosis.

BACKGROUND: Cognitive decline remains difficult to predict as structural brain damage cannot fully explain the extensive heterogeneity found between MS patients. OBJECTIVE: To investigate whether functional brain network organization measured with magnetoencephalography (MEG) predicts cognitive decline in MS patients after 5 years and to explore its value beyond structural pathology. METHODS: Resting-state MEG recordings, structural MRI, and neuropsychological assessments were analyzed of 146 MS patients, and 100 patients had a 5-year follow-up neuropsychological assessment. Network properties of the minimum spanning tree (i.e. backbone of the functional brain network) indicating network integration and overload were related to baseline and longitudinal cognition, correcting for structural damage. RESULTS: A more integrated beta band network (i.e. smaller diameter) and a less integrated delta band network (i.e. lower leaf fraction) predicted cognitive decline after 5 years (Radj2=15%), independent of structural damage. Cross-sectional analyses showed that a less integrated network (e.g. lower tree hierarchy) related to worse cognition, independent of frequency band. CONCLUSIONS: The level of functional brain network integration was an independent predictive marker of cognitive decline, in addition to the severity of structural damage. This work thereby indicates the promise of MEG-derived network measures in predicting disease progression in MS.

Mapping functional brain networks from the structural connectome: Relating the series expansion and eigenmode approaches.

Functional brain networks are shaped and constrained by the underlying structural network. However, functional networks are not merely a one-to-one reflection of the structural network. Several theories have been put forward to understand the relationship between structural and functional networks. However, it remains unclear how these theories can be unified. Two existing recent theories state that 1) functional networks can be explained by all possible walks in the structural network, which we will refer to as the series expansion approach, and 2) functional networks can be explained by a weighted combination of the eigenmodes of the structural network, the so-called eigenmode approach. To elucidate the unique or common explanatory power of these approaches to estimate functional networks from the structural network, we analysed the relationship between these two existing views. Using linear algebra, we first show that the eigenmode approach can be written in terms of the series expansion approach, i.e., walks on the structural network associated with different hop counts correspond to different weightings of the eigenvectors of this network. Second, we provide explicit expressions for the coefficients for both the eigenmode and series expansion approach. These theoretical results were verified by empirical data from Diffusion Tensor Imaging (DTI) and functional Magnetic Resonance Imaging (fMRI), demonstrating a strong correlation between the mappings based on both approaches. Third, we analytically and empirically demonstrate that the fit of the eigenmode approach to measured functional data is always at least as good as the fit of the series expansion approach, and that errors in the structural data lead to large errors of the estimated coefficients for the series expansion approach. Therefore, we argue that the eigenmode approach should be preferred over the series expansion approach. Results hold for eigenmodes of the weighted adjacency matrices as well as eigenmodes of the graph Laplacian. â€‹Taken together, these results provide an important step towards unification of existing theories regarding the structure-function relationships in brain networks.

Postoperative oscillatory brain activity as an add-on prognostic marker in diffuse glioma.

INTRODUCTION: Progression-free survival (PFS) in glioma patients varies widely, even when stratifying for known predictors (i.e. age, molecular tumor subtype, presence of epilepsy, tumor grade and Karnofsky performance status). Neuronal activity has been shown to accelerate tumor growth in an animal model, suggesting that brain activity may be valuable as a PFS predictor. We investigated whether postoperative oscillatory brain activity, assessed by resting-state magnetoencephalography is of additional value when predicting PFS in glioma patients. METHODS: We included 27 patients with grade II-IV gliomas. Each patient's oscillatory brain activity was estimated by calculating broadband power (0.5-48 Hz) in 56 epochs of 3.27 s and averaged over 78 cortical regions of the Automated Anatomical Labeling atlas. Cox proportional hazard analysis was performed to test the predictive value of broadband power towards PFS, adjusting for known predictors by backward elimination. RESULTS: Higher broadband power predicted shorter PFS after adjusting for known prognostic factors (n = 27; HR 2.56 (95% confidence interval (CI) 1.15-5.70); p = 0.022). Post-hoc univariate analysis showed that higher broadband power also predicted shorter overall survival (OS; n = 38; HR 1.88 (95% CI 1.00-3.54); p = 0.038). CONCLUSIONS: Our findings suggest that postoperative broadband power is of additional value in predicting PFS beyond already known predictors.

Increased segregation of functional networks in developing brains.

A growing literature conceptualises typical brain development from a network perspective. However, largely due to technical and methodological challenges inherent in paediatric functional neuroimaging, there remains an important gap in our knowledge regarding the typical development of functional brain networks in "preschool" childhood (i.e., children younger than 6 years of age). In this study, we recorded brain oscillatory activity using age-appropriate magnetoencephalography in 24 children, including 14 preschool children aged from 4 to 6 years and 10 school children aged from 7 to 12 years. We compared the topology of the resting-state brain networks in these children, estimated using minimum spanning tree (MST) constructed from phase synchrony between beamformer-reconstructed time-series, with that of 24 adults. Our results show that during childhood the MST topology shifts from a star-like (centralised) toward a more line-like (de-centralised) configuration, indicating the functional brain networks become increasingly segregated. In addition, the increasing global network segregation is frequency-independent and accompanied by decreases in centrality (or connectedness) of cortical regions with age, especially in areas of the default mode network. We propose a heuristic MST model of "network space", which posits a clear developmental trajectory for the emergence of complex brain networks. Our results not only revealed topological reorganisation of functional networks across multiple temporal and spatial scales in childhood, but also fill a gap in the literature regarding neurophysiological mechanisms of functional brain maturation during the preschool years of childhood.

Tracking dynamic brain networks using high temporal resolution MEG measures of functional connectivity.

Fluctuations in functional interactions between brain regions typically occur at the millisecond time scale. Conventional connectivity metrics are not adequately time-resolved to detect such fast fluctuations in functional connectivity. At the same time, attempts to use conventional metrics in a time-resolved manner usually come with the selection of sliding windows of fixed arbitrary length. In the current work, we evaluated the use of high temporal resolution metrics of functional connectivity in conjunction with non-negative tensor factorisation to detect fast fluctuations in connectivity and temporally evolving subnetworks. To this end, we used the phase difference derivative, wavelet coherence, and we also introduced a new metric, the instantaneous amplitude correlation. In order to deal with the inherently noisy nature of magnetoencephalography data and large datasets, we make use of recurrence plots and we used pair-wise orthogonalisation to avoid spurious estimates of functional connectivity due to signal leakage. Firstly, metrics were evaluated in the context of dynamically coupled neural mass models in the presence and absence of delays and also compared to conventional static metrics with fixed sliding windows. Simulations showed that these high temporal resolution metrics outperformed conventional static connectivity metrics. Secondly, the sensitivity of the metrics to fluctuations in connectivity was analysed in post-movement beta rebound magnetoencephalography data, which showed time locked sensorimotor subnetworks that modulated with the post-movement beta rebound. Finally, sensitivity of the metrics was evaluated in resting-state magnetoencephalography, showing similar spatial patterns across metrics, thereby indicating the robustness of the current analysis. The current methods can be applied in cognitive experiments that involve fast modulations in connectivity in relation to cognition. In addition, these methods could also be used as input to temporal graph analysis to further characterise the rapid fluctuation in brain network topology.

A systematic review of MEG-based studies in Parkinson's disease: The motor system and beyond.

Parkinson's disease (PD) is accompanied by functional changes throughout the brain, including changes in the electromagnetic activity recorded with magnetoencephalography (MEG). An integrated overview of these changes, its relationship with clinical symptoms, and the influence of treatment is currently missing. Therefore, we systematically reviewed the MEG studies that have examined oscillatory activity and functional connectivity in the PD-affected brain. The available articles could be separated into motor network-focused and whole-brain focused studies. Motor network studies revealed PD-related changes in beta band (13-30 Hz) neurophysiological activity within and between several of its components, although it remains elusive to what extent these changes underlie clinical motor symptoms. In whole-brain studies PD-related oscillatory slowing and decrease in functional connectivity correlated with cognitive decline and less strongly with other markers of disease progression. Both approaches offer a different perspective on PD-specific disease mechanisms and could therefore complement each other. Combining the merits of both approaches will improve the setup and interpretation of future studies, which is essential for a better understanding of the disease process itself and the pathophysiological mechanisms underlying specific PD symptoms, as well as for the potential to use MEG in clinical care.

Resting-state MEG measurement of functional activation as a biomarker for cognitive decline in MS.

BACKGROUND: Neurophysiological measures of brain function, such as magnetoencephalography (MEG), are widely used in clinical neurology and have strong relations with cognitive impairment and dementia but are still underdeveloped in multiple sclerosis (MS). OBJECTIVES: To demonstrate the value of clinically applicable MEG-measures in evaluating cognitive impairment in MS. METHODS: In eyes-closed resting-state, MEG data of 83 MS patients and 34 healthy controls (HCs) peak frequencies and relative power of six canonical frequency bands for 78 cortical and 10 deep gray matter (DGM) areas were calculated. Linear regression models, correcting for age, gender, and education, assessed the relation between cognitive performance and MEG biomarkers. RESULTS: Increased alpha1 and theta power was strongly associated with impaired cognition in patients, which differed between cognitively impaired (CI) patients and HCs in bilateral parietotemporal cortices. CI patients had a lower peak frequency than HCs. Oscillatory slowing was also widespread in the DGM, most pronounced in the thalamus. CONCLUSION: There is a clinically relevant slowing of neuronal activity in MS patients in parietotemporal cortical areas and the thalamus, strongly related to cognitive impairment. These measures hold promise for the application of resting-state MEG as a biomarker for cognitive disturbances in MS in a clinical setting.

Direction of information flow in large-scale resting-state networks is frequency-dependent.

Normal brain function requires interactions between spatially separated, and functionally specialized, macroscopic regions, yet the directionality of these interactions in large-scale functional networks is unknown. Magnetoencephalography was used to determine the directionality of these interactions, where directionality was inferred from time series of beamformer-reconstructed estimates of neuronal activation, using a recently proposed measure of phase transfer entropy. We observed well-organized posterior-to-anterior patterns of information flow in the higher-frequency bands (alpha1, alpha2, and beta band), dominated by regions in the visual cortex and posterior default mode network. Opposite patterns of anterior-to-posterior flow were found in the theta band, involving mainly regions in the frontal lobe that were sending information to a more distributed network. Many strong information senders in the theta band were also frequent receivers in the alpha2 band, and vice versa. Our results provide evidence that large-scale resting-state patterns of information flow in the human brain form frequency-dependent reentry loops that are dominated by flow from parieto-occipital cortex to integrative frontal areas in the higher-frequency bands, which is mirrored by a theta band anterior-to-posterior flow.

Interactions between frontal and posterior oscillatory dynamics support adjustment of stimulus processing during reinforcement learning.

Reinforcement learning (RL) in humans is subserved by a network of striatal and frontal brain areas. The electrophysiological signatures of feedback evaluation are increasingly well understood, but how those signatures relate to the use of feedback to guide subsequent behavioral adjustment remains unclear. One mechanism for post-feedback behavioral optimization is the modulation of sensory processing. We used source-reconstructed MEG to test whether feedback affects the interactions between sources of oscillatory activity in the learning network and task-relevant stimulus-processing areas. Participants performed a probabilistic RL task in which they learned associations between colored faces and response buttons using trial-and-error feedback. Delta-band (2-4 Hz) and theta-band (4-8 Hz) power in multiple frontal regions were sensitive to feedback valence. Low and high beta-band power (12-20 and 20-30 Hz) in occipital, parietal, and temporal regions differentiated between color and face information. Consistent with our hypothesis, single-trial power-power correlations between frontal and posterior-sensory areas were modulated by the interaction between feedback valence and the relevant stimulus characteristic (color versus identity). These results suggest that long-range oscillatory coupling supports post-feedback updating of stimulus processing.

Comparing multilayer brain networks between groups: Introducing graph metrics and recommendations.

There is an increasing awareness of the advantages of multi-modal neuroimaging. Networks obtained from different modalities are usually treated in isolation, which is however contradictory to accumulating evidence that these networks show non-trivial interdependencies. Even networks obtained from a single modality, such as frequency-band specific functional networks measured from magnetoencephalography (MEG) are often treated independently. Here, we discuss how a multilayer network framework allows for integration of multiple networks into a single network description and how graph metrics can be applied to quantify multilayer network organisation for group comparison. We analyse how well-known biases for single layer networks, such as effects of group differences in link density and/or average connectivity, influence multilayer networks, and we compare four schemes that aim to correct for such biases: the minimum spanning tree (MST), effective graph resistance cost minimisation, efficiency cost optimisation (ECO) and a normalisation scheme based on singular value decomposition (SVD). These schemes can be applied to the layers independently or to the multilayer network as a whole. For correction applied to whole multilayer networks, only the SVD showed sufficient bias correction. For correction applied to individual layers, three schemes (ECO, MST, SVD) could correct for biases. By using generative models as well as empirical MEG and functional magnetic resonance imaging (fMRI) data, we further demonstrated that all schemes were sensitive to identify network topology when the original networks were perturbed. In conclusion, uncorrected multilayer network analysis leads to biases. These biases may differ between centres and studies and could consequently lead to unreproducible results in a similar manner as for single layer networks. We therefore recommend using correction schemes prior to multilayer network analysis for group comparisons.

Consistency of magnetoencephalographic functional connectivity and network reconstruction using a template versus native MRI for co-registration.

INTRODUCTION: Studies using functional connectivity and network analyses based on magnetoencephalography (MEG) with source localization are rapidly emerging in neuroscientific literature. However, these analyses currently depend on the availability of costly and sometimes burdensome individual MR scans for co-registration. We evaluated the consistency of these measures when using a template MRI, instead of native MRI, for the analysis of functional connectivity and network topology. METHODS: Seventeen healthy participants underwent resting-state eyes-closed MEG and anatomical MRI. These data were projected into source space using an atlas-based peak voxel and a centroid beamforming approach either using (1) participants' native MRIs or (2) the Montreal Neurological Institute's template. For both methods, time series were reconstructed from 78 cortical atlas regions. Relative power was determined in six classical frequency bands per region and globally averaged. Functional connectivity (phase lag index) between each pair of regions was calculated. The adjacency matrices were then used to reconstruct functional networks, of which regional and global metrics were determined. Intraclass correlation coefficients were calculated and Bland-Altman plots were made to quantify the consistency and potential bias of the use of template versus native MRI. RESULTS: Co-registration with the template yielded largely consistent relative power, connectivity, and network estimates compared to native MRI. DISCUSSION: These findings indicate that there is no (systematic) bias or inconsistency between template and native MRI co-registration of MEG. They open up possibilities for retrospective and prospective analyses to MEG datasets in the general population that have no native MRIs available. Hum Brain Mapp, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc. Hum Brain Mapp 39:104-119, 2018. © 2017 Wiley Periodicals, Inc.