Brain Struct Funct DOI 10.1007/s00429-014-0936-0
ORIGINAL ARTICLE
The connectome of the basal ganglia Oliver Schmitt • Peter Eipert • Richard Kettlitz Felix Leßmann • Andreas Wree
•
Received: 27 February 2014 / Accepted: 30 October 2014 Ó Springer-Verlag Berlin Heidelberg 2014
Abstract 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 O. Schmitt (&) P. Eipert R. Kettlitz F. Leßmann A. Wree Department of Anatomy, University of Rostock, Rostock, Germany e-mail:
[email protected] 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 graphtheoretical 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. Keywords Connectome Connectomics Basal ganglia Caudate putamen Striatum Substantia nigra Neuroontology Digital atlasing Tract tracing Multiscale Network analysis Graph analysis Abbreviations A All (all inputs and outputs) AD Average degree Ac Accumbens nucleus AC Amygdaloid complex AGl Lateral agranular prefrontal cortex AGm Medial agranular prefrontal cortex Aut Authoritativeness AW Average weight BC Betweenness centrality BG Basal ganglia C Circle CC Cluster-coefficient CE Closeness centrality
123
Brain Struct Funct
chain CL CM CNS CPu DG Dic Dii Dis DNN Doc Doi Dos DR EC Ent HIPP Hub I in INN L LGP LHb MDL MDM MDS MGP MRF O out PC PCA PF Pir PL PRC Pub RADin RADout Rec Rel Sic Sii Sis SG SNC SNR Soc Soi Sos
123
Chain pattern of a motif Centrolateral thalamic nucleus Central medial thalamic nucleus Central nervous system Caudate putamen Degree Direct input from contralateral Direct input from ipsilateral Direct input from ipsi- and contralateral Direct neighbor network Direct output to contralateral Direct output to ipsilateral Direct output to ipsi- and contralateral Dorsal raphe nucleus Eigenvector centrality Entorhinal cortex Hippocampus Hubness In (input to a region; used in tables) Symmetric input connection to a central node of a motif Indirect neighbor network Laterality Lateral globus pallidus Lateral habenular nucleus Mediodorsal thalamic nucleus lateral part Mediodorsal thalamic nucleus medial part Metric multidimensional scaling Medial globus pallidus Mesencephalic reticular formation Out (Output of region; used in tables only) Symmetric output connection from a central node of a motif Paracentral thalamic nucleus Principal component analysis Parafascicular thalamic nucleus Piriform cortex Path length Page rank centrality Number of articles Radiality of the input Radiality of the output Reciprocal Reliability Subtree input from contralateral Subtree input from ipsilateral Subtree input from ipsi- and contralateral Subgraph centrality Substantia nigra compact part Substantia nigra reticular part Subtree output to contralateral Subtree output to ipsilateral Subtree output to ipsi- and contralateral
SP SPN STh VA VL VM VTA
Length of shortest path Spiny neurons of the CPu Subthalamic nucleus Ventro anterior thalamic nucleus Ventrolateral thalamic nucleus Ventromedial thalamic nucleus Ventral tegmental area A10
Introduction The basal ganglia (BG) are composed of the four major nuclei neostriatum, globus pallidus (GP), substantia nigra (SN) and subthalamic nucleus (STh). The canonical basal ganglia system consists of the BG and of those cortical areas that are connected with BG (Gerfen and Bolam 2010). In rodents the neostriatum is a single nucleus, the caudate putamen (CPu) which is divided by the internal capsule into the caudate nucleus and putamen in higher vertebrates. The GP consists of an external segment or lateral globus pallidus (LGP) and the internal segment or medial globus pallidus (MGP). The MGP in rodents is also called the entopeduncular nucleus. The SN consists of the pars compacta (SNc) and the pars reticulata (SNr). These nuclear complexes are related to the dorsal aspects of the BG. The ventral aspects contain the accumbens nucleus (Ac), the ventral pallidum and the medial aspects of the STh and SN. The dorsal aspect of the basal ganglia has motor and associative functions, whereas the ventral aspect is related to limbic functions. Additional motor components that are strongly associated with the dorsal aspect of the BG are the frontal cortex especially the primary motor cortex (lateral agranular prefrontal cortex) and subnuclei of the thalamus (Gerfen 2004). Principal afferents of the BG arise from the ipsi- and contralateral cerebral cortex (layer 5 glutamatergic neurons), intralaminar thalamic nuclei such as the centromedian and parafascicular nucleus, the dorsal raphe nucleus (DR) and the amygdaloid complex (AC). The most extensive input to the BG is glutamatergic from the cerebral cortex and thalamus. These projections terminate in the heads of spines of the GABAergic spiny neurons (SPN) of the CPu. The latter are the principal projection neurons innervating the SN and GP. These SPN possess a dense local axon collateral system that innervates other SPN and interneurons of the CPu. At the level of macroconnectomes and macrocircuits all three major pathways of the BG are represented by the connectome data in this study. The direct pathway is an excitatory pathway which has the same source region (CPu) and BGoutput region (STh) as the inhibitory indirect pathway. The macrocircuitry of the BG is dominated by these two
Brain Struct Funct
principal pathways (direct, indirect) (Alexander and Crutcher 1990; DeLong 1990; Smith et al. 1998) and a hyperdirect pathway (Kita and Kitai 1994) which will be investigated in the introductory part of the ‘‘Results’’. These pathways along with the cerebral cortex and thalamus form three loops (cortical-basal ganglia loops). It is known that the GABAergic neurons in the LGP influence the activity at every level of the BG by their extensive axon collaterals (Tepper et al. 2007). Single cell labeling studies (Kita and Kitai 1994; Bevan et al. 1998; Sato et al. 2000) have indicated that LGP projection neurons have axonal collaterals to the STh and BG output nuclei. Additionally, up to one-third of LGP neurons project to the CPu (Kita and Kitai 1994; Bevan et al. 1998) and to its GABAergic interneurons which innervate the spiny projection neurons (Kita 1993; Koos and Tepper 1999; Bennett and Bolam 1994). Besides the major inputs to and outputs from the BG nearly all areas of the neocortex are interconnected with regions of the BG (Gerfen and Bolam 2010). In terms of functional segregation, frontal cortical areas and intralaminar thalamic nuclei are involved in planning and execution of movement behavior. Outputs of the BG to the intermediate layers of the superior colliculus are involved in the generation of eye and head movements and output to the pedunculopontine nucleus is necessary for orienting movements of the body. A functional and hodological grouping of subregions the BG into four categories has been proposed (Parent 1986; Alexander and Crutcher 1990; DeLong 1990; Smith et al. 1998; Gerfen and Bolam 2010). The input nuclei are the CPu and STh which receive cortical inputs. The output nuclei are the MGP and SNR which project to the thalamus, midbrain and brainstem. The LGP is a relay nucleus connecting the input and output nuclei. The modulator nucleus is the SNC. The major objectives of this study are 1. 2. 3. 4.
the collation of connectivity data, a weighted and directed analysis of the BG connectome, the analysis at two levels of detail and the intrinsic, extrinsic, and differential analysis
which are specified in the following. In addition to this general organization of the BG abundant information of interconnections and densities of connections between topographic, functional, cytoarchitectonic and chemoarchitectonic subdivisions of BG regions (intrinsic BG connections) is available. More specifically, afferent, efferent intrinsic and extrinsic connections of the BG regions are described in 789 reports (Lanciego and Wouterlood 2011) (objective 1). These connectivity data can be collated for a meta-study as shown by others (Bakker et al. 2012; Sporns et al. 2000b; Stephan et al. 2001b; Press et al. 2001) to
analyze and characterize the resulting network in terms of graph theory (Sporns 2002, 2011; Rubinov and Sporns 2010; Young et al. 1994, 1996; Young 1992b). One particular objective of this study is a weighted and directed analysis of the BG connectome with regard to strengths of connections in combination with ipsi- and contralateral connectivity (objective 2). Connectivity data sets of parts of nervous systems have been developed and analyzed by several groups (Felleman and Essen 1991; Young 1992a, 1993; Sporns et al. 2000a, 2002, 2004; Sporns and Ko¨tter 2004; Honey et al. 2007; Modha and Singh 2010; Sugar et al. 2011), however, these data do not provide a complete collection of BG connections. Most of the connections in these data sets are compiled into neocortical connectomes within the framework of meta-studies of tract-tracing publications. The sources of information in this study originate exclusively from peerreviewed tract-tracing publications where anterograde and/ or retrograde tracers were applied (lesion studies were not evaluated). Only those publications were considered that describe connectivity in juvenile or adult rats. Connections between neuronal regions are handled with the help of a consistent neuroontology (Schmitt and Eipert 2012) that can be updated frequently in parallel to the fast progress of the identification of connections in tract-tracing publications. In most tract-tracing studies that describe new details of afferents and efferents, modifications of nomenclatures and overlap of regions occur which must be interpreted and integrated into connectome projects (Schmitt and Eipert 2012). Therefore, it is necessary to propose refinements of subdivisions of regions that are documented in stereotaxic atlases (Paxinos and Watson 2004, 2009; Swanson 2003). Such subdivisions can be presented in acyclic graphs or hierarchies. This approach allows to expand and aggregate regions and to apply such a selection of super- and subregions for defining a connectome and partial connectomes, respectively, without loosing non-standard terminologies of publications. This technique was applied here to define a core connectome of the basal ganglia (BG1) and a more detailed connectome that includes especially limbic connections (BG2) (objective 3) (the composition of the BG1and BG2-networks is described in the ‘‘Materials and methods’’ as well as in the ‘‘Conditional connectomics’’). The BG-connectome was analyzed in three contexts (objective 4): (1) embedded in regions that directly connect to it (Vlachos et al. 2012) (connectome embedding), to measure extrinsic connectivity; (2) stand-alone, to measure intrinsic connectivity (3) differential, to compare the role of local regions in the other two contexts (differential connectomics). A further aim of this investigation was to collate the connectivity data (Burns and Cheng 2006) of all peer-
123
Brain Struct Funct
reviewed non-viral tract tracing publications (Bohland et al. 2009a) in which connections of regions of the BG were described. Viral tract tracing was not considered in this work because some neurotrophic viruses propagate transsynaptically in a time-dependent manner. The density and reliability of collated connections should be as large as possible. Such a BG data source has been proposed recently (Looi et al. 2013), but was not available, so far. Now, this BG connectome is available within a database (http://neuroviisas.med.uni-rostock.de/connectome/index. php) and can be downloaded for sharing (Van Horn and Gazzaniga 2013; Koslow 1997, 2005; Ko¨tter 2003), (https://neuinfo.org-/datasharing/index.shtm). Based on these tract-tracing data, the intrinsic, extrinsic, bilateral and weighted connectivity of the BG connectome was characterized in terms of global and local parameters, as well as multivariate statistics to determine regions of specific importance of intrinsic and extrinsic BG connectomes of the rat brain. By applying modularity analysis it was intended to investigate a possible correlation of functionally similar regions and their connectivity patterns. A population-based large-scale simulation of a thalamocortical network of the whole human cortex (Izhikevich and Edelman 2008) has provided evidence for the propagation of an excitatory wave from the visual cortex to frontal areas following a spike injection into the occipital areas. Because macroscopic connectivity in this model was primarily based on tractographic data of DTI-measurements, one assumption of this model was a complete reciprocity of all macroscopic interconnections of cortical regions. The hypothesis that reciprocal connections in tract tracing studies are highly specifically distributed has been investigated. The result of this analysis of reciprocal connections is of particular interest for modeling connectivity at macroscopic levels based on DTI-measurements and tractographic analysis.
Materials and methods In a meta-analysis of 4,513 peer reviewed tract-tracing publications describing non-viral tract-tracing experiments in juvenile and adult rats we found 3,265 articles that provide connectivity data that are useful for a mesoscale connectome analysis (Schmitt et al. 2012b; Burns and Cheng 2006; Burns et al. 2008a, b; Bota and Swanson 2007b, 2010). Nearly 800 of these 3,265 publications are considered as the data core because they contain connectivity data of particular regions of the basal ganglia. The inquiry of these publications was performed through PubMed (http://www.ncbi.nlm.nih.gov/pubmed) using search expressions (Schmitt et al. 2012b). Among the 3,265 articles numerous studies were found that applied tract tracing
123
techniques for verifying and confirming experiments in adult laboratory rats as controls. Also these studies were taken into consideration, however, only those connections of them were collated which were observed in control rats. The sources, targets, semiquantitative weights (strengths of connections), type of tracer, directions of transport and features of somata (transmitters, enzymes, receptors) were transferred to tables. The raw data were verified by at least two experts in neuroanatomy. The verified data tables were imported into neuroVIISAS. This software is a platformindependent generic framework (Schmitt and Eipert 2012) developed in JAVATM for advanced connectome analysis (Bullmore and Sporns 2009; Stam and Reijneveld 2007; Albert and Barabasi 2002; Sporns et al. 2000a, 2004), visual analytics (Ware 2013; Archambault et al. 2008; Wong et al. 2006a, b) and nervous system specific 3Dmultiscale visualizations using stereotaxic coordinate systems (Hjornevik et al. 2007; Gustafson et al. 2004, 2007; Ju et al. 2006; Bjaalie 2002). The advantage of neuroVIISAS is that different modalities like ontologies, digital atlasing, connectomics and neuron population simulation (Feng et al. 2005; Burns et al. 2006; Gewaltig and Diesmann 2007; Moore et al. 2007; Niggemann et al. 2008; Lee et al. 2008; Gouws et al. 2009) are integrated into one flexible framework (Martone et al. 2004). The neuroVIISAS installation package can be downloaded from (http://neu roviisas.med.uni-rostock.de). All computations and visualizations were realized exclusively with the neuroVIISAS framework (chord diagrams were generated using the circular layout engine CIRCOS (http://circos.ca) through the CIRCOS–neuroVIISAS-interface). In most reports, connections are described in the form of ordinal categories (1, 2, 3, 4), symbols (*, **, ***) or expressions like: sparse, moderate, dense. Occasionally, authors have performed a stereologic analysis of retrogradely marked perikarya or a densitometric quantification of axonal terminals. There exist no standardization or normalization between these ordinal categories and estimations of the density of connections between independent studies or different reports. Obviously, connection weights are not discrete categories (Hilgetag and Grant 2000) and extend over several orders of dimensions (Markov et al. 2011, 2014). The latter determined a lognormal distribution of the fraction of neurons following injections of different tract tracers. Lognormal distributions of corticocortical connectivity was also found by Oh et al. (2014), ErcseyRavasz et al. (2013) and Wang et al. (2012). Therefore, ordinal categories were transformed to a logarithmic scale. This was realized in neuroVIISAS by a transformation table to approximate the distribution of the connection strengths logarithmically (Appendix Fig. 10). We would like to point out that these estimated connections strengths should be considered as a hypothetical estimate.
Brain Struct Funct
In general, connections are described in tract tracing studies in terms of ipsi- and contralateral courses rather than by providing exact stereotaxic coordinates (The sagittal axes in the rat atlas of Paxinos and Watson (2007) have the same values for the left and right side. In the rat atlas of Swanson (2003) one hemisphere is considered, only.) or an explicit indication of a left and/or right side injection into a specific region. In neuroVIISAS it is possible to manage injections and labelled sites either with explicit laterality, or with a symmetric representation where the labelled sites are either ipsi- or contralateral. To judge the reliability R of the data of connections, a reliability estimation of connectivity data is introduced. Because the connection data and not the measurement of connections itself constitute R it should be considered as an observation score. If this observation score is large then the probability that the connection really exists is large. An important assumption for calculating R is, that at least most of the tract tracing reports of a set of regions of interest have been evaluated and connectivity data are available for computation. The observation score of a connection is estimated by adding reliability weights. Reliability weights (Table 1) are defined for Table 1 Reliability weights used for estimating the reliability parameter R Variable
Case
Value
t
a/r
0.25
t
r
0.5
t
a
0.5
t
r?a/r
0.7
t
a?a/r
0.7
t
a?r
1.0
t w
a?r?a/r 3:0 unknown
1.0 0.7
w
2:0 fibers of passage
0.0
w
1:0 not clear
0.8
w
0:5 exists
w
0.0 not present
w
0.5 very light
1.0
w
1.0 light/sparse
1.0
w
1.5 light/moderate
1.0
w
2.0 moderate/dense
1.0
w
2.5 moderate/strong
1.0
w
3.0 strong
1.0
w
4.0 very strong
1.0
0.9 -1.0
Value reliability weight of connection strength, t variable of reliability weight for transport directions of tracers, w variable of reliability weight for strengths of connections, a anterograde tracer transport, r retrograde tracer transport, a/r bidirectional tracer transport, a?r?a/ r means that a connection has been proved by an anterograde, a retrograde and a bidirectional transported tracer
–
–
types of tracer transport directions (anterograde, retrograde) (t weight of the transport direction of a tracer) and the weight or strength of a connection (w weight of the connection strength).
The reliability R is defined as follows: X X R¼ wþ þ tþ þ w t
ð1Þ
Identical connections (connections between identical regions) may be described severalfold in a particular report in different animals or ‘‘cases’’ and/or in different reports. Since all these data are available in our database we are able to calculate sums of w of such case-based and noncase-based experimental observations. However, if a connection has been proved only by a retrograde or anterograde method then a smaller t ¼ 0:5 is assigned instead of a prove by a anterograde and retrograde method in two independent experiments (t ¼ 1). A bilateral transport (a/r) of a tracer is considered also in this weighting scheme (Table 1). The variables have subscripts and þ which indicate an observation of an existing connection (wþ , tþ ) or an explicit description that a connection does not exist (w , t ). Hence, explicitly not existent connections are weighted by negative reliability weights. The number of P observations of a specific connections is added ( w). However, the value of t is determined by identifying different directions of tracer transports within all connections that were added. For example, a specific connection that was observed in 10 different tract tracing experiments by applying anterogradely transported tracers obtains a w ¼ 10 and a t ¼ 0:5 (R ¼ 10:5). If the connection that was observed 10 times was found using 7 anterograde and 3 retrograde tracers then t ¼ 1 (R ¼ 11). And if the connection is observed 9 times (wþ ¼ 9) with 7 retrograde and 2 anterograde tracers (tþ ¼ 1) and in addition this specific connection has not been found in one experiment (w ¼ 1) with an anterograde tracer (t ¼ 0:5) then R ¼ 8:5. If there are many descriptions in publications that a specific connection does not exist, it gets a strong negative observation score. Finally, these scores are presented in a reliability matrix. We must emphasize that the observation score is defined for the data of a connection between a specific source or efferent region and a target or afferent region. So, it is not a score for connections rather than for the underlying data of the connections. Some connections and regions may receive more attention in tract tracing research so that frequently investigated or popular regions may have higher observation scores. This aspect is explicated in ‘‘Discussion’’. Regions of the nervous system of this particular project of the rat connectome are organized hierarchically (Schmitt and Eipert 2012), in graph theory such a hierarchy is called
123
Brain Struct Funct
a tree. An example of a hierarchy is shown in the Appendix Fig. 11. The root region of this hierarchy can be the nervous system of the rat. If a region is divided into smaller parts, these subregions are the childs of this parent region, which could be further subdivided. A region and all its descendents give rise to a so called subtree. Regions without descendents are called leafs. With regard to this concept the terms direct input or direct output of a region sum up all those connections that were described to end at (afferents), respectively, start from (efferents) this specific region. The terms subtree input or subtree output specify all those connections that end in, respectively, start from any region of the subtree of a given region including the region itself. The BG1 network consists of 10 regions (SNC, SNR, VA, VL, LGP, MGP, CPu, STh, AGl, AGm). These core regions of the BG are involved in the processing of motor functions. There exist several additional regions which are strongly interconnected with the BG1 regions, however, they are merely indirectly connected with regard to the principal BG pathways (direct, indirect, hyperdirect). To obtain a more comprehensive network, conditions for the selection of further regions can be applied (conditional connectomics). These selected regions belong to specific thalamic nuclei (lateral habenula among other regions), mescencephalic nuclei like the ventral tegmental area as well as to the limbic system like the ventral striatum (ventral CPu, accumbens nucleus). Further regions like the piriform cortex, hippocampus, entorhinal cortex and mesencephalic reticular formation (MRF) were selected in combination with their thalamic relay nuclei. In the following we will refer to this network as the BG2 network. The BG2 network is more extensive and contains with regard to the BG1 regions 15 further regions (VTA, MRF, CM, CL, PC, PF, MDL, MDM, VM, LHb, AC, Ac, HIPP, Ent, Pir). These two network configurations of the BG can be considered in a wider sense as a motoric BG connectome (BG1) and an associative BG connectome (BG2). The 10 regions of the BG1 connectome are localized on levels 9–14 (most on level 12) and the regions of the BG2 connectome on levels 8–14 of the hierarchy. The complex spatial arrangement of regions (Fig. 1) and their connections (Fig. 2) becomes visible by applying interactive 3Dvisualization of neuroVIISAS. The BG1 network was embedded in a surrounding network of regions that all have at least one input and one output connection to a region of the BG1 network. All regions of the surrounding network are located at a comparable hierarchical level like the BG1 region. The regions or parts of the selected regions of the BG1 or BG2 network do not overlap with or are not included in other regions (Bohland et al. 2009b). A further hypothesis to be tested in this study is, that specific distributions of reciprocal connections could
123
characterize the three principal pathways of the BG. Such specific distributions were analyzed by a recently developed pathway analysis of the BG-connectome data. Moreover, it was tested if additional routes from a source to a target region may exist in the directed BG connectome, a phenomenon which may entail a routing problem in structural connectomics. To determine pathways from source to target regions a sequential search of interconnected regions (pathway analysis) was applied in neuroVIISAS. The search-technique allows to add individual regions that should be passed in pathways from the provided source to the target regions. The result is shown in a table that contains weights, reliability and reciprocal connections within the detected pathways. To share the connectivity data of the BG, the rat connectome database can directly be queried through a web interface (http://neuroviisas.med.uni-rostock.de/connectome/ index.php).
Results Principal pathways The well-known principal pathways of the BG can be identified in the connectome. The adjacency matrix of the core regions of the BG is shown in Fig. 3a (BG1 network). Each of the following pathways can be derived from this matrix. The direct excitatory pathway [AGl, AGm] ! [CPu] ! [SNR, MGP] ! [VA, VL] ! [AGl, AGm] of the BG can be reconstructed using the adjacency matrix of principal regions (DeLong 1990; Alexander and Crutcher 1990; Albin et al. 1989) (Fig. 4a). The indirect inhibitory pathway [AGl, AGm] ! [CPu] ! [LGP] ! [STh] ! [SNR, MGP] ! [VA, VL] ! [AGl, AGm] of the BG can also be built using the same selection of regions (DeLong 1990; Alexander and Crutcher 1990; Albin et al. 1989) (Fig. 4b). The hyperdirect inhibitory pathway [AGl, AGm] ! [STh] ! [SNR, MGP] ! [VA, VL] ! [AGl, AGm] of the BG can be built by means of the same selection of regions (Bosch et al. 2012; Nambu et al. 2002) (Fig. 4c), too. The number of tract-tracing experiments that describe a certain connection is indicated by the number above each edge. Hence, the connection of the SNC to the CPu is described in 44 tract-tracing publications (Fig. 4a–c). The adjacency matrix of the ipsi- and contralateral BG network suggests that there could be additional pathways from AGl over the CPu back to AGl. Therefore, a pathway analysis was applied to detect additional routes or parallel pathways in the BG network (Table 2). If the maximum pathlength is
Brain Struct Funct Fig. 1 The regions of the BG2 network in stereotaxic space. a The stereotaxic space of the regions of the BG2 network. This view shows the region from dorsal, respectively, from top. Cortical regions are shown in red, subcortical parts in light red, thalamic regions are light magenta and SN are presented in dark magenta. The same color scheme was used for the coloring of abbreviations of regions in matrices. b Abbreviations of the regions of one hemisphere. c View from rostral directly onto the accumbens nucleus and ventral striatum. d View from caudal where midbrain structures and the turquoise MRF appear
restricted to 4 (because the minimum length of the principal pathways is 4), 653 ipsilateral pathways exist from AGl to AGl (circular pathway) in the BG2 network. If the CPu is determined as the first region in a pathway from the AGl back to the AGl, then 74 ipsilateral pathways (8 of these pathways have less than 4 edges) are found. The latter were sorted with regard to the reliability value (see ‘‘Reliability of connections’’). 6 pathways have reliability values that are larger than 10 and all of them have the SNR as the third region within the pathway. Most variability was found for thalamic regions (VM, PF, VL, PC, CM, CL) as output region to AGl. The direct pathway of this group has a reliability of 12. The indirect pathway may use MDL and PF as thalamic output regions to AGl or the VTA as an efferent region to AGl (reliability values are lying between 4 to 5). The hyperdirect pathway has a reliability of 7.5 and may use different output regions (MDL, PF, VM, VTA) to the AGl, too. Besides variations of BG-output nuclei of the thalamus, there are principal changes of input nuclei within these pathways, e.g. AGl has projections to AGm, CL, MDL, MGP, MRF, PC, PF, SNR, VL and VM. The sequence of connections of directly connected regions within a pathway were examined with regard to reciprocal connections. Then the patterns of reciprocity of connections within all principal pathways of the BG (only those regions were considered that were selected in the
BG2-network) were analyzed. In the first part, MGP was selected as a connecting region from the BG to thalamic and mesencephalic regions that have connections to AGl. In the second part SNR was selected as a connecting region from the BG to thalamic and mesencephalic regions that have connections to AGl. Interestingly, it was found in the first part of the analysis that PF is the sole region with reciprocal connections to MGP and AGl within the output of the BG2 network in all three principal BG-pathways. There were no reciprocal connections reported from VA, VL, VM and MDL to MGP. However, MGP seems to have a reciprocal connection with VTA, however, VTA is nonreciprocally connected with AGl. In the second part of this reciprocal pattern analysis with regard to SNR it was found that in addition to PF also VM has reciprocal connections to SNR and to AGl in all three principal BG pathways. Again VTA has a reciprocal connection to SNR and the non-reciprocal connection to VTA in this SNR BG-output pathway. Contralateral connectivity The core regions of the basal ganglia show larger connection strengths in ipsi- than in contralateral connections. The AGl and AGm possess more contralateral connections than the core regions of the BG (Fig. 3b). Reciprocal
123
Brain Struct Funct Fig. 2 The weighted connections and color-coded weighted DGAll parameter of the bilateral BG2 network. a The ipsilateral connections of the BG2 network are shown using a 3-axis expansion technique of neuroVIISAS. b A unilateral view with labels of regions. c Bilateral regions in a comparable view like in (a) with labels
connections of the BG core regions occur exclusively ipsilaterally, however, two reciprocal connections of the contralateral AGl with VL and CPu were identified. Ipsiand contralateral connections are visualized in 3D in Fig. 2. Stronger and weaker ipsi- and contralateral connectivity of AGl is clearly visible in the low-resolution BG1-network visualization using the chord diagrams of CIRCOS (Appendix Fig. 12). However, the CIRCOS visualization is not suitable for the high-resolution BG2 network (Appendix Fig. 13). Reliability of connections The observation scores of the BG network are presented in Fig. 3c. The most reliable connection in the BG network is the ipsilateral SNC ! CPu projection (43.6) followed by
123
the ipsilateral SNR ! CPu (24.1) projection. The largest reliability of a contralateral connection is observed for SNC ! CPu (5). Further, large reliabilities exist for ipsilateral efferences of STh to SNR (14) and MGP (10). Global and local connectome analysis of the BG2 network The ipsilateral BG2 network consists of 25 regions and 235 connections and the adjacency matrix has a line density of 39.167 %. The connections of thalamic regions to mesencephalic and diencephalic regions of the adjacency matrix are sparse, a phenomenon that is well presented by the adjacency matrix (Fig. 3d) as well as in the distance matrix that show relatively larger distances from these regions to other regions of the BG network (Appendix Fig. 14). The
Brain Struct Funct
a
b
c
d
Fig. 3 Matrix representations of the BG connectome. a Adjacency matrix of the core regions of the motoric basal ganglia. Numbers and gray values are indicating the number of experiments that describe connections. Source regions are shown in the rows and target regions in the columns. b The bilateral connectivity of the core regions of the motoric basal ganglia. The weighted adjacency matrix displays the mean weights of connections. Weights are ranked as shown in the color table. Altogether, the ipsilateral connectivity is denser than the contralateral. Relative strong connections exist between contralateral cortical regions, only. c The reliability of the weighted connectivity of the core regions of the motoric basal ganglia. The reliability of the weighted bilateral connections is calculated using the number of
experiments that describe connections (if different experiments within a particular publication occur than all experiments are considered) and variability of weights. The maximum is 43.6 for SNC ! CPu and the minimum is 0.1. d The adjacency matrix of the ipsilateral BG2 network with regions that are not primarily involved in motoric processing (BG2). The intrinsic connectivity of thalamic regions (magenta) is very weak. The regions of the midbrain and MRF have abundant connections to most other part of the BG2 network. Especially SNR has strong connections to nearly all thalamic regions in this selection. The amygdalar complex and the CPu have comparable number of output connections. The variability of weights of the bilateral BG2 network is shown in the Appendix Fig. 17
central connectional position of the CPu and the substantial contralateral connections have been visualized using an organic planar graph visualization (Fig. 4d). The parafascicular thalamic region has a conspicuously dense connectivity to other regions. This is confirmed by computing the communicability (Estrada and Hatano 2008) (Appendix Fig. 15). Communicability is a measure for the number of paths from a source to a target region. Here, PF has significantly larger values than other thalamic regions. Furthermore, the LGP has relatively small communicability values. A related matrix is the generalized topological overlap matrix GTOM (Yip and Horvath 2007). It is a
measure of pairwise interconnectedness that is proportional to the number of neighbors that a pair of nodes share in common. Interestingly, VTA and MRF possess a noticeable amount of common neighbors resulting in relative large GTOM values (Appendix Fig. 16). The global analysis of the BG2 network reveals a network architecture close to a small-world network because the small worldness is 1.258. The average cluster coefficient of 0.661 is larger than in the Erdo¨s Re´nyi, Watts– Strogatz or Barabasi–Albert randomizations. The scale-free property (power law distribution of edges) of the BG2 network can be asserted because with 2.2, the deviation D
123
Brain Struct Funct Fig. 4 The three principal pathways of the motoric BG. a Direct pathway. b Indirect pathway. c Hyperdirect pathway. The levels of the regions within the hierarchy (the hierarchical organization of regions is shown in Appendix Fig. 11) are shown at the left side of the graphs. The number of publications of specific parts of a pathway are displayed by the numbers on the connections and the maximal weight by the same colors as shown in Fig. 3b. Thick lines indicate reciprocal connections. d The bilateral BG2 network with weighted edges in a planar graph visualization with an organic edge layout. The CPu is localized around the center of regions of each hemisphere. _L region of the ipsilateral hemisphere, _R region of the contralateral hemisphere
a
b
c
d
from the power function is relatively small. The mean scale-free property as well as mean global values of the rewiring randomization resemble the real BG2 network values more closely, than those of other randomizations (Appendix Fig. 18). The weighted modularity analysis (Hagmann et al. 2008) distinguishes 5 modules (Fig. 5). In module 1 most regions are strongly interconnected with the CPu. Module 2 contains the primary motor cortex and important core regions of the BG (SNR, STh, LGP) as well as specific thalamic regions and MRF. Module 3 and 5 consist of regions that have memory functions and are involved in emotional processing. Module 4 contains MGP and VL (VA has light connections only). Based on the
123
weighted modularity analysis modules 1, 2 and 4 contain regions that are associated with motoric function and modules 3 and 5 those regions involved in memory and emotional processing. Hence, the computed modules partly reflects the functional separation of regions. Some parameters of the local-network analysis of the BG2 network are summarized in Table 6 (a complete list of all local parameters of this analysis is shown in Appendix Table 15). After sorting the regions by DGall , it turns out that AC has most inputs and outputs (42), followed by the CPu (41) and VTA (40). If the analysis is done on the weighted rather than binarized data, the top three regions are CPu, VTA and SNR (Appendix Table 7). The
Brain Struct Funct Table 2 Pathway analysis starting with projections from the AGl to the CPu Path
PL
Rec
AW
Pub
Rel
AGl ! CPu ! SNR ! VM ! AGl
4
1
3
13
13.13
AGl ! CPu ! SNR ! PF ! AGl
4
1
3
12
12.95
AGl ! CPu ! SNR ! VL ! AGl
4
0
3
11.75
12.02
AGl ! CPu ! SNR ! PC ! AGl
4
0
3
11.25
11.73
AGl ! CPu ! SNR ! CM ! AGl
4
0
3
11
11.48
AGl ! CPu ! SNR ! CL ! AGl
4
0
3
10.75
10.98
AGl ! CPu ! SNR ! MDL ! AGl
4
0
2.75
10
9.75
AGl ! CPu ! SNR ! VTA ! AGl
4
0
2.25
9.25
9.75
AGl ! CPu ! SNR ! VA ! AGl
4
0
2.25
9
9.30
AGl ! CPu ! SNC ! CM ! AGl
4
0
3
7.25
7.63
AGl ! CPu ! SNC ! CL ! AGl AGl ! CPu ! MGP ! VL ! AGl
4 4
0 0
3 3
7.5 8.5
7.53 7.50
AGl ! CPu ! SNC ! VTA ! AGl
4
0
2.75
6.75
7.13
AGl ! CPu ! SNC ! AGm ! AGl
4
0
2.38
7
7.03
AGl ! CPu ! SNC ! PC ! AGl
4
0
2.75
7
7.03
AGl ! CPu ! SNC ! PF ! AGl
4
1
2.38
6.75
7.03
AGl ! CPu ! MGP ! PF ! AGl
4
1
2.75
7.75
6.97
AGl ! CPu ! SNC ! VM ! AGl
4
0
2.38
7.75
6.97
AGl ! CPu ! SNC ! MDL ! AGl
4
0
2
6.25
6.32
AGl ! CPu ! Ac ! VTA ! AGl
4
0
2.25
5.5
6.05
AGl ! CPu ! VL ! AGl
3
1
2.67
5.33
5.93
AGl ! CPu ! MGP ! VM ! AGl
4
0
2.5
7
5.85
AGl ! CPu ! STh ! PF ! AGl
4
1
2.25
4.75
5.78
AGl ! CPu ! AGl
2
1
2.5
5
5.55
AGl ! CPu ! MRF ! CL ! AGl
4
0
2
4.75
5.55
AGl ! CPu ! AGm ! VL ! AGl
4
1
2.5
4.75
5.32
AGl ! CPu ! VTA ! Ent ! AGl AGl ! CPu ! VM ! AGl
4 3
0 1
2.38 2.67
5.25 4.67
5.32 5.17
AGl ! CPu ! MRF ! PF ! AGl
4
1
2
4.25
5.07
AGl ! CPu ! MGP ! MDL ! AGl
4
0
2.5
6.25
4.88
AGl ! CPu ! MGP ! VTA ! AGl
4
0
2
5.75
4.88
AGl ! CPu ! AGm ! AGl
3
1
2.33
4.33
4.87
AGl ! CPu ! PC ! AGl
3
1
2.33
4.33
4.87
AGl ! CPu ! PF ! AGl
3
1
2.33
4.33
4.87
AGl ! CPu ! VTA ! AGl
3
0
2.33
5
4.83
AGl ! CPu ! MGP ! VA ! AGl
4
0
2
5.75
4.80
AGl ! CPu ! MRF ! VL ! AGl
4
0
1.25
4.25
4.78
AGl ! CPu ! VTA ! CL ! AGl
4
1
2.25
4.75
4.75
AGl ! CPu ! VTA ! AGm ! AGl
4
1
2.25
5
4.72
AGl ! CPu ! VTA ! VM ! AGl
4
1
2.5
4.75
4.72
AGl ! CPu ! AC ! VL ! AGl
4
0
2.25
4.25
4.70
AGl ! CPu ! STh ! MDL ! AGl AGl ! CPu ! STh ! VTA ! AGl
4 4
0 0
1.75 1.75
4 4.25
4.70 4.61
AGl ! CPu ! MRF ! PC ! AGl
4
0
2
3.75
4.57
AGl ! CPu ! MRF ! VM ! AGl
4
1
1.75
3.75
4.57
AGl ! CPu ! CM ! AGl
3
0
2.33
4
4.53
AGl ! CPu ! AGm ! CL ! AGl
4
1
2.5
4
4.53
AGl ! CPu ! AGm ! VM ! AGl
4
1
2.5
4
4.50
123
Brain Struct Funct Table 2 continued Path
PL
Rec
AW
Pub
Rel
AGl ! CPu ! AGm ! PC ! AGl
4
1
2.25
3.75
4.28
AGl ! CPu ! PF ! AGm ! AGl
4
1
2
3.75
4.28
AGl ! CPu ! CM ! AGm ! AGl AGl ! CPu ! PC ! AGm ! AGl
4 4
1 1
2 2.25
4 4
4.25 4.25
AGl ! CPu ! PF ! VM ! AGl
4
0
2.25
3.75
4.25
AGl ! CPu ! VTA ! CM ! AGl
4
0
2.25
4.75
4.15
AGl ! CPu ! MRF ! CM ! AGl
4
0
2
3.5
4.13
AGl ! CPu ! VTA ! PF ! AGl
4
1
2.5
4.75
4.13
AGl ! CPu ! MDL ! AGm ! AGl
4
0
1.75
4
4.05
AGl ! CPu ! AGm ! CM ! AGl
4
0
2.25
3.5
4.03
AGl ! CPu ! AGm ! PF ! AGl
4
1
2.25
3.5
4.03
AGl ! CPu ! VL ! AGm ! AGl
4
1
2.5
3.5
4.03
AGl ! CPu ! VM ! AGm ! AGl
4
1
2.5
3.75
4.03
AGl ! CPu ! PF ! VTA ! AGl
4
0
2
3.5
4.00
AGl ! CPu ! MDL ! AGl
3
0
1.33
3.67
3.93
AGl ! CPu ! MRF ! AGm ! AGl
4
0
1.5
3.5
3.88
AGl ! CPu ! MRF ! MDL ! AGl
4
0
1.75
3
3.88
AGl ! CPu ! VTA ! MDL ! AGl AGl ! CPu ! AGm ! MDL ! AGl
4 4
1 1
1.75 2.25
4 3.5
3.80 3.78
AGl ! CPu ! AGm ! VTA ! AGl
4
0
2
3.25
3.72
AGl ! CPu ! MRF ! VTA ! AGl
4
0
1.5
3.25
3.72
AGl ! CPu ! AC ! CM ! AGl
4
0
2
3.25
3.71
AGl ! CPu ! AC ! VTA ! AGl
4
0
1.5
3.75
3.61
AGl ! CPu ! VM ! VTA ! AGl
4
0
2.13
3
3.53
AGl ! CPu ! PC ! VTA ! AGl
4
0
1.75
3
3.50
AGl ! CPu ! MDL ! VTA ! AGl
4
0
1.25
3
3.33
PL path length, Rec reciprocal, AW average weight, Pub number of articles, Rel reliability
importance of a region in this network is estimated by the Shapley-value (Keinan et al. 2004; Ko¨tter et al. 2007). The lowest Shapley-value (largest importance) was calculated for CPu (0:309), then the AGl (0:277) and AGm (0:243) also have smaller values than other regions. The betweenness centrality (BC) (Freeman 1977) is a coefficient that divides the number of shortest paths between regions passing a particular region by the number of all shortest paths. It provides a useful measure for the interconnectedness and embedding of a particular region in the network (Vlachos et al. 2012). Again the CPu is most important and has the largest BC of 0.092, followed by the AC (0.074) and the AGm (0.069). The BC has been visualized in 3D in the bilateral BG network by assigning the BC values with a color code to each region (Fig. 6). The same ranking of CPu (1), AC (0.997) and AGm (0.817) holds for the authoritativeness as well (Kleinberg 1999). The local parameters have been visualized in parallel coordinates (Appendix Fig. 19) after sorting by their similarities. The presentation of the high-dimensional parameter space clearly indicates the central role of the CPu in
123
the BG2 network. Also the cortical input areas AGl and AGm are lying close to the CPu-curve in this visualization. A further finding concerns the laterality (LA ) of connections (Appendix Tables 8, 9). The AGm and AGl have a laterality of 0.62 and 0.607 (ipsilateral connections divided by DGall ). A smaller laterality value was found for MRF (0.552). Interestingly, the core regions of the BG all have lateralities larger than 0.7. This phenomenon is also found when considering reciprocal laterality (LR ). Rich-club analysis and knotty-centers A rich club is a set of high-degree nodes (rich nodes) which are connected with each other even more frequently than it would be expected by their high degree. Nodes with numerous links are much more likely to form well interconnected subgraphs (clubs) than low degree nodes. The rich-club phenomenon can be computed by the rich-club coefficient / (Zhou and Mondragn 2004). /ðkÞ measures the fraction of edges actually connecting nodes with degree [ k out of the maximum number of possible connections
Brain Struct Funct
(Colizza et al. 2006). Rich-club analysis has been applied successfully to detect structural linkages of hubs of the cerebral cortex (Shanahan et al. 2013; Harriger et al. 2012; van den Heuvel and Sporns 2011). The diagram of richclub coefficients in dependence on node degrees /ðkÞ of the bilateral BG2 network allows to detect a rich-club
phenomenon of nodes with degrees k 44 (Fig. 7). These nodes have a / 0:8. The feeder nodes are those nodes of the BG2-network that have direct connections to the richclub nodes. This connectivity is shown in Appendix Fig. 20. It visualizes the dense connectivity of rich nodes and a substantial connectivity between rich-club nodes and feeder nodes. CPu, AGl, AGm and SNR are members of the rich club. Interestingly, the non-canonical BG regions AC, VTA and MRF are also rich-club nodes. The knottycentrality (KC) quantifies the connectivity within a topologically central connective core (Shanahan and Wildie 2012). The knotty-center consists of the following regions found in the ipsi- and contralateral hemispheres (Fig. 21): VTA, MRF, CPu, AGl, AGm, AC. The thalamic output regions do not belong to the knotty-set. However, knottyset members MRF and VTA are massively connected with cortical as well as subcortical hub regions. Motif analysis
Fig. 5 The weighted-modularity analysis of the ipsilateral BG2 network. Five modules have more weighted connections within than among the modules. Thick lines indicate large weights
Motif analysis has been applied by cortical and subcortical connectomes in different species to determine such subgraphs that are significantly more abundant or considerably less frequent than in suitable randomizations (Sporns and Ko¨tter 2004; Sporns 2006; Schreiber and Schwo¨bbermeyer 2004, 2005). The motifs of the ipsi- (Fig. 8a) and bilateral (Fig. 8b) BG2 networks show characteristic distribution patterns. The frequency of 13 3-node motifs and 3 complex 5-node motifs was determined. The 5-node motifs consist of 2 reciprocal edges and two non-reciprocal edges which
Fig. 6 The weighted connections and color-coded weighted betweenness (BC) parameter of the bilateral BG2network. The color code of connections weights is the same as in Fig. 3d. The weighted betweenness centrality BC of the regions has been integrated into this 3D-view. The CPu belongs to the class of maximalweighted BC values indicating its connectional importance in the BG2-network AGm and Pir have smaller values. The large amount of contralateral projections is also visible
123
Brain Struct Funct Fig. 7 The rich-club coefficients / of the ipsi and contralateral BG2 network in dependence of the ranks k of the regions. 1,000 rewiring randomizations were applied and mean rich-club coefficients were computed. A secondary yaxis in green shows the coefficient /ðkÞ=rand/ðkÞ
are oriented in 3 ways: chain-like (chain), centripetal (in) to and centrifugal (out) from the center vertex to the vertices with reciprocal edges. The frequency of motifs in the real ipsilateral BG2 network can be compared with the frequencies of motifs in 1,000 rewiring randomizations with the same number of nodes and edges. It turns out that the motifs 3-09 and 3-13 are significantly more abundant in the real BG2 network than in randomizations (Fig. 8a). In the bilateral BG2 network 3-04, 3-06, 3-09, 3-13 motifs are more abundant than in 1,000 rewiring simulations. All these 3-node motifs have at least one reciprocal edge. In contrast to the ipsilateral network, the 5-node motifs are very abundant (Fig. 8b). The contribution of regions to motifs is shown in Table 3 and Appendix Table 10. AC, CPu and VTA are most often elements of the fully reciprocal motif 3-13 (Table 3). The 5-node motifs only appear sporadically. The MGP is most often (326) an element of the centripetal 5-node motif and 333 times an element in the chain 5-node motif. The centrifugal 5-node motif is built by AGl in the majority of cases. In the bilateral BG2 network MRF participates most often in the full reciprocal motif 3-13. Multivariate analysis Modularity analysis provides a classification of regions based on their number of connections in modules and in between modules. The result of a weighted modularity analysis has been described before (Fig. 5). Metric multidimensional scaling (MDS) arranges regions that have similar connectivity patterns closer together than to regions that have more different connections. The result of a MDS in a normalized and quadratic diagram, respectively, is shown in Appendix Fig. 24 for the ipsilateral and in Appendix Fig. 25 for the bilateral BG2 network. Most of the thalamic regions are located within the lower left quadrant whereas LHb and MDM are located in the upper left quadrant. Cortical motor regions, MGP, CPu, SNR and
123
MRF are in the lower right quadrant. SNC and VTA are located close together in the upper right quadrant. Ent, Pir and HIPP which also share functional features are relatively closely arranged in the upper left quadrant. In the bilateral BG2-network cortical motor, SNC, SNR, VTA, CPu and AC are located close together in the lower quadrants whereas MRF appears as an outlier in the lower quadrants. All other regions are located in the upper quadrants. The principal component analysis was applied to the DGAll , CluCAll , CluC2 , AvgDGnb , VCDG , Loc as proposed by (Echtermeyer et al. 2011) (Appendix Table 13, the definitions of these parameters are summarized in the Appendix) to determine those two dimensions that have the largest influence in the bilateral BG2 network (Fig. 9a, b). It turns out that 5 parameters of the first component (first row of Appendix Table 13) have similar absolute values and determine the distribution of regions in the first dimension. The second dimension (second row of Appendix Table 13) is determined by the largest absolute value by CluC2 . This means, that the variance within the y-axis of the PCA-plane is dominated by the variance of CluC2 . The gray values of the PCA-plane diagram (Fig. 9a) indicate the densities of regions within a ring. The lightest gray value in the middle covers most regions. The STh was selected and the connectional relations of STh in the BG2 network are shown in the circular layout in Fig. 9b. In the inner ring 24 directly connected neighbors are located and in the outer ring 25 indirect neighbors are shown. The number of indirect and direct neighbors is well balanced and typical for all regions within the center of the PCA plane. This result is important since the target of deep brain stimulation is the STh. The CPu shows another pattern of direct–indirect-neighbor connectivity. The number of direct neighbors is much larger than the number of indirect neighbors (Fig. 9c). Interestingly, the number of indirect neighbors of VL (smaller first component than within the CPu) is larger than the number of direct neighbors (Fig. 9d).
Brain Struct Funct
a
b Fig. 8 Overview of the distributions of motifs in the ipsi and contralateral BG2 network. a The frequencies of motifs in the ipsilateral BG2 network. 13 directed 3-node motifs and 3 directed 5-node motifs have been searched: syOut (‘‘Out’’ with regard to the central node of the motif), syIn and syCh. The frequency of a motif in the real network is indicated by the blue point. The BG2 network was randomized using a rewiring procedure. Each randomization contains exactly the same number of regions and connections. The rewiring randomization was repeated 1,000 times, after each randomization the frequencies of motifs were determined, indicated by a small black dot. The red bar indicates the variation of frequency within the standard
syOut
syIn
syCh
syOut
syCh
syIn
deviation. The motifs are sorted by their z values. Hence, motifs that have larger frequencies in the real network appear on the right side of the diagram. Most motifs that have reciprocal edges are more abundant in the real network than in the 1,000 rewiring simulations. b The frequency of motifs in the contralateral BG2 network. The complex motifs syOut, syCh and syIn are much more often represented in the contralateral BG2 network than in the ipsilateral case. The differences of frequencies of motifs in the real BG2 network and in rewiring BG2 networks are stronger in the contralateral than in the ipsilateral case
123
Brain Struct Funct Table 3 Frequency of motifs of the bilateral BG2 network Region
3-01
3-02
3-03
3-04
3-05
MRF_L
121
MRF_R
109
27
2
228
26
28
3
220
23
AC_L
32
17
15
184
AC_R
32
16
15
VTA_L
41
28
11
AGm_L
22
36
VTA_R
41
31
AGm_R
23
CPu_R
3-06
3-07
3-08
3-09
3-10
3-11
3-12
3-13
chain
in
out
36
3
115
98
33
20
164
125
7
0
0
42
2
112
103
31
18
157
125
2
0
0
16
56
2
54
134
19
16
154
114
72
130
8
191
16
56
2
53
162
21
16
157
111
68
136
8
166
10
48
0
54
110
25
15
139
95
78
54
25
24
145
16
112
0
25
151
21
26
105
93
188
273
82
13
162
10
51
0
51
104
24
15
131
90
113
60
25
42
28
143
18
120
1
25
140
22
27
102
88
212
254
77
30
29
30
103
33
176
2
44
94
35
41
127
80
135
42
45
CPu_L
31
28
33
103
33
172
2
41
93
37
39
131
80
114
51
58
AGl_L AGl_R
8 8
7 8
5 5
76 73
1 1
135 136
1 1
2 2
198 184
15 15
4 4
80 72
78 76
249 294
104 108
131 173
PF_L
25
24
15
150
9
65
0
31
100
8
7
67
62
145
160
33
PF_R
22
20
15
145
8
64
0
29
102
7
5
63
62
136
184
30
STh_L
20
45
14
144
6
60
3
23
81
18
13
60
60
265
161
61
STh_R
17
42
12
141
6
63
2
21
83
16
11
55
60
243
173
69
SNR_L
77
55
39
142
13
91
5
63
54
33
13
58
56
193
191
63
SNR_R
78
56
39
141
13
90
5
63
53
33
13
58
56
190
204
51
SNC_R
58
31
20
140
29
77
1
55
41
23
25
68
55
168
153
96
SNC_L
58
30
20
141
29
78
1
55
42
23
25
68
55
151
135
107
MGP_R
34
50
23
105
12
67
2
34
62
17
6
52
49
333
326
105
MGP_L
34
48
23
108
12
67
2
34
63
17
6
52
49
329
299
100
Ac_L
32
45
20
140
19
51
2
17
44
28
16
68
33
201
165
75
LHb_R
33
57
14
124
20
67
6
30
48
23
7
60
31
251
155
111
LHb_L
33
57
14
127
20
66
6
30
48
23
7
60
31
242
149
93
Ac_R
35
52
24
137
21
53
2
20
40
28
15
63
30
192
163
87
VM_R VM_L
51 51
30 29
13 13
99 101
26 26
72 73
2 2
31 31
83 83
11 11
6 6
56 56
27 27
149 141
68 68
34 34
CM_R
70
39
17
110
24
66
4
27
88
16
10
49
21
108
54
20
CM_L
67
37
19
104
24
64
4
28
85
21
12
51
21
101
55
19
Ent_L
25
43
29
76
10
79
3
16
67
14
12
51
19
281
158
112
PC_R
52
20
12
77
2
54
0
6
98
7
5
29
14
120
74
51
PC_L
52
20
12
77
2
53
0
6
99
7
5
29
14
103
73
56
HIPP_L
37
58
33
82
24
117
2
23
35
21
22
62
13
258
120
135
HIPP_R
37
53
32
91
27
112
4
21
39
23
22
59
10
277
189
143
Pir_L
56
47
18
88
29
71
3
44
37
22
16
46
10
249
110
100
CL_R
81
43
18
75
8
53
3
19
59
21
9
30
8
96
51
63
CL_L
82
44
18
75
8
52
3
19
61
21
9
30
8
91
48
67
MDL_R
76
37
8
103
18
24
0
29
62
4
2
35
7
84
145
39
MDL_L
76
38
8
103
18
23
0
29
63
4
2
35
7
71
136
39
Pir_R
59
55
21
85
32
71
4
42
37
24
13
44
6
250
91
96
LGP_L LGP_R
46 44
30 31
7 7
89 87
6 6
24 24
2 2
12 12
19 19
8 8
3 3
25 25
6 6
194 183
193 204
29 26
VL_R
105
28
23
113
22
49
0
46
45
7
8
26
6
62
191
30
VL_L
105
29
23
113
22
48
0
46
46
7
8
26
6
59
181
35
MDM_L
77
19
2
108
6
18
0
20
31
1
8
19
5
159
140
48
MDM_R
78
20
2
105
6
12
0
20
32
1
7
15
2
91
116
17
123
Brain Struct Funct Table 3 continued Region
3-01
3-02
3-03
3-04
3-05
3-06
3-07
3-08
3-09
3-10
3-11
3-12
3-13
chain
in
out
Ent_R
4
19
4
28
0
13
0
0
28
1
0
2
1
57
63
34
VA_L
30
22
17
58
1
36
1
2
0
6
1
0
0
0
67
0
VA_R
30
22
17
58
1
35
1
2
0
6
1
0
0
0
70
0
Regions are sorted using motif 3-13 frequencies chain symmetricCHAIN, in symmetricIN, out symmetricOUT (see legend of Fig. 8)
a
b
c
d
Fig. 9 Principal component analysis of the ipsi and contralateral BG network. a PCA plane of the ipsi- and contralateral BG2 network. The regions of the BG2 network are distributed in the PCA-plane (defined by the first (abscissa) and second component (ordinate) of the PCA) as strongly overlapping pairs. The gray values within a ring present the density of regions in the ring. Most regions are located in the innermost ring. This presentation of regions in the PCA plane need to be considered together with the values of the local network parameters (Appendix Table 13) and their contribution to the first (abscissa) and second component (ordinate). b STh was selected from the bilateral BG2 network. The direct (first) neighbors of STh are located on a ring around STh. The connections of these regions are
dark blue. The indirect (second) neighbors of STh are on the external ring. Connections between second neighbors are light blue and connections between second and first neighbors are yellow. The number of direct and indirect neighbors is similar as well as the number of connections inside the sets of neighbors. c CPu was selected from the bilateral BG2 network. The number of first neighbors of the CPu is much larger than the number of second neighbors. Also the density of connections between second neighbors is much more different from the knotty-center presentation of STh (Appendix Fig. 21). d VL was selected from the bilateral BG2 network. The number of indirect neighbors is much larger than the number of direct neighbors of the VL
123
Brain Struct Funct Table 4 Vulnerability of regions in the ipsilateral BG2 network Region Caudate putamen
Vulnerability 2.2792
Amygdaloid complex
2.1273
Ventral tegmental area A10
1.7475
Mesencephalic reticular formation
1.5196
Medial agranular prefrontal cortex
1.4057
Substantia nigra compact part
0.8360
Lateral agranular prefrontal cortex
0.7600
Substantia nigra reticular part
0.7220
Subthalamic nucleus
0.6461
Medial globus pallidus
0.6081
Accumbens nucleus Parafascicular thalamic nucleus
0.3802 0.2663
Lateral habenular nucleus
-0.3034
Piriform cortex
-0.3414
Entorhinal cortex
-0.4173
Ventromedial thalamic nucleus
-0.5313
Hippocampus
-0.6072
Central medial thalamic nucleus
-0.6452
Mediodorsal thalamic nucleus lateral part
-0.9870
Centrolateral thalamic nucleus
-1.1010
Lateral globus pallidus
-1.1010
Ventrolateral thalamic nucleus
-1.2529
Paracentral thalamic nucleus
-1.3289
Mediodorsal thalamic nucleus medial part
-1.3668
Ventro anterior thalamic nucleus
-2.5062
Differential connectomics: vulnerability Vulnerability analysis provides an estimate of the effect of removing a region from a network. If the CPu is removed, then the mean distance between all resting regions increases by 2.279 % and if AC is removed then the mean distance increases by 2.127 %. A sorted list of vulnerabilities of regions is shown in Table 4. The number of connections or the degree of regions is positively correlated with most local network measures especially the centrality parameters (Appendix Table 12). The number of connections of the ipsilateral BG2 network was visualized using the K-core analysis where regions are sorted with regard to their number of connections and those regions which have the smallest number of connections are removed until the regions with the largest number of connections remain in the center of the K-core visualization (Appendix Fig. 22). This has been performed for the ipsilateral BG2 and the bilateral BG2 network (Appendix Fig. 23). It turns out that CPu, AGl, AGm, VTA, MRF and PF are the core regions of the BG2 network. The MDS of the bilateral BG2 network is relatively symmetric. In this case, all thalamic regions of the ipsilateral hemisphere are
123
arranged in the upper left quadrant. Ent, Pir and HIPP are in between these thalamic regions. SNR, SNC, CPu, AGl and AGm are located in the lower left and lower right quadrant, respectively, for the contralateral regions. If MDS is performed on the weighted connectivity matrix, then the regions are moving closer together in direction of the midline of the MDS diagram (not shown). Connectomics of an embedded BG-network An embedding of the BG1 network into a network of regions which have direct input and output connections to the BG1 regions was performed to determine the importance and ranks of BG1 regions in the context of extrinsic direct and indirect neighbors. The results of the analysis are shown in Appendix Table 14. Changes of ranks of regions were found when comparing the BG1 network with the network of all regions that are directly connected to the BG1 network. In terms of BC, the CPu, AGl and STh are covering the first three ranks. In the direct neighbor network (DNN) the AGl, AGm and LGP have ranks 1, 2 and 3. In the indirect neighbor network (INN) CPu, SNC and AGl have ranks 1, 2 and 3. With regard to EC, the SNC, CPu and MGP are on the first three ranks in the BG1 network, in the DNN and INN these are SNC, SNR and MGP. Hence, with regard to EC there is not so much change in the ranks as compared with the BC parameter. Interestingly, the ranking of regions by means of PRC, SC, Hub, Aut, RADout and RADin in the BG1 network, DNN and INN is similar: CPu, AGl, STh, MGP, SNC or SNR are covering the first three ranks. A further interesting finding is that VL is considered as an important output region of the BG, however, the region rankings of the embedded BG1-network indicate variable ranks of 4 or 7, only. Furthermore, the CPu changes its rank in the DNN when considering the RADin parameter from 1 in the BG1-network to 7. Global extrinsic connectivity of the BG2 network The global extrinsic connectivity of the BG2 network has been determined ipsilaterally (Table 5) and bilaterally (Appendix Table 11). Extrinsic connections are those where at least the source or target is not a member of the set of regions of BG2. The CL has most direct ipsi- and contralateral inputs (284) in the BG2 network. The AGl has most direct outputs followed by MRF and AGm. However, AGm has the most ipsilateral outputs (516). When considering all connections of subregions of a region of BG2, then the MRF has most inputs (78,691) followed by AC (5,424) and CPu (1,536). The ranks of regions changes when regarding their subtree outputs: AC is 4,903, for MRF it is 1,893 and for AGl it is 1,375. In the bilateral
Brain Struct Funct Table 5 Global extrinsic connections of regions of the ipsilateral BG2 network Region
Dic
Dii
Dis
Doc
Doi
Dos
Sic
Sii
Sis
Soc
Soi
Sos
Centrolateral_thalamic_nucleus_L
398
446
844
10
68
78
413
470
883
10
103
Accumbens_nucleus_L
277
378
655
5
74
79
373
829
1,202
17
304
321
Central_medial_thalamic_nucleus_L
278
322
600
18
70
88
288
340
628
18
184
202
Paracentral_thalamic_nucleus_L
242
309
551
10
81
91
245
413
658
10
89
99
Mesencephalic_reticular_formation_L
226
311
537
235
334
569
3,418
4,273
7,691
1,383
3,510
4,893
Amygdaloid_complex_L
148
375
523
15
71
86
951
4,473
5,424
933
3,970
4,903
Ventrolateral_thalamic_nucleus_L
154
344
498
12
52
64
175
384
559
13
56
69
82
311
393
31
79
110
479
1,057
1,536
35
172
207
Parafascicular_thalamic_nucleus_L
124
213
337
42
136
178
160
377
537
127
291
418
Ventromedial_thalamic_nucleus_L
91
244
335
37
54
91
128
356
484
44
161
205
Ventral_tegmental_area_A10_L Medial_agranular_prefrontal_cortex_L
78 70
238 131
316 201
138 130
231 416
369 546
127 112
476 212
603 324
188 203
440 584
628 787
Hippocampus_L
37
137
174
27
46
73
263
1,030
1,293
228
631
859
Entorhinal_cortex_L
35
137
172
36
63
99
166
900
1,066
80
546
626
Lateral_habenular_nucleus_L
30
128
158
44
99
143
70
226
296
52
129
181
Subthalamic_nucleus_L
25
124
149
31
83
114
38
181
219
42
127
169
Substantia_nigra_compact_part_L
32
101
133
48
126
174
47
136
183
58
181
239
Lateral_agranular_prefrontal_cortex_L
30
95
125
284
381
665
89
436
525
574
801
1,375
Piriform_cortex_L
17
103
120
27
53
80
86
300
386
73
191
264
Mediodorsal_thalamic_nucleus_medial_part_L
15
91
106
1
14
15
15
91
106
1
14
15
Mediodorsal_thalamic_nucleus_lateral_part_L
22
82
104
2
18
20
22
83
105
2
18
20
Substantia_nigra_reticular_part_L
20
73
93
50
133
183
20
86
106
74
200
274 91
Caudate_putamen_L
113
Medial_globus_pallidus_L
15
53
68
17
59
76
16
58
74
18
73
Lateral_globus_pallidus_L
10
22
32
0
4
4
10
22
32
0
5
5
Ventro_anterior_thalamic_nucleus_L
12
12
24
0
10
10
15
14
29
0
10
10
The regions were sorted by the direct input sum (Dis) Dic direct input from contralateral, Dii direct input from ipsilateral, Dis direct input from ipsi- and contralateral, Doc direct output to contralateral, Doi direct output to ipsilateral, Dos direct output to ipsi- and contralateral, Sic subtree input from contralateral, Sii subtree input from ipsilateral, Sis subtree input from ipsi- and contralateral, Soc subtree output to contralateral, Soi subtree output to ipsilateral, Sos subtree output to ipsi- and contralateral
BG2 network (Appendix Table 11) CL has most inputs (838) followed by Ac (645) and CM (587). The AGl has most outputs (640) followed by MRF (561) and AGm (524). Most subtree (for details see ‘‘Materials and methods’’) inputs are going to MRF (7,589), AC (5,291) and HIPP (1,358). The output from subtrees is dominated by AC (4,821) followed by MRF (4,700) and AGl (1,307). Local extrinsic connectivity of the BG2 network In this section extrinsic regions are determined that have abundant connections to the intrinsic regions of the BG2 network. Furthermore, extrinsic regions were determined that receive many connections from intrinsic regions. To identify these extrinsic regions which are densely interconnected with intrinsic regions, a search through the whole connectome was performed to filter out those regions which are located around hierarchical level 9 and
possess direct ipsilateral inputs or outputs to the intrinsic regions. Those regions that have no input or no output were removed from this selection. After this conditional expansion of the network, the interactive hierarchical table method of neuroVIISAS was used to determine extrinsic regions that have abundant inputs to or get outputs from intrinsic regions. It turns out that the medulla oblongata provides most inputs to the CPu, followed by the lateral thalamic nuclear group, medial geniculate nucleus, ventral thalamus and intralaminar nuclei of the thalamus. Most outputs of the CPu are going to intrinsic nuclei. However, the lateral thalamic group and the substantia innominata nucleus basalis complex receive a few outputs of the CPu. Abundant inputs to the SNC derive from the medulla oblongata, tuberomammillary nucleus and the pedunculopontine tegmental nucleus. Most external outputs from SNC are going to the ventral thalamus. The medial global pallidus does not have many extrinsic outputs nor does it
123
Brain Struct Funct
receive many extrinsic inputs. The subthalamic nucleus receives abundant extrinsic inputs from the pedunculopontine tegmental nucleus and the ventral thalamus. The sole outputs of MGP are going to the olfactory tubercle. The accumbens nucleus receives many inputs from mesocortical regions, the medulla oblongata as well as from hypothalamic zones. Most extrinsic outputs are going to hypothalamic zones, the olfactory tubercle and the medulla oblongata. The VL receives abundant connections from the cerebellum and hypothalamic zone. VL projects extrinsically to the medulla oblongata and mesocortical regions.
Discussion The two levels (BG1, BG2) of the mesoscale connectome of the BG were investigated with regard to intrinsic, extrinsic, weighted, ipsi- and contralateral configurations. Connectome embedding and differential connectomics were applied to determine the importance of regions (ranking) in terms of their connectivity in the network. The routing problem and reciprocity patterns The three principal pathways (direct, indirect, hyperdirect) were revealed in the BG connectome. By applying pathway analysis, many further possible pathways from the AGl over the CPu back to AGl were detected in the ipsilateral BG2 network. So far, alternative pathways in the orofacial system (Mascaro et al. 2005) and the pedunculopontine nucleus of the interface between the BG and brainstem nuclei (Martinez-Gonzalez et al. 2013) have been described. Which pathway or set of pathways are used predominantly in terms of biosignal processing and function cannot be worked out by network analytical methods. Maybe this routing problem can be resolved by considering the coding features of spike patterns (intrinsically bursting, regularspiking, fast-spiking) or inhibitory and excitatory functions of local circuits (Tiesinga et al. 2008; Steriade 2004). In an ongoing meta-study, the transsynaptic pathways after virus (especially pseudorabies) injections are being collated and will be compared with the pathway results of monosynaptic tracing studies. Certain reciprocal and non-reciprocal patterns between MGP and SNR to thalamic and further on to the AGl were detected. Continuous reciprocal links of VM and PF between BG output nuclei and AGl were determined. However, VTA is reciprocally connected with SNR and MGP but non-reciprocally with AGl. Rat connectome studies There are not so many publications available which report graph theoretic analyzes of partial connectomes of the rat
123
nervous system. A major connectional analysis of some parts of the rat central nervous system was made (Burns 1997) and a sophisticated web-page (BAMS2 ) for connectivity retrieval of the rat CNS is available (http://bran cusi.usc.edu). Connectome studies focusing on specific parts or functional systems of the rat CNS were performed in terms of the nucleus of the solitary tract (Palombi et al. 2006), the hippocampus (Burns and Young 2000), the thalamocortical connectivity (Costa and Sporns 2006), the brainstem reticular formation (Humphries et al. 2006) and the retrosplenial cortex (Sugar et al. 2011). Moreover, there exist an increasing amount of meta-studies of connectomes of many other organisms from Caenorhabditis elegans (Towlson et al. 2013) up to primates (Ko¨tter 2004; Stephan et al. 2001a; Passingham et al. 2002). The unique feature of the rat project within neuroVIISAS is (1) the preservation of hierarchical terminologies in the form of an advanced neuroontology (Bota and Swanson 2007a), (2) the explicit collation of ipsi- and contralateral connections, (3) the extraction of collateral projections, (4) transneuronal and transsynaptic connectivity, (5) the multimodal integration of digital atlases, (6) and the integration of populationbased simulations using NEST (Gewaltig and Diesmann 2007). Moreover, the rat connectome project encompasses about 90 % of the available tract-tracing literature of the rat nervous system indexed in PubMed (http://www.ncbi. nlm.nih.gov/pubmed). A list of all reports is available on the webpage http://neuroviisas.med.uni-rostock.de/refer ences.html. In the future, the databases Web of Science, Scopus and Google Scholar will be queried for searching those reports that are not indexed in PubMed. Nevertheless, it is unknown how many connections at the scale of nuclei exist (ground truth of connections) and how dynamic changes of connectivity strengths influences variability of the estimates of connectivity strengths. In particular, the question of completeness becomes evident with regard to contralateral connectivity because in most tract-tracing studies contralateral connections are not systematically investigated. This means that especially contralateral connectivity should be considered as underestimated. In the macaque Barbas et al. (2005) provide evidence that ipsiand contralateral connectivity depend on topological proximity and the structural type of linked areas. Furthermore, they found that projections of both hemispheres are highly correlated with the structural architecture of cortical areas and that the contralateral connectivity mirrors ipsilateral connectivity while being weaker. This suggests that contralateral connectivity brings out a specific laterality pattern which should be investigated further. Hence, the successive accumulation of contralateral connections in the rat connectome project should elucidate the many aspects of contralateral connectivity at the cortical and subcortical level of the central nervous system.
Brain Struct Funct
Reliability of connections
Network granularity
It should be emphasized that some connections and regions may receive more attention in tract-tracing research so that frequently investigated or popular regions may have higher observation scores. This means, that the observation score introduced here should be taken with caution. There could be different reasons like disease (Parkinson syndrome, Huntington-disease and so forth), function (whisker cortex) or accessibility (size of source and target regions in tract tracing experiments, e.g., CPu) why connections or regions are investigated more often than others. Furthermore, most experimental tract-tracing studies typically have been hypothesis driven and restricted to the analysis of one or a few brain regions. The number of tract tracing studies that aim to map connections of a larger functional or structural system of the CNS is rather limited. However, the mouse connectome project at the Allen Institute (http://www. mouseconnectome.org) performs a high-throughput mapping of the interconnectivity of the whole mouse brain based on the tract-tracing technique (Oh et al. 2014). Nevertheless it is useful to have the reliability score available, so that one can filter the connectivity matrix to only contain results that have been observed multiple times. It is reasonable to consider specific reliability weights of connections between two regions in a particular report and in independent reports. In such a way, it could be possible to judge inherent observer biases of connections. The analysis of observation scores shows large values for most connections of the BG network data. Especially, the ipsilateral connections have larger reliability values than the contralateral connections. Obviously, the ipsilateral connectivity of the bilateral BG1 network is stronger than the contralateral connectivity. Furthermore, the laterality of cortical regions is larger than those of the subcortical regions. In the chord diagram of the BG1 network the strong connectivity of SNC to CPu attracts attention. However, in the vulnerability analysis the removal of the SNC, comparable with its neuro-degenerative loss of dopaminergic neurons in the SNC in Parkinson disease, shows an unexpected 6th position in the vulnerability ranking, only. This is in agreement with the modularity role of the SNC, indicating that the SNC is not a complex integrator or hub-region (hub = 0.672) (Sporns et al. 2007). This is also confirmed by the relatively large Shapley-value of 0.008 in the bilateral BG2 network (Small or negative Shapley-values indicate an increase of importance of nodes in a network. Because the Shapley-value is a sum of the differences of connected components, negative values occur. The formal definition can be found in the Appendix.).
The BG2 network contains functionally important limbic regions and specific thalamic, mesencephalic and brainstem regions. In particular, the brainstem and mesencephalic regions show a very dense connectivity to most BG2 regions. However, diencephalic regions of the thalamic nuclei have dense connections to cortical areas and specific subcortical regions of the BG. The larger granularity of the BG2 network emphasizes the impressive contralateral connectivity as shown in the planar network diagram (Fig. 4). This means that the laterality should be considered in dependence of the network granularity. The global connectivity of the BG2 network indicates that it has a small world (Sporns 2006; Bassett and Bullmore 2006) and scale-free structure (Stam and Reijneveld 2007). However, it remains to be seen whether the conclusion will hold when the complete rat connectome is considered (LimaMendez and van Helden 2009). Connectome visualization and thalamic regions The BG2 network was also visualized in 3D-stereotaxic coordinate space (Paxinos and Watson 2007) in neuroVIISAS (Fig. 1) in combination with local-network parameters and weighted connections (Fig. 2). Before computing the surfaces of regions, the contours of all regions of the stereotaxic atlas were manually edited and related to the ontology (Schmitt and Eipert 2012). The surfaces were determined by the marching-cube algorithm of VTK (http://www.vtk.org/) in neuroVIISAS. In the BG2network the CPu has a strong ipsilateral connectivity whereas cortical regions can be distinguished by many contralateral connections. The thalamic PF region is a particular region which has many direct (adjacency matrix) and indirect (communicability matrix) connections to most other regions of the BG1- and BG2-network. In contrast, VA appears to be a region with marginal connectivity only. Modularity analysis and MDS The modular grouping that was found by modularity analysis and MDS indicates a relation of the connectional architecture with functional entities. The regions of the ipsilateral and contralateral network were also grouped by applying metric multidimensional scaling (MDS). Those regions that have similar connections are located close together in the MDS-diagram (Appendix Fig. 24) and those which have different connections are located far away from each other. In the ipsilateral BG2 network nearly all thalamic regions are located in one quadrant of the MDS-
123
Brain Struct Funct
diagram and the motor cortical regions, CPu, MGP and SNR are located together with MRF in another quadrant. A trend towards such a thalamic—basal ganglia—grouping was also found in the bilateral BG2 network (Appendix Fig. 25). Local parameters of the basal ganglia connectome Most local parameters of the ipsilateral BG2 network show relatively high correlations with the DGAll parameter. Especially the centrality parameters build a group with similar numeric trends in the parallel coordinate visualization (Appendix Fig. 19). Comparable patterns of correlations of local network parameters were also described by Bounova and de Weck (2012). However, not for each region does the correlation of DGAll (Appendix Table 12) and another parameter like the Shapley-value and EC-value show comparable trends. The CPu has rank 2 with regard to DGAll , however, for EC it has rank 3 (0.92) and for the Shapley-value it has rank 1 (0:31) (Table 6). AGl has a Shapley-value of 0:28 and rank 2, whereas, the EC-value is 0.685 and has a rank of 11, only. The DGAll of AGl is 31 which equals rank 8. AGm shows comparable values: its Table 6 Local parameters of the ipsilateral BG2 network
The regions are sorted using the degree all. Abbreviations and parameter definitions are described in the Appendix. A complete list of all local parameters is shown in the Appendix Table 15 A all, DG degree, I in, O out
123
Shapley-value is 0:24 (rank 3) and the EC-value is 0.693 (rank 10), however, its DGAll is 34 with rank 5, only. This means that the total number of connections of a particular region does not fully determine its importance. Furthermore, connectional aspects like connections to other important regions (authoritativeness) and local clustering influence its importance on the connectional level. The weighted-ipsilateral BG2 network provides other ranks, yet the CPu still belongs to the group of regions with highest ranks. The VTA and SNR have ranks 2 and 3 with regard to DGAll . The CPu and VTA stay on ranks 1 and 2 when ECvalues are considered. However, SNR has only rank 7 and MRF gets up to rank 3. The Shapley-value cannot be computed for weighted networks. Hence, some numeric similarities can be observed when connectional roles of particular regions are considered in a non-weighted and weighted BG2 network. A complex region at midbrain and brainstem levels like the MRF with many large weights may strongly change the ranks of importance of regions in a particular functional system like the BG. This means important regions of an isolated functional system or subnetwork like the BG2 that is embedded in a significantly more complex network may change their significances of
Region
REC
DGA
DGO
DGI
BC
Shapley
Hub
Aut
AC
20
42
22
20
0.0736
-0.16
0.96
1
CPu
18
41
20
21
0.0923
-0.31
0.89
1
VTA
17
40
23
17
0.051
-0.1
1
0.89
MRF
15
38
23
15
0.0287
-0.01
1
0.87
AGm
15
34
17
17
0.0694
-0.24
0.72
0.82
SNC
11
32
21
11
0.0177
0.03
0.93
0.67
SNR
11
32
19
13
0.0279
-0.04
0.81
0.77
AGl
14
31
15
16
0.0487
-0.28
0.67
0.76
MGP
11
30
17
13
0.0239
-0.02
0.74
0.77
STh
12
30
17
13
0.0175
0.02
0.78
0.78
PF
11
28
17
11
0.0073
0.06
0.78
0.73
Ac
10
28
13
15
0.0214
0.01
0.58
0.82
LHb
8
22
11
11
0.0057
0.07
0.51
0.71
Pir
5
22
8
14
0.0105
0.07
0.36
0.8
VM Ent
7 8
21 21
9 9
12 12
0.0034 0.0086
0.11 0.05
0.45 0.42
0.77 0.65
CM
5
20
10
10
0.0042
0.12
0.48
0.66
HIPP
6
20
11
9
0.0048
0.11
0.52
0.55
MDL
4
17
4
13
0.0025
0.19
0.2
0.8
CL
4
16
8
8
0.0017
0.14
0.38
0.55
LGP
4
16
5
11
0.0017
0.16
0.23
0.69
PC
5
14
6
8
0.0009
0.15
0.31
0.56
MDM
4
14
4
10
0.0013
0.2
0.19
0.59
VL
3
14
3
11
0.0008
0.25
0.15
0.7
VA
0
5
2
3
0.0001
0.42
0.1
0.18
Brain Struct Funct
regions heavily. This was proved by embedding the lesscomplex BG1 network (contains all BG-core regions) into (1) the extrinsic direct neighborhood network and (2) into the indirect neighborhood network. It was found that some regions-ranks of the direct and indirect networks are constant in terms of EC whereby SNC has rank 1 in this comparison. For the parameters PRC, SC, Hub, Aut, RADout and RADin region-ranks change a little, however, the most important regions stay in the upper ranks. The low ranking of VL does not reflect its significance with regard to its BG output role. This finding is in line with the high variability of regions which have direct input to the AGl in the pathway analysis and supports the hypothesis that the projections from BG-core regions are relatively divergent to different thalamic nuclei. In the non-weighted bilateral BG2 network the MRF, AGm, AGl and VL have the largest DGAll values. However, the Shapley-value of the CPu is conspicuously small (0:37) as well as the value for the AGl (0:37), followed by the AC (0:35). Again, another context like the bilateral selection of BG2 regions changes the importance and ranks of particular regions. With this direct confirmation of the change of importance of brain regions in different contexts of networks it is emphasized that connectomeanalysis depends strongly on the selection of intrinsic regions and extrinsic regions of particular functional subsystems of nervous systems. Motif-analysis A further local feature is the frequency by which regions contribute to particular motifs (Fig. 8a, b; Table 3 and Appendix Table 10). This was investigated in the ipsilateral and contralateral BG2 network. The most impressing result of this analysis was the significantly higher frequency of some motifs (3-12, 3-13) with reciprocal edges. A significantly higher number of complex 5-node motifs was found in the BG2 networks. Those regions, that have high degrees are more frequently found to build reciprocal motifs such as AC, CPu, VTA and MRF in the ipsilateral BG2 network and MRF, AC, VTA and AGm in the bilateral BG2 network. The frequencies in Table 3 of the bilateral BG2-network were sorted with regard to the fully reciprocal motif 3-13. However, regions with high abundance in motifs 3-12 and 3-13 do not have high abundance in the complex 5-node motifs. It is striking that regions with many contralateral connections like AGl and AGm contribute significantly more often to the symmetricINmotif. Surprisingly, the frequency of the MGP-contribution to the symmetricIN-motif is maximal which means that MGP is most often involved in building the symmetricINmotif.
Principal component analysis Principal component analysis (PCA) was performed to determine significant parameters out of a set of 6 local parameters that reflect local connectivity (Echtermeyer et al. 2011). In the PCA-plane the thalamic regions have a relatively small component 1 and a larger component 2, while limbic regions tend to be localized in the center of the PCA-plane and BG-core regions appear to have a relatively larger component 1 (Fig. 9). Again, this distribution of regions can be interpreted as a grouping with functional meaning. Within the PCA-analysis it was found that the CPu has a very large amount of direct connections to most of the BG2 regions (Fig. 9c). Knotty centers and rich-club analysis The knotty-center in the bilateral BG2 network contains the cortical regions AGl and AGm, followed by the CPu and AC as well as MRF and VTA (Appendix Fig. 21). By comparing the knotty-center with the rich-club regions (Appendix Fig. 20) of the bilateral BG2 network, the SNR is the sole region which is a rich-club member, however, not in the knotty-center. Some of these regions were also found in the K-core analysis. This analysis has assigned in the ipsilateral BG2 network the CPu, STh, MGP, AC, PF, MRF, VTA, SNR and SNC to the inner core and in the bilateral BG2 network AGm was assigned to the interior core, however, AGl, SNR, SNC and STh to the first shell. Partly, these groupings provide comparable results for CPu, AC, MRF and VTA, however, some differences for SNC, SNR, AGl, AGm, MGP and STh were found as well. Representation and interpretation of ipsi- and contralateral connections So far, no evidence has been reported of an asymmetry of structural connectivity in the rat brain. However, Chida and Toyosawa (1994) reported a functional asymmetry of the rat hippocampus using small animal EEG. Only, a few studies exist that describe asymmetries of olfactoric processing (Riddle and Purves 1995) and the volume of whisker representing fields of the primary somatosensory cortex in rats (Parthasarathy and Bhalla 2013). In the peripheral nervous system of the rat a significant asymmetry of connections exists, especially for the innervation of impar organs of the abdomen. To allow the representation of peripheral and central connections for extensive connectivity analysis of the complete nervous system of an organism it is mandatory to realize side-dependent connectivity mapping. Finally, the well-known structural as well as functional asymmetry of connections of the human
123
Brain Struct Funct
CNS support the concept of a generic framework allowing left- and right-side representations of connectivity. Organization of connections in a hierarchy of regions A further important challenge of the approach described before by Schmitt and Eipert (2012) is the hierarchical organization of subdivisions of regions in the rat connectivity project in neuroVIISAS. The concept of hierarchy is based on spatial scales and semantic relations. The concept is extended by defining variants of the hierarchy to map connections in terms of functional subdivisions. This approach allows the consistent representation of connections of very fine details such as terminal fields as well as connections of functional subsystems like, e.g., the sensory and motoric barrel cortices. Such a scale-based organization of connectivity data has several advantages with regard to multiscale connectivity analysis and connectome based population simulations. By controlling the granularity of a connectome through spatial scales it is also possible to filter and compare connectomes at coarse scales where most connectivity data have larger reliability values. This scaling-technique can be used in combination with functional and chemoarchitectonic data (data are not presented here) which were collated along with the connection data. Such an exceptional database of experimental raw data could be the basis for more realistic large scale simulations of the rat physiom [http://neuroviisas.med.unirostock.de/ratphysiome.html (Bai et al. 2006)]. One important perspective is to use these tract tracing based connectome data of the BG and to combine it with functional attributes at the microcircuit level (Diesmann et al. 1999) to define more realistic and population-based models for simulation studies (Gewaltig and Diesmann 2007). As shown before by (Schmitt and Eipert 2012) this step can be performed also in neuroVIISAS using a hybrid neuron model (Thibeault and Srinivasa 2013) in the NEST simulation engine (Gewaltig and Diesmann 2007).
connectivity reflects, at last partly, functionality. We have worked out a problem by pathway analysis termed the routing problem because multiple pathways from a particular source to a target were found in the data. So far, it had not been possible by network analytical methods to determine a pathway between a source and a target node that is the principal route or a restricted set of functional routes in vivo. Now a comprehensive data set is available and can be shared by the neuroscience community to be used for simulation studies in computational neuroscience as well as in connectomics and neuroinformatics. Acknowledgments The authors thank Klaus-Peter Schmitz (Department of Biomedical Engineering, University of Rostock) for the support of the neuroVIISAS project. We thank Frauke Winzer, Susanne Lehmann, Antje Schu¨mann, Jennifer Meinhardt, Ann-Christin Klu¨nker for their faithful work on extending the database and mappings. All work was supported by the Faculty of Mathematics and Natural Sciences and of the Faculty of Medicine of the University of Rostock.
Appendix A The appendix contains in the first part figures, in the second part tables and the third part formal definitions of matrices, graph-theoretical parameters and randomization models. See Figs. 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26.
Summary The connections of the basal ganglia and important regions directly connected with the basal ganglia were collated in a standardized and systematic way from all available articles which have documented results of tract-tracing experiments in normal juvenile and adult rats. The connectome was analyzed with regard to laterality in a ipsilateral and contralateral way as well as with regard to connection weights in a weighted and non-weighted configuration. Functional groups like thalamic, limbic, basal ganglia core regions and a motor associated region can be deduced from the connectional patterns by multivariate methods. Hence,
123
Fig. 10 Ordinal categories of connection strengths are plotted on the abscissa and the hypothetical estimates of connection densities shown on the left ordinate which scales the blue straight line. The right ordinate is linearly scaled for the red graph. The largest slope lies between the semiquantitative values 3.5 and 4. This logarithmic transformation has been used for adapting the ordinal categories to a distribution of connection densities that was observed by others as well [Hilgetag and Grant (2000); Markov et al. (2011, 2014); Oh et al. (2014); Ercsey-Ravasz et al. (2013); Wang et al. (2012)]. The function for transforming ordinal categories of connection strengths to the logarithmic distributed connection densities is f ðxÞ ¼ elnð10Þx4lnð10Þ
Brain Struct Funct
Fig. 11 The regions of the BG2 network are organized in a hierarchy. It is presented as a triangle hierarchy (top) and a region hierarchy (left). The sequence of regions of the hierarchy (from left to right) is
the same as that in the adjacency matrix. All regions of the connectomes that were investigated are organized in a hierarchy
Appendix B
disjunctive region that contains neuron perikarya (sources of physiological action potential) and/or axonal terminals (targets of physiological action potentials). Set of indexed nodes: The set of all indices of nodes is
Tables 7, 8, 9, 10, 11, 12, 13, 14, 15.
Appendix C: Formal definitions of parameters, matrices and concepts The definitions of expressions, parameters, matrices and simulation models (random graph models) used in this article are summarized in the following. More detailed descriptions, algorithms and proofs are provided elsewhere (Newman 2010; Rubinov and Sporns 2010; Newman et al. 2006; Ko¨tter 2007; Brandes and Erlebach 2005; Ko¨tter 2003). Basic definitions Node, vertex: The smallest subunit of a network. With regard to connectomes a node is a circumscribed or
N ¼ f1; 2; 3; . . .; ng
ð2Þ
Number of nodes: The number of nodes (regions, vertices) is n ¼ jNj
ð3Þ
Edge: A directed edge ði; jÞ 2 N N is the line that connects vertices i and j with source i and target j. The set of directed edges E is E NN
ð4Þ
Edges: The number of edges (connections, links) is ¼ jEj
ð5Þ
Set of edges:
123
Brain Struct Funct
Fig. 12 Chord diagram visualization of the bilateral BG1 network using CIRCOS. The diagram consists of 4 nested rings. The outer ring of arcs indicates the relative frequency of inputs and outputs to a particular region. The second ring is the input arc ring and the third ring is the output ring of arcs. The inner ring indicates the absolute
numbers of connections followed by a thin layer of arcs with an input part and output part for each region. The extensive connectivity of the CPu with SNC and SNR is visualized prominently. Stronger contralateral connections can be easily found by the thicker red connections between the CPu and contralateral AGl
L ¼ fði; jÞ 2 Eji 6¼ jg
Weighted matrix: The weighted matrix W is
ð6Þ
W ¼ ðwij Þni;j¼1
The set of all not self-referencing edges is ‘ ¼ jLj
ð7Þ
Graph: G ¼ ðN; EÞ
ð8Þ
Adjacency matrix: The adjacency matrix (connectivity matrix) A is A¼
ðaij Þni;j¼1
123
where aij ¼
1
if ði; jÞ 2 E
0
else
ð9Þ
ð10Þ
whereas wij is the weight of the edge ði; jÞ that connects i and j. 0 wij 1. Path: A sequence of vertices (v1 ; . . .vk ) is a path from (v1 to vk ) if 8i 2 f1; . . .; k 1g : ðvi ; viþ1 Þ 2 E. The length of a path v1 ; . . .; vk is k 1. Distance matrix: The distance matrix D is D ¼ ðdij Þni;j¼1 where
ð11Þ
Brain Struct Funct
Fig. 13 Chord diagram visualization of the bilateral BG2 network using CIRCOS. The amount of connections is relatively large. The strong connections of the CPu and AC are diverging into many different regions
dij ¼ dði; jÞ length of the shortest path from i to j; ¼ 1;
if such a path exists else
ð12Þ Generalized topology matrices (GTOM): Let Nm ðiÞ be the m step neighborhood of node i: Nm ðiÞ ¼ fj 6¼ ij minfdði; jÞ; dðj; iÞg mg
ð13Þ
Nmout ðiÞ ¼ fj 6¼ ijdði; jÞ mg
ð14Þ
Nmin ðiÞ ¼ fj 6¼ ijdðj; iÞ mg
ð15Þ
then the GTOM-matrix of step m is defined as
n m GTOMðmÞ ¼ ðgm ij Þi;j¼1 with gij 8 jNm ðiÞ \ Nm ðjÞj þ aij aji > > > > < minfjNm ðiÞj; jNm ðjÞjg þ 1 aij aji ¼ > > > > : 1
if i 6¼ j if i ¼ j
ð16Þ out The definitions of GTOMin ðmÞ and GTOMðmÞ with the directed in out m step neighborhoods Nm ðiÞ and Nm ðiÞ are analog. Degree all (degall , DGAll ): Self-references of nodes are not considered for all three degree measures. degðiÞ ¼ degall ðiÞ
123
Brain Struct Funct
Fig. 14 The distance matrix of the ipsilateral BG2 network. The distances between regions correspond to the smallest number of connections that are necessary to connect pairs of regions. The thalamic paths are in most cases longer (lighter gray values) than paths of midbrain regions or MRF
Fig. 15 The communicability matrix of the ipsilateral BG2 network. The number of shortest paths through a pair of nodes is presented by the communicability matrix (Estrada and Hatano 2008). Light gray values are indicating large communicability values. The CPu, AC, STh, PF, MRF and VTA are contributing to shortest paths of the BG2network
degðiÞ ¼
n X j¼1 j6¼i
Degree out:
123
aij þ aji
Fig. 16 The generalized topology overlap matrix (GTOM) of the ipsilateral BG2 network. The pairwise interconnectedness in relation to the number of neighbors that a pair of nodes share in common is presented by the GTOM. Large values indicate many similarities of neighbors and connections of a pair of nodes. E.g., the pair CL and CPu has a relatively large GTOM value like the pair STh and VTA
Fig. 17 The variability of connection strengths of the bilateral BG2 network. Large-standard deviations of connection weights are indicated by light gray values. The variability of connection weights of the SNC has relatively smaller values in comparison to VTA, HIPP and AGl
ð17Þ degout ðiÞ ¼
n X j¼1 j6¼i
aij
ð18Þ
Brain Struct Funct
Fig. 18 The global parameters of the ipsilateral BG2 network. The first row (beginning with ‘‘Nodes’’) presents the general parameters of the BG2 network. The line density is 52.333 %. If each regions would be connected with each other region then the line density is 100 %. Cy number of shortest cycles and CyC is the cycle coefficient, AvgHD is the average hierarchical level of the regions in the BG2 network. Six-different randomizations have been performed. Each
Degree in: degin ðiÞ ¼
n X
Latout ðiÞ ¼ aji
ð19Þ
j¼1 j6¼i
Reciprocal edge count RecðiÞ: RecðiÞ is the number of reciprocal edges adjacent to a node i. X RecðiÞ ¼ aij aji ð20Þ j2N i6¼j
degIPSI out ðiÞ degout ðiÞ
ð23Þ
Laterality of the reciprocal edge count LatRec ðiÞ: The laterality of the reciprocal edge count is the fraction of ipsilateral reciprocal edges. For a node i 2 N IPSI N. X 1 LatRec ðiÞ ¼ aij aji ð24Þ RecðiÞ IPSI j2N i6¼j
Neighborhoods: Out-neighbors of i: IPSI
N be a subset of nodes and Laterality: Let N degIPSI ðiÞ the degree of the node i in the subset N IPSI . Then the lateralities are defined as follows: degIPSI all ðiÞ degall ðiÞ
ð21Þ
degIPSI in ðiÞ Latin ðiÞ ¼ degin ðiÞ
ð22Þ
Latall ðiÞ ¼
randomization was repeated 1,000 times using exactly the same number of regions and connections as in the real BG2 network. The global network parameters of the real BG network are shown in the first column followed by mean values of the randomizations. The rewiring randomization is showing the smallest differences with regard to global parameters when comparing with the real network
Niout ¼ fj 2 Nnfigjaij ¼ 1g
ð25Þ
In-neighbors of i: Niin ¼ fj 2 Nnfigjaji ¼ 1g
ð26Þ
All neighbors of i: Ni ¼ Niout [ Niin
ð27Þ
123
Fig. 19 The local parameters of the ipsilateral BG2 network in parallel coordinates. Parallel coordinate visualization provides an overview of a high-dimensional parameter space. Parameters are sorted by similarity. The CPu, AGl and AGm are highlighted by dashed lines. The CPu has almost largest or smallest values with regard to centrality measures. For abbreviations of parameters see Appendix
Brain Struct Funct
123
Niþ ¼ Ni [ fig
ð28Þ
Network parameters Communicability matrix G: Gpq ¼
1 ðAk Þ X pq
k!
k¼0
¼ ðeA Þpq
ð29Þ
Modularity measure: Let M ¼ fM1 ; . . .; Mm g be a partition of N. Mi is a group, module or cluster of vertices. With 1 X ðajk þ akj Þ; ei ¼ ð30Þ ‘ j;k2M i
j\k
the fraction of edges that fall within group Mi N and 1X X ðajk þ akj Þ; ai ¼ ð31Þ 2‘ j2M k2Nnfjg the fraction of ends of edges that are attached to vertices in group Mi , the Modularity Q¼
m X
ðei a2i Þ;
ð32Þ
i¼1
whereas a2i is the fraction of edges that would connect vertices within group Mi if they were connected at random. A large modularity implies that the fraction of edges that fall within groups is larger than expected in the random case. The partition is generated by a ‘‘Greedy’’ optimization algorithm. Starting with a partition where every single node has its own group, stepwise those two groups are joined that increase Q most. The algorithm ends if there are no more such groups. The weighted case is similar, only the aij are replaced by wij and ‘ is replaced by the sum of the edge weights X ! wij ‘w ¼ ð33Þ i;j2N i6¼j
The method of Newman and Girvan (2004) was used. Global efficiency GE: X 1 1 GE ¼ ð34Þ nðn 1Þ i;j2N dðijÞ i6¼j
! ! GE and GE w analog with d! ði; jÞ and d w ði; jÞ Directed global efficiency: X 1 1 GE! ¼ nðn 1Þ i;j2N d! ðijÞ !
i6¼j
Harmonic mean HM:
ð35Þ
Brain Struct Funct
Fig. 20 The rich-club regions of the ipsi and contralateral BG2 network have a rich-club coefficient / [ 0:8 respectively a rank of 44: The direct neighbors or feeders are shown in the upper arc. Their interconnections are presented by blue edges. The connection of
HM ¼
1 GE
ð36Þ
The directed and weighted versions use the directed and weighted global efficiencies. Local efficiency: The local efficiency indicates how strong neighbors of nodes are interconnected. For each node i the inverse lengths of the shortest paths of the neighbors of i that are passing i are added. The local efficiency is this sum divided by the maximal possible sum of paths between neighbors that are connected with i. The efficiency of the network (global efficiency) is the average local efficiency of all nodes. Directed local efficiency: P 1 j;k2Ni djk ðNi Þ X 1 j6¼k ! LE ¼ ð37Þ n i2N ni ðni 1Þ
rich nodes with feeder nodes is shown by yellow edges and the interconnectivity of rich nodes by black edges. The indices _L and _R are indicating the sides of the hemispheres
Weighted directed local efficiency: ! P 1 w j;k2Ni djk ðNi Þ X ! 1 j6 ¼ k LE w ¼ n i2N ni ðni 1Þ
ð38Þ
ni [ 1
! whereby ni ¼ jNi j and djk ðNi Þ, respectively, djkw ðNi Þ is the length of the shortest path between j and k that contains only neighbors of i. Directed assortativity coefficient r ! : P r
!
¼ 1 2
P
ði;jÞ2L
degout ðiÞ degin ðjÞ 4‘1
2 ði;jÞ2L ðdegout ðiÞ
hP
þ degin ðjÞ2 Þ 4‘1
ði;jÞ2L ðdegout ðiÞ
hP
þ degin ðjÞÞ
ði;jÞ2L ðdegout ðiÞ
i2
þ degin ðjÞÞ
i2
ð39Þ ! Directed and weighted assortativity coefficient r w :
ni [ 1
i2 w w w ðdeg ðiÞ þ deg ðjÞÞ ij ! out in ði;jÞ2L rw ¼ P hP i2 2 2 w w w w 1 1 w ðdeg ðiÞ þ deg ðjÞ Þ w ðdeg ðiÞ þ deg ðjÞÞ ij ij out in out in ði;jÞ2L ði;jÞ2L 2 4‘ P
w w 1 ði;jÞ2L wij ðdegout ðiÞ degin ðjÞÞ 4‘
hP
ð40Þ
123
Brain Struct Funct
Fig. 21 The regions of the ipsi- and contralateral BG2 network that belong to the knotty-center in an arc diagram visualization. The knotty-center regions are on the lower arc. Direct neighbors of knotty centers are shown in the upper arc. Their interconnections are
presented by blue edges. The connections of the knotty-center with their direct neighbors are shown by yellow edges and the interconnectivity of rich-nodes by black edges. The indices _L and _R are indicating the sides of the hemispheres
Fig. 22 The K-core analysis of the ipsilateral BG2 network. The connections of CPu are highlighted. In the outer circles are those regions which have the smallest number of connections and in the center are regions that have the largest number of connections
Fig. 23 The K-core analysis of the ipsi- and contralateral BG2 network. The connections of the CPu are highlighted
123
Brain Struct Funct
The correlation of the degrees of nodes that are connected: 1 r 1. Large positive values imply that nodes are mainly connected to nodes with similar degrees. Large negative values imply that nodes with a large degree are mainly to nodes that have a small degree. If r 0 there is no relation detectable. Average path length = characteristic path length ðdÞ: With P ¼ fði; jÞ 2 N Njdði; jÞ\1g, the set of paths. 8 X > < 1 dði; jÞ; P 6¼ £ ð41Þ d ¼ jPj ði;jÞ2P > : 0; P¼£ In the weighted case the distances dði; jÞ are replaced by the weighted distances. Average directed degree deg: deg ¼
Fig. 24 Metric multidimensional scaling (MDS) applied to the ipsilateral BG2 network. Regions that are adjacent to each other have more similarities with regard to their connectivity than with regions that are further away from each other. The 8 thalamic regions that are located within the lower-left quadrant are relatively similar with regard to connectivity. AGm and AGl in the lower-right quadrant have also a large similarity. Interestingly, SNR is adjacent to MRF, however, the second major output node of the BG is MGP which is not to far from SNR. Limbic regions like Ac, Ent, HIPP and Pir are located in the upper left quadrant
2‘ n
ð42Þ
Heterogeneity VC: Coefficient of variation (VC) of the degreeall parameter. X 1 rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ðdegall ðiÞ degÞ2 HVC ¼ ð43Þ deg i2N If HVC ¼ 0, all nodes have the same degree. The larger HVC the more diverse are the node degrees. In the weighted case the versions of the degrees are used. The heterogeneity measure of Estrada (2010) was not implemented because it is not defined for directed and weighted graphs. Line density Ld: Ld ¼
‘ n ðn 1Þ
ð44Þ
Without self-referencing edges. Rich-club coefficient /ðkÞ: With Nk ¼ fi 2 NjdegðiÞ [ kg and Ek ¼ fði; jÞ 2 Nk Nk jði; jÞ 2 Eg Gk ¼ ðNk ; Ek Þ is the subgraph of G ¼ ðN; EÞ containing all vertices with a degree greater than k. The rich-club coefficient /ðkÞ of a graph G is defined as the line density of the subgraph Gk : /ðkÞ ¼ LdðGk Þ
ð45Þ
Diameter Diam: Diam ¼ maxfdði; jÞjdði; jÞ\1g
ð46Þ
Katz index: The Katz index (Katz status index, Katz centrality) is a measure for the direct and indirect input of a node (Foster et al. 2001). CKatz ðiÞ ¼ Fig. 25 MDS applied to the bilateral BG2 network. The major output regions SNR and MGP are located in the upper- and lower- left quadrant, however, not too far from each other. Most limbic and all thalamic regions are located in the upper-left quadrant (left hemisphere) or in the upper-right quadrant (right hemisphere)
1 X X
ak ðAk Þji
ð47Þ
k¼1 j¼1
The attenuation factor a has to be smaller than the reciprocal of the absolute value of the largest eigenvalue of A. For a better readability and comparability of the results, in
123
Brain Struct Funct Fig. 26 Counts of cyclic connections or loops in the bilateral BG2 network. The logarithm (Log) of the frequency of contributions of regions to a cycle of size 2–7 is presented by the y-axis. Every second region (all regions of the left hemisphere) are labeled. For a cycle from MRF back to MRF passing 4 regions within the BG2 network there exist 8,255 possible paths or 41,191,973 paths when 7 regions are allowed
neuroVIISAS the Katz centrality is multiplied by the mean in ðiÞ of all nodes with CKatz ðiÞ [ 0. Hence, of the quotient Cdeg Katz ðiÞ
the values lie in the same range as the indegrees. Number of triangles: X ðaij þ aji Þðaik þ aki Þðajk þ akj Þ t! ðiÞ ¼ j;k2Nnfig j\k
ð48Þ
The maximum number of possible triangles that can be deviated from a complete reciprocal triangle is 8. Weighted number of triangles: X 1 1 1 1 1 1 ! t w ðiÞ ¼ ðw3ij þ w3ji Þðw3ik þ w3ki Þðw3jk þ w3kj Þ ð49Þ j;k2Nnfig j\k
Instead of the sum of triangles (t! ðiÞ) the sum of geometric means of edge weights of each triangle is calculated. The 1
following example provides ðwij wjk wik Þ3 as the summand:
Directed transitivity: The general definition of transitivity (T) is the sum of number of triangles around all nodes divided by the maximum possible sum of triangles around all nodes. P t! ðiÞ ! T ¼ P i2N ð50Þ i2N tmax ðiÞ Directed and weighted transitivity: ! P ! t w ðiÞ w i2N T ¼P i2N tmax ðiÞ
ð51Þ
whereby tmax ðiÞ ¼ degðiÞ ðdegðiÞ 1Þ 2 recðiÞ with degðiÞ ¼ number of adjacent edges of i and recðiÞ ¼ number of reciprocal edges of i (the two directions of one reciprocal edge are considered as one reciprocal edge). The degree deg and the reciprocity rec are defined as: X degðiÞ ¼ aij þ aji ð52Þ j2Nnfig
recðiÞ ¼
X
aij aji
j2Nnfig
For the directed and weighted case: 1 wij [ 0 aij ¼ 0 else
123
ð53Þ
ð54Þ
Brain Struct Funct Table 7 Local parameters of the weighted-ipsilateral BG2 network Region
DGA
DGO
DGI
CDC
Katz
SPC
Triag
CyclC
EccO
CPu
26.25
12.25
14.0
0.53333
12.17519
0.16667
503.6196
0.66462
375
VTA
20.875
12.5
8.375
0.4012
8.04604
0.25
436.14384
0.55527
0.41667
SNR
19.125
11.0
8.125
0.42484
8.51141
0.29167
323.26377
0.66477
0.41667
MRF
18.875
12.375
AC
18.0
9.625
6.5
0.34437
6.58671
0.25
399.3736
0.56189
0.5
8.375
0.46528
7.82869
0.25
402.70699
0.49854
0.45
SNC
17.375
11.125
6.25
0.35971
6.90902
0.25
317.85196
0.59338
0.5
STh
16.875
9.0
7.875
0.46667
7.93142
0.25
313.32162
0.64057
0.41667
AGm
16.25
9.75
6.5
0.4
6.16135
375
246.8733
0.51577
0.41667
MGP
15.5
8.75
6.75
0.43548
6.93033
0.33333
270.12611
0.59842
0.41667
AGl
14.5
7.625
6.875
0.47414
6.36241
0.16667
237.70384
0.56222
0.5
PF
13.875
8.5
5.375
0.38739
5.68001
375
263.27964
0.5863
0.5
Ac
12.75
6.25
6.5
0.5098
6.29047
0.25
227.32308
0.55463
0.5
Pir
11.0
4.5
6.5
0.59091
6.51555
0.25
157.8558
0.56038
0.54167
VM
10.375
4.375
6.0
0.57831
6.2192
0.33333
154.06266
0.5957
0.5
HIPP CM
10.25 10.0
5.25 4.75
5.0 5.25
0.4878 525
4.75053 5.51932
0.25 375
133.31931 131.01033
0.57456 0.59741
0.5 0.5
LHb
10.0
4.875
5.125
0.5125
5.42794
0.25
164.73628
0.50748
0.5
VL
9.0
1.5
7.5
0.83333
7.67315
0.2
81.79831
0.73374
625
MDL
8.75
1.5
7.25
0.82857
7.35711
0.41667
101.08281
0.58321
625
LGP
8.75
2.625
6.125
0.7
6.23351
0.41667
106.87098
0.65648
0.54167
CL
8.5
4.0
4.5
0.52941
4.73061
0.41667
96.16541
0.62436
0.5
Ent
8.375
3.25
5.125
0.61194
4.82634
0.25
113.68525
0.43629
0.66667
PC
7.0
2.5
4.5
0.64286
4.89828
0.41667
71.94521
0.60614
0.54167
MDM
6.0
1.5
4.5
0.75
4.34144
0.58333
56.44901
0.52978
625
VA
2.5
1.0
1.5
0.6
1.51929
625
7.97699
0.70278
625
Lev
Loc
Region
EccI
CCO
CCI
CCA
CCT
CC2
ADnb
VCDG
CPu
0.41667
0.29079
0.26399
0.25321
0.31398
0
12.54348
366
0.3729
0.96524
VTA
0.5
0.27075
0.30556
0.27075
0.27191
0
12.92935
0.37537
0.2589
0.98554
SNR
0.5
0.25804
0.41506
0.25952
0.33326
0.4375
12.70238
0.43298
0.23739
0.8122
MRF
0.5
0.2747
0.37619
0.2747
0.29024
0
13.0163
0.38224
0.21062
0.98553
AC
0.5
0.28436
0.30362
0.28436
0.23942
0
13.26136
0.37947
0.17976
0.93529
SNC
0.5
0.28185
0.4625
0.28185
0.32768
125
13.2381
0.39005
0.1655
0.86018
STh
0.45
0.33686
0.40545
0.32557
0.37036
0.03333
14.38889
0.34253
0.10538
0.74536
AGm
0.5
0.23989
0.27574
0.25
0.22607
375
12.54605
0.45852
0.17302
0.68563
MGP
0.54167
0.3102
0.42468
0.29496
0.31854
0.23125
13.50658
0.40979
0.1091
0.76064
AGl
0.5
0.3125
0.27552
0.29412
0.26353
0.24405
13.53676
0.43892
0.08385
0.64739
PF
0.5
0.34835
0.45455
0.34835
0.35869
0.04167
14.90441
0.32782
-0.00814
0.68479
Ac Pir
0.5 0.5
0.37099 0.32639
0.34286 0.36882
0.31127 0.33869
0.30886 0.31955
0.3 0.20238
13.70139 14.35294
0.4002 0.36632
0.00466 -0.09458
0.72482 0.68825
VM
0.5
0.41667
0.42898
0.38805
0.37946
0.11528
15.75893
0.29619
-0.18057
0.54
HIPP
625
0.31439
0.41319
0.33333
0.32676
0.20139
14.39167
0.37484
-0.12721
0.56576
CM
0.5
0.31389
0.42778
0.34345
0.35408
0.0625
15.2
0.32213
-0.17637
0.52686
LHb
625
375
0.46023
0.3875
0.3362
0.16493
15.325
0.31079
-0.18284
0.62755
VL
0.41667
625
0.40568
0.40568
0.46476
0.11298
17.03409
0.25543
-0.28858
0.38215
MDL
0.5
0.34375
0.40064
0.40064
0.38289
0.07273
16.52885
0.26345
-0.28663
0.4676
LGP
0.54167
0.64375
0.43523
0.45265
0.46065
0.16761
16.5
0.28107
-0.28201
0.49548
CL
0.54167
0.41295
0.49107
0.45265
0.41451
0.14489
16.75
0.25151
-0.30743
0.48319
123
Brain Struct Funct Table 7 continued Region
EccI
CCO
CCI
CCA
CCT
CC2
ADnb
VCDG
Lev
Loc
Ent
0.58333
0.34201
0.32576
0.35897
0.2814
0.18977
15.03846
0.34071
-0.24859
0.47373
PC
0.5
0.40417
0.47098
0.4566
0.41829
0.13929
18.22222
0.20137
-0.43369
0.3041
MDM VA
0.58333 0.66667
375 625
0.38611 0.16667
0.38611 0.4625
0.32442 0.39885
0.23214 0.22807
15.2875 18,325
0.26926 0.23199
-0.41216 -0.75126
0.34497 0.14264
Region
CEO
CEI
BC
EC
SC
PRC
FC
Stress
CPu
4.20438
4.4376
0.23583
1
243.76129
1
0.50716
182
1.41218
VTA
3.94521
3.25424
0.05577
0.91751
154.20322
0.58706
0.39402
45
1.86591
SNR
3.86577
3.65482
0.0939
0.69998
130.14345
0.53169
0.33929
72
0.64769
MRF
3.86577
2.88
0.01208
897
129.7987
0.41633
0.31682
10
0.78248
AC SNC
3.37632 3.91837
3.13043 3.40828
0.04026 0.02159
0.68099 0.80884
107.60925 116.2092
0.64706 0.38099
0.38934 0.2869
30 26
0.65307 0.74325
STh
3.81457
3.42043
0.0237
0.77207
125.52561
0.4873
0.38685
23
0.26545
AGm
3.72574
2.83744
0.07555
0.67544
88.77745
0.70423
0.38213
64
1.49691
MGP
3.34884
3.18232
0.0226
0.59287
86.04644
0.43381
0.31835
22
-0.02124
AGl
3.40828
2.83744
0.01736
0.66608
85.19876
0.64031
0.38948
17
1.02062
PF
3.40828
2.8263
0.00645
0.72779
99.67888
0.32626
0.31893
10
0.45657
Ac
3.21788
2.83744
0.00543
0.56331
77.21747
0.43082
0.27614
7
-0.18659
Pir
2.96907
2.97214
0.01585
0.35111
45.73839
0.45033
0.21232
14
0.74636
VM
2.89447
3.22148
0.00218
0.3923
55.50867
0.36658
0.27266
5
-0.40825
HIPP
2.66667
2.67907
0.01479
0.35601
34.50464
0.35927
0.24167
10
0.37318
CM
3.16484
2.86567
0.00494
0.39959
45.3308
0.30888
0.22619
7
-0.64639
LHb
3
2.71698
0.00453
0.42068
54.47565
0.29981
0.25625
5
-0.65307
VL
2.26772
3.64557
0.00595
0.17041
24.99282
0.46127
0.16023
5
-1.12268
MDL
2.16867
3.23595
0.00177
0.14331
22.56384
0.43209
0.14423
4
-1.49273
LGP CL
2.67907 2.78261
3.08021 2.8263
0 0.00203
0.26767 0.35679
38.75035 35.96419
0.35154 0.27445
0.20786 0.24242
0 6
-0.5946 -0.69016
Ent
2.37037
2.65438
0.00242
0.26855
26.31024
0.36863
0.25
2
PC
2.65438
2.75598
0.00112
0.26166
28.23121
0.28222
0.2934
4
-1.12268
MDM
2.18347
2.75862
0.00192
0.10016
9.94491
0.28595
0.16944
4
-1.30614
VA
2.19847
2.20015
0
0.11114
5.36512
0.1381
0.15
0
-2.21914
Region
ZI
ZA
PCO
PCI
PCA
RadO
RadI
ZO
0
CenO
CenI
CPu
1.63169
1.63651
0.64973
0.64541
0.64907
1.42882
1.44132
3
10
VTA SNR
1.45048 0.46465
1.89849 0.59929
0.6392 0.66219
0.63711 0.61917
0.63839 0.6512
1.41319 1.40799
1.35938 1.39306
-5 -4
-16 -10
MRF
-0.42875
0.49136
0.65177
0.65385
0.65471
1.40799
1.31944
-5
-19
0.72524
0.73019
0.66352
0.66251
0.66397
1.37049
1.34722
-11
-17
AC SNC
-0.31337
0.23626
0.64083
0.5952
0.62771
1.41146
1.37326
-3
-11
STh
0.65916
0.49556
0.63735
0.59965
0.63473
1.40451
1.37431
-8
-15
AGm
0.77174
1.44677
0.65976
0.66198
0.66225
1.39826
1.31424
-9
-17
MGP
0.46465
0.23626
0.66245
0.55967
0.63541
1.36806
1.35243
-11
-16
AGl
1.37199
1.25569
0.59446
0.62281
0.6088
1.37326
1.31424
-11
-18
PF
-0.79964
-0.17863
Ac
0.48349
Pir
-0.96699
VM
0.1715
HIPP
0.48349
123
0.63149
0.64467
0.64216
1.37326
1.31285
-12
-18
0
0.5632
0.63018
609
1.3559
1.31424
-18
-17
0.29208
0.27778
0.64793
0.58264
1.32986
1.33021
-20
-16
-0.27298
0.65959
0.63281
0.6538
1.32118
1.35625
-18
-19
0.43811
0.52608
0.54
0.55205
1.29167
1.2934
-16
-15
Brain Struct Funct Table 7 continued Region
ZI
ZA
PCO
PCI
PCA
RadO
RadI
CenO
CenI
CM
-1.62924
-1.03731
0.65374
0.63492
0.66375
1.35069
1.31771
-16
-18
LHb
-1.69223
-1.02226
0.65746
0.56871
0.62625
1.33333
1.29861
-16
-18
VL MDL
0.77174 -0.48349
-0.65514 -1.31434
0.44444 0.5
0.65778 0.65161
0.64815 0.66286
1.22569 1.20556
1.39236 1.35764
-23 -22
-12 -14 -19
0.65916
0.02881
0
0.46647
0.35714
1.2934
1.34201
-25
CL
LGP
-0.89689
-0.85283
0.53125
0.54938
0.55536
1.30729
1.31285
-21
-22
Ent
0.96699
0.29208
0.27219
0.51636
0.44108
1.24479
1.28993
-21
-16
PC
-1.02899
-1.22839
0.66
0.64815
0.65561
1.28993
1.30382
-22
-21
MDM
-0.96699
-1.31434
0.44444
0.62346
0.60069
1.20868
1.30417
-24
-19
VA
-1.86942
-2.20123
0.5
0.44444
0.48
1.21181
1.21215
-25
-25
The regions are sorted using the weighted degree all. Abbreviations and parameter definitions are described in the Supporting Information A all, AD average degree, C circle, CC cluster-coefficient, CE closeness centrality, DG degree, I in, O out, SP length of shortest path
Table 8 Local parameters of the bilateral BG2 network Region
HD
REC
DGA
DGO
DGI
CDC
LA
LO
LI
LR
Katz
SPC
MRF_R
12
24
67
41
26
0.39
0.55
0.54
0.58
0.63
25.16
1
MRF_L
12
24
68
42
26
0.38
0.56
0.55
0.58
0.63
25.2
1
AGm_R
14
21
55
27
28
0.51
0.58
0.59
0.57
0.67
26.34
1
AGl_R
14
22
48
24
24
0.5
0.6
0.58
0.63
0.59
23.25
1
AGm_L
14
22
56
28
28
0.5
0.61
0.61
0.61
0.68
26.39
1
AGl_L
14
23
50
25
25
0.5
0.62
0.6
0.64
0.61
23.96
1
VL_L
12
4
22
5
17
0.77
0.64
0.6
0.65
0.75
18.12
1
VL_R
12
4
22
5
17
0.77
0.64
0.6
0.65
0.75
18.08
1
HIPP_L
10
8
31
18
13
0.42
0.65
0.61
0.69
0.75
11.15
1
HIPP_R
10
8
31
18
13
0.42
0.65
0.61
0.69
0.75
11.11
1
AC_R
8
26
64
37
27
0.42
0.66
0.59
0.74
0.77
24.82
1
CM_L
11
8
30
13
17
0.57
0.67
0.77
0.59
0.63
19.12
2
LHb_L
10
10
33
16
17
0.52
0.67
0.69
0.65
0.8
16.61
1
CM_R
11
8
30
13
17
0.57
0.67
0.77
0.59
0.63
19.08
2
LHb_R AC_L
10 8
10 25
33 62
16 36
17 26
0.52 0.42
0.67 0.68
0.69 0.61
0.65 0.77
0.8 0.8
16.56 24.76
1 1
CPu_R
11
21
58
23
35
0.6
0.69
0.87
0.57
0.86
31.66
1
PC_L
11
7
20
9
11
0.55
0.7
0.67
0.73
0.71
13.3
2
VM_L
11
9
30
12
18
0.6
0.7
0.75
0.67
0.78
19.37
2
PC_R
11
7
20
9
11
0.55
0.7
0.67
0.73
0.71
13.27
2
VM_R
11
9
30
12
18
0.6
0.7
0.75
0.67
0.78
19.33
2
VTA_R
11
21
54
33
21
0.39
0.7
0.67
0.76
0.76
20.72
1
CPu_L
11
21
58
23
35
0.6
0.71
0.87
0.6
0.86
31.67
1
PF_R
11
15
38
21
17
0.45
0.71
0.76
0.65
0.73
17.29
2
VTA_L
11
22
56
34
22
0.39
0.71
0.68
0.77
0.77
21.41
1
Pir_R
10
6
28
12
16
0.57
0.71
0.58
0.81
0.67
16.27
1
PF_L
11
15
39
22
17
0.44
0.72
0.77
0.65
0.73
17.33
2
SNR_L
12
13
44
26
18
0.41
0.73
0.73
0.72
0.85
17.68
2
CL_L
11
5
22
10
12
0.55
0.73
0.8
0.67
0.8
13.02
2
SNR_R CL_R
12 11
13 5
44 22
26 10
18 12
0.41 0.55
0.73 0.73
0.73 0.8
0.72 0.67
0.85 0.8
17.64 12.98
2 2
123
Brain Struct Funct Table 8 continued Region
HD
REC
Ent_L
11
Pir_L
10
SNC_L SNC_R
12 12
Ac_R
11
10
34
14
20
0.59
0.76
0.86
0.7
0.9
18.73
1
STh_R
9
14
37
19
18
0.49
0.78
0.84
0.72
0.86
18.09
1
STh_L
9
14
38
20
18
0.47
0.79
0.85
0.72
0.86
18.14
1
Ac_L
11
11
35
15
20
0.57
0.8
0.87
0.75
0.91
18.75
1
MDL_L
12
5
21
5
16
0.76
0.81
0.8
0.81
0.8
16.92
2
MDL_R
12
5
21
5
16
0.76
0.81
0.8
0.81
0.8
16.84
2
VA_L
12
0
6
2
4
0.67
0.83
1
0.75
0
4.23
3
MGP_L
12
12
36
22
14
0.39
0.83
0.77
0.93
0.92
VA_R
12
0
6
2
4
0.67
0.83
1
0.75
0
MGP_R
12
12
36
22
14
0.39
0.83
0.77
0.93
0.92
13.95
1
Ent_R
12
2
6
3
3
0.5
0.83
0.67
1
1
2.76
1
MDM_R
12
3
14
3
11
0.79
0.86
1
0.82
1
11.32
2
MDM_L
12
4
16
4
12
0.75
0.88
1
0.83
1
12.1
2
LGP_L LGP_R
12 12
4 4
16 16
5 5
11 11
0.69 0.69
1 1
1 1
1 1
1 1
10.96 10.91
2 2
Region
Triag
CyclC
EccO
EccI
CCO
CCI
CCA
CCT
CC2
ADnb
VCDG
MRF_R
1,975
0.07
2
3
0.39
0.56
0.39
0.45
0.1
36.4
0.39
MRF_L
2,021
0.07
2
3
0.39
0.56
0.38
0.45
0.05
36.23
0.39
AGm_R
1,279
0.07
2
2
0.43
0.45
0.39
0.44
0.34
36.29
0.48
940
0.07
2
2
0.43
0.42
0.42
0.42
0.38
36.58
0.44
AGl_R AGm_L AGl_L VL_L
DGA
DGO
DGI
CDC
LA
LO
LI
LR
Katz
SPC
9
28
14
14
0.5
0.75
0.64
0.86
0.89
13.85
1
7
29
12
17
0.59
0.76
0.67
0.82
0.71
17.02
1
13 13
42 42
27 27
15 15
0.36 0.36
0.76 0.76
0.78 0.78
0.73 0.73
0.85 0.85
15.25 15.24
1 1
14.03
1
4.22
3
1,324
0.07
2
2
0.42
0.45
0.39
0.44
0.34
36.41
0.47
988 296
0.07 0.09
2 3
3 2
0.42 0.85
0.41 0.64
0.41 0.64
0.41 0.65
0.39 0.18
36.26 50.28
0.44 0.23
VL_R
296
0.09
3
2
0.85
0.64
0.64
0.65
0.18
50.17
0.23
HIPP_L
510
0.08
2
3
0.58
0.55
0.52
0.56
0.25
42.3
0.41
HIPP_R
479
0.08
2
3
0.52
0.54
0.52
0.52
0.25
42.17
0.41
AC_R
1,714
0.07
2
3
0.4
0.48
0.4
0.43
0.31
36.79
0.41
CM_L
522
0.09
3
2
0.61
0.65
0.57
0.61
0.18
46
0.31
LHb_L
634
0.1
3
3
0.59
0.64
0.55
0.61
0.26
42.52
0.32
CM_R
498
0.09
2
2
0.53
0.65
0.54
0.58
0.19
44.82
0.36
LHb_R
634
0.1
3
3
0.59
0.64
0.55
0.61
0.27
42.39
0.32 0.39
AC_L
1,724
0.08
2
3
0.42
0.53
0.42
0.46
0.26
37.76
CPu_R
1,423
0.07
2
2
0.5
0.42
0.41
0.44
0.28
37.43
0.4
PC_L
266
0.12
3
3
0.71
0.75
0.72
0.73
0.24
53.77
0.17
VM_L
564
0.1
3
2
0.69
0.67
0.62
0.66
0.18
47.33
0.28
PC_R
266
0.12
2
3
0.71
0.75
0.72
0.73
0.23
53.69
0.18
VM_R VTA_R
564 1,434
0.1 0.08
2 2
2 3
0.69 0.46
0.67 0.59
0.62 0.46
0.66 0.51
0.19 0.2
47.14 39.82
0.28 0.38
CPu_L
1,433
0.07
2
2
0.5
0.42
0.41
0.44
0.28
37.49
0.4 0.32
PF_R
0.1
2
3
0.58
0.65
0.56
0.61
0.22
43.91
1,514
0.08
2
3
0.46
0.56
0.46
0.5
0.2
39.59
0.37
Pir_R
418
0.08
3
3
0.49
0.66
0.56
0.56
0.28
42.36
0.36
PF_L
865
0.09
3
3
0.56
0.65
0.54
0.6
0.22
43.38
0.32
VTA_L
123
838
Brain Struct Funct Table 8 continued Region
Triag
CyclC
EccO
SNR_L
916
0.08
3
3
0.39
0.64
0.39
0.49
0.37
36.23
0.44
CL_L
293
0.1
3
3
0.64
0.68
0.63
0.65
0.26
46.71
0.27
SNR_R CL_R
916 293
0.08 0.1
3 3
3 3
0.39 0.64
0.64 0.68
0.39 0.63
0.49 0.65
0.39 0.27
36.1 46.41
0.44 0.27
Ent_L
453
0.09
2
3
0.6
0.63
0.58
0.61
0.27
44.11
0.34
Pir_L
460
0.08
3
3
0.51
0.65
0.58
0.58
0.29
42.64
0.36
SNC_L
948
0.09
3
3
0.47
0.75
0.48
0.56
0.29
39.76
0.41
SNC_R
948
0.09
2
3
0.47
0.75
0.48
0.56
0.3
39.59
0.41
EccI
CCO
CCI
CCA
CCT
CC2
ADnb
VCDG
Ac_R
641
0.1
3
3
0.63
0.56
0.51
0.58
0.37
39.21
0.41
STh_R
804
0.1
2
3
0.61
0.63
0.55
0.62
0.25
42.83
0.34
STh_L
837
0.1
3
3
0.6
0.63
0.54
0.61
0.25
42.42
0.35
Ac_L
679
0.1
3
3
0.6
0.57
0.51
0.58
0.37
39.38
0.4
MDL_L
284
0.12
3
3
0.75
0.69
0.69
0.69
0.25
48.31
0.25
MDL_R
284
0.12
3
3
0.75
0.69
0.69
0.69
0.25
48.19
0.26
20
0.15
3
3
1
0.33
0.67
0.67
0.34
49.83
0.15
VA_L MGP_L
728
0.1
3
3
0.49
0.78
0.48
0.59
0.39
37.79
0.42
20
0.15
3
3
1
0.33
0.67
0.67
0.35
49.5
0.16
MGP_R Ent_R
728 18
0.1 0.12
2 3
3 3
0.49 0.83
0.78 0.67
0.48 0.58
0.59 0.69
0.39 0.42
37.63 39
0.42 0.37
MDM_R
138
0.11
3
3
0.83
0.79
0.79
0.78
0.29
49.91
0.28
MDM_L
180
0.11
3
3
0.92
0.77
0.77
0.78
0.29
48.5
0.3
LGP_L
202
0.15
3
3
1
0.82
0.83
0.87
0.38
44.17
0.3
LGP_R
202
0.15
3
3
1
0.82
0.83
0.87
0.36
43.83
0.31
VA_R
The regions are sorted using the degree all. Abbreviations and parameter definitions are described in the Supporting Information A all, AD average degree, C circle, CC cluster-coefficient, DG degree, I input, L laterality, O out, O output, R reciprocal, SP length of shortest path
Table 9 Local parameters of the bilateral BG2 network Region
Lev
Loc
CCO
CCI
BC
EC
SC
PRC
FC
Stress
Shapley
MRF_R
0.32
0.88
0.86
0.67
0.04
0.98
25,442,898.34
0.58
0.58
735
-0.21
MRF_L
0.33
0.91
0.88
0.67
0.04
1
26,026,701.64
0.58
0.56
763
-0.21
AGm_R
0.25
0.62
0.69
0.7
0.04
0.64
17,358,891.41
0.85
0.66
649
-0.22
AGl_R
0.18
0.49
0.66
0.66
0.04
0.55
13,243,875.48
0.89
0.85
516
-0.37
AGm_L
0.26
0.63
0.7
0.7
0.05
0.66
18,011,872.88
0.83
0.68
675
-0.22 -0.4
AGl_L
0.2
0.5
0.67
0.66
0.04
0.57
14,180,107.63
0.92
0.86
552
VL_L
-0.37
0.31
0.51
0.6
0
0.15
2,727,904.72
0.41
0.26
26
0.21
VL_R
-0.37
0.32
0.51
0.6
0
0.15
2,700,497.46
0.41
0.26
26
0.21
HIPP_L
-0.1
0.43
0.61
0.54
0.01
0.51
5,787,286.52
0.34
0.45
181
0.07
HIPP_R
-0.09
0.43
0.61
0.54
0.01
0.49
5,588,129.56
0.34
0.45
217
0.03
AC_R
0.3
0.76
0.8
0.67
0.07
0.89
22,779,869.63
0.78
0.69
942
-0.35
CM_L
-0.18
0.39
0.57
0.6
0.01
0.34
6,998,306.99
0.43
0.46
129
0.09
LHb_L
-0.1
0.45
0.59
0.6
0.01
0.48
8,298,279.27
0.39
0.52
141
0.03
CM_R
-0.15
0.38
0.58
0.6
0.01
0.32
6,632,895.27
0.43
0.46
230
0.09
LHb_R AC_L
-0.1 0.27
0.46 0.75
0.59 0.79
0.6 0.66
0.01 0.04
0.48 0.89
8,215,504.87 22,732,362.73
0.4 0.73
0.52 0.68
138 698
0.03 -0.25
CPu_R
0.25
0.74
0.65
0.78
0.05
0.55
18,102,852.8
1
0.59
726
-0.37
123
Brain Struct Funct Table 9 continued Region
Lev
Loc
CCO
CCI
BC
EC
SC
PRC
FC
Stress
Shapley
PC_L
-0.45
0.23
0.54
0.55
0
0.28
3,976,306.75
0.31
0.59
34
0.16
VM_L
-0.2
0.4
0.56
0.61
0
0.35
7,402,875.66
0.43
0.49
83
0.1
PC_R VM_R
-0.45 -0.2
0.23 0.4
0.55 0.57
0.56 0.61
0 0
0.28 0.35
3,930,742.69 733,2016.86
0.32 0.44
0.59 0.49
35 82
0.16 0.1
VTA_R
0.18
0.68
0.75
0.62
0.03
0.88
18,725,331.26
0.55
0.64
429
-0.18
CPu_L
0.25
0.74
0.65
0.78
0.05
0.56
18,415,083.57
0.99
0.59
714
-0.37
-0.04
0.44
0.64
0.59
0.01
0.59
10,995,770.08
0.38
0.68
195
-0.01
0.2
0.7
0.77
0.63
0.03
0.9
19,880,271.57
0.57
0.65
485
-0.2
Pir_R
-0.16
0.44
0.56
0.58
0.01
0.29
4,746,984.35
0.45
0.4
123
0.09
PF_L
-0.02
0.46
0.64
0.58
0.01
0.61
11,479,685.57
0.38
0.65
205
-0.02 -0.03
PF_R VTA_L
SNR_L CL_L SNR_R
0.14
0.54
0.67
0.59
0.02
0.58
11,064,009.13
0.45
0.49
336
-0.33
0.31
0.54
0.56
0
0.27
3,797,240.62
0.3
0.42
67
0.17
0.14
0.54
0.67
0.6
0.02
0.57
10,923,968.28
0.45
0.49
347
-0.03
CL_R
-0.33
0.31
0.54
0.56
0
0.27
3,721,042.61
0.3
0.42
67
0.17
Ent_L
-0.19
0.37
0.58
0.57
0.01
0.4
5,650,397.26
0.42
0.55
125
0.06
Pir_L
-0.15
0.45
0.56
0.59
0.01
0.29
5,000,332.71
0.46
0.43
116
0.09
SNC_L
0.07
0.58
0.68
0.58
0.01
0.68
10,703,603.42
0.37
0.48
260
0
SNC_R Ac_R
0.07 -0.03
0.59 0.48
0.69 0.57
0.58 0.61
0.01 0.01
0.67 0.4
10,568,914.69 7,650,218.92
0.38 0.48
0.48 0.49
266 150
0 0.02
STh_R
-0.04
0.46
0.62
0.6
0.01
0.56
10,461,940.53
0.46
0.65
252
-0.02
STh_L
-0.02
0.48
0.62
0.6
0.01
0.59
10,996,004.71
0.46
0.63
263
-0.03
Ac_L
-0.02
0.49
0.58
0.61
0.01
0.43
8,161,457.86
0.48
0.52
166
MDL_L
-0.37
0.31
0.52
0.59
0
0.17
3,135,070.85
0.41
0.31
31
0.2
MDL_R
-0.37
0.31
0.52
0.59
0
0.17
3,078,873.46
0.41
0.31
32
0.2
VA_L
-0.78
0.09
0.43
0.48
0
0.05
249,113.28
0.15
0.27
1
0.47
0.02
0.49
0.64
0.56
0.01
0.57
8,355,735.62
0.38
0.54
217
0.02
-0.78
0.09
0.43
0.48
0
0.05
241,958.15
0.15
0.27
1
0.47
MGP_L VA_R MGP_R
0
0.03
0.49
0.64
0.56
0.01
0.57
8,210,767.81
0.38
0.54
229
0.02
Ent_R
-0.71
0.11
0.47
0.44
0
0.1
268,736.61
0.13
0.58
4
0.31
MDM_R
-0.54
0.22
0.48
0.54
0
0.09
1,061,603.37
0.27
0.27
8
0.27
MDM_L
-0.48
0.25
0.49
0.54
0
0.11
1,396,820.87
0.3
0.33
11
0.23
LGP_L
-0.44
0.3
0.46
0.52
0
0.14
1,710,007.94
0.29
0.39
6
0.21
LGP_R
-0.44
0.3
0.46
0.53
0
0.14
1,670,074.84
0.3
0.39
6
0.21
Region
ZO
ZI
ZA
PCO
PCI
PCA
RadO
RadI
CenO
CenI
Hub
Aut
MRF_R
1.3
0.64
1.15
0.65
0.64
0.65
2.84
2.51
-1
-10
0.98
0.87
MRF_L
1.4
0.71
1.25
0.66
0.65
0.66
2.86
2.51
1
-10
1
0.88
AGm_R
0.53
1
0.75
0.62
0.64
0.63
2.55
2.57
-15
-7
0.67
0.81
AGl_R
0.53
1.36
0.88
0.53
0.52
0.53
2.49
2.49
-18
-11
0.61
0.59
AGm_L
0.59
1.09
0.82
0.64
0.65
0.65
2.57
2.57
-14
-7
0.69
0.81
AGl_L
0.79
1.47
1.11
0.56
0.59
0.58
2.51
2.49
-17
-11
0.63
0.61
-1.24 -1.38
-0.42 -0.43
-1.03 -1.12
0.48 0.48
0.64 0.6
0.62 0.58
2.02 2.02
2.35 2.35
-39 -39
-18 -18
0.15 0.15
0.65 0.65
HIPP_L
0.72
0.52
0.72
0.62
0.52
0.59
2.37
2.14
-24
-25
0.49
0.4
HIPP_R
0.38
0.52
0.49
0.64
0.47
0.59
2.37
2.16
-24
-24
0.47
0.4
AC_R
1.42
1.5
1.64
0.64
0.55
0.61
2.76
2.51
-5
-10
0.89
0.78
CM_L
-0.63
-0.42
-0.6
0.64
0.62
0.64
2.24
2.35
-31
-19
0.35
0.65
LHb_L
-0.23
-0.8
-0.46
0.62
0.65
0.64
2.31
2.33
-27
-19
0.41
0.59
VL_L VL_R
123
Brain Struct Funct Table 9 continued Region
ZO
ZI
ZA
PCO
PCI
PCA
RadO
RadI
CenO
CenI
Hub
Aut
CM_R
-0.8
-0.43
-0.72
0.64
0.58
0.62
2.27
2.35
-30
-19
0.34
0.65
LHb_R
-0.23
-0.79
-0.46
0.59
0.62
0.62
2.31
2.33
-27
-19
0.41
0.59
AC_L CPu_R
1.07 1.3
1.01 2.08
1.18 1.68
0.65 0.42
0.58 0.64
0.63 0.58
2.73 2.47
2.49 2.71
-6 -19
-11 0
0.88 0.54
0.78 1
PC_L
-1.24
-1.18
-1.31
0.67
0.6
0.65
2.16
2.18
-35
-27
0.27
0.45
VM_L
-0.63
0.34
-0.32
0.63
0.59
0.6
2.22
2.37
-32
-18
0.35
0.67
PC_R
-1.19
-1.15
-1.26
0.64
0.51
0.59
2.18
2.2
-34
-26
0.27
0.45
VM_R
-0.61
0.28
-0.32
0.57
0.54
0.55
2.24
2.37
-31
-18
0.35
0.67
VTA_R
1.3
1
1.28
0.62
0.54
0.59
2.67
2.39
-9
-16
0.81
0.67
CPu_L
1.4
2.23
1.82
0.47
0.65
0.61
2.47
2.71
-19
0
0.54
1
PF_R
0.53
-0.08
0.34
0.54
0.55
0.55
2.43
2.31
-22
-21
0.54
0.64
VTA_L
0.72
-0.94
0.04
0.62
0.57
0.6
2.69
2.41
-8
-15
0.83
0.7
Pir_R
0.38
-0.94
-0.19
0.49
0.51
0.56
2.2
2.29
-32
-21
0.3
0.54
PF_L
0.59
-0.04
0.39
0.6
0.58
0.6
2.43
2.29
-22
-22
0.56
0.65
SNR_L
0.99
0.34
0.82
0.6
0.59
0.6
2.51
2.31
-18
-21
0.64
0.61
-0.74
0.56
0.54
0.56
2.16
2.2
-35
-26
0.26
0.48
CL_L SNR_R
-0.63
-0.8
0.91
0.28
0.75
0.57
0.55
0.56
2.51
2.33
-18
-20
0.63
0.6
CL_R Ent_L
-0.61 0.72
-0.79 0.52
-0.72 0.72
0.46 0.52
0.44 0.49
0.46 0.52
2.16 2.29
2.2 2.24
-34 -28
-26 -23
0.25 0.37
0.47 0.44
Pir_L
0.38
0.04
0.26
0.49
0.55
0.56
2.2
2.31
-32
-20
0.3
0.57
SNC_L
1.2
-0.42
0.68
0.59
0.6
0.6
2.53
2.27
-16
-22
0.67
0.54
SNC_R
1.11
-0.43
0.61
0.57
0.55
0.56
2.55
2.27
-15
-22
0.67
0.54
-1.02
0.52
-0.42
0.5
0.54
0.53
2.24
2.37
-29
-17
0.32
0.65
0.34
0.28
0.34
0.53
0.55
0.54
2.39
2.35
-23
-18
0.47
0.63
Ac_R STh_R STh_L
0.38
0.34
0.39
0.59
0.59
0.59
2.39
2.35
-23
-18
0.49
0.63
Ac_L
-0.32
0.52
0.04
0.55
0.58
0.57
2.27
2.37
-29
-17
0.34
0.65
MDL_L
-1.45
-0.04
-1.03
0.64
0.59
0.61
2.06
2.31
-40
-21
0.13
0.59
MDL_R
-1.38
-0.08
-0.99
0.56
0.54
0.54
2.08
2.31
-38
-21
0.13
0.58
VA_L
-1.45
-2.31
-1.88
0
0.38
0.28
1.69
1.92
-47
-37
0.06
0.14
0.52
0.44
0.5
2.43
2.2
-21
-25
0.53
0.47
0
0.38
0.28
1.67
1.92
-47
-37
0.06
0.14
MGP_L
0.99
0.34
0.82
-1.57
-2.59
-2.06
0.91
0.28
0.75
0.48
0.36
0.44
2.45
2.22
-20
-24
0.53
0.47
Ent_R MDM_R
-1.37 -1.72
-2.4 -0.94
-2.03 -1.57
0 0.44
0.44 0.58
0.28 0.56
1.86 1.92
1.71 2.14
-44 -42
-39 -27
0.06 0.08
0.1 0.43
MDM_L
-1.37
0.04
-0.88
0.38
0.54
0.51
1.94
2.16
-41
-26
0.1
0.46
LGP_L
-0.84
-0.42
-0.74
0
0.4
0.3
1.84
2.08
-42
-32
0.13
0.37
LGP_R
-0.99
-0.43
-0.86
0
0.3
0.22
1.82
2.1
-42
-31
0.12
0.36
VA_R MGP_R
The regions are sorted using the degree all. Abbreviations and parameter definitions are described in the Supporting Information A all, AD average degree, C circle, CC cluster-coefficient, DG degree, I in, L laterality, O out, SP length of shortest path
123
Brain Struct Funct Table 10 Frequency of motifs of the ipsilateral BG2 network Region
3-01
3-02
AC&1
0
1
27
0
CPu&1
5
2
20
2
VTA&10
0
1
40
2
MRF&23
0
0
44
2
STh&3
9
3
41
AGm&1
5
3
37
SNR&25
9
6
42
PF&8
3
6
45
SNC&26
2
0
43
MGP&12
17
2
42
5
AGl&3 Ac&5
2 14
2 7
26 34
0 7
LHb&4
10
1
47
5
12
2
13
21
8
1
27
26
2
5
VM&17
9
3
24
8
13
1
14
36
3
1
34
18
2
1
Ent&12
7
3
26
0
24
0
5
34
5
8
32
15
1
6
HIPP&12
8
6
28
8
29
0
12
22
6
8
29
11
2
2
CM&24
11
10
26
10
16
3
15
27
12
4
30
10
0
0
PC&22
4
4
30
0
11
0
3
47
4
1
15
9
3
2
LGP&18
9
4
52
6
6
2
12
11
8
3
25
6
3
5
13
7
27
15
15
2
28
11
14
11
32
6
3
2
MDL&28
9
3
36
11
8
0
22
28
3
2
25
5
0
3
MDM&26
8
1
38
1
11
0
13
22
1
7
13
5
1
3
CL&24
12
5
25
5
14
1
7
28
13
5
23
5
0
0
VL&31
7
2
43
4
7
0
22
23
1
4
15
3
0
1
VA&13
7
5
30
0
25
1
1
0
5
1
0
0
0
7
Pir&20
3-03
3-04
3-05
3-06
3-07
3-08
3-09
3-10
3-11
3-12
3-13
chain
5
0
12
30
0
9
53
1
38
12
0
65
83
0
0
3
73
68
0
1
0
28
34
0
8
3
66
65
0
0
2
0
40
17
8
3
54
65
0
0
4 2
8
0
22
0
16
29
8
6
38
48
0
0
6
63
5
11
36
46
3
2
16
7
3
32
16
20
1
33
45
0
3
5 10
9
0
17
27
6
2
34
45
0
0
14
0
34
10
10
9
39
44
0
0
11
1
21
21
10
3
38
42
5
8
30 15
1 1
0 8
57 24
8 19
3 5
41 41
39 32
0 0
0 0
in
out
Regions are sorted using motif 3-13 frequencies chain symmetricCHAIN, in symmetricIN, out symmetricOUT (see legend of Fig. 8)
Table 11 Global extrinsic connections of regions of the bilateral BG2 network Region
Dic
Dii
Dis
Doc
Doi
Dos
Centrolateral_thalamic_nucleus_L
392
446
838
7
68
75
403
470
873
7
103
Centrolateral_thalamic_nucleus_R
392
446
838
7
68
75
403
470
873
7
103
110
Accumbens_nucleus_R
265
380
645
4
74
78
308
853
1,161
11
309
320
Accumbens_nucleus_L
265
378
643
4
74
78
318
829
1,147
11
304
315
Central_medial_thalamic_nucleus_L
265
322
587
12
70
82
270
340
610
12
184
196
Central_medial_thalamic_nucleus_R
265
322
587
12
72
84
270
340
610
12
188
200
Paracentral_thalamic_nucleus_L
235
309
544
6
81
87
238
413
651
6
89
95
Paracentral_thalamic_nucleus_R
235
309
544
6
81
87
238
413
651
6
89
95
Mesencephalic_reticular_formation_L
222
311
533
227
334
561
3,316
4,273
7,589
1,188
3,510
4,698
Mesencephalic_reticular_formation_R
222
311
533
227
334
561
3,316
4,273
7,589
1,188
3,512
4,700
Amygdaloid_complex_R Amygdaloid_complex_L
144 144
378 375
522 519
4 4
74 71
78 75
632 633
4,659 4,473
5,291 5,106
525 528
4,296 3,970
4,821 4,498
Ventrolateral_thalamic_nucleus_L
144
344
488
8
52
60
160
384
544
9
56
65
Ventrolateral_thalamic_nucleus_R
144
344
488
8
52
60
160
384
544
9
56
65
Caudate_putamen_L
54
311
365
21
79
100
279
1,057
1,336
22
172
194
Caudate_putamen_R
54
311
365
21
79
100
273
1,064
1,337
22
172
194
123
Sic
Sii
Sis
Soc
Soi
Sos 110
Brain Struct Funct Table 11 continued Region
Dic
Dii
Dis
Doc
Doi
Dos
Sic
Sii
Sis
Soc
Soi
Sos
Ventromedial_thalamic_nucleus_L
82
244
326
33
54
87
105
356
461
40
161
201
Ventromedial_thalamic_nucleus_R
82
244
326
33
54
87
105
356
461
40
161
201
Parafascicular_thalamic_nucleus_L Parafascicular_thalamic_nucleus_R
111 111
213 213
324 324
35 35
136 137
171 172
140 140
377 377
517 517
112 112
291 292
403 404
Ventral_tegmental_area_A10_L
60
238
298
99
231
330
96
476
572
124
440
564
Ventral_tegmental_area_A10_R
60
238
298
99
235
334
96
476
572
124
446
570
Medial_agranular_prefrontal_cortex_R
50
138
188
106
418
524
85
219
304
159
586
745
Hippocampus_R
20
162
182
8
47
55
138
1,220
1,358
130
649
779
Medial_agranular_prefrontal_cortex_L
50
131
181
106
416
522
85
212
297
159
584
743
Entorhinal_cortex_L
27
137
164
14
63
77
154
900
1,054
21
546
567
Hippocampus_L
25
137
162
8
46
54
154
1,030
1,184
130
631
761
Lateral_habenular_nucleus_L
24
128
152
37
99
136
55
226
281
40
129
169
Lateral_habenular_nucleus_R
24
128
152
37
99
136
55
226
281
40
129
169
Subthalamic_nucleus_L
19
124
143
23
83
106
29
181
210
34
127
161
Subthalamic_nucleus_R
19
124
143
23
83
106
29
181
210
34
129
163
Substantia_nigra_compact_part_L
24
101
125
30
126
156
37
136
173
33
181
214
Substantia_nigra_compact_part_R
24
101
125
30
126
156
37
136
173
33
181
214
Piriform_cortex_R Lateral_agranular_prefrontal_cortex_R
11 22
112 99
123 121
11 257
55 383
66 640
34 71
309 440
343 511
37 504
227 803
264 1,307
Lateral_agranular_prefrontal_cortex_L
22
95
117
257
381
638
71
436
507
504
801
1,305
Piriform_cortex_L
11
103
114
11
53
64
34
300
334
38
191
229
Mediodorsal_thalamic_nucleus_medial_part_L
11
91
102
0
14
14
11
91
102
0
14
14
Mediodorsal_thalamic_nucleus_medial_part_R
11
91
102
0
16
16
11
91
102
0
16
16
Mediodorsal_thalamic_nucleus_lateral_part_L
17
82
99
0
18
18
17
83
100
0
18
18
Mediodorsal_thalamic_nucleus_lateral_part_R
17
82
99
0
18
18
17
83
100
0
18
18
Substantia_nigra_reticular_part_L
14
73
87
34
133
167
14
86
100
45
200
245
Substantia_nigra_reticular_part_R
14
73
87
34
133
167
14
86
100
45
200
245
Medial_globus_pallidus_L
12
53
65
11
59
70
13
58
71
11
73
84
Medial_globus_pallidus_R
12
53
65
11
59
70
13
58
71
11
73
84
Entorhinal_cortex_R
3
40
43
0
5
5
9
68
77
0
22
22
Lateral_globus_pallidus_L
4
22
26
0
4
4
4
22
26
0
5
5
Lateral_globus_pallidus_R
4
22
26
0
4
4
4
22
26
0
5
5
12 12
12 12
24 24
0 0
10 10
10 10
14 14
14 14
28 28
0 0
10 10
10 10
Ventro_anterior_thalamic_nucleus_L Ventro_anterior_thalamic_nucleus_R
The regions were sorted by the direct input sum (Dis) Dic direct input from contralateral, Dii direct input from ipsilateral, Dis direct input from ipsi- and contralateral, Doc direct output to contralateral, Doi direct output to ipsilateral, Dos direct output to ipsi- and contralateral, Sic subtree input from contralateral, Sii subtree input from ipsilateral, Sis subtree input from ipsi- and contralateral, Soc subtree output to contralateral, Soi subtree output to ipsilateral, Sos subtree output to ipsi- and contralateral
123
123 0.971
0.829 -0.558 0.957
-0.062
-0.090
-0.006
0.181
-0.009
0.030
-0.043
-0.152
0.002 0.147
-0.090
0.080
DGA
DGO
DGI
CDC
LA
LO
LI
LR
Katz LC
Triag
CyclC
0.015 0.024
CCI BC
0.114
-0.017
-0.012
-0.139
-0.139
PRC
FC
Stress
ZO
ZI
-0.172
-0.110
CCO
-0.084
-0.127
Loc
EC
-0.059
Lev
SC
0.106
VCDG
-0.138
0.285
-0.018
CCT
-0.093
-0.012
CCA
CC2
-0.086
CCI
ADnb
-0.833
0.094
CCO
0.822
0.822
0.944
0.814
0.834
0.969
0.915
0.801 0.904
0.922
0.907
0.926
0.686
-0.760
-0.827
-0.445
-0.808
-0.146
EccI
-0.672
0.019
-0.289
EccO
-0.741
0.047
-0.433
-0.487
-0.573
-0.624
0.855
0.930
1.000
-0.018
1.000
-0.018
REC
HD
HD
REC
Region
0.863
0.863
0.941
0.697
0.811
0.978
0.936
0.833 0.878
0.953
0.973
0.969
0.719
-0.780
-0.219
-0.863
-0.867
-0.402
-0.859
-0.166
-0.631
-0.810
0.981
0.860 -0.589
0.076
-0.486
-0.511
-0.615
-0.636
0.884
0.955
1.000
0.971
-0.062
DGA
0.896
0.896
0.881
0.693
0.629
0.950
0.988
0.636 0.786
0.989
0.935
0.926
0.709
-0.787
-0.199
-0.824
-0.862
-0.320
-0.867
0.031
-0.616
-0.748
0.957
0.678 -0.547
0.036
-0.375
-0.529
-0.560
-0.803
0.705
1.000
0.955
0.930
-0.090
DGO
0.648
0.648
0.858
0.573
0.945
0.839
0.679
0.987 0.860
0.718
0.851
0.856
0.600
-0.622
-0.209
-0.763
-0.711
-0.456
-0.685
-0.446
-0.535
-0.756
0.834
0.987 -0.545
0.125
-0.570
-0.386
-0.586
-0.253
1.000
0.705
0.884
0.855
-0.006
DGI
Table 12 Correlation coefficients of local parameters of the bilateral BG2 network
-0.807
-0.807
-0.532
-0.657
-0.217
-0.620
-0.821
-0.179 -0.418
-0.767
-0.616
-0.677
-0.625
0.673
0.009
0.622
0.729
0.107
0.776
-0.278
0.476
0.555
-0.627
-0.221 0.442
-0.048
0.148
0.477
0.377
1.000
-0.253
-0.803
-0.636
-0.624
0.181
CDC
-0.423
-0.423
-0.577
-0.450
-0.548
-0.593
-0.545
-0.603 -0.540
-0.616
-0.536
-0.576
-0.305
0.285
0.503
0.735
0.653
0.491
0.649
0.416
0.493
0.779
-0.565
-0.612 0.412
0.281
0.830
0.814
1.000
0.377
-0.586
-0.560
-0.615
-0.573
-0.009
LA
-0.448
-0.448
-0.452
-0.469
-0.388
-0.480
-0.520
-0.356 -0.413
-0.587
-0.460
-0.505
-0.373
0.340
0.339
0.603
0.565
0.213
0.634
0.156
0.433
0.700
-0.464
-0.405 0.589
0.052
0.373
1.000
0.814
0.477
-0.386
-0.529
-0.511
-0.487
0.030
LO
-0.246
-0.246
-0.446
-0.280
-0.493
-0.464
-0.365
-0.630 -0.417
-0.408
-0.412
-0.447
-0.121
0.114
0.506
0.576
0.488
0.523
0.440
0.526
0.330
0.578
-0.436
-0.606 0.106
0.365
1.000
0.373
0.830
0.148
-0.570
-0.375
-0.486
-0.433
-0.043
LI
0.118
0.118
-0.044
0.094
0.038
-0.004
0.071
0.076 -0.054
0.052
0.151
0.182
0.347
-0.195
0.131
0.193
0.078
0.652
-0.103
0.110
0.064
-0.106
0.033
0.117 -0.404
1.000
0.365
0.052
0.281
-0.048
0.125
0.036
0.076
0.047
-0.152
LR
0.598
0.598
0.812
0.559
0.920
0.816
0.655
0.988 0.807
0.703
0.823
0.831
0.529
-0.541
-0.292
-0.721
-0.667
-0.389
-0.675
-0.494
-0.511
-0.747
0.808
1.000 -0.496
0.117
-0.606
-0.405
-0.612
-0.221
0.987
0.678
0.860
0.829
0.002
Katz
Brain Struct Funct
-0.211
-0.043
-0.194
-0.135
-0.019
-0.103
-0.016
-0.096
LC
PCI
PCA
RadO
RadI
CenO
CenI
Hub
Aut
Region
0.589
0.106
-0.404
-0.496
1.000 -0.537
0.680
0.401
0.071
0.595
0.114
0.596
0.564
0.026
LO
LI
LR
Katz
LC Triag
CyclC
EccO
EccI
CCO
CCI
CCA
CCT
CC2
-0.682
0.412
LA
VCDG
0.442
CDC
0.623
0.808
-0.545
DGI
ADnb
0.033
-0.547
DGO
0.957
0.642
-0.747
-0.273
-0.810
-0.822
-0.401
-0.786
-0.107
-0.617
-0.725
-0.537 1.000
-0.436
-0.464
-0.565
-0.627
0.834
0.957
0.981
-0.558
-0.589
REC
-0.090
Triag
0.925
0.820
0.934
0.790
0.912
0.521
0.481
0.517
DGA
0.147
-0.109
HD
0.747
-0.013
ZA
PCO 0.861
REC
HD
Region
Table 12 continued
-0.722
0.647
0.228
0.884
0.853
0.341
0.906
0.284
0.514
1.000
0.680 -0.725
-0.747
-0.106
0.578
0.700
0.779
0.555
-0.756
-0.748
-0.810
-0.741
0.080
CyclC
0.949
0.857
0.962
0.833
0.949
0.601
0.545
0.584
0.897
0.765
DGA
-0.446
0.422
0.216
0.580
0.534
0.382
0.538
0.210
1.000
0.514
0.401 -0.617
-0.511
0.064
0.330
0.433
0.493
0.476
-0.535
-0.616
-0.631
-0.672
0.019
EccO
0.997
0.673
0.995
0.648
0.967
0.518
0.409
0.554
0.863
0.618
DGO
-0.074
-0.010
0.204
0.224
0.160
0.286
0.051
1.000
0.210
0.284
0.071 -0.107
-0.494
0.110
0.526
0.156
0.416
-0.278
-0.446
0.031
-0.166
-0.146
-0.289
EccI
0.695
0.987
0.731
0.969
0.743
0.619
0.657
0.521
0.783
0.855
DGI
-0.779
0.757
0.105
0.877
0.898
0.260
1.000
0.051
0.538
0.906
0.595 -0.786
-0.675
-0.103
0.440
0.634
0.649
0.776
-0.685
-0.867
-0.859
-0.808
0.094
CCO
-0.813
-0.231
-0.795
-0.200
-0.788
-0.290
-0.094
-0.433
-0.666
-0.278
CDC
-0.208
0.345
-0.057
0.658
0.542
1.000
0.260
0.286
0.382
0.341
0.114 -0.401
-0.389
0.652
0.523
0.213
0.491
0.107
-0.456
-0.320
-0.402
-0.445
-0.086
CCI
-0.590
-0.656
-0.594
-0.619
-0.642
-0.747
-0.684
-0.693
-0.449
-0.402
LA
-0.796
0.856
-0.008
0.966
1.000
0.542
0.898
0.160
0.534
0.853
0.596 -0.822
-0.667
0.078
0.488
0.565
0.653
0.729
-0.711
-0.862
-0.867
-0.827
-0.012
CCA
-0.559
-0.438
-0.562
-0.389
-0.610
-0.586
-0.429
-0.605
-0.419
-0.299
LO
-0.721
0.769
0.057
1.000
0.966
0.658
0.877
0.224
0.580
0.884
0.564 -0.810
-0.721
0.193
0.576
0.603
0.735
0.622
-0.763
-0.824
-0.863
-0.833
-0.018
CCT
-0.398
-0.637
-0.400
-0.639
-0.446
-0.666
-0.722
-0.565
-0.314
-0.352
LI
0.293
-0.297
1.000
0.057
-0.008
-0.057
0.105
0.204
0.216
0.228
0.026 -0.273
-0.292
0.131
0.506
0.339
0.503
0.009
-0.209
-0.199
-0.219
-0.138
0.285
CC2
0.021
0.109
0.076
0.108
0.134
0.134
0.101
0.167
0.198
0.309
LR
-0.930
1.000
-0.297
0.769
0.856
0.345
0.757
-0.010
0.422
0.647
0.623 -0.747
-0.541
-0.195
0.114
0.340
0.285
0.673
-0.622
-0.787
-0.780
-0.760
-0.093
ADnb
0.674
0.987
0.708
0.986
0.735
0.673
0.700
0.581
0.735
0.822
Katz
Brain Struct Funct
123
123
0.709
0.600
-0.625
-0.305
-0.373
-0.121
0.347
0.529
-0.682
CDC
LA
LO
LI
LR
Katz
LC
VCDG
Region
DGI
-0.538
Aut
DGO
-0.569
Hub
0.719
-0.590
CenI
DGA
-0.487
CenO
0.686
-0.605
RadI
0.106
-0.426
RadO
REC
-0.397 -0.373
PCI PCA
HD
-0.618
-0.480
FC
PCO
-0.527
PRC
-0.537
-0.516
SC
-0.598
-0.570
EC
ZA
-0.489
BC
ZI
-0.486
CCI
-0.511
-0.556
CCO
-0.598
-0.622
Loc
ZO
-0.645
Lev
Stress
LC
Region
Table 12 continued
-0.645
0.831
0.182
-0.447
-0.505
-0.576
-0.677
0.856
0.926
0.969
0.926
-0.059
Lev
0.947
0.797
0.953
0.767
0.919
0.531
0.502 0.504
0.840
0.682
0.828
0.828
0.940
0.652
0.742
0.996
0.942
0.864
0.779
0.959
0.968
0.913
Triag
-0.622
0.823
0.151
-0.412
-0.460
-0.536
-0.616
0.851
0.935
0.973
0.907
-0.127
Loc
-0.757
-0.794
-0.786
-0.749
-0.839
-0.762
-0.732 -0.654
-0.783
-0.679
-0.753
-0.753
-0.723
-0.558
-0.711
-0.720
-0.722
-0.679
-0.728
-0.771
-0.774
-0.849
CyclC
-0.556
0.703
0.052
-0.408
-0.587
-0.616
-0.767
0.718
0.989
0.953
0.922
-0.110
CCO
-0.630
-0.536
-0.648
-0.494
-0.633
-0.303
-0.359 -0.161
-0.615
-0.544
-0.583
-0.583
-0.658
-0.604
-0.572
-0.634
-0.619
-0.639
-0.516
-0.634
-0.564
-0.588
EccO
-0.486
0.988
0.076
-0.630
-0.356
-0.603
-0.179
0.987
0.636
0.833
0.801
0.015
CCI
0.013
-0.493
0.001
-0.480
-0.034
-0.265
-0.136 -0.392
-0.088
-0.322
0.046
0.046
-0.204
-0.062
-0.470
-0.125
0.067
-0.267
-0.519
0.015
-0.107
-0.153
EccI
-0.489
0.807
-0.054
-0.417
-0.413
-0.540
-0.418
0.860
0.786
0.878
0.904
0.024
BC
-0.872
-0.693
-0.890
-0.673
-0.932
-0.669
-0.731 -0.493
-0.840
-0.632
-0.860
-0.860
-0.753
-0.684
-0.644
-0.783
-0.849
-0.682
-0.641
-0.869
-0.823
-0.912
CCO
-0.570
0.655
0.071
-0.365
-0.520
-0.545
-0.821
0.679
0.988
0.936
0.915
-0.172
EC
-0.313
-0.442
-0.300
-0.368
-0.242
-0.081
0.010 -0.197
-0.307
-0.299
-0.272
-0.272
-0.557
-0.330
-0.563
-0.440
-0.280
-0.616
-0.456
-0.289
-0.300
-0.310
CCI
-0.516
0.816
-0.004
-0.464
-0.480
-0.593
-0.620
0.839
0.950
0.978
0.969
-0.084
SC
-0.855
-0.701
-0.865
-0.639
-0.862
-0.496
-0.506 -0.448
-0.817
-0.613
-0.837
-0.837
-0.830
-0.653
-0.692
-0.820
-0.830
-0.776
-0.659
-0.839
-0.818
-0.882
CCA
-0.527
0.920
0.038
-0.493
-0.388
-0.548
-0.217
0.945
0.629
0.811
0.834
0.114
PRC
-0.822
-0.763
-0.832
-0.703
-0.832
-0.547
-0.530 -0.490
-0.803
-0.644
-0.799
-0.799
-0.858
-0.634
-0.774
-0.818
-0.783
-0.830
-0.726
-0.813
-0.797
-0.857
CCT
-0.480
0.559
0.094
-0.280
-0.469
-0.450
-0.657
0.573
0.693
0.697
0.814
-0.017
FC
-0.228
-0.279
-0.226
-0.308
-0.267
-0.549
-0.483 -0.495
-0.043
-0.041
-0.039
-0.039
-0.102
0.081
-0.005
-0.274
-0.248
0.005
-0.254
-0.297
-0.267
-0.121
CC2
-0.511
0.812
-0.044
-0.446
-0.452
-0.577
-0.532
0.858
0.881
0.941
0.944
-0.012
Stress
-0.759
-0.555
-0.774
-0.488
-0.731
-0.178
-0.203 -0.180
-0.800
-0.617
-0.814
-0.814
-0.765
-0.621
-0.626
-0.731
-0.745
-0.727
-0.525
-0.731
-0.777
-0.824
ADnb
Brain Struct Funct
-0.121
-0.824
-0.446
-0.074
-0.779
-0.208
-0.796
-0.721
0.293
-0.930
1.000 0.816
0.713
0.665
0.509
0.649
0.673
0.623
0.610
0.580
0.675
0.768
0.768
0.680
0.794
0.346
0.244 0.307
0.714
0.500
0.721
0.563
0.692
EccO
EccI
CCO
CCI
CCA
CCT
CC2
ADnb
VCDG Lev
Loc
CCO
CCI
BC
EC
SC
PRC
FC
Stress
ZO
ZI
ZA
PCO
PCI
PCA RadO
RadI
CenO
CenI
Hub
Aut
0.913
0.923
0.837
0.944
0.823
0.952
0.517 0.618
0.635
0.933
0.804
0.892
0.892
0.866
0.700
0.786
0.903
0.913
0.804
0.801
0.915
0.949
0.816 1.000
-0.857
-0.882
-0.310
-0.912
-0.153
-0.588
-0.849
0.642
-0.722
Lev
CyclC
VCDG
Triag
Region
Table 12 continued
0.923
0.821
0.941
0.801
0.923
0.514 0.572
0.545
0.885
0.735
0.862
0.862
0.890
0.570
0.733
0.949
0.922
0.809
0.795
0.936
1.000
0.713 0.949
-0.777
-0.267
-0.797
-0.818
-0.300
-0.823
-0.107
-0.564
-0.774
0.968
Loc
0.989
0.699
0.992
0.675
0.977
0.467 0.590
0.624
0.839
0.621
0.860
0.860
0.874
0.671
0.627
0.952
0.980
0.774
0.656
1.000
0.936
0.665 0.915
-0.731
-0.297
-0.813
-0.839
-0.289
-0.869
0.015
-0.634
-0.771
0.959
CCO
0.629
0.989
0.667
0.982
0.689
0.685 0.639
0.532
0.725
0.821
0.584
0.584
0.809
0.524
0.935
0.786
0.610
0.817
1.000
0.656
0.795
0.509 0.801
-0.525
-0.254
-0.726
-0.659
-0.456
-0.641
-0.519
-0.516
-0.728
0.779
CCI
0.770
0.812
0.783
0.757
0.741
0.412 0.395
0.330
0.777
0.736
0.708
0.708
0.978
0.668
0.901
0.880
0.734
1.000
0.817
0.774
0.809
0.649 0.804
-0.727
0.005
-0.830
-0.776
-0.616
-0.682
-0.267
-0.639
-0.679
0.864
BC
0.989
0.654
0.989
0.632
0.968
0.413 0.541
0.589
0.848
0.600
0.884
0.884
0.837
0.705
0.586
0.935
1.000
0.734
0.610
0.980
0.922
0.673 0.913
-0.745
-0.248
-0.783
-0.830
-0.280
-0.849
0.067
-0.619
-0.722
0.942
EC
0.941
0.803
0.946
0.771
0.913
0.511 0.534
0.506
0.823
0.676
0.808
0.808
0.948
0.695
0.762
1.000
0.935
0.880
0.786
0.952
0.949
0.623 0.903
-0.731
-0.274
-0.818
-0.820
-0.440
-0.783
-0.125
-0.634
-0.720
0.996
SC
0.622
0.927
0.654
0.895
0.655
0.525 0.486
0.403
0.747
0.837
0.610
0.610
0.850
0.652
1.000
0.762
0.586
0.901
0.935
0.627
0.733
0.610 0.786
-0.626
-0.005
-0.774
-0.692
-0.563
-0.644
-0.470
-0.572
-0.711
0.742
PRC
0.701
0.548
0.706
0.510
0.714
0.281 0.336
0.407
0.655
0.547
0.645
0.645
0.673
1.000
0.652
0.695
0.705
0.668
0.524
0.671
0.570
0.580 0.700
-0.621
0.081
-0.634
-0.653
-0.330
-0.684
-0.062
-0.604
-0.558
0.652
FC
0.867
0.813
0.877
0.763
0.835
0.449 0.451
0.402
0.826
0.725
0.784
0.784
1.000
0.673
0.850
0.948
0.837
0.978
0.809
0.874
0.890
0.675 0.866
-0.765
-0.102
-0.858
-0.830
-0.557
-0.753
-0.204
-0.658
-0.723
0.940
Stress
Brain Struct Funct
123
123
-0.583
0.046
-0.860
-0.272
-0.837
-0.799
-0.039
-0.814
0.768
0.892
0.862
0.860
0.584 0.708
0.884
0.808
0.610
0.645
0.784
0.000
1.000
EccO
EccI
CCO
CCI
CCA
CCT
CC2
ADnb
VCDG
Lev
Loc
CCO
CCI BC
EC
SC
PRC
FC
Stress
ZO
ZI
0.598 -0.598
0.828
0.118
LR
Katz LC
-0.753
-0.246
LI
CyclC
0.309
-0.448
Triag
-0.352
-0.423
LO
0.671
1.000
0.725
0.547
0.837
0.676
0.600
0.821 0.736
0.621
0.735
0.804
0.680
-0.617
-0.041
-0.644
-0.613
-0.299
-0.632
-0.322
-0.544
-0.679
0.682
0.822 -0.537
-0.299
-0.402
-0.278
0.855
0.618
LA
0.896
DGO
0.765
0.648
0.863
DGA
0.747
-0.013
-0.807
0.822
REC
CDC
-0.139
HD
ZI
DGI
ZO
Region
Table 12 continued
0.953
0.671
0.826
0.655
0.747
0.823
0.848
0.725 0.777
0.839
0.885
0.933
0.794
-0.800
-0.043
-0.803
-0.817
-0.307
-0.840
-0.088
-0.615
-0.783
0.840
0.735 -0.618
0.198
-0.314
-0.419
-0.449
-0.666
0.783
0.863
0.897
0.861
-0.109
ZA
0.445
0.953
0.402
0.407
0.403
0.506
0.589
0.532 0.330
0.624
0.545
0.635
0.346
-0.203
-0.483
-0.530
-0.506
0.010
-0.731
-0.136
-0.359
-0.732
0.502
0.581 -0.397
0.167
-0.565
-0.605
-0.693
-0.433
0.521
0.554
0.584
0.517
-0.211
PCO
0.276
0.445
0.449
0.281
0.525
0.511
0.413
0.685 0.412
0.467
0.514
0.517
0.244
-0.180
-0.495
-0.490
-0.448
-0.197
-0.493
-0.392
-0.161
-0.654
0.504
0.700 -0.373
0.101
-0.722
-0.429
-0.684
-0.094
0.657
0.409
0.545
0.481
-0.043
PCI
0.406
0.276
0.451
0.336
0.486
0.534
0.541
0.639 0.395
0.590
0.572
0.618
0.307
-0.178
-0.549
-0.547
-0.496
-0.081
-0.669
-0.265
-0.303
-0.762
0.531
0.673 -0.426
0.134
-0.666
-0.586
-0.747
-0.290
0.619
0.518
0.601
0.521
-0.194
PCA
0.883
0.406
0.835
0.714
0.655
0.913
0.968
0.689 0.741
0.977
0.923
0.952
0.714
-0.731
-0.267
-0.832
-0.862
-0.242
-0.932
-0.034
-0.633
-0.839
0.919
0.735 -0.605
0.134
-0.446
-0.610
-0.642
-0.788
0.743
0.967
0.949
0.912
-0.135
RadO
0.584
0.883
0.763
0.510
0.895
0.771
0.632
0.982 0.757
0.675
0.801
0.823
0.500
-0.488
-0.308
-0.703
-0.639
-0.368
-0.673
-0.480
-0.494
-0.749
0.767
0.986 -0.487
0.108
-0.639
-0.389
-0.619
-0.200
0.969
0.648
0.833
0.790
-0.019
RadI
0.899
0.584
0.877
0.706
0.654
0.946
0.989
0.667 0.783
0.992
0.941
0.944
0.721
-0.774
-0.226
-0.832
-0.865
-0.300
-0.890
0.001
-0.648
-0.786
0.953
0.708 -0.590
0.076
-0.400
-0.562
-0.594
-0.795
0.731
0.995
0.962
0.934
-0.103
CenO
0.610
0.899
0.813
0.548
0.927
0.803
0.654
0.989 0.812
0.699
0.821
0.837
0.563
-0.555
-0.279
-0.763
-0.701
-0.442
-0.693
-0.493
-0.536
-0.794
0.797
0.987 -0.569
0.109
-0.637
-0.438
-0.656
-0.231
0.987
0.673
0.857
0.820
-0.016
CenI
0.894
0.610
0.867
0.701
0.622
0.941
0.989
0.629 0.770
0.989
0.923
0.923
0.692
-0.759
-0.228
-0.822
-0.855
-0.313
-0.872
0.013
-0.630
-0.757
0.947
0.674 -0.538
0.021
-0.398
-0.559
-0.590
-0.813
0.695
0.997
0.949
0.925
-0.096
Hub
0.560
0.894
0.736
0.441
0.803
0.785
0.646
0.950 0.700
0.690
0.822
0.798
0.459
-0.460
-0.430
-0.643
-0.610
-0.268
-0.642
-0.460
-0.434
-0.727
0.788
0.969 -0.462
0.179
-0.630
-0.368
-0.597
-0.204
0.944
0.653
0.826
0.759
-0.094
Aut
Brain Struct Funct
0.954
0.649
0.669
1.000
Cluster coefficient (triangle based): The triangle based cluster coefficient (Fagiolo 2007) of a node n is the number of triangles around n divided by the maximum possible number. In this version of the cluster coefficient reciprocal edges to a neighbor of a node n can affect the cluster coefficient of node n. In the other version only edges between neighbors of n have an influence to the cluster coefficient of node n. t! ðiÞ tmax ðiÞ ! ! t w ðiÞ w CT ¼ tmax ðiÞ
0.710
0.995
0.984
0.648
0.669
0.685 0.995 1.000 0.686
1.000
0.964 0.648 0.686 1.000
0.710
0.724 0.972 0.985 0.723
0.984
0.736 0.545 0.572 0.715
0.739
0.760 0.421 0.445 0.706
0.699
0.635 0.586 0.612 0.635
0.720
0.688 0.859
0.603
0.765 0.613
0.879 0.730
0.749
0.655 0.825
0.836
Aut CenO RadI
CenI
Hub
Brain Struct Funct
CT! ¼
ð55Þ ð56Þ
Cluster coefficient: Number of edges between the neighbors of a node divided by the maximum possible number. C ! ðiÞ refers to all neighbors of i. X 1 ajk C ! ðiÞ ¼ ð57Þ jNi j ðjNi j 1Þ j;k2N i
! ðiÞ refers to the out-neighbors of i. Cout X 1 ! Cout ðiÞ ¼ out ajk out jNi j ðjNi j 1Þ j;k2N out
0.972
0.739
0.985
0.723
1.000
0.673
0.519
0.724
0.876
0.677
RadO
j6¼k
ð58Þ
i
0.545
! Cin ðiÞ refers to the in-neighbors of i. X 1 ! Cin ðiÞ ¼ in ajk jNi j ðjNiin j 1Þ in j;k2Ni j6¼k
ð59Þ
In the weighted case the aij are replaced by the wij . Average cluster coefficient: n 1X C! ¼ C! ð60Þ n i¼1 i Most local parameters show high correlations with DGA
0.421 Aut
0.894
0.613
0.859
0.586
0.699 0.720 Hub
0.610
0.836
0.749
0.603
0.572 0.445 CenI
0.899
0.655
0.879
0.612
0.715 0.706 CenO
0.584
0.825
0.730
0.635
0.673 0.519
1.000 0.880
0.876
0.724
0.496
0.677 0.883
0.406 RadO
RadI
0.479
0.942
0.880 1.000 PCA
0.276
0.421
0.355
0.704
0.942 0.704 PCI
0.445
0.491
0.505
1.000
0.479 0.355
0.496 0.421
0.505 1.000
0.491 1.000
0.862 0.953
0.671 ZA
PCO
0.862
PCI ZO Region
Table 12 continued
ZI
ZA
PCO
PCA
j6¼k
and n ! 1X ! Cw ¼ Ciw n i¼1
ð61Þ
Small worldness S: C Crand
S¼
d
ð62Þ
drand
Centrality: Pn degmax degðiÞ n degmax 2 ‘ ¼ CD ¼ i¼1 ðn 1Þ ðn 2Þ ðn 1Þ ðn 2Þ
ð63Þ
This centrality (degree centrality) is defined for an undirected network based on undirected degrees. A directed or
123
Brain Struct Funct Table 13 The local parameters of the bilateral BG2 network were applied to the PCA Component
DGAll
CluCAll
CluC2
AvgDGnb
VCDG
Loc
Share (%)
1
0.440
-0.439
0.154
-0.457
0.438
0.434
74.978
2
-0.280
0.101
0.860
-0.166
0.216
-0.312
19.582
3
-0.040
0.818
-0.148
-0.182
0.434
0.292
2.705
4
-0.425
-0.296
-0.417
-0.049
0.626
-0.404
1.875
5
-0.521
-0.056
-0.133
-0.713
-0.408
0.180
0.621
6
-0.523
-0.191
0.151
0.468
0.117
0.659
0.24
The first two components determine the axes of the PCA-plane Share The share of variance of the data represented by the component
weighted version is not available yet. For the calculation the directed network is transferred to an undirected one. Circle length LC: dði; iÞ; dði; iÞ\1 LCðiÞ ¼ ð64Þ 0; dði; iÞ ¼ 1
! 1 X degwall ðjÞ degNB w ðiÞ ¼ jNi j j2Ni
ð70Þ
Variation coefficient of neighbor degree: qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi P 2 1 ðdeg ðjÞ degNBðiÞÞ all j2N jNi j i
ð71Þ
Eccentricity out: Eccentricity out, the output eccentricity of the vertex i is the maximum distance from i to any vertex.
VCðiÞ ¼
Eccout ðiÞ ¼ maxfdði; jÞjj 2 Ng
The weighted case is analogue. Locality index of node i ðLocðiÞÞ: The locality index of node i is the fraction of edges adjacent to nodes in Niþ whose source and target lie in Niþ . P P k2Niþ ajk j2N þ
ð65Þ
Eccentricity in: Eccentricity in, the input eccentricity of the vertex i is the maximum distance from i to any vertex. Eccin ðiÞ ¼ maxfdðj; iÞjj 2 Ng
ð66Þ
Cluster-coefficient of second neighbors: The cluster-coefficient of second neighbors (Hierarchical directed cluster coefficient of second (indirect) neighbors) C2 ðiÞ is the number of edges between the 2nd neighbors of node i, divided by the maximum possible number of edges. In the weighted case it is the sum of weights of the edges between the 2nd neighbors of node i, divided by the maximum possible sum. With ![ N2 ðiÞ ¼ Nj Niþ ; ð67Þ
degNBðiÞ
i
LocðiÞ ¼ P
j2Niþ
CCOUT ðiÞ ¼ P ð68Þ
In the weighted case the aij are replaced by wij . Average neighbor degree: The non-weighted average neighbor degree NBðiÞ of node i is 1 X degNBðiÞ ¼ degall ðjÞ ð69Þ j2Ni jNi j Weighted average neighbor degree: The weighted average neighbor degree NBðiÞ of node i is
123
k2N k6¼j
ð72Þ
ajk
The weighted case is analogue. A value of 0 means that the node is isolated. The larger the value, the less edges connect the neighborhood of i to outside node. The maximum of one is reached if the neighborhood of i is not connected to outside nodes. Closeness centrality out CCout ðiÞ: The closeness centrality out with indices of nodes from which node i can be reached (RNOUT ðiÞ)RNOUT ðiÞ ¼ fj 2 Nnfigjdði; jÞ\1g
j2Ni
the set of second neighbors of node i is: 8 1 < ; if jN2 ðiÞj [ 1 C2 ðiÞ ¼ jN2 ðiÞj ðjN2 ðiÞj 1Þ : 0; otherwise
P k6¼j
jRNOUT ðiÞj j2RNOUT ðiÞ dði; jÞ
ð73Þ
Closeness centrality in CCin ðiÞ: The closeness centrality in with indices of nodes which can be reached from node i (RNIN ðiÞ)RNIN ðiÞ ¼ fj 2 Nnfigjdðj; iÞ\1g CCIN ðiÞ ¼ P
jRNIN ðiÞj j2RNIN ðiÞ dðj; iÞ
ð74Þ
Betweenness centrality BC: BCðiÞ ¼
1 ðn 1Þðn 2Þ
X j;k2Nnfig
qj;k ðiÞ qj;k
ð75Þ
0.60000
1.00000
0.69500
0.71600
0.62300
MGP
CPu
STh
AGl
AGm
PRC
0.41700
0.04282
AGm
LGP
0.06481
AGl
0.62400
0.06366
STh
0.35400
0.17824
CPu
VL
0.04282
MGP
VA
0.01736 0.00000
VL LGP
0.51000 0.59000
PRC
0.00694
VA
SNC SNR
0.02778
0.03472
SNR
0.35100
0.44400
0.59000
0.58900
0.28800
0.20600
0.34100
0.11000
0.46800 0.32300
0.03206
0.02035
0.02062
0.00850
0.00697 0.02064
0.00114
0.00011
0.00002
0.02083
SNC
Direct BC
BG1 BC
Region
0.17400
0.21900
0.27900
0.29300
0.14200
0.10300
0.17000
0.05400
0.23000 0.16000
PRC
0.00290
0.00610
0.00590
0.00880
0.00200
0.00030 0.00002
0.00001
0.00510
0.00710
Indirect BC
VA
LGP
SNC
SNR
MGP
AGm
VL
STh
CPu AGl
PRC
LGP
VA
VL
SNC
SNR
MGP AGm
STh
AGl
CPu
BG1 BC
Table 14 The local parameters of the embedded BG1 network
VA
LGP
MGP
SNR
VL
AGm
AGl
SNC
STh CPu
PRC
SNC
SNR
VA
VL
MGP
CPu STh
LGP
AGm
AGl
Direct BC
VA
LGP
MGP
SNR
VL
AGm
AGl
SNC
CPu STh
PRC
VA
LGP
VL
MGP
AGm
STh SNR
AGl
SNC
CPu
Indirect BC
0.39200
0.47700
0.79600
1.00000
0.63900
0.42600
0.31000
0.11800
0.59200 0.55700
SC
0.64900
0.86500
0.92800
0.98200
0.94400
0.39400 0.72100
0.29200
0.71700
1.00000
BG1 EC
0.54100
0.65600
0.82500
0.50400
0.46700
0.14000
0.33100
0.02600
1.00000 0.70500
SC
0.43700
0.40000
0.39900
0.48000
0.48200
0.23700 0.16300
0.06400
0.57900
0.68000
Direct EC
0.55020
0.46600
0.64800
0.81660
0.50320
0.13510
0.33170
0.02520
1.00000 0.70280
SC
0.39000
0.35300
0.35300
0.42500
0.42800
0.21100 0.13900
0.05500
0.51400
0.60200
Indirect EC
VA
VL
AGm
LGP
AGl
SNR
SNC
MGP
CPu STh
SC
VA
VL
AGm
SNR
LGP
STh AGl
MGP
CPu
SNC
BG1 EC
VA
LGP
VL
MGP
CPu
AGm
AGl
SNR
SNC STh
SC
VA
LGP
VL
STh
AGl
CPu AGm
MGP
SNR
SNC
Direct EC
VA
LGP
VL
AGl
MGP
AGm
STh
SNR
SNC CPu
SC
VA
LGP
VL
STh
AGl
CPu AGm
MGP
SNR
SNC
Indirect EC
Brain Struct Funct
123
123
0.62000
2.66667
2.22222
2.33333
2.55556
2.77778
2.88889
2.77778
2.66667
2.55550
SNR
VA
VL
LGP
MGP
CPu
STh
AGl
AGm
3.04255
3.00000
3.08511
3.09574
3.04255
2.79787
2.85106
2.61702
3.24468
3.24468
RADout
0.45600
0.41500
0.38500
0.50200
0.50800
0.17000
0.23100
0.06700
0.67600
0.74700
Hub
3.07580
3.01080
3.00720
3.07580
3.09390
2.72560
2.86280
2.54870
3.14800
3.15520
RADout
0.37400
0.33500
0.31000
0.40600
0.40900
0.13300
0.18700
0.05300
0.54300
0.60400
Hub
VA
VL
AGm
LGP
AGl
SNR
STh
MGP
SNC
CPu
RADout
VA
VL
AGm
LGP
SNR
AGl
STh
MGP
SNC
CPu
Hub
VA
LGP
VL
AGl
AGm
MGP
STh
CPu
SNR
SNC
RADout
VA
LGP
VL
STh
AGl
AGm
CPu
MGP
SNR
SNC
Hub
VA
LGP
VL
STh
AGl
AGm
CPu
MGP
SNR
SNC
RADout
VA
LGP
VL
STh
AGl
AGm
CPu
MGP
SNR
SNC
Hub
2.55600
2.55600
2.77800
3.00000
2.66700
2.44400
2.66700
2.33300
2.66700
2.55600
RADin
0.63900
0.48800
0.92200
1.00000
0.82100
0.58600
0.80600
0.36200
0.83700
0.67300
Aut
3.01100
3.20200
3.19100
2.97900
3.01100
2.89400
3.04300
2.74500
3.01100
3.11700
RADin
0.49300
0.42900
0.66400
0.68800
0.42900
0.33300
0.58000
0.18600
0.48800
0.63700
Aut
3.03250
3.02890
3.09390
3.11550
3.00000
2.95670
3.05050
2.80140
3.01440
3.07580
RADin
0.43900
0.38300
0.57600
0.59700
0.37600
0.29600
0.50900
0.16200
0.42600
0.55600
Aut
VA
LGP
AGm
AGl
SNC
MGP
VL
SNR
STh
CPu
RADin
VA
AGl
LGP
AGm
SNC
VL
MGP
SNR
STh
CPu
Aut
VA
LGP
CPu
AGm
MGP
SNR
VL
SNC
STh
AGl
RADin
VA
LGP
AGl
MGP
SNR
AGm
VL
SNC
STh
CPu
Aut
VA
LGP
MGP
SNR
AGl
AGm
VL
SNC
STh
CPu
RADin
VA
LGP
MGP
AGl
SNR
AGm
VL
SNC
STh
CPu
Aut
EC eigenvector centrality, PRC page rank centrality, Hub hubness, Aut authoritativeness, RADout radiality of the output, RADin radiality of the input
The BG1-network was embedded in a network that contains all extrinsic regions that have direct input and output (direct) or indirect input and output (indirect) connections to the regions of the BG1-network. The italic short names are presenting the rank of a particular region with regard to a descending sorted specific local parameter (maximum rank 1 for a region; smallest value: rank 10). For example the CPu has the rank 1 concerning the betweenness centrality (BC) in the BG1-network. In the extrinsic direct neighbor network it has the rank 4 and in the indirect neighbor network it has rank 1. The subgraph centralities (SC) were normed
2.77778
SNC
RADout
0.87100
AGm
0.89800
MGP
AGl
0.73700
LGP
1.00000
0.36900
VL
0.87400
0.25800
VA
STh
0.79400
SNR
CPu
0.97200
SNC
Hub
Table 14 continued
Brain Struct Funct
11
11
11
10
12
11
12
11
12
12
12
VM
Ent
CM
HIPP
MDL
CL
LGF
PC
MDM
VL
VA
MGP
10
12
AGI
10
14
SNR
Pir
12
SNC
LHb
12
AGm
9
14
MRF
11 11
12
VTA
PF Ac
11
CPu
STh
8
11
AC
HD
Region
0
3
4
6
4
4
4
6
5
8
7
5
8
11 10
12
11
14
11
11
15
15
17
18
20
REC
5
14
14
14
16
16
17
20
20
21
21
22
22
28 28
30
30
31
32
32
34
38
40
41
42
DGA
Table 15 Local parameters of the ipsilateral BG2 network
2
3
4
6
5
8
4
11
10
9
9
8
11
17 13
17
17
15
19
21
17
23
23
20
22
DGC
3
11
10
8
11
8
13
9
10
12
12
14
11
11 15
13
13
16
13
11
17
15
17
21
20
DGI
0.6
0.79
0.71
0.57
0.69
0.5
0.76
0.45
0.5
0.57
0.57
0.64
0.5
0.39 0.54
0.43
0.43
0.52
0.41
0.34
0.5
3.39
0.43
0.51
0.48
CDC
3.04
11.47
10.08
8.94
10.97
8.45
13.37
9.2
10.73
11.85
12.37
13.67
11.51
11.69 14.02
13.19
12.92
14.76
13.17
11.33
15.75
14.86
16.05
18.79
18.82
Katz
3
1
2
2
2
2
2
1
2
1
2
1
1
2 1
1
1
1
2
1
1
1
1
1
1
SPC
15
142
135
148
202
188
205
264
275
284
325
299
367
551 492
600
562
499
603
624
558
840
864
886
950
Triag
3.17
0.12
0.11
0.14
0.15
0.13
0.13
0.1
0.11
0.11
0.12
0.09
0.13
0.12 0.11
0.11
0.11
0.09
0.1
0.1
0.08
0.1
0.09
0.09
0.09
CyclC
3
3
2
2
2
2
2
2
2
2
2
3
2
2 2
2
2
2
2
2
2
2
2
2
2
EccO
Brain Struct Funct
123
EccI
2
2
2
2
2
2
2
2
2
2
2 2
2
2
2
2
2
3
2
2
2
2
3
2
3
Region
AC
CPu
VTA
MRF
AGm
SNC
SNR
AGI
MGP
STh
PF Ac
LHb
Pir
VM
Ent
CM
HIPP
MDL
CL
LGF
PC
MDM
VL
VA
Table 15 continued
123
1
1
0.92
0.83
1
0.75
0.92
0.68
0.68
0.71
0.78
0.61
0.75
0.68 0.71
0.67
0.58
0.57
0.5
0.54
0.48
0.54
0.53
0.58
0.54
CCO
0.33
0.81
0.78
0.88
0.82
0.88
0.76
0.81
0.84
0.69
0.84
0.71
0.9
0.86 0.64
0.78
0.79
0.54
0.81
0.89
0.54
0.73
0.61
0.53
0.59
CCI
0.75
0.81
0.78
0.85
0.83
0.82
0.76
0.67
0.69
0.72
0.74
0.67
0.75
0.68 0.59
0.65
0.56
0.56
0.5
0.54
0.49
0.54
0.53
0.51
0.54
CCA
0.75
0.81
0.78
0.86
0.87
0.81
0.78
0.72
0.74
0.7
0.8
0.66
0.82
0.75 0.67
0.71
0.66
0.55
0.62
0.64
0.51
0.61
0.57
0.55
0.56
CCT
0.47
0.27
0.41
0.27
0.36
0.32
0.17
0.37
0.13
0.35
0.26
0.36
0.33
0.07 0.57
0.07
0.65
0.52
0.83
0.5
0.75
0
0
0
0
CC2
33.6
32.82
30.9
35.33
31.42
31.67
31.85
28.79
30
29.85
30.43
28.24
30.07
29.18 26.56
28.33
25.95
26.06
24.81
25.67
24.53
25.43
25.35
24.91
25.68
ADnb
0.12
0.22
0.24
0.13
0.25
0.2
0.22
0.33
0.26
0.29
0.27
0.3
0.26
0.28 0.35
0.3
0.39
0.41
0.41
0.36
0.42
0.34
0.34
0.37
0.32
VCDG
-0.74
-0.39
-0.36
-0.43
-0.3
-0.31
-0.29
-0.15
-0.18
-0.15
-0.16
-0.1
-0.14
0 0.06
0.05
0.11
0.13
0.16
0.14
0.2
0.22
0.24
0.27
0.26
Lev
0.13
0.39
0.35
0.29
0.47
0.46
0.46
0.52
0.54
0.5
0.54
0.69
0.58
0.68 0.72
0.74
0.75
0.64
0.79
0.85
0.69
0.99
0.99
0.96
0.94
Loc
0.5
0.52
0.55
0.57
0.56
0.6
0.55
0.65
0.63
0.62
0.62
0.59
0.65
0.77 0.69
0.77
0.77
0.73
0.83
0.89
0.77
0.96
0.96
0.86
0.92
CEO
0.52
0.65
0.62
0.6
0.65
0.6
0.69
0.59
0.63
0.67
0.67
0.71
0.65
0.65 0.73
0.69
0.69
0.75
0.69
0.65
0.77
0.73
0.77
0.89
0.86
CEI
Brain Struct Funct
1.35
0.95
SNC
AGI
1.35
AGm
0.32
0.32
MRF
SNR
0.55
VTA
ZA
Region
1.35
0.0001
VA
1.43
0.0008
VL
CPu
0.0013
MDM
AC
0.0017
0.0009
0.0017
CL
LGP
0.0025
MDL
PC
0.0042
0.0048
CM
HIPP
0.63
0.66
0.67
0.62
0.67
0.67
0.6
0.67
PCO
0.11
0.17
0.17
0.34
0.29
0.41
0.23
0.49
0.33
0.6
0.68
0.45
0.0214
Ac
0.82
0.86
0.45 0.44
0.0073
PF
0.0034 0.0086
3.0175
STh
0.8
0.69
VM Ent
0.0239
MGP
0.0057
0.0487
AGI
0.77
0.0105
0.0279
SNR
0.92
0.69
LHb
0.0177
SNC
1
Pir
0.0694
AGm
1
0.0510
0.0287
VTA
MRF
0.92
0.0923
CPu
0.97
0.0739
AC
EC
BC
Region
Table 15 continued
0.63
0.6
0.51
0.65
0.66
0.66
0.66
0.67
PCI
1,668.26
8,992.33
8,278.07
15,064.76
16,085.3
17,062.22
15,199.19
20,626.62
23,714.09
27,703.34 23,826.42
20,854.48
32,238.35
44,068.22
47,202.32
52,625.08
47,865.52
47,501.37
50,872.73
48,687.55
50,691.03
69,052.81
74,470.26
80,557.24
84,751.47
SC
0.63
0.66
0.65
0.64
0.67
0.66
0.66
0.67
PCA
0.18
0.45
0.44
0.35
0.41
0.33
0.53
0.4
0.41
0.46 0.57
0.6
0.43
0.6
0.43
0.54
0.53
0.86
0.54
0.45
0.9
0.58
0.74
1
0.94
PRC
2.63
2.79
2.88
2.71
2.96
2.96
2.83
2.92
RadO
0.3
0.27
0.4
0.6
0.39
0.45
0.31
0.51
0.45
0.55 0.64
0.39
0.62
0.6
0.65
0.68
0.61
0.83
0.56
0.52
0.8
0.65
0.74
0.79
0.91
FC
2.67
2.54
2.46
2.71
2.63
2.71
2.88
2.83
RadI
1
5
7
6
6
10
12
20
22
17 30
35
19
67
35
61
71
118
82
62
141
113
154
208
187
Stress
-8
-5
-2
-6
0
0
-3
-1
CenO
0.42
0.25
0.2
0.15
0.16
0.14
0.19
0.11
0.12
0.11 0.05
0.07
0.07
0.01
0.06
0.02
-0.02
-0.28
-0.04
0.03
-0.24
-0.01
-0.1
-0.31
-0.16
Shapley
-5
-9
-10
-4
-6
-4
1
-1
CenI
-1.24
-0.75
-1.32
-0.75
-2.06
-0.75
-1.88
0.9
-0.27
-0.27 0.35
0.35
-1.15
-0.21
0.69
0.69
0.69
1.19
-0.23
0.89
1.67
0.69
0.9
1.19
0.9
ZO
0.67
0.81
0.93
0.72
1
1
0.89
0.96
Hub
-1.51
-0.24
0.73
-0.66
-1.89
-0.66
-1.6
-0.44
-0.66
-0.24 0.73
-0.44
-0.38
-0.44
-0.38
-0.38
1.13
1.46
1.13
1.13
1.03
-0.38
-0.44
1.46
1.9
ZI
0.76
0.77
0.67
0.82
0.87
0.89
1
1
Aut
Brain Struct Funct
123
123
0.55
-2.09
-0.71
-2.21
-0.71
MDL
CL
LGP
PC
0
0
0.38
0.5
0.32
0.66
0.5
0.51
0.64
0.41
0.63
0.66 0.56
0.44
0.66
0.54
0.66
0.64
0.59
0.64
0.49
0.62
0.63
0.61
0.66
0.56
0.6 0.66
0.64
0.6
PCI
0.48
0.62
0.52
0.64
0.59
0.63
0.66
0.56
0.67
0.57
0.65
0.62
0.6
0.64 0.63
0.65
0.64
PCA
2
2.08
2.17
2.25
2.21
2.33
2.17
2.46
2.42
2.38
2.38
2.29
2.46
2.71 2.54
2.71
2.71
RadO
2.08
2.46
2.38
2.33
2.46
2.33
2.54
2.29
2.42
2.5
2.5
2.58
2.46
2.46 2.63
2.54
2.54
RadI
-23
-20
-20
-18
-19
-16
-20
-12
-14
-14
-15
-16
-12
-7 -10
-6
-6
CenO
-19
-10
-12
-14
-11
-14
-9
-13
-12
-9
-10
-7
-10
-11 -6
-8
-8
CenI
A all, AD average degree, C circle, CC cluster-coefficient, CE closeness centrality, DG degree, I in, O out, SP length of shortest path
The regions are sorted using the degree all. Abbreviations and parameter definitions are described in the Supporting Information
-1.4
-0.48
CM
HIPP
VA
0.55
-0.77
-0.25
VM
Ent
-0.48
0.11
MDM
0.64
-0.95
LHb
Pir
VL
0.44
0.32 -0.33
PF Ac
0.65
0.32
0.63
0.95
PCO
STh
ZA
MGP
Region
Table 15 continued
0.1
0.15
0.19
0.31
0.23
0.38
0.2
0.52
0.48
0.32
0.45
0.36
0.51
0.78 0.58
0.78
0.74
Hub
0.18
0.7
0.59
0.56
0.69
0.55
0.8
0.55
0.66
0.65
0.77
0.8
0.71
0.73 0.82
0.78
0.77
Aut
Brain Struct Funct
Brain Struct Funct
where qj;k is the number of shortest paths from j to k and qj;k ðiÞ is the number of shortest paths from j to k that pass through i. The directed and weighted definitions are the same. ^ [ 1 be a subset of Knotty centrality KC: Let N^ N, jNj N. Then P ^ ¼ LdðNÞ ^ Pi2N^ BCðiÞ KCðNÞ ð76Þ i2N BCðiÞ ^ of the subgraph N. ^ The knotty with the line density LdðNÞ center of a graph G is a subset NKC of nodes with ^ KCðNKC Þ ¼ min fKCðNÞg ^ NN ^ [1 jNj
ð77Þ
KCðNKC Þ is called the knotty-centerdness of the graph G. The knotty-centrality of a node i is defined as KCðiÞ ¼
1;
i 2 NKC
0;
else
ð78Þ
Stress S: SðiÞ ¼
X
qj;k ðiÞ
j;k2Nnfig
ð79Þ
The directed and weighted definitions are the same. Central point distance CPD: n 1 X BCmax BCðiÞ CPD ¼ ð80Þ n 1 i¼1 BCmax where BCmax ¼ maxi2N fBCðiÞg is the maximum Betweenness centrality. The directed and weighted versions use the directed and weighted Betweenness centralities. Participation coefficient: The partition M¼ fM1 ; . . .Mm g is generated as described in the definition of modularity. X deg ði; Mj Þ2 x ! PCx ðiÞ ¼ 1 ð81Þ degx ðiÞ Mj 2M with x 2 fin; out; allg and X degin ði; Mj Þ ¼ aki k2Mj nfig
(Number of edges from vertices of Mj to i). X aik degout ði; Mj Þ ¼ k2Mj nfig
(Number of edges from i to vertices of Mj ). X ðaik þ aki Þ degall ði; Mj Þ ¼ k2Mj nfig
ð82Þ
(Number of edges between i and vertices of Mj ). X degw ði; Mj Þ2 ! w x PCx ðiÞ ¼ 1 degwx ðiÞ Mj 2M
ð85Þ
with the same x and weighted definitions of degrees. One has 0 PCðiÞ 1. If PCðiÞ ¼ 1, the node i has no edges (in, out, all). If PCðiÞ ¼ 0 all edges (in, out all) come from, go to or stay in the same cluster. The larger PCðiÞ the more clusters are involved in the edges of node i. Z score/within module degree: Let Mi be the module containing node i. degx ði; Mi Þ x 2 fin; out; allg is defined in the participation coefficient. 1 X degx ðMi Þ ¼ degx ðj; Mi Þ ð86Þ jM j i j2Mi
is the mean and vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi !2ffi u u 1 X rdegx ðMi Þ ¼ t degx ðj; Mi Þ degx ðMi Þ jMi j j2Mi
ð87Þ
the standard deviation of the within module Mi degree distribution. Then the Z score is defined as Zx! ðiÞ ¼
degx ði; Mi Þ degx ðMi Þ rdegx ðMi Þ
ð88Þ
and analogous ! degwx ði; Mi Þ degw x ðMi Þ Zxw ðiÞ ¼ rdegwx ðMi Þ
ð89Þ
with the weighted versions of the mean and standard deviation. A value above one or below minus one implies that a node has significantly more or less edges from, to or from and to nodes in its cluster than the average node in its cluster has. Eigenvector centrality: The eigenvector centrality ECðiÞ is the i-th component of the eigenvector with the largest corresponding eigenvalue of the adjacency matrix resp. weight matrix. Shapley rating /: The Shapley rating is a measure that provides information about the loss of connectivity following the removal of a node. X ^ SRðiÞ ¼ ðjSCCðN^ [ figÞj jSCCðNÞjÞ ^ NNnfig
ð83Þ
ð84Þ
ð90Þ ^ 1Þ! jNj! ^ ðn jNj n! ^ where SCCðNÞ is the set of strongly connected components ^ The smaller the value is, the more important is the of N. node in the sense of connectivity of the graph. Because of the exponential number of subsets, this parameter can be approximated for large networks, only.
123
Brain Struct Funct
Radiality: The radiality of a node Rad is a measure of the distance of a node to all other nodes. Nodes that have a small radiality have larger distances to other nodes than those with a greater radiality. Input radiality Radin : The input radiality of a node Radin is X 1 Diam þ 1 dðj; iÞ Radin ðiÞ ¼ ð91Þ n 1 j2N dðj;iÞ\1
In the weighted case the weighted distances are used. Output radiality Radout : The output radiality of a node Radout is X 1 Radout ðiÞ ¼ Diam þ 1 dði; jÞ ð92Þ n 1 j2N dði;jÞ\1
In the weighted case the weighted distances are used. Centroid value Cen: With gout ði; jÞ ¼ jfk 2 Njdði; kÞ\dðj; kÞ\1gj and gin ði; jÞ ¼ jfk 2 Njdðk; iÞ\dðk; jÞ\1gj which are the number of nodes closer to node i than to node j with regard to In- and Outdistance, the centroid value is defined in the following. Output centroid value Cenout : Cenout ðiÞ ¼ minfgout ði; jÞ gout ðj; iÞjj 2 Nnf igg
ð93Þ
FCðiÞ ¼
X 1 aji aik jNi j ðjNi j 1Þ j;k2N
In the weighted case we define the flow coefficient as the sum of weights of paths of length 2 between neighbors of a node i that pass node i divided by the maximum possible sum. X X ! 1 FC w ðiÞ ¼ ðwji þ wik Þ 2 jNi j ðjNi j 1Þ j2Ni k2N nfjg
Cenin ðiÞ ¼ minfgin ði; jÞ gin ðj; iÞjj 2 Nnf igg
ð94Þ
In the weighted case the weighted distances are used. A value \0 implies, that there exists a node that is closer to most other nodes. A value 0 implies, that this node is most central in the network. A value ¼0 implies, that there are more than one most central nodes. Page rank centrality PRC: PRCðiÞ ¼ ri where r is the solution of the linear system 0 1 1 1 B .. C T ð95Þ ðI a A BÞ r ¼ ð1 aÞ @ . A n 1 with the damping factor a ¼ 0:85, the identity matrix I and the diagonal matrix B, whereby 8 1 < ; degout ðiÞ [ 0 bii ¼ degout ðiÞ ð96Þ : 0; otherwise In the weighted case the weight matrix W is used instead of ! A and the weighted version deg w ðiÞ of the outdegree. out
Flow coefficient FC: Number of paths of length 2 between neighbors of a node i that pass node i divided by the maximum possible numbers of sub paths.
123
i
wji [ 0
wik [ 0
ð98Þ Average flow coefficient FC: Number of paths of length 2 between neighbors of a node i that pass node i divided by the maximum possible numbers of sub paths. FC ¼
n 1X FCðiÞ n i¼1
ð99Þ
n ! ! 1X FC w ðiÞ FC w ¼ n i¼1
ð100Þ
Subgraph centrality SC: SCðiÞ ¼
1 X ðAk Þ k¼0
Input centroid value Cenin :
ð97Þ
i
j6¼k
ii
ð101Þ
k!
1 X ! ðW k Þii SC w ðiÞ ¼ k! k¼0
ð102Þ
The subgraph centrality of the network is the average subgraph centrality of its nodes. SC ¼
n 1X SCðiÞ n i¼1
ð103Þ
n ! 1X ! SC w ¼ SC w ðiÞ n i¼1
ð104Þ
Undirected cyclic coefficient CyclC: The undirected cyclic coefficient as published by (Kim et al. 2005) Cyclic topology in complex networks. X 2 1 CyclCðiÞ ¼ jNi j ðjNi j 1Þ ðj;kÞ2N N 2 þ disti ðj; kÞ i
i
j6¼k
ð105Þ With 8 length of the shortest path from j to k > > > < that does not contains i; disti ðj; kÞ ¼ > if such a path exists > > : 1; otherwise ð106Þ
Brain Struct Funct
Directed cyclic coefficient CyclC! : A publication about the directed cyclic coefficient is unknown. The directed cyclic coefficient is implemented here as follows: CyclC! ðiÞ ¼
jNiout j
1 jNiout \ Niin j 1 2 þ disti ðj; kÞ in
jNiin j
X
ð107Þ
ðj;kÞ2Niout Ni j6¼k
! Directed weighted cyclic coefficient CyclC w : A publication about the directed weighted cyclic coefficient is unknown. The directed weighted cyclic coefficient is implemented here as follows: ! CyclC w ðiÞ ¼
1 jNiout j jNiin j jNiout \ Niin j X 1 wij þ wki þ distwi ðj; kÞ out in
ð108Þ
ðj;kÞ2Ni Ni j6¼k
with distwi ðj; kÞ is the weighted version of disti ðj; kÞ with the weighted path length. Cyclic network coefficient CyclC! : The cyclic coefficient of the network is the average cyclic coefficient of its nodes: n 1X CyclC! ¼ CyclC! ðiÞ ð109Þ n i¼1 Hubness and authoritativeness: A hub is a node that points to many authorities and an authority is a node that has numerous input connections from many hubs (Kleinberg 1999; Sporns et al. 2007). The hubness HubðiÞ of a node i is: X HubðiÞ ¼ AuthðjÞ ð110Þ out j2Ni
with the authoritativeness AuthðiÞ X AuthðiÞ ¼ HubðjÞ j2Niin
ð111Þ
An iterative algorithm is used to calculate a fixed point of these equations. Vulnerability V: The vulnerability V is the maximum relative decrease of the global efficiency removing a single node. GE GEðiÞ V ¼ max ð112Þ i2N GE where GEðiÞ is the global efficiency of the graph (Nnfig; fðj; kÞ 2 Ejj 6¼ i 6¼ kg) that originates by removal of node i and all edges adjacent to i. The weighted version is analog using the weighted global efficiencies.
Random models The following random graph models are compared to the real network of the intrinsic amygdala connectivity. By comparing the average path length and the cluster coefficient of the models with the real network it is feasible to determine a model that is most similar to the real network. Erdo¨s Re´nyi graph: ð113Þ
Gðn; pÞ
where n is the number of vertices and p is the probability that an edge ði; jÞ exists, for all i; j. The degree distribution of the Erdo¨s Re´nyi random graph is binomial in terms of n1 k PðdegðvÞ ¼ kÞ ¼ p ð1 pÞn1k ð114Þ k Watts–Strogatz graph: The small-world model of Watts– Strogatz is a random graph generation model that provides graphs with small-world properties. The network (initially it has a non-random lattice structure) is build by linking each node to its hki closest neighbors using a rewiring probability p. Hence, an edge has the probability p that it will be rewired as a random edge. The number of rewired links can be estimated by: pE ¼ pNhki=2
ð115Þ
Barabasi–Albert graph: The Barabasi–Albert graph is used to generate preferential attachments between nodes. The probability pi that the new node is connected to node i is ki pi ¼ P
j kj
ð116Þ
The degree distribution of a Barabasi–Albert network is scale free following the power law distribution of the form: PðkÞ k3
ð117Þ
Eipert graph: The modified Eipert model (EN: Eipert network) is based on the Barabasi–Albert graph. However, the algorithm starts at a fixed number of nodes and edges are added iteratively. Ozik–Hunt–Ott graph: The Ozik–Hunt–Ott model (OHO) (Ozik et al. 2004) is a small-world randomization approach that was modified for directed networks and a fixed number of edges. The OHO-model uses a growing mechanism in which all connections are made locally to topographical nearby regions. Rewiring graph: The rewiring-models connects each target of an edge of a network to another target node. Power law: PðkÞ ¼ a kc
ð118Þ
D is the deviation (error) of an empirical distribution of degrees from the power law function. A small D value
123
Brain Struct Funct
means that the empirical distribution is similar with the power law function.
References Albert R, Barabasi AL (2002) Statistical mechanics of complex networks. Rev Mod Phys 74:47–97 Albin RL, Young AB, Penney JB (1989) The functional anatomy of basal ganglia disorders. Trends Neurosci 12:366–375 Alexander GE, Crutcher MD (1990) Functional architecture of basal ganglia circuits: neural substrates of parallel processing. Trends Neurosci 13:266–271 Archambault D, Munzner T, Auber D (2008) Grouseflocks: steerable exploration of graph hierarchy space. IEEE Trans Vis Comput Graph 14:900–913 Bai X, Yu L, Liu Q, Zhang J, Li A, Han D, Luo Q, Gong H (2006) A high-resolution anatomical rat atlas. J Anat 209:707–708 Bakker R, Wachtler T, Diesmann M (2012) CoCoMac 2.0 and the future of tract-tracing databases. Front Neuroinform 6:1–6 Barbas H, Hilgetag CC, Saha S, Dermon CR, Suski JL (2005) Parallel organization of contralateral and ipsilateral prefrontal cortical projections in the rhesus monkey. BMC Neurosci 6:32 Bassett DS, Bullmore E (2006) Small-world brain networks. Neuroscientist 12:512–523 Bennett BD, Bolam JP (1994) Synaptic input and output of parvalbumin-immunoreactive neurons in the neostriatum of the rat. Neuroscience 62:707–719 Bevan MD, Booth PA, Eaton SA, Bolam JP (1998) Selective innervation of neostriatal interneurons by a subclass of neuron in the globus pallidus of the rat. J Neurosci 18:9438–9452 Bjaalie JG (2002) Localization in the brain: new solutions emerging. Nat Rev Neurosci 3:322–325 Bohland J, Bokil H, Allen C, Mitra P (2009a) The brain atlas concordance problem: quantitative comparison of anatomical parcellations. PloS One 4:e7200 Bohland JW, Wu C, Barbas H, Bokil H, Bota M, Breiter HC, Cline HT, Doyle JC, Freed PJ, Greenspan RJ, Haber SN, Hawrylycz M, Herrera DG, Hilgetag CC, Huang JZ, Jones A, Jones EG, Karten HJ, Kleinfeld D, Ko¨tter R, Lester HA, Lin JM, Mensh BD, Mikula S, Panksepp J, Price JL, Safdieh J, Saper CB, Schiff ND, Schmahmann JD, Stillman BW, Svoboda K, Swanson LW, Toga AW, Van Essen DC, Watson JD, Mitra PP (2009b) A proposal for a coordinated effort for the determination of brainwide neuroanatomical connectivity in model organisms at a mesoscopic scale. PLoS Comput Biol 5:1–9 Bosch C, Mailly P, Degos B, Deniau JM, Venance L (2012) Preservation of the hyperdirect pathway of basal ganglia in a rodent brain slice. Neuroscience 215:31–41 Bota M, Swanson L (2007a) The neuron classification problem. Brain Res Rev 56:79–88 Bota M, Swanson L (2007b) Online workbenches for neural network connections. J Comp Neurol 500:807–814 Bota M, Swanson L (2010) Collating and curating neuroanatomical nomenclatures: principles and use of the brain architecture knowledge management system (BAMS). Front Neuroinf 4:1–16 Bounova G, de Weck O (2012) Overview of metrics and their correlation patterns for multiple-metric topology analysis on heterogeneous graph ensembles. Phys Rev E 85:1–11 Brandes U, Erlebach T (2005) Network analysis. Springer, Berlin Bullmore E, Sporns O (2009) Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci 10:186–198
123
Burns G (1997) Neural connectivity of the rat: theory, methods and applications. PhD thesis, University of Oxford Burns G, Cheng WC (2006) Tools for knowledge acquisition within the neuroscholar system and their application to anatomical tract-tracing data. J Biomed Discov Collab 1:10–16 Burns G, Young M (2000) Analysis of the connectional organization of neural systems associated with the hippocampus in rats. Phil Trans R Soc Lond Ser B 355:55–70 Burns G, Cheng WC, Thompson R, Swanson L (2006) The NeuARt II system: a viewing tool for neuroanatomical data based on published neuroanatomical atlases. BMC Bioinf 7:531–550 Burns G, Cheng WC, Thompson R, Swanson L (2008a) The NeuARt II system: a viewing tool for neuroanatomical data based on published neuroanatomical atlases. LNCS 5151:9–18 Burns G, Feng D, Hovy E (2008b) Intelligent approaches to mining the primary research literature: techniques, systems, and examples. In: Computational intelligence in medical informatics. Studies in computational intelligence, vol 85. Springer, Berlin, pp 17–50 Chida Y, Toyosawa K (1994) Laterality of theta waves recorded from the bilateral hippocampi in freely moving male rats. J Vet Med Sci 56:407–410 Colizza V, Flammini A, Serrano MA, Vespignani A (2006) Detecting rich-club ordering in complex networks. Nat Phys 2:110–115 Costa LdF, Sporns O (2006) Correlating thalamocortical connectivity and activity. Appl Phy Lett 89:1–3 DeLong MR (1990) Primate models of movement disorders of basal ganglia origin. Trends Neurosci 13:281–285 Diesmann M, Gewaltig MO, Aertsen A (1999) Stable propagation of synchronous spiking in cortical neural networks. Nature 402:529 Echtermeyer C, Costa L, Rodrigues F, Kaiser M (2011) Automatic network fingerprinting through single-node motifs. Plos One 6:1–9 Ercsey-Ravasz M, Markov NT, Lamy C, Van Essen DC, Knoblauch K, Toroczkai Z, Kennedy H (2013) A predictive network model of cerebral cortical connectivity based on a distance rule. Neuron 80:184–197 Estrada E (2010) Quantifying network heterogeneity. Phys Rev E 82(066):102 Estrada E, Hatano N (2008) Communicability in complex networks. Phys Rev E 77(036):111 Fagiolo G (2007) Clustering in complex directed networks. Phys Rev E 76(026):107 Felleman D, Essen DV (1991) Distributed hierarchical processing in the primate cerebral cortex. Cereb Cortex 1:1–47 Feng G, Burton N, Hill B, Davidson D, Kerwin J, Scott M, Lindsay S, Baldock R (2005) JAtlasView: a Java atlas-viewer for browsing biomedical 3D images and atlases. BMC Bioinform 6:1–7 Foster KC, Muth SQ, Potterat JJ, Rothenberg RB (2001) A faster Katz status score algorithm. Comput Math Organ Theory 7:275–285 Freeman LC (1977) A set of measures of centrality based on betweenness. Sociometry 40:35–41 Gerfen C (2004) Basal ganglia. In: The rat nervous system. Elsevier Academic Press, San Diego, pp 455–508 Gerfen C, Bolam J (2010) The neuroanatomical organization of the basal ganglia. In: Handbook of basal ganglia structure and function: a decade of progress. Elsevier Academic Press, San Diego, pp 3–28 Gewaltig MO, Diesmann M (2007) Nest (neural simulation tool). Scholarpedia 2:1430 Gouws A, Woods W, Millman R, Morland A, Green G (2009) DataViewer3D: an open-source, cross-platform multi-modal neuroimaging data visualization tool. Front Neuroinf 3:1–18 Gustafson C, Tretiak O, Bertrand L, Nissanov J (2004) Design and implementation of software for assembly and browsing of 3D brain atlases. Comput Methods Programs Biomed 74:53–61
Brain Struct Funct Gustafson C, Bug W, Nissanov J (2007) Neuroterrain—a client– server system for browsing 3d biomedical image data sets. BMC Bioinform 8 Hagmann P, Cammoun L, Gigandet X, Meuli R, Honey CJ, Wedeen VJ, Sporns O (2008) Mapping the structural core of human cerebral cortex. PLoS Biol 6:1–15 Harriger L, van den Heuvel MP, Sporns O (2012) Rich club organization of macaque cerebral cortex and its role in network communication. PLoS One 7(e46):497 Hilgetag CC, Grant S (2000) Uniformity, specificity and variability of corticocortical connectivity. Philos Trans R Soc Lond B Biol Sci 355:7–20 Hjornevik T, Leergaard T, Darine D, Moldestad O, Dale A, Willoch F, Bjaalie J (2007) Three-dimensional atlas system for mouse and rat brain imaging data. Front Neuroinf 1:1–11 Honey C, Ko¨tter R, Breakspear M, Sporns O (2007) Network structure of cerebral cortex shapes functional connectivity on multiple time scales. Proc Natl Acad Sci USA 104:10240–10245 Humphries M, Gurney K, Prescott T (2006) The brainstem reticular formation is a small-world, not scale-free, network. Proc R Soc B 273:503–511 Izhikevich EM, Edelman GM (2008) Large-scale model of mammalian thalamocortical systems. Proc Natl Acad Sci USA 105:3593–3598 Ju T, Warren J, Carson J, Bello M, Kakadiaris I, Chiu W, Thaller C, Eichele G (2006) 3d volume reconstruction of a mouse brain from histological sections using warp filtering. J Neurosci Methods 156:84–100 Keinan A, Hilgetag CC, Meilijson I, Ruppin E (2004) Causal localization of neural function: the shapley value method. Neurocomputing 58–60:215–222 Kita H (1993) GABAergic circuits of the striatum. Prog Brain Res 99:51–72 Kita H, Kitai ST (1994) The morphology of globus pallidus projection neurons in the rat: an intracellular staining study. Brain Res 636:308–319 Kleinberg JM (1999) Authoritative sources in a hyperlinked environment. J ACM 46:604–632 Koos T, Tepper JM (1999) Inhibitory control of neostriatal projection neurons by GABAergic interneurons. Nat Neurosci 2:467–472 Koslow S (1997) Neuroinformatic: an overview of the human brain project. Lawrence Erlbaum Associates, New Jersey Koslow S (2005) Databasing the brain: from data to knowledge. Wiley, New York Ko¨tter R (2003) Neuroscience databases: a practical guide. Kluwer Academic Publishers, Dordrecht Ko¨tter R (2004) Online retrieval, processing, and visualization of primate connectivity data from the CoCoMac database. Neuroinformatics 2:127–144 Ko¨tter R (2007) Anatomical concepts of brain connectivity. Springer, Berlin Ko¨tter R, Reid AT, Krumnack A, Wanke E, Sporns O (2007) Shapley ratings in brain networks. Front Neuroinf 1:1–9 Lanciego JL, Wouterlood FG (2011) A half century of experimental neuroanatomical tracing. J Chem Neuroanat 42:157–183 Lee J, Munch K, Carlis J, Pardo J (2008) Internet image viewer (iiv). BMC Med Imaging 29:1–20 Lima-Mendez G, van Helden J (2009) The powerful law of the power law and other myths in network biology. Mol BioSyst 5:1482–1493 Looi JC, Walterfang M, Nilsson C, Power BD, van Westen D, Velakoulis D, Wahlund LO, Thompson PM (2013) The subcortical connectome: hubs, spokes and the space between—a vision for further research in neurodegenerative disease. Aust N Z J Psychiatry EPub:1–4
Markov NT, Misery P, Falchier A, Lamy C, Vezoli J, Quilodran R, Gariel MA, Giroud P, Ercsey-Ravasz M, Pilaz LJ, Huissoud C, Barone P, Dehay C, Toroczkai Z, Van Essen DC, Kennedy H, Knoblauch K (2011) Weight consistency specifies regularities of macaque cortical networks. Cereb Cortex 21:1254–1272 Markov NT, Vezoli J, Chameau P, Falchier A, Quilodran R, Huissoud C, Lamy C, Misery P, Giroud P, Ullman S, Barone P, Dehay C, Knoblauch K, Kennedy H (2014) Anatomy of hierarchy: feedforward and feedback pathways in macaque visual cortex. J Comp Neurol 522(1):225–259 Martinez-Gonzalez C, van Andel J, Bolam JP, Mena-Segovia J (2013) Divergent motor projections from the pedunculopontine nucleus are differentially regulated in Parkinsonism. Brain Struct Funct 218 (epub ahead of print) Martone ME, Gupta A, Ellisman MH (2004) E-neuroscience: challenges and triumphs in integrating distributed data from molecules to brains. Nat Neurosci 7:467–472 Mascaro MB, Bittencourt JC, Casatti CA, Elias CF (2005) Alternative pathways for catecholamine action in oral motor control. Neurosci Lett 386:34–39 Modha D, Singh R (2010) Network architecture of the long-distance pathways in the macaque brain. PNAS 107:13485–13490 Moore E, Poliakov A, Lincoln P, Brinkley J (2007) MindSeer: a portable and extensible tool for visualization of structural and functional neuroimaging data. BMC Bioinform 8(389):1–12 Nambu A, Tokuno H, Takada M (2002) Functional significance of the cortico-subthalamo-pallidal ’hyperdirect’ pathway. Neurosci Res 43:111–117 Newman M (2010) Networks: an introduction. Oxford University Press Inc., New York Newman M, Girvan M (2004) Finding and evaluating community structure in networks. Phys Rev 69:1–15 Newman M, Barabasi AL, Watts DJ (eds) (2006) The structure and dynamics of networks. Princeton University Press, New Jersey Niggemann J, Gebert A, Schulz S (2008) Model formulation: modeling functional neuroanatomy for an anatomy information system. JAMIA 15:671–678 Oh SW, Harris JA, Ng L, Winslow B, Cain N, Mihalas S, Wang Q, Lau C, Kuan L, Henry AM, Mortrud MT, Ouellette B, Nguyen TN, Sorensen SA, Slaughterbeck CR, Wakeman W, Li Y, Feng D, Ho A, Nicholas E, Hirokawa KE, Bohn P, Joines KM, Peng H, Hawrylycz MJ, Phillips JW, Hohmann JG, Wohnoutka P, Gerfen CR, Koch C, Bernard A, Dang C, Jones AR, Zeng H (2014) A mesoscale connectome of the mouse brain. Nature 508:207–214 Ozik J, Hunt B, Ott E (2004) Growing networks with geographical attachment preference: emergence of small worlds. Phys Rev E 69(026):108 Palombi O, Shin JW, Watson C, Paxinos G (2006) Neuroanatomical affiliation visualization-interface system. Neuroinformatics 4:299–317 Parent A (1986) Comparative neurobiology of the basal ganglia. Wiley, New York Parthasarathy K, Bhalla US (2013) Laterality and symmetry in rat olfactory behavior and in physiology of olfactory input. J Neurosci 33:5750–5760 Passingham RE, Stephan KE, Kotter R (2002) The anatomical basis of functional localization in the cortex. Nat Rev Neurosci 3:606–616 Paxinos G, Watson C (2004) The rat brain in stereotaxic coordinates, 5th edn. Elsevier Academic Press, London Paxinos G, Watson C (2007) The rat brain in stereotaxic coordinates, 6th edn. Elsevier Academic Press, Amsterdam Paxinos G, Watson C (2009) BrainNavigator. Academic Press Inc., San Diego
123
Brain Struct Funct Press W, Olshausen B, Van Essen D (2001) A graphical anatomical database of neural connections. Phil Trans R Soc Lond B 356:1147–1157 Riddle DR, Purves D (1995) Individual variation and lateral asymmetry of the rat primary somatosensory cortex. J Neurosci 15:4184–4195 Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52:1059–1069 Sato F, Lavallee P, Levesque M, Parent A (2000) Single-axon tracing study of neurons of the external segment of the globus pallidus in primate. J Comp Neurol 417:17–31 Schmitt O, Eipert P (2012) Neuroviisas: approaching multiscale simulation of the rat connectome. Neuroinformatics 10:243–267 Schmitt O, Eipert P, Philipp K, Kettlitz R, Fuellen G, Wree A (2012) The intrinsic connectome of the rat amygdala. Front Neural Circuits 6:81 Schreiber F, Schwo¨bbermeyer H (2004) Towards motif detection in networks: frequency concepts and flexible search. In: Proceedings of the international workshop on network tools and applications in biology NETTAB04, pp 91–102 Schreiber F, Schwo¨bbermeyer H (2005) Mavisto: a tool for the exploration of network motifs. Bioinformatics 21:3572–3574 Shanahan M, Wildie M (2012) Knotty-centrality: finding the connective core of a complex network. PLoS One 7(e36):579 Shanahan M, Bingman VP, Shimizu T, Wild M, Gu¨ntu¨rku¨n O (2013) Large-scale network organization in the avian forebrain: a connectivity matrix and theoretical analysis. Front Comp Neurosci 7(89):1–17 Smith Y, Bevan MD, Shink E, Bolam JP (1998) Microcircuitry of the direct and indirect pathways of the basal ganglia. Neuroscience 86:353–387 Sporns O (2002) Graph theory methods for the analysis of neural connectivity patterns. In: Ko¨tter R (ed) Neuroscience databases. Kluwer Academic, Boston, pp 171–185 Sporns O (2006) Small-world connectivity, motif composition, and complexity of fractal neuronal connections. Biosystems 85:55–64 Sporns O (2011) Networks of the brain. The MIT Press, Cambridge Sporns O, Ko¨tter R (2004) Motifs in brain networks. PLoS Biol 2:1910–1918 Sporns O, Tononi G, Edelman GM (2000a) Connectivity and complexity: the relationship between neuroanatomy and brain dynamics. Neural Netw 13:909–922 Sporns O, Tononi G, Edelman GM (2000b) Theoretical neuroanatomy: relating anatomical and functional connectivity in graphs and cortical connection matrices. Cereb Cortex 10:127–141 Sporns O, Tononi G, Edelman GM (2002) Theoretical neuroanatomy and the connectivity of the cerebral cortex. Behav Brain Res 135:69–74 Sporns O, Chialvo DR, Kaiser M, Hilgetag CC (2004) Organization, development and function of complex brain networks. Trends Cogn Sci 8:418–425 Sporns O, Honey CJ, Ko¨tter R (2007) Identification and classification of hubs in brain networks. PLoS One 2:1–14 Stam CJ, Reijneveld JC (2007) Graph theoretical analysis of complex networks in the brain. Nonlinear Biomed Phys 1:1–19 Stephan K, Kamper L, Bozkurt A, Burns G, Young M, Ko¨tter R (2001a) Advanced database methodology for the collation of
123
connectivity data on the macaque brain (CoCoMac). Philos Trans R Soc Lond B Biol Sci 356:1159–1186 Stephan KE, Kamper L, Bozkurt A, Burns GAPC, Young MP, Ko¨tter R (2001b) Advanced database methodology for the Collation of Connectivity data on the Macaque brain (CoCoMac). Philos Trans R Soc Lond Ser B 356:1159–1186 Steriade M (2004) Neocortical cell classes are flexible entities. Nat Rev Neurosci 5:121–134 Sugar J, Witter M, van Strien N, Cappaert N (2011) The retrosplenial cortex: intrinsic connectivity and connections with the (para)hippocampal region in the rat. An interactive connectome. Front Neuroinf 5:1–13 Swanson L (2003) Brain maps: structure of the rat brain, vol 3. Elsevier, Amsterdam Tepper JM, Abercrombie ED, Bolam JP (2007) Basal ganglia macrocircuits. Prog Brain Res 160:3–7 Thibeault CM, Srinivasa N (2013) Using a hybrid neuron in physiologically inspired models of the basal ganglia. Front Comput Neurosci 7:88. (eCollection) Tiesinga P, Fellous JM, Sejnowski TJ (2008) Regulation of spike timing in visual cortical circuits. Nat Rev Neurosci 9:97–107 Towlson EK, Ve´rtes PE, Ahnert SE, Schafer WR, Bullmore ET (2013) The rich club of the C. elegans neuronal connectome. J Neurosci 33:6380–6387 van den Heuvel MP, Sporns O (2011) Rich-club organization of the human connectome. J Neurosci 31:15775–15786 Van Horn JD, Gazzaniga MS (2013) Why share data? Lessons learned from the fMRIDC. Neuroimage 82:677–682 Vlachos I, Aertsen A, Kumar A (2012) Beyond statistical significance: implications of network structure on neuronal activity. PLoS Comput Biol 8(e1002):311 Wang Q, Sporns O, Burkhalter A (2012) Network analysis of corticocortical connections reveals ventral and dorsal processing streams in mouse visual cortex. J Neurosci 32:4386–4399 Ware C (2013) Information visualization: perception for design. Morgan Kaufmann Publishers Inc., San Francisco Wong PC, Foote H, Chin G, Mackey P, Perrine K (2006a) Graph signatures for visual analytics. IEEE Trans Vis Comput Graph 12:1399–1413 Wong PC, Foote H, Mackey P, Perrine K, Chin G (2006b) Generating graphs for visual analytics through interactive sketching. IEEE Trans Vis Comput Graph 12:1386–1398 Yip AM, Horvath S (2007) Gene network interconnectedness and the generalized topological overlap measure. BMC Bioinform 8:22 Young M (1992a) Objective analysis of the topological organization of the primate cortical visual system. Nature 358:152–155 Young MP (1992b) Objective analysis of the topological organization of the primate cortical visual system. Nature 358:152–155 Young M (1993) The organization of neural systems in the primate cerebral cortex. Proc Biol Sci 252:13–18 Young M, Scannell J, Burns G, Blakemore C (1994) Analysis of connectivity: neural systems in the cerebral cortex. Rev Neurosci 5:227–250 Young M, Scannell J, Burns G (1996) The analysis of cortical connectivity. Springer, Berlin Zhou S, Mondragn RJ (2004) The rich-club phenomenon in the internet topology. IEEE Commun Lett 8:180–182