Search

Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Blood-level oxygenation-dependent (BOLD) functional magnetic resonance imaging (fMRI) is the most common modality to study functional connectivity in the human brain. Most research to date has focused on connectivity between pairs of brain regions. However, attention has recently turned towards connectivity involving more than two regions, that is, higher-order connectivity. It is not yet clear how higher-order connectivity can best be quantified. The measures that are currently in use cannot distinguish between pairwise (i.e., second-order) and higher-order connectivity. We show that genuine higher-order connectivity can be quantified by using multivariate cumulants. We explore the use of multivariate cumulants for quantifying higher-order connectivity and the performance of block bootstrapping for statistical inference. In particular, we formulate a generative model for fMRI signals exhibiting higher-order connectivity and use it to assess bias, standard errors, and detection probabilities. Application to resting-state fMRI data from the Human Connectome Project demonstrates that spontaneous fMRI signals are organized into higher-order networks that are distinct from second-order resting-state networks. Application to a clinical cohort of patients with multiple sclerosis further demonstrates that cumulants can be used to classify disease groups and explain behavioral variability. Hence, we present a novel framework to reliably estimate genuine higher-order connectivity in fMRI data which can be used for constructing hyperedges, and finally, which can readily be applied to fMRI data from populations with neuropsychiatric disease or cognitive neuroscientific experiments., (© 2024 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

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

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., (© 2022 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

Longitudinal analyses of magnetoencephalography (MEG) data are essential for a full understanding of the pathophysiology of brain diseases and the development of brain activity over time. However, time-dependent factors, such as the recording environment and the type of MEG recording system may affect such longitudinal analyses. We hypothesized that, using source-space analysis, hardware and software differences between two recordings systems may be overcome, with the aim of finding consistent neurophysiological results. We studied eight healthy subjects who underwent three consecutive MEG recordings over 7 years, using two different MEG recordings systems; a 151-channel VSM-CTF system for the first two time points and a 306-channel Elekta Vectorview system for the third time point. We assessed the within (longitudinal) and between-subject (cross-sectional) consistency of power spectra and functional connectivity matrices. Consistency of within-subject spectral power and functional connectivity matrices was good and was not significantly different when using different MEG recording systems as compared to using the same system. Importantly, we confirmed that within-subject consistency values were higher than between-subject values. We demonstrated consistent neurophysiological findings in healthy subjects over a time span of seven years, despite using data recorded on different MEG systems and different implementations of the analysis pipeline., (© 2021. The Author(s).)

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., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

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.

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

Functional networks obtained from magnetoencephalography (MEG) from different frequency bands show distinct spatial patterns. It remains to be elucidated how distinct spatial patterns in MEG networks emerge given a single underlying structural network. Recent work has suggested that the eigenmodes of the structural network might serve as a basis set for functional network patterns in the case of functional MRI. Here, we take this notion further in the context of frequency band specific MEG networks. We show that a selected set of eigenmodes of the structural network can predict different frequency band specific networks in the resting state, ranging from delta (1-4 Hz) to the high gamma band (40-70 Hz). These predictions outperform predictions based from surrogate data, suggesting a genuine relationship between eigenmodes of the structural network and frequency specific MEG networks. We then show that the relevant set of eigenmodes can be excited in a network of neural mass models using linear stability analysis only by including delays. Excitation of an eigenmode in this context refers to a dynamic instability of a network steady state to a spatial pattern with a corresponding coherent temporal oscillation. Simulations verify the results from linear stability analysis and suggest that theta, alpha and beta band networks emerge very near to the bifurcation. The delta and gamma bands in the resting state emerges further away from the bifurcation. These results show for the first time how delayed interactions can excite the relevant set of eigenmodes that give rise to frequency specific functional connectivity patterns., (Copyright © 2018 Elsevier Inc. All rights reserved.)

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., (Copyright © 2017 Elsevier Inc. All rights reserved.)

One of the challenges of brain network analysis is to directly compare network organization between subjects, irrespective of the number or strength of connections. In this study, we used minimum spanning tree (MST; a unique, acyclic subnetwork with a fixed number of connections) analysis to characterize the human brain network to create an empirical reference network. Such a reference network could be used as a null model of connections that form the backbone structure of the human brain. We analyzed the MST in three diffusion-weighted imaging datasets of healthy adults. The MST of the group mean connectivity matrix was used as the empirical null-model. The MST of individual subjects matched this reference MST for a mean 58%-88% of connections, depending on the analysis pipeline. Hub nodes in the MST matched with previously reported locations of hub regions, including the so-called rich club nodes (a subset of high-degree, highly interconnected nodes). Although most brain network studies have focused primarily on cortical connections, cortical-subcortical connections were consistently present in the MST across subjects. Brain network efficiency was higher when these connections were included in the analysis, suggesting that these tracts may be utilized as the major neural communication routes. Finally, we confirmed that MST characteristics index the effects of brain aging. We conclude that the MST provides an elegant and straightforward approach to analyze structural brain networks, and to test network topological features of individual subjects in comparison to empirical null models., (© 2018 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.)

To understand the heterogeneity of functional connectivity results reported in the literature, we analyzed the separate effects of grey and white matter damage on functional connectivity and networks in multiple sclerosis. For this, we employed a biophysical thalamo-cortical model consisting of interconnected cortical and thalamic neuronal populations, informed and amended by empirical diffusion MRI tractography data, to simulate functional data that mimic neurophysiological signals. Grey matter degeneration was simulated by decreasing within population connections and white matter degeneration by lowering between population connections, based on lesion predilection sites in multiple sclerosis. For all simulations, functional connectivity and functional network organization are quantified by phase synchronization and network integration, respectively. Modeling results showed that both cortical and thalamic grey matter damage induced a global increase in functional connectivity, whereas white matter damage induced an initially increased connectivity followed by a global decrease. Both white and especially grey matter damage, however, induced a decrease in network integration. These empirically informed simulations show that specific topology and timing of structural damage are nontrivial aspects in explaining functional abnormalities in MS. Insufficient attention to these aspects likely explains contradictory findings in multiple sclerosis functional imaging studies so far., (© 2018 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.)

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., (Copyright © 2017 Elsevier Inc. All rights reserved.)

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., (© 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.)

Resective neurosurgery carries the risk of postoperative cognitive deterioration. The concept of 'hub (over)load', caused by (over)use of the most important brain regions, has been theoretically postulated in relation to symptomatology and neurological disease course, but lacks experimental confirmation. We investigated functional hub load and postsurgical cognitive deterioration in patients undergoing lesion resection. Patients (n = 28) underwent resting-state magnetoencephalography and neuropsychological assessments preoperatively and 1-year after lesion resection. We calculated stationary hub load score (SHub) indicating to what extent brain regions linked different subsystems; high SHub indicates larger processing pressure on hub regions. Dynamic hub load score (DHub) assessed its variability over time; low values, particularly in combination with high SHub values, indicate increased load, because of consistently high usage of hub regions. Hypothetically, increased SHub and decreased DHub relate to hub overload and thus poorer/deteriorating cognition. Between time points, deteriorating verbal memory performance correlated with decreasing upper alpha DHub. Moreover, preoperatively low DHub values accurately predicted declining verbal memory performance. In summary, dynamic hub load relates to cognitive functioning in patients undergoing lesion resection: postoperative cognitive decline can be tracked and even predicted using dynamic hub load, suggesting it may be used as a prognostic marker for tailored treatment planning., Competing Interests: M.K. has received consultancy fees from Hoffmann-La Roche, outside the submitted work. E.W.S.C., E.v.D., P.T., P.C.d.W.H., J.C.B., J.J.G.G., J.C.R., C.J.S. and L.D. report no conflicts of interest.

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

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

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.

Communication between brain regions is still insufficiently understood. Applying concepts from network science has shown to be successful in gaining insight in the functioning of the brain. Recent work has implicated that especially shortest paths in the structural brain network seem to play a major role in the communication within the brain. So far, for the functional brain network, only the average length of the shortest paths has been analyzed. In this article, we propose to construct the union of shortest path trees (USPT) as a new topology for the functional brain network. The minimum spanning tree, which has been successful in a lot of recent studies to comprise important features of the functional brain network, is always included in the USPT. After interpreting the link weights of the functional brain network as communication probabilities, the USPT of this network can be uniquely defined. Using data from magnetoencephalography, we applied the USPT as a method to find differences in the network topology of multiple sclerosis patients and healthy controls. The new concept of the USPT of the functional brain network also allows interesting interpretations and may represent the highways of the brain.

Introduction: Gray matter (GM) atrophy is common in multiple sclerosis (MS), but the relationship with white matter (WM) pathology is largely unknown. Some studies found a co-occurrence in specific systems, but a regional analysis across the brain in different clinical phenotypes is necessary to further understand the disease mechanism underlying GM atrophy in MS. Therefore, we investigated the association between regional GM atrophy and pathology in anatomically connected WM tracts., Methods: Conventional and diffusion tensor imaging was performed at 3T in 208 patients with long-standing MS and 60 healthy controls. Deep and cortical GM regions were segmented and quantified, and both lesion volumes and average normal appearing WM fractional anisotropy of their associated tracts were derived using an atlas obtained by probabilistic tractography in the controls. Linear regression was then performed to quantify the amount of regional GM atrophy that can be explained by WM pathology in the connected tract., Results: MS patients showed extensive deep and cortical GM atrophy. Cortical atrophy was particularly present in frontal and temporal regions. Pathology in connected WM tracts statistically explained both regional deep and cortical GM atrophy in relapsing-remitting (RR) patients, but only deep GM atrophy in secondary-progressive (SP) patients., Conclusion: In RRMS patients, both deep and cortical GM atrophy were associated with pathology in connected WM tracts. In SPMS patients, only regional deep GM atrophy could be explained by pathology in connected WM tracts. This suggests that in SPMS patients cortical GM atrophy and WM damage are (at least partly) independent disease processes., (© 2015 Wiley Periodicals, Inc.)

Both gray matter atrophy and disruption of functional networks are important predictors for physical disability and cognitive impairment in multiple sclerosis (MS), yet their relationship is poorly understood. Graph theory provides a modality invariant framework to analyze patterns of gray matter morphology and functional coactivation. We investigated, how gray matter and functional networks were affected within the same MS sample and examined their interrelationship. Magnetic resonance imaging and magnetoencephalography (MEG) were performed in 102 MS patients and 42 healthy controls. Gray matter networks were computed at the group-level based on cortical thickness correlations between 78 regions across subjects. MEG functional networks were computed at the subject level based on the phase-lag index between time-series of regions in source-space. In MS patients, we found a more regular network organization for structural covariance networks and for functional networks in the theta band, whereas we found a more random network organization for functional networks in the alpha2 band. Correlation analysis revealed a positive association between covariation in thickness and functional connectivity in especially the theta band in MS patients, and these results could not be explained by simple regional gray matter thickness measurements. This study is a first multimodal graph analysis in a sample of MS patients, and our results suggest that a disruption of gray matter network topology is important to understand alterations in functional connectivity in MS as regional gray matter fails to take into account the inherent connectivity structure of the brain., (© 2014 Wiley Periodicals, Inc.)

Background: The majority of patients with multiple sclerosis (MS) present with spinal cord pathology. Spinal cord atrophy is thought to be a marker of disease severity, but in long-disease duration its relation to brain pathology and clinical disability is largely unknown., Objective: Our aim was to investigate mean upper cervical cord area (MUCCA) in patients with long-standing MS and assess its relation to brain magnetic resonance imaging (MRI) measures and clinical disability., Methods: MUCCA was measured in 196 MS patients and 55 healthy controls using 3DT1-weighted cervical images obtained at 3T MRI. Clinical disability was measured using the Expanded Disability Status Scale (EDSS), Nine-Hole-Peg test (9-HPT), and 25 feet Timed Walk Test (TWT). Stepwise linear regression was performed to assess the association between MUCCA and MRI measures, and between MUCCA and clinical disability., Results: MUCCA was smaller (mean 11.7%) in MS patients compared with healthy controls (72.56±9.82 and 82.24±7.80 mm2 respectively; p<0.001), most prominently in male patients. MUCCA was associated with normalized brain volume, and number of cervical cord lesions. MUCCA was independently associated with EDSS, TWT, and 9-HPT., Conclusion: MUCCA was reduced in MS patients compared with healthy controls. It provides a relevant marker for clinical disability in long-standing disease, independent of other MRI measures., (© The Author(s), 2014.)

Purpose: To identify the measures of focal and diffuse white matter (WM) abnormalities that are related to whole-brain, deep, and cortical gray matter (GM) atrophy in long-standing multiple sclerosis (MS)., Materials and Methods: The institutional review board approved the study; all subjects gave written informed consent. Magnetic resonance (MR) imaging was performed at 3 T in 208 patients with MS of long-standing duration (disease duration ≥ 10 years) and in 60 healthy control subjects. Normalized GM volume (NGMV), normalized WM volume (NWMV), normalized deep GM volume (NDGMV), cortical thickness, and normalized lesion volume (NLV) were quantified. Tissue integrity of normal-appearing WM (NAWM) and lesions was measured by using diffusion-tensor MR imaging. Multivariate associations between measures of GM atrophy and WM abnormalities were assessed in the patient group by using multiple linear regression., Results: NGMV, NDGMV, and cortical thickness were reduced in patients with MS (all P < .001). The final model for NGMV consisted of NWMV, NLV, and patient age and sex (adjusted R(2) = 0.58, P < .001). NWMV, NLV, and patient sex were the explanatory variables for NDGMV (adjusted R(2) = 0.75, P < .001). The model for cortical thickness consisted of fractional anisotropy of NAWM, NLV, and patient age and sex (adjusted R(2) = 0.32, P < .001). The relationship between GM atrophy and WM abnormalities was weaker in primary and secondary progressive disease than in relapsing-remitting disease., Conclusion: Whole-brain and deep GM atrophy were particularly explained by WM atrophy and lesion volume, while cortical atrophy was associated with NAWM integrity loss. The weaker relationship between GM atrophy and WM abnormalities in patients with progressive disease might indicate a more independent neurodegenerative disease process in these patients.

The EEG is an accessible tool for detecting encephalopathy, which usually manifests as delirium and sometimes as coma. Several disturbances have been described in the EEG of patients with encephalopathy, including diffuse slowing and periodic discharges. The pathophysiology of these EEG alterations, however, is poorly understood. This article shows that simulating activity of large populations of neurons, using neural mass models and neural network analysis, may increase our understanding of EEG disturbances in encephalopathy. We provide a brief introduction on the concepts of neural mass modeling and graph theoretical network analysis, and insights from this approach in previous work on neurologic disease, with a focus on encephalopathy. Finally, we speculate how anatomically coupled neural mass modeling combined with network analysis could provide new insights in pathophysiology of encephalopathy.

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. • Sliding window approaches for dynamic connectivity come with selection of fixed and arbitrary window lengths. • We evaluate the use of high temporal resolution metrics of connectivity for electrophysiological data. • We evaluate two existing measures: the phase difference derivative and the wavelet coherence. • We introduce one new high temporal resolution metric: the instantaneous amplitude correlation. • All metrics can detect genuine fluctuations in dynamic connectivity as was shown in simulations, task- and resting-state data. [ABSTRACT FROM AUTHOR]

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