Data and model considerations for estimating time-varying functional connectivity in fMRI.

Functional connectivity (FC) in the brain has been shown to exhibit subtle but reliable modulations within a session. One way of estimating time-varying FC is by using state-based models that describe fMRI time series as temporal sequences of states, each with an associated, characteristic pattern of FC. However, the estimation of these models from data sometimes fails to capture changes in a meaningful way, such that the model estimation assigns entire sessions (or the largest part of them) to a single state, therefore failing to capture within-session state modulations effectively; we refer to this phenomenon as the model becoming static, or model stasis. Here, we aim to quantify how the nature of the data and the choice of model parameters affect the model's ability to detect temporal changes in FC using both simulated fMRI time courses and resting state fMRI data. We show that large between-subject FC differences can overwhelm subtler within-session modulations, causing the model to become static. Further, the choice of parcellation can also affect the model's ability to detect temporal changes. We finally show that the model often becomes static when the number of free parameters per state that need to be estimated is high and the number of observations available for this estimation is low in comparison. Based on these findings, we derive a set of practical recommendations for time-varying FC studies, in terms of preprocessing, parcellation and complexity of the model.

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.

Transient spectral events in resting state MEG predict individual task responses.

Even in response to simple tasks such as hand movement, human brain activity shows remarkable inter-subject variability. Recently, it has been shown that individual spatial variability in fMRI task responses can be predicted from measurements collected at rest; suggesting that the spatial variability is a stable feature, inherent to the individual's brain. However, it is not clear if this is also true for individual variability in the spatio-spectral content of oscillatory brain activity. Here, we show using MEG (N â€‹= â€‹89) that we can predict the spatial and spectral content of an individual's task response using features estimated from the individual's resting MEG data. This works by learning when transient spectral 'bursts' or events in the resting state tend to reoccur in the task responses. We applied our method to motor, working memory and language comprehension tasks. All task conditions were predicted significantly above chance. Finally, we found a systematic relationship between genetic similarity (e.g. unrelated subjects vs. twins) and predictability. Our approach can predict individual differences in brain activity and suggests a link between transient spectral events in task and rest that can be captured at the level of individuals.

A dynamic system of brain networks revealed by fast transient EEG fluctuations and their fMRI correlates.

Resting state brain activity has become a significant area of investigation in human neuroimaging. An important approach for understanding the dynamics of neuronal activity in the resting state is to use complementary imaging modalities. Electrophysiological recordings can access fast temporal dynamics, while functional magnetic resonance imaging (fMRI) studies reveal detailed spatial patterns. However, the relationship between these two measures is not fully established. In this study, we used simultaneously recorded electroencephalography (EEG) and fMRI, along with Hidden Markov Modelling, to investigate how network dynamics at fast sub-second time-scales, accessible with EEG, link to the slower time-scales and higher spatial detail of fMRI. We found that the fMRI correlates of fast transient EEG dynamic networks show highly reproducible spatial patterns, and that their spatial organization exhibits strong similarity with traditional fMRI resting state networks maps. This further demonstrates the potential of electrophysiology as a tool for understanding the fast network dynamics that underlie fMRI resting state networks.

Discovery of key whole-brain transitions and dynamics during human wakefulness and non-REM sleep.

The modern understanding of sleep is based on the classification of sleep into stages defined by their electroencephalography (EEG) signatures, but the underlying brain dynamics remain unclear. Here we aimed to move significantly beyond the current state-of-the-art description of sleep, and in particular to characterise the spatiotemporal complexity of whole-brain networks and state transitions during sleep. In order to obtain the most unbiased estimate of how whole-brain network states evolve through the human sleep cycle, we used a Markovian data-driven analysis of continuous neuroimaging data from 57 healthy participants falling asleep during simultaneous functional magnetic resonance imaging (fMRI) and EEG. This Hidden Markov Model (HMM) facilitated discovery of the dynamic choreography between different whole-brain networks across the wake-non-REM sleep cycle. Notably, our results reveal key trajectories to switch within and between EEG-based sleep stages, while highlighting the heterogeneities of stage N1 sleep and wakefulness before and after sleep.

A new feature extraction method for signal classification applied to cord dorsum potential detection.

In the spinal cord of the anesthetized cat, spontaneous cord dorsum potentials (CDPs) appear synchronously along the lumbo-sacral segments. These CDPs have different shapes and magnitudes. Previous work has indicated that some CDPs appear to be specially associated with the activation of spinal pathways that lead to primary afferent depolarization and presynaptic inhibition. Visual detection and classification of these CDPs provides relevant information on the functional organization of the neural networks involved in the control of sensory information and allows the characterization of the changes produced by acute nerve and spinal lesions. We now present a novel feature extraction approach for signal classification, applied to CDP detection. The method is based on an intuitive procedure. We first remove by convolution the noise from the CDPs recorded in each given spinal segment. Then, we assign a coefficient for each main local maximum of the signal using its amplitude and distance to the most important maximum of the signal. These coefficients will be the input for the subsequent classification algorithm. In particular, we employ gradient boosting classification trees. This combination of approaches allows a faster and more accurate discrimination of CDPs than is obtained by other methods.