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 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.

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.

Prediction of regional functional impairment following experimental stroke via connectome analysis

Recent advances in functional connectivity suggest that shared neuronal activation patterns define brain networks linking anatomically separate brain regions. We sought to investigate how cortical stroke disrupts multiple brain regions in processing spatial information. We conducted a connectome investigation at the mesoscale-level using the neuroVIISAS-framework, enabling the analysis of directed and weighted connectivity in bilateral hemispheres of cortical and subcortical brain regions. We found that spatial-exploration induced brain activation mapped by Fos, a proxy of neuronal activity, was differentially affected by stroke in a region-specific manner. The extent of hypoactivation following spatial exploration is inversely correlated with the spatial distance between the region of interest and region damaged by stroke, in particular within the parietal association and the primary somatosensory cortex, suggesting that the closer a region is to a stroke lesion, the more it would be affected during functional activation. Connectome modelling with 43 network parameters failed to reliably predict regions of hypoactivation in stroke rats exploring a novel environment, despite a modest correlation found for the centrality and hubness parameters in the home-caged animals. Further investigation in the inhibitory versus excitatory neuronal networks and microcircuit connectivity is warranted to improve the accuracy of predictability in post-stroke functional impairment.

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.

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.

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.

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.