Reaction-diffusion models in weighted and directed connectomes.

Connectomes represent comprehensive descriptions of neural connections in a nervous system to better understand and model central brain function and peripheral processing of afferent and efferent neural signals. Connectomes can be considered as a distinctive and necessary structural component alongside glial, vascular, neurochemical, and metabolic networks of the nervous systems of higher organisms that are required for the control of body functions and interaction with the environment. They are carriers of functional phenomena such as planning behavior and cognition, which are based on the processing of highly dynamic neural signaling patterns. In this study, we examine more detailed connectomes with edge weighting and orientation properties, in which reciprocal neuronal connections are also considered. Diffusion processes are a further necessary condition for generating dynamic bioelectric patterns in connectomes. Based on our precise connectome data, we investigate different diffusion-reaction models to study the propagation of dynamic concentration patterns in control and lesioned connectomes. Therefore, differential equations for modeling diffusion were combined with well-known reaction terms to allow the use of connection weights, connectivity orientation and spatial distances. Three reaction-diffusion systems Gray-Scott, Gierer-Meinhardt and Mimura-Murray were investigated. For this purpose, implicit solvers were implemented in a numerically stable reaction-diffusion system within the framework of neuroVIISAS. The implemented reaction-diffusion systems were applied to a subconnectome which shapes the mechanosensitive pathway that is strongly affected in the multiple sclerosis demyelination disease. It was found that demyelination modeling by connectivity weight modulation changes the oscillations of the target region, i.e. the primary somatosensory cortex, of the mechanosensitive pathway. In conclusion, a new application of reaction-diffusion systems to weighted and directed connectomes has been realized. Because the implementation was realized in the neuroVIISAS framework many possibilities for the study of dynamic reaction-diffusion processes in empirical connectomes as well as specific randomized network models are available now.

Connectome verification: inter-rater and connection reliability of tract-tracing-based intrinsic hypothalamic connectivity.

MOTIVATION: Structural connectomics supports understanding aspects of neuronal dynamics and brain functions. Conducting metastudies of tract-tracing publications is one option to generate connectome databases by collating neuronal connectivity data. Meanwhile, it is a common practice that the neuronal connections and their attributes of such retrospective data collations are extracted from tract-tracing publications manually by experts. As the description of tract-tracing results is often not clear-cut and the documentation of interregional connections is not standardized, the extraction of connectivity data from tract-tracing publications could be complex. This might entail that different experts interpret such non-standardized descriptions of neuronal connections from the same publication in variable ways. Hitherto, no investigation is available that determines the variability of extracted connectivity information from original tract-tracing publications. A relatively large variability of connectivity information could produce significant misconstructions of adjacency matrices with faults in network and graph analyzes. The objective of this study is to investigate the inter-rater and inter-observation variability of tract-tracing-based documentations of neuronal connections. To demonstrate the variability of neuronal connections, data of 16 publications which describe neuronal connections of subregions of the hypothalamus have been assessed by way of example. RESULTS: A workflow is proposed that allows detecting variability of connectivity at different steps of data processing in connectome metastudies. Variability between three blinded experts was found by comparing the connection information in a sample of 16 publications that describe tract-tracing-based neuronal connections in the hypothalamus. Furthermore, observation scores, matrix visualizations of discrepant connections and weight variations in adjacency matrices are analyzed. AVAILABILITY: The resulting data and software are available at http://neuroviisas.med.uni-rostock.de/neuroviisas.shtml.

Modeling and simulation for prediction of multiple sclerosis progression.

In light of extensive work that has created a wide range of techniques for predicting the course of multiple sclerosis (MS) disease, this paper attempts to provide an overview of these approaches and put forth an alternative way to predict the disease progression. For this purpose, the existing methods for estimating and predicting the course of the disease have been categorized into clinical, radiological, biological, and computational or artificial intelligence-based markers. Weighing the weaknesses and strengths of these prognostic groups is a profound method that is yet in need and works directly at the level of diseased connectivity. Therefore, we propose using the computational models in combination with established connectomes as a predictive tool for MS disease trajectories. The fundamental conduction-based Hodgkin-Huxley model emerged as promising from examining these studies. The advantage of the Hodgkin-Huxley model is that certain properties of connectomes, such as neuronal connection weights, spatial distances, and adjustments of signal transmission rates, can be taken into account. It is precisely these properties that are particularly altered in MS and that have strong implications for processing, transmission, and interactions of neuronal signaling patterns. The Hodgkin-Huxley (HH) equations as a point-neuron model are used for signal propagation inside a small network. The objective is to change the conduction parameter of the neuron model, replicate the changes in myelin properties in MS and observe the dynamics of the signal propagation across the network. The model is initially validated for different lengths, conduction values, and connection weights through three nodal connections. Later, these individual factors are incorporated into a small network and simulated to mimic the condition of MS. The signal propagation pattern is observed after inducing changes in conduction parameters at certain nodes in the network and compared against a control model pattern obtained before the changes are applied to the network. The signal propagation pattern varies as expected by adapting to the input conditions. Similarly, when the model is applied to a connectome, the pattern changes could give an insight into disease progression. This approach has opened up a new path to explore the progression of the disease in MS. The work is in its preliminary state, but with a future vision to apply this method in a connectome, providing a better clinical tool.

Connectome-based prediction of functional impairment in experimental stroke models.

Experimental rat models of stroke and hemorrhage are important tools to investigate cerebrovascular disease pathophysiology mechanisms, yet how significant patterns of functional impairment induced in various models of stroke are related to changes in connectivity at the level of neuronal populations and mesoscopic parcellations of rat brains remain unresolved. To address this gap in knowledge, we employed two middle cerebral artery occlusion models and one intracerebral hemorrhage model with variant extent and location of neuronal dysfunction. Motor and spatial memory function was assessed and the level of hippocampal activation via Fos immunohistochemistry. Contribution of connectivity change to functional impairment was analyzed for connection similarities, graph distances and spatial distances as well as the importance of regions in terms of network architecture based on the neuroVIISAS rat connectome. We found that functional impairment correlated with not only the extent but also the locations of the injury among the models. In addition, via coactivation analysis in dynamic rat brain models, we found that lesioned regions led to stronger coactivations with motor function and spatial learning regions than with other unaffected regions of the connectome. Dynamic modeling with the weighted bilateral connectome detected changes in signal propagation in the remote hippocampus in all 3 stroke types, predicting the extent of hippocampal hypoactivation and impairment in spatial learning and memory function. Our study provides a comprehensive analytical framework in predictive identification of remote regions not directly altered by stroke events and their functional implication.

The brainstem connectome database.

Connectivity data of the nervous system and subdivisions, such as the brainstem, cerebral cortex and subcortical nuclei, are necessary to understand connectional structures, predict effects of connectional disorders and simulate network dynamics. For that purpose, a database was built and analyzed which comprises all known directed and weighted connections within the rat brainstem. A longterm metastudy of original research publications describing tract tracing results form the foundation of the brainstem connectome (BC) database which can be analyzed directly in the framework neuroVIISAS. The BC database can be accessed directly by connectivity tables, a web-based tool and the framework. Analysis of global and local network properties, a motif analysis, and a community analysis of the brainstem connectome provides insight into its network organization. For example, we found that BC is a scale-free network with a small-world connectivity. The Louvain modularity and weighted stochastic block matching resulted in partially matching of functions and connectivity. BC modeling was performed to demonstrate signal propagation through the somatosensory pathway which is affected in Multiple sclerosis.

Towards Differential Connectomics with NeuroVIISAS.

The comparison of connectomes is an essential step to identify changes in structural and functional neuronal networks. However, the connectomes themselves as well as the comparisons of connectomes could be manifold. In most applications, comparisons of connectomes are applied to specific sets of data. In many studies collections of scripts are applied optimized for certain species (non-generic approaches) or diseases (control versus disease group connectomes). These collections of scripts have a limited functionality which do not support functional and topographic mappings of connectomes (hemispherical asymmetries, peripheral nervous system). The platform-independent and generic neuroVIISAS framework is built to circumvent limitations that come with variants of nomenclatures, connectivity lists and connectional hierarchies as well as restrictions to structural connectome analyses. A new analytical module is introduced into the framework to compare different types of connectomes and different representations of the same connectome within a unique software environment. As an example a differential analysis of the partial connectome of the laboratory rat that is based on virus tract tracing with the same regions of non-virus tract tracing has been performed. A relatively large connectional coherence between the two different techniques was found. However, some detected connections are described by virus tract-tracing only.

The connectome of the basal ganglia.

The basal ganglia of the laboratory rat consist of a few core regions that are specifically interconnected by efferents and afferents of the central nervous system. In nearly 800 reports of tract-tracing investigations the connectivity of the basal ganglia is documented. The readout of connectivity data and the collation of all the connections of these reports in a database allows to generate a connectome. The collation, curation and analysis of such a huge amount of connectivity data is a great challenge and has not been performed before (Bohland et al. PloS One 4:e7200, 2009) in large connectomics projects based on meta-analysis of tract-tracing studies. Here, the basal ganglia connectome of the rat has been generated and analyzed using the consistent cross-platform and generic framework neuroVIISAS. Several advances of this connectome meta-study have been made: the collation of laterality data, the network-analysis of connectivity strengths and the assignment of regions to a hierarchically organized terminology. The basal ganglia connectome offers differences in contralateral connectivity of motoric regions in contrast to other regions. A modularity analysis of the weighted and directed connectome produced a specific grouping of regions. This result indicates a correlation of structural and functional subsystems. As a new finding, significant reciprocal connections of specific network motifs in this connectome were detected. All three principal basal ganglia pathways (direct, indirect, hyperdirect) could be determined in the connectome. By identifying these pathways it was found that there exist many further equivalent pathways possessing the same length and mean connectivity weight as the principal pathways. Based on the connectome data it is unknown why an excitation pattern may prefer principal rather than other equivalent pathways. In addition to these new findings the local graph-theoretical features of regions of the connectome have been determined. By performing graph theoretical analyses it turns out that beside the caudate putamen further regions like the mesencephalic reticular formation, amygdaloid complex and ventral tegmental area are important nodes in the basal ganglia connectome. The connectome data of this meta-study of tract-tracing reports of the basal ganglia are available for further network studies, the integration into neocortical connectomes and further extensive investigations of the basal ganglia dynamics in population simulations.

neuroVIISAS: approaching multiscale simulation of the rat connectome.

neuroVIISAS is a generic platform which allows the integration of neuroontologies, mapping functions for brain atlas development, and connectivity data administration; all of which are required for the analysis of structurally and neurobiologically realistic simulations of networks. What makes neuroVIISAS unique is the ability to integrate neuroontologies, image stacks, mappings, visualizations, analyzes and simulations to use them for modelling and simulations. Based on the analysis of over 2020 tracing studies, atlas terminologies and registered histological stacks of images, neuroVIISAS permits the definition of neurobiologically realistic networks that are transferred to the simulation engine NEST. The analysis on a local and global level, the visualization of connectivity data and the results of simulations offer new possibilities to study structural and functional relationships of neural networks. This paper describes the major components and techniques of how to analyse, visualize and simulate with neuroVIISAS shown on a model network at a coarse CNS level (106 regions, 1566 connections) out of 13681 regions and 134043 connections of the left and right part of the CNS. This network of major components of the left and right hemisphere has small-world properties of the Watts-Strogatz model. Furthermore, synchronized subpopulations, oscillations of rate distributions and a time shift of population activities of the left and right hemisphere were observed in the neurocomputational simulations. In summary, a generic platform has been developed that realizes data-analysis-visualization integration for the exploration of network dynamics on multiple levels.

The intrinsic connectome of the rat amygdala.

The connectomes of nervous systems or parts there of are becoming important subjects of study as the amount of connectivity data increases. Because most tract-tracing studies are performed on the rat, we conducted a comprehensive analysis of the amygdala connectome of this species resulting in a meta-study. The data were imported into the neuroVIISAS system, where regions of the connectome are organized in a controlled ontology and network analysis can be performed. A weighted digraph represents the bilateral intrinsic (connections of regions of the amygdala) and extrinsic (connections of regions of the amygdala to non-amygdaloid regions) connectome of the amygdala. Its structure as well as its local and global network parameters depend on the arrangement of neuronal entities in the ontology. The intrinsic amygdala connectome is a small-world and scale-free network. The anterior cortical nucleus (72 in- and out-going edges), the posterior nucleus (45), and the anterior basomedial nucleus (44) are the nuclear regions that posses most in- and outdegrees. The posterior nucleus turns out to be the most important nucleus of the intrinsic amygdala network since its Shapley rate is minimal. Within the intrinsic amygdala, regions were determined that are essential for network integrity. These regions are important for behavioral (processing of emotions and motivation) and functional (memory) performances of the amygdala as reported in other studies.

Orexinergic innervation of the extended amygdala and basal ganglia in the rat.

The orexinergic system interacts with several functional states of emotions, stress, hunger, wakefulness and behavioral arousal through four pathways originating in the lateral hypothalamus (LH). Hundreds of orexinergic efferents have been described by tracing studies and direct immunohistochemistry of orexin in the forebrain, olfactory regions, hippocampus, amygdala, septum, basal ganglia, thalamus, hypothalamus, brain stem and spinal cord. Most of these tracing studies investigated the whole orexinergic projection to all regions of the intracranial part of the CNS. To identify the orexinergic efferents at the subnuclear level of resolution, we focussed on the orexinergic target in the amygdala, which is substantially involved in the LH output and contributes mostly to the functional outcome of the orexinergic system and the basal ganglia. Immunohistochemical identification of axonal orexin A and orexin B in male adult rats has been performed on serial sections. In the extended amygdala many new orexinergic targets were found in the anterior amygdaloid area (dense), anterior cortical nucleus (moderate), amygdalostriatal transition region (moderate), basolateral regions (moderate), basomedial nucleus (moderate), several bed nucleus of the stria terminals regions (few to dense), central amygdaloid subdivisions (dense), posteromedial cortical nucleus (moderate) and medial amygdaloid subnuclei (dense). Furthermore, the entopeduncular nucleus has been newly identified as another target for orexinergic fibers with a high density. These results suggest that subdivisions and subnuclei of the extended amygdala are specific targets of the orexinergic system.