The importance of the whole: Topological data analysis for the network neuroscientist

Juvenile Myoclonic Epilepsy Imaging Endophenotypes and Relationship With Cognition and Resting-State EEG

Structural neuroimaging studies of patients with Juvenile Myoclonic Epilepsy (JME) typically present two findings: 1-volume reduction of subcortical gray matter structures, and 2-abnormalities of cortical thickness. The general trend has been to observe increased cortical thickness primarily in medial frontal regions, but heterogeneity across studies is common, including reports of decreased cortical thickness. These differences have not been explained. The cohort of patients investigated here originates from the Juvenile Myoclonic Epilepsy Connectome Project, which included comprehensive neuropsychological testing, 3 T MRI, and high-density 256-channel EEG. 64 JME patients aged 12-25 and 41 age and sex-matched healthy controls were included. Data-driven approaches were used to compare cortical thickness and subcortical volumes between the JME and control participants. After differences were identified, supervised machine learning was used to confirm their classification power. K-means clustering was used to generate imaging endophenotypes, which were then correlated with cognition, EEG frequency band lagged coherence from resting state high-density EEG, and white and grey matter based spatial statistics from diffusion imaging. The volumes of subcortical gray matter structures, particularly the thalamus and the motor-associated thalamic nuclei (ventral anterior), were found to be smaller in JME. In addition, the right hemisphere (primarily) sulcal pre-motor cortex was abnormally thicker in an age-dependent manner in JME with an asymmetry in the pre-motor cortical findings. These results suggested that for some patients JME may be an asymmetric disease, at least at the cortical level. Cluster analysis revealed three discrete imaging endophenotypes (left, right, symmetric). Clinically, the groups were not substantially different except for cognition, where left hemisphere disease was linked with a lower performance on a general cognitive factor ("g"). HD-EEG demonstrated statistically significant differences between imaging endophenotypes. Tract-based spatial statistics showed significant changes between endophenotypes as well. The left dominant disease group exhibited diffuse white matter changes. JME patients present with heterogeneity in underlying imaging endophenotypes that are defined by the presence and laterality of asymmetric abnormality at the level of the pre-motor sulcal cortex; these endophenotypes are linked to orderly relationships with cognition, EEG, and white matter pathology. The relationship of JME's adolescent onset, age-dependent cortical thickness loss, and seizure upon awakening all suggest that synaptic pruning may be a key element in the pathogenesis of JME. Individualized treatment approaches for neuromodulation are needed to target the most relevant cortical and subcortical structures as well as develop disease-modifying and neuroprotective strategies.

Heterogeneous and higher-order cortical connectivity undergirds efficient, robust, and reliable neural codes

We hypothesized that the heterogeneous architecture of biological neural networks provides a substrate to regulate the well-known tradeoff between robustness and efficiency, thereby allowing different subpopulations of the same network to optimize for different objectives. To distinguish between subpopulations, we developed a metric based on the mathematical theory of simplicial complexes that captures the complexity of their connectivity by contrasting its higher-order structure to a random control and confirmed its relevance in several openly available connectomes. Using a biologically detailed cortical model and an electron microscopic dataset, we showed that subpopulations with low simplicial complexity exhibit efficient activity. Conversely, subpopulations of high simplicial complexity play a supporting role in boosting the reliability of the network as a whole, softening the robustness-efficiency tradeoff. Crucially, we found that both types of subpopulations can and do coexist within a single connectome in biological neural networks, due to the heterogeneity of their connectivity.

Higher-order connectomics of human brain function reveals local topological signatures of task decoding, individual identification, and behavior

Traditional models of human brain activity often represent it as a network of pairwise interactions between brain regions. Going beyond this limitation, recent approaches have been proposed to infer higher-order interactions from temporal brain signals involving three or more regions. However, to this day it remains unclear whether methods based on inferred higher-order interactions outperform traditional pairwise ones for the analysis of fMRI data. To address this question, we conducted a comprehensive analysis using fMRI time series of 100 unrelated subjects from the Human Connectome Project. We show that higher-order approaches greatly enhance our ability to decode dynamically between various tasks, to improve the individual identification of unimodal and transmodal functional subsystems, and to strengthen significantly the associations between brain activity and behavior. Overall, our approach sheds new light on the higher-order organization of fMRI time series, improving the characterization of dynamic group dependencies in rest and tasks, and revealing a vast space of unexplored structures within human functional brain data, which may remain hidden when using traditional pairwise approaches.

Specific inhibition and disinhibition in the higher-order structure of a cortical connectome

Neurons are thought to act as parts of assemblies with strong internal excitatory connectivity. Conversely, inhibition is often reduced to blanket inhibition with no targeting specificity. We analyzed the structure of excitation and inhibition in the MICrONS $mm^{3}$ dataset, an electron microscopic reconstruction of a piece of cortical tissue. We found that excitation was structured around a feed-forward flow in large non-random neuron motifs with a structure of information flow from a small number of sources to a larger number of potential targets. Inhibitory neurons connected with neurons in specific sequential positions of these motifs, implementing targeted and symmetrical competition between them. None of these trends are detectable in only pairwise connectivity, demonstrating that inhibition is structured by these large motifs. While descriptions of inhibition in cortical circuits range from non-specific blanket-inhibition to targeted, our results describe a form of targeting specificity existing in the higher-order structure of the connectome. These findings have important implications for the role of inhibition in learning and synaptic plasticity.

Kernel machine tests of association using extrinsic and intrinsic cluster evaluation metrics

Modeling the network topology of the human brain within the mesoscale has become an increasing focus within the neuroscientific community due to its variation across diverse cognitive processes, in the presence of neuropsychiatric disease or injury, and over the lifespan. Much research has been done on the creation of algorithms to detect these mesoscopic structures, called communities or modules, but less has been done to conduct inference on these structures. The literature on analysis of these community detection algorithms has focused on comparing them within the same subject. These approaches, however, either do not accomodate a more general association between community structure and an outcome or cannot accommodate additional covariates that may confound the association of interest. We propose a semiparametric kernel machine regression model for either a continuous or binary outcome, where covariate effects are modeled parametrically and brain connectivity measures are measured nonparametrically. By incorporating notions of similarity between network community structures into a kernel distance function, the high-dimensional feature space of brain networks, defined on input pairs, can be generalized to non-linear spaces, allowing for a wider class of distance-based algorithms. We evaluate our proposed methodology on both simulated and real datasets.

Functional loops: Monitoring functional organization of deep neural networks using algebraic topology

Various topological methods have emerged in recent years to investigate the inner workings of deep neural networks (DNNs) based on the structural and weight information. However, their effectiveness is restricted due to the stratified structure and volatile weight information. In this study, we explore the relationship between functional organizations and network performance using algebraic topology. Our results indicate that functional loops reveal functional interaction patterns of multiple neurons in DNNs. We also propose functional persistence as a measure of functional complexity and develop an early stopping criterion that achieves competitive results without requiring a validation set.

Topology-based Clustering of Functional Brain Networks in an Alzheimer's Disease Cohort

Alzheimer's disease is a progressive neurodegenerative disease with many identifying biomarkers for diagnosis. However, whole-brain phenomena, particularly in functional MRI modalities, are not fully understood nor characterized. Here we employ the novel application of topological data analysis (TDA)-based methods of persistent homology to functional brain networks from ADNI-3 cohort to perform a subtyping experiment using unsupervised clustering techniques. We then investigate variations in QT-PAD challenge features across the identified clusters. Using a Wasserstein distance kernel with a variety of clustering algorithms, we found that the 0th-homology Wasserstein distance kernel and spectral clustering yielded clusters with significant differences in whole brain and medial temporal lobe (MTL) volume, thus demonstrating an intrinsic link between whole brain functional topology and brain morphometric structure. These findings demonstrate the importance of MTL in functional connectivity and the efficacy of using TDA-based machine learning methods in network neuroscience and neurodegenerative disease subtyping.

From pixels to connections: exploring in vitro neuron reconstruction software for network graph generation

Digital reconstruction has been instrumental in deciphering how in vitro neuron architecture shapes information flow. Emerging approaches reconstruct neural systems as networks with the aim of understanding their organization through graph theory. Computational tools dedicated to this objective build models of nodes and edges based on key cellular features such as somata, axons, and dendrites. Fully automatic implementations of these tools are readily available, but they may also be purpose-built from specialized algorithms in the form of multi-step pipelines. Here we review software tools informing the construction of network models, spanning from noise reduction and segmentation to full network reconstruction. The scope and core specifications of each tool are explicitly defined to assist bench scientists in selecting the most suitable option for their microscopy dataset. Existing tools provide a foundation for complete network reconstruction, however more progress is needed in establishing morphological bases for directed/weighted connectivity and in software validation.

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.

Topological structure of population activity in mouse visual cortex encodes densely sampled stimulus rotations

The primary visual cortex is one of the most well understood regions supporting the processing involved in sensory computation. Following the popularization of high-density neural recordings, it has been observed that the activity of large neural populations is often constrained to low dimensional manifolds. In this work, we quantify the structure of such neural manifolds in the visual cortex. We do this by analyzing publicly available two-photon optical recordings of mouse primary visual cortex in response to visual stimuli with a densely sampled rotation angle. Using a geodesic metric along with persistent homology, we discover that population activity in response to such stimuli generates a circular manifold, encoding the angle of rotation. Furthermore, we observe that this circular manifold is expressed differently in subpopulations of neurons with differing orientation and direction selectivity. Finally, we discuss some of the obstacles to reliably retrieving the truthful topology generated by a neural population.

Topological and geometric analysis of cell states in single-cell transcriptomic data

Single-cell RNA sequencing (scRNA-seq) enables dissecting cellular heterogeneity in tissues, resulting in numerous biological discoveries. Various computational methods have been devised to delineate cell types by clustering scRNA-seq data, where clusters are often annotated using prior knowledge of marker genes. In addition to identifying pure cell types, several methods have been developed to identify cells undergoing state transitions, which often rely on prior clustering results. The present computational approaches predominantly investigate the local and first-order structures of scRNA-seq data using graph representations, while scRNA-seq data frequently display complex high-dimensional structures. Here, we introduce scGeom, a tool that exploits the multiscale and multidimensional structures in scRNA-seq data by analyzing the geometry and topology through curvature and persistent homology of both cell and gene networks. We demonstrate the utility of these structural features to reflect biological properties and functions in several applications, where we show that curvatures and topological signatures of cell and gene networks can help indicate transition cells and the differentiation potential of cells. We also illustrate that structural characteristics can improve the classification of cell types.

Psilocybin enhances insightfulness in meditation: a perspective on the global topology of brain imaging during meditation

In this study, for the first time, we explored a dataset of functional magnetic resonance images collected during focused attention and open monitoring meditation before and after a five-day psilocybin-assisted meditation retreat using a recently established approach, based on the Mapper algorithm from topological data analysis. After generating subject-specific maps for two groups (psilocybin vs. placebo, 18 subjects/group) of experienced meditators, organizational principles were uncovered using graph topological tools, including the optimal transport (OT) distance, a geometrically rich measure of similarity between brain activity patterns. This revealed characteristics of the topology (i.e. shape) in space (i.e. abstract space of voxels) and time dimension of whole-brain activity patterns during different styles of meditation and psilocybin-induced alterations. Most interestingly, we found that (psilocybin-induced) positive derealization, which fosters insightfulness specifically when accompanied by enhanced open-monitoring meditation, was linked to the OT distance between open-monitoring and resting state. Our findings suggest that enhanced meta-awareness through meditation practice in experienced meditators combined with potential psilocybin-induced positive alterations in perception mediate insightfulness. Together, these findings provide a novel perspective on meditation and psychedelics that may reveal potential novel brain markers for positive synergistic effects between mindfulness practices and psilocybin.

Distinct serum GDNF coupling with brain structural and functional changes underlies cognitive status in Parkinson's disease

AIM

Aberrations in brain connections are implicated in the pathogenesis of Parkinson's disease (PD). We previously demonstrated that Glial cell-derived neurotrophic factor (GDNF) reduction is associated with cognition decline. Nonetheless, it is elusive if the pattern of brain topological connectivity differed across PD with divergent serum GDNF levels, and the accompanying profile of cognitive deficits has yet to be determined.

METHODS

We collected data on the participants' cognition, demographics, and serum GDNF levels. Participants underwent 3.0T magnetic resonance imaging, and we assessed the degree centrality, brain network topology, and cortical thickness of the healthy control (HC) (n = 25), PD-high-GDNF (n = 19), and PD-low-GDNF (n = 19) groups using graph-theoretic measures of resting-state functional MRI to reveal how much brain connectivity varies and its clinical correlates, as well as to determine factors predicting the cognitive status in PD.

RESULTS

The results show different network properties between groups. Degree centrality abnormalities were found in the right inferior frontal gyrus and right parietal lobe postcentral gyrus, linked with cognition scores. The two aberrant clusters serve as a potentially powerful signal for determining whether a patient has PD and the patient's cognition level after integrating with GDNF, duration, and dopamine dosage. Moreover, we found a significant positive relationship between the thickness of the left caudal middle frontal lobe and a plethora of cognitive domains. Further discriminant analysis revealed that the cortical thickness of this region could distinguish PD patients from healthy controls. The mental state evaluation will also be more precise when paired with GDNF and duration.

CONCLUSION

Our findings reveal that the topological features of brain networks and cortical thickness are altered in PD patients with cognitive deficits. The above change, accompanied by the serum GDNF, may have merit as a diagnosis marker for PD and, arguably, cognition status.

Topological features of spike trains in recurrent spiking neural networks that are trained to generate spatiotemporal patterns

In this study, we focus on training recurrent spiking neural networks to generate spatiotemporal patterns in the form of closed two-dimensional trajectories. Spike trains in the trained networks are examined in terms of their dissimilarity using the Victor-Purpura distance. We apply algebraic topology methods to the matrices obtained by rank-ordering the entries of the distance matrices, specifically calculating the persistence barcodes and Betti curves. By comparing the features of different types of output patterns, we uncover the complex relations between low-dimensional target signals and the underlying multidimensional spike trains.

Mild cognitive impairment prediction and cognitive score regression in the elderly using EEG topological data analysis and machine learning with awareness assessed in affective reminiscent paradigm

Introduction

The main objective of this study is to evaluate working memory and determine EEG biomarkers that can assist in the field of health neuroscience. Our ultimate goal is to utilize this approach to predict the early signs of mild cognitive impairment (MCI) in healthy elderly individuals, which could potentially lead to dementia. The advancements in health neuroscience research have revealed that affective reminiscence stimulation is an effective method for developing EEG-based neuro-biomarkers that can detect the signs of MCI.

Methods

We use topological data analysis (TDA) on multivariate EEG data to extract features that can be used for unsupervised clustering, subsequent machine learning-based classification, and cognitive score regression. We perform EEG experiments to evaluate conscious awareness in affective reminiscent photography settings.

Results

We use EEG and interior photography to distinguish between healthy cognitive aging and MCI. Our clustering UMAP and random forest application accurately predict MCI stage and MoCA scores.

Discussion

Our team has successfully implemented TDA feature extraction, MCI classification, and an initial regression of MoCA scores. However, our study has certain limitations due to a small sample size of only 23 participants and an unbalanced class distribution. To enhance the accuracy and validity of our results, future research should focus on expanding the sample size, ensuring gender balance, and extending the study to a cross-cultural context.

Topological Data Analysis Captures Task-Driven fMRI Profiles in Individual Participants: A Classification Pipeline Based on Persistence

BOLD-based fMRI is the most widely used method for studying brain function. The BOLD signal while valuable, is beset with unique vulnerabilities. The most notable of these is the modest signal to noise ratio, and the relatively low temporal and spatial resolution. However, the high dimensional complexity of the BOLD signal also presents unique opportunities for functional discovery. Topological Data Analyses (TDA), a branch of mathematics optimized to search for specific classes of structure within high dimensional data may provide particularly valuable applications. In this investigation, we acquired fMRI data in the anterior cingulate cortex (ACC) using a basic motor control paradigm. Then, for each participant and each of three task conditions, fMRI signals in the ACC were summarized using two methods: a) TDA based methods of persistent homology and persistence landscapes and b) non-TDA based methods using a standard vectorization scheme. Finally, using machine learning (with support vector classifiers), classification accuracy of TDA and non-TDA vectorized data was tested across participants. In each participant, TDA-based classification out-performed the non-TDA based counterpart, suggesting that our TDA analytic pipeline better characterized task- and condition-induced structure in fMRI data in the ACC. Our results emphasize the value of TDA in characterizing task- and condition-induced structure in regional fMRI signals. In addition to providing our analytical tools for other users to emulate, we also discuss the unique role that TDA-based methods can play in the study of individual differences in the structure of functional brain signals in the healthy and the clinical brain.

Efficient, continual, and generalized learning in the brain - neural mechanism of Mental Schema 2.0

There has been tremendous progress in artificial neural networks (ANNs) over the past decade; however, the gap between ANNs and the biological brain as a learning device remains large. With the goal of closing this gap, this paper reviews learning mechanisms in the brain by focusing on three important issues in ANN research: efficiency, continuity, and generalization. We first discuss the method by which the brain utilizes a variety of self-organizing mechanisms to maximize learning efficiency, with a focus on the role of spontaneous activity of the brain in shaping synaptic connections to facilitate spatiotemporal learning and numerical processing. Then, we examined the neuronal mechanisms that enable lifelong continual learning, with a focus on memory replay during sleep and its implementation in brain-inspired ANNs. Finally, we explored the method by which the brain generalizes learned knowledge in new situations, particularly from the mathematical generalization perspective of topology. Besides a systematic comparison in learning mechanisms between the brain and ANNs, we propose "Mental Schema 2.0," a new computational property underlying the brain's unique learning ability that can be implemented in ANNs.

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.

Topology recapitulates morphogenesis of neuronal dendrites

Branching allows neurons to make synaptic contacts with large numbers of other neurons, facilitating the high connectivity of nervous systems. Neuronal arbors have geometric properties such as branch lengths and diameters that are optimal in that they maximize signaling speeds while minimizing construction costs. In this work, we asked whether neuronal arbors have topological properties that may also optimize their growth or function. We discovered that for a wide range of invertebrate and vertebrate neurons the distributions of their subtree sizes follow power laws, implying that they are scale invariant. The power-law exponent distinguishes different neuronal cell types. Postsynaptic spines and branchlets perturb scale invariance. Through simulations, we show that the subtree-size distribution depends on the symmetry of the branching rules governing arbor growth and that optimal morphologies are scale invariant. Thus, the subtree-size distribution is a topological property that recapitulates the functional morphology of dendrites.

Topological data analysis of the synchronization of a network of Rössler chaotic electronic oscillators

Synchronization study allows a better understanding of the exchange of information among systems. In this work, we study experimental data recorded from a set of Rössler-like chaotic electronic oscillators arranged in a complex network, where the interactions between the oscillators are given in terms of a connectivity matrix, and their intensity is controlled by a global coupling parameter. We use the zero and one persistent homology groups to characterize the point clouds obtained from the signals recorded in pairs of oscillators. We show that the normalized persistent entropy (NPE) allows us to characterize the effective coupling between pairs of oscillators because it tends to increase with the coupling strength and to decrease with the distance between the oscillators. We also observed that pairs of oscillators that have similar degrees and are nearest neighbors tend to have higher NPE values than pairs with different degrees. However, large variability is found in the NPE values. Comparing the NPE behavior with that of the phase-locking value (PLV, commonly used to evaluate the synchronization of phase oscillators), we find that for large enough coupling, PLV only displays a monotonic increase, while NPE shows a richer behavior that captures variations in the behavior of the oscillators. This is due to the fact that PLV only captures coupling-induced phase changes, while NPE also captures amplitude changes. Moreover, when we consider the same network but with Kuramoto phase oscillators, we also find that NPE captures the transition to synchronization (as it increases with the coupling strength), and it also decreases with the distance between the oscillators. Therefore, we propose NPE as a data analysis technique to try to differentiate pairs of oscillators that have strong effective coupling because they are first or near neighbors, from those that have weaker coupling because they are distant neighbors.

Machine Learning Interpretability Methods to Characterize Brain Network Dynamics in Epilepsy

The rapid adoption of machine learning (ML) algorithms in a wide range of biomedical applications has highlighted issues of trust and the lack of understanding regarding the results generated by ML algorithms. Recent studies have focused on developing interpretable ML models and establish guidelines for transparency and ethical use, ensuring the responsible integration of machine learning in healthcare. In this study, we demonstrate the effectiveness of ML interpretability methods to provide important insights into the dynamics of brain network interactions in epilepsy, a serious neurological disorder affecting more than 60 million persons worldwide. Using high-resolution intracranial electroencephalogram (EEG) recordings from a cohort of 16 patients, we developed high accuracy ML models to categorize these brain activity recordings into either seizure or non-seizure classes followed by a more complex task of delineating the different stages of seizure progression to different parts of the brain as a multi-class classification task. We applied three distinct types of interpretability methods to the high-accuracy ML models to gain an understanding of the relative contributions of different categories of brain interaction patterns, including multi-focii interactions, which play an important role in distinguishing between different states of the brain. The results of this study demonstrate for the first time that post-hoc interpretability methods enable us to understand why ML algorithms generate a given set of results and how variations in value of input values affect the accuracy of the ML algorithms. In particular, we show in this study that interpretability methods can be used to identify brain regions and interaction patterns that have a significant impact on seizure events. The results of this study highlight the importance of the integrated implementation of ML algorithms together with interpretability methods in aberrant brain network studies and the wider domain of biomedical research.

Tree representations of brain structural connectivity via persistent homology

The brain structural connectome is generated by a collection of white matter fiber bundles constructed from diffusion weighted MRI (dMRI), acting as highways for neural activity. There has been abundant interest in studying how the structural connectome varies across individuals in relation to their traits, ranging from age and gender to neuropsychiatric outcomes. After applying tractography to dMRI to get white matter fiber bundles, a key question is how to represent the brain connectome to facilitate statistical analyses relating connectomes to traits. The current standard divides the brain into regions of interest (ROIs), and then relies on an adjacency matrix (AM) representation. Each cell in the AM is a measure of connectivity, e.g., number of fiber curves, between a pair of ROIs. Although the AM representation is intuitive, a disadvantage is the high-dimensionality due to the large number of cells in the matrix. This article proposes a simpler tree representation of the brain connectome, which is motivated by ideas in computational topology and takes topological and biological information on the cortical surface into consideration. We demonstrate that our tree representation preserves useful information and interpretability, while reducing dimensionality to improve statistical and computational efficiency. Applications to data from the Human Connectome Project (HCP) are considered and code is provided for reproducing our analyses.

A topological classifier to characterize brain states: When shape matters more than variance

Despite the remarkable accuracies attained by machine learning classifiers to separate complex datasets in a supervised fashion, most of their operation falls short to provide an informed intuition about the structure of data, and, what is more important, about the phenomena being characterized by the given datasets. By contrast, topological data analysis (TDA) is devoted to study the shape of data clouds by means of persistence descriptors and provides a quantitative characterization of specific topological features of the dataset under scrutiny. Here we introduce a novel TDA-based classifier that works on the principle of assessing quantifiable changes on topological metrics caused by the addition of new input to a subset of data. We used this classifier with a high-dimensional electro-encephalographic (EEG) dataset recorded from eleven participants during a previous decision-making experiment in which three motivational states were induced through a manipulation of social pressure. We calculated silhouettes from persistence diagrams associated with each motivated state with a ready-made band-pass filtered version of these signals, and classified unlabeled signals according to their impact on each reference silhouette. Our results show that in addition to providing accuracies within the range of those of a nearest neighbour classifier, the TDA classifier provides formal intuition of the structure of the dataset as well as an estimate of its intrinsic dimension. Towards this end, we incorporated variance-based dimensionality reduction methods to our dataset and found that in most cases the accuracy of our TDA classifier remains essentially invariant beyond a certain dimension.

Persistent homology analysis of type 2 diabetes genome-wide association studies in protein-protein interaction networks

Genome-wide association studies (GWAS) involving increasing sample sizes have identified hundreds of genetic variants associated with complex diseases, such as type 2 diabetes (T2D); however, it is unclear how GWAS hits form unique topological structures in protein-protein interaction (PPI) networks. Using persistent homology, this study explores the evolution and persistence of the topological features of T2D GWAS hits in the PPI network with increasing p-value thresholds. We define an n-dimensional persistent disease module as a higher-order generalization of the largest connected component (LCC). The 0-dimensional persistent T2D disease module is the LCC of the T2D GWAS hits, which is significantly detected in the PPI network (196 nodes and 235 edges, P< 0.05). In the 1-dimensional homology group analysis, all 18 1-dimensional holes (loops) of the T2D GWAS hits persist over all p-value thresholds. The 1-dimensional persistent T2D disease module comprising these 18 persistent 1-dimensional holes is significantly larger than that expected by chance (59 nodes and 83 edges, P< 0.001), indicating a significant topological structure in the PPI network. Our computational topology framework potentially possesses broad applicability to other complex phenotypes in identifying topological features that play an important role in disease pathobiology.

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.

Simple topological task-based functional connectivity features predict longitudinal behavioral change of fluid reasoning in the RANN cohort

Recent attention has been given to topological data analysis (TDA), and more specifically persistent homology (PH), to identify the underlying shape of brain network connectivity beyond simple edge pairings by computing connective components across different connectivity thresholds (see Sizemore et al., 2019). In the present study, we applied PH to task-based functional connectivity, computing 0-dimension Betti (B0) curves and calculating the area under these curves (AUC); AUC indicates how quickly a single connected component is formed across correlation filtration thresholds, with lower values interpreted as potentially analogous to lower whole-brain system segregation (e.g., Gracia-Tabuenca et al., 2020). One hundred sixty-three participants from the Reference Ability Neural Network (RANN) longitudinal lifespan cohort (age 20-80 years) were tested in-scanner at baseline and five-year follow-up on a battery of tests comprising four domains of cognition (i.e., Stern et al., 2014). We tested for 1.) age-related change in the AUC of the B0 curve over time, 2.) the predictive utility of AUC in accounting for longitudinal change in behavioral performance and 3.) compared system segregation to the PH approach. Results demonstrated longitudinal age-related decreases in AUC for Fluid Reasoning, with these decreases predicting longitudinal declines in cognition, even after controlling for demographic and brain integrity factors; moreover, change in AUC partially mediated the effect of age on change in cognitive performance. System segregation also significantly decreased with age in three of the four cognitive domains but did not predict change in cognition. These results argue for greater application of TDA to the study of aging.

Manifold learning for fMRI time-varying functional connectivity

Whole-brain functional connectivity (FC) measured with functional MRI (fMRI) evolves over time in meaningful ways at temporal scales going from years (e.g., development) to seconds [e.g., within-scan time-varying FC (tvFC)]. Yet, our ability to explore tvFC is severely constrained by its large dimensionality (several thousands). To overcome this difficulty, researchers often seek to generate low dimensional representations (e.g., 2D and 3D scatter plots) hoping those will retain important aspects of the data (e.g., relationships to behavior and disease progression). Limited prior empirical work suggests that manifold learning techniques (MLTs)-namely those seeking to infer a low dimensional non-linear surface (i.e., the manifold) where most of the data lies-are good candidates for accomplishing this task. Here we explore this possibility in detail. First, we discuss why one should expect tvFC data to lie on a low dimensional manifold. Second, we estimate what is the intrinsic dimension (ID; i.e., minimum number of latent dimensions) of tvFC data manifolds. Third, we describe the inner workings of three state-of-the-art MLTs: Laplacian Eigenmaps (LEs), T-distributed Stochastic Neighbor Embedding (T-SNE), and Uniform Manifold Approximation and Projection (UMAP). For each method, we empirically evaluate its ability to generate neuro-biologically meaningful representations of tvFC data, as well as their robustness against hyper-parameter selection. Our results show that tvFC data has an ID that ranges between 4 and 26, and that ID varies significantly between rest and task states. We also show how all three methods can effectively capture subject identity and task being performed: UMAP and T-SNE can capture these two levels of detail concurrently, but LE could only capture one at a time. We observed substantial variability in embedding quality across MLTs, and within-MLT as a function of hyper-parameter selection. To help alleviate this issue, we provide heuristics that can inform future studies. Finally, we also demonstrate the importance of feature normalization when combining data across subjects and the role that temporal autocorrelation plays in the application of MLTs to tvFC data. Overall, we conclude that while MLTs can be useful to generate summary views of labeled tvFC data, their application to unlabeled data such as resting-state remains challenging.

Mining Correlation between Fluid Intelligence and Whole-brain Large Scale Structural Connectivity

Exploring the neural basis of intelligence and the corresponding associations with brain network has been an active area of research in network neuroscience. Up to now, the majority of explorations mining human intelligence in brain connectomics leverages whole-brain functional connectivity patterns. In this study, structural connectivity patterns are instead used to explore relationships between brain connectivity and different behavioral/cognitive measures such as fluid intelligence. Specifically, we conduct a study using the 397 unrelated subjects from Human Connectome Project (Young Adults) dataset to estimate individual level structural connectivity matrices. We show that topological features, as quantified by our proposed measurements: Average Persistence (AP) and Persistent Entropy (PE), has statistically significant associations with different behavioral/cognitive measures. We also perform a parallel study using traditional graph-theoretical measures, provided by Brain Connectivity Toolbox, as benchmarks for our study. Our findings indicate that individual's structural connectivity indeed offers reliable predictive power of different behavioral/cognitive measures, including but not limited to fluid intelligence. Our results suggest that structural connectomes provide complementary insights (compared to using functional connectomes) in predicting human intelligence and warrants future studies on human intelligence and/or other behavioral/cognitive measures involving multi-modal approach.

Topological insights into the neural basis of flexible behavior

It is widely accepted that there is an inextricable link between neural computations, biological mechanisms, and behavior, but it is challenging to simultaneously relate all three. Here, we show that topological data analysis (TDA) provides an important bridge between these approaches to studying how brains mediate behavior. We demonstrate that cognitive processes change the topological description of the shared activity of populations of visual neurons. These topological changes constrain and distinguish between competing mechanistic models, are connected to subjects' performance on a visual change detection task, and, via a link with network control theory, reveal a tradeoff between improving sensitivity to subtle visual stimulus changes and increasing the chance that the subject will stray off task. These connections provide a blueprint for using TDA to uncover the biological and computational mechanisms by which cognition affects behavior in health and disease.

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.

Multivariate information theory uncovers synergistic subsystems of the human cerebral cortex

One of the most well-established tools for modeling the brain is the functional connectivity network, which is constructed from pairs of interacting brain regions. While powerful, the network model is limited by the restriction that only pairwise dependencies are considered and potentially higher-order structures are missed. Here, we explore how multivariate information theory reveals higher-order dependencies in the human brain. We begin with a mathematical analysis of the O-information, showing analytically and numerically how it is related to previously established information theoretic measures of complexity. We then apply the O-information to brain data, showing that synergistic subsystems are widespread in the human brain. Highly synergistic subsystems typically sit between canonical functional networks, and may serve an integrative role. We then use simulated annealing to find maximally synergistic subsystems, finding that such systems typically comprise ≈10 brain regions, recruited from multiple canonical brain systems. Though ubiquitous, highly synergistic subsystems are invisible when considering pairwise functional connectivity, suggesting that higher-order dependencies form a kind of shadow structure that has been unrecognized by established network-based analyses. We assert that higher-order interactions in the brain represent an under-explored space that, accessible with tools of multivariate information theory, may offer novel scientific insights.

Development of the Functional Connectome Topology in Adolescence: Evidence from Topological Data Analysis

Adolescence is a crucial developmental period in terms of behavior and mental health. Therefore, understanding how the brain develops during this stage is a fundamental challenge for neuroscience. Recent studies have modeled the brain as a network or connectome, mainly applying measures from graph theory, showing a change in its functional organization, such as an increase in its segregation and integration. Topological Data Analysis (TDA) complements such modeling by extracting high-dimensional features across the whole range of connectivity values instead of exploring a fixed set of connections. This study inquires into the developmental trajectories of such properties using a longitudinal sample of typically developing human participants (N = 98; 53/45 female/male; 6.7-18.1 years), applying TDA to their functional connectomes. In addition, we explore the effect of puberty on individual developmental trajectories. Results showed that the adolescent brain has a more distributed topology structure compared with random networks but is more densely connected at the local level. Furthermore, developmental effects showed nonlinear trajectories for the topology of the whole brain and fronto-parietal networks, with an inflection point and increasing trajectories after puberty onset. These results add to the insights into the development of the functional organization of the adolescent brain.

Manifold Learning for fMRI time-varying FC

Whole-brain functional connectivity ( FC ) measured with functional MRI (fMRI) evolve over time in meaningful ways at temporal scales going from years (e.g., development) to seconds (e.g., within-scan time-varying FC ( tvFC )). Yet, our ability to explore tvFC is severely constrained by its large dimensionality (several thousands). To overcome this difficulty, researchers seek to generate low dimensional representations (e.g., 2D and 3D scatter plots) expected to retain its most informative aspects (e.g., relationships to behavior, disease progression). Limited prior empirical work suggests that manifold learning techniques ( MLTs )-namely those seeking to infer a low dimensional non-linear surface (i.e., the manifold) where most of the data lies-are good candidates for accomplishing this task. Here we explore this possibility in detail. First, we discuss why one should expect tv FC data to lie on a low dimensional manifold. Second, we estimate what is the intrinsic dimension (i.e., minimum number of latent dimensions; ID ) of tvFC data manifolds. Third, we describe the inner workings of three state-of-the-art MLTs : Laplacian Eigenmaps ( LE ), T-distributed Stochastic Neighbor Embedding ( T-SNE ), and Uniform Manifold Approximation and Projection ( UMAP ). For each method, we empirically evaluate its ability to generate neuro-biologically meaningful representations of tvFC data, as well as their robustness against hyper-parameter selection. Our results show that tvFC data has an ID that ranges between 4 and 26, and that ID varies significantly between rest and task states. We also show how all three methods can effectively capture subject identity and task being performed: UMAP and T-SNE can capture these two levels of detail concurrently, but L E could only capture one at a time. We observed substantial variability in embedding quality across MLTs , and within- MLT as a function of hyper-parameter selection. To help alleviate this issue, we provide heuristics that can inform future studies. Finally, we also demonstrate the importance of feature normalization when combining data across subjects and the role that temporal autocorrelation plays in the application of MLTs to tvFC data. Overall, we conclude that while MLTs can be useful to generate summary views of labeled tvFC data, their application to unlabeled data such as resting-state remains challenging.

A Topological Loss Function for Deep-Learning Based Image Segmentation Using Persistent Homology

We introduce a method for training neural networks to perform image or volume segmentation in which prior knowledge about the topology of the segmented object can be explicitly provided and then incorporated into the training process. By using the differentiable properties of persistent homology, a concept used in topological data analysis, we can specify the desired topology of segmented objects in terms of their Betti numbers and then drive the proposed segmentations to contain the specified topological features. Importantly this process does not require any ground-truth labels, just prior knowledge of the topology of the structure being segmented. We demonstrate our approach in four experiments: one on MNIST image denoising and digit recognition, one on left ventricular myocardium segmentation from magnetic resonance imaging data from the UK Biobank, one on the ACDC public challenge dataset and one on placenta segmentation from 3-D ultrasound. We find that embedding explicit prior knowledge in neural network segmentation tasks is most beneficial when the segmentation task is especially challenging and that it can be used in either a semi-supervised or post-processing context to extract a useful training gradient from images without pixelwise labels.

Inference of hyperedges and overlapping communities in hypergraphs

Hypergraphs, encoding structured interactions among any number of system units, have recently proven a successful tool to describe many real-world biological and social networks. Here we propose a framework based on statistical inference to characterize the structural organization of hypergraphs. The method allows to infer missing hyperedges of any size in a principled way, and to jointly detect overlapping communities in presence of higher-order interactions. Furthermore, our model has an efficient numerical implementation, and it runs faster than dyadic algorithms on pairwise records projected from higher-order data. We apply our method to a variety of real-world systems, showing strong performance in hyperedge prediction tasks, detecting communities well aligned with the information carried by interactions, and robustness against addition of noisy hyperedges. Our approach illustrates the fundamental advantages of a hypergraph probabilistic model when modeling relational systems with higher-order interactions.

Geometry of spiking patterns in early visual cortex: a topological data analytic approach

In the brain, spiking patterns live in a high-dimensional space of neurons and time. Thus, determining the intrinsic structure of this space presents a theoretical and experimental challenge. To address this challenge, we introduce a new framework for applying topological data analysis (TDA) to spike train data and use it to determine the geometry of spiking patterns in the visual cortex. Key to our approach is a parametrized family of distances based on the timing of spikes that quantifies the dissimilarity between neuronal responses. We applied TDA to visually driven single-unit and multiple single-unit spiking activity in macaque V1 and V2. TDA across timescales reveals a common geometry for spiking patterns in V1 and V2 which, among simple models, is most similar to that of a low-dimensional space endowed with Euclidean or hyperbolic geometry with modest curvature. Remarkably, the inferred geometry depends on timescale and is clearest for the timescales that are important for encoding contrast, orientation and spatial correlations.

From calcium imaging to graph topology

Systems neuroscience is facing an ever-growing mountain of data. Recent advances in protein engineering and microscopy have together led to a paradigm shift in neuroscience; using fluorescence, we can now image the activity of every neuron through the whole brain of behaving animals. Even in larger organisms, the number of neurons that we can record simultaneously is increasing exponentially with time. This increase in the dimensionality of the data is being met with an explosion of computational and mathematical methods, each using disparate terminology, distinct approaches, and diverse mathematical concepts. Here we collect, organize, and explain multiple data analysis techniques that have been, or could be, applied to whole-brain imaging, using larval zebrafish as an example model. We begin with methods such as linear regression that are designed to detect relations between two variables. Next, we progress through network science and applied topological methods, which focus on the patterns of relations among many variables. Finally, we highlight the potential of generative models that could provide testable hypotheses on wiring rules and network progression through time, or disease progression. While we use examples of imaging from larval zebrafish, these approaches are suitable for any population-scale neural network modeling, and indeed, to applications beyond systems neuroscience. Computational approaches from network science and applied topology are not limited to larval zebrafish, or even to systems neuroscience, and we therefore conclude with a discussion of how such methods can be applied to diverse problems across the biological sciences.

Brain Complexities and Anesthesia: Their Meaning and Measurement

A complex system is often associated with emergence of new phenomena from the interactions between the system's components. General anesthesia reduces brain complexity and so inhibits the emergence of consciousness. An understanding of complexity is necessary for the interpretation of brain monitoring algorithms. Complexity indices capture the "difficulty" of understanding brain activity over time and/or space. Complexity-entropy plots reveal the types of complexity indices and their balance of randomness and structure. Lempel-Ziv complexity is a common index of temporal complexity for single-channel electroencephalogram containing both power spectral and nonlinear effects, revealed by phase-randomized surrogate data. Computing spatial complexities involves forming a connectivity matrix and calculating the complexity of connectivity patterns. Spatiotemporal complexity can be estimated in multiple ways including temporal or spatial concatenation, estimation of state switching, or integrated information. This article illustrates the concept and application of various complexities by providing working examples; a website with interactive demonstrations has also been created.

Topological Data Analysis Helps to Improve Accuracy of Deep Learning Models for Fake News Detection Trained on Very Small Training Sets

Topological data analysis has recently found applications in various areas of science, such as computer vision and understanding of protein folding. However, applications of topological data analysis to natural language processing remain under-researched. This study applies topological data analysis to a particular natural language processing task: fake news detection. We have found that deep learning models are more accurate in this task than topological data analysis. However, assembling a deep learning model with topological data analysis significantly improves the model’s accuracy if the available training set is very small.

Multiscale topology characterizes dynamic tumor vascular networks

Advances in imaging techniques enable high-resolution three-dimensional (3D) visualization of vascular networks over time and reveal abnormal structural features such as twists and loops, and their quantification is an active area of research. Here, we showcase how topological data analysis, the mathematical field that studies the "shape" of data, can characterize the geometric, spatial, and temporal organization of vascular networks. We propose two topological lenses to study vasculature, which capture inherent multiscale features and vessel connectivity, and surpass the single-scale analysis of existing methods. We analyze images collected using intravital and ultramicroscopy modalities and quantify spatiotemporal variation of twists, loops, and avascular regions (voids) in 3D vascular networks. This topological approach validates and quantifies known qualitative trends such as dynamic changes in tortuosity and loops in response to antibodies that modulate vessel sprouting; furthermore, it quantifies the effect of radiotherapy on vessel architecture.

Multiscale topology characterizes dynamic tumor vascular networks

Advances in imaging techniques enable high-resolution three-dimensional (3D) visualization of vascular networks over time and reveal abnormal structural features such as twists and loops, and their quantification is an active area of research. Here, we showcase how topological data analysis, the mathematical field that studies the "shape" of data, can characterize the geometric, spatial, and temporal organization of vascular networks. We propose two topological lenses to study vasculature, which capture inherent multiscale features and vessel connectivity, and surpass the single-scale analysis of existing methods. We analyze images collected using intravital and ultramicroscopy modalities and quantify spatiotemporal variation of twists, loops, and avascular regions (voids) in 3D vascular networks. This topological approach validates and quantifies known qualitative trends such as dynamic changes in tortuosity and loops in response to antibodies that modulate vessel sprouting; furthermore, it quantifies the effect of radiotherapy on vessel architecture.

Topological data analysis in biomedicine: A review

Significant technological advances made in recent years have shepherded a dramatic increase in utilization of digital technologies for biomedicine- everything from the widespread use of electronic health records to improved medical imaging capabilities and the rising ubiquity of genomic sequencing contribute to a "digitization" of biomedical research and clinical care. With this shift toward computerized tools comes a dramatic increase in the amount of available data, and current tools for data analysis capable of extracting meaningful knowledge from this wealth of information have yet to catch up. This article seeks to provide an overview of emerging mathematical methods with the potential to improve the abilities of clinicians and researchers to analyze biomedical data, but may be hindered from doing so by a lack of conceptual accessibility and awareness in the life sciences research community. In particular, we focus on topological data analysis (TDA), a set of methods grounded in the mathematical field of algebraic topology that seeks to describe and harness features related to the "shape" of data. We aim to make such techniques more approachable to non-mathematicians by providing a conceptual discussion of their theoretical foundations followed by a survey of their published applications to scientific research. Finally, we discuss the limitations of these methods and suggest potential avenues for future work integrating mathematical tools into clinical care and biomedical informatics.

Analysis and Prevention and Control System of Domino Accident Risk Data in Chemical Parks Based on Topological Neural Network

A topologically based neural network algorithm is used to conduct an in-depth study and analysis of domino accident risk data in chemical parks, and this is used to construct a prevention and control system applied to the safety prediction of chemical parks. Firstly, the operating model of the flue gas turbine is expanded and analyzed according to the basic theory of topology, and the object element model is constructed to determine the feature vector and potential risk level. Then, the idea of differential evolution is introduced into the topological neural network to solve the problem that the learning rate and weighting coefficients are difficult to determine, and then the complete DE-ENN algorithm is proposed and tested with the UCI standard data set to verify the effectiveness of the algorithm. Finally, the algorithm is applied to the potential risk identification of the smoke machine operation model, and the experimental results show that the method not only has a simple structure, short running time, and high prediction accuracy but also has excellent generalization ability. For the inherent risk and domino effect risk of chemical equipment in chemical fiber enterprises, the accident risk assessment method based on the protection layer analysis method is proposed; combined with the probability of domino accident and personnel vulnerability model based on the comprehensive analysis of the research results of the allowable risk standard, the allowable risk standard applicable to chemical fiber production enterprises in China is proposed. Given the potential accident risk characteristics of chemical fiber production enterprises, the calculation method of firefighting demand and firefighting capacity of firefighting system is given; the index system of firefighting system emergency response capacity assessment is constructed from three aspects of firefighting system integrity, reliability, and effectiveness, and the assessment model and grade classification standard of firefighting emergency response capacity of chemical fiber production enterprises are determined.

A hands-on tutorial on network and topological neuroscience

The brain is an extraordinarily complex system that facilitates the optimal integration of information from different regions to execute its functions. With the recent advances in technology, researchers can now collect enormous amounts of data from the brain using neuroimaging at different scales and from numerous modalities. With that comes the need for sophisticated tools for analysis. The field of network neuroscience has been trying to tackle these challenges, and graph theory has been one of its essential branches through the investigation of brain networks. Recently, topological data analysis has gained more attention as an alternative framework by providing a set of metrics that go beyond pairwise connections and offer improved robustness against noise. In this hands-on tutorial, our goal is to provide the computational tools to explore neuroimaging data using these frameworks and to facilitate their accessibility, data visualisation, and comprehension for newcomers to the field. We will start by giving a concise (and by no means complete) overview of the field to introduce the two frameworks and then explain how to compute both well-established and newer metrics on resting-state functional magnetic resonance imaging. We use an open-source language (Python) and provide an accompanying publicly available Jupyter Notebook that uses the 1000 Functional Connectomes Project dataset. Moreover, we would like to highlight one part of our notebook dedicated to the realistic visualisation of high order interactions in brain networks. This pipeline provides three-dimensional (3-D) plots of pairwise and higher-order interactions projected in a brain atlas, a new feature tailor-made for network neuroscience.

Brain-wide visual habituation networks in wild type and fmr1 zebrafish

Habituation is a form of learning during which animals stop responding to repetitive stimuli, and deficits in habituation are characteristic of several psychiatric disorders. Due to technical challenges, the brain-wide networks mediating habituation are poorly understood. Here we report brain-wide calcium imaging during larval zebrafish habituation to repeated visual looming stimuli. We show that different functional categories of loom-sensitive neurons are located in characteristic locations throughout the brain, and that both the functional properties of their networks and the resulting behavior can be modulated by stimulus saliency and timing. Using graph theory, we identify a visual circuit that habituates minimally, a moderately habituating midbrain population proposed to mediate the sensorimotor transformation, and downstream circuit elements responsible for higher order representations and the delivery of behavior. Zebrafish larvae carrying a mutation in the fmr1 gene have a systematic shift toward sustained premotor activity in this network, and show slower behavioral habituation.

Geodesic-based distance reveals nonlinear topological features in neural activity from mouse visual cortex

An increasingly popular approach to the analysis of neural data is to treat activity patterns as being constrained to and sampled from a manifold, which can be characterized by its topology. The persistent homology method identifies the type and number of holes in the manifold, thereby yielding functional information about the coding and dynamic properties of the underlying neural network. In this work, we give examples of highly nonlinear manifolds in which the persistent homology algorithm fails when it uses the Euclidean distance because it does not always yield a good approximation to the true distance distribution of a point cloud sampled from a manifold. To deal with this issue, we instead estimate the geodesic distance which is a better approximation of the true distance distribution and can therefore be used to successfully identify highly nonlinear features with persistent homology. To document the utility of the method, we utilize a toy model comprised of a circular manifold, built from orthogonal sinusoidal coordinate functions and show how the choice of metric determines the performance of the persistent homology algorithm. Furthermore, we explore the robustness of the method across different manifold properties, like the number of samples, curvature and amount of added noise. We point out strategies for interpreting its results as well as some possible pitfalls of its application. Subsequently, we apply this analysis to neural data coming from the Visual Coding-Neuropixels dataset recorded at the Allen Institute in mouse visual cortex in response to stimulation with drifting gratings. We find that different manifolds with a non-trivial topology can be seen across regions and stimulus properties. Finally, we interpret how these changes in manifold topology along with stimulus parameters and cortical region inform how the brain performs visual computation.

Edges in brain networks: Contributions to models of structure and function

Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.