Dynamic BOLD functional connectivity in humans and its electrophysiological correlates

Neural oscillations subserve many human perceptual and cognitive operations. Accordingly, brain functional connectivity is not static in time, but fluctuates dynamically following the synchronization and desynchronization of neural populations. This dynamic functional connectivity has recently been demonstrated in spontaneous fluctuations of the Blood Oxygen Level-Dependent (BOLD) signal, measured with functional Magnetic Resonance Imaging (fMRI). We analyzed temporal fluctuations in BOLD connectivity and their electrophysiological correlates, by means of long (≈50 min) joint electroencephalographic (EEG) and fMRI recordings obtained from two populations: 15 awake subjects and 13 subjects undergoing vigilance transitions. We identified positive and negative correlations between EEG spectral power (extracted from electrodes covering different scalp regions) and fMRI BOLD connectivity in a network of 90 cortical and subcortical regions (with millimeter spatial resolution). In particular, increased alpha (8-12 Hz) and beta (15-30 Hz) power were related to decreased functional connectivity, whereas gamma (30-60 Hz) power correlated positively with BOLD connectivity between specific brain regions. These patterns were altered for subjects undergoing vigilance changes, with slower oscillations being correlated with functional connectivity increases. Dynamic BOLD functional connectivity was reflected in the fluctuations of graph theoretical indices of network structure, with changes in frontal and central alpha power correlating with average path length. Our results strongly suggest that fluctuations of BOLD functional connectivity have a neurophysiological origin. Positive correlations with gamma can be interpreted as facilitating increased BOLD connectivity needed to integrate brain regions for cognitive performance. Negative correlations with alpha suggest a temporary functional weakening of local and long-range connectivity, associated with an idling state

Altered fMRI Connectivity Dynamics in Temporal Lobe Epilepsy Might Explain Seizure Semiology

Temporal lobe epilepsy (TLE) can be conceptualized as a network disease. The network can be characterized by inter-regional functional connectivity, i.e., blood oxygen level-dependent (BOLD) signal correlations between any two regions. However, functional connectivity is not constant over time, thus computing correlation at a given time and then at some later time could give different results (non-stationarity). We hypothesized (1) that non-stationarities can be induced by epilepsy (e.g., interictal epileptic activity) increasing local signal variance and that (2) these transient events contribute to fluctuations in connectivity leading to pathological functioning, i.e., TLE semiology. We analyzed fMRI data from 27 patients with TLE and 22 healthy controls focusing on EEG-confirmed wake epochs only to protect against sleep-induced connectivity changes. Testing hypothesis (1), we identified brain regions where the BOLD signal variance was significantly greater in TLE than in controls: the temporal pole – including the hippocampus. Taking the latter as the seed region and testing hypothesis (2), we calculated the time-varying inter-regional correlation values (dynamic functional connectivity) to other brain regions and found greater connectivity variance in the TLE than the control group mainly in the precuneus, the supplementary and sensorimotor, and the frontal cortices. We conclude that the highest BOLD signal variance in the hippocampi is highly suggestive of a specific epilepsy-related effect. The altered connectivity dynamics in TLE patients might help to explain the hallmark semiological features of dyscognitive seizures including impaired consciousness (precuneus, frontal cortex), sensory disturbance, and motor automatisms (sensorimotor cortices, supplementary motor cortex). Accounting for the non-stationarity and state-dependence of functional connectivity are a prerequisite in the search for potential connectivity-derived biomarkers in TLE

EEG Microstate Sequences From Different Clustering Algorithms Are Information-Theoretically Invariant

We analyse statistical and information-theoretical properties of EEG microstate sequences, as seen through the lens of five different clustering algorithms. Microstate sequences are computed for n = 20 resting state EEG recordings during wakeful rest. The input for all clustering algorithms is the set of EEG topographic maps obtained at local maxima of the spatial variance. This data set is processed by two classical microstate clustering algorithms (1) atomize and agglomerate hierarchical clustering (AAHC) and (2) a modified K-means algorithm, as well as by (3) K-medoids, (4) principal component analysis (PCA) and (5) fast independent component analysis (Fast-ICA). Using this technique, EEG topographies can be substituted with microstate labels by competitive fitting based on spatial correlation, resulting in a symbolic, non-metric time series, the microstate sequence. Microstate topographies and symbolic time series are further analyzed statistically, including static and dynamic properties. Static properties, which do not contain information about temporal dependencies of the microstate sequence include the maximum similarity of microstate maps within and between the tested clustering algorithms, the global explained variance and the Shannon entropy of the microstate sequences. Dynamic properties are sensitive to temporal correlations between the symbols and include the mixing time of the microstate transition matrix, the entropy rate of the microstate sequences and the location of the first local maximum of the autoinformation function. We also test the Markov property of microstate sequences, the time stationarity of the transition matrix and detect periodicities by means of time-lagged mutual information. Finally, possible long-range correlations of microstate sequences are assessed via Hurst exponent estimation. We find that while static properties partially reflect properties of the clustering algorithms, information-theoretical quantities are largely invariant with respect to the clustering method used. As each clustering algorithm has its own profile of computational speed, ease of implementation, determinism vs. stochasticity and theoretical underpinnings, our results convey a positive message concerning the free choice of method and the comparability of results obtained from different algorithms. The invariance of these quantities implies that the tested properties are algorithm-independent, inherent features of resting state EEG derived microstate sequences

Causal hierarchy within the thalamo-cortical network in spike and wave discharges

Background: Generalised spike wave (GSW) discharges are the electroencephalographic (EEG) hallmark of absence seizures, clinically characterised by a transitory interruption of ongoing activities and impaired consciousness, occurring during states of reduced awareness. Several theories have been proposed to explain the pathophysiology of GSW discharges and the role of thalamus and cortex as generators. In this work we extend the existing theories by hypothesizing a role for the precuneus, a brain region neglected in previous works on GSW generation but already known to be linked to consciousness and awareness. We analysed fMRI data using dynamic causal modelling (DCM) to investigate the effective connectivity between precuneus, thalamus and prefrontal cortex in patients with GSW discharges. Methodology and Principal Findings: We analysed fMRI data from seven patients affected by Idiopathic Generalized Epilepsy (IGE) with frequent GSW discharges and significant GSW-correlated haemodynamic signal changes in the thalamus, the prefrontal cortex and the precuneus. Using DCM we assessed their effective connectivity, i.e. which region drives another region. Three dynamic causal models were constructed: GSW was modelled as autonomous input to the thalamus (model A), ventromedial prefrontal cortex (model B), and precuneus (model C). Bayesian model comparison revealed Model C (GSW as autonomous input to precuneus), to be the best in 5 patients while model A prevailed in two cases. At the group level model C dominated and at the population-level the p value of model C was ∼1. Conclusion: Our results provide strong evidence that activity in the precuneus gates GSW discharges in the thalamo-(fronto) cortical network. This study is the first demonstration of a causal link between haemodynamic changes in the precuneus - an index of awareness - and the occurrence of pathological discharges in epilepsy. © 2009 Vaudano et al

The voxel-wise functional connectome can be efficiently derived from co-activations in a sparse spatio-temporal point-process

Large efforts are currently under way to systematically map functional connectivity between all pairs of millimeter-scale brain regions based on large neuroimaging databases. The exploratory unraveling of this "functional connectome" based on functional Magnetic Resonance Imaging (fMRI) can benefit from a better understanding of the contributors to resting state functional connectivity. In this work, we introduce a sparse representation of fMRI data in the form of a discrete point-process encoding high-amplitude events in the blood oxygenation level-dependent (BOLD) signal and we show it contains sufficient information for the estimation of functional connectivity between all pairs of voxels. We validate this method by replicating results obtained with standard whole-brain voxel-wise linear correlation matrices in two datasets. In the first one (n = 71), we study the changes in node strength (a measure of network centrality) during deep sleep. The second is a large database (n = 1147) of subjects in which we look at the age-related reorganization of the voxel-wise network of functional connections. In both cases it is shown that the proposed method compares well with standard techniques, despite requiring only data on the order of 1% of the original BOLD signal time series. Furthermore, we establish that the point-process approach does not reduce (and in one case increases) classification accuracy compared to standard linear correlations. Our results show how large fMRI datasets can be drastically simplified to include only the timings of large-amplitude events, while still allowing the recovery of all pair-wise interactions between voxels. The practical importance of this dimensionality reduction is manifest in the increasing number of collaborative efforts aiming to study large cohorts of healthy subjects as well as patients suffering from brain disease. Our method also suggests that the electrophysiological signals underlying the dynamics of fMRI time series consist of all-or-none temporally localized events, analogous to the avalanches of neural activity observed in recordings of local field potentials (LFP), an observation of potentially high neurobiological relevance.Fil: Tagliazucchi, Enzo. Christian Albrechts Universitat Zu Kiel.; Alemania. University Frankfurt am Main; AlemaniaFil: Siniatchkin, Michael. Christian Albrechts Universitat Zu Kiel.; AlemaniaFil: Laufs, Helmut. University Frankfurt am Main; Alemania. University Hospital Schleswig Holstein; AlemaniaFil: Chialvo, Dante Renato. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Universidad Nacional de San Martin. Escuela de Ciencia y Tecnologia. Centro de Estudios Multidisciplinarios En Sistemas Complejos y Ciencias del Cerebro.; Argentin

Connectivity dynamics from wakefulness to sleep

Interest in time-resolved connectivity in fMRI has grown rapidly in recent years. The most widely used technique for studying connectivity changes over time utilizes a sliding windows approach. There has been some debate about the utility of shorter versus longer windows, the use of fixed versus adaptive windows, as well as whether observed resting state dynamics during wakefulness may be predominantly due to changes in sleep state and subject head motion. In this work we use an independent component analysis (ICA)-based pipeline applied to concurrent EEG/fMRI data collected during wakefulness and various sleep stages and show: 1) connectivity states obtained from clustering sliding windowed correlations of resting state functional network time courses well classify the sleep states obtained from EEG data, 2) using shorter sliding windows instead of longer non-overlapping windows improves the ability to capture transition dynamics even at windows as short as 30 ​s, 3) motion appears to be mostly associated with one of the states rather than spread across all of them 4) a fixed tapered sliding window approach outperforms an adaptive dynamic conditional correlation approach, and 5) consistent with prior EEG/fMRI work, we identify evidence of multiple states within the wakeful condition which are able to be classified with high accuracy. Classification of wakeful only states suggest the presence of time-varying changes in connectivity in fMRI data beyond sleep state or motion. Results also inform about advantageous technical choices, and the identification of different clusters within wakefulness that are separable suggest further studies in this direction.Fil: Damaraju, Eswar. Instituto Tecnológico de Georgia; Estados UnidosFil: Tagliazucchi, Enzo Rodolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Laufs, Helmut. Goethe Universitat Frankfurt; AlemaniaFil: Calhoun, Vince D.. Instituto Tecnológico de Georgia; Estados Unido

Data-driven discovery of canonical large-scale brain dynamics

Human behavior and cognitive function correlate with complex patterns of spatio-temporal brain dynamics, which can be simulated using computational models with different degrees of biophysical realism. We used a data-driven optimization algorithm to determine and classify the types of local dynamics that enable the reproduction of different observables derived from functional magnetic resonance recordings. The phase space analysis of the resulting equations revealed a predominance of stable spiral attractors, which optimized the similarity to the empirical data in terms of the synchronization, metastability, and functional connectivity dynamics. For stable limit cycles, departures from harmonic oscillations improved the fit in terms of functional connectivity dynamics. Eigenvalue analyses showed that proximity to a bifurcation improved the accuracy of the simulation for wakefulness, whereas deep sleep was associated with increased stability. Our results provide testable predictions that constrain the landscape of suitable biophysical models, while supporting noise-driven dynamics close to a bifurcation as a canonical mechanism underlying the complex fluctuations that characterize endogenous brain activity.Fil: Piccinini, Juan Ignacio. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Deco, Gustavo. Universitat Pompeu Fabra; EspañaFil: Kringelbach, Morten. University of Oxford; Reino UnidoFil: Laufs, Helmut. University of Kiel; AlemaniaFil: Sanz Perl Hernandez, Yonatan. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Tagliazucchi, Enzo Rodolfo. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentin

FMRI spectral signatures of sleep

Sleep can be distinguished from wake by changes in brain electrical activity, typically assessed using electroencephalography (EEG). The hallmark of nonrapid-eye-movement (NREM) sleep is the shift from high-frequency, low-amplitude wake EEG to low-frequency, high-amplitude sleep EEG dominated by spindles and slow waves. Here we identified signatures of sleep in brain hemodynamic activity, using simultaneous functional MRI (fMRI) and EEG. We found that, at the transition from wake to sleep, fMRI blood oxygen level–dependent (BOLD) activity evolved from a mixed-frequency pattern to one dominated by two distinct oscillations: a low-frequency (0.1 Hz) prominent in deep sleep and correlated with the occurrence of slow waves. The two oscillations were both detectable across the brain but exhibited distinct spatiotemporal patterns. During the falling-asleep process, the low-frequency oscillation first appeared in the thalamus, then the posterior cortex, and lastly the frontal cortex, while the high-frequency oscillation first appeared in the midbrain, then the frontal cortex, and lastly the posterior cortex. During the waking-up process, both oscillations disappeared first from the thalamus, then the frontal cortex, and lastly the posterior cortex. The BOLD oscillations provide local signatures of spindle and slow wave activity. They may be employed to monitor the regional occurrence of sleep or wakefulness, track which regions are the first to fall asleep or wake up at the wake–sleep transitions, and investigate local homeostatic sleep processes