Virtual neural network-guided optimization of non-invasive brain stimulation in Alzheimer's disease.

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation technique with potential for counteracting disrupted brain network activity in Alzheimer's disease (AD) to improve cognition. However, the results of tDCS studies in AD have been variable due to different methodological choices such as electrode placement. To address this, a virtual brain network model of AD was used to explore tDCS optimization. We compared a large, representative set of virtual tDCS intervention setups, to identify the theoretically optimized tDCS electrode positions for restoring functional network features disrupted in AD. We simulated 20 tDCS setups using a computational dynamic network model of 78 neural masses coupled according to human structural topology. AD network damage was simulated using an activity-dependent degeneration algorithm. Current flow modeling was used to estimate tDCS-targeted cortical regions for different electrode positions, and excitability of the pyramidal neurons of the corresponding neural masses was modulated to simulate tDCS. Outcome measures were relative power spectral density (alpha bands, 8-10 Hz and 10-13 Hz), total spectral power, posterior alpha peak frequency, and connectivity measures phase lag index (PLI) and amplitude envelope correlation (AEC). Virtual tDCS performance varied, with optimized strategies improving all outcome measures, while others caused further deterioration. The best performing setup involved right parietal anodal stimulation, with a contralateral supraorbital cathode. A clear correlation between the network role of stimulated regions and tDCS success was not observed. This modeling-informed approach can guide and perhaps accelerate tDCS therapy development and enhance our understanding of tDCS effects. Follow-up studies will compare the general predictions to personalized virtual models and validate them with tDCS-magnetoencephalography (MEG) in a clinical AD patient cohort.

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.

Combining advanced MRI and EEG techniques better explains long-term motor outcome after very preterm birth.

BACKGROUND: Preterm born children are at high risk for adverse motor neurodevelopment. The aim of this study was to establish the relationship between motor outcome and advanced magnetic resonance imaging (MRI) and electroencephalography (EEG) measures. METHODS: In a prospective cohort study of 64 very preterm born children, the motor outcome was assessed at 9.83 (SD 0.70) years. Volumetric MRI, diffusion tensor imaging (DTI), and EEG were acquired at 10.85 (SD 0.49) years. We investigated associations between motor outcome and brain volumes (white matter, deep gray matter, cerebellum, and ventricles), white matter integrity (fractional anisotropy and mean, axial and radial diffusivity), and brain activity (upper alpha (A2) functional connectivity and relative A2 power). The independence of associations with motor outcome was investigated with a final model. For each technique, the measure with the strongest association was selected to avoid multicollinearity. RESULTS: Ventricular volume, radial diffusivity, mean diffusivity, relative A2 power, and A2 functional connectivity were significantly correlated to motor outcome. The final model showed that ventricular volume and relative A2 power were independently associated with motor outcome (B = -9.42 × 10-5, p = 0.027 and B = 28.9, p = 0.007, respectively). CONCLUSIONS: This study suggests that a lasting interplay exists between brain structure and function that might underlie motor outcome at school age. IMPACT: This is the first study that investigates the relationships between motor outcome and brain volumes, DTI, and brain function in preterm born children at school age. Ventricular volume and relative upper alpha power on EEG have an independent relation with motor outcome in preterm born children at school age. This suggests that there is a lasting interplay between structure and function that underlies adverse motor outcome.

Neonatal quantitative electroencephalography and long-term outcomes: a systematic review.

AIM: To evaluate quantitative electroencephalogram (EEG) measures as predictors of long-term neurodevelopmental outcome in infants with a postconceptional age below 46 weeks, including typically developing infants born at term, infants with heterogeneous underlying pathologies, and infants born preterm. METHOD: A comprehensive search was performed using PubMed, Embase, and Web of Science from study inception up to 8th January 2021. Studies that examined associations between neonatal quantitative EEG measures, based on conventional and amplitude-integrated EEG, and standardized neurodevelopmental outcomes at 2 years of age or older were reviewed. Significant associations between neonatal quantitative EEG and long-term outcome measures were grouped into one or more of the following categories: cognitive outcome; motor outcome; composite scores; and other standardized outcome assessments. RESULTS: Twenty-four out of 1740 studies were included. Multiple studies showed that conventional EEG-based absolute power in the delta, theta, alpha, and beta frequency bands and conventional and amplitude-integrated EEG-related amplitudes were positively associated with favourable long-term outcome across several domains, including cognition and motor performance. Furthermore, a lower presence of discontinuous background pattern was also associated with favourable outcomes. However, interpretation of the results is limited by heterogeneity in study design and populations. INTERPRETATION: Neonatal quantitative EEG measures may be used as prognostic biomarkers to identify those infants who will develop long-term difficulties and who might benefit from early interventions.

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.

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.

Modern network science of neurological disorders.

Modern network science has revealed fundamental aspects of normal brain-network organization, such as small-world and scale-free patterns, hierarchical modularity, hubs and rich clubs. The next challenge is to use this knowledge to gain a better understanding of brain disease. Recent developments in the application of network science to conditions such as Alzheimer's disease, multiple sclerosis, traumatic brain injury and epilepsy have challenged the classical concept of neurological disorders being either 'local' or 'global', and have pointed to the overload and failure of hubs as a possible final common pathway in neurological disorders.

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.

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.