Discovering common change-point patterns in functional connectivity across subjects

Uncovering shape signatures of resting-state functional connectivity by geometric deep learning on Riemannian manifold

Functional neural activities manifest geometric patterns, as evidenced by the evolving network topology of functional connectivities (FC) even in the resting state. In this work, we propose a novel manifold‐based geometric neural network for functional brain networks (called “Geo‐Net4Net” for short) to learn the intrinsic low‐dimensional feature representations of resting‐state brain networks on the Riemannian manifold. This tool allows us to answer the scientific question of how the spontaneous fluctuation of FC supports behavior and cognition. We deploy a set of positive maps and rectified linear unit (ReLU) layers to uncover the intrinsic low‐dimensional feature representations of functional brain networks on the Riemannian manifold taking advantage of the symmetric positive‐definite (SPD) form of the correlation matrices. Due to the lack of well‐defined ground truth in the resting state, existing learning‐based methods are limited to unsupervised methodologies. To go beyond this boundary, we propose to self‐supervise the feature representation learning of resting‐state functional networks by leveraging the task‐based counterparts occurring before and after the underlying resting state. With this extra heuristic, our Geo‐Net4Net allows us to establish a more reasonable understanding of resting‐state FCs by capturing the geometric patterns (aka. spectral/shape signature) associated with resting states on the Riemannian manifold. We have conducted extensive experiments on both simulated data and task‐based functional resonance magnetic imaging (fMRI) data from the Human Connectome Project (HCP) database, where our Geo‐Net4Net not only achieves more accurate change detection results than other state‐of‐the‐art counterpart methods but also yields ubiquitous geometric patterns that manifest putative insights into brain function.

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.

Cross-covariance isolate detect: A new change-point method for estimating dynamic functional connectivity

Evidence of the non stationary behavior of functional connectivity (FC) networks has been observed in task based functional magnetic resonance imaging (fMRI) experiments and even prominently in resting state fMRI data. This has led to the development of several new statistical methods for estimating this time-varying connectivity, with the majority of the methods utilizing a sliding window approach. While computationally feasible, the sliding window approach has several limitations. In this paper, we circumvent the sliding window, by introducing a statistical method that finds change-points in FC networks where the number and location of change-points are unknown a priori. The new method, called cross-covariance isolate detect (CCID), detects multiple change-points in the second-order (cross-covariance or network) structure of multivariate, possibly high-dimensional time series. CCID allows for change-point detection in the presence of frequent changes of possibly small magnitudes, can assign change-points to one or multiple brain regions, and is computationally fast. In addition, CCID is particularly suited to task based data, where the subject alternates between task and rest, as it firstly attempts isolation of each of the change-points within subintervals, and secondly their detection therein. Furthermore, we also propose a new information criterion for CCID to identify the change-points. We apply CCID to several simulated data sets and to task based and resting state fMRI data and compare it to recent change-point methods. CCID may also be applicable to electroencephalography (EEG), magentoencephalography (MEG) and electrocorticography (ECoG) data. Similar to other biological networks, understanding the complex network organization and functional dynamics of the brain can lead to profound clinical implications. Finally, the R package ccid implementing the method from the paper is available from CRAN.

Dual stream neural networks for brain signal classification

Objective. The primary objective of this work is to develop a neural nework classifier for arbitrary collections of functional neuroimaging signals to be used in brain–computer interfaces (BCIs). Approach. We propose a dual stream neural network (DSNN) for the classification problem. The first stream is an end-to-end classifier taking raw time-dependent signals as input and generating feature identification signatures from them. The second stream enhances the identified features from the first stream by adjoining a dynamic functional connectivity matrix aimed at incorporating nuanced multi-channel information during specified BCI tasks. Main results. The proposed DSNN classifier is benchmarked against three publicly available datasets, where the classifier demonstrates performance comparable to, or better than the state-of-art in each instance. An information theoretic examination of the trained network is also performed, utilizing various tools, to demonstrate how to glean interpretive insight into how the hidden layers of the network parse the underlying biological signals. Significance . The resulting DSNN is a subject-independent classifier that works for any collection of 1D functional neuroimaging signals, with the option of integrating domain specific information in the design.

Learning Common Harmonic Waves on Stiefel Manifold - A New Mathematical Approach for Brain Network Analyses

Converging evidence shows that disease-relevant brain alterations do not appear in random brain locations, instead, their spatial patterns follow large-scale brain networks. In this context, a powerful network analysis approach with a mathematical foundation is indispensable to understand the mechanisms of neuropathological events as they spread through the brain. Indeed, the topology of each brain network is governed by its native harmonic waves, which are a set of orthogonal bases derived from the Eigen-system of the underlying Laplacian matrix. To that end, we propose a novel connectome harmonic analysis framework that provides enhanced mathematical insights by detecting frequency-based alterations relevant to brain disorders. The backbone of our framework is a novel manifold algebra appropriate for inference across harmonic waves. This algebra overcomes the limitations of using classic Euclidean operations on irregular data structures. The individual harmonic differences are measured by a set of common harmonic waves learned from a population of individual Eigen-systems, where each native Eigen-system is regarded as a sample drawn from the Stiefel manifold. Specifically, a manifold optimization scheme is tailored to find the common harmonic waves, which reside at the center of the Stiefel manifold. To that end, the common harmonic waves constitute a new set of neurobiological bases to understand disease progression. Each harmonic wave exhibits a unique propagation pattern of neuropathological burden spreading across brain networks. The statistical power of our novel connectome harmonic analysis approach is evaluated by identifying frequency-based alterations relevant to Alzheimer's disease, where our learning-based manifold approach discovers more significant and reproducible network dysfunction patterns than Euclidean methods.