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

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

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

Validation of non-REM sleep stage decoding from resting state fMRI using linear support vector machines

A growing body of literature suggests that changes in consciousness are reflected in specific connectivity patterns of the brain as obtained from resting state fMRI (rs-fMRI). As simultaneous electroencephalography (EEG) is often unavailable, decoding of potentially confounding sleep patterns from rs-fMRI itself might be useful and improve data interpretation. Linear support vector machine classifiers were trained on combined rs-fMRI/EEG recordings from 25 subjects to separate wakefulness (S0) from non-rapid eye movement (NREM) sleep stages 1 (S1), 2 (S2), slow wave sleep (SW) and all three sleep stages combined (SX). Classifier performance was quantified by a leave-one-subject-out cross-validation (LOSO-CV) and on an independent validation dataset comprising 19 subjects. Results demonstrated excellent performance with areas under the receiver operating characteristics curve (AUCs) close to 1.0 for the discrimination of sleep from wakefulness (S0|SX), S0|S1, S0|S2 and S0|SW, and good to excellent performance for the classification between sleep stages (S1|S2:~0.9; S1|SW:~1.0; S2|SW:~0.8). Application windows of fMRI data from about 70 s were found as minimum to provide reliable classifications. Discrimination patterns pointed to subcortical-cortical connectivity and within-occipital lobe reorganization of connectivity as strongest carriers of discriminative information. In conclusion, we report that functional connectivity analysis allows valid classification of NREM sleep stages

ProductGraphSleepNet: Sleep Staging using Product Spatio-Temporal Graph Learning with Attentive Temporal Aggregation

The classification of sleep stages plays a crucial role in understanding and diagnosing sleep pathophysiology. Sleep stage scoring relies heavily on visual inspection by an expert that is time consuming and subjective procedure. Recently, deep learning neural network approaches have been leveraged to develop a generalized automated sleep staging and account for shifts in distributions that may be caused by inherent inter/intra-subject variability, heterogeneity across datasets, and different recording environments. However, these networks ignore the connections among brain regions, and disregard the sequential connections between temporally adjacent sleep epochs. To address these issues, this work proposes an adaptive product graph learning-based graph convolutional network, named ProductGraphSleepNet, for learning joint spatio-temporal graphs along with a bidirectional gated recurrent unit and a modified graph attention network to capture the attentive dynamics of sleep stage transitions. Evaluation on two public databases: the Montreal Archive of Sleep Studies (MASS) SS3; and the SleepEDF, which contain full night polysomnography recordings of 62 and 20 healthy subjects, respectively, demonstrates performance comparable to the state-of-the-art (Accuracy: 0.867;0.838, F1-score: 0.818;0.774 and Kappa: 0.802;0.775, on each database respectively). More importantly, the proposed network makes it possible for clinicians to comprehend and interpret the learned connectivity graphs for sleep stages.Comment: 9 pages, 6 figure

Metabolic connectivity of resting-state networks in alpha synucleinopathies, from prodromal to dementia phase

Previous evidence suggests that the derangement of large-scale brain networks reflects structural, molecular, and functional mechanisms underlying neurodegenerative diseases. Although the alterations of multiple large-scale brain networks in Parkinson’s disease (PD) and Dementia with Lewy Bodies (DLB) are reported, a comprehensive study on connectivity reconfiguration starting from the preclinical phase is still lacking. We aimed to investigate shared and disease-specific changes in the large-scale networks across the Lewy Bodies (LB) disorders spectrum using a brain metabolic connectivity approach. We included 30 patients with isolated REM sleep behavior disorder (iRBD), 28 with stable PD, 30 with DLB, and 30 healthy controls for comparison. We applied seed-based interregional correlation analyses (IRCA) to evaluate the metabolic connectivity in the large-scale resting-state networks, as assessed by [18F]FDG-PET, in each clinical group compared to controls. We assessed metabolic connectivity changes by applying the IRCA and specific connectivity metrics, such as the weighted and unweighted Dice similarity coefficients (DC), for the topographical similarities. All the investigated large-scale brain resting-state networks showed metabolic connectivity alterations, supporting the widespread involvement of brain connectivity within the alpha-synuclein spectrum. Connectivity alterations were already evident in iRBD, severely affecting the posterior default mode, attentive and limbic networks. Strong similarities emerged in iRBD and DLB that showed comparable connectivity alterations in most large-scale networks, particularly in the posterior default mode and attentive networks. Contrarily, PD showed the main connectivity alterations limited to motor and somatosensory networks. The present findings reveal that metabolic connectivity alterations in the large-scale networks are already present in the early iRBD phase, resembling the DLB metabolic connectivity changes. This suggests and confirms iRBD as a risk condition for progression to the severe LB disease phenotype. Of note, the neurobiology of stable PD supports its more benign phenotype

Spatiotemporal Modeling of Multivariate Signals With Graph Neural Networks and Structured State Space Models

Multivariate signals are prevalent in various domains, such as healthcare, transportation systems, and space sciences. Modeling spatiotemporal dependencies in multivariate signals is challenging due to (1) long-range temporal dependencies and (2) complex spatial correlations between sensors. To address these challenges, we propose representing multivariate signals as graphs and introduce GraphS4mer, a general graph neural network (GNN) architecture that captures both spatial and temporal dependencies in multivariate signals. Specifically, (1) we leverage Structured State Spaces model (S4), a state-of-the-art sequence model, to capture long-term temporal dependencies and (2) we propose a graph structure learning layer in GraphS4mer to learn dynamically evolving graph structures in the data. We evaluate our proposed model on three distinct tasks and show that GraphS4mer consistently improves over existing models, including (1) seizure detection from electroencephalography signals, outperforming a previous GNN with self-supervised pretraining by 3.1 points in AUROC; (2) sleep staging from polysomnography signals, a 4.1 points improvement in macro-F1 score compared to existing sleep staging models; and (3) traffic forecasting, reducing MAE by 8.8% compared to existing GNNs and by 1.4% compared to Transformer-based models