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

Classification of 3-node restricted excitatory-inhibitory networks

We classify connected 3-node restricted excitatory-inhibitory networks, extending our previous paper (Aguiar et al., 2024). We assume that there are two node-types and two arrow-types, excitatory and inhibitory; all excitatory arrows are identical and all inhibitory arrows are identical; and excitatory (resp. inhibitory) nodes can only output excitatory (resp. inhibitory) arrows. The classification is performed under the following two network perspectives: ODE-equivalence and minimality; and valence ≤2. The results of this and the previous work constitute a first step towards analysing dynamics and bifurcations of excitatory-inhibitory networks and have potential applications to biological network models.

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

Understanding the dynamical behavior of complex systems from their underlying network architectures is a long-standing question in complexity theory. Therefore, many metrics have been devised to extract network features like motifs, centrality, and modularity measures. It has previously been proposed that network symmetries are of particular importance since they are expected to underly the synchronization of a system's units, which is ubiquitously observed in nervous system activity patterns. However, perfectly symmetrical structures are difficult to assess in noisy measurements of biological systems, like neuronal connectomes. Here, we devise a principled method to infer network symmetries from combined connectome and neuronal activity data. Using nervous system-wide population activity recordings of the C.elegans backward locomotor system, we infer structures in the connectome called fibration symmetries, which can explain which group of neurons synchronize their activity. Our analysis suggests functional building blocks in the animal's motor periphery, providing new testable hypotheses on how descending interneuron circuits communicate with the motor periphery to control behavior. Our approach opens a new door to exploring the structure-function relations in other complex systems, like the nervous systems of larger animals.

Fibration symmetry-breaking supports functional transitions in a brain network engaged in language

In his book 'A Beautiful Question', physicist Frank Wilczek argues that symmetry is 'nature's deep design,' governing the behavior of the universe, from the smallest particles to the largest structures. While symmetry is a cornerstone of physics, it has not yet been found widespread applicability to describe biological systems, particularly the human brain. In this context, we study the human brain network engaged in language and explore the relationship between the structural connectivity (connectome or structural network) and the emergent synchronization of the mesoscopic regions of interest (functional network). We explain this relationship through a different kind of symmetry than physical symmetry, derived from the categorical notion of Grothendieck fibrations. This introduces a new understanding of the human brain by proposing a local symmetry theory of the connectome, which accounts for how the structure of the brain's network determines its coherent activity. Among the allowed patterns of structural connectivity, synchronization elicits different symmetry subsets according to the functional engagement of the brain. We show that the resting state is a particular realization of the cerebral synchronization pattern characterized by a fibration symmetry that is broken in the transition from rest to language. Our findings suggest that the brain's network symmetry at the local level determines its coherent function, and we can understand this relationship from theoretical principles.

Dynamics and bifurcations in genetic circuits with fibration symmetries

Circuit building blocks of gene regulatory networks (GRN) have been identified through the fibration symmetries of the underlying biological graph. Here, we analyse analytically six of these circuits that occur as functional and synchronous building blocks in these networks. Of these, the lock-on, toggle switch, Smolen oscillator, feed-forward fibre and Fibonacci fibre circuits occur in living organisms, notably Escherichia coli; the sixth, the repressilator, is a synthetic GRN. We consider synchronous steady states determined by a fibration symmetry (or balanced colouring) and determine analytic conditions for local bifurcation from such states, which can in principle be either steady-state or Hopf bifurcations. We identify conditions that characterize the first bifurcation, the only one that can be stable near the bifurcation point. We model the state of each gene in terms of two variables: mRNA and protein concentration. We consider all possible 'admissible' models-those compatible with the network structure-and then specialize these general results to simple models based on Hill functions and linear degradation. The results systematically classify using graph symmetries the complexity and dynamics of these circuits, which are relevant to understand the functionality of natural and synthetic cells.

Classification of 2-node excitatory-inhibitory networks

We classify connected 2-node excitatory-inhibitory networks under various conditions. We assume that, as well as for connections, there are two distinct node-types, excitatory and inhibitory. In our classification we consider four different types of excitatory-inhibitory networks: restricted, partially restricted, unrestricted and completely unrestricted. For each type we give two different classifications. Using results on ODE-equivalence and minimality, we classify the ODE-classes and present a minimal representative for each ODE-class. We also classify all the networks with valence ≤2. These classifications are up to renumbering of nodes and the interchange of 'excitatory' and 'inhibitory' on nodes and arrows. These classifications constitute a first step towards analysing dynamics and bifurcations of excitatory-inhibitory networks. The results have potential applications to biological network models, especially neuronal networks, gene regulatory networks, and synthetic gene networks.

Fibration symmetry-breaking supports functional transitions in a brain network engaged in language

In his book 'A Beautiful Question' 1, physicist Frank Wilczek argues that symmetry is 'nature's deep design,' governing the behavior of the universe, from the smallest particles to the largest structures 1-4. While symmetry is a cornerstone of physics, it has not yet been found widespread applicability to describe biological systems 5, particularly the human brain. In this context, we study the human brain network engaged in language and explore the relationship between the structural connectivity (connectome or structural network) and the emergent synchronization of the mesoscopic regions of interest (functional network). We explain this relationship through a different kind of symmetry than physical symmetry, derived from the categorical notion of Grothendieck fibrations 6. This introduces a new understanding of the human brain by proposing a local symmetry theory of the connectome, which accounts for how the structure of the brain's network determines its coherent activity. Among the allowed patterns of structural connectivity, synchronization elicits different symmetry subsets according to the functional engagement of the brain. We show that the resting state is a particular realization of the cerebral synchronization pattern characterized by a fibration symmetry that is broken 7 in the transition from rest to language. Our findings suggest that the brain's network symmetry at the local level determines its coherent function, and we can understand this relationship from theoretical principles.

Fibration symmetries and cluster synchronization in the Caenorhabditis elegans connectome

Capturing how the Caenorhabditis elegans connectome structure gives rise to its neuron functionality remains unclear. It is through fiber symmetries found in its neuronal connectivity that synchronization of a group of neurons can be determined. To understand these we investigate graph symmetries and search for such in the symmetrized versions of the forward and backward locomotive sub-networks of the Caenorhabditi elegans worm neuron network. The use of ordinarily differential equations simulations admissible to these graphs are used to validate the predictions of these fiber symmetries and are compared to the more restrictive orbit symmetries. Additionally fibration symmetries are used to decompose these graphs into their fundamental building blocks which reveal units formed by nested loops or multilayered fibers. It is found that fiber symmetries of the connectome can accurately predict neuronal synchronization even under not idealized connectivity as long as the dynamics are within stable regimes of simulations.

Fibration symmetries and cluster synchronization in the Caenorhabditis elegans connectome

Capturing how the Caenorhabditis elegans connectome structure gives rise to its neuron functionality remains unclear. It is through fiber symmetries found in its neuronal connectivity that synchronization of a group of neurons can be determined. To understand these we investigate graph symmetries and search for such in the symmetrized versions of the forward and backward locomotive sub-networks of the Caenorhabditi elegans worm neuron network. The use of ordinarily differential equations simulations admissible to these graphs are used to validate the predictions of these fiber symmetries and are compared to the more restrictive orbit symmetries. Additionally fibration symmetries are used to decompose these graphs into their fundamental building blocks which reveal units formed by nested loops or multilayered fibers. It is found that fiber symmetries of the connectome can accurately predict neuronal synchronization even under not idealized connectivity as long as the dynamics are within stable regimes of simulations.

Breaking reflection symmetry: evolving long dynamical cycles in Boolean systems

In interacting dynamical systems, specific local interaction rules for system components give rise to diverse and complex global dynamics. Long dynamical cycles are a key feature of many natural interacting systems, especially in biology. Examples of dynamical cycles range from circadian rhythms regulating sleep to cell cycles regulating reproductive behavior. Despite the crucial role of cycles in nature, the properties of network structure that give rise to cycles still need to be better understood. Here, we use a Boolean interaction network model to study the relationships between network structure and cyclic dynamics. We identify particular structural motifs that support cycles, and other motifs that suppress them. More generally, we show that the presence of dynamical reflection symmetry in the interaction network enhances cyclic behavior. In simulating an artificial evolutionary process, we find that motifs that break reflection symmetry are discarded. We further show that dynamical reflection symmetries are over-represented in Boolean models of natural biological systems. Altogether, our results demonstrate a link between symmetry and functionality for interacting dynamical systems, and they provide evidence for symmetry's causal role in evolving dynamical functionality.

Fast algorithm to identify minimal patterns of synchrony through fibration symmetries in large directed networks

Recent studies have revealed the interplay between the structure of network circuits with fibration symmetries and the functionality of biological networks within which they have been identified. The presence of these symmetries in complex networks predicts the phenomenon of cluster synchronization, which produces patterns of a synchronized group of nodes. Here, we present a fast, and memory efficient, algorithm to identify fibration symmetries in networks. The algorithm is particularly suitable for large networks since it has a runtime of complexity O(Mlog⁡N) and requires O(M+N) of memory resources, where N and M are the number of nodes and edges in the network, respectively. The algorithm is a modification of the so-called refinement paradigm to identify circuits that are symmetrical to information flow (i.e., fibers) by finding the coarsest refinement partition over the network. Finally, we show that the algorithm provides an optimal procedure for identifying fibers, overcoming current approaches used in the literature.

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

BACKGROUND

Gene regulatory networks coordinate the expression of genes across physiological states and ensure a synchronized expression of genes in cellular subsystems, critical for the coherent functioning of cells. Here we address the question whether it is possible to predict gene synchronization from network structure alone. We have recently shown that synchronized gene expression can be predicted from symmetries in the gene regulatory networks described by the concept of symmetry fibrations. We showed that symmetry fibrations partition the genes into groups called fibers based on the symmetries of their 'input trees', the set of paths in the network through which signals can reach a gene. In idealized dynamic gene expression models, all genes in a fiber are perfectly synchronized, while less idealized models-with gene input functions differencing between genes-predict symmetry breaking and desynchronization.

RESULTS

To study the functional role of gene fibers and to test whether some of the fiber-induced coexpression remains in reality, we analyze gene fibrations for the gene regulatory networks of E. coli and B. subtilis and confront them with expression data. We find approximate gene coexpression patterns consistent with symmetry fibrations with idealized gene expression dynamics. This shows that network structure alone provides useful information about gene synchronization, and suggest that gene input functions within fibers may be further streamlined by evolutionary pressures to realize a coexpression of genes.

CONCLUSIONS

Thus, gene fibrations provide a sound conceptual tool to describe tunable coexpression induced by network topology and shaped by mechanistic details of gene expression.