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

•

Time-varying FC models sometimes fail to detect temporal changes in fMRI data.

•

Between-subject and within-session FC variability affect model stasis.

•

The choice of parcellation affects model stasis in real fMRI data.

•

The number of observations and free parameters per state critically affect model stasis.

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.

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.

Sleep is composed of a number of different stages, each associated with a different pattern of brain activity. Here, using a data-driven Hidden Markov Model (HMM) of fMRI data, the authors discover a more complex set of neural activity states underlying the conventional stages of non-REM 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.

Accelerating forward genetics for cell wall deconstruction

The elucidation of the genes involved in cell wall synthesis and assembly remains one of the biggest challenges of cell wall biology. Although traditional genetic approaches, using simple yet elegant screens, have identified components of the cell wall, many unknowns remain. Exhausting the genetic toolbox by performing sensitized screens, adopting chemical genetics or combining these with improved cell wall imaging, hold the promise of new gene discovery and function. With the recent introduction of next-generation sequencing technologies, it is now possible to quickly and efficiently map and clone genes of interest in record time. The combination of a classical genetics approach and cutting edge technology will propel cell wall biology in plants forward into the future.

Next-generation mapping of Arabidopsis genes

Next-generation genomic sequencing technologies have made it possible to directly map mutations responsible for phenotypes of interest via direct sequencing. However, most mapping strategies proposed to date require some prior genetic analysis, which can be very time-consuming even in genetically tractable organisms. Here we present a de novo method for rapidly and robustly mapping the physical location of EMS mutations by sequencing a small pooled F₂ population. This method, called Next Generation Mapping (NGM), uses a chastity statistic to quantify the relative contribution of the parental mutant and mapping lines to each SNP in the pooled F₂ population. It then uses this information to objectively localize the candidate mutation based on its exclusive segregation with the mutant parental line. A user-friendly, web-based tool for performing NGM analysis is available at http://bar.utoronto.ca/NGM. We used NGM to identify three genes involved in cell-wall biology in Arabidopsis thaliana, and, in a power analysis, demonstrate success in test mappings using as few as ten F₂ lines and a single channel of Illumina Genome Analyzer data. This strategy can easily be applied to other model organisms, and we expect that it will also have utility in crops and any other eukaryote with a completed genome sequence.