BOLD fMRI Correlation Reflects Frequency-Specific Neuronal Correlation

Active suicidal ideation associated with dysfunction in default mode network using resting-state EEG and functional MRI - Findings from the T-RAD Study

Suicide in youth and young adults is a serious public health problem. However, the biological mechanisms of suicidal ideation (SI) remain poorly understood. The primary goal of these analyses was to identify the connectome profile of suicidal ideation using resting state electroencephalography (EEG). We evaluated the neurocircuitry of SI in a sample of youths and young adults (aged 10-26 years, n = 111) with current or past diagnoses of either a depressive disorder or bipolar disorder who were enrolled in the Texas Resilience Against Depression Study (T-RAD). Neurocircuitry was analyzed using orthogonalized power envelope connectivity computed from resting state EEG. Suicidal ideation was assessed with the 3-item Suicidal Thoughts factor of the Concise Health Risk Tracking self-report scale. The statistical pipeline involved dimension reduction using principal component analysis, and the association of neuroimaging data with SI using regularized canonical correlation analysis. From the original 111 participants and the correlation matrix of 4950 EEG connectivity pairs in each band (alpha, beta, theta), dimension reduction generated 1305 EEG connectivity pairs in the theta band, 2337 EEG pairs in the alpha band, and 914 EEG connectivity pairs in the beta band. Overall, SI was consistently involved with dysfunction of the default mode network (DMN). This report provides preliminary evidence of DMN dysfunction associated with active suicidal ideation in adolescents. Using EEG using power envelopes to compute connectivity moves us closer to using neurocircuit dysfunction in the clinical setting to identify suicidal ideation.

Spectral-switching analysis reveals real-time neuronal network representations of concurrent spontaneous naturalistic behaviors in human brain

Despite abundant evidence of functional networks in the human brain, their neuronal underpinnings, and relationships to real-time behavior have been challenging to resolve. Analyzing brain-wide intracranial-EEG recordings with video monitoring, acquired in awake subjects during clinical epilepsy evaluation, we discovered the tendency of each brain region to switch back and forth between 2 distinct power spectral densities (PSDs 2-55Hz). We further recognized that this 'spectral switching' occurs synchronously between distant sites, even between regions with differing baseline PSDs, revealing long-range functional networks that would be obscured in analysis of individual frequency bands. Moreover, the real-time PSD-switching dynamics of specific networks exhibited striking alignment with activities such as conversation and hand movements, revealing a multi-threaded functional network representation of concurrent naturalistic behaviors. Network structures and their relationships to behaviors were stable across days, but were altered during N3 sleep. Our results provide a new framework for understanding real-time, brain-wide neural-network dynamics.

Mobile neuroimaging: What we have learned about the neural control of human walking, with an emphasis on EEG-based research

Our understanding of the neural control of human walking has changed significantly over the last twenty years and mobile brain imaging methods have contributed substantially to current knowledge. High-density electroencephalography (EEG) has the advantages of being lightweight and mobile while providing temporal resolution of brain changes within a gait cycle. Advances in EEG hardware and processing methods have led to a proliferation of research on the neural control of locomotion in neurologically intact adults. We provide a narrative review of the advantages and disadvantages of different mobile brain imaging methods, then summarize findings from mobile EEG studies quantifying electrocortical activity during human walking. Contrary to historical views on the neural control of locomotion, recent studies highlight the widespread involvement of many areas, such as the anterior cingulate, posterior parietal, prefrontal, premotor, sensorimotor, supplementary motor, and occipital cortices, that show active fluctuations in electrical power during walking. The electrocortical activity changes with speed, stability, perturbations, and gait adaptation. We end with a discussion on the next steps in mobile EEG research.

Altered correlation of concurrently recorded EEG-fMRI connectomes in temporal lobe epilepsy

Whole-brain functional connectivity networks (connectomes) have been characterized at different scales in humans using EEG and fMRI. Multimodal epileptic networks have also been investigated, but the relationship between EEG and fMRI defined networks on a whole-brain scale is unclear. A unified multimodal connectome description, mapping healthy and pathological networks would close this knowledge gap. Here, we characterize the spatial correlation between the EEG and fMRI connectomes in right and left temporal lobe epilepsy (rTLE/lTLE). From two centers, we acquired resting-state concurrent EEG-fMRI of 35 healthy controls and 34 TLE patients. EEG-fMRI data was projected into the Desikan brain atlas, and functional connectomes from both modalities were correlated. EEG and fMRI connectomes were moderately correlated. This correlation was increased in rTLE when compared to controls for EEG-delta/theta/alpha/beta. Conversely, multimodal correlation in lTLE was decreased in respect to controls for EEG-beta. While the alteration was global in rTLE, in lTLE it was locally linked to the default mode network. The increased multimodal correlation in rTLE and decreased correlation in lTLE suggests a modality-specific lateralized differential reorganization in TLE, which needs to be considered when comparing results from different modalities. Each modality provides distinct information, highlighting the benefit of multimodal assessment in epilepsy.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and the underlying patterns of neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the amygdala of two male macaques. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5 Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-modal approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

Rhythmicity of neuronal oscillations delineates their cortical and spectral architecture

Neuronal oscillations are commonly analyzed with power spectral methods that quantify signal amplitude, but not rhythmicity or 'oscillatoriness' per se. Here we introduce a new approach, the phase-autocorrelation function (pACF), for the direct quantification of rhythmicity. We applied pACF to human intracerebral stereoelectroencephalography (SEEG) and magnetoencephalography (MEG) data and uncovered a spectrally and anatomically fine-grained cortical architecture in the rhythmicity of single- and multi-frequency neuronal oscillations. Evidencing the functional significance of rhythmicity, we found it to be a prerequisite for long-range synchronization in resting-state networks and to be dynamically modulated during event-related processing. We also extended the pACF approach to measure 'burstiness' of oscillatory processes and characterized regions with stable and bursty oscillations. These findings show that rhythmicity is double-dissociable from amplitude and constitutes a functionally relevant and dynamic characteristic of neuronal oscillations.

Data-driven MEG analysis to extract fMRI resting-state networks

The electrophysiological basis of resting-state networks (RSN) is still under debate. In particular, no principled mechanism has been determined that is capable of explaining all RSN equally well. While magnetoencephalography (MEG) and electroencephalography are the methods of choice to determine the electrophysiological basis of RSN, no standard analysis pipeline of RSN yet exists. In this article, we compare the two main existing data-driven analysis strategies for extracting RSNs from MEG data and introduce a third approach. The first approach uses phase-amplitude coupling to determine the RSN. The second approach extracts RSN through an independent component analysis of the Hilbert envelope in different frequency bands, while the third new approach uses a singular value decomposition instead. To evaluate these approaches, we compare the MEG-RSN to the functional magnetic resonance imaging (fMRI)-RSN from the same subjects. Overall, it was possible to extract RSN with MEG using all three techniques, which matched the group-specific fMRI-RSN. Interestingly the new approach based on SVD yielded significantly higher correspondence to five out of seven fMRI-RSN than the two existing approaches. Importantly, with this approach, all networks-except for the visual network-had the highest correspondence to the fMRI networks within one frequency band. Thereby we provide further insights into the electrophysiological underpinnings of the fMRI-RSNs. This knowledge will be important for the analysis of the electrophysiological connectome.

Cognitive abilities are associated with rapid dynamics of electrophysiological connectome states

Time-varying changes in whole-brain connectivity patterns, or connectome state dynamics, hold significant implications for cognition. However, connectome dynamics at fast (> 1Hz) timescales highly relevant to cognition are poorly understood due to the dominance of inherently slow fMRI in connectome studies. Here, we investigated the behavioral significance of rapid electrophysiological connectome dynamics using source-localized EEG connectomes during resting-state (N=926, 473 females). We focused on dynamic connectome features pertinent to individual differences, specifically those with established heritability: Fractional Occupancy (i.e., the overall duration spent in each recurrent connectome state) in beta and gamma bands, and Transition Probability (i.e., the frequency of state switches) in theta, alpha, beta, and gamma bands. Canonical correlation analysis found a significant relationship between the heritable phenotypes of sub-second connectome dynamics and cognition. Specifically, principal components of Transition Probabilities in alpha (followed by theta and gamma bands) and a cognitive factor representing visuospatial processing (followed by verbal and auditory working memory) most notably contributed to the relationship. We conclude that the specific order in which rapid connectome states are sequenced shapes individuals' cognitive abilities and traits. Such sub-second connectome dynamics may inform about behavioral function and dysfunction and serve as endophenotypes for cognitive abilities.

Rapid dynamics of electrophysiological connectome states are heritable

Time-varying changes in whole-brain connectivity patterns, or connectome state dynamics, are a prominent feature of brain activity with broad functional implications. While infra-slow (<0.1Hz) connectome dynamics have been extensively studied with fMRI, rapid dynamics highly relevant for cognition are poorly understood. Here, we asked whether rapid electrophysiological connectome dynamics constitute subject-specific brain traits and to what extent they are under genetic influence. Using source-localized EEG connectomes during resting-state (N=928, 473 females), we quantified heritability of multivariate (multi-state) features describing temporal or spatial characteristics of connectome dynamics. States switched rapidly every ~60-500ms. Temporal features were heritable, particularly, Fractional Occupancy (in theta, alpha, beta, and gamma bands) and Transition Probability (in theta, alpha, and gamma bands), representing the duration spent in each state and the frequency of state switches, respectively. Genetic effects explained a substantial proportion of phenotypic variance of these features: Fractional Occupancy in beta (44.3%) and gamma (39.8%) bands and Transition Probability in theta (38.4%), alpha (63.3%), beta (22.6%), and gamma (40%) bands. However, we found no evidence for heritability of spatial features, specifically states' Modularity and connectivity pattern. We conclude that genetic effects strongly shape individuals' connectome dynamics at rapid timescales, specifically states' overall occurrence and sequencing.

Abstract perceptual choice signals during action-linked decisions in the human brain

Humans can make abstract choices independent of motor actions. However, in laboratory tasks, choices are typically reported with an associated action. Consequentially, knowledge about the neural representation of abstract choices is sparse, and choices are often thought to evolve as motor intentions. Here, we show that in the human brain, perceptual choices are represented in an abstract, motor-independent manner, even when they are directly linked to an action. We measured MEG signals while participants made choices with known or unknown motor response mapping. Using multivariate decoding, we quantified stimulus, perceptual choice, and motor response information with distinct cortical distributions. Choice representations were invariant to whether the response mapping was known during stimulus presentation, and they occupied a distinct representational space from motor signals. As expected from an internal decision variable, they were informed by the stimuli, and their strength predicted decision confidence and accuracy. Our results demonstrate abstract neural choice signals that generalize to action-linked decisions, suggesting a general role of an abstract choice stage in human decision-making.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the macaque amygdala and activated them with a highly selective and potent DREADD agonist, deschloroclozapine. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Interestingly, activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-disciplinary approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

Gender may be related to the side of the motor syndrome and cognition in idiopathic Parkinson's disease

BACKGROUND

and Sex and cognitive profile may be related to the laterality of motor symptoms in idiopathic Parkinson's disease.

INTRODUCTION

Parkinson's disease (PD) is well recognised as an inherently asymmetric disease with unilateral onset of motor symptoms. The laterality of motor symptoms may be linked to sex, clinical and demographic variables, and neuropsychological disorders. However, the available data are inconsistent. This study aimed to explore the potential association between the laterality of motor symptoms and clinical and demographic variables and deficits in specific cognitive domains.

MATERIAL AND METHODS

We retrospectively recruited 97 participants with idiopathic PD without dementia; 60 presented motor symptoms on the left side and 37 on the right side. Both groups were comparable in terms of age, age at disease onset, disease duration, and severity of the neurological deficits according to the Unified Parkinson's Disease Rating Scale and the Hoehn and Yahr scale.

RESULTS

Participants with left-side motor symptoms scored lower on the Schwab and England Activities of Daily Living scale. Our sample included more men than women (67% vs. 33%). Both sexes were not equally represented in the 2 groups: there were significantly more men than women in the group of patients with left-side motor symptoms (77% vs. 23%), whereas the percentages of men and women in the group of patients with right-side motor symptoms were similar (51% vs. 49%). Both groups performed similarly in all neuropsychological tasks, but women, independently of laterality, performed better than men in the naming task.

CONCLUSION

We found a clear prevalence of men in the group of patients with left-side motor symptoms; this group also scored lower on the Schwab and England Scale. Female sex was predictive of better performance in the naming task. Sex should always be considered in disorders that cause asymmetric involvement of the brain, such as PD.

Mapping cross-species connectome atlas of human and macaque striatum

Cross-species connectome atlas (CCA) that can provide connectionally homogeneous and homologous brain nodes is essential and customized for cross-species neuroscience. However, existing CCAs were flawed in design and coarse-grained in results. In this study, a normative mapping framework of CCA was proposed and applied on human and macaque striatum. Specifically, all striatal voxels in the 2 species were mixed together and classified based on their represented and characterized feature of within-striatum resting-state functional connectivity, which was shared between the species. Six pairs of striatal parcels in these species were delineated in both hemispheres. Furthermore, this striatal parcellation was demonstrated by the best-matched whole-brain functional and structural connectivity between interspecies corresponding subregions. Besides, detailed interspecies differences in whole-brain multimodal connectivities and involved brain functions of these subregions were described to flesh out this CCA of striatum. In particular, this flexible and scalable mapping framework enables reliable construction of CCA of the whole brain, which would enable reliable findings in future cross-species research and advance our understandings into how the human brain works.

Dissociation between phase and power correlation networks in the human brain is driven by co-occurrent bursts

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.

The neural correlates of amplitude of low-frequency fluctuation: a multimodal resting-state MEG and fMRI-EEG study

The amplitude of low-frequency fluctuation (ALFF) describes the regional intensity of spontaneous blood-oxygen-level-dependent signal in resting-state functional magnetic resonance imaging (fMRI). How the fMRI-ALFF relates to the amplitude in electrophysiological signals remains unclear. We here aimed to investigate the neural correlates of fMRI-ALFF by comparing the spatial difference of amplitude between the eyes-closed (EC) and eyes-open (EO) states from fMRI and magnetoencephalography (MEG), respectively. By synthesizing MEG signal into amplitude-based envelope time course, we first investigated 2 types of amplitude in MEG, meaning the amplitude of neural activities from delta to gamma (i.e. MEG-amplitude) and the amplitude of their low-frequency modulation at the fMRI range (i.e. MEG-ALFF). We observed that the MEG-ALFF in EC was increased at parietal sensors, ranging from alpha to beta; whereas the MEG-amplitude in EC was increased at the occipital sensors in alpha. Source-level analysis revealed that the increased MEG-ALFF in the sensorimotor cortex overlapped with the most reliable EC-EO differences observed in fMRI at slow-3 (0.073-0.198 Hz), and these differences were more significant after global mean standardization. Taken together, our results support that (i) the amplitude at 2 timescales in MEG reflect distinct physiological information and that (ii) the fMRI-ALFF may relate to the ALFF in neural activity.

Multistage classification identifies altered cortical phase- and amplitude-coupling in Multiple Sclerosis

Distinguishing groups of subjects or experimental conditions in a high-dimensional feature space is a common goal in modern neuroimaging studies. Successful classification depends on the selection of relevant features as not every neuronal signal component or parameter is informative about the research question at hand. Here, we developed a novel unsupervised multistage analysis approach that combines dimensionality reduction, bootstrap aggregating and multivariate classification to select relevant neuronal features. We tested the approach by identifying changes of brain-wide electrophysiological coupling in Multiple Sclerosis. Multiple Sclerosis is a demyelinating disease of the central nervous system that can result in cognitive decline and physical disability. However, related changes in large-scale brain interactions remain poorly understood and corresponding non-invasive biomarkers are sparse. We thus compared brain-wide phase- and amplitude-coupling of frequency specific neuronal activity in relapsing-remitting Multiple Sclerosis patients (n = 17) and healthy controls (n = 17) using magnetoencephalography. Changes in this dataset included both, increased and decreased phase- and amplitude-coupling in wide-spread, bilateral neuronal networks across a broad range of frequencies. These changes allowed to successfully classify patients and controls with an accuracy of 84%. Furthermore, classification confidence predicted behavioral scores of disease severity. In sum, our results unravel systematic changes of large-scale phase- and amplitude coupling in Multiple Sclerosis. Furthermore, our results establish a new analysis approach to efficiently contrast high-dimensional neuroimaging data between experimental groups or conditions.

The relationship between electrophysiological and hemodynamic measures of neural activity varies across picture naming tasks: A multimodal magnetoencephalography-functional magnetic resonance imaging study

Different neuroimaging methods can yield different views of task-dependent neural engagement. Studies examining the relationship between electromagnetic and hemodynamic measures have revealed correlated patterns across brain regions but the role of the applied stimulation or experimental tasks in these correlation patterns is still poorly understood. Here, we evaluated the across-tasks variability of MEG-fMRI relationship using data recorded during three distinct naming tasks (naming objects and actions from action images, and objects from object images), from the same set of participants. Our results demonstrate that the MEG-fMRI correlation pattern varies according to the performed task, and that this variability shows distinct spectral profiles across brain regions. Notably, analysis of the MEG data alone did not reveal modulations across the examined tasks in the time-frequency windows emerging from the MEG-fMRI correlation analysis. Our results suggest that the electromagnetic-hemodynamic correlation could serve as a more sensitive proxy for task-dependent neural engagement in cognitive tasks than isolated within-modality measures.

Genetic polymorphisms in COMT and BDNF influence synchronization dynamics of human neuronal oscillations

Neuronal oscillations, their inter-areal synchronization, and scale-free dynamics constitute fundamental mechanisms for cognition by regulating communication in neuronal networks. These oscillatory dynamics have large inter-individual variability that is partly heritable. We hypothesized that this variability could be partially explained by genetic polymorphisms in neuromodulatory genes. We recorded resting-state magnetoencephalography (MEG) from 82 healthy participants and investigated whether oscillation dynamics were influenced by genetic polymorphisms in catechol-O-methyltransferase (COMT) Val158Met and brain-derived neurotrophic factor (BDNF) Val66Met. Both COMT and BDNF polymorphisms influenced local oscillation amplitudes and their long-range temporal correlations (LRTCs), while only BDNF polymorphism affected the strength of large-scale synchronization. Our findings demonstrate that COMT and BDNF genetic polymorphisms contribute to inter-individual variability in neuronal oscillation dynamics. Comparison of these results to computational modeling of near-critical synchronization dynamics further suggested that COMT and BDNF polymorphisms influenced local oscillations by modulating the excitation-inhibition balance according to the brain criticality framework.

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.

Resting state EEG power spectrum and functional connectivity in autism: a cross-sectional analysis

BACKGROUND

Understanding the development of the neuronal circuitry underlying autism spectrum disorder (ASD) is critical to shed light into its etiology and for the development of treatment options. Resting state EEG provides a window into spontaneous local and long-range neuronal synchronization and has been investigated in many ASD studies, but results are inconsistent. Unbiased investigation in large and comprehensive samples focusing on replicability is needed.

METHODS

We quantified resting state EEG alpha peak metrics, power spectrum (PS, 2-32 Hz) and functional connectivity (FC) in 411 children, adolescents and adults (n = 212 ASD, n = 199 neurotypicals [NT], all with IQ > 75). We performed analyses in source-space using individual head models derived from the participants' MRIs. We tested for differences in mean and variance between the ASD and NT groups for both PS and FC using linear mixed effects models accounting for age, sex, IQ and site effects. Then, we used machine learning to assess whether a multivariate combination of EEG features could better separate ASD and NT participants. All analyses were embedded within a train-validation approach (70%-30% split).

RESULTS

In the training dataset, we found an interaction between age and group for the reactivity to eye opening (p = .042 uncorrected), and a significant but weak multivariate ASD vs. NT classification performance for PS and FC (sensitivity 0.52-0.62, specificity 0.59-0.73). None of these findings replicated significantly in the validation dataset, although the effect size in the validation dataset overlapped with the prediction interval from the training dataset.

LIMITATIONS

The statistical power to detect weak effects-of the magnitude of those found in the training dataset-in the validation dataset is small, and we cannot fully conclude on the reproducibility of the training dataset's effects.

CONCLUSIONS

This suggests that PS and FC values in ASD and NT have a strong overlap, and that differences between both groups (in both mean and variance) have, at best, a small effect size. Larger studies would be needed to investigate and replicate such potential effects.

Spectral Fingerprints of Cortical Neuromodulation

Pupil size has been established as a versatile marker of noradrenergic and cholinergic neuromodulation, which has profound effects on neuronal processing, cognition, and behavior. However, little is known about the cortical control and effects of pupil-linked neuromodulation. Here, we show that pupil dynamics are tightly coupled to temporally, spectrally, and spatially specific modulations of local and large-scale cortical population activity in the human brain. We quantified the dynamics of band-limited cortical population activity in resting human subjects using magnetoencephalography and investigated how neural dynamics were linked to simultaneously recorded pupil dynamics. Our results show that pupil-linked neuromodulation does not merely affect cortical population activity in a stereotypical fashion. Instead, we identified three frontal, precentral, and occipitoparietal networks, in which local population activity with distinct spectral profiles in the theta, beta, and alpha bands temporally preceded and followed changes in pupil size. Furthermore, we found that amplitude coupling at ∼16 Hz in a large-scale frontoparietal network predicted pupil dynamics. Our results unravel network-specific spectral fingerprints of cortical neuromodulation in the human brain that likely reflect both the causes and effects of neuromodulation.SIGNIFICANCE STATEMENT Brain function is constantly affected by modulatory neurotransmitters. Pupil size has been established as a versatile marker of noradrenergic and cholinergic neuromodulation. However, because the cortical correlates of pupil dynamics are largely unknown, fundamental questions remain unresolved. Which cortical networks control pupil-linked neuromodulation? Does neuromodulation affect cortical activity in a stereotypical or region-specific fashion? To address this, we quantified the dynamics of cortical population activity in human subjects using magnetoencephalography. We found that pupil dynamics are coupled to highly specific modulations of local and large-scale cortical activity in the human brain. We identified four cortical networks with distinct spectral profiles that temporally predicted and followed pupil size dynamics. These effects likely reflect both the cortical control and effect of neuromodulation.

Functional Connectivity of the Brain Across Rodents and Humans

Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.

Harmoni: a Method for Eliminating Spurious Interactions due to the Harmonic Components in Neuronal Data

Cross-frequency synchronization (CFS) has been proposed as a mechanism for integrating spatially and spectrally distributed information in the brain. However, investigating CFS in Magneto- and Electroencephalography (MEG/EEG) is hampered by the presence of spurious neuronal interactions due to the non-sinusoidal waveshape of brain oscillations. Such waveshape gives rise to the presence of oscillatory harmonics mimicking genuine neuronal oscillations. Until recently, however, there has been no methodology for removing these harmonics from neuronal data. In order to address this long-standing challenge, we introduce a novel method (called HARMOnic miNImization - Harmoni) that removes the signal components which can be harmonics of a non-sinusoidal signal. Harmoni's working principle is based on the presence of CFS between harmonic components and the fundamental component of a non-sinusoidal signal. We extensively tested Harmoni in realistic EEG simulations. The simulated couplings between the source signals represented genuine and spurious CFS and within-frequency phase synchronization. Using diverse evaluation criteria, including ROC analyses, we showed that the within- and cross-frequency spurious interactions are suppressed significantly, while the genuine activities are not affected. Additionally, we applied Harmoni to real resting-state EEG data revealing intricate remote connectivity patterns which are usually masked by the spurious connections. Given the ubiquity of non-sinusoidal neuronal oscillations in electrophysiological recordings, Harmoni is expected to facilitate novel insights into genuine neuronal interactions in various research fields, and can also serve as a steppingstone towards the development of further signal processing methods aiming at refining within- and cross-frequency synchronization in electrophysiological recordings.

Connectomics of human electrophysiology

We present both a scientific overview and conceptual positions concerning the challenges and assets of electrophysiological measurements in the search for the nature and functions of the human connectome. We discuss how the field has been inspired by findings and approaches from functional magnetic resonance imaging (fMRI) and informed by a small number of significant multimodal empirical studies, which show that the canonical networks that are commonplace in fMRI are in fact rooted in electrophysiological processes. This review is also an opportunity to produce a brief, up-to-date critical survey of current data modalities and analytical methods available for deriving both static and dynamic connectomes from electrophysiology. We review hurdles that challenge the significance and impact of current electrophysiology connectome research. We then encourage the field to take a leap of faith and embrace the wealth of electrophysiological signals, despite their apparent, disconcerting complexity. Our position is that electrophysiology connectomics is poised to inform testable mechanistic models of information integration in hierarchical brain networks, constructed from observable oscillatory and aperiodic signal components and their polyrhythmic interactions.

Brain functional and effective connectivity based on electroencephalography recordings: A review

Functional connectivity and effective connectivity of the human brain, representing statistical dependence and directed information flow between cortical regions, significantly contribute to the study of the intrinsic brain network and its functional mechanism. Many recent studies on electroencephalography (EEG) have been focusing on modeling and estimating brain connectivity due to increasing evidence that it can help better understand various brain neurological conditions. However, there is a lack of a comprehensive updated review on studies of EEG-based brain connectivity, particularly on visualization options and associated machine learning applications, aiming to translate those techniques into useful clinical tools. This article reviews EEG-based functional and effective connectivity studies undertaken over the last few years, in terms of estimation, visualization, and applications associated with machine learning classifiers. Methods are explored and discussed from various dimensions, such as either linear or nonlinear, parametric or nonparametric, time-based, and frequency-based or time-frequency-based. Then it is followed by a novel review of brain connectivity visualization methods, grouped by Heat Map, data statistics, and Head Map, aiming to explore the variation of connectivity across different brain regions. Finally, the current challenges of related research and a roadmap for future related research are presented.

Dynamic primitives of brain network interaction Special Issue "Advances in Mapping the Connectome"

What dynamic processes underly functional brain networks? Functional connectivity (FC) and functional connectivity dynamics (FCD) are used to represent the patterns and dynamics of functional brain networks. FC(D) is related to the synchrony of brain activity: when brain areas oscillate in a coordinated manner this yields a high correlation between their signal time series. To explain the processes underlying FC(D) we review how synchronized oscillations emerge from coupled neural populations in brain network models (BNMs). From detailed spiking networks to more abstract population models, there is strong support for the idea that the brain operates near critical instabilities that give rise to multistable or metastable dynamics that in turn lead to the intermittently synchronized slow oscillations underlying FC(D). We explore further consequences from these fundamental mechanisms and how they fit with reality. We conclude by highlighting the need for integrative brain models that connect separate mechanisms across levels of description and spatiotemporal scales and link them with cognitive function.

Deficient Interhemispheric Connectivity Underlies Movement Irregularities in Parkinson's Disease

BACKGROUND

Movement execution is impaired in patients with Parkinson's disease. Evolving neurodegeneration leads to altered connectivity between distinct regions of the brain and altered activity at interconnected areas. How connectivity alterations influence complex movements like drawing spirals in Parkinson's disease patients remains largely unexplored.

OBJECTIVE

We investigated whether deteriorations in interregional connectivity relate to impaired execution of drawing.

METHODS

Twenty-nine patients and 31 age-matched healthy control participants drew spirals with both hands on a digital graphics tablet, and the regularity of drawing execution was evaluated by sample entropy. We recorded resting-state fMRI and task-related EEG, and calculated the time-resolved partial directed coherence to estimate effective connectivity for both imaging modalities to determine the extent and directionality of interregional interactions.

RESULTS

Movement performance in Parkinson's disease patients was characterized by increased sample entropy, corresponding to enhanced irregularities in task execution. Effective connectivity between the motor cortices of both hemispheres, derived from resting-state fMRI, was significantly reduced in Parkinson's disease patients in comparison to controls. The connectivity strength in the nondominant to dominant hemisphere direction in both modalities was inversely correlated with irregularities during drawing, but not with the clinical state.

CONCLUSION

Our findings suggest that interhemispheric connections are affected both at rest and during drawing movements by Parkinson's disease. This provides novel evidence that disruptions of interhemispheric information exchange play a pivotal role for impairments of complex movement execution in Parkinson's disease patients.

Frequency-Specific Blood Oxygen Level Dependent Oscillations Associated With Pain Relief From Ankle Acupuncture in Patients With Chronic Low Back Pain

Objective: Recent advances in brain imaging have deepened our knowledge of the neural activity in distinct brain areas associated with acupuncture analgesia. However, there has not been conclusive research into the frequency-specific resting-state functional changes associated with acupuncture analgesia in patients with chronic pain. Here, we aimed to characterize changes across multiple frequencies of resting-state cortical activity associated with ankle acupuncture stimulation (AAS) in patients with chronic low back pain (CLBP) and healthy controls.

Methods: Twenty seven patients with CLBP and Twenty five age- and gender-matched healthy volunteers were enrolled in the study. Participants received tactile sham acupuncture (TSA) and AAS, respectively. The whole-brain amplitude of low-frequency fluctuation (ALFF) in the range 0.01–0.25 Hz was assessed for changes associated with each intervention. Further, a visual analog scale (VAS) was used to collect subjective measures of pain intensity in patients. Linear mixed-effect modeling (LME) was used to examine the mean ALFF values of AAS and TSA between patients and healthy controls.

Results: The ALFF was modulated in the default mode network (an increase in the medial prefrontal cortex, and a decrease in the cerebellum/posterior ingulate/parahippocampus, P < 0.01, corrected) in both patients and controls. Decreased ALFF in the bilateral insular was frequency-dependent. Modulations in the cerebellum and right insular were significantly correlated with VAS pain score after AAS (P < 0.01).

Conclusion: Hence, frequency-specific resting-state activity in the cerebellum and insular was correlated to AAS analgesia. Our frequency-specific analysis of ALFF may provide novel insights related to pain relief from acupuncture.

Cortical correlation structure of aperiodic neuronal population activity

Electrophysiological population signals contain oscillatory and non-oscillatory aperiodic (1/frequency-like) components. So far research has largely focused on oscillatory activity, and only recently, interest in aperiodic population activity has gained momentum. Accordingly, while the cortical correlation structure of oscillatory population activity has been characterized, little is known about the correlation of aperiodic neuronal activity. To address this, we investigated aperiodic neuronal population activity in the human brain using resting-state magnetoencephalography (MEG). We combined source-analysis, signal orthogonalization and irregular-resampling auto-spectral analysis (IRASA) to systematically characterize the cortical distribution and correlation of aperiodic neuronal activity. We found that aperiodic population activity is robustly correlated across the cortex and that this correlation is spatially well structured. Furthermore, we found that the cortical correlation structure of aperiodic activity is similar but distinct from the correlation structure of oscillatory neuronal activity. Anterior cortical regions showed the strongest differences between oscillatory and aperiodic correlation patterns. Our results suggest that correlations of aperiodic population activity serve as robust markers of cortical network interactions. Furthermore, our results show that aperiodic and oscillatory signal components provide non-redundant information about large-scale neuronal correlations. This may reflect at least partly distinct neuronal mechanisms underlying and reflected by oscillatory and aperiodic neuronal population activity.

Resting state network connectivity is attenuated by fMRI acoustic noise

INTRODUCTION

During the past decades there has been an increasing interest in tracking brain network fluctuations in health and disease by means of resting state functional magnetic resonance imaging (rs-fMRI). Rs-fMRI however does not provide the ideal environmental setting, as participants are continuously exposed to noise generated by MRI coils during acquisition of Echo Planar Imaging (EPI). We investigated the effect of EPI noise on resting state activity and connectivity using magnetoencephalography (MEG), by reproducing the acoustic characteristics of rs-fMRI environment during the recordings. As compared to fMRI, MEG has little sensitivity to brain activity generated in deep brain structures, but has the advantage to capture both the dynamic of cortical magnetic oscillations with high temporal resolution and the slow magnetic fluctuations highly correlated with BOLD signal.

METHODS

Thirty healthy subjects were enrolled in a counterbalanced design study including three conditions: a) silent resting state (Silence), b) resting state upon EPI noise (fMRI), and c) resting state upon white noise (White). White noise was employed to test the specificity of fMRI noise effect. The amplitude envelope correlation (AEC) in alpha band measured the connectivity of seven Resting State Networks (RSN) of interest (default mode network, dorsal attention network, language, left and right auditory and left and right sensory-motor). Vigilance dynamic was estimated from power spectral activity.

RESULTS

fMRI and White acoustic noise consistently reduced connectivity of cortical networks. The effects were widespread, but noise and network specificities were also present. For fMRI noise, decreased connectivity was found in the right auditory and sensory-motor networks. Progressive increase of slow theta-delta activity related to drowsiness was found in all conditions, but was significantly higher for fMRI . Theta-delta significantly and positively correlated with variations of cortical connectivity.

DISCUSSION

rs-fMRI connectivity is biased by unavoidable environmental factors during scanning, which warrant more careful control and improved experimental designs. MEG is free from acoustic noise and allows a sensitive estimation of resting state connectivity in cortical areas. Although underutilized, MEG could overcome issues related to noise during fMRI, in particular when investigation of motor and auditory networks is needed.

Brain rhythms define distinct interaction networks with differential dependence on anatomy

Cognitive functions are subserved by rhythmic neuronal synchronization across widely distributed brain areas. In 105 area pairs, we investigated functional connectivity (FC) through coherence, power correlation, and Granger causality (GC) in the theta, beta, high-beta, and gamma rhythms. Between rhythms, spatial FC patterns were largely independent. Thus, the rhythms defined distinct interaction networks. Importantly, networks of coherence and GC were not explained by the spatial distributions of the strengths of the rhythms. Those networks, particularly the GC networks, contained clear modules, with typically one dominant rhythm per module. To understand how this distinctiveness and modularity arises on a common anatomical backbone, we correlated, across 91 area pairs, the metrics of functional interaction with those of anatomical projection strength. Anatomy was primarily related to coherence and GC, with the largest effect sizes for GC. The correlation differed markedly between rhythms, being less pronounced for the beta and strongest for the gamma rhythm.

Dynamic large-scale connectivity of intrinsic cortical oscillations supports adaptive listening in challenging conditions

In multi-talker situations, individuals adapt behaviorally to this listening challenge mostly with ease, but how do brain neural networks shape this adaptation? We here establish a long-sought link between large-scale neural communications in electrophysiology and behavioral success in the control of attention in difficult listening situations. In an age-varying sample of N = 154 individuals, we find that connectivity between intrinsic neural oscillations extracted from source-reconstructed electroencephalography is regulated according to the listener's goal during a challenging dual-talker task. These dynamics occur as spatially organized modulations in power-envelope correlations of alpha and low-beta neural oscillations during approximately 2-s intervals most critical for listening behavior relative to resting-state baseline. First, left frontoparietal low-beta connectivity (16 to 24 Hz) increased during anticipation and processing of a spatial-attention cue before speech presentation. Second, posterior alpha connectivity (7 to 11 Hz) decreased during comprehension of competing speech, particularly around target-word presentation. Connectivity dynamics of these networks were predictive of individual differences in the speed and accuracy of target-word identification, respectively, but proved unconfounded by changes in neural oscillatory activity strength. Successful adaptation to a listening challenge thus latches onto two distinct yet complementary neural systems: a beta-tuned frontoparietal network enabling the flexible adaptation to attentive listening state and an alpha-tuned posterior network supporting attention to speech.

Principles and open questions in functional brain network reconstruction

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network representation involves often covert theoretical assumptions and methodological choices which affect the way networks are reconstructed from experimental data, and ultimately the resulting network properties and their interpretation. Here, we review some fundamental conceptual underpinnings and technical issues associated with brain network reconstruction, and discuss how their mutual influence concurs in clarifying the organization of brain function.

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.

Impaired intrinsic functional connectivity among medial temporal lobe and sub-regions related to memory deficits in intracranial dural arteriovenous fistula

PURPOSE

The functional changes concerning memory deficits in dural arteriovenous fistula (dAVF) brain are inadequately understood. This study aimed to understand the functional connectivity alterations of brain regions widely affirmed for explicit and implicit memory functions in dAVF patients (DP) and look into the frequency effects of the altered functional networks.

METHODS

Resting-state functional magnetic resonance imaging (rsfMRI) analysis was done in the memory-associated regions of 30 DP and 30 healthy controls (HC). Frequency decomposition was used to determine potential frequency-dependent functional connectivity changes. They underwent neuropsychological tests and were correlated with changes in memory networks compared with HC.

RESULTS

The results showed weaker functional connectivity among the medial temporal lobe and sub-regions in DP suggestive of dysfunction of explicit and implicit memory functions, which corroborated with the positive correlation between memory scores and hippocampal-parahippocampal connectivity of DP, along with a significant group difference of lower memory and cognitive performance in DP assessed by neuropsychological tests. A frequency-dependent study of the altered rsFC revealed lower functional connectivity strength and impaired neural coupling manifested at some sub-band frequencies indicative of disturbed cortical rhythm in DP.

CONCLUSION

This pilot study gives insights into significant intrinsic functional connectivity changes in the memory regions of the dAVF brain. The results may have clinical implications in the choice of interventional management of dAVF and can impact clinical decision making for realizable prevention of progressive memory impairment and irreversible brain damage in such patients.

Behavioural relevance of spontaneous, transient brain network interactions in fMRI

How spontaneously fluctuating functional magnetic resonance imaging (fMRI) signals in different brain regions relate to behaviour has been an open question for decades. Correlations in these signals, known as functional connectivity, can be averaged over several minutes of data to provide a stable representation of the functional network architecture for an individual. However, associations between these stable features and behavioural traits have been shown to be dominated by individual differences in anatomy. Here, using kernel learning tools, we propose methods to assess and compare the relation between time-varying functional connectivity, time-averaged functional connectivity, structural brain data, and non-imaging subject behavioural traits. We applied these methods to Human Connectome Project resting-state fMRI data to show that time-varying fMRI functional connectivity, detected at time-scales of a few seconds, has associations with some behavioural traits that are not dominated by anatomy. Despite time-averaged functional connectivity accounting for the largest proportion of variability in the fMRI signal between individuals, we found that some aspects of intelligence could only be explained by time-varying functional connectivity. The finding that time-varying fMRI functional connectivity has a unique relationship to population behavioural variability suggests that it might reflect transient neuronal communication fluctuating around a stable neural architecture.

Gender may be related to the side of the motor syndrome and cognition in idiopathic Parkinson's disease

BACKGROUND

and Sex and cognitive profile may be related to the laterality of motor symptoms in idiopathic Parkinson's disease.

INTRODUCTION

Parkinson's disease (PD) is well recognised as an inherently asymmetric disease with unilateral onset of motor symptoms. The laterality of motor symptoms may be linked to sex, clinical and demographic variables, and neuropsychological disorders. However, the available data are inconsistent. This study aimed to explore the potential association between the laterality of motor symptoms and clinical and demographic variables and deficits in specific cognitive domains.

MATERIAL AND METHODS

We retrospectively recruited 97 participants with idiopathic PD without dementia; 60 presented motor symptoms on the left side and 37 on the right side. Both groups were comparable in terms of age, age at disease onset, disease duration, and severity of the neurological deficits according to the Unified Parkinson's Disease Rating Scale and the Hoehn and Yahr scale.

RESULTS

Participants with left-side motor symptoms scored lower on the Schwab and England Activities of Daily Living scale. Our sample included more men than women (67% vs. 33%). Both sexes were not equally represented in the 2 groups: there were significantly more men than women in the group of patients with left-side motor symptoms (77% vs. 23%), whereas the percentages of men and women in the group of patients with right-side motor symptoms were similar (51% vs. 49%). Both groups performed similarly in all neuropsychological tasks, but women, independently of laterality, performed better than men in the naming task.

CONCLUSION

We found a clear prevalence of men in the group of patients with left-side motor symptoms; this group also scored lower on the Schwab and England Scale. Female sex was predictive of better performance in the naming task. Sex should always be considered in disorders that cause asymmetric involvement of the brain, such as PD.

The relationship between EEG and fMRI connectomes is reproducible across simultaneous EEG-fMRI studies from 1.5T to 7T

Both electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are non-invasive methods that show complementary aspects of human brain activity. Despite measuring different proxies of brain activity, both the measured blood-oxygenation (fMRI) and neurophysiological recordings (EEG) are indirectly coupled. The electrophysiological and BOLD signal can map the underlying functional connectivity structure at the whole brain scale at different timescales. Previous work demonstrated a moderate but significant correlation between resting-state functional connectivity of both modalities, however there is a wide range of technical setups to measure simultaneous EEG-fMRI and the reliability of those measures between different setups remains unknown. This is true notably with respect to different magnetic field strengths (low and high field) and different spatial sampling of EEG (medium to high-density electrode coverage). Here, we investigated the reproducibility of the bimodal EEG-fMRI functional connectome in the most comprehensive resting-state simultaneous EEG-fMRI dataset compiled to date including a total of 72 subjects from four different imaging centers. Data was acquired from 1.5T, 3T and 7T scanners with simultaneously recorded EEG using 64 or 256 electrodes. We demonstrate that the whole-brain monomodal connectivity reproducibly correlates across different datasets and that a moderate crossmodal correlation between EEG and fMRI connectivity of r ≈ 0.3 can be reproducibly extracted in low- and high-field scanners. The crossmodal correlation was strongest in the EEG-β frequency band but exists across all frequency bands. Both homotopic and within intrinsic connectivity network (ICN) connections contributed the most to the crossmodal relationship. This study confirms, using a considerably diverse range of recording setups, that simultaneous EEG-fMRI offers a consistent estimate of multimodal functional connectomes in healthy subjects that are dominantly linked through a functional core of ICNs across spanning across the different timescales measured by EEG and fMRI. This opens new avenues for estimating the dynamics of brain function and provides a better understanding of interactions between EEG and fMRI measures. This observed level of reproducibility also defines a baseline for the study of alterations of this coupling in pathological conditions and their role as potential clinical markers.

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.

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

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.

Dissociating harmonic and non-harmonic phase-amplitude coupling in the human brain

Phase-amplitude coupling (PAC) has been hypothesized to coordinate cross-frequency interactions of neuronal activity in the brain. However, little is known about the distribution of PAC across the human brain and the frequencies involved. Furthermore, it remains unclear to what extent PAC may reflect spurious cross-frequency coupling induced by physiological artifacts or rhythmic non-sinusoidal signals with higher harmonics. Here, we combined MEG, source-reconstruction and different measures of cross-frequency coupling to systematically characterize local PAC across the resting human brain. We show that cross-frequency measures of phase-amplitude, phase-phase, and amplitude-amplitude coupling are all sensitive to signals with higher harmonics. In conjunction, these measures allow to distinguish harmonic and non-harmonic PAC. Based on these insights, we found no evidence for non-harmonic local PAC in resting-state MEG. Instead, we found cortically and spectrally wide-spread PAC driven by harmonic signals. Furthermore, we show how physiological artifacts and spectral leakage cause spurious PAC across wide frequency ranges. Our results clarify how different measures of cross-frequency interactions can be combined to characterize PAC, and cast doubt on the presence of prominent non-harmonic phase-amplitude coupling in human resting-state MEG.

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.

Oscillation-Based Connectivity Architecture Is Dominated by an Intrinsic Spatial Organization, Not Cognitive State or Frequency

Functional connectivity of neural oscillations (oscillation-based FC) is thought to afford dynamic information exchange across task-relevant neural ensembles. Although oscillation-based FC is classically defined relative to a prestimulus baseline, giving rise to rapid, context-dependent changes in individual connections, studies of distributed spatial patterns show that oscillation-based FC is omnipresent, occurring even in the absence of explicit cognitive demands. Thus, the issue of whether oscillation-based FC is primarily shaped by cognitive state or is intrinsic in nature remains open. Accordingly, we sought to reconcile these observations by interrogating the ECoG recordings of 18 presurgical human patients (8 females) for state dependence of oscillation-based FC in five canonical frequency bands across an array of six task states. FC analysis of phase and amplitude coupling revealed a highly similar, largely state-invariant (i.e., intrinsic) spatial component across cognitive states. This spatial organization was shared across all frequency bands. Crucially, however, each band also exhibited temporally independent FC dynamics capable of supporting frequency-specific information exchange. In conclusion, the spatial organization of oscillation-based FC is largely stable over cognitive states (i.e., primarily intrinsic in nature) and shared across frequency bands. Together, our findings converge with previous observations of spatially invariant patterns of FC derived from extremely slow and aperiodic fluctuations in fMRI signals. Our observations indicate that "background" FC should be accounted for in conceptual frameworks of oscillation-based FC targeting task-related changes.SIGNIFICANCE STATEMENT A fundamental property of neural activity is that it is periodic, enabling functional connectivity (FC) between distant regions through coupling of their oscillations. According to task-based studies, such oscillation-based FC is rapid and malleable to meet cognitive task demands. Studying distributed FC patterns instead of FC in a few individual connections, we found that oscillation-based FC is largely stable across various cognitive states and shares a common layout across oscillation frequencies. This stable spatial organization of FC in fast oscillatory brain signals parallels the known stability of fMRI-based intrinsic FC architecture. Despite the observed spatial state and frequency invariance, FC of individual connections was temporally independent between frequency bands, suggesting a putative mechanism for malleable frequency-specific FC to support cognitive tasks.

Functional connectivity of EEG is subject-specific, associated with phenotype, and different from fMRI

A variety of psychiatric, behavioral and cognitive phenotypes have been linked to brain ''functional connectivity'' -- the pattern of correlation observed between different brain regions. Most commonly assessed using functional magnetic resonance imaging (fMRI), here, we investigate the connectivity-phenotype associations with functional connectivity measured with electroencephalography (EEG), using phase-coupling. We analyzed data from the publicly available Healthy Brain Network Biobank. This database compiles a growing sample of children and adolescents, currently encompassing 1657 individuals. Among a variety of assessment instruments we focus on ten phenotypic and additional demographic measures that capture most of the variance in this sample. The largest effect sizes are found for age and sex for both fMRI and EEG. We replicate previous findings of an association of Intelligence Quotient (IQ) and Attention Deficit Hyperactivity Disorder (ADHD) with the pattern of fMRI functional connectivity. We also find an association with socioeconomic status, anxiety and the Child Behavior Checklist Score. For EEG we find a significant connectivity-phenotype relationship with IQ. The actual spatial patterns of functional connectivity are quite different between fMRI and source-space EEG. However, within EEG we observe clusters of functional connectivity that are consistent across frequency bands. Additionally we analyzed reproducibility of functional connectivity. We compare connectivity obtained with different tasks, including resting state, a video and a visual flicker task. For both EEG and fMRI the variation between tasks was smaller than the variability observed between subjects. We also found an increase of reliability with increasing frequency of the EEG, and increased sampling duration. We conclude that, while the patterns of functional connectivity are distinct between fMRI and phase-coupling of EEG, they are nonetheless similar in their robustness to the task, and similar in that idiosyncratic patterns of connectivity predict individual phenotypes.

Multi-timescale hybrid components of the functional brain connectome: A bimodal EEG-fMRI decomposition

Concurrent electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) bridge brain connectivity across timescales. During concurrent EEG-fMRI resting-state recordings, whole-brain functional connectivity (FC) strength is spatially correlated across modalities. However, cross-modal investigations have commonly remained correlational, and joint analysis of EEG-fMRI connectivity is largely unexplored. Here we investigated if there exist (spatially) independent FC networks linked between modalities. We applied the recently proposed hybrid connectivity independent component analysis (connICA) framework to two concurrent EEG-fMRI resting-state datasets (total 40 subjects). Two robust components were found across both datasets. The first component has a uniformly distributed EEG frequency fingerprint linked mainly to intrinsic connectivity networks (ICNs) in both modalities. Conversely, the second component is sensitive to different EEG frequencies and is primarily linked to intra-ICN connectivity in fMRI but to inter-ICN connectivity in EEG. The first hybrid component suggests that connectivity dynamics within well-known ICNs span timescales, from millisecond range in all canonical frequencies of FC_EEG to second range of FC_fMRI. Conversely, the second component additionally exposes linked but spatially divergent neuronal processing at the two timescales. This work reveals the existence of joint spatially independent components, suggesting that parts of resting-state connectivity are co-expressed in a linked manner across EEG and fMRI over individuals.

Functional connectivity is governed by a whole-brain organization measurable over multiple timescales by functional magnetic resonance imaging (fMRI) and electroencephalography (EEG). The relationship across the whole-brain organization captured at the different timescales of EEG and fMRI is largely unknown. Using concurrent EEG-fMRI, we identified spatially independent components consisting of brain connectivity patterns that co-occur in EEG and fMRI over subjects. We observed a component with similar connectivity organization across EEG and fMRI as well as a component with divergent connectivity. The former component governed all EEG frequencies while the latter was modulated by frequency. These findings show that part of functional connectivity organizes in a common spatial layout over several timescales, while a spatially independent part is modulated by frequency-specific information.

Magnetoencephalography resting state connectivity patterns as indicatives of surgical outcome in epilepsy patients

OBJECTIVE

Focal epilepsy is a disorder affecting several brain networks; however, epilepsy surgery usually targets a restricted region, the so-called epileptic focus. There is a growing interest in embedding resting state (RS) connectivity analysis into pre-surgical workup.

APPROACH

In this retrospective study, we analyzed Magnetoencephalography (MEG) long-range RS functional connectivity patterns in patients with drug-resistant focal epilepsy. MEG recorded prior to surgery from seven seizure-free (Engel Ia) and five non seizure-free (Engel III or IV) patients were analyzed (minimum 2-years post-surgical follow-up). MEG segments without any detectable epileptic activity were source localized using wavelet-based Maximum Entropy on the Mean method. Amplitude envelope correlation in the theta (4-8 Hz), alpha (8-13 Hz), and beta (13-26 Hz) bands were used for assessing connectivity.

MAIN RESULTS

For seizure-free patients, we found an isolated epileptic network characterized by weaker connections between the brain region where interictal epileptic discharges (IED) are generated and the rest of the cortex, when compared to connectivity between the corresponding contralateral homologous region and the rest of the cortex. Contrarily, non seizure-free patients exhibited a widespread RS epileptic network characterized by stronger connectivity between the IED generator and the rest of the cortex, in comparison to the contralateral region and the cortex. Differences between the two seizure outcome groups concerned mainly distant long-range connections and were found in the alpha-band.

SIGNIFICANCE

Importantly, these connectivity patterns suggest specific mechanisms describing the underlying organization of the epileptic network and were detectable at the individual patient level, supporting the prospect use of MEG connectivity patterns in epilepsy to predict post-surgical seizure outcome.

Concurrent EEG- and fMRI-derived functional connectomes exhibit linked dynamics

Long-range connectivity has become the most studied feature of human functional Magnetic Resonance Imaging (fMRI), yet the spatial and temporal relationship between its whole-brain dynamics and electrophysiological connectivity remains largely unknown. FMRI-derived functional connectivity exhibits spatial reconfigurations or time-varying dynamics at infraslow (<0.1Hz) speeds. Conversely, electrophysiological connectivity is based on cross-region coupling of fast oscillations (~1-100Hz). It is unclear whether such fast oscillation-based coupling varies at infraslow speeds, temporally coinciding with infraslow dynamics across the fMRI-based connectome. If so, does the association of fMRI-derived and electrophysiological dynamics spatially vary over the connectome across the functionally distinct electrophysiological oscillation bands? In two concurrent electroencephalography (EEG)-fMRI resting-state datasets, oscillation-based coherence in all canonical bands (delta through gamma) indeed reconfigured at infraslow speeds in tandem with fMRI-derived connectivity changes in corresponding region-pairs. Interestingly, irrespective of EEG frequency-band the cross-modal tie of connectivity dynamics comprised a large proportion of connections distributed across the entire connectome. However, there were frequency-specific differences in the relative strength of the cross-modal association. This association was strongest in visual to somatomotor connections for slower EEG-bands, and in connections involving the Default Mode Network for faster EEG-bands. Methodologically, the findings imply that neural connectivity dynamics can be reliably measured by fMRI despite heavy susceptibility to noise, and by EEG despite shortcomings of source reconstruction. Biologically, the findings provide evidence that contrast with known territories of oscillation power, oscillation coupling in all bands slowly reconfigures in a highly distributed manner across the whole-brain connectome.

Combining magnetoencephalography with magnetic resonance imaging enhances learning of surrogate-biomarkers

Electrophysiological methods, that is M/EEG, provide unique views into brain health. Yet, when building predictive models from brain data, it is often unclear how electrophysiology should be combined with other neuroimaging methods. Information can be redundant, useful common representations of multimodal data may not be obvious and multimodal data collection can be medically contraindicated, which reduces applicability. Here, we propose a multimodal model to robustly combine MEG, MRI and fMRI for prediction. We focus on age prediction as a surrogate biomarker in 674 subjects from the Cam-CAN dataset. Strikingly, MEG, fMRI and MRI showed additive effects supporting distinct brain-behavior associations. Moreover, the contribution of MEG was best explained by cortical power spectra between 8 and 30 Hz. Finally, we demonstrate that the model preserves benefits of stacking when some data is missing. The proposed framework, hence, enables multimodal learning for a wide range of biomarkers from diverse types of brain signals.

Genuine cross-frequency coupling networks in human resting-state electrophysiological recordings

Phase synchronization of neuronal oscillations in specific frequency bands coordinates anatomically distributed neuronal processing and communication. Typically, oscillations and synchronization take place concurrently in many distinct frequencies, which serve separate computational roles in cognitive functions. While within-frequency phase synchronization has been studied extensively, less is known about the mechanisms that govern neuronal processing distributed across frequencies and brain regions. Such integration of processing between frequencies could be achieved via cross-frequency coupling (CFC), either by phase–amplitude coupling (PAC) or by n:m-cross–frequency phase synchrony (CFS). So far, studies have mostly focused on local CFC in individual brain regions, whereas the presence and functional organization of CFC between brain areas have remained largely unknown. We posit that interareal CFC may be essential for large-scale coordination of neuronal activity and investigate here whether genuine CFC networks are present in human resting-state (RS) brain activity. To assess the functional organization of CFC networks, we identified brain-wide CFC networks at mesoscale resolution from stereoelectroencephalography (SEEG) and at macroscale resolution from source-reconstructed magnetoencephalography (MEG) data. We developed a novel, to our knowledge, graph-theoretical method to distinguish genuine CFC from spurious CFC that may arise from nonsinusoidal signals ubiquitous in neuronal activity. We show that genuine interareal CFC is present in human RS activity in both SEEG and MEG data. Both CFS and PAC networks coupled theta and alpha oscillations with higher frequencies in large-scale networks connecting anterior and posterior brain regions. CFS and PAC networks had distinct spectral patterns and opposing distribution of low- and high-frequency network hubs, implying that they constitute distinct CFC mechanisms. The strength of CFS networks was also predictive of cognitive performance in a separate neuropsychological assessment. In conclusion, these results provide evidence for interareal CFS and PAC being 2 distinct mechanisms for coupling oscillations across frequencies in large-scale brain networks.

Genuine interareal cross-frequency coupling (CFC) can be identified from human resting state activity using magnetoencephalography, stereoelectroencephalography, and novel network approaches. CFC couples slow theta and alpha oscillations to faster oscillations across brain regions.