Mixtures of large-scale dynamic functional brain network modes.

Accurate temporal modelling of functional brain networks is essential in the quest for understanding how such networks facilitate cognition. Researchers are beginning to adopt time-varying analyses for electrophysiological data that capture highly dynamic processes on the order of milliseconds. Typically, these approaches, such as clustering of functional connectivity profiles and Hidden Markov Modelling (HMM), assume mutual exclusivity of networks over time. Whilst a powerful constraint, this assumption may be compromising the ability of these approaches to describe the data effectively. Here, we propose a new generative model for functional connectivity as a time-varying linear mixture of spatially distributed statistical "modes". The temporal evolution of this mixture is governed by a recurrent neural network, which enables the model to generate data with a rich temporal structure. We use a Bayesian framework known as amortised variational inference to learn model parameters from observed data. We call the approach DyNeMo (for Dynamic Network Modes), and show using simulations it outperforms the HMM when the assumption of mutual exclusivity is violated. In resting-state MEG, DyNeMo reveals a mixture of modes that activate on fast time scales of 100-150 ms, which is similar to state lifetimes found using an HMM. In task MEG data, DyNeMo finds modes with plausible, task-dependent evoked responses without any knowledge of the task timings. Overall, DyNeMo provides decompositions that are an approximate remapping of the HMM's while showing improvements in overall explanatory power. However, the magnitude of the improvements suggests that the HMM's assumption of mutual exclusivity can be reasonable in practice. Nonetheless, DyNeMo provides a flexible framework for implementing and assessing future modelling developments.

osl-dynamics, a toolbox for modeling fast dynamic brain activity.

Neural activity contains rich spatiotemporal structure that corresponds to cognition. This includes oscillatory bursting and dynamic activity that span across networks of brain regions, all of which can occur on timescales of tens of milliseconds. While these processes can be accessed through brain recordings and imaging, modeling them presents methodological challenges due to their fast and transient nature. Furthermore, the exact timing and duration of interesting cognitive events are often a priori unknown. Here, we present the OHBA Software Library Dynamics Toolbox (osl-dynamics), a Python-based package that can identify and describe recurrent dynamics in functional neuroimaging data on timescales as fast as tens of milliseconds. At its core are machine learning generative models that are able to adapt to the data and learn the timing, as well as the spatial and spectral characteristics, of brain activity with few assumptions. osl-dynamics incorporates state-of-the-art approaches that can be, and have been, used to elucidate brain dynamics in a wide range of data types, including magneto/electroencephalography, functional magnetic resonance imaging, invasive local field potential recordings, and electrocorticography. It also provides novel summary measures of brain dynamics that can be used to inform our understanding of cognition, behavior, and disease. We hope osl-dynamics will further our understanding of brain function, through its ability to enhance the modeling of fast dynamic processes.

The cortical neurophysiological signature of amyotrophic lateral sclerosis.

The progressive loss of motor function characteristic of amyotrophic lateral sclerosis is associated with widespread cortical pathology extending beyond primary motor regions. Increasing muscle weakness reflects a dynamic, variably compensated brain network disorder. In the quest for biomarkers to accelerate therapeutic assessment, the high temporal resolution of magnetoencephalography is uniquely able to non-invasively capture micro-magnetic fields generated by neuronal activity across the entire cortex simultaneously. This study examined task-free magnetoencephalography to characterize the cortical oscillatory signature of amyotrophic lateral sclerosis for having potential as a pharmacodynamic biomarker. Eight to ten minutes of magnetoencephalography in the task-free, eyes-open state was recorded in amyotrophic lateral sclerosis (n = 36) and healthy age-matched controls (n = 51), followed by a structural MRI scan for co-registration. Extracted magnetoencephalography metrics from the delta, theta, alpha, beta, low-gamma, high-gamma frequency bands included oscillatory power (regional activity), 1/f exponent (complexity) and amplitude envelope correlation (connectivity). Groups were compared using a permutation-based general linear model with correction for multiple comparisons and confounders. To test whether the extracted metrics could predict disease severity, a random forest regression model was trained and evaluated using nested leave-one-out cross-validation. Amyotrophic lateral sclerosis was characterized by reduced sensorimotor beta band and increased high-gamma band power. Within the premotor cortex, increased disability was associated with a reduced 1/f exponent. Increased disability was more widely associated with increased global connectivity in the delta, theta and high-gamma bands. Intra-hemispherically, increased disability scores were particularly associated with increases in temporal connectivity and inter-hemispherically with increases in frontal and occipital connectivity. The random forest model achieved a coefficient of determination (R2) of 0.24. The combined reduction in cortical sensorimotor beta and rise in gamma power is compatible with the established hypothesis of loss of inhibitory, GABAergic interneuronal circuits in pathogenesis. A lower 1/f exponent potentially reflects a more excitable cortex and a pathology unique to amyotrophic lateral sclerosis when considered with the findings published in other neurodegenerative disorders. Power and complexity changes corroborate with the results from paired-pulse transcranial magnetic stimulation. Increased magnetoencephalography connectivity in worsening disability is thought to represent compensatory responses to a failing motor system. Restoration of cortical beta and gamma band power has significant potential to be tested in an experimental medicine setting. Magnetoencephalography-based measures have potential as sensitive outcome measures of therapeutic benefit in drug trials and may have a wider diagnostic value with further study, including as predictive markers in asymptomatic carriers of disease-causing genetic variants.

Multi-dynamic modelling reveals strongly time-varying resting fMRI correlations.

The activity of functional brain networks is responsible for the emergence of time-varying cognition and behaviour. Accordingly, time-varying correlations (Functional Connectivity) in resting fMRI have been shown to be predictive of behavioural traits, and psychiatric and neurological conditions. Typically, methods that measure time varying Functional Connectivity (FC), such as sliding windows approaches, do not separately model when changes occur in the mean activity levels from when changes occur in the FC, therefore conflating these two distinct types of modulation. We show that this can bias the estimation of time-varying FC to appear more stable over time than it actually is. Here, we propose an alternative approach that models changes in the mean brain activity and in the FC as being able to occur at different times to each other. We refer to this method as the Multi-dynamic Adversarial Generator Encoder (MAGE) model, which includes a model of the network dynamics that captures long-range time dependencies, and is estimated on fMRI data using principles of Generative Adversarial Networks. We evaluated the approach across several simulation studies and resting fMRI data from the Human Connectome Project (1003 subjects), as well as from UK Biobank (13301 subjects). Importantly, we find that separating fluctuations in the mean activity levels from those in the FC reveals much stronger changes in FC over time, and is a better predictor of individual behavioural variability.