Critical dynamics and interictal epileptiform discharge: a comparative analysis with respect to tracking seizure risk cycles.

Epilepsy is characterized by recurrent, unprovoked seizures. Accurate prediction of seizure occurrence has long been a clinical goal since this would allow to optimize patient treatment, prevent injuries due to seizures, and alleviate the patient burden of unpredictability. Advances in implantable electroencephalographic (EEG) devices, allowing for long-term interictal EEG recordings, have facilitated major progress in this field. Recently, it has been discovered that interictal brain activity demonstrates circadian and multi-dien cycles that are strongly aligned, or phase locked, with seizure risk. Thus, cyclical brain activity patterns have been used to forecast seizures. However, in the effort to develop a clinically useful EEG based seizure forecasting system, challenges remain. Firstly, multiple EEG features demonstrate cyclical patterns, but it remains unclear which feature is best suited for predicting seizures. Secondly, the technology for long-term EEG recording is currently limited in both spatial and temporal sampling resolution. In this study, we compare five established EEG metrics:synchrony, spatial correlation, temporal correlation, signal variance which have been motivated from critical dynamics theory, and interictal epileptiform discharge (IED) which are a traditional marker of seizure propensity. We assess their effectiveness in detecting 24-h and seizure cycles as well as their robustness under spatial and temporal subsampling. Analyzing intracranial EEG data from 23 patients, we report that all examined features exhibit 24-h cycles. Spatial correlation, signal variance, and synchrony showed the highest phase locking with seizures, while IED rates were the lowest. Notably, spatial and temporal correlation were also found to be highly correlated to each other, as were signal variance and IED-suggesting some features may reflect similar aspects of cortical dynamics, whereas others provide complementary information. All features proved robust under subsampling, indicating that the dynamic properties of interictal activity evolve slowly and are not confined to specific brain regions. Our results may aid future translational research by assisting in design and testing of EEG based seizure forecasting systems.

A deep learning approach to estimating initial conditions of Brain Network Models in reference to measured fMRI data.

Introduction: Brain Network Models (BNMs) are mathematical models that simulate the activity of the entire brain. These models use neural mass models to represent local activity in different brain regions that interact with each other via a global structural network. Researchers have been interested in using these models to explain measured brain activity, particularly resting state functional magnetic resonance imaging (rs-fMRI). BNMs have shown to produce similar properties as measured data computed over longer periods of time such as average functional connectivity (FC), but it is unclear how well simulated trajectories compare to empirical trajectories on a timepoint-by-timepoint basis. During task fMRI, the relevant processes pertaining to task occur over the time frame of the hemodynamic response function, and thus it is important to understand how BNMs capture these dynamics over these short periods. Methods: To test the nature of BNMs' short-term trajectories, we used a deep learning technique called Neural ODE to simulate short trajectories from estimated initial conditions based on observed fMRI measurements. To compare to previous methods, we solved for the parameterization of a specific BNM, the Firing Rate Model, using these short-term trajectories as a metric. Results: Our results show an agreement between parameterization of using previous long-term metrics with the novel short term metrics exists if also considering other factors such as the sensitivity in accuracy with relative to changes in structural connectivity, and the presence of noise. Discussion: Therefore, we conclude that there is evidence that by using Neural ODE, BNMs can be simulated in a meaningful way when comparing against measured data trajectories, although future studies are necessary to establish how BNM activity relate to behavioral variables or to faster neural processes during this time period.

Comparison of Resting-State Functional MRI Methods for Characterizing Brain Dynamics.

Resting-state functional MRI (fMRI) exhibits time-varying patterns of functional connectivity. Several different analysis approaches have been developed for examining these resting-state dynamics including sliding window connectivity (SWC), phase synchrony (PS), co-activation pattern (CAP), and quasi-periodic patterns (QPP). Each of these approaches can be used to generate patterns of activity or inter-areal coordination which vary across time. The individual frames can then be clustered to produce temporal groupings commonly referred to as "brain states." Several recent publications have investigated brain state alterations in clinical populations, typically using a single method for quantifying frame-wise functional connectivity. This study directly compares the results of k-means clustering in conjunction with three of these resting-state dynamics methods (SWC, CAP, and PS) and quantifies the brain state dynamics across several metrics using high resolution data from the human connectome project. Additionally, these three dynamics methods are compared by examining how the brain state characterizations vary during the repeated sequences of brain states identified by a fourth dynamic analysis method, QPP. The results indicate that the SWC, PS, and CAP methods differ in the clusters and trajectories they produce. A clear illustration of these differences is given by how each one results in a very different clustering profile for the 24s sequences explicitly identified by the QPP algorithm. PS clustering is sensitive to QPPs with the mid-point of most QPP sequences grouped into the same single cluster. CAPs are also highly sensitive to QPPs, separating each phase of the QPP sequences into different sets of clusters. SWC (60s window) is less sensitive to QPPs. While the QPPs are slightly more likely to occur during specific SWC clusters, the SWC clustering does not vary during the 24s QPP sequences, the goal of this work is to improve both the practical and theoretical understanding of different resting-state dynamics methods, thereby enabling investigators to better conceptualize and implement these tools for characterizing functional brain networks.

Resting Brain Fluctuations Are Intrinsically Coupled to Visual Response Dynamics.

How do intrinsic brain dynamics interact with processing of external sensory stimuli? We sought new insights using functional magnetic resonance imaging to track spatiotemporal activity patterns at the whole brain level in lightly anesthetized mice, during both resting conditions and visual stimulation trials. Our results provide evidence that quasiperiodic patterns (QPPs) are the most prominent component of mouse resting brain dynamics. These QPPs captured the temporal alignment of anticorrelation between the default mode (DMN)- and task-positive (TPN)-like networks, with global brain fluctuations, and activity in neuromodulatory nuclei of the reticular formation. Specifically, the phase of QPPs prior to stimulation could significantly stratify subsequent visual response magnitude, suggesting QPPs relate to brain state fluctuations. This is the first observation in mice that dynamics of the DMN- and TPN-like networks, and particularly their anticorrelation, capture a brain state dynamic that affects sensory processing. Interestingly, QPPs also displayed transient onset response properties during visual stimulation, which covaried with deactivations in the reticular formation. We conclude that QPPs appear to capture a brain state fluctuation that may be orchestrated through neuromodulation. Our findings provide new frontiers to understand the neural processes that shape functional brain states and modulate sensory input processing.

Brain network constraints and recurrent neural networks reproduce unique trajectories and state transitions seen over the span of minutes in resting-state fMRI.

Large-scale patterns of spontaneous whole-brain activity seen in resting-state functional magnetic resonance imaging (rs-fMRI) are in part believed to arise from neural populations interacting through the structural network (Honey, Kötter, Breakspear, & Sporns, 2007). Generative models that simulate this network activity, called brain network models (BNM), are able to reproduce global averaged properties of empirical rs-fMRI activity such as functional connectivity (FC) but perform poorly in reproducing unique trajectories and state transitions that are observed over the span of minutes in whole-brain data (Cabral, Kringelbach, & Deco, 2017; Kashyap & Keilholz, 2019). The manuscript demonstrates that by using recurrent neural networks, it can fit the BNM in a novel way to the rs-fMRI data and predict large amounts of variance between subsequent measures of rs-fMRI data. Simulated data also contain unique repeating trajectories observed in rs-fMRI, called quasiperiodic patterns (QPP), that span 20 s and complex state transitions observed using k-means analysis on windowed FC matrices (Allen et al., 2012; Majeed et al., 2011). Our approach is able to estimate the manifold of rs-fMRI dynamics by training on generating subsequent time points, and it can simulate complex resting-state trajectories better than the traditional generative approaches.

Author Correction: Quasi-Periodic Patterns of Neural Activity improve Classification of Alzheimer's Disease in Mice.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Dynamic properties of simulated brain network models and empirical resting-state data.

Brain network models (BNMs) have become a promising theoretical framework for simulating signals that are representative of whole-brain activity such as resting-state fMRI. However, it has been difficult to compare the complex brain activity obtained from simulations to empirical data. Previous studies have used simple metrics to characterize coordination between regions such as functional connectivity. We extend this by applying various different dynamic analysis tools that are currently used to understand empirical resting-state fMRI (rs-fMRI) to the simulated data. We show that certain properties correspond to the structural connectivity input that is shared between the models, and certain dynamic properties relate more to the mathematical description of the brain network model. We conclude that the dynamic properties that explicitly examine patterns of signal as a function of time rather than spatial coordination between different brain regions in the rs-fMRI signal seem to provide the largest contrasts between different BNMs and the unknown empirical dynamical system. Our results will be useful in constraining and developing more realistic simulations of whole-brain activity.

Quasi-periodic patterns contribute to functional connectivity in the brain.

Functional connectivity is widely used to study the coordination of activity between brain regions over time. Functional connectivity in the default mode and task positive networks is particularly important for normal brain function. However, the processes that give rise to functional connectivity in the brain are not fully understood. It has been postulated that low-frequency neural activity plays a key role in establishing the functional architecture of the brain. Quasi-periodic patterns (QPPs) are a reliably observable form of low-frequency neural activity that involve the default mode and task positive networks. Here, QPPs from resting-state and working memory task-performing individuals were acquired. The spatiotemporal pattern, strength, and frequency of the QPPs between the two groups were compared and the contribution of QPPs to functional connectivity in the brain was measured. In task-performing individuals, the spatiotemporal pattern of the QPP changes, particularly in task-relevant regions, and the QPP tends to occur with greater strength and frequency. Differences in the QPPs between the two groups could partially account for the variance in functional connectivity between resting-state and task-performing individuals. The QPPs contribute strongly to connectivity in the default mode and task positive networks and to the strength of anti-correlation seen between the two networks. Many of the connections affected by QPPs are also disrupted during several neurological disorders. These findings contribute to understanding the dynamic neural processes that give rise to functional connectivity in the brain and how they may be disrupted during disease.

Quasi-Periodic Patterns of Neural Activity improve Classification of Alzheimer's Disease in Mice.

Resting state (rs)fMRI allows measurement of brain functional connectivity and has identified default mode (DMN) and task positive (TPN) network disruptions as promising biomarkers for Alzheimer's disease (AD). Quasi-periodic patterns (QPPs) of neural activity describe recurring spatiotemporal patterns that display DMN with TPN anti-correlation. We reasoned that QPPs could provide new insights into AD network dysfunction and improve disease diagnosis. We therefore used rsfMRI to investigate QPPs in old TG2576 mice, a model of amyloidosis, and age-matched controls. Multiple QPPs were determined and compared across groups. Using linear regression, we removed their contribution from the functional scans and assessed how they reflected functional connectivity. Lastly, we used elastic net regression to determine if QPPs improved disease classification. We present three prominent findings: (1) Compared to controls, TG2576 mice were marked by opposing neural dynamics in which DMN areas were anti-correlated and displayed diminished anti-correlation with the TPN. (2) QPPs reflected lowered DMN functional connectivity in TG2576 mice and revealed significantly decreased DMN-TPN anti-correlations. (3) QPP-derived measures significantly improved classification compared to conventional functional connectivity measures. Altogether, our findings provide insight into the neural dynamics of aberrant network connectivity in AD and indicate that QPPs might serve as a translational diagnostic tool.

Instantaneous brain dynamics mapped to a continuous state space.

Measures of whole-brain activity, from techniques such as functional Magnetic Resonance Imaging, provide a means to observe the brain's dynamical operations. However, interpretation of whole-brain dynamics has been stymied by the inherently high-dimensional structure of brain activity. The present research addresses this challenge through a series of scale transformations in the spectral, spatial, and relational domains. Instantaneous multispectral dynamics are first developed from input data via a wavelet filter bank. Voxel-level signals are then projected onto a representative set of spatially independent components. The correlation distance over the instantaneous wavelet-ICA state vectors is a graph that may be embedded onto a lower-dimensional space to assist the interpretation of state-space dynamics. Applying this procedure to a large sample of resting-state and task-active data (acquired through the Human Connectome Project), we segment the empirical state space into a continuum of stimulus-dependent brain states. Upon observing the local neighborhood of brain-states adopted subsequent to each stimulus, we may conclude that resting brain activity includes brain states that are, at times, similar to those adopted during tasks, but that are at other times distinct from task-active brain states. As task-active brain states often populate a local neighborhood, back-projection of segments of the dynamical state space onto the brain's surface reveals the patterns of brain activity that support many experimentally-defined states.