Abstract:Volume electron microscopy together with computer-based image analysis are yielding neural circuit diagrams of ever larger regions of the brain. These datasets are usually represented in a cell-to-cell connectivity matrix and contain important information about prevalent circuit motifs allowing to directly test various theories on the computation in that brain structure. Of particular interest are the detection of cell assemblies and the quantification of feedback, which can profoundly change circuit properties. While the ordering of cells along the rows and columns doesn't change the connectivity, it can make special connectivity patterns recognizable. For example, ordering the cells along the flow of information, feedback and feedforward connections are segregated above and below the main matrix diagonal, respectively. Different algorithms are used to renumber matrices such as to minimize a given cost function, but either their performance becomes unsatisfying at a given size of the circuit or the CPU time needed to compute them scales in an unfavorable way with increasing number of neurons. Based on previous ideas, I describe an algorithm which is effective in matrix reordering with respect to both its performance as well as to its scaling in computing time. Rather than trying to reorder the matrix in discrete steps, the algorithm transiently relaxes the integer program by assigning a real-valued parameter to each cell describing its location on a continuous axis ('smooth-index') and finds the parameter set that minimizes the cost. I find that the smooth-index algorithm outperforms all algorithms I compared it to, including those based on topological sorting. Connectomic data provide researchers with neural circuit diagrams of ever larger regions of the brain. These datasets are usually represented in a cell-to-cell connectivity matrix and contain important information about prevalent circuit motifs. Such motifs, however, only become visible if the connectivity matrix is reordered appropriately. For example, ordering the cells along the flow of information, feedback and feedforward connections are segregated above and below the main matrix diagonal, respectively. While most previous approaches rely on topological sorting, my method relaxes the discrete vertex indices to real numbers ('smooth-index') along independent parameter axes and defines a differentiable cost function, thus, allowing gradient-based algorithms to find a minimum. The parameter set at this minimum is then re-discretized to reorder the connectivity matrix accordingly. I find my method to scale favorably with the circuit size and to outperform all algorithms I compared it to.

Fast algorithm to identify cluster synchrony through fibration symmetries in large information-processing networks

Quasifibrations of Graphs to Find Symmetries in Biological Networks

Fibration symmetries uncover the building blocks of biological networks

A functorial Dowker theorem and persistent homology of asymmetric networks

Fibration symmetries and cluster synchronization in the Caenorhabditis elegans connectome

Predicting synchronized gene coexpression patterns from fibration symmetries in gene regulatory networks in bacteria

General Optimization Framework for Accurate and Efficient Reconstruction of Symmetric Complex Networks from Dynamical Data.

Endowing networks with desired symmetries and modular behavior

Sifting out communities in large sparse networks

Connectivity Matrix Seriation via Relaxation

Efficient Topology Reconstruction Via Machine Learning Based Traffic Patterns Recognition in Optically Interconnected Computing System.

Symmetries and synchronization from whole-neural activity in C. elegans connectome: Integration of functional and structural networks

Non-Negative Matrix Factorizations for Multiplex Network Analysis

Clustering Large Attributed Information Networks: an Efficient Incremental Computing Approach.

Dynamics and bifurcations in genetic circuits with fibration symmetries

Circuits with broken fibration symmetries perform core logic computations in biological networks

Hypergraph reconstruction from dynamics

FALCON or how to compute measures time efficiently on dynamically evolving dense complex networks?

Identifying Representative Network Motifs for Inferring Higher-order Structure of Biological Networks

Symmetries and synchronization from whole-neural activity in {\it C. elegans} connectome: Integration of functional and structural networks