-
Topological Analysis of Seizure-Induced Changes in Brain Hierarchy Through Effective Connectivity
Authors:
Anass B. El-Yaagoubi,
Moo K. Chung,
Hernando Ombao
Abstract:
Traditional Topological Data Analysis (TDA) methods, such as Persistent Homology (PH), rely on distance measures (e.g., cross-correlation, partial correlation, coherence, and partial coherence) that are symmetric by definition. While useful for studying topological patterns in functional brain connectivity, the main limitation of these methods is their inability to capture the directional dynamics…
▽ More
Traditional Topological Data Analysis (TDA) methods, such as Persistent Homology (PH), rely on distance measures (e.g., cross-correlation, partial correlation, coherence, and partial coherence) that are symmetric by definition. While useful for studying topological patterns in functional brain connectivity, the main limitation of these methods is their inability to capture the directional dynamics - which is crucial for understanding effective brain connectivity. We propose the Causality-Based Topological Ranking (CBTR) method, which integrates Causal Inference (CI) to assess effective brain connectivity with Hodge Decomposition (HD) to rank brain regions based on their mutual influence. Our simulations confirm that the CBTR method accurately and consistently identifies hierarchical structures in multivariate time series data. Moreover, this method effectively identifies brain regions showing the most significant interaction changes with other regions during seizures using electroencephalogram (EEG) data. These results provide novel insights into the brain's hierarchical organization and illuminate the impact of seizures on its dynamics.
△ Less
Submitted 18 July, 2024;
originally announced July 2024.
-
Topological Embedding of Human Brain Networks with Applications to Dynamics of Temporal Lobe Epilepsy
Authors:
Moo K. Chung,
Ji Bi Che,
Veena A. Nair,
Camille Garcia Ramos,
Jedidiah Ray Mathis,
Vivek Prabhakaran,
Elizabeth Meyerand,
Bruce P. Hermann,
Jeffrey R. Binder,
Aaron F. Struck
Abstract:
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 model…
▽ More
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.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Spectral Topological Data Analysis of Brain Signals
Authors:
Anass B. El-Yaagoubi,
Shuhao Jiao,
Moo K. Chung,
Hernando Ombao
Abstract:
Topological data analysis (TDA) has become a powerful approach over the last twenty years, mainly due to its ability to capture the shape and the geometry inherent in the data. Persistence homology, which is a particular tool in TDA, has been demonstrated to be successful in analyzing functional brain connectivity. One limitation of standard approaches is that they use arbitrarily chosen threshold…
▽ More
Topological data analysis (TDA) has become a powerful approach over the last twenty years, mainly due to its ability to capture the shape and the geometry inherent in the data. Persistence homology, which is a particular tool in TDA, has been demonstrated to be successful in analyzing functional brain connectivity. One limitation of standard approaches is that they use arbitrarily chosen threshold values for analyzing connectivity matrices. To overcome this weakness, TDA provides a filtration of the weighted brain network across a range of threshold values. However, current analyses of the topological structure of functional brain connectivity primarily rely on overly simplistic connectivity measures, such as the Pearson orrelation. These measures do not provide information about the specific oscillators that drive dependence within the brain network. Here, we develop a frequency-specific approach that utilizes coherence, a measure of dependence in the spectral domain, to evaluate the functional connectivity of the brain. Our approach, the spectral TDA (STDA), has the ability to capture more nuanced and detailed information about the underlying brain networks. The proposed STDA method leads to a novel topological summary, the spectral landscape, which is a 2D-generalization of the persistence landscape. Using the novel spectral landscape, we analyze the EEG brain connectivity of patients with attention deficit hyperactivity disorder (ADHD) and shed light on the frequency-specific differences in the topology of brain connectivity between the controls and ADHD patients.
△ Less
Submitted 1 December, 2023;
originally announced January 2024.
-
Sulcal Pattern Matching with the Wasserstein Distance
Authors:
Zijian Chen,
Soumya Das,
Moo K. Chung
Abstract:
We present the unified computational framework for modeling the sulcal patterns of human brain obtained from the magnetic resonance images. The Wasserstein distance is used to align the sulcal patterns nonlinearly. These patterns are topologically different across subjects making the pattern matching a challenge. We work out the mathematical details and develop the gradient descent algorithms for…
▽ More
We present the unified computational framework for modeling the sulcal patterns of human brain obtained from the magnetic resonance images. The Wasserstein distance is used to align the sulcal patterns nonlinearly. These patterns are topologically different across subjects making the pattern matching a challenge. We work out the mathematical details and develop the gradient descent algorithms for estimating the deformation field. We further quantify the image registration performance. This method is applied in identifying the differences between male and female sulcal patterns.
△ Less
Submitted 1 July, 2023;
originally announced July 2023.
-
Altered Topological Structure of the Brain White Matter in Maltreated Children through Topological Data Analysis
Authors:
Moo K. Chung,
Tahmineh Azizi,
Jamie L. Hanson,
Andrew L. Alexander,
Richard J. Davidson,
Seth D. Pollak
Abstract:
Childhood maltreatment may adversely affect brain development and consequently influence behavioral, emotional, and psychological patterns during adulthood. In this study, we propose an analytical pipeline for modeling the altered topological structure of brain white matter in maltreated and typically developing children. We perform topological data analysis (TDA) to assess the alteration in the g…
▽ More
Childhood maltreatment may adversely affect brain development and consequently influence behavioral, emotional, and psychological patterns during adulthood. In this study, we propose an analytical pipeline for modeling the altered topological structure of brain white matter in maltreated and typically developing children. We perform topological data analysis (TDA) to assess the alteration in the global topology of the brain white-matter structural covariance network among children. We use persistent homology, an algebraic technique in TDA, to analyze topological features in the brain covariance networks constructed from structural magnetic resonance imaging (MRI) and diffusion tensor imaging (DTI). We develop a novel framework for statistical inference based on the Wasserstein distance to assess the significance of the observed topological differences. Using these methods in comparing maltreated children to a typically developing control group, we find that maltreatment may increase homogeneity in white matter structures and thus induce higher correlations in the structural covariance; this is reflected in the topological profile. Our findings strongly suggest that TDA can be a valuable framework to model altered topological structures of the brain. The MATLAB codes and processed data used in this study can be found at https://github.com/laplcebeltrami/maltreated.
△ Less
Submitted 14 November, 2023; v1 submitted 12 April, 2023;
originally announced April 2023.
-
Unified Topological Inference for Brain Networks in Temporal Lobe Epilepsy Using the Wasserstein Distance
Authors:
Moo K. Chung,
Camille Garcia Ramos,
Felipe Branco De Paiva,
Jedidiah Mathis,
Vivek Prabharakaren,
Veena A. Nair,
Elizabeth Meyerand,
Bruce P. Hermann,
Jeffrey R. Binder,
Aaron F. Struck
Abstract:
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.…
▽ More
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.
△ Less
Submitted 20 September, 2023; v1 submitted 13 February, 2023;
originally announced February 2023.
-
Hodge-Decomposition of Brain Networks
Authors:
D. Vijay Anand,
Moo K. Chung
Abstract:
We analyze brain networks by decomposing them into three orthogonal components: gradient, curl, and harmonic flows, through the Hodge decomposition, a technique advantageous for capturing complex topological features. A Wasserstein distance based topological inference is developed to determine the statistical significance of each component. The Hodge decomposition is applied to human brain network…
▽ More
We analyze brain networks by decomposing them into three orthogonal components: gradient, curl, and harmonic flows, through the Hodge decomposition, a technique advantageous for capturing complex topological features. A Wasserstein distance based topological inference is developed to determine the statistical significance of each component. The Hodge decomposition is applied to human brain networks obtained from a resting-state fMRI study. Our results indicate statistically significant differences in the topological features between male and female brain networks.
△ Less
Submitted 1 April, 2024; v1 submitted 18 November, 2022;
originally announced November 2022.
-
Dynamic Topological Data Analysis of Functional Human Brain Networks
Authors:
Moo K. Chung,
Soumya Das,
Hernando Ombao
Abstract:
Developing reliable methods to discriminate different transient brain states that change over time is a key neuroscientific challenge in brain imaging studies. Topological data analysis (TDA), a novel framework based on algebraic topology, can handle such a challenge. However, existing TDA has been somewhat limited to capturing the static summary of dynamically changing brain networks. We propose…
▽ More
Developing reliable methods to discriminate different transient brain states that change over time is a key neuroscientific challenge in brain imaging studies. Topological data analysis (TDA), a novel framework based on algebraic topology, can handle such a challenge. However, existing TDA has been somewhat limited to capturing the static summary of dynamically changing brain networks. We propose a novel dynamic-TDA framework that builds persistent homology over a time series of brain networks. We construct a Wasserstein distance based inference procedure to discriminate between time series of networks. The method is applied to the resting-state functional magnetic resonance images of human brain. We demonstrate that our proposed dynamic-TDA approach can distinctly discriminate between the topological patterns of male and female brain networks. MATLAB code for implementing this method is available at https://github.com/laplcebeltrami/PH-STAT.
△ Less
Submitted 18 December, 2023; v1 submitted 17 October, 2022;
originally announced October 2022.
-
Topological Data Analysis for Multivariate Time Series Data
Authors:
Anass El Yaagoubi Bourakna,
Moo K. Chung,
Hernando Ombao
Abstract:
Over the last two decades, topological data analysis (TDA) has emerged as a very powerful data analytic approach which can deal with various data modalities of varying complexities. One of the most commonly used tools in TDA is persistent homology (PH) which can extract topological properties from data at various scales. Our aim in this article is to introduce TDA concepts to a statistical audienc…
▽ More
Over the last two decades, topological data analysis (TDA) has emerged as a very powerful data analytic approach which can deal with various data modalities of varying complexities. One of the most commonly used tools in TDA is persistent homology (PH) which can extract topological properties from data at various scales. Our aim in this article is to introduce TDA concepts to a statistical audience and provide an approach to analyze multivariate time series data. The application focus will be on multivariate brain signals and brain connectivity networks. Finally, the paper concludes with an overview of some open problems and potential application of TDA to modeling directionality in a brain network as well as the casting of TDA in the context of mixed effects models to capture variations in the topological properties of data collected from multiple subjects
△ Less
Submitted 28 April, 2022;
originally announced April 2022.
-
Embedding of Functional Human Brain Networks on a Sphere
Authors:
Moo K. Chung,
Zijian Chen
Abstract:
Human brain activity is often measured using the blood-oxygen-level dependent (BOLD) signals obtained through functional magnetic resonance imaging (fMRI). The strength of connectivity between brain regions is then measured as a Pearson correlation matrix. As the number of brain regions increases, the dimension of matrix increases. It becomes extremely cumbersome to even visualize and quantify suc…
▽ More
Human brain activity is often measured using the blood-oxygen-level dependent (BOLD) signals obtained through functional magnetic resonance imaging (fMRI). The strength of connectivity between brain regions is then measured as a Pearson correlation matrix. As the number of brain regions increases, the dimension of matrix increases. It becomes extremely cumbersome to even visualize and quantify such weighted complete networks. To remedy the problem, we propose to embed brain networks onto a sphere, which is a Riemannian manifold with constant positive curvature. The Matlab code for the spherical embedding is given in https://github.com/laplcebeltrami/sphericalMDS.
△ Less
Submitted 19 May, 2022; v1 submitted 7 April, 2022;
originally announced April 2022.
-
Topological Data Analysis of Human Brain Networks Through Order Statistics
Authors:
Soumya Das,
D. Vijay Anand,
Moo K. Chung
Abstract:
Understanding the common topological characteristics of the human brain network across a population is central to understanding brain functions. The abstraction of human connectome as a graph has been pivotal in gaining insights on the topological properties of the brain network. The development of group-level statistical inference procedures in brain graphs while accounting for the heterogeneity…
▽ More
Understanding the common topological characteristics of the human brain network across a population is central to understanding brain functions. The abstraction of human connectome as a graph has been pivotal in gaining insights on the topological properties of the brain network. The development of group-level statistical inference procedures in brain graphs while accounting for the heterogeneity and randomness still remains a difficult task. In this study, we develop a robust statistical framework based on persistent homology using the order statistics for analyzing brain networks. The use of order statistics greatly simplifies the computation of the persistent barcodes. We validate the proposed methods using comprehensive simulation studies and subsequently apply to the resting-state functional magnetic resonance images. We found a statistically significant topological difference between the male and female brain networks.
△ Less
Submitted 13 October, 2022; v1 submitted 5 April, 2022;
originally announced April 2022.
-
Statistical Analysis on Brain Surfaces
Authors:
Moo K. Chung,
Jamie L. Hanson,
Seth D. Pollak
Abstract:
In this paper, we review widely used statistical analysis frameworks for data defined along cortical and subcortical surfaces that have been developed in last two decades. The cerebral cortex has the topology of a 2D highly convoluted sheet. For data obtained along curved non-Euclidean surfaces, traditional statistical analysis and smoothing techniques based on the Euclidean metric structure are i…
▽ More
In this paper, we review widely used statistical analysis frameworks for data defined along cortical and subcortical surfaces that have been developed in last two decades. The cerebral cortex has the topology of a 2D highly convoluted sheet. For data obtained along curved non-Euclidean surfaces, traditional statistical analysis and smoothing techniques based on the Euclidean metric structure are inefficient. To increase the signal-to-noise ratio (SNR) and to boost the sensitivity of the analysis, it is necessary to smooth out noisy surface data. However, this requires smoothing data on curved cortical manifolds and assigning smoothing weights based on the geodesic distance along the surface. Thus, many cortical surface data analysis frameworks are differential geometric in nature. The smoothed surface data is then treated as smooth random fields and statistical inferences can be performed within Keith Worsley's random field theory. The methods described in this paper are illustrated with the hippocampus surface data set. Using this case study, we will determine if there is an effect of family income on the growth of hippocampus in children in detail. There are a total of 124 children and 82 of them have repeat magnetic resonance images (MRI) two years later.
△ Less
Submitted 13 March, 2022;
originally announced March 2022.
-
Persistent Homological State-Space Estimation of Functional Human Brain Networks at Rest
Authors:
Moo K. Chung,
Shih-Gu Huang,
Ian C. Carroll,
Vince D. Calhoun,
H. Hill Goldsmith
Abstract:
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…
▽ More
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. MATLAB code for the method is available at https://github.com/laplcebeltrami/PH-STAT.
△ Less
Submitted 16 April, 2024; v1 submitted 31 December, 2021;
originally announced January 2022.
-
Hodge-Laplacian of Brain Networks
Authors:
D. Vijay Anand,
Moo K. Chung
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 January, 2023; v1 submitted 15 October, 2021;
originally announced October 2021.
-
Spatial mapping of protein composition and tissue organization: a primer for multiplexed antibody-based imaging
Authors:
John W. Hickey,
Elizabeth K. Neumann,
Andrea J. Radtke,
Jeannie M. Camarillo,
Rebecca T. Beuschel,
Alexandre Albanese,
Elizabeth McDonough,
Julia Hatler,
Anne E. Wiblin,
Jeremy Fisher,
Josh Croteau,
Eliza C. Small,
Anup Sood,
Richard M. Caprioli,
R. Michael Angelo,
Garry P. Nolan,
Kwanghun Chung,
Stephen M. Hewitt,
Ronald N. Germain,
Jeffrey M. Spraggins,
Emma Lundberg,
Michael P. Snyder,
Neil L. Kelleher,
Sinem K. Saka
Abstract:
Tissues and organs are composed of distinct cell types that must operate in concert to perform physiological functions. Efforts to create high-dimensional biomarker catalogs of these cells are largely based on transcriptomic single-cell approaches that lack the spatial context required to understand critical cellular communication and correlated structural organization. To probe in situ biology wi…
▽ More
Tissues and organs are composed of distinct cell types that must operate in concert to perform physiological functions. Efforts to create high-dimensional biomarker catalogs of these cells are largely based on transcriptomic single-cell approaches that lack the spatial context required to understand critical cellular communication and correlated structural organization. To probe in situ biology with sufficient coverage depth, several multiplexed protein imaging methods have recently been developed. Though these antibody-based technologies differ in strategy and mode of immunolabeling and detection tags, they commonly utilize antibodies directed against protein biomarkers to provide detailed spatial and functional maps of complex tissues. As these promising antibody-based multiplexing approaches become more widely adopted, new frameworks and considerations are critical for training future users, generating molecular tools, validating antibody panels, and harmonizing datasets. In this perspective, we provide essential resources and key considerations for obtaining robust and reproducible multiplexed antibody-based imaging data compiling specialized knowledge from domain experts and technology developers.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
Lattice Paths for Persistent Diagrams
Authors:
Moo K. Chung,
Hernando Ombao
Abstract:
Persistent homology has undergone significant development in recent years. However, one outstanding challenge is to build a coherent statistical inference procedure on persistent diagrams. In this paper, we first present a new lattice path representation for persistent diagrams. We then develop a new exact statistical inference procedure for lattice paths via combinatorial enumerations. The lattic…
▽ More
Persistent homology has undergone significant development in recent years. However, one outstanding challenge is to build a coherent statistical inference procedure on persistent diagrams. In this paper, we first present a new lattice path representation for persistent diagrams. We then develop a new exact statistical inference procedure for lattice paths via combinatorial enumerations. The lattice path method is applied to the topological characterization of the protein structures of the COVID-19 virus. We demonstrate that there are topological changes during the conformational change of spike proteins.
△ Less
Submitted 30 July, 2021; v1 submitted 1 May, 2021;
originally announced May 2021.
-
Graph Theory in Brain Networks
Authors:
Moo K. Chung
Abstract:
Recent developments in graph theoretic analysis of complex networks have led to deeper understanding of brain networks. Many complex networks show similar macroscopic behaviors despite differences in the microscopic details. Probably two most often observed characteristics of complex networks are scale-free and small-world properties. In this paper, we will explore whether brain networks follow sc…
▽ More
Recent developments in graph theoretic analysis of complex networks have led to deeper understanding of brain networks. Many complex networks show similar macroscopic behaviors despite differences in the microscopic details. Probably two most often observed characteristics of complex networks are scale-free and small-world properties. In this paper, we will explore whether brain networks follow scale-free and small-worldness among other graph theory properties.
△ Less
Submitted 9 March, 2021;
originally announced March 2021.
-
Reviews: Topological Distances and Losses for Brain Networks
Authors:
Moo K. Chung,
Alexander Smith,
Gary Shiu
Abstract:
Almost all statistical and machine learning methods in analyzing brain networks rely on distances and loss functions, which are mostly Euclidean or matrix norms. The Euclidean or matrix distances may fail to capture underlying subtle topological differences in brain networks. Further, Euclidean distances are sensitive to outliers. A few extreme edge weights may severely affect the distance. Thus i…
▽ More
Almost all statistical and machine learning methods in analyzing brain networks rely on distances and loss functions, which are mostly Euclidean or matrix norms. The Euclidean or matrix distances may fail to capture underlying subtle topological differences in brain networks. Further, Euclidean distances are sensitive to outliers. A few extreme edge weights may severely affect the distance. Thus it is necessary to use distances and loss functions that recognize topology of data. In this review paper, we survey various topological distance and loss functions from topological data analysis (TDA) and persistent homology that can be used in brain network analysis more effectively. Although there are many recent brain imaging studies that are based on TDA methods, possibly due to the lack of method awareness, TDA has not taken as the mainstream tool in brain imaging field yet. The main purpose of this paper is provide the relevant technical survey of these powerful tools that are immediately applicable to brain network data.
△ Less
Submitted 17 February, 2021;
originally announced February 2021.
-
Topological Learning for Brain Networks
Authors:
Tananun Songdechakraiwut,
Moo K. Chung
Abstract:
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 exte…
▽ More
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.
△ Less
Submitted 26 January, 2023; v1 submitted 25 November, 2020;
originally announced December 2020.
-
Active Learning Pipeline for Brain Mapping in a High Performance Computing Environment
Authors:
Adam Michaleas,
Lars A. Gjesteby,
Michael Snyder,
David Chavez,
Meagan Ash,
Matthew A. Melton,
Damon G. Lamb,
Sara N. Burke,
Kevin J. Otto,
Lee Kamentsky,
Webster Guan,
Kwanghun Chung,
Laura J. Brattain
Abstract:
This paper describes a scalable active learning pipeline prototype for large-scale brain mapping that leverages high performance computing power. It enables high-throughput evaluation of algorithm results, which, after human review, are used for iterative machine learning model training. Image processing and machine learning are performed in a batch layer. Benchmark testing of image processing usi…
▽ More
This paper describes a scalable active learning pipeline prototype for large-scale brain mapping that leverages high performance computing power. It enables high-throughput evaluation of algorithm results, which, after human review, are used for iterative machine learning model training. Image processing and machine learning are performed in a batch layer. Benchmark testing of image processing using pMATLAB shows that a 100$\times$ increase in throughput (10,000%) can be achieved while total processing time only increases by 9% on Xeon-G6 CPUs and by 22% on Xeon-E5 CPUs, indicating robust scalability. The images and algorithm results are provided through a serving layer to a browser-based user interface for interactive review. This pipeline has the potential to greatly reduce the manual annotation burden and improve the overall performance of machine learning-based brain mapping.
△ Less
Submitted 25 June, 2020;
originally announced June 2020.
-
Statistical Analysis of Dynamic Functional Brain Networks in Twins
Authors:
Moo K. Chung,
Shih-Gu Huang,
Tananun Songdechakraiwut,
Ian C. Carroll,
H. Hill Goldsmith
Abstract:
Recent studies have shown that functional brain brainwork is dynamic even during rest. A common approach to modeling the brain network in whole brain resting-state fMRI is to compute the correlation between anatomical regions via sliding windows. However, the direct use of the sample correlation matrices is not reliable due to the image acquisition, processing noises and the use of discrete window…
▽ More
Recent studies have shown that functional brain brainwork is dynamic even during rest. A common approach to modeling the brain network in whole brain resting-state fMRI is to compute the correlation between anatomical regions via sliding windows. However, the direct use of the sample correlation matrices is not reliable due to the image acquisition, processing noises and the use of discrete windows that often introduce spurious high-frequency fluctuations and the zig-zag pattern in the estimated time-varying correlation measures. To address the problem and obtain more robust correlation estimates, we propose the heat kernel based dynamic correlations. We demonstrate that the proposed heat kernel method can smooth out the unwanted high-frequency fluctuations in correlation estimations and achieve higher accuracy in identifying dynamically changing distinct states. The method is further used in determining if such dynamic state change is genetically heritable using a large-scale twin study. Various methodological challenges for analyzing paired twin dynamic networks are addressed.
△ Less
Submitted 11 October, 2020; v1 submitted 6 November, 2019;
originally announced November 2019.
-
Statistical Model for Dynamically-Changing Correlation Matrices with Application to Brain Connectivity
Authors:
Shih-Gu Huang,
S. Balqis Samdin,
Chee-Ming Ting,
Hernando Ombao,
Moo K. Chung
Abstract:
Background: Recent studies have indicated that functional connectivity is dynamic even during rest. A common approach to modeling the dynamic functional connectivity in whole-brain resting-state fMRI is to compute the correlation between anatomical regions via sliding time windows. However, the direct use of the sample correlation matrices is not reliable due to the image acquisition and processin…
▽ More
Background: Recent studies have indicated that functional connectivity is dynamic even during rest. A common approach to modeling the dynamic functional connectivity in whole-brain resting-state fMRI is to compute the correlation between anatomical regions via sliding time windows. However, the direct use of the sample correlation matrices is not reliable due to the image acquisition and processing noises in resting-sate fMRI.
New method: To overcome these limitations, we propose a new statistical model that smooths out the noise by exploiting the geometric structure of correlation matrices. The dynamic correlation matrix is modeled as a linear combination of symmetric positive-definite matrices combined with cosine series representation. The resulting smoothed dynamic correlation matrices are clustered into disjoint brain connectivity states using the k-means clustering algorithm.
Results: The proposed model preserves the geometric structure of underlying physiological dynamic correlation, eliminates unwanted noise in connectivity and obtains more accurate state spaces. The difference in the estimated dynamic connectivity states between males and females is identified.
Comparison with existing methods: We demonstrate that the proposed statistical model has less rapid state changes caused by noise and improves the accuracy in identifying and discriminating different states.
Conclusions: We propose a new regression model on dynamically changing correlation matrices that provides better performance over existing windowed correlation and is more reliable for the modeling of dynamic connectivity.
△ Less
Submitted 3 November, 2019; v1 submitted 25 December, 2018;
originally announced December 2018.
-
Harmonic holes as the submodules of brain network and network dissimilarity
Authors:
Hyekyoung Lee,
Moo K. Chung,
Hongyoon Choi,
Hyejin Kang,
Seunggyun Ha,
Yu Kyeong Kim,
Dong Soo Lee
Abstract:
Persistent homology has been applied to brain network analysis for finding the shape of brain networks across multiple thresholds. In the persistent homology, the shape of networks is often quantified by the sequence of $k$-dimensional holes and Betti numbers.The Betti numbers are more widely used than holes themselves in topological brain network analysis. However, the holes show the local connec…
▽ More
Persistent homology has been applied to brain network analysis for finding the shape of brain networks across multiple thresholds. In the persistent homology, the shape of networks is often quantified by the sequence of $k$-dimensional holes and Betti numbers.The Betti numbers are more widely used than holes themselves in topological brain network analysis. However, the holes show the local connectivity of networks, and they can be very informative features in analysis. In this study, we propose a new method of measuring network differences based on the dissimilarity measure of harmonic holes (HHs). The HHs, which represent the substructure of brain networks, are extracted by the Hodge Laplacian of brain networks. We also find the most contributed HHs to the network difference based on the HH dissimilarity. We applied our proposed method to clustering the networks of 4 groups, normal control (NC), stable and progressive mild cognitive impairment (sMCI and pMCI), and Alzheimer's disease (AD). The results showed that the clustering performance of the proposed method was better than that of network distances based on only the global change of topology.
△ Less
Submitted 11 November, 2018;
originally announced November 2018.
-
Topological Brain Network Distances
Authors:
Moo K. Chung,
Hyekyoung Lee,
Andrey Gritsenko,
Alex DiChristofano,
Dustin Pluta,
Hernando Ombao,
Victor Solo
Abstract:
Existing brain network distances are often based on matrix norms. The element-wise differences in the existing matrix norms may fail to capture underlying topological differences. Further, matrix norms are sensitive to outliers. A major disadvantage to element-wise distance calculations is that it could be severely affected even by a small number of extreme edge weights. Thus it is necessary to de…
▽ More
Existing brain network distances are often based on matrix norms. The element-wise differences in the existing matrix norms may fail to capture underlying topological differences. Further, matrix norms are sensitive to outliers. A major disadvantage to element-wise distance calculations is that it could be severely affected even by a small number of extreme edge weights. Thus it is necessary to develop network distances that recognize topology. In this paper, we provide a survey of bottleneck, Gromov-Hausdorff (GH) and Kolmogorov-Smirnov (KS) distances that are adapted for brain networks, and compare them against matrix-norm based network distances. Bottleneck and GH-distances are often used in persistent homology. However, they were rarely utilized to measure similarity between brain networks. KS-distance is recently introduced to measure the similarity between networks across different filtration values. The performance analysis was conducted using the random network simulations with the ground truths. Using a twin imaging study, which provides biological ground truth, we demonstrate that the KS distance has the ability to determine heritability.
△ Less
Submitted 9 September, 2018;
originally announced September 2018.
-
Exact Combinatorial Inference for Brain Images
Authors:
Moo K. Chung,
Zhan Luo,
Alex D. Leow,
Andrew L. Alexander,
Richard J. Davidson,
H. Hill Goldsmith
Abstract:
The permutation test is known as the exact test procedure in statistics. However, often it is not exact in practice and only an approximate method since only a small fraction of every possible permutation is generated. Even for a small sample size, it often requires to generate tens of thousands permutations, which can be a serious computational bottleneck. In this paper, we propose a novel combin…
▽ More
The permutation test is known as the exact test procedure in statistics. However, often it is not exact in practice and only an approximate method since only a small fraction of every possible permutation is generated. Even for a small sample size, it often requires to generate tens of thousands permutations, which can be a serious computational bottleneck. In this paper, we propose a novel combinatorial inference procedure that enumerates all possible permutations combinatorially without any resampling. The proposed method is validated against the standard permutation test in simulation studies with the ground truth. The method is further applied in twin DTI study in determining the genetic contribution of the minimum spanning tree of the structural brain connectivity.
△ Less
Submitted 8 July, 2018;
originally announced July 2018.
-
Statistical Challenges of Big Brain Network Data
Authors:
Moo K. Chung
Abstract:
We explore the main characteristics of big brain network data that offer unique statistical challenges. The brain networks are biologically expected to be both sparse and hierarchical. Such unique characterizations put specific topological constraints onto statistical approaches and models we can use effectively. We explore the limitations of the current models used in the field and offer alternat…
▽ More
We explore the main characteristics of big brain network data that offer unique statistical challenges. The brain networks are biologically expected to be both sparse and hierarchical. Such unique characterizations put specific topological constraints onto statistical approaches and models we can use effectively. We explore the limitations of the current models used in the field and offer alternative approaches and explain new challenges.
△ Less
Submitted 23 December, 2017; v1 submitted 27 September, 2017;
originally announced September 2017.
-
Topological Distances between Networks and Its Application to Brain Imaging
Authors:
Hyekyoung Lee,
Zhiwei Ma,
Yuan Wang,
Moo K. Chung
Abstract:
This paper surveys various distance measures for networks and graphs that were introduced in persistent homology. The scope of the paper is limited to network distances that were actually used in brain networks but the methods can be easily adapted to any weighted graph in other fields. The network version of Gromov-Hausdorff, bottleneck, kernel distances are introduced. We also introduce a recent…
▽ More
This paper surveys various distance measures for networks and graphs that were introduced in persistent homology. The scope of the paper is limited to network distances that were actually used in brain networks but the methods can be easily adapted to any weighted graph in other fields. The network version of Gromov-Hausdorff, bottleneck, kernel distances are introduced. We also introduce a recently developed KS-test like distance based on monotonic topology features such as the zeroth Betti number. Numerous toy examples and the result of applying many different distances to the brain networks of different clinical status and populations are given.
△ Less
Submitted 11 July, 2017; v1 submitted 16 January, 2017;
originally announced January 2017.
-
The Importance of Being Negative: A serious treatment of non-trivial edges in brain functional connectome
Authors:
Liang Zhan,
Lisanne M. Jenkins,
Ouri E. Wolfson,
Johnson J. GadElkarim,
Kevin Nocito,
Paul M. Thompson,
Olusola A. Ajilore,
Moo K. Chung,
Alex D. Leow
Abstract:
Understanding the modularity of fMRI-derived brain networks or connectomes can inform the study of brain function organization. However, fMRI connectomes additionally involve negative edges, which are not rigorously accounted for by existing approaches to modularity that either ignores or arbitrarily weight these connections. Furthermore, most Q maximization-based modularity algorithms yield varia…
▽ More
Understanding the modularity of fMRI-derived brain networks or connectomes can inform the study of brain function organization. However, fMRI connectomes additionally involve negative edges, which are not rigorously accounted for by existing approaches to modularity that either ignores or arbitrarily weight these connections. Furthermore, most Q maximization-based modularity algorithms yield variable results with suboptimal reproducibility. Here we present an alternative, reproducible approach that exploits how frequent the BOLD-signal correlation between two nodes is negative. We validated this novel probability-based modularity approach on two independent publicly-available resting-state connectome dataset (the Human Connectome Project and the 1000 Functional Connectomes) and demonstrated that negative correlations alone are sufficient in understanding resting-state modularity. In fact, this approach a) permits a dual formulation, leading to equivalent solutions regardless of whether one considers positive or negative edges; b) is theoretically linked to the Ising model defined on the connectome, thus yielding modularity result that maximizes data likelihood. We additionally were able to detect sex differences in modularity that the most widely utilized methods did not. Results confirmed the superiority of our approach in that: a) correlations with the highest probability of being negative are consistently placed between modules, b) due to the equivalent dual forms, no arbitrary weighting factor is required to balance the influence between negative and positive correlations, unlike existing Q maximization-based modularity approaches. As datasets like HCP become widely available for analysis by the neuroscience community at large, appropriate computational tools to understand the neurobiological information of negative edges in fMRI connectomes are increasingly important.
△ Less
Submitted 5 June, 2017; v1 submitted 6 September, 2016;
originally announced September 2016.
-
Mapping Heritability of Large-Scale Brain Networks with a Billion Connections {\em via} Persistent Homology
Authors:
Moo K. Chung,
Victoria Vilalta-Gil,
Paul J. Rathouz,
Benjamin B. Lahey,
David H. Zald
Abstract:
In many human brain network studies, we do not have sufficient number (n) of images relative to the number (p) of voxels due to the prohibitively expensive cost of scanning enough subjects. Thus, brain network models usually suffer the small-n large-p problem. Such a problem is often remedied by sparse network models, which are usually solved numerically by optimizing L1-penalties. Unfortunately,…
▽ More
In many human brain network studies, we do not have sufficient number (n) of images relative to the number (p) of voxels due to the prohibitively expensive cost of scanning enough subjects. Thus, brain network models usually suffer the small-n large-p problem. Such a problem is often remedied by sparse network models, which are usually solved numerically by optimizing L1-penalties. Unfortunately, due to the computational bottleneck associated with optimizing L1-penalties, it is not practical to apply such methods to construct large-scale brain networks at the voxel-level. In this paper, we propose a new scalable sparse network model using cross-correlations that bypass the computational bottleneck. Our model can build sparse brain networks at the voxel level with p > 25000. Instead of using a single sparse parameter that may not be optimal in other studies and datasets, the computational speed gain enables us to analyze the collection of networks at every possible sparse parameter in a coherent mathematical framework via persistent homology. The method is subsequently applied in determining the extent of heritability on a functional brain network at the voxel-level for the first time using twin fMRI.
△ Less
Submitted 29 June, 2016; v1 submitted 15 September, 2015;
originally announced September 2015.
-
Integrated multimodal network approach to PET and MRI based on multidimensional persistent homology
Authors:
Hyekyoung Lee,
Hyejin Kang,
Moo K. Chung,
Seonhee Lim,
Bung-Nyun Kim,
Dong Soo Lee
Abstract:
Finding the underlying relationships among multiple imaging modalities in a coherent fashion is one of challenging problems in the multimodal analysis. In this study, we propose a novel multimodal network approach based on multidi- mensional persistent homology. In this extension of the previous threshold-free method of persistent homology, we visualize and discriminate the topological change of i…
▽ More
Finding the underlying relationships among multiple imaging modalities in a coherent fashion is one of challenging problems in the multimodal analysis. In this study, we propose a novel multimodal network approach based on multidi- mensional persistent homology. In this extension of the previous threshold-free method of persistent homology, we visualize and discriminate the topological change of integrated brain networks by varying not only threshold but also mixing ratios between two different imaging modalities. Moreover, we also pro- pose an integration method for multimodal networks, called one-dimensional projection, with a specific mixing ratio between modalities. We applied the proposed methods to PET and MRI data from 21 autism spectrum disorder (ASD) children and 10 pediatric control subjects. From the results, we found that the brain networks of ASD children and controls differ significantly, with ASD showing asymmetrical changes of connected structures between PET and MRI. The integrated MRI and PET networks showed that ASD children had weaker connections than controls within the visual cortex, between dorsal and ventral parts of the temporal pole, between frontal and parietal regions, and between the left perisylvian and other brain regions. These results provide a multidimensional homological understanding of disease-related PET and MRI networks that discloses the network association with ASD.
△ Less
Submitted 24 February, 2016; v1 submitted 16 October, 2014;
originally announced October 2014.
-
Generalized epidemic process on modular networks
Authors:
Kihong Chung,
Yongjoo Baek,
Daniel Kim,
Meesoon Ha,
Hawoong Jeong
Abstract:
Social reinforcement and modular structure are two salient features observed in the spreading of behavior through social contacts. In order to investigate the interplay between these two features, we study the generalized epidemic process on modular networks with equal-sized finite communities and adjustable modularity. Using the analytical approach originally applied to clique-based random networ…
▽ More
Social reinforcement and modular structure are two salient features observed in the spreading of behavior through social contacts. In order to investigate the interplay between these two features, we study the generalized epidemic process on modular networks with equal-sized finite communities and adjustable modularity. Using the analytical approach originally applied to clique-based random networks, we show that the system exhibits a bond-percolation type continuous phase transition for weak social reinforcement, whereas a discontinuous phase transition occurs for sufficiently strong social reinforcement. Our findings are numerically verified using the finite-size scaling analysis and the crossings of the bimodality coefficient.
△ Less
Submitted 22 May, 2014; v1 submitted 29 November, 2013;
originally announced December 2013.
-
The Open Connectome Project Data Cluster: Scalable Analysis and Vision for High-Throughput Neuroscience
Authors:
Randal Burns,
William Gray Roncal,
Dean Kleissas,
Kunal Lillaney,
Priya Manavalan,
Eric Perlman,
Daniel R. Berger,
Davi D. Bock,
Kwanghun Chung,
Logan Grosenick,
Narayanan Kasthuri,
Nicholas C. Weiler,
Karl Deisseroth,
Michael Kazhdan,
Jeff Lichtman,
R. Clay Reid,
Stephen J. Smith,
Alexander S. Szalay,
Joshua T. Vogelstein,
R. Jacob Vogelstein
Abstract:
We describe a scalable database cluster for the spatial analysis and annotation of high-throughput brain imaging data, initially for 3-d electron microscopy image stacks, but for time-series and multi-channel data as well. The system was designed primarily for workloads that build connectomes---neural connectivity maps of the brain---using the parallel execution of computer vision algorithms on hi…
▽ More
We describe a scalable database cluster for the spatial analysis and annotation of high-throughput brain imaging data, initially for 3-d electron microscopy image stacks, but for time-series and multi-channel data as well. The system was designed primarily for workloads that build connectomes---neural connectivity maps of the brain---using the parallel execution of computer vision algorithms on high-performance compute clusters. These services and open-science data sets are publicly available at http://openconnecto.me.
The system design inherits much from NoSQL scale-out and data-intensive computing architectures. We distribute data to cluster nodes by partitioning a spatial index. We direct I/O to different systems---reads to parallel disk arrays and writes to solid-state storage---to avoid I/O interference and maximize throughput. All programming interfaces are RESTful Web services, which are simple and stateless, improving scalability and usability. We include a performance evaluation of the production system, highlighting the effectiveness of spatial data organization.
△ Less
Submitted 18 June, 2013; v1 submitted 14 June, 2013;
originally announced June 2013.