Multimodal description of whole brain connectivity: A comparison of resting state MEG, fMRI, and DWI

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Brain identification of IBS patients based on GBDT and multiple imaging techniques

The brain biomarker of irritable bowel syndrome (IBS) patients is still lacking. The study aims to explore a new technology studying the brain alterations of IBS patients based on multi-source brain data. In the study, a decision-level fusion method based on gradient boosting decision tree (GBDT) was proposed. Next, 100 healthy subjects were used to validate the effectiveness of the method. Finally, the identification of brain alterations and the pain evaluation in IBS patients were carried out by the fusion method based on the resting-state fMRI and DWI for 46 patients and 46 controls selected randomly from 100 healthy subjects. The results showed that the method can achieve good classification between IBS patients and controls (accuracy = 95%) and pain evaluation of IBS patients (mean absolute error = 0.1977). Moreover, both the gain-based and the permutation-based evaluation instead of statistical analysis showed that left cingulum bundle contributed most significantly to the classification, and right precuneus contributed most significantly to the evaluation of abdominal pain intensity in the IBS patients. The differences seem to suggest a probable but unexplored separation about the central regions between the identification and progression of IBS. This finding may provide one new thought and technology for brain alteration related to IBS.

Connectivity of neural signals to the primary motor area during preparatory periods for movement following external and internal cues

PURPOSE

We investigated the connectivity of neural signals from movement-related cortical areas to the primary motor area (M1) in the hemisphere contralateral to the movement side during the period of movement-related magnetic fields before movement.

MATERIALS AND METHODS

Participants were 13 healthy adults, and nerual signals were recorded using magnetoencephalography. Spontaneous extension of the right wrist was performed at the participant's own pace and following a visual cue in internal (IC) and external (EC) cue tasks. The connectivity of neural signals to M1 from each movement-related motor area was assessed by Granger causality analysis (GCA). The GCA was performed on the neural activity elicited in a frequency band between 7.8 and 46.9 Hz during the pre-movement periods, which occurred durng the readiness field (RF) and the negative slope prime (NSp). F-values, as connectivity values obtained by GCA, were compared between the EC and IC cue tasks.

RESULTS

For NSp periods, the connectivity of neural signals from the left superior frontal area (SF-L) to M1 was dominant in the IC task, whereas that from the left superior parietal area (SP-L) to M1 was dominant in the EC task. The F value in the GCA from SP-L to M1 was greater in the EC task during RF than in the IC task during equivalent periods.

CONSLUSIONS

In the present study, there were differences in the connectivity of neural signals to M1 between IC and EC tasks. The present results suggested that the pattern of pre-movement neural activity that resulted in a movement was not uniform but differed between movement tasks just before the movement.

Similarity between structural and proxy estimates of brain connectivity

Functional magnetic resonance and diffusion weighted imaging have so far made a major contribution to delineation of the brain connectome at the macroscale. While functional connectivity (FC) was shown to be related to structural connectivity (SC) to a certain degree, their spatial overlap is unknown. Even less clear are relations of SC with estimates of connectivity from inter-subject covariance of regional F18-fluorodeoxyglucose uptake (FDGcov) and grey matter volume (GMVcov). Here, we asked to what extent SC underlies three proxy estimates of brain connectivity: FC, FDGcov and GMVcov. Simultaneous PET/MR acquisitions were performed in 56 healthy middle-aged individuals. Similarity between four networks was assessed using Spearman correlation and convergence ratio (CR), a measure of spatial overlap. Spearman correlation coefficient was 0.27 for SC-FC, 0.40 for SC-FDGcov, and 0.15 for SC-GMVcov. Mean CRs were 51% for SC-FC, 48% for SC-FDGcov, and 37% for SC-GMVcov. These results proved to be reproducible and robust against image processing steps. In sum, we found a relevant similarity of SC with FC and FDGcov, while GMVcov consistently showed the weakest similarity. These findings indicate that white matter tracts underlie FDGcov to a similar degree as FC, supporting FDGcov as estimate of functional brain connectivity.

A Model to Study Time Lagged Interactions, Source Connectivity and Source Activities Using Multi-channel EEG

A computational model to examine time lagged interactions; identify number of interacting pairs of neuronal sources; and determine source activities from multi-channel EEG measurements is described. It is based on the imaginary part of the cross spectrum between the EEG channels. The imaginary part of the cross spectrum between the EEG channels provides the most suitable property that reflects the presence of interacting sources. The model assumes that not all sources are activated simultaneously and that there is a time lag amongst some of them. A new analytical expression derived for the imaginary part of cross spectrum between channels shows that it is different from the zero lag case. A method is then proposed to identify time lag interactions, by studying its variation as a function of frequency. Assuming pair wise interaction between sources, the model shows that simultaneous diagonalization at different frequencies of symmetric matrices formed by multiplying the anti-symmetric matrix of the imaginary part of cross spectrum with its transpose will provide information on the number of interacting source pairs as a function of frequency. The matrix that simultaneously diagonalizes all the symmetric matrices is identified as the mixing matrix. This can be used to obtain the source activities.

Integrating multimodal and multiscale connectivity blueprints of the human cerebral cortex in health and disease

The brain is composed of disparate neural populations that communicate and interact with one another. Although fiber bundles, similarities in molecular architecture, and synchronized neural activity all reflect how brain regions potentially interact with one another, a comprehensive study of how all these interregional relationships jointly reflect brain structure and function remains missing. Here, we systematically integrate 7 multimodal, multiscale types of interregional similarity ("connectivity modes") derived from gene expression, neurotransmitter receptor density, cellular morphology, glucose metabolism, haemodynamic activity, and electrophysiology in humans. We first show that for all connectivity modes, feature similarity decreases with distance and increases when regions are structurally connected. Next, we show that connectivity modes exhibit unique and diverse connection patterns, hub profiles, spatial gradients, and modular organization. Throughout, we observe a consistent primacy of molecular connectivity modes-namely correlated gene expression and receptor similarity-that map onto multiple phenomena, including the rich club and patterns of abnormal cortical thickness across 13 neurological, psychiatric, and neurodevelopmental disorders. Finally, to construct a single multimodal wiring map of the human cortex, we fuse all 7 connectivity modes and show that the fused network maps onto major organizational features of the cortex including structural connectivity, intrinsic functional networks, and cytoarchitectonic classes. Altogether, this work contributes to the integrative study of interregional relationships in the human cerebral cortex.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Structural basis of envelope and phase intrinsic coupling modes in the cerebral cortex

Intrinsic coupling modes (ICMs) can be observed in ongoing brain activity at multiple spatial and temporal scales. Two families of ICMs can be distinguished: phase and envelope ICMs. The principles that shape these ICMs remain partly elusive, in particular their relation to the underlying brain structure. Here we explored structure-function relationships in the ferret brain between ICMs quantified from ongoing brain activity recorded with chronically implanted micro-ECoG arrays and structural connectivity (SC) obtained from high-resolution diffusion MRI tractography. Large-scale computational models were used to explore the ability to predict both types of ICMs. Importantly, all investigations were conducted with ICM measures that are sensitive or insensitive to volume conduction effects. The results show that both types of ICMs are significantly related to SC, except for phase ICMs when using measures removing zero-lag coupling. The correlation between SC and ICMs increases with increasing frequency which is accompanied by reduced delays. Computational models produced results that were highly dependent on the specific parameter settings. The most consistent predictions were derived from measures solely based on SC. Overall, the results demonstrate that patterns of cortical functional coupling as reflected in both phase and envelope ICMs are both related, albeit to different degrees, to the underlying structural connectivity in the cerebral cortex.

Functional brain network properties of healthy full-term newborns quantified by scalp and source-reconstructed EEG

OBJECTIVE

Identifying the functional brain network properties of term low-risk newborns using high-density EEG (HD-EEG) and comparing these properties with those of established functional magnetic resonance image (fMRI) - based networks.

METHODS

HD-EEG was collected from 113 low-risk term newborns before delivery hospital discharge and within 72 hours of birth. Functional brain networks were reconstructed using coherence at the scalp and source levels in delta, theta, alpha, beta, and gamma frequency bands. These networks were characterized for the global and local network architecture.

RESULTS

Source-level networks in all the frequency bands identified the presence of the efficient small world (small-world propensity (SWP) > 0.6) architecture with four distinct modules linked by hub regions and rich-club (coefficient > 1) topology. The modular regions included primary, association, limbic, paralimbic, and subcortical regions, which have been demonstrated in fMRI studies. In contrast, scalp-level networks did not display consistent small world architecture (SWP < 0.6), and also identified only 2-3 modules in each frequency band.The modular regions of the scalp-network primarily included frontal and occipital regions.

CONCLUSIONS

Our findings show that EEG sources in low-risk newborns corroborate fMRI-based connectivity results.

SIGNIFICANCE

EEG source analysis characterizes functional connectivity at the bedside of low-risk newborn infants soon after birth.

Assessing the Repeatability of Multi-Frequency Multi-Layer Brain Network Topologies Across Alternative Researcher's Choice Paths

There is a growing interest in the neuroscience community on the advantages of multilayer functional brain networks. Researchers usually treated different frequencies separately at distinct functional brain networks. However, there is strong evidence that these networks share complementary information while their interdependencies could reveal novel findings. For this purpose, neuroscientists adopt multilayer networks, which can be described mathematically as an extension of trivial single-layer networks. Multilayer networks have become popular in neuroscience due to their advantage to integrate different sources of information. Here, Ι will focus on the multi-frequency multilayer functional connectivity analysis on resting-state fMRI (rs-fMRI) recordings. However, constructing a multilayer network depends on selecting multiple pre-processing steps that can affect the final network topology. Here, I analyzed the rs-fMRI dataset from a single human performing scanning over a period of 18 months (84 scans in total), and the rs-fMRI dataset containing 25 subjects with 3 repeat scans. I focused on assessing the reproducibility of multi-frequency multilayer topologies exploring the effect of two filtering methods for extracting frequencies from BOLD activity, three connectivity estimators, with or without a topological filtering scheme, and two spatial scales. Finally, I untangled specific combinations of researchers' choices that yield consistently brain networks with repeatable topologies, giving me the chance to recommend best practices over consistent topologies.

Multimodal multilayer network centrality relates to executive functioning

Executive functioning (EF) is a higher order cognitive process that is thought to depend on a network organization facilitating integration across subnetworks, in the context of which the central role of the fronto-parietal network (FPN) has been described across imaging and neurophysiological modalities. However, the potentially complementary unimodal information on the relevance of the FPN for EF has not yet been integrated. We employ a multilayer framework to allow for integration of different modalities into one 'network of networks.' We used diffusion MRI, resting-state functional MRI, MEG, and neuropsychological data obtained from 33 healthy adults to construct modality-specific single-layer networks as well as a single multilayer network per participant. We computed single-layer and multilayer eigenvector centrality of the FPN as a measure of integration in this network and examined their associations with EF. We found that higher multilayer FPN centrality, but not single-layer FPN centrality, was related to better EF. We did not find a statistically significant change in explained variance in EF when using the multilayer approach as compared to the single-layer measures. Overall, our results show the importance of FPN integration for EF and underline the promise of the multilayer framework toward better understanding cognitive functioning.

Topological reorganization of brain network might contribute to the resilience of cognitive functioning in mildly disabled relapsing remitting multiple sclerosis

Multiple sclerosis (MS) is an inflammatory and demyelinating disease which leads to impairment in several functional systems including cognition. Alteration of brain networks is linked to disability and its progression. However, results are mostly cross-sectional and yet contradictory as putative adaptive and maladaptive mechanisms were found. Here, we aimed to explore longitudinal reorganization of brain networks over 2-years by combining diffusion tensor imaging (DTI), resting-state functional MRI (fMRI), magnetoencephalography (MEG), and a comprehensive neuropsychological-battery. In 37 relapsing-remitting MS (RRMS) and 39 healthy-controls, cognition remained stable over-time. We reconstructed network models based on the three modalities and analyzed connectivity in relation to the hierarchical topology and functional subnetworks. Network models were compared across modalities and in their association with cognition using linear-mixed-effect-regression models. Loss of hub connectivity and global reduction was observed on a structural level over-years (p < .010), which was similar for functional MEG-networks but not for fMRI-networks. Structural hub connectivity increased in controls (p = .044), suggesting a physiological mechanism of healthy aging. Despite a general loss in structural connectivity in RRMS, hub connectivity was preserved (p = .002) over-time in default-mode-network (DMN). MEG-networks were similar to DTI and weakly correlated with fMRI in MS (p < .050). Lower structural (β between .23-.33) and both lower (β between .40-.59) and higher functional connectivity (β = -.54) in DMN was associated with poorer performance in attention and memory in RRMS (p < .001). MEG-networks involved no association with cognition. Here, cognitive stability despite ongoing neurodegeneration might indicate a resilience mechanism of DMN hubs mimicking a physiological reorganization observed in healthy aging.

Structural connectome constrained graphical lasso for MEG partial coherence

Structural connectivity provides the backbone for communication between neural populations. Since axonal transmission occurs on a millisecond time scale, measures of M/EEG functional connectivity sensitive to phase synchronization, such as coherence, are expected to reflect structural connectivity. We develop a model of MEG functional connectivity whose edges are constrained by the structural connectome. The edge strengths are defined by partial coherence, a measure of conditional dependence. We build a new method-the adaptive graphical lasso (AGL)-to fit the partial coherence to perform inference on the hypothesis that the structural connectome is reflected in MEG functional connectivity. In simulations, we demonstrate that the structural connectivity's influence on the partial coherence can be inferred using the AGL. Further, we show that fitting the partial coherence is superior to alternative methods at recovering the structural connectome, even after the source localization estimates required to map MEG from sensors to the cortex. Finally, we show how partial coherence can be used to explore how distinct parts of the structural connectome contribute to MEG functional connectivity in different frequency bands. Partial coherence offers better estimates of the strength of direct functional connections and consequently a potentially better estimate of network structure.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic–haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

What is the relationship between functional network architectures inferred from electromagnetic and haemodynamic data? This study shows that superposition of electromagnetic networks at canonical frequency bands manifests as highly structured patterns of haemodynamic functional connectivity in the human brain, reflecting systematic variation in cytoarchitecture.

Reproducibility of graph measures derived from resting-state MEG functional connectivity metrics in sensor and source spaces

Prior studies have used graph analysis of resting‐state magnetoencephalography (MEG) to characterize abnormal brain networks in neurological disorders. However, a present challenge for researchers is the lack of guidance on which network construction strategies to employ. The reproducibility of graph measures is important for their use as clinical biomarkers. Furthermore, global graph measures should ideally not depend on whether the analysis was performed in the sensor or source space. Therefore, MEG data of the 89 healthy subjects of the Human Connectome Project were used to investigate test–retest reliability and sensor versus source association of global graph measures. Atlas‐based beamforming was used for source reconstruction, and functional connectivity (FC) was estimated for both sensor and source signals in six frequency bands using the debiased weighted phase lag index (dwPLI), amplitude envelope correlation (AEC), and leakage‐corrected AEC. Reliability was examined over multiple network density levels achieved with proportional weight and orthogonal minimum spanning tree thresholding. At a 100% density, graph measures for most FC metrics and frequency bands had fair to excellent reliability and significant sensor versus source association. The greatest reliability and sensor versus source association was obtained when using amplitude metrics. Reliability was similar between sensor and source spaces when using amplitude metrics but greater for the source than the sensor space in higher frequency bands when using the dwPLI. These results suggest that graph measures are useful biomarkers, particularly for investigating functional networks based on amplitude synchrony.

Investigating the spectral features of the brain meso-scale structure at rest

Recent studies provide novel insights into the meso-scale organization of the brain, highlighting the co-occurrence of different structures: classic assortative (modular), disassortative, and core-periphery. However, the spectral properties of the brain meso-scale remain mostly unexplored. To fill this knowledge gap, we investigated how the meso-scale structure is organized across the frequency domain. We analyzed the resting state activity of healthy participants with source-localized high-density electroencephalography signals. Then, we inferred the community structure using weighted stochastic block-model (WSBM) to capture the landscape of meso-scale structures across the frequency domain. We found that different meso-scale modalities co-exist and are diversely organized over the frequency spectrum. Specifically, we found a core-periphery structure dominance, but we also highlighted a selective increase of disassortativity in the low frequency bands (<8 Hz), and of assortativity in the high frequency band (30-50 Hz). We further described other features of the meso-scale organization by identifying those brain regions which, at the same time, (a) exhibited the highest degree of assortativity, disassortativity, and core-peripheriness (i.e., participation) and (b) were consistently assigned to the same community, irrespective from the granularity imposed by WSBM (i.e., granularity-invariance). In conclusion, we observed that the brain spontaneous activity shows frequency-specific meso-scale organization, which may support spatially distributed and local information processing.

Effects of lactate and carbon monoxide interactions on neuroprotection and neuropreservation

Lactate, historically considered a waste product of anerobic metabolism, is a metabolite in whole-body metabolism needed for normal central nervous system (CNS) functions and a potent signaling molecule and hormone in the CNS. Neuronal activity signals normally induce its formation primarily in astrocytes and production is dependent on anerobic and aerobic metabolisms. Functions are dependent on normal dynamic, expansive, and evolving CNS functions. Levels can change under normal physiologic conditions and with CNS pathology. A readily combusted fuel that is sshuttled throughout the body, lactate is used as an energy source and is needed for CNS hemostasis, plasticity, memory, and excitability. Diffusion beyond the neuron active zone impacts activity of neurons and astrocytes in other areas of the brain. Barriergenesis, function of the blood-brain barrier, and buffering between oxidative metabolism and glycolysis and brain metabolism are affected by lactate. Important to neuroprotection, presence or absence is associated with L-lactate and heme oxygenase/carbon monoxide (a gasotransmitter) neuroprotective systems. Effects of carbon monoxide on L-lactate affect neuroprotection - interactions of the gasotransmitter with L-lactate are important to CNS stability, which will be reviewed in this article.

Static versus Functional PET: Making Sense of Metabolic Connectivity

Recently, Jamadar et al. (2021, Metabolic and hemodynamic resting-state connectivity of the human brain: a high-temporal resolution simultaneous BOLD-fMRI and FDG-fPET multimodality study. Cereb Cortex. 31(6), 2855-2867) compared the patterns of brain connectivity or covariance as obtained from 3 neuroimaging measures: 1) functional connectivity estimated from temporal correlations in the functional magnetic resonance imaging blood oxygen level-dependent signal, metabolic connectivity estimated, 2) from temporal correlations in 16-s frames of dynamic [18F]-fluorodeoxyglucose-positron emission tomography (FDG-PET), which they designate as functional FDG-PET (fPET), and 3) from intersubject correlations in static FDG-PET images (sPET). Here, we discuss a number of fundamental issues raised by the Jamadar study. These include the choice of terminology, the interpretation of cross-modal findings, the issue of group- to single-subject level inferences, and the meaning of metabolic connectivity as a biomarker. We applaud the methodological approach taken by the authors, but wish to present an alternative perspective on their findings. In particular, we argue that sPET and fPET can both provide valuable information about brain connectivity. Certainly, resolving this conundrum calls for further experimental and theoretical efforts to advance the developing framework of PET-based brain connectivity indices.

Exploring MEG brain fingerprints: Evaluation, pitfalls, and interpretations

Individual characterization of subjects based on their functional connectome (FC), termed "FC fingerprinting", has become a highly sought-after goal in contemporary neuroscience research. Recent functional magnetic resonance imaging (fMRI) studies have demonstrated unique characterization and accurate identification of individuals as an accomplished task. However, FC fingerprinting in magnetoencephalography (MEG) data is still widely unexplored. Here, we study resting-state MEG data from the Human Connectome Project to assess the MEG FC fingerprinting and its relationship with several factors including amplitude- and phase-coupling functional connectivity measures, spatial leakage correction, frequency bands, and behavioral significance. To this end, we first employ two identification scoring methods, differential identifiability and success rate, to provide quantitative fingerprint scores for each FC measurement. Secondly, we explore the edgewise and nodal MEG fingerprinting patterns across the different frequency bands (delta, theta, alpha, beta, and gamma). Finally, we investigate the cross-modality fingerprinting patterns obtained from MEG and fMRI recordings from the same subjects. We assess the behavioral significance of FC across connectivity measures and imaging modalities using partial least square correlation analyses. Our results suggest that fingerprinting performance is heavily dependent on the functional connectivity measure, frequency band, identification scoring method, and spatial leakage correction. We report higher MEG fingerprinting performances in phase-coupling methods, central frequency bands (alpha and beta), and in the visual, frontoparietal, dorsal-attention, and default-mode networks. Furthermore, cross-modality comparisons reveal a certain degree of spatial concordance in fingerprinting patterns between the MEG and fMRI data, especially in the visual system. Finally, the multivariate correlation analyses show that MEG connectomes have strong behavioral significance, which however depends on the considered connectivity measure and temporal scale. This comprehensive, albeit preliminary investigation of MEG connectome test-retest identifiability offers a first characterization of MEG fingerprinting in relation to different methodological and electrophysiological factors and contributes to the understanding of fingerprinting cross-modal relationships. We hope that this first investigation will contribute to setting the grounds for MEG connectome identification.

Interactions Between Aging and Alzheimer's Disease on Structural Brain Networks

Normative aging and Alzheimer's disease (AD) propagation alter anatomical connections among brain parcels. However, the interaction between the trajectories of age- and AD-linked alterations in the topology of the structural brain network is not well understood. In this study, diffusion-weighted magnetic resonance imaging (MRI) datasets of 139 subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database were used to document their structural brain networks. The 139 participants consist of 45 normal controls (NCs), 37 with early mild cognitive impairment (EMCI), 27 with late mild cognitive impairment (LMCI), and 30 AD patients. All subjects were further divided into three subgroups based on their age (56-65, 66-75, and 71-85 years). After the structural connectivity networks were built using anatomically-constrained deterministic tractography, their global and nodal topological properties were estimated, including network efficiency, characteristic path length, transitivity, modularity coefficient, clustering coefficient, and betweenness. Statistical analyses were then performed on these metrics using linear regression, and one- and two-way ANOVA testing to examine group differences and interactions between aging and AD propagation. No significant interactions were found between aging and AD propagation in the global topological metrics (network efficiency, characteristic path length, transitivity, and modularity coefficient). However, nodal metrics (clustering coefficient and betweenness centrality) of some cortical parcels exhibited significant interactions between aging and AD propagation, with affected parcels including left superior temporal, right pars triangularis, and right precentral. The results collectively confirm the age-related deterioration of structural networks in MCI and AD patients, providing novel insight into the cross effects of aging and AD disorder on brain structural networks. Some early symptoms of AD may also be due to age-associated anatomic vulnerability interacting with early anatomic changes associated with AD.

Predicting MEG resting-state functional connectivity from microstructural information

Understanding how human brain microstructure influences functional connectivity is an important endeavor. In this work, magnetic resonance imaging data from 90 healthy participants were used to calculate structural connectivity matrices using the streamline count, fractional anisotropy, radial diffusivity, and a myelin measure (derived from multicomponent relaxometry) to assign connection strength. Unweighted binarized structural connectivity matrices were also constructed. Magnetoencephalography resting-state data from those participants were used to calculate functional connectivity matrices, via correlations of the Hilbert envelopes of beamformer time series in the delta, theta, alpha, and beta frequency bands. Nonnegative matrix factorization was performed to identify the components of the functional connectivity. Shortest path length and search-information analyses of the structural connectomes were used to predict functional connectivity patterns for each participant. The microstructure-informed algorithms predicted the components of the functional connectivity more accurately than they predicted the total functional connectivity. This provides a methodology to understand functional mechanisms better. The shortest path length algorithm exhibited the highest prediction accuracy. Of the weights of the structural connectivity matrices, the streamline count and the myelin measure gave the most accurate predictions, while the fractional anisotropy performed poorly. Overall, different structural metrics paint very different pictures of the structural connectome and its relationship to functional connectivity.

Brain Connectivity Studies on Structure-Function Relationships: A Short Survey with an Emphasis on Machine Learning

This short survey reviews the recent literature on the relationship between the brain structure and its functional dynamics. Imaging techniques such as diffusion tensor imaging (DTI) make it possible to reconstruct axonal fiber tracks and describe the structural connectivity (SC) between brain regions. By measuring fluctuations in neuronal activity, functional magnetic resonance imaging (fMRI) provides insights into the dynamics within this structural network. One key for a better understanding of brain mechanisms is to investigate how these fast dynamics emerge on a relatively stable structural backbone. So far, computational simulations and methods from graph theory have been mainly used for modeling this relationship. Machine learning techniques have already been established in neuroimaging for identifying functionally independent brain networks and classifying pathological brain states. This survey focuses on methods from machine learning, which contribute to our understanding of functional interactions between brain regions and their relation to the underlying anatomical substrate.

Understanding neural flexibility from a multifaceted definition

Flexibility is a hallmark of human intelligence. Emerging studies have proposed several flexibility measurements at the level of individual regions, to produce a brain map of neural flexibility. However, flexibility is usually inferred from separate components of brain activity (i.e., intrinsic/task-evoked), and different definitions are used. Moreover, recent studies have argued that neural processing may be more than a task-driven and intrinsic dichotomy. Therefore, the understanding to neural flexibility is still incomplete. To address this issue, we propose a multifaceted definition of neural flexibility according to three key features: broad cognitive engagement, distributed connectivity, and adaptive connectome dynamics. For these three features, we first review the advances in computational approaches, their functional relevance, and their potential pitfalls. We then suggest a set of metrics that can help us assign a flexibility rating to each region. Subsequently, we present an emergent probabilistic view for further understanding the functional operation of individual regions in the unified framework of intrinsic and task-driven states. Finally, we highlight several areas related to the multifaceted definition of neural flexibility for future research. This review not only strengthens our understanding of flexible human brain, but also suggests that the measure of neural flexibility could bridge the gap between understanding intrinsic and task-driven brain function dynamics.

Multi-band MEG signatures of BOLD connectivity reorganization during visuospatial attention

The functional architecture of the resting brain, as measured with the blood oxygenation level-dependent functional connectivity (BOLD-FC), is slightly modified during task performance. In previous work, we reported behaviorally relevant BOLD-FC modulations between visual and dorsal attention regions when subjects performed a visuospatial attention task as compared to central fixation (Spadone et al., 2015). Here we use magnetoencephalography (MEG) in the same group of subjects to identify the electrophysiological correlates of the BOLD-FC modulation found in our previous work. While BOLD-FC topography, separately at rest and during visual attention, corresponded to neuromagnetic Band-Limited Power (BLP) correlation in the alpha and beta bands (8-30 Hz), BOLD-FC modulations evoked by performing the visual attention task (Spadone et al. 2015) did not match any specific oscillatory band BLP modulation. Conversely, following the application of an orthogonal spatial decomposition that identifies common inter-subject co-variations, we found that attention-rest BOLD-FC modulations were recapitulated by multi-spectral BLP-FC components. Notably, individual variability of alpha connectivity between Frontal Eye Fields and visual occipital regions, jointly with decreased interaction in the Visual network, correlated with visual discrimination accuracy. In summary, task-rest BOLD connectivity modulations match multi-spectral MEG BLP connectivity.

Integrated and segregated frequency architecture of the human brain network

The frequency of brain activity modulates the relationship between the brain and human behavior. Insufficient understanding of frequency-specific features may thus lead to inconsistent explanations of human behavior. However, to date, the frequency-specific features of the human brain functional network at the whole-brain level remain poorly understood. Here, we used resting-state fMRI data and graph-theory analyses to investigate the frequency-specific characteristics of fMRI signals in 12 frequency bands (frequency range 0.01-0.7 Hz) in 75 healthy participants. We found that brain regions with higher level and more complex functions had a more variable functional connectivity pattern but engaged less in higher frequency ranges. Moreover, brain regions that engaged in fewer frequency bands played more integrated roles (i.e., higher network participation coefficient and lower within-module degree) in the functional network, whereas regions that engaged in broader frequency ranges exhibited more segregated functions (i.e., lower network participation coefficient and higher within-module degree). Finally, behavioral analyses revealed that regional frequency variability was associated with a spectrum of behavioral functions from sensorimotor functions to complex cognitive and social functions. Taken together, our results showed that segregated functions are executed in wide frequency ranges, whereas integrated functions are executed mainly in lower frequency ranges. These frequency-specific features of brain networks provided crucial insights into the frequency mechanism of fMRI signals, suggesting that signals in higher frequency ranges should be considered for their relation to cognitive functions.

Brain dysconnectivity relates to disability and cognitive impairment in multiple sclerosis

The pathophysiology of cognitive dysfunction in multiple sclerosis (MS) is still unclear. This magnetoencephalography (MEG) study investigates the impact of MS on brain resting-state functional connectivity (rsFC) and its relationship to disability and cognitive impairment. We investigated rsFC based on power envelope correlation within and between different frequency bands, in a large cohort of participants consisting of 99 MS patients and 47 healthy subjects. Correlations were investigated between rsFC and outcomes on disability, disease duration and 7 neuropsychological scores within each group, while stringently correcting for multiple comparisons and possible confounding factors. Specific dysconnections correlating with MS-induced physical disability and disease duration were found within the sensorimotor and language networks, respectively. Global network-level reductions in within- and cross-network rsFC were observed in the default-mode network. Healthy subjects and patients significantly differed in their scores on cognitive fatigue and verbal fluency. Healthy subjects and patients showed different correlation patterns between rsFC and cognitive fatigue or verbal fluency, both of which involved a shift in patients from the posterior default-mode network to the language network. Introducing electrophysiological rsFC in a regression model of verbal fluency and cognitive fatigue in MS patients significantly increased the explained variance compared to a regression limited to structural MRI markers (relative thalamic volume and lesion load). This MEG study demonstrates that MS induces distinct changes in the resting-state functional brain architecture that relate to disability, disease duration and specific cognitive functioning alterations. It highlights the potential value of electrophysiological intrinsic rsFC for monitoring the cognitive impairment in patients with MS.

Tau pathology in early Alzheimer's disease is linked to selective disruptions in neurophysiological network dynamics

Understanding the role of Tau protein aggregation in the pathogenesis of Alzheimer's disease is critical for the development of new Tau-based therapeutic strategies to slow or prevent dementia. We tested the hypothesis that Tau pathology is associated with functional organization of widespread neurophysiological networks. We used electro-magnetoencephalography with [¹⁸F]AV-1451 PET scanning to quantify Tau-dependent network changes. Using a graph theoretical approach to brain connectivity, we quantified nodal measures of functional segregation, centrality, and the efficiency of information transfer and tested them against levels of [¹⁸F]AV-1451. Higher Tau burden in early Alzheimer's disease was associated with a shift away from the optimal small-world organization and a more fragmented network in the beta and gamma bands, whereby parieto-occipital areas were disconnected from the anterior parts of the network. Similarly, higher Tau burden was associated with decreases in both local and global efficiency, especially in the gamma band. The results support the translational development of neurophysiological "signatures" of Alzheimer's disease, to understand disease mechanisms in humans and facilitate experimental medicine studies.

Understanding the dynamics of biological and neural oscillator networks through exact mean-field reductions: a review

Many biological and neural systems can be seen as networks of interacting periodic processes. Importantly, their functionality, i.e., whether these networks can perform their function or not, depends on the emerging collective dynamics of the network. Synchrony of oscillations is one of the most prominent examples of such collective behavior and has been associated both with function and dysfunction. Understanding how network structure and interactions, as well as the microscopic properties of individual units, shape the emerging collective dynamics is critical to find factors that lead to malfunction. However, many biological systems such as the brain consist of a large number of dynamical units. Hence, their analysis has either relied on simplified heuristic models on a coarse scale, or the analysis comes at a huge computational cost. Here we review recently introduced approaches, known as the Ott-Antonsen and Watanabe-Strogatz reductions, allowing one to simplify the analysis by bridging small and large scales. Thus, reduced model equations are obtained that exactly describe the collective dynamics for each subpopulation in the oscillator network via few collective variables only. The resulting equations are next-generation models: Rather than being heuristic, they exactly link microscopic and macroscopic descriptions and therefore accurately capture microscopic properties of the underlying system. At the same time, they are sufficiently simple to analyze without great computational effort. In the last decade, these reduction methods have become instrumental in understanding how network structure and interactions shape the collective dynamics and the emergence of synchrony. We review this progress based on concrete examples and outline possible limitations. Finally, we discuss how linking the reduced models with experimental data can guide the way towards the development of new treatment approaches, for example, for neurological disease.

Intrinsic connectome organization across temporal scales: New insights from cross-modal approaches

The discovery of a stable, whole-brain functional connectivity organization that is largely independent of external events has drastically extended our view of human brain function. However, this discovery has been primarily based on functional magnetic resonance imaging (fMRI). The role of this whole-brain organization in fast oscillation-based connectivity as measured, for example, by electroencephalography (EEG) and magnetoencephalography (MEG) is only beginning to emerge. Here, we review studies of intrinsic connectivity and its whole-brain organization in EEG, MEG, and intracranial electrophysiology with a particular focus on direct comparisons to connectome studies in fMRI. Synthesizing this literature, we conclude that irrespective of temporal scale over four orders of magnitude, intrinsic neurophysiological connectivity shows spatial similarity to the connectivity organization commonly observed in fMRI. A shared structural connectivity basis and cross-frequency coupling are possible mechanisms contributing to this similarity. Acknowledging that a stable whole-brain organization governs long-range coupling across all timescales of neural processing motivates researchers to take "baseline" intrinsic connectivity into account when investigating brain-behavior associations, and further encourages more widespread exploration of functional connectomics approaches beyond fMRI by using EEG and MEG modalities.

Maladaptive alterations of resting state cortical network in Tinnitus: A directed functional connectivity analysis of a larger MEG data set

The present study used resting state MEG whole-head recordings to identify how chronic tonal tinnitus relates to altered functional connectivity of brain's intrinsic cortical networks. Resting state MEG activity of 40 chronic tinnitus patients and 40 matched human controls was compared identifying significant alterations in intrinsic networks of the tinnitus population. Directed functional connectivity of the resting brain, at a whole cortex level, was estimated by means of a statistical comparison of the estimated phase Transfer Entropy (pTE) between the time-series of cortical activations, as reconstructed by LORETA. As pTE identifies the direction of the information flow, a detailed analysis of the connectivity differences between tinnitus patients and controls was possible. Results indicate that the group of tinnitus patients show increased connectivity from right dorsal prefrontal to right medial temporal areas. Our results go beyond previous findings by indicating that the role of the left para-hippocampal area is dictated by a modulation from dmPFC; a region that is part of the dorsal attention network (DAN), as well as implicated in the regulation of emotional processing. Additionally, this whole cortex analysis showed a crucial role of the left inferior parietal cortex, which modulated the activity of the right superior temporal gyrus, providing new hypotheses for the role of this area within the context of current tinnitus models. Overall, these maladaptive alterations of the structure of intrinsic cortical networks show a decrease in efficiency and small worldness of the resting state network of tinnitus patients, which is correlated to tinnitus distress.

New Horizons in Early Dementia Diagnosis: Can Cerebellar Stimulation Untangle the Knot?

Differentiating Mild Cognitive Impairment (MCI) from dementia and estimating the risk of MCI-to-dementia conversion (MDC) are challenging tasks. Thus, objective tools are mandatory to get early diagnosis and prognosis. About that, there is a growing interest on the role of cerebellum-cerebrum connectivity (CCC). The aim of this study was to differentiate patients with an early diagnosis of dementia and MCI depending on the effects of a transcranial magnetic stimulation protocol (intermittent theta-burst stimulation -iTBS) delivered on the cerebellum able to modify cortico-cortical connectivity. Indeed, the risk of MDC is related to the response to iTBS, being higher in non-responder individuals. All patients with MCI, but eight (labelled as MCI-), showed preserved iTBS aftereffect. Contrariwise, none of the patients with dementia showed iTBS aftereffects. None of the patients showed EEG aftereffects following a sham TBS protocol. Five among the MCI- patients converted to dementia at 6-month follow-up. Our data suggest that cerebellar stimulation by means of iTBS may support the differential diagnosis between MCI and dementia and potentially identify the individuals with MCI who may be at risk of MDC. These findings may help clinicians to adopt a better prevention/follow-up strategy in such patients.

The homotopic connectivity of the functional brain: a meta-analytic approach

Homotopic connectivity (HC) is the connectivity between mirror areas of the brain hemispheres. It can exhibit a marked and functionally relevant spatial variability, and can be perturbed by several pathological conditions. The voxel-mirrored homotopic connectivity (VMHC) is a technique devised to enquire this pattern of brain organization, based on resting state functional connectivity. Since functional connectivity can be revealed also in a meta-analytical fashion using co-activations, here we propose to calculate the meta-analytic homotopic connectivity (MHC) as the meta-analytic counterpart of the VMHC. The comparison between the two techniques reveals their general similarity, but also highlights regional differences associated with how HC varies from task to rest. Two main differences were found from rest to task: (i) regions known to be characterized by global hubness are more similar than regions displaying local hubness; and (ii) medial areas are characterized by a higher degree of homotopic connectivity, while lateral areas appear to decrease their degree of homotopic connectivity during task performance. These findings show that MHC can be an insightful tool to study how the hemispheres functionally interact during task and rest conditions.

The Importance of the Validation of M/EEG With Current Biomarkers in Alzheimer's Disease

Current biomarkers used in research and in clinical practice in Alzheimer's Disease (AD) are the analysis of cerebral spinal fluid (CSF) to detect levels of Aβ42 and phosphorylated-tau, amyloid and FDG-PET, and MRI volumetry. Some of these procedures are still invasive for patients or expensive. Electroencephalography (EEG) and Magnetoencephalography (MEG) are two non-invasive techniques able to detect the early synaptic dysfunction and track the course of the disease. However, in spite of its added value they are not part of the standard of care in clinical practice in dementia. In this paper we review what these neurophysiological techniques can add to the early diagnosis of AD, whether results in both modalities are related to each other or not, as well as the need of its validation against current biomarkers. We discuss their potential implications for the better understanding of the pathophysiological mechanisms of the disease as well as the need of performing simultaneous M/EEG recordings to better understand discrepancies between these two techniques. Finally, more studies are needed studying M/EEG with amyloid and Tau biomarkers.

How do spatially distinct frequency specific MEG networks emerge from one underlying structural connectome? The role of the structural eigenmodes

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.

Linking structural and effective brain connectivity: structurally informed Parametric Empirical Bayes (si-PEB)

Despite the potential for better understanding functional neuroanatomy, the complex relationship between neuroimaging measures of brain structure and function has confounded integrative, multimodal analyses of brain connectivity. This is particularly true for task-related effective connectivity, which describes the causal influences between neuronal populations. Here, we assess whether measures of structural connectivity may usefully inform estimates of effective connectivity in larger scale brain networks. To this end, we introduce an integrative approach, capitalising on two recent statistical advances: Parametric Empirical Bayes, which provides group-level estimates of effective connectivity, and Bayesian model reduction, which enables rapid comparison of competing models. Crucially, we show that structural priors derived from high angular resolution diffusion imaging on a dynamic causal model of a 12-region network-based on functional MRI data from the same subjects-substantially improve model evidence (posterior probability 1.00). This provides definitive evidence that structural and effective connectivity depend upon each other in mediating distributed, large-scale interactions in the brain. Furthermore, this work offers novel perspectives for understanding normal brain architecture and its disintegration in clinical conditions.

Sensor Level Functional Connectivity Topography Comparison Between Different References Based EEG and MEG

Sensor-level functional connectivity topography (sFCT) contributes significantly to our understanding of brain networks. sFCT can be constructed using either electroencephalography (EEG) or magnetoencephalography (MEG). Here, we compared sFCT within the EEG modality and between EEG and MEG modalities. We first used simulations to look at how different EEG references-including the Reference Electrode Standardization Technique (REST), average reference (AR), linked mastoids (LM), and left mastoid references (LR)-affect EEG-based sFCT. The results showed that REST decreased the reference effects on scalp EEG recordings, making REST-based sFCT closer to the ground truth (sFCT based on ideal recordings). For the inter-modality simulation comparisons, we compared each type of EEG-sFCT with MEG-sFCT using three metrics to quantize the differences: Relative Error (RE), Overlap Rate (OR), and Hamming Distance (HD). When two sFCTs are similar, RE and HD are low, while OR is high. Results showed that among all reference schemes, EEG-and MEG-sFCT were most similar when the EEG was REST-based and the EEG and MEG were recorded simultaneously. Next, we analyzed simultaneously recorded MEG and EEG data from publicly available face-recognition experiments using a similar procedure as in the simulations. The results showed (1) if MEG-sFCT is the standard, REST-and LM-based sFCT provided results closer to this standard in the terms of HD; (2) REST-based sFCT and MEG-sFCT had the highest similarity in terms of RE; (3) REST-based sFCT had the most overlapping edges with MEG-sFCT in terms of OR. This study thus provides new insights into the effect of different reference schemes on sFCT and the similarity between MEG and EEG in terms of sFCT.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

SIGNIFICANCE

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

How Different EEG References Influence Sensor Level Functional Connectivity Graphs

Highlights: Hamming Distance is applied to distinguish the difference of functional connectivity networkThe orientations of sources are testified to influence the scalp Functional Connectivity Graph (FCG) from different references significantlyREST, the reference electrode standardization technique, is proved to have an overall stable and excellent performance in variable situations. The choice of an electroencephalograph (EEG) reference is a practical issue for the study of brain functional connectivity. To study how EEG reference influence functional connectivity estimation (FCE), this study compares the differences of FCE resulting from the different references such as REST (the reference electrode standardization technique), average reference (AR), linked mastoids (LM), and left mastoid references (LR). Simulations involve two parts. One is based on 300 dipolar pairs, which are located on the superficial cortex with a radial source direction. The other part is based on 20 dipolar pairs. In each pair, the dipoles have various orientation combinations. The relative error (RE) and Hamming distance (HD) between functional connectivity matrices of ideal recordings and that of recordings obtained with different references, are metrics to compare the differences of the scalp functional connectivity graph (FCG) derived from those two kinds of recordings. Lower RE and HD values imply more similarity between the two FCGs. Using the ideal recording (IR) as a standard, the results show that AR, LM and LR perform well only in specific conditions, i.e., AR performs stable when there is no upward component in sources' orientation. LR achieves desirable results when the sources' locations are away from left ear. LM achieves an indistinct difference with IR, i.e., when the distribution of source locations is symmetric along the line linking the two ears. However, REST not only achieves excellent performance for superficial and radial dipolar sources, but also achieves a stable and robust performance with variable source locations and orientations. Benefitting from the stable and robust performance of REST vs. other reference methods, REST might best recover the real FCG of EEG. Thus, REST based FCG may be a good candidate to compare the FCG of EEG based on different references from different labs.

The impact of MEG source reconstruction method on source-space connectivity estimation: A comparison between minimum-norm solution and beamforming

Despite numerous important contributions, the investigation of brain connectivity with magnetoencephalography (MEG) still faces multiple challenges. One critical aspect of source-level connectivity, largely overlooked in the literature, is the putative effect of the choice of the inverse method on the subsequent cortico-cortical coupling analysis. We set out to investigate the impact of three inverse methods on source coherence detection using simulated MEG data. To this end, thousands of randomly located pairs of sources were created. Several parameters were manipulated, including inter- and intra-source correlation strength, source size and spatial configuration. The simulated pairs of sources were then used to generate sensor-level MEG measurements at varying signal-to-noise ratios (SNR). Next, the source level power and coherence maps were calculated using three methods (a) L2-Minimum-Norm Estimate (MNE), (b) Linearly Constrained Minimum Variance (LCMV) beamforming, and (c) Dynamic Imaging of Coherent Sources (DICS) beamforming. The performances of the methods were evaluated using Receiver Operating Characteristic (ROC) curves. The results indicate that beamformers perform better than MNE for coherence reconstructions if the interacting cortical sources consist of point-like sources. On the other hand, MNE provides better connectivity estimation than beamformers, if the interacting sources are simulated as extended cortical patches, where each patch consists of dipoles with identical time series (high intra-patch coherence). However, the performance of the beamformers for interacting patches improves substantially if each patch of active cortex is simulated with only partly coherent time series (partial intra-patch coherence). These results demonstrate that the choice of the inverse method impacts the results of MEG source-space coherence analysis, and that the optimal choice of the inverse solution depends on the spatial and synchronization profile of the interacting cortical sources. The insights revealed here can guide method selection and help improve data interpretation regarding MEG connectivity estimation.

Source Localization of Brain States Associated with Canonical Neuroimaging Postures

Cognitive neuroscientists rarely consider the influence that body position exerts on brain activity; yet, postural variation holds important implications for the acquisition and interpretation of neuroimaging data. Whereas participants in most behavioral and EEG experiments sit upright, many prominent brain imaging techniques (e.g., fMRI) require participants to lie supine. Here we demonstrate that physical comportment profoundly alters baseline brain activity as measured by magnetoencephalography (MEG)-an imaging modality that permits multipostural acquisition. We collected resting-state MEG data from 12 healthy participants in three postures (lying supine, reclining at 45°, and sitting upright). Source-modeling analysis revealed a broadly distributed influence of posture on resting brain function. Sitting upright versus lying supine was associated with greater high-frequency (i.e., beta and gamma) activity in widespread parieto-occipital cortex. Moreover, sitting upright and reclining postures correlated with dampened activity in prefrontal regions across a range of bandwidths (i.e., from alpha to low gamma). The observed effects were large, with a mean Cohen's d of 0.95 ( SD = 0.23). In addition to neural activity, physiological parameters such as muscle tension and eye blinks may have contributed to these posture-dependent changes in brain signal. Regardless of the underlying mechanisms, however, the present results have important implications for the acquisition and interpretation of multimodal imaging data (e.g., studies combining fMRI or PET with EEG or MEG). More broadly, our findings indicate that generalizing results-from supine neuroimaging measurements to erect positions typical of ecological human behavior-would call for considering the influence that posture wields on brain dynamics.

Optimal phase synchronization in networks of phase-coherent chaotic oscillators

We investigate the existence of an optimal interplay between the natural frequencies of a group of chaotic oscillators and the topological properties of the network they are embedded in. We identify the conditions for achieving phase synchronization in the most effective way, i.e., with the lowest possible coupling strength. Specifically, we show by means of numerical and experimental results that it is possible to define a synchrony alignment function J(ω,L) linking the natural frequencies ω_i of a set of non-identical phase-coherent chaotic oscillators with the topology of the Laplacian matrix L, the latter accounting for the specific organization of the network of interactions between oscillators. We use the classical Rössler system to show that the synchrony alignment function obtained for phase oscillators can be extended to phase-coherent chaotic systems. Finally, we carry out a series of experiments with nonlinear electronic circuits to show the robustness of the theoretical predictions despite the intrinsic noise and parameter mismatch of the electronic components.

A Hitchhiker's Guide to Functional Magnetic Resonance Imaging

Functional Magnetic Resonance Imaging (fMRI) studies have become increasingly popular both with clinicians and researchers as they are capable of providing unique insights into brain functions. However, multiple technical considerations (ranging from specifics of paradigm design to imaging artifacts, complex protocol definition, and multitude of processing and methods of analysis, as well as intrinsic methodological limitations) must be considered and addressed in order to optimize fMRI analysis and to arrive at the most accurate and grounded interpretation of the data. In practice, the researcher/clinician must choose, from many available options, the most suitable software tool for each stage of the fMRI analysis pipeline. Herein we provide a straightforward guide designed to address, for each of the major stages, the techniques, and tools involved in the process. We have developed this guide both to help those new to the technique to overcome the most critical difficulties in its use, as well as to serve as a resource for the neuroimaging community.