Organizations: National Key Laboratory for Multimedia Information Processing, School of Computer Science, Peking University, Beijing, China · School of Electronics Engineering and Computer Science, Peking University, Beijing, China · School of Computer Science, Beijing University of Posts and Telecommunications, Beijing, China · College of Engineering, Peking University, Beijing, China · Yuanpei College, Peking University, Beijing, China
Connectomes, graph-level maps of neurons and their synaptic connections, provide a structural basis for understanding how brain circuits support function and computation. However, mapping connectome structure to computation remains difficult because these graphs are high-dimensional, sparse, and sensitive to local structural variation. Existing approaches often depend on hand-crafted structural descriptors or task-specific predictors, which limits their ability to represent connectomes in a form that is both generative and functionally meaningful. We propose a conditional generative latent framework that encodes connectome graphs into a compact structural space while using available node-level conditions to guide reconstruction and generation. From this space, the model can reconstruct observed connectivity with a mean edge-reconstruction AUC up to 0.910 and generate new candidate connectomes, enabling a unified analysis of graph structure and computational behavior. Using connectome-derived graphs as recurrent computational substrates, we found that the learned latent space captures functional variation across reservoir-computing experiments, with cross-validated R2 values up to approximately 0.87. Interpretability analysis further revealed task-specific structural mechanisms: in our examples, memory performance is associated with reciprocal recurrent connectivity, whereas prediction and classification are more strongly associated with spectral properties of the recurrent network. These findings suggest an AI-for-science approach to linking neural connectivity to computation and provide a generative and interpretable basis for studying how distinct structural mechanisms shape computational capacity.
Figures & tables
Figure 1: Connectome-to-function latent space. Observed connectomes are encoded into a generative latent space where coordinates predict reservoir performance across copy memory, Lorenz system prediction, and sequential MNIST. Arrows indicate task-specific functional directions, and bottom values report cross-validated R2 .
Figure 2: Framework overview. (a) The conditional VAE encodes directed local circuits sampled from connectomes with node coordinates and cell type into a 32-dimensional structural latent coordinate, then decodes edge probabilities under node-level conditions. (b) Decoded graphs define fixed Dale-signed reservoirs evaluated on copy memory, Lorenz system prediction, and sequential MNIST; a predictor maps latent coordinates zi to functional scores Fi .
Figure 3: Connectome reconstruction fidelity and latent topological organization. (a) Representative binary adjacency matrices for the original connectome, the naive VAE baseline, and the full CVAE reconstruction. (b) PCA projection of the learned 32-dimensional latent vectors from 498 connectome samples, colored by clustering coefficient and mean degree. The smooth variation of these graph-level descriptors over the PC1–PC2 plane indicates that the latent space preserves continuous macroscopic topological variation.
Model
AUC ↑
Deg. ↓
Eff. ↓
Clust. ↓
Assort. ↓
Mod. ↓
Trans. ↓
Louv. ↓
Naive VAE baseline
0.722
0.084
0.288
0.820
1.681
0.326
0.733
0.392
Full CVAE
0.910
0.198
0.314
0.380
1.210
0.090
0.163
0.119
Table 1: Reconstruction comparison between the full CVAE and naive VAE baseline. AUC measures edge-level reconstruction. Deg., Eff., Clust., Assort., Mod., Trans., and Louv. denote mean degree, global efficiency, clustering coefficient, assortativity, modularity, transitivity, and directed Louvain modularity relative errors ( Blondel et al. 2008 ) .
Figure 4: Predicted versus measured reservoir performance from 32-dimensional CVAE latent coordinates for copy memory, Lorenz system prediction, and sequential MNIST.
Task
Mediator
rα,S
rS,F
βS∣α
βα
βα∣S
Red.
Copy
E/I recip.
0.787
0.698
0.267
0.757
0.547
0.278
Copy
CW5
0.527
0.508
0.151
0.757
0.678
0.105
Lorenz
Non-normality
-0.808
0.699
0.330
-0.723
-0.456
0.369
Lorenz
Top singular value
-0.785
0.649
0.213
-0.723
-0.556
0.231
sMNIST
Stable rank
0.904
0.912
0.319
0.945
0.657
0.305
sMNIST
Top-1 energy
-0.802
-0.825
-0.190
0.945
0.793
0.161
Table 2: Candidate partial mediators for task-specific latent functional directions. rα,S and rS,F are Pearson correlations for the two marginal associations. βS∣α is the mediator coefficient in the controlled model F∼α+S . βα is the coefficient of α in the total-effect model F∼α , and βα∣S is the coefficient of α after adding the mediator in F∼α+S . Red. is the fractional direct-effect reduction. CW5 denotes normalized length-5 closed walks.
Task
Linear
Ridge
SVM
XGB
GPR
Copy
0.624
0.625
-2.477
0.626
0.651
Lorenz
0.434
0.436
0.445
0.396
0.459
sMNIST
0.872
0.872
0.821
0.844
0.873
Table 3: Latent-to-function regression performance across reservoir-computing tasks. Entries report 5-fold cross-validated R2 from the 32-dimensional CVAE latent coordinates to the sample-level task metric. Higher values are better for all entries. XGB and GPR denote XGBoost and Gaussian process regression, respectively.
The aim of this work is to examine the connectomes of Caenorhabditis elegans through a computational lens using the reservoir computing framework. Connectomes are mappings of biological neural networks; C. elegans is the first organism for which physical connectomes covering the whole nervous system have been published. The connectomes of C. elegans used in this paper have been derived at different ages of the organism and are based on three different ways of measuring inter-cellular connections. They have, with minimal preprocessing, been implemented as reservoirs in the form of echo state networks, which are recurrent neural networks. In reservoir computing, the reservoir itself is not trained, rather the output of the reservoir is passed to a comparatively small read-out module in which training takes place. Training and testing is conducted in different neuro-inspired tasks, with the aim of using these tasks as a benchmark for the connectomes. This process has been repeated with different configurations of the reservoir and equally sized but randomized null models have been used for comparison. The results show that the biological wiring and a bio-informed configuration of input and output nodes of the reservoirs do not necessarily lead to better performance. Contrarily, the randomized null models are often outperforming the original connectomes on the chosen benchmarks. At the same time it becomes clear that the results depend a lot on the configuration of the reservoir and the way the connectome has been derived from the organism. Connectomes from different ages may produce varying outcome, without a clear trend becoming visible.
Felix S. Reimers, Ola Huse Ramstad, Aliaksandr Hubin +1
Department of Computer Science and Communication, Østfold University College, Halden, 1757, Østfold, Norway. · Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, ˚As, 1432, Akershus, Norway. · Department of Mathematics, University of Oslo, Oslo, 0316, Oslo, Norway.
Understanding how the architecture of neural networks shapes the computations they carry is a central challenge in neuroscience and machine learning. While specific circuit architectures have been linked to particular network computations and theoretical bounds on expressivity of broad classes of networks have been found, we are still missing general principles connecting the structure of finite networks to their computational capabilities. Here, we characterize the computational abilities of recurrent neural networks as a function of their connectivity by training a large collection of different networks to compute a large set of Boolean functions. For small networks, we constructed the complete ``catalogs'' of network-function performance, which revealed that computational capacity varies widely across architectures and that most networks show poor performance, and most functions are hard to compute. However, we show that having local 2- and 3-cycles in a network strongly enhances its computational ability, and networks with such cycles are often the minimal architectures that can solve particular functions. We further show that a small set of structural statistics accurately predict networks' performance. Extending our analysis to large networks showed that typical networks fail even to approximate a randomly selected function. Surprisingly, adding a small number of sparsely connected biologically-inspired interneurons to the network dramatically increases computational capacity. As in small networks, adding short cycles improved networks' capacity, outperforming acyclic or reachability-matched controls. Thus, our results identify local cycles as design principles linking neural connectivity to computational power, and offer a general framework to explore structure-function relations in computing networks.
Tom Talpir, Elad Schneidman
Department of Brain Sciences, Weizmann Institute of Science, Rehovot, Israel
Reservoir computing exploits the fixed dynamics of a recurrent network for temporal processing, requiring only a trained linear readout. Biological neural connectomes, shaped by millions of years of evolution, may encode computational structure beyond what random reservoirs provide, yet whether that structure can be further enhanced by principled optimisation remains an open question. We address it by applying four gradient-free, bio-inspired optimisers (Particle Swarm Optimisation, Differential Evolution, Grey Wolf Optimiser, and Whale Optimisation Algorithm) to the edge weights of connectome-based echo-state networks across six species spanning six orders of magnitude in neural complexity: C. elegans (279 neurons), Drosophila (49 nodes), mouse (112), rat (73), macaque (29 regions, continuous FLNe synaptic strengths), and human structural MRI connectivity (83 parcels). Each connectome is evaluated on four canonical reservoir computing benchmarks: Memory Capacity (MC), Lorenz attractor prediction, NARMA-10 system identification, and Mackey-Glass chaotic time-series prediction. All four optimisers consistently outperform unoptimised biological baselines across every task and species when initialised from biological weights. WOA achieves the largest gains on every task: up to a 17x MC improvement (C. elegans: 1.39 to 23.91) and up to 89% NRMSE reduction (Mackey-Glass, human), corresponding to an average 214% improvement across all species and tasks. Crucially, random initialisation on the same topology reliably underperforms biology, establishing biological weight values as an essential inductive bias that topology alone cannot recover. These results position bio-inspired, biologically-initialised optimisation as a principled and broadly effective strategy for connectome reservoir computing across the animal kingdom.
Anmol Guragain, Savvas Kakalis, Juan Ignacio Godino-Llorente
ETSI de Telecomunicación, Universidad Politécnica de Madrid, Madrid, Spain