Unsupervised spatiotemporal analysis of fMRI data using graph-based visualizations of self-organizing maps

The degenerate coding of psychometric profiles through functional connectivity archetypes

Introduction

Degeneracy in the brain-behavior code refers to the brain's ability to utilize different neural configurations to support similar functions, reflecting its adaptability and robustness. This study aims to explore degeneracy by investigating the non-linear associations between psychometric profiles and resting-state functional connectivity (RSFC).

Methods

The study analyzed RSFC data from 500 subjects to uncover the underlying neural configurations associated with various psychometric outcomes. Self-organized maps (SOM), a type of unsupervised machine learning algorithm, were employed to cluster the RSFC data. And identify distinct archetypal connectivity profiles characterized by unique within- and between-network connectivity patterns.

Results

The clustering analysis using SOM revealed several distinct archetypal connectivity profiles within the RSFC data. Each archetype exhibited unique connectivity patterns that correlated with various cognitive, physical, and socioemotional outcomes. Notably, the interaction between different SOM dimensions was significantly associated with specific psychometric profiles.

Discussion

This study underscores the complexity of brain-behavior interactions and the brain's capacity for degeneracy, where different neural configurations can lead to similar behavioral outcomes. These findings highlight the existence of multiple brain architectures capable of producing similar behavioral outcomes, illustrating the concept of neural degeneracy, and advance our understanding of neural degeneracy and its implications for cognitive and emotional health.

Brain functional network modeling and analysis based on fMRI: a systematic review

In recent years, the number of patients with neurodegenerative diseases (i.e., Alzheimer's disease, Parkinson's disease, mild cognitive impairment) and mental disorders (i.e., depression, anxiety and schizophrenia) have increased dramatically. Researchers have found that complex network analysis can reveal the topology of brain functional networks, such as small-world, scale-free, etc. In the study of brain diseases, it has been found that these topologies have undergoed abnormal changes in different degrees. Therefore, the research of brain functional networks can not only provide a new perspective for understanding the pathological mechanism of neurological and psychiatric diseases, but also provide assistance for the early diagnosis. Focusing on the study of human brain functional networks, this paper reviews the research results in recent years. First, this paper introduces the background of the study of brain functional networks under complex network theory and the important role of topological properties in the study of brain diseases. Second, the paper describes how to construct a brain functional network using neural image data. Third, the common methods of functional network analysis, including network structure analysis and disease classification, are introduced. Fourth, the role of brain functional networks in pathological study, analysis and diagnosis of brain functional diseases is studied. Finally, the paper summarizes the existing studies of brain functional networks and points out the problems and future research directions.

Density-based clustering of static and dynamic functional MRI connectivity features obtained from subjects with cognitive impairment

Various machine-learning classification techniques have been employed previously to classify brain states in healthy and disease populations using functional magnetic resonance imaging (fMRI). These methods generally use supervised classifiers that are sensitive to outliers and require labeling of training data to generate a predictive model. Density-based clustering, which overcomes these issues, is a popular unsupervised learning approach whose utility for high-dimensional neuroimaging data has not been previously evaluated. Its advantages include insensitivity to outliers and ability to work with unlabeled data. Unlike the popular k-means clustering, the number of clusters need not be specified. In this study, we compare the performance of two popular density-based clustering methods, DBSCAN and OPTICS, in accurately identifying individuals with three stages of cognitive impairment, including Alzheimer’s disease. We used static and dynamic functional connectivity features for clustering, which captures the strength and temporal variation of brain connectivity respectively. To assess the robustness of clustering to noise/outliers, we propose a novel method called recursive-clustering using additive-noise (R-CLAN). Results demonstrated that both clustering algorithms were effective, although OPTICS with dynamic connectivity features outperformed in terms of cluster purity (95.46%) and robustness to noise/outliers. This study demonstrates that density-based clustering can accurately and robustly identify diagnostic classes in an unsupervised way using brain connectivity.

A decomposition model and voxel selection framework for fMRI analysis to predict neural response of visual stimuli

This paper presents an fMRI signal analysis methodology using geometric mean curve decomposition (GMCD) and mutual information-based voxel selection framework. Previously, the fMRI signal analysis has been conducted using empirical mean curve decomposition (EMCD) model and voxel selection on raw fMRI signal. The erstwhile methodology loses frequency component, while the latter methodology suffers from signal redundancy. Both challenges are addressed by our methodology in which the frequency component is considered by decomposing the raw fMRI signal using geometric mean rather than arithmetic mean and the voxels are selected from EMCD signal using GMCD components, rather than raw fMRI signal. The proposed methodologies are adopted for predicting the neural response. Experimentations are conducted in the openly available fMRI data of six subjects, and comparisons are made with existing decomposition models and voxel selection frameworks. Subsequently, the effect of degree of selected voxels and the selection constraints are analyzed. The comparative results and the analysis demonstrate the superiority and the reliability of the proposed methodology.

Variance decomposition for single-subject task-based fMRI activity estimates across many sessions

Here we report an exploratory within-subject variance decomposition analysis conducted on a task-based fMRI dataset with an unusually large number of repeated measures (i.e., 500 trials in each of three different subjects) distributed across 100 functional scans and 9 to 10 different sessions. Within-subject variance was segregated into four primary components: variance across-sessions, variance across-runs within a session, variance across-blocks within a run, and residual measurement/modeling error. Our results reveal inhomogeneous and distinct spatial distributions of these variance components across significantly active voxels in grey matter. Measurement error is dominant across the whole brain. Detailed evaluation of the remaining three components shows that across-session variance is the second largest contributor to total variance in occipital cortex, while across-runs variance is the second dominant source for the rest of the brain. Network-specific analysis revealed that across-block variance contributes more to total variance in higher-order cognitive networks than in somatosensory cortex. Moreover, in some higher-order cognitive networks across-block variance can exceed across-session variance. These results help us better understand the temporal (i.e., across blocks, runs and sessions) and spatial distributions (i.e., across different networks) of within-subject natural variability in estimates of task responses in fMRI. They also suggest that different brain regions will show different natural levels of test-retest reliability even in the absence of residual artifacts and sufficiently high contrast-to-noise measurements. Further confirmation with a larger sample of subjects and other tasks is necessary to ensure generality of these results.

Experimental Validation of Dynamic Granger Causality for Inferring Stimulus-Evoked Sub-100 ms Timing Differences from fMRI

Decoding the sequential flow of events in the human brain non-invasively is critical for gaining a mechanistic understanding of brain function. In this study, we propose a method based on dynamic Granger causality analysis to measure timing differences in brain responses from fMRI. We experimentally validate this method by detecting sub-100 ms timing differences in fMRI responses obtained from bilateral visual cortex using fast sampling, ultra-high field and an event-related visual hemifield paradigm with known timing difference between the hemifields. Classical Granger causality was previously shown to be able to detect sub-100 ms timing differences in the visual cortex. Since classical Granger causality does not differentiate between spontaneous and stimulus-evoked responses, dynamic Granger causality has been proposed as an alternative, thereby necessitating its experimental validation. In addition to detecting timing differences as low as 28 ms using dynamic Granger causality, the significance of the inference from our method increased with increasing delay both in simulations and experimental data. Therefore, it provides a methodology for understanding mental chronometry from fMRI in a data-driven way.

Learning from data: recognizing glaucomatous defect patterns and detecting progression from visual field measurements

A hierarchical approach to learn from visual field data was adopted to identify glaucomatous visual field defect patterns and to detect glaucomatous progression. The analysis pipeline included three stages, namely, clustering, glaucoma boundary limit detection, and glaucoma progression detection testing. First, cross-sectional visual field tests collected from each subject were clustered using a mixture of Gaussians and model parameters were estimated using expectation maximization. The visual field clusters were further estimated to recognize glaucomatous visual field defect patterns by decomposing each cluster into several axes. The glaucoma visual field defect patterns along each axis then were identified. To derive a definition of progression, the longitudinal visual fields of stable glaucoma eyes on the abnormal cluster axes were projected and the slope was approximated using linear regression (LR) to determine the confidence limit of each axis. For glaucoma progression detection, the longitudinal visual fields of each eye on the abnormal cluster axes were projected and the slope was approximated by LR. Progression was assigned if the progression rate was greater than the boundary limit of the stable eyes; otherwise, stability was assumed. The proposed method was compared to a recently developed progression detection method and to clinically available glaucoma progression detection software. The clinical accuracy of the proposed pipeline was as good as or better than the currently available methods.

Recognizing patterns of visual field loss using unsupervised machine learning

Glaucoma is a potentially blinding optic neuropathy that results in a decrease in visual sensitivity. Visual field abnormalities (decreased visual sensitivity on psychophysical tests) are the primary means of glaucoma diagnosis. One form of visual field testing is Frequency Doubling Technology (FDT) that tests sensitivity at 52 points within the visual field. Like other psychophysical tests used in clinical practice, FDT results yield specific patterns of defect indicative of the disease. We used Gaussian Mixture Model with Expectation Maximization (GEM), (EM is used to estimate the model parameters) to automatically separate FDT data into clusters of normal and abnormal eyes. Principal component analysis (PCA) was used to decompose each cluster into different axes (patterns). FDT measurements were obtained from 1,190 eyes with normal FDT results and 786 eyes with abnormal (i.e., glaucomatous) FDT results, recruited from a university-based, longitudinal, multi-center, clinical study on glaucoma. The GEM input was the 52-point FDT threshold sensitivities for all eyes. The optimal GEM model separated the FDT fields into 3 clusters. Cluster 1 contained 94% normal fields (94% specificity) and clusters 2 and 3 combined, contained 77% abnormal fields (77% sensitivity). For clusters 1, 2 and 3 the optimal number of PCA-identified axes were 2, 2 and 5, respectively. GEM with PCA successfully separated FDT fields from healthy and glaucoma eyes and identified familiar glaucomatous patterns of loss.