For several years it has been argued that neural synchronisation is crucial for cognition. The idea that synchronised temporal patterns between different neural groups carries information above and beyond the isolated activity of these groups has inspired a shift in focus in the field of functional neuroimaging. Specifically, investigation into the activation elicited within certain regions by some stimulus or task has, in part, given way to analysis of patterns of co-activation or functional connectivity between distal regions. Recently, the functional connectivity community has been looking beyond the assumptions of stationarity that earlier work was based on, and has introduced methods to incorporate temporal dynamics into the analysis of connectivity. In particular, non-invasive electrophysiological data (magnetoencephalography/electroencephalography (MEG/EEG)), which provides direct measurement of whole-brain activity and rich temporal information, offers an exceptional window into such (potentially fast) brain dynamics. In this review, we discuss challenges, solutions, and a collection of analysis tools that have been developed in recent years to facilitate the investigation of dynamic functional connectivity using these imaging modalities. Further, we discuss the applications of these approaches in the study of cognition and neuropsychiatric disorders. Finally, we review some existing developments that, by using realistic computational models, pursue a deeper understanding of the underlying causes of non-stationary connectivity.

No abstract

Network connectivity is an integral feature of human brain function, and characterising its maturational trajectory is a critical step towards understanding healthy and atypical neurodevelopment. Here, we used magnetoencephalography (MEG) to investigate both stationary (i.e. time averaged) and rapidly modulating (dynamic) electrophysiological connectivity, in participants aged from mid-childhood to early adulthood (youngest participant 9 years old; oldest participant 25 years old). Stationary functional connectivity (measured via inter-regional coordination of neural oscillations) increased with age in the alpha and beta frequency bands, particularly in bilateral parietal and temporo-parietal connections. Our dynamic analysis (also applied to alpha/beta oscillations) revealed the spatiotemporal signatures of 8 dynamic networks; these modulate on a ∼100 ms time scale, and temporal stability in attentional networks was found to increase with age. Significant overlap was found between age-modulated dynamic networks and inter-regional oscillatory coordination, implying that altered network dynamics underlie age related changes in functional connectivity. Our results provide novel insights into brain network electrophysiology, and lay a foundation for future work in childhood disorders.

Fluctuations in functional interactions between brain regions typically occur at the millisecond time scale. Conventional connectivity metrics are not adequately time-resolved to detect such fast fluctuations in functional connectivity. At the same time, attempts to use conventional metrics in a time-resolved manner usually come with the selection of sliding windows of fixed arbitrary length. In the current work, we evaluated the use of high temporal resolution metrics of functional connectivity in conjunction with non-negative tensor factorisation to detect fast fluctuations in connectivity and temporally evolving subnetworks. To this end, we used the phase difference derivative and wavelet coherence, and we also introduced a new metric, the instantaneous amplitude correlation. In order to deal with the inherently noisy nature of magnetoencephalography data and large datasets, we make use of recurrence plots and we used pair-wise orthogonalisation to avoid spurious estimates of functional connectivity due to signal leakage. Firstly, metrics were evaluated in the context of dynamically coupled neural mass models in the presence and absence of delays and also compared to conventional static metrics with fixed sliding windows. Simulations showed that these high temporal resolution metrics outperformed conventional static connectivity metrics. Secondly, the sensitivity of the metrics to fluctuations in connectivity was analysed in post-movement beta rebound magnetoencephalography data, which showed time locked sensorimotor subnetworks that modulated with the beta rebound. Finally, sensitivity of the metrics was evaluated in resting-state magnetoencephalography, showing similar spatial patterns across metrics, thereby indicating the robustness of the current analysis. The current methods can be applied in cognitive experiments that involve fast modulations in connectivity in relation to cognition. In addition, these methods could also be used as input to temporal graph analysis to further characterise the rapid fluctuation in brain network topology.

The study of functional connectivity using magnetoencephalography (MEG) is an expanding area of neuroimaging, and adds an extra dimension to the more common assessments made using fMRI. The importance of such metrics is growing, with recent demonstrations of their utility in clinical research, however previous reports suggest that whilst group level resting state connectivity is robust, single session recordings lack repeatability. Such robustness is critical if MEG measures in individual subjects are to prove clinically valuable. In the present paper, we test how practical aspects of experimental design affect the intra-subject repeatability of MEG findings; specifically we assess the effect of coregistration method and data recording duration. We show that the use of a foam head-cast, which is known to improve coregistration accuracy, increased significantly the between session repeatability of both beamformer reconstruction and connectivity estimation. We also show that recording duration is a critical parameter, with large improvements in repeatability apparent when using ten minute, compared to five minute recordings. Further analyses suggest that the origin of this latter effect is not underpinned by technical aspects of source reconstruction, but rather by a genuine effect of brain state; short recordings are simply inefficient at capturing the canonical MEG network in a single subject. Our results provide important insights on experimental design and will prove valuable for future MEG connectivity studies.3

Despite advances in the field of dynamic connectivity, fixed sliding window approaches for the detection of fluctuations in functional connectivity are still widely used. The use of conventional connectivity metrics in conjunction with a fixed sliding window comes with the arbitrariness of the chosen window lengths. In this paper we use multivariate autoregressive and neural mass models with a priori defined ground truths to systematically analyze the sensitivity of conventional metrics in combination with different window lengths to detect genuine fluctuations in connectivity for various underlying state durations. Metrics of interest are the coherence, imaginary coherence, phase lag index, phase locking value and the amplitude envelope correlation. We performed analysis for two nodes and at the network level. We demonstrate that these metrics show indeed higher variability for genuine temporal fluctuations in connectivity compared to a static connectivity state superimposed by noise. Overall, the error of the connectivity estimates themselves decreases for longer state durations (order of seconds), while correlations of the connectivity fluctuations with the ground truth was higher for longer state durations. In general, metrics, in combination with a sliding window, perform poorly for very short state durations. Increasing the SNR of the system only leads to a moderate improvement. In addition, at the network level, only longer window widths were sufficient to detect plausible resting state networks that matched the underlying ground truth, especially for the phase locking value, amplitude envelope correlation and coherence. The length of these longer window widths did not necessarily correspond to the underlying state durations. For short window widths resting state network connectivity patterns could not be retrieved. We conclude that fixed sliding window approaches for connectivity can detect modulations of connectivity, but mostly if the underlying dynamics operate on moderate to slow timescales. In practice, this can be a drawback, as state durations can vary significantly in empirical data.

No abstract

Despite its omnipresence in everyday interactions and its importance for mental health, mood and its neuronal underpinnings are poorly understood. Computational models can help identify parameters affecting self-reported mood during mood induction tasks. Here, we test if computationally modeled dynamics of self-reported mood during monetary gambling can be used to identify trial-by-trial variations in neuronal activity. To this end, we shifted mood in healthy (N = 24) and depressed (N = 30) adolescents by delivering individually tailored reward prediction errors while recording magnetoencephalography (MEG) data. Following a pre-registered analysis, we hypothesize that the expectation component of mood would be predictive of beta-gamma oscillatory power (25–40 Hz). We also hypothesize that trial variations in the source localized responses to reward feedback would be predicted by mood and by its reward prediction error component. Through our multilevel statistical analysis, we found confirmatory evidence that beta-gamma power is positively related to reward expectation during mood shifts, with localized sources in the posterior cingulate cortex. We also confirmed reward prediction error to be predictive of trial-level variations in the response of the paracentral lobule. To our knowledge, this is the first study to harness computational models of mood to relate mood fluctuations to variations in neural oscillations with MEG.