Statistical model for dynamically-changing correlation matrices with application to brain connectivity

Riemannian manifold-based disentangled representation learning for multi-site functional connectivity analysis

Functional connectivity (FC), derived from resting-state functional magnetic resonance imaging (rs-fMRI), has been widely used to characterize brain abnormalities in disorders. FC is usually defined as a correlation matrix that is a symmetric positive definite (SPD) matrix lying on the Riemannian manifold. Recently, a number of learning-based methods have been proposed for FC analysis, while the geometric properties of Riemannian manifold have not yet been fully explored in previous studies. Also, most existing methods are designed to target one imaging site of fMRI data, which may result in limited training data for learning reliable and robust models. In this paper, we propose a novel Riemannian Manifold-based Disentangled Representation Learning (RM-DRL) framework which is capable of learning invariant representations from fMRI data across multiple sites for brain disorder diagnosis. In RM-DRL, we first employ an SPD-based encoder module to learn a latent unified representation of FC from different sites, which can preserve the Riemannian geometry of the SPD matrices. In latent space, a disentangled representation module is then designed to split the learned features into domain-specific and domain-invariant parts, respectively. Finally, a decoder module is introduced to ensure that sufficient information can be preserved during disentanglement learning. These designs allow us to introduce four types of training objectives to improve the disentanglement learning. Our RM-DRL method is evaluated on the public multi-site ABIDE dataset, showing superior performance compared with several state-of-the-art methods.

Deep learning models reveal the link between dynamic brain connectivity patterns and states of consciousness

Decoding states of consciousness from brain activity is a central challenge in neuroscience. Dynamic functional connectivity (dFC) allows the study of short-term temporal changes in functional connectivity (FC) between distributed brain areas. By clustering dFC matrices from resting-state fMRI, we previously described "brain patterns" that underlie different functional configurations of the brain at rest. The networks associated with these patterns have been extensively analyzed. However, the overall dynamic organization and how it relates to consciousness remains unclear. We hypothesized that deep learning networks would help to model this relationship. Recent studies have used low-dimensional variational autoencoders (VAE) to learn meaningful representations that can help explaining consciousness. Here, we investigated the complexity of selecting such a generative model to study brain dynamics, and extended the available methods for latent space characterization and modeling. Therefore, our contributions are threefold. First, compared with probabilistic principal component analysis and sparse VAE, we showed that the selected low-dimensional VAE exhibits balanced performance in reconstructing dFCs and classifying brain patterns. We then explored the organization of the obtained low-dimensional dFC latent representations. We showed how these representations stratify the dynamic organization of the brain patterns as well as the experimental conditions. Finally, we proposed to delve into the proposed brain computational model. We first applied a receptive field analysis to identify preferred directions in the latent space to move from one brain pattern to another. Then, an ablation study was achieved where we virtually inactivated specific brain areas. We demonstrated the model's efficiency in summarizing consciousness-specific information encoded in key inter-areal connections, as described in the global neuronal workspace theory of consciousness. The proposed framework advocates the possibility of developing an interpretable computational brain model of interest for disorders of consciousness, paving the way for a dynamic diagnostic support tool.

Topological Embedding of Human Brain Networks with Applications to Dynamics of Temporal Lobe Epilepsy

We introduce a novel, data-driven topological data analysis (TDA) approach for embedding brain networks into a lower-dimensional space in quantifying the dynamics of temporal lobe epilepsy (TLE) obtained from resting-state functional magnetic resonance imaging (rs-fMRI). This embedding facilitates the orthogonal projection of 0D and 1D topological features, allowing for the visualization and modeling of the dynamics of functional human brain networks in a resting state. We then quantify the topological disparities between networks to determine the coordinates for embedding. This framework enables us to conduct a coherent statistical inference within the embedded space. Our results indicate that brain network topology in TLE patients exhibits increased rigidity in 0D topology but more rapid flections compared to that of normal controls in 1D topology.

Topological state-space estimation of functional human brain networks

We introduce an innovative, data-driven topological data analysis (TDA) technique for estimating the state spaces of dynamically changing functional human brain networks at rest. Our method utilizes the Wasserstein distance to measure topological differences, enabling the clustering of brain networks into distinct topological states. This technique outperforms the commonly used k-means clustering in identifying brain network state spaces by effectively incorporating the temporal dynamics of the data without the need for explicit model specification. We further investigate the genetic underpinnings of these topological features using a twin study design, examining the heritability of such state changes. Our findings suggest that the topology of brain networks, particularly in their dynamic state changes, may hold significant hidden genetic information.

Unified topological inference for brain networks in temporal lobe epilepsy using the Wasserstein distance

Persistent homology offers a powerful tool for extracting hidden topological signals from brain networks. It captures the evolution of topological structures across multiple scales, known as filtrations, thereby revealing topological features that persist over these scales. These features are summarized in persistence diagrams, and their dissimilarity is quantified using the Wasserstein distance. However, the Wasserstein distance does not follow a known distribution, posing challenges for the application of existing parametric statistical models. To tackle this issue, we introduce a unified topological inference framework centered on the Wasserstein distance. Our approach has no explicit model and distributional assumptions. The inference is performed in a completely data driven fashion. We apply this method to resting-state functional magnetic resonance images (rs-fMRI) of temporal lobe epilepsy patients collected from two different sites: the University of Wisconsin-Madison and the Medical College of Wisconsin. Importantly, our topological method is robust to variations due to sex and image acquisition, obviating the need to account for these variables as nuisance covariates. We successfully localize the brain regions that contribute the most to topological differences. A MATLAB package used for all analyses in this study is available at https://github.com/laplcebeltrami/PH-STAT.

Unified Topological Inference for Brain Networks in Temporal Lobe Epilepsy Using the Wasserstein Distance

Persistent homology offers a powerful tool for extracting hidden topological signals from brain networks. It captures the evolution of topological structures across multiple scales, known as filtrations, thereby revealing topological features that persist over these scales. These features are summarized in persistence diagrams, and their dissimilarity is quantified using the Wasserstein distance. However, the Wasserstein distance does not follow a known distribution, posing challenges for the application of existing parametric statistical models. To tackle this issue, we introduce a unified topological inference framework centered on the Wasserstein distance. Our approach has no explicit model and distributional assumptions. The inference is performed in a completely data driven fashion. We apply this method to resting-state functional magnetic resonance images (rs-fMRI) of temporal lobe epilepsy patients collected from two different sites: the University of Wisconsin-Madison and the Medical College of Wisconsin. Importantly, our topological method is robust to variations due to sex and image acquisition, obviating the need to account for these variables as nuisance covariates. We successfully localize the brain regions that contribute the most to topological differences. A MATLAB package used for all analyses in this study is available at https://github.com/laplcebeltrami/PH-STAT.

Topological analysis of functional connectivity in Parkinson's disease

Parkinson's disease (PD) is a clinically heterogeneous disorder, which mainly affects patients' motor and non-motor function. Functional connectivity was preliminary explored and studied through resting state functional magnetic resonance imaging (rsfMRI). Through the topological analysis of 54 PD scans and 31 age-matched normal controls (NC) in the Neurocon dataset, leveraging on rsfMRI data, the brain functional connection and the Vietoris-Rips (VR) complex were constructed. The barcodes of the complex were calculated to reflect the changes of functional connectivity neural circuits (FCNC) in brain network. The 0-dimensional Betti number β0 means the number of connected branches in VR complex. The average number of connected branches in PD group was greater than that in NC group when the threshold δ ≤ 0.7. Two-sample Mann-Whitney U test and false discovery rate (FDR) correction were used for statistical analysis to investigate the FCNC changes between PD and NC groups. In PD group, under threshold of 0.7, the number of FCNC involved was significantly differences and these brain regions include the Cuneus_R, Lingual_R, Fusiform_R and Heschl_R. There are also significant differences in brain regions in the Frontal_Inf_Orb_R and Pallidum_R, when the threshold increased to 0.8 and 0.9 (p < 0.05). In addition, when the length of FCNC was medium, there was a significant statistical difference between the PD group and the NC group in the Neurocon dataset and the Parkinson's Progression Markers Initiative (PPMI) dataset. Topological analysis based on rsfMRI data may provide comprehensive information about the changes of FCNC and may provide an alternative for clinical differential diagnosis.

Hodge Laplacian of Brain Networks

The closed loops or cycles in a brain network embeds higher order signal transmission paths, which provide fundamental insights into the functioning of the brain. In this work, we propose an efficient algorithm for systematic identification and modeling of cycles using persistent homology and the Hodge Laplacian. Various statistical inference procedures on cycles are developed. We validate the our methods on simulations and apply to brain networks obtained through the resting state functional magnetic resonance imaging. The computer codes for the Hodge Laplacian are given in https://github.com/laplcebeltrami/hodge.

TOPOLOGICAL LEARNING FOR BRAIN NETWORKS

This paper proposes a novel topological learning framework that integrates networks of different sizes and topology through persistent homology. Such challenging task is made possible through the introduction of a computationally efficient topological loss. The use of the proposed loss bypasses the intrinsic computational bottleneck associated with matching networks. We validate the method in extensive statistical simulations to assess its effectiveness when discriminating networks with different topology. The method is further demonstrated in a twin brain imaging study where we determine if brain networks are genetically heritable. The challenge here is due to the difficulty of overlaying the topologically different functional brain networks obtained from resting-state functional MRI onto the template structural brain network obtained through diffusion MRI.

Learning Brain Functional Networks with Latent Temporal Dependency for MCI Identification

AbstractResting-state functional magnetic resonance imaging (rs-fMRI) has become a popular non-invasive way of diagnosing neurological disorders or their early stages by probing functional connectivity between different brain regions of interest (ROIs) across subjects. In the past decades, researchers have proposed many methods to estimate brain functional networks (BFNs) based on blood-oxygen-level-dependent (BOLD) signals captured by rs-fMRI. However, most of the existing methods estimate BFNs under the assumption that signals are independently sampled, which ignores the temporal dependency and sequential order of different time points (or volumes). To address this problem, in this paper, we first propose a novel BFN estimation model by introducing a latent variable to control the sequence of volumes for encoding the temporal dependency and sequential information of signals into the estimated BFNs. Then, we develop an efficient learning algorithm to solve the proposed model by the alternating optimization scheme. To verify the effectiveness of the proposed method, the estimated BFNs are used to identify subjects with mild cognitive impairment (MCIs) from normal controls (NCs). Experimental results show that our method outperforms the baseline methods in the sense of classification performance.