Impact Factor 3.005

The Frontiers in Neuroscience journal series is the 1st most cited in Neurosciences

Original Research ARTICLE

Front. Neural Circuits, 11 December 2012 |

The intrinsic connectome of the rat amygdala

  • 1Department of Anatomy, University of Rostock, Rostock, Germany
  • 2Institute for Biostatistics and Informatics in Medicine and Ageing Research, University of Rostock, Rostock, Germany

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

1. Introduction

A connectome is a level-dependent representation of connections between biological entities. Levels could be molecular, cellular (micro level), cohesive structural and/or functional ensembles of cells (meso level), or distinct groups of cohesive ensembles (macro level). With regard to the nervous system of the rat, we investigated data of a partial connectome of the central nervous system, the connectome of the amygdala (amygdaloid complex), and its extrinsic efferent targets as well as extrinsic afferent sources at the meso level of nuclei and their subdivisions. The intrinsic and extrinsic amygdala data are part of an evaluation of about 2100 publications conducted for this study and contain exclusively tract-tracing data of the rat nervous system (peripheral nervous system, spinal cord, brain). One rationale of this study is to uncover the connectivity of the intrinsic amygdala network of the rat based on high-resolution tract-tracing data. So far, the connectome of the intrinsic amygdala network has not been investigated in terms of network measures, graph theoretical quantities, and multivariate statistics. After making these connectome data accessible through a weighted connectivity matrix, specific questions concerning motifs, reciprocity, and most strongly connected regions upon others can be addressed.

Connectivity data sets of parts of nervous systems have been developed and analyzed by several groups (Felleman and Essen, 1991; Young, 1992, 1993; Scannell et al., 1999; Sporns et al., 2000, 2002; Sporns and Kötter, 2004; Sporns and Zwi, 2004; Honey et al., 2007; Modha and Singh, 2010; Sugar et al., 2011). Most of these connections are compiled into by meta-studies of tract-tracing publications. The sources of information in this study originate exclusively from peer-reviewed 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 adult rats. In our case, 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.

Partial connectomes of the rat nervous system have been elaborated for the retrosplenial cortex (Sugar et al., 2011), nucleus of the solitary tract (Palombi et al., 2006), reticular formation (Humphries et al., 2006), hippocampus (Burns and Young, 2000), and the visual system (Burns, 1997), for an overview of connectome studies we refer to Sporns (2011), Schmitt and Eipert (2012). A connectome of the amygdala of the rat is reported here for the first time, even though, reviews of the hodology of the amygdala are numerous (Swanson and Petrovich, 1998; Pitkänen, 2000; Sah et al., 2003; de Olmos et al., 2004). In contrast to low-resolution, non-directed, and non-weighted connectomes derived from tractographic analysis of in vivo DTI measurements (Essen et al., 2012), the type of connectome described in the following is a high-resolution, directed, and weighted network where nodes are organized in a neuroontology containing functional, neurophysiological, and molecular biological information.

Nodes of the network correspond to regions which are distinguishable in terms of cytoarchitecture. These regions can be related to each other with regard to hierarchical subdivisions. Since many possibilities exist to build such hierarchies of subdivisions of the amygdaloid complex that influence the analysis of its connectome, an overview will be provided in the following. The German translation Mandelkern of the term amygdala has been introduced by Burdach (1819–1826). A first topographic description of the amygdala was published by Meynert (1867). A nomenclature of subdivisions was introduced by Johnston (1923) and further developed by Pitkänen (2000) and de Olmos et al. (2004). Lists of terms of amygdala nuclei are contained in the stereotaxic atlases of Paxinos and Watson (1986) and Swanson (1992). A comprehensive comparison of nomenclatures was published by Price (1981) and Pitkänen (2000). The rat amygdaloid complex (amygdalar complex, amygdaloid body), amygdala in short, is a heterogeneous gray complex of 13 larger nuclear and cortical regions. This multinuclear complex is located in the depth of the anteromedial temporal lobe ventral to the lentiform nucleus (Figure 1). The topographical criteria proposed by Brockhaus (1938) to subdivide the amygdala into a superficial and a deep nuclear group has been adopted here as well as by other authors (Pitkänen, 2000; de Olmos et al., 2004) and will be described in detail in the next subsection.


Figure 1. Visualization of hierarchical regions of the rat central nervous system using neuroVIISAS. The subdivision of regions of the amygdala is based on the stereotactic atlas of Paxinos and Watson (2007). (A) Location of the amygdala (arrow) in the central nervous system of the rat. (B) The view from the midline into the right hemisphere shows the amygdala (magenta structures). Cortical regions are labeled with red tones [same orientation as shown in (A)]. This scheme of colors is consistent with the color mapping of the complete rat nervous system connectome where further colors are already assigned to other regions (see Materials and Methods). (C) The view from ventral (bottom) shows the rostrocaudal extension of the amygdala. The rough surfaces are caused by the intersection distance, serial sectioning, and small subregions of the amygdala.

In compliance with the structural diversity of the amygdala the various nuclear groups and subdivisions also differ functionally. Major outputs to functional systems generating experimental, motor, mnemonic, autonomic, or endocrine responses originate largely in different nuclei (Russchen, 1986). Four functional systems can be related to nuclear groups of the amygdala: the accessory olfactory, the main olfactory, the autonomic, and the frontotemporal system. The amygdala plays a crucial role in conditioned fear, anxiety, and attention. Especially the lateral, basolateral, and central nuclei of the amygdala build a functionally unified system necessary for the acquisition and expression of conditioned or instrumental fear learning and memory modulation (Davis, 1992). For example, the electric stimulation of the amygdala elicits a pattern of behaviors that mimic natural or conditioned fear (Rosen and Davis, 1990) and anxiolytic effects can be induced by applying benzodiazepine receptor antagonists in the basolateral nucleus (Hart et al., 2010). Distinct groups of neurons respond to (1) primary (unlearned) reinforcers like taste, (2) visual stimuli which previously have been paired with primary reinforcers like taste (positively reinforcing effects produced by novel stimuli independent of previous association of visual stimuli with primary reinforcement), and (3) responds to faces. Substantial evidence from both infra-human and human subject studies suggest that the sympathetic nervous system and the amygdala are crucial for the modulation of long-term memory storage for emotionally arousing events. The central and basolateral amygdaloid nuclei are involved in conditioned taste aversion, a unique type of innately predisposed (prepared) learning, in which the subject associates a taste with malaise over long delays (Yamamoto and Ueji, 2011). The consolidation of long-term explicit/declarative memory can be blocked after lesioning the basolateral amygdaloid nucleus (Gale et al., 2004).

Dysfunctions of the amygdala in humans have been reported in different neurological and psychiatric diseases like Klüver-Bucy-syndrome, Urbach-Wiethe-syndrome, temporal lobe epilepsy, Alzheimer disease, dementing disorders, and major depressive disorder. Generally, changes in emotionality and modulation of memory are attributed to lesions of the amygdala. A deficit in identifying fear in facial expressions of emotion while other expressions are intact was described by Adolphs et al. (1994) and Young et al. (1995). The importance of the amygdala for fear conditioning as well as promotion of effective memories (e.g., aversive conditioning) by arousing situations have been found in rodents (Chen et al., 2011, 2012). Neurotoxic lesion of the amygdala alters affective responses (Meunier et al., 1999). After lesion of the basolateral amygdaloid nucleus it was shown that the antidepressant fluoxetine has a positive effect on hippocampal cell survival (Castro et al., 2010). This emphasizes the fact that the amygdala may modulate the antidepressant action in hippocampal neurogenesis and its relation to depression-like behaviors. Furthermore, the basolateral amygdaloid nucleus seems to play an important role in the representation of the sensory features of motivationally significant events. Dwyer and Killcross (2006) demonstrated that animals with lesions of the basolateral amygdaloid nucleus can represent the sensory aspects of neutral events but not the sensory aspects of motivational events.

The connectional role of these regions which have been distinguished functionally can be determined by investigating the connectome of the amygdala in terms of graph theoretic and multivariate analysis. This contribution aims to establish and elucidate the connectome of the rat amygdala with a focus on intrinsic connectivity and emphasis on central, medial, lateral, and basal nuclear complexes.

2. Materials and Methods

Organization, visualization, and analysis of connectome data have been performed in neuroVIISAS (Schmitt and Eipert, 2012), a platform-independent generic framework developed in JAVA™. The neuroVIISAS installation package can be downloaded from All computations and visualization have been realized within the neuroVIISAS framework. The following pipeline describes the connectome analysis performed in this study:

1. Search, selection, and evaluation of tract-tracing publications

2. Defining the network of interest by selecting regions

3. Computing matrix representations

4. Global parameter computation

5. Local parameter computation

6. Motif analysis

7. Multivariate analysis

8. Vulnerability analysis

9. Extrinsic connectivity estimation

2.1. Nomenclature

The superficial cortical-like nuclear group consists of the anterior, posterolateral, posteromedial cortical amygdaloid nuclei, the amygdalohippocampal and amygdalopiriform transition area, the nucleus of the lateral olfactory tract, the nucleus of the accessory olfactory tract, the medial amygdaloid nucleus, and the anterior amygdaloid area. The nuclear groups that are located in the deep zones of the amygdaloid body are the ventral basolateral nucleus, the central amygdaloid nucleus, and the laterobasal nuclear complex as well as three divisions, the lateral, basolateral and basomedial nuclei, intercalated masses, amygdalostriatal zone, granular and parvicellular interfascicular islands, and intramedullary griseum. Additional groups are the extended amygdala (central and medial sublenticular extended amygdala, paraseptal lateral, and medial bed nuclei of the stria terminalis) and the unclassified cell group (part of “Other amygdala areas” in Table 1). This grouping of major nuclei and primary subdivisions is based on the work of de Olmos et al. (2004) and it is visualized in a hierarchy exported from neuroVIISAS (Figure 2). The hierarchical nomenclature of de Olmos et al. (2004) is based on cyto- and fibroarchitectonical criteria, ontogenetic, histo- and immunocytochemical, and hodological data. However, in the tract-tracing literature of the amygdala, much smaller subparts of nuclei and zones are described. In order to preserve this high-resolution spatial information, further subdivisions based on the nomenclature of de Olmos et al. (2004) have been taken into account. The de Olmos et al. (2004) nomenclature has been related to further (Paxinos and Watson, 1986; Swanson, 1992; Pitkänen, 2000; Price, 1981) terminologies (Table 1). In Table 1 major nuclear groups such as LA, Bnc, AB, Me, Ce, and AHi possess several subdivisions and are located at the 4th level of the amygdala hierarchy. The basal nucleus possesses a maximum of 89 subdivisions that were found in the amygdala tract-tracing literature. In the following, we refer only to an extended version of the (de Olmos et al., 2004) nomenclature obtained by subdividing the regions mentioned in tract-tracing studies.


Table 1. Comparison of amygdala nomenclatures.


Figure 2. A hierarchical nomenclature of subregions of the rat amygdala based on (de Olmos et al., 2004) has been implemented for the rat connectome project in neuroVIISAS. Only the top levels of the hierarchy are shown. Different intensities of magenta indicate different subregions.

2.2. Data Collection

Data retrieval was based on publications reporting tract-tracing studies of intrinsic and extrinsic amygdala connections using PubMed1, describing:

1. tract-tracing experiments using anterograde, retrograde, or bidirectional (where location of labeled perikaryons and axonal terminals are explicitly described) transport of substances

2. juvenile or adult rats

3. healthy, genetically un-altered and untreated control rats

4. results in English, German, or French

The search expression [amyg* AND rat AND brain AND (pathway* OR projection* OR afferent* OR efferent* OR connect*) AND (trace* OR tracing OR retrog* OR anterog*)] NOT (virus OR viral) yielded 487 papers, which were added to the connectome bibliography references.bib in JabRef2. This bibliographic information is also available online at the neuroVIISAS webpage. However, the PubMed search as described consider titles and abstracts of publications only. More connectivity data of the amygdala have been collected by evaluating tract-tracing publications of cortical areas and basal ganglia. At present, the whole rat connectome project contains data from 2100 tract-tracing publications. Papers not conforming with the criteria were marked in references.bib in order to exclude them automatically from search results. Intrinsic connections of the amygdala are described in 81, extrinsic input to the amygdala was found in 332, and extrinsic output from the amygdala in 337 publications.

The connectional information was gathered from these publications by the following criteria:

1. anterogradely labeled terminals and retrogradely marked perikaryons are clearly related to regions

2. locations of injection sites are described unambiguously

3. injection sites do not overlap with adjacent regions or fiber bundles

4. lesion and transsynaptic-tracing studies are excluded

The following data were retrieved from the publications and integrated in the rat connectome project of neuroVIISAS:

1. ipsilateral, contralateral, and bilateral connections

2. semiquantitative information that is often designated as the “weight” of a connection

3. tract-tracing substance

4. transport direction of tracer

5. case or individual rat where a connection has been observed

6. bibliographic details (stored in bibtex format) related to each connection

7. links to the pdf-documents of the tract-tracing publications

All information about connections is stored as objects in a zip-compressed project file that can be exported in spreadsheet format as well as in the owl ontology format of Protégé (Cimino and Zhu, 2006; Zhang et al., 2006; Poliakov et al., 2007; Musen et al., 2009). Before connectivity data were imported into neuroVIISAS, correctness was checked by at least two of the authors. The results of independent retrograde and anterograde experiments were integrated, such that the source of a connection, the target, its weight, and further connectivity attributes are available in a network representation that can be analyzed and visualized. All amygdala connections were added to the rat connectome data based on the hierarchy as shown in Figure 2. The interactive navigation, analysis, and visualization functions of neuroVIISAS allow a consistent and complete multi-level inspection of connectomes.

The subdivision of de Olmos et al. (2004) was used as an initial hierarchical nomenclature (Figure 2). This subdivision of regions of the amygdala is chiefly based upon cytoarchitectonic differences. However, it was extended based on descriptive subdivisions introduced in tract-tracing studies. The hierarchy consists of nodes (neuroanatomic regions) and relations whereby the relations in terms of an ontology belong to the type “is part.” Within the hierarchy, a root node (supranuclear division of the amygdala), nodes without further subdivisions (leafs), and nodes with further subdivisions exist. The sequence of hierarchically predefined leafs will be applied to matrix representations of connections and computations of further matrices to allow comparison. So far, 8 levels of subdivision of the amygdaloid complex are available in the connectome. Three regions have been mentioned in the literature that are possessing connections at the 8th level of subdivision of the amygdala: posterior basolateral nucleus parvicellular subdivision caudal part (Groenewegen et al., 1997), posterior basomedial nucleus medial half in the caudal part, and posterior basomedial nucleus lateral half in the rostral part (Kishi et al., 2006). The left- and right-hemispheric amygdala are connected by contralateral projections. The intrinsic connections of the left and right amygdala are the same. Since we have not found a lateralization of the amygdala connectivity, results of the left amygdala are presented.

The amygdala connectome is a partial connectome within the total connectome of the whole rat nervous system. Different ranges of colors are assigned to basic parts of the central nervous system, e.g., diencephalon, midbrain, cerebral cortex, oblongate medulla. Within this scheme of colors the amygdala subregions are assigned to shades of magenta to differentiate them from non-amygdala adjacent regions. The same colors of regions were consistently used in the hierarchy, 3D-visualization (Figure 1), matrices (Figure 4), and principal component analysis (Figure 9).

3. Results

In addition to the definitions, descriptions and interpretations of parameters and statistical procedures that are used in the results subsection, the supplement contains a complete list of definitions.

3.1. Global Parameters of the Intrinsic Amygdala Network

The hierarchy of regions of the left amygdala consists of 362 nodes. 279 of them do not posses further subregions. For 162 of these 279 regions of the left amygdala no intrinsic connections are documented. These 162 regions and, if needed, their parent regions were successively removed in order to build a network where each node has at least one connection. The resulting network contains 132 regions that are forming a single component connected by 665 edges (Figure 3A). A multi-level representation of the adjacency matrix is shown in Figure 4A. The final condensed intrinsic network contains only regions which have at least one input and one output (Figures 4B,C). This network has the benefit of being well-defined because regions where no connectivity data are available do not contribute to transmission of information in the network and should not influence global and local network parameters, limiting the value of any comparison with randomized networks. The condensed network contains 49 regions that are connected by 464 edges (Figure 3B). Both figures show the global network parameters and results of simulations [Erdös Rényi, Watts-Strogatz, Barabasi-Albert, Eipert (Eipert network), modified (directed graph) OHO (Ozik et al., 2004), and rewiring network generators (for details, see Schmitt and Eipert, 2012)] with the same number of nodes and edges as the real intrinsic left amygdala network. Simulations are used in order to compare the quantitative features of the real network with averaged parameters of different types of random networks that are simulated several times (see below). In the condensed network, the mean number of edges per node is 18.939 (average valency) and the line density (connections of the network divided by the number of all possible connections) 19.728%. The heterogeneity (standard deviation of valency divided by the average value) of 0.68 suggests that each node has a sum of inputs and outputs in the same range (Estrada, 2010). On average each node can be reached by using 2.478 edges (average path length). The average, number of links incident upon a node (centrality) is 0.438. A measure of how complete the neighborhood of a node is, is given by the average cluster coefficient (Rubinov and Sporns, 2010) and turns out be 0.505. Both global parameters are larger than the mean of 1000 Erdös Rényi simulations with the same number of nodes and edges hinting at a specific wiring structure. However, these parameters are more similar for the EN, modified OHO, and rewiring networks. In a small-world network most nodes are not neighbors of one another, however, most of them can be reached from every other node by a small number of edges. The small-worldness parameter (cluster coefficient of real and Erdös Rényi networks divided by average path lengths of real and Erdös Rényi networks) of 1.99 is the largest, compared to simulated networks (Humphries et al., 2006; Humphries and Gurney, 2008). Hence, the intrinsic amygdala network has a small-world property. Interestingly, the average path length of 2.478 is the largest compared to all simulations. A small average path length and a large average cluster coefficient hint at a small-world network, consistent with the small-worldness parameter. The condensed amygdala network has a relatively small delta error (mean deviation between the data and the approximation of the power law distribution multiplied by 100) of 1.4 in comparison with a power law distribution, pointing to a scale-free property. The modularity Newman (2004) of 0.36 is relatively small, indicating dense connections between modules of the intrinsic amygdala network (Figures 3C,D).


Figure 3. Global analysis of the left intrinsic amygdala network applied to (A) all regions (132 regions and 665 edges) and (B) 49 regions that have at least one input and one output that are connected by 464 edges (condensed intrinsic amygdala network). The selected regions of the network are visualized by hierarchies of triangles (top panels). The bottom panels display the global network statistics. Both networks are simulated 1000 times by 6 different network generators and global parameters of 1000 simulations were averaged. The result of a modularity analysis is shown in (C) (all regions) and (D) (regions with at least one input and output). (C) Within the 5 modules the number of connections is larger than in between the modules. (D) The number of modules is reduced to 3 if only those regions are considered that have at least one input and one output (90° rotated view).


Figure 4. Hierarchical connectomes allow the selection of different levels of subdivisions of regions. (A) This multi-level representation is shown in the panel with gray background and applied to the left amygdala with 132 regions from level 3 to level 8. Sources or efferent regions are arranged in rows and target or afferent regions are arranged in columns. The sequence of regions is predefined by the sequence of leafs of the defined region hierarchy (see Materials and Methods). (B) Triangle representation of all regions that have at least one input and one output at level 8 and (C) the corresponding adjacency matrix describing 49 regions and 464 connections.

3.2. Network

The rows of the adjacency matrix indicate sources, perikaryon locations or efferents, and the columns are indicating targets, axonal terminal locations, or afferents. The adjacency matrix shows a cluster of dense interconnections of the subregions of the superficial amygdaloid cortex (superficial cortical-like nuclear group; Figure 5A). The different subdivisions of the amygdaloid nucleus are subregions of the superficial amygdaloid cortex and possess numerous outputs (103) to anterior basomedial and accessory basal nuclei. The cluster of strong local connectivity within subregions of the superficial amygdaloid cortex can be found again in the distance matrix where most of these connections turn out to use only one edge (Figure 5B) which means that most of these regions are directly interconnected. The ventromedial part of the lateral nucleus, the bed nucleus of the stria terminalis lateral division ventral part, and the bed nucleus of the stria terminalis medial division posterolateral part attract attention due to many long distances indicated by dark gray values in the correspondent rows. The bed nucleus of the stria terminalis medial division posterolateral part and the ventromedial part of the lateral nucleus show low communicability values (Estrada and Hatano, 2008; only few shortest paths between a pair of connected regions exist), whereby the bed nucleus of the stria terminalis lateral division ventral part has similar values as neighboring nuclei (Figure 5B). The subregions of the superficial amygdaloid cortex are forming a cluster with large communicability values, hence, many shortest paths between pairs of connected regions exist. The matrix of connectivity in degree matchings (Sporns, 2002) shows a new set of regions that have quite similar input patterns of connectivity. These subregions belong to the bed nucleus of the stria terminalis and are located within the upper left part of the matrix indicated by a yellow square (Figure 5F). The subregions of the bed nucleus of the stria terminalis have lower connectivity output matching values (Sporns, 2002; Figure 5E) than input values which indicates more similarity of input connections to than output connections from BST-subregions. In the connectivity input and output matching matrix, the larger values of the subregions of the bed nucleus of the stria terminalis are highlighted again and the second cluster of the superficial amygdaloid nucleus is also visible (Figure 5D). The large input connectivity matching means that pairs of regions have similar afferents. Because the same regions have low output connectivity matching values, processed information of them is distributed divergently to target regions.


Figure 5. Connectivity matrices of the condensed intrinsic left amygdala network. (A) Adjacency matrix with highlighted regions (transparent magenta) of the superficial amygdaloid cortex. Ten connectivity weights are color-coded. (B) Distance matrix. Gray levels 1–7 indicate the smallest number of edges between regions. Light gray values indicate short distances. (C) Communicability matrix. The gray scale codes values between 0.0553 and 7.915E10. Large values indicate (e.g., ACo) that many short shortest paths between a pair of regions exist. (D) Connectivity matching matrix for inputs and outputs. The gray level scale is coding values between 0 and 1. Large values indicate similar input and output connections of two regions. (E) Connectivity matching matrix for inputs. Red shades are coding similar inputs of two regions. (F) Connectivity matching matrix for outputs. Green shades are coding the similar outputs of two regions. The yellow square indicates the regions of the bed nucleus of the stria terminalis.

3.3. Local Network Parameters

Local network analysis using the same set of regions as in the global network analysis reveals the impact of single regions with regard to connectivity and flow of information. Most of the 30 local network parameters calculated by neuroVIISAS (Table 2) are strongly correlated with the degree all parameter [the degree of a node (vertex) is the number of incoming (afferent) and outgoing (efferent) edges (connections)]. If incoming connections are differentiated from outgoing connections, then the terms indegrees (in) and outdegrees (out) are used (Sporns 2002). The anterior cortical amygdaloid nucleus has the largest degree all (55) followed by the posterior amygdaloid nucleus (41) and the bed nucleus of the accessory olfactory tract (33). However, the Shapley rate (a low Shapley rate indicates a large impact of a node in a network; Kötter et al., 2007) of the posterior amygdaloid nucleus is smaller (−2.532) than that of the anterior cortical amygdaloid nucleus (−0.759). Interestingly, the Shapley rate of the supracapsular bed nucleus of the stria terminalis lateral part is the second smallest (−0.91), although only having 13 inputs and outputs. The Shapley rates are visualized in the 3D-atlas using a multiaxes extension (Figure 6).


Table 2. Local parameters of the condensed amygdala network sorted by the degree all parameter.


Figure 6. Visualization of Shapley rates in 3D-expansion view of the unilateral and bilateral amygdala. (A) View from dorsal. The left- and right-hemispheric amygdala are visible through the transparent pars cranialis of the central nervous system. (B) Unilateral amygdala with regions that have the lowest Shapley rates (see text) and color-coded weights of connections (color-codes of connections are the same as in Figure 5). Color-coded Shapley rates are assigned to regions. (C) Bilateral amygdala with ipsilateral and contralateral connections.

Most cyclic pathways from the anterior cortical amygdaloid nucleus back to it are possible via 21 different 2-edge-cycles, followed by the intercalated nuclei (16 different 2-edge-cycles) and the accessory basal nucleus parvicellular division, periamygdaloid complex sulcal division, and accessory basal nucleus magnocellular division (15 different 2-edge-cycles).

As described in the introduction, the lateral (LA), the basolateral (BL), and the central amygdaloid nucleus (CE) are frequently analyzed in behavioral, neurophysiologic, and pharmacologic studies. In addition, the medial amygdaloid nucleus (ME) is considered as an important entity (Aggleton 2000; Pitkänen 2000) of amygdaloid circuits and has been analyzed in the context of LA, BL, and CE. In the hierarchy that is used here, the LA and BL are child nodes of the laterobasal nuclear complex, CE is a subregion of the central extended amygdala, and ME is a subregion of the medial extended amygdala. They are not directly visible because some of their subregions are expanded. However, by reducing some branches of the hierarchy, they become leafs of the hierarchy (26 regions, 350 connections; not shown here) and local parameters of these regions can be computed. The CE has the smallest Shapley rate and most inputs and outputs (indegree: 22, outdegree: 22) indicate its large impact of the intrinsic amygdala connectome. ME (outdegree: 20, indegree: 21) has the second lowest Shapley rate of −0.076. The Shapley rate of BL is −0.003 (outdegree: 14, indegree: 15) and that of LA is −0.029 (outdegree: 17, indegree: 17). With regard to the other 24 regions, ME has the 3rd, BL the 9th, and LA the 10th smallest Shapley value. Hence, their intrinsic impact of BL and LA are considerably smaller than those of CE and ME. In the network where CE, ME, LA, and BL are leafs, most cyclic pathways are passing CE (2 node cycles: 19 passages, 3-node passages: 254) and ME (2 node cycles: 18 passages, 3-node passages: 247; LA 2 node cycles: 15, 3-node cycles: 198; BL 2 node cycles: 10, 3-node cycles: 121).

The sum of extrinsic contralateral and ipsilateral input of CE is 506 followed by BL (286), LA (230), and ME (216). The extrinsic output of CE is 380, of BL is 174, of LA is 109, and ME has the 2nd largest: 207. These four regions receive most extrinsic inputs among all regions of the 26 node hierarchy.

3.4. Motifs and Circuits

Motif analysis for directed 3-node and 4-node motifs revealed a significantly more frequent expression of completely reciprocal motifs (all edges of a motif are reciprocal) than in a random network (Figure 7A). A maximum of 13 directed motifs that consist of 3-nodes and 199 that consist of 4 nodes can be constructed. The frequency of all motifs were determined by a paralleled (using n cores in parallel) isomorphism search (Schmitt 2012). The connections were randomized 1000 times (1000 networks with the same number of nodes and edges as the condensed intrinsic amygdala network were computed) by the rewiring simulation, the motif frequencies were determined and visualized by black dots in the motif-diagram. Blue dots that are above or beyond the black point clouds indicate a significant larger or lower frequency of occurrence in the real amygdala network in comparison to rewired networks. The first 6 motifs can be considered as variations of convergent and divergent motifs. These are all less frequent in the real than in the rewired networks. The cyclic motifs (3-07, 3-10) are also less frequent in the real amygdala network (Table 3). The completely reciprocal motif (3-13) is much more abundant in the real than in rewired networks (Figure 7A, arrow). The 4-node subgraph analysis shows similar results as the 3-node subgraph analysis (Figure 7B). Variations of divergent and convergent 4-node motifs (data not shown) as well as the completely reciprocal motifs are much more frequent in the rewired than in the real amygdala network. In Figure 7B only those motifs are shown for which p = 0 (p = 0 indicates motifs that are significant more frequent in a real network, p = 1 indicates motifs that are significant more frequent in a random network; Table 4). In the motif 4-116 the motif 3-13 is included and it occurs most frequently.


Figure 7. Motif analysis of the condensed intrinsic amygdala network. Frequency of motifs in the intrinsic amygdala network is indicated by a blue dot. Frequencies of motifs in rewiring randomizations are indicated by black points. Red bars show the standard deviation of motif frequencies. (A) 3-Node directed motif analysis of 13 subgraphs. The abundant fully reciprocal motif (e.g., 3-13, where 3 is the number of nodes of this motif and 13 is the identification number of a this motif) is indicated by an arrow and magnified. (B) 4-Node directed motif analysis. Those subgraphs that were significantly more frequent in the real network are shown.


Table 3. All 13 3-node motifs and resampling calculations based on 1000 rewiring randomizations of 49 regions and 464 connections.


Table 4. 13 Significant 4-node motifs from a total of 199 4-node motifs and 1000 rewiring randomizations of 49 regions and 464 connections.

The particular regions that are part of the completely reciprocal 3-13 motif are the anterior cortical amygdaloid nucleus (90 times), the intercalated nuclei of the amygdala (90), the periamygdaloid complex sulcal division (84), and accessory basal nuclei (83). The central division of the sublenticular extended amygdala (3), the ventral basolateral nucleus (2), and the posterior basolateral nucleus lateral part (2) occur in the cyclic 3-07 motif with almost no reciprocal connections. The frequencies of regions in the completely reciprocal 4-199 motif is similar to the case of motif 3-13: anterior cortical amygdaloid nucleus (258), intercalated nuclei of the amygdala (268), periamygdaloid complex sulcal division (254), and accessory basal nuclei (254). In contrast to the motif 3-13 region, participation frequencies of the intercalated nuclei of the amygdala are higher than the anterior cortical amygdaloid nucleus for motif 4-199.

The four regions LA, BL, CE, and ME have been analyzed with regard to local parameters and extrinsic connectivity (see above). The motif-search revealed that ME is part of the 3-13 motif in 90 cases, CE (86), LA (69), and BL (19). Interestingly, BL is most frequently a member in the 3-04 motif with one reciprocal edge and a single directed edge (CE: 31, LA: 23, ME: 20). The motif 3-09 with two reciprocal edges has another frequency distribution with regard to these regions (CE: 47, ME: 40, LA: 32, BL: 25).

LA receives sensory inputs (Kim and Jung, 2006; Ciocchi et al., 2010; Haubensak et al., 2010; Tye and Deisseroth, 2012). Some of these sensory inputs originate in the somatosensory, auditory, and visual cortical regions. Typically LA projects to the basal nuclear complex which projects to CE (Ciocchi et al., 2010; Haubensak et al., 2010; Tye and Deisseroth, 2012). In terms of behavior this pathway is interpreted as the fear conditioning pathway (Kim and Jung, 2006; Ciocchi et al., 2010; Haubensak et al., 2010; Tye and Deisseroth, 2012). The rat connectome was analyzed by a constrained search using the pathway analysis approach of neuroVIISAS to detect pathways through LA → Bnc → CE. CE seems to be necessary for the fear response. As shown in Figure 8, the pathways from somatosensory, visual, and auditory regions to LA, from LA to Bnc and from Bnc to the output complex, CE can be reconstructed using a transhierarchical search strategy. The subregions of each of the three complexes are densely interconnected by differentially weighted connections.


Figure 8. Pathway analysis of somatosensory, visual and auditory regions of the cerebral cortex and the lateral, basal, and central nuclear complexes of the amygdala. (A) Somatosensory, (B) visual, and (C) auditory cortical regions are defined as source regions of the projection to the lateral (LA) nuclear complex of the amygdala. All subregions of the LA are taken into account as putative target regions. From these LA-target regions connections to the basal nuclear complex (including BL) where searched. From the basal nuclear complex, connections were chosen which project to CE subregions. The colors of lines corresponds to the weight code in Figure 4. Thick lines indicate reciprocal connections and the numbers indicate the number of studies that have documented a particular connection. In the left part are the source regions displayed, then the LA subtree regions are following, then the basal nuclear complex subtree regions, and on the right the subregions of the CE subtree. The numbers at the left part of the pathway-diagrams are indicating hierarchical levels.

3.5. Principal Component Analysis

Principal component analysis (PCA) of the condensed intrinsic amygdala network was applied to detect groups of regions with similar patterns of connectivity (Echtermeyer et al., 2011; Schmitt and Eipert, 2012). The degree all (DGall; all refers to the sum of input and output connections), average degree of neighbors (AvgDGnb; Rubinov and Sporns, 2010), cluster coefficient all (CluCall; Rubinov and Sporns, 2010), cluster coefficient of 2nd neighbors (Echtermeyer et al., 2011; CluC2), variation coefficient of neighbor degree (VCDG1; Echtermeyer et al., 2011), and locality (Loc; da Costa et al., 2006) were used as 6 dimensions (feature vector) for PCA because these local parameters quantify connectivity relations of direct and indirect neighbors. Here, PCA maps feature vectors of all nodes to a plane to detect differences between nodes with regard to their feature expression. A Parzen-window (Echtermeyer et al., 2011) was used to realize a probability density function in the two-dimensional PCA-plane. The contribution of each parameter to the components is shown in Table 5. The parameters locality (Loc), cluster coefficient all (CluCall), and average neighbor degree (AvgDGnb) are the largest among all 6 variables and have the biggest influence on the x-axis (component 1). The absolute values of parameters degree all (DGall) and the variation coefficient of neighbor degree (VCDG) have the biggest influence on the y-axis (component 2; Figure 9A). The subventricular nucleus (SV) is located at the upper left corner of the PCA-plane. SV has only three direct neighbors which have a similar amount of connections to 2nd neighbors of SV like the 2nd neighbors among themselves (1st weak connectivity). The anterior cortical amygdaloid nucleus (ACo) has a large degree all and many connections between 1st neighbors, however, sparse connections from 1st to 2nd neighbors and among 2nd neighbors (1st NB focused connectivity). The amygdalohippocampal area medial part (AHimp) has many neighbors which are very densely connected among themselves and more densely than the 2nd neighbors among themselves (1st dense connectivity). A similar connectional pattern can be found for the regions mapped around AHimp: Pacs, Bmc, AHilp, I, MeC, PACm, MeRo, ABm, MeCVABp. The posteromedial cortical nucleus (PMCo) is located almost in the center of the PCA-plane. It has a medium number of connections to direct neighbors and the 2nd neighbors are strongly connected among themselves. In relation to the 1st neighbors, the 2nd neighbors of PMCo are connected more densely (1st balanced mid connectivity). Around PMCo is a second larger group of regions that join the latter connectivity feature of PMCo. The bed nuclei of the stria terminalis are among this PMCo-group. The amygdalostriatal transition area rostrocaudal part (AStrrc) is located at the bottom of the PCA-plane and has only two 1st neighbors. The 2nd neighbors of AStrrc are strongly connected among themselves (2nd NB focused connectivity).


Table 5. A feature vector consists of the parameters degree all (DGall), average degree (AvgDGnb), cluster coefficient all (CluCall), cluster coefficient of 2nd neighbors (CluC2), variation coefficient of neighbor degree (VCDG), and locality (Loc; for a definition of parameters, see Echtermeyer et al., 2011).


Figure 9. PCA analysis of the condensed intrinsic amygdala network. (A) Probability density diagram of the 2-dimensional plane of PCA showing the 49 mapped regions. Component axes are indicated and correspond to Table 5. Yellow letters refer to the circle-diagrams. All circle-diagrams have a center circle that correspond to a region of interest of the PCA-plane in (A). The inner-ring of circles are the 1st neighbors and the outer-ring of circles are the 2nd neighbors of the center region. (B) Subventricular nucleus (SV). (C) Anterior cortical nucleus (ACo). (D) Amygdalohippocampal area medial part (AHimp). (E) Posteromedial cortical nucleus (PMCo). (F) Amygdalostriatal transition area rostrocaudal part.

3.6. Vulnerability of the Intrinsic Amygdala Network

Vulnerability analysis in neuroVIISAS was applied to identify nodes and connections that contribute significantly to the integrity of the intrinsic amygdala network. The analysis has been performed by iteratively removing nodes (significance list Table 6) or edges (vulnerability matrix Figure 10) from the condensed intrinsic amygdala network, followed by calculating the closeness (inverse of the sum of distances to all other nodes) of remaining nodes. The relative change in percent of the closeness turns out to be a meaningful quantity for the significance of removed nodes. If the relative average closeness is large (average path length increases) after removing a node, then this node possesses a large significance for the connectivity of the network. A negative relative average change of closeness (average path length decreases) may occur if the removed node is relatively far away from most other nodes. In Table 6 the posterior amygdaloid nucleus (9.434%), the anterior cortical amygdaloid nucleus (9.227%), and the anterior basomedial nucleus (2.695%) have largest significances and if these nodes are removed the average path length increases much more than removing other nodes. The relative average change of closeness after removing connections is visualized in the vulnerability matrix (Figure 10). The cluster of regions in the lower right quadrant of the vulnerability matrix shows similar values and belongs to the superficial amygdaloid complex.


Table 6. Significance of regions in percent based on the change of closeness after removing a particular region from the network.


Figure 10. Vulnerability matrix of the condensed intrinsic amygdala network. Values indicate the decrease of closeness after removing a particular connection. Removal of the connection from the posterior basolateral nucleus lateral part to the supracapsular bed nucleus of the stria terminalis lateral part has the largest significance of 2.78% (lightest gray value).

3.7. Extrinsic Connectivity of the Amygdala

The condensed intrinsic network of an unilateral amygdala has ipsilateral and contralateral inputs and outputs to regions that do not belong to the intrinsic amygdala connectome. Direct connections start at leafs (regions that are not split into subregions or where subregions are not opened in the user-defined hierarchy presentation) of subtrees and were considered to calculate extrinsic inputs and outputs. The anterior cortical amygdaloid nucleus has the largest direct ipsilateral input (101) and the second largest direct contralateral input (11). The direct ipsilateral output consists of 138 and the direct contralateral output consists of 19 connections (rank 5). The basal nucleus intermediate division has the most directed contralateral outputs (37), yet, only 19 ipsilateral outputs (rank 29). The anterior cortical amygdaloid nucleus turns out to be a part of the condensed intrinsic amygdala network that has most ipsilateral inputs and outputs. Table 7 documents that ipsilateral and contralateral connectivity is not positively correlated for each region, e.g., the interstitial nucleus of the posterior limb of the anterior commissure medial part.


Table 7. Extrinsic bilateral inputs and outputs of all regions of the condensed intrinsic amygdala network.

4. Discussion

Analysis and visualization of neuronal connections with regard to networks and connectomes have been pushed forward in the last 30 years (MacDonald, 1983; Felleman and Essen, 1991; Tononi et al., 1994; Young et al., 1994). Specific connectivity analysis based on meta-studies have been performed on the nucleus of the solitary tract of the rat (Palombi et al., 2006), hippocampus of the rat (Burns and Young, 2000), thalamocortical connectivity (Scannell et al., 1999; da Costa et al., 2006), brainstem reticular formation (Humphries et al., 2006), and retrosplenial cortex (Sugar et al., 2011).

In an ongoing meta-study of tract-tracing publications of the rat nervous system, the connectivity data of 2100 articles are collected in neuroVIISAS, a generic platform for digital atlasing, connectivity analysis and visualization (Schmitt et al., 2012), and population based simulations. In all tract-tracing studies of the amygdala, groups of regions were described as subdivisions and parts of regions. Therefore, regions of the nervous system are arranged in a neuroontology within neuroVIISAS. Because most tract-tracing studies have been performed in the rat, these rat-specific articles have been evaluated.

The total number of peer-reviewed publications of tract-tracing studies of the rat, without considering other organisms, can be determined by convenient expressions that filter the PubMed database (see text footnote 1). 4528 Of such articles can be found in September 2012. 491 Articles of these 4528 mention amyg* (amygdaloideum, amygdaloid, etc.) in their abstracts. Hence, they can be considered as the core of the meta-analysis of tract-tracing based connections of the amygdala. In addition, the bibliographies of review articles and monographs (de Olmos et al., 2004) have been evaluated. In other projects, the connectomes of further functional systems (basal ganglia, cerebellum, cerebral cortex, hippocampus, spinal cord, peripheral nervous system) are analyzed. In many of these tract-tracing publications that do not use the term amyg* in their abstract, input and output connections to and from the amygdala are described and have been captured in the whole rat connectome project. Interestingly, the rat is the organism in which most peer-reviewed tract-tracing studies are available.

Comparable connectome studies of the rat are considering neurons (Arenas et al., 2008) or several areas (whole brain and/or whole CNS; Bohland et al., 2009). The painstaking work of connectivity evaluation of Burns (1997) is based on the terminology of Swanson and Petrovich (1998), ranging at a comparable granularity as the terminology of Paxinos and Watson (2007). The systematic evaluation and databasing of connectivities of the rat brain was proposed and extensively investigated by Burns (1997). This work has been refined by Burns and Cheng (2006), Burns et al. (2008), and Bota and Swanson (2008). So far, a connectome analysis based on a meta-study of tract-tracing publications of the rat amygdala has not been performed.

A principal difficulty of extracting source-target-weight data from tract-tracing publications is the definition of interlinked source and target regions in different studies. This is aggravated by the fact that ambiguities of region definitions may occur in different contexts:

1. spatial ambiguities

2. meaning and definition of region terms (ambiguity of definition)

3. location of a region term in a hierarchical terminology (ambiguity of classification)

4. comparative or interspecies alignment of terms (e.g., application of Brodmann area terms of human cerebral cortex parcelation to the rat cerebral cortex; ambiguity of homology)

These particular terminological contexts are the cause for the correspondence or concordance problem (Stephan et al., 2000; Bezgin et al., 2009; Bohland et al., 2009). An adaptive solution of the problem with regard to the development of neuroanatomical knowledge of regions of nervous systems of different species is not available. The objective relation transformation (ORT; Stephan et al., 2000) reduces the pattern of spatial overlap between region pairs to a discrete set of relational attributes containing identity, subset, superset, and partial overlapping. To realize ORT, a documentation of spatial features (e.g., extension, neighboring regions, super regions, subregions, development) of regions in tract-tracing publications is necessary. Unfortunately, the objective of tract-tracing studies is a description of sources, targets, and pathways visualized by axonal transport of tracing substances. Hence, a definition of regions and the comparison of them with other studies is out of scope in most cases (Pitkänen, 2000; Jones et al., 2005). Data that are necessary to apply ORT consistently are either missing in most tract-tracing studies or they must be inferred (indirect availability). The latter may also contribute to imponderables and variability in region definitions.

In this work, extracted source-target-weight data are directly derived from the original text to reduce inconsistencies as much as possible. This has the advantage of reproducibility of sources and targets in a connectome. Interpretations of how much a region partially overlaps with regions in other publications or if subregions are assigned to different superregions are avoided to prevent errors introduced by misinterpretations. Furthermore, necessary data that allow correct interpretation of terminological overlap and subdivision problems are often missing in publications. Hence, source, target, and weight data are extracted with a minimum of further interpretation to accumulate them continuously during evaluation of new tract-tracing publications. By using accumulated data, it is possible to generate different selections of nodes from the whole neuroontology of the rat nervous system. In neuroVIISAS, several mechanisms to assemble regions of the complete neuroontology have been implemented to generate problem-dependent adjacency matrices. A coarse selection of superregions (e.g., thalamus, hypothalamus, sensomotoric cortex, cerebellum) for connectivity analysis provides more validity with regard to region definitions (comparable with tract-tracing data judgment in ORT), however, less important and relevant connectivity information is the result. Nevertheless, the approach of accumulating raw data from original publications in a consistent framework avoids the introduction of errors while solving overlap and classification problems if necessary data are not available for the problem at hand.

Beside the concordance problem, there are further factors that constrain the informative value of tract-tracing based connectome construction:

1. different dynamics of tracer transport

2. selectivity of tracer transport (unidirectional, bidirectional)

3. identification of axonal terminals and fibers of passages

4. tracer uptake by damaged fibers of passage

5. technical details of tracer application (volume, velocity, iontophoresis)

6. distance of sources and targets from the location of tracer application and survival time

7. spread of tracer in adjacent regions (especially in stereotaxic pressure injections)

8. sensitivity of tracer and immunohistochemical visualization

9. semiquantitative and qualitative interpretation of connection strength (weight) is not normalized between different studies and most criteria for weights are not based on stereological and densitometric evaluations

10. investigator bias (interpretation of descriptive representations of tract-tracing results)

The amygdala connectome has been analyzed using neuroVIISAS by several global and network measures, matrices, motifs, and multivariate techniques. However, other methods are available to investigate the structures and features of complex networks (Kötter and Stephan, 2003; Tononi and Sporns, 2003; Sporns and Kötter, 2004; Sporns and Zwi, 2004; Bassett and Bullmore, 2006; Nagyessy et al., 2006; da Costa et al., 2007; Kötter et al., 2007; Sporns et al., 2007). More specifically, neural complexity (Tononi et al., 1994), community structure (Newman and Girvan, 2004), connectivity descriptors (Kammer and Täubig, 2004; Goryczka and Arodz, 2006), functional topology analysis (Blinder et al., 2005), efficiency and cost of networks (Achard and Bullmore, 2007), node conformity of factorizable networks (Dong and Horvath, 2007), lesion analysis, and transfer entropy (Schreiber, 2000; Honey and Sporns, 2008) as well as spectral analysis of networks (Baltz and Kliemann, 2004) are advanced methods to reveal structures and functions of networks. In addition, multivariate methods of non-metric multidimensional scaling and Procrustes R2 statistics (Goodhill et al., 1995; Burns and Young, 2000), hierarchical optimization (Hilgetag et al., 1996; Burns and Young, 2000), optimal hierarchical orderings (Hilgetag et al., 2000), and canonical variable analysis (da Costa et al., 2007) are options to further investigate results that are obtained by the principal component analysis as found here.

By selecting those regions where inputs and outputs within the intrinsic amygdala network are reported, a condensed intrinsic amygdala network was generated. This network is a scale-free and, to some degree, small-world network. The adjacency matrix and complexer matrices exhibit a cluster of connections (Figure 5) within the subregions of the superficial amygdaloid cortex which can also be observed at different levels of the region hierarchy (Figure 4). The medial amygdaloid nucleus has the largest number of efferents and afferents within this cluster, described here for the first time. The dense intraamygdaloid connection of the medial amygdaloid nucleus and the central amygdaloid nucleus that have also been reported in the review of Pitkänen (2000) can be confirmed here (subsection 4.2. Network, Figure 5). These densely interlinked regions of the superficial amygdaloid cortex belong to olfaction, autonomic hypothalamus and prefrontal cortex cognition.

The CE, ME, BL, and LA regions are considered as hodologic and functional entities of the amygdala (Aggleton, 2000; Pitkänen, 2000). Based on the subdivision of de Olmos et al. (2004), the resulting connectome was modified to obtain CE, ME, BL, and LA regions as connectional nodes. CE and ME have comparable quantitative features with regard to local parameters (indegrees, outdegrees, Shapley rates) and cycle counts. This indicates that these regions are very important with regard to integrity of the intrinsic amygdala connectome. In addition, their frequency of participation for the complete reciprocal motif 3-13 is larger in comparison to other regions. In contrast, BL and LA have lower Shapley rates and they do not participate as often as CE and ME to 2 cycle and 3 cycle pathways within the intrinsic connectome of the amygdala. However, BL receives the second most extrinsic inputs from non-amygdala regions of the ipsi- and contralateral hemisphere. From the connectional perspective this corresponds with the point of view that BL and LA are entries of sensory information (Kim and Jung, 2006). Furthermore, BL participates most often in the 3-04 motif with one reciprocal and a single non-reciprocal edge. The reciprocal edge can be interpreted in terms of a local regulatory function (positive or negative feedback) and the connection to a node that does not project to the reciprocally connected nodes could be regarded as a regulatory output of this motif. To summarize a connectional role of BL, we can consider it as an extrinsic receiver that may regulate extrinsic inputs and then distribute them in the intrinsic connectome of the amygdala. ME possesses second most extrinsic outputs and a relatively small Shapley rate which could emphasize the connectional role as a sender to extrinsic targets and a strong intrinsic network entity. However, CE receives and sends most extrinsic inputs and outputs and has the smallest Shapley rate which points to a more integrative function in comparison to BL and ME. These interpretations are speculative and need to be verified. Nevertheless, they are based on a huge number of observations in tract-tracing studies.

The sensory inputs from the somatosensory, visual, and auditory cortex to the LA have been confirmed in the connectome data. Furthermore, we can agree with dense interconnections between LA and basal complex regions. Reciprocal connections between LA and Bnc (BL is a subregion of Bnc) were also found. Dense interconnections with reciprocal projections between Bnc and CE are found as well. Hence, the pathway of fear conditioning (Kim and Jung, 2006; Ciocchi et al., 2010; Haubensak et al., 2010; Tye and Deisseroth, 2012) from different sensory regions can be reconstructed. However, it turns out that there are a lot more interconnections between specific subregions and reciprocal projections are known than considered in Haubensak et al. (2010), Ciocchi et al. (2010), Tye and Deisseroth (2012), and Kim and Jung (2006).

The Shapley rate is an indicator for the importance of a particular region in the network. Regions which posses largest importance with regard to the Shapley rate are the posterior amygdaloid nucleus, the anterior cortical amygdaloid nucleus, and the supracapsular bed nucleus of the stria terminalis lateral part. However, the latter has only a small degree all of 13 connections. Hence, Shapley rates should always be judged in comparison with other local parameters. Also, the Shapley rates of the bilateral condensed intrinsic amygdala network indicate that these regions are as important as those in the unilateral network.

The complete reciprocal 3-node and 4-node motifs occur significantly more often in the real condensed amygdala network than in rewiring simulations, a finding reported here for the first time. However, circular motifs without reciprocal connections are less abundant in the real network. Nodes that possess large centrality values and importance for the network structure are always involved in the completely reciprocal motifs. This could indicate a strong regulatory role and self control by connectional dense feedbacks.

By removing the posterior amygdaloid nucleus, the anterior cortical amygdaloid nucleus, and the anterior basomedial nucleus, strongest changes of relative closeness have been observed. However, the anterior basomedial nucleus does not have such strong centrality features like the regions with lowest Shapley rates, indicating that vulnerability analysis may provide additional information regarding the role of regions that has not been revealed by quantities of local network analysis.

In the PCA we found that the anterior cortical amygdaloid nucleus has many 1st neighbors, however, connectivity between 1st and 2nd neighbors is weak. In contrast, the amygdalohippocampal area medial part has many 1st neighbors and they are also strongly interconnected. The posteromedial cortical amygdaloid nucleus has numerous 1st neighbors which are strongly interconnected. In addition, the 2nd neighbors of PMCo are densely interconnected, too. These features of local connectivity of PMCo in relation to its smallest Shapley rate of the intrinsic amygdala connectome indicate its important role in the network.

Finally, it should be emphasized that the connectional information accumulated from the evaluation of a large amount of tract-tracing publications may only reflect partial trends in connectivity detection of the tract-tracing community. The real and complete connectome of the rat amygdala might have features that cannot be detected by analysis of all published tract-tracing data. However, by also including publications in this work which are not focused on tracer injections in the amygdala, validity of the basic network structure is expected to become larger. Further limitations of neural connectome analysis are the spatiotemporal dynamics of connectivity at the synaptic level.

In conclusion, the quantitative analysis of the amygdala connectome allows to identify regions of the superficial amygdaloid complex that are densely interlinked and may be important for the internal regulation of information processing of the intrinsic amygdala network. The amygdala network has small-word and scale-free properties. It contains complete reciprocal motifs significantly more frequently than in randomizations. The posterior amygdaloid nucleus, anterior cortical amygdaloid nucleus, and bed nucleus of the accessory olfactory tract turned out to send and receive most connections within the condensed amygdaloid network.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


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, and Hannah Ormanns for their work on extending the database and mappings. All work was supported by the Faculty of Mathematics and Natural Sciences and the Faculty of Medicine of the University of Rostock.


AA, AAA, anterior amygdaloid area; AAD, dorsal anterior amygdaloid area; AB, accessory basal nucleus; AB, accessory basal nucleus modified from de Olmos et al. (2004); ABm, accessory basal nucleus magnocellular division; ABmc, accessory basal nucleus magnocellular division; ABp, accessory basal nucleus parvicellular division; ABpc, accessory basal nucleus parvicellular division; ACo, layers 1-2 = anterior cortical nucleus; ACo, anterior cortical amygdaloid nucleus; AHA, amygdalohippocampal area; AHAl, amygdalohippocampal area lateral division; AHAm, amygdalohippocampal area medial division; AHi, amygdalohippocampal area; AHiAL, amygdalohippocampal area anterolateral; AHilp, Amygdalohippocampal area lateral part; AHilp, amygdalohippocampal area lateral part; AHimp, amygdalohippocampal area medial part; AHimp, amygdalohippocampal area medial part; AHiPM, amygdalohippocampal area, posteromedial; APir PM, amygdalopiriform transition area posteromedial part; AStrrc, amygdalostriatal transition area rostrocaudal part; B, basal nucleus; BA, bed nucleus accessory olfactory tract; BAOT, bed nucleus of the accessory olfactory tract; BAOT*, bed nucleus of accessory olfactory tract; Bi, basal nucleus intermediate division; BL, BLA, basolateral nucleus; BLA*, anterior basolateral nucleus; BLAa, basolateral nucleus, anterior part; BLAp, basolateral nucleus, posterior part; BLPLP, posterior basolateral nucleus lateral part; BLV, ventral basolateral nucleus; BM, BMA, basomedial nucleus; BMA, layer 3 = Anterior basomedial nucleus; BMA, anterior basomedial nucleus; BMAa, layer 3 = Basomedial nucleus, anterior part; BMAp, basomedial nucleus, posterior part; Bmc, basal nucleus magnocellular division; Bmc, basal nucleus magnocellular part; BMP, posterior basomedial nucleus; BMP, posterior basomedial nucleus; Bnc, basal nucleus modified from de Olmos et al. (2004); Bnci, basal nucleus intermediate division; Bpc, BPc, basal nucleus parvicellular division; BSTAD, bed nucleus of the stria terminalis anterior dorsal area; BSTAV, bed nucleus of the stria terminalis anterior ventral area; BSTdm, bed nucleus of the stria terminalis anterior division dorsomedial nucleus; BSTif, bed nucleus of the stria terminalis interfascicular nucleus; BSTLD, bed nucleus of the stria terminalis lateral division dorsal part; BSTLP, bed nucleus of the stria terminalis lateral division posterior part; BSTLV, bed nucleus of the stria terminalis lateral division ventral part; BSTmc, bed nucleus of the stria terminalis magnocellular nucleus; BSTMPL, bed nucleus of the stria terminalis medial division posterolateral part; BSTpr, bed nucleus of the stria terminalis principal nucleus; BSTSl, supracapsular bed nucleus of the stria terminalis lateral part; BSTt, bed nucleus of the stria terminalis transverse nucleus; CE, Ce, CEA, central nucleus; CEl, central nucleus lateral division; CEAl, central nucleus, lateral part; CEAc, central nucleus, capsular part; CEAm, central nucleus, medial part; Cec, CEc, central nucleus capsular division; CeCD, central amygdaloid nucleus caudal division; CeI, central amygdaloid nucleus intermediate division; CeI, CEi, central nucleus intermediate division; CeL, central nucleus, lateral division; CeLC, central nucleus, lateral division, capsular; CeLCn, central nucleus, lateral division, central; CEm, central nucleus medial division; CeM, central nucleus, medial division; Co, cortical nucleus; COA, cortical nucleus of the amygala; COa, anterior cortical nucleus; COAa, layers 1-2 = Cortical nucleus, anterior part; COAp, posterior amygdaloid nucleus; COApl, cortical nucleus, posterior part, lateral zone; COApm, cortical nucleus, posterior part, medial zone; COApm, posterior amygdaloid nucleus medial part; COp, posterior cortical nucleus; EN, Eipert network; Fu, bed nucleus of the stria terminalis fusiform part; I, intercalated nuclei, Intercalated nuclei of the amygdala; IA, intercalated nuclei; IM, intercalated masses; i.n., included in; IPACL, interstitial nucleus of the posterior limb of the anterior commissure lateral part; IPACM, interstitial nucleus of the posterior limb of the anterior commissure medial part; I*, intercalated masses; L, La, lateral nucleus; LAM, lateral nucleus medial part; LaVL, lateral nucleus ventrolateral nucleus; LaVM, ventromedial lateral nucleus; LaVM, ventromedial part of the lateral nucleus; Ldl, LaDL, lateral nucleus dorsolateral division; Level, level of subdivision (first level is the supranuclear division of the amygdala); Lm, lateral nucleus medial division; LOT, nucleus of the lateral olfactory tract; LOTL2, nucleus of the lateral olfactory tract layer 2; LOTL3, nucleus of the lateral olfactory tract layer 3; Lvl, lateral nucleus ventrolateral division; M, medial nucleus; Mc, medial nucleus caudal division; Mcd, dorsal part of central division of the medial nucleus; Mcv, ventral part of central division of the medial nucleus; Me, MEA, medial nucleus; MEAad, medial nucleus, anterodorsal part; MEAav, medial nucleus, anteroventral part; MeAD, anterodorsal medial nucleus; MeAD, medial amygdaloid nucleus anterodorsal part; MEApd, medial nucleus, posterodorsal part; MEApv, medial nucleus, posteroventral part; MeAV, anteroventral medial nucleus; MeC, medial amygdaloid nucleus caudal part; MeCd, medial nucleus central division; MeCD, medial amygdaloid nucleus central dorsal part; MeCV, medial amygdaloid nucleus central ventral part; MePD, medial amygdaloid nucleus posterodorsal part; MePD, posterodorsal medial nucleus; MePV, medial amygdaloid nucleus posteroventral part; MePV, posteroventral medial nucleus; MeRo, medial amygdaloid nucleus rostral part; Mr, medial nucleus rostral division; NLOT, nucleus of the lateral olfactory tract; PA, posterior nucleus amygdala; PAC, periamygdaloid cortex; PACm, periamygdaloid complex medial division; PACm, periamygdaloid cortex, medial division; PACs, periamygdaloid complex sulcal division; PACs, periamygdaloid cortex, sulcal division; PLCo, posterolateral cortical nucleus; PMCo, posteromedial cortical nucleus; SLEAc, central division of sublenticular extended amygdala; SV, Subventricular nucleus; w.n.s., with no subdivisions.



Achard, S., and Bullmore, E. (2007). Efficiency and cost of economical brain functional networks. PLoS Comput. Biol. 3:e17. doi:10.1371/journal.pcbi.0030017

CrossRef Full Text

Adolphs, R., Tranel, D., Damasio, H., and Damasio, A. (1994). Impaired recognition of emotion in facial expressions following bilateral damage to the human amygdala. Nature 372, 669–672.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Aggleton, J. (2000). “The amygdala – what’s happened in the last decade?” in The Amygdala: A Functional Analysis, ed. J. P. Aggleton (New York: Oxford University Press), 2–29.

Arenas, A., Fernández, A., and Gómez, S. (2008). A complex network approach to the determination of functional groups in the neural system of C. elegans. Lect. Notes Comput. Sci. 5151, 9–18.

CrossRef Full Text

Baltz, A., and Kliemann, L. (2004). “Spectral analysis,” in Network Analysis. Lecture Notes in Computer Science, Vol. 3418, eds U. Brandes and T. Erlebach (Berlin: Springer), 373–416.

Bassett, D. S., and Bullmore, E. (2006). Small-world brain networks. Neuroscientist 12, 512–523.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bezgin, G., Reid, A. T., Schubert, D., and Kötter, R. (2009). Matching spatial with ontological brain regions using java tools for visualization, database access, and integrated data analysis. Neuroinformatics 7, 7–22.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Blinder, P., Baruchi, I., Volman, V., Levine, H., Baranes, D., and Jacob, E. (2005). Functional topology classification of biological computing networks. Nat. Comput. 4, 339–361.

CrossRef Full Text

Bohland, J., Bokil, H., Allen, C., and Mitra, P. (2009). The brain atlas concordance problem: quantitative comparison of anatomical parcellations. PLoS ONE 4:e7200. doi:10.1371/journal.pone.0007200

CrossRef Full Text

Bota, M., and Swanson, L. (2008). Bams neuroanatomical ontology: design and implementation. Front. Neuroinform. 2:2. doi:10.3389/neuro.11.002.2008

CrossRef Full Text

Brockhaus, H. (1938). Zur normalen und pathologischen Anatomie des Mandelkerngebietes. J. Psychol. Neurol. 49, 1–136.

Burdach, K. (1819–1826). Vom Baue und Leben des Gehirns. Leipzig: Dyk’sche Buchhandlung.

Burns, G. (1997). Neural Connectivity of the Rat: Theory, Methods and Applications. Ph.D. thesis, University of Oxford, Oxford.

Burns, G., and Cheng, W.-C. (2006). Tools for knowledge acquisition within the neuroscholar system and their application to anatomical tract-tracing data. J. Biomed. Discov. Collab. 1, 10–16.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Burns, G., Feng, D., and Hovy, E. (2008). “Intelligent approaches to mining the primary research literature: Techniques, systems, and examples,” in Computational Intelligence in Medical Informatics. Studies in Computational Intelligence, Vol. 85, eds A. Kelemen, Y. Liang, and A. V. Vasilakos (Berlin: Springer), 17–50.

Burns, G., and Young, M. (2000). Analysis of the connectional organization of neural systems associated with the hippocampus in rats. Philos. Trans. R. Soc. Lond. B Biol. Sci. 355, 55–70.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Castro, J., Varea, E., Márquez, C., Cordero, M., Poirier, G., and Sandi, C. (2010). Role of the amygdala in antidepressant effects on hippocampal cell proliferation and survival and on depression-like behavior in the rat. PLoS ONE 5:e8618. doi:10.1371/journal.pone.0008618

CrossRef Full Text

Chen, G.-D., Manohar, S., and Salvi, R. (2012). Amygdala hyperactivity and tonotopic shift after salicylate exposure. Brain Res. 1485, 63–76.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chen, L., Liu, J., Zhang, Q., Feng, J., Gui, Z., Ali, U., et al. (2011). Alterations of emotion, cognition and firing activity of the basolateral nucleus of the amygdala after partial bilateral lesions of the nigrostriatal pathway in rats. Brain Res. Bull. 85, 329–338.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cimino, J., and Zhu, X. (2006). The practical impact of ontologies on biomedical informatics. Yearb. Med. Inform. 2006, 200–211.

Ciocchi, S., Herry, C., Grenier, F., Wolff, S. B., Letzkus, J. J., Vlachos, I., et al. (2010). Encoding of conditioned fear in central amygdala inhibitory circuits. Nature 468, 277–282.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

da Costa, L. F., Rodrigues, F. A., Travieso, G., and Boas, P. V. (2007). Characterization of complex networks: a survey of measurements. Adv. Phys. 56, 167–242.

CrossRef Full Text

da Costa, L. F., and Sporns, O. (2006). Correlating thalamocortical connectivity and activity. Appl. Phys. Lett. 89, 1–3.

Davis, M. (1992). The role of the amygdala in fear and anxiety. Annu. Rev. Neurosci. 15, 353–375.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

de Olmos, J. S., Beltramino, C. A., and Alheid, G. (2004). “Amygdala and extended amygdala of the rat: a cytoarchitectonical fibroarchitectonical, and chemoarchitectonical survey,” in The Rat Nervous System, ed. G. Paxinos (San Diego: Elsevier Academic Press), 509–603.

Dong, J., and Horvath, S. (2007). Understanding network concepts in modules. BMC Syst. Biol. 1:24. doi:10.1186/1752-0509-1-24

CrossRef Full Text

Dwyer, D., and Killcross, S. (2006). Lesions of the basolateral amygdala disrupt conditioning based on the retrieved representations of motivationally significant events. J. Neurosci. 2, 8305–8309.

CrossRef Full Text

Echtermeyer, C., da, F., Costa, L., Rodrigues, F., and Kaiser, M. (2011). Automatic network fingerprinting through single-node motifs. PLoS ONE 6:e15765. doi:10.1371/journal.pone.0015765

CrossRef Full Text

Essen, D. V., Ugurbil, K., Auerbach, E., Barch, D., Behrens, T., Bucholz, R., et al. (2012). The human connectome project: a data acquisition perspective. Neuroimage 62, 2222–2231.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Estrada, E. (2010). Quantifying network heterogeneity. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 82, 066102.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Estrada, E., and Hatano, N. (2008). Communicability in complex networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 77, 036111.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fagiolo, G. (2007). Clustering in complex directed networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 76, 026107.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Felleman, D., and Essen, D. V. (1991). Distributed hierarchical processing in the primate cerebral cortex. Cereb. Cortex 1, 1–47.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gale, G., Anagnostaras, S., Godsil, B., Mitchell, S., Nozawa, T., Sage, J., et al. (2004). Role of the basolateral amygdala in the storage of fear memories across the adult lifetime of rats. J. Neurosci. 24, 3810–3815.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Goodhill, G., Simmen, M., and Willshaw, D. (1995). An evaluation of the use of multidimensional scaling for understanding brain connectivity. Philos. Trans. R. Soc. Lond. B Biol. Sci. 348, 265–280.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Goryczka, S., and Arodz, T. (2006). Complex-network-based methodology for analysis of biomedical data. BioAlg. MedSys. 3, 19–26.

Groenewegen, H., Wright, C., and Uylings, H. (1997). The anatomical relationships of the prefrontal cortex with limbic structures and the basal ganglia. J. Psychopharmacol. (Oxford) 11, 99–106.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hart, G., Harris, J. A., and Westbrook, R. F. (2010). Systemic or intra-amygdala infusion of the benzodiazepine, midazolam, impairs learning, but facilitates re-learning to inhibit fear responses in extinction. Learn. Mem. 17, 210–220.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Haubensak, W., Kunwar, P. S., Cai, H., Ciocchi, S., Wall, N. R., Ponnusamy, R., et al. (2010). Genetic dissection of an amygdala microcircuit that gates conditioned fear. Nature 468, 270–276.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hilgetag, C. C., O’Neill, M. A., and Young, M. P. (1996). Indeterminate organization of the visual system. Science 271, 776–777.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hilgetag, C. C., O’Neill, M. A., and Young, M. P. (2000). Hierarchical organization of macaque and cat cortical sensory systems explored with a novel network processor. Philos. Trans. R. Soc. Lond. B Biol. Sci. 355, 71–89.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Honey, C., Kötter, R., Breakspear, M., and Sporns, O. (2007). Network structure of cerebral cortex shapes functional connectivity on multiple time scales. Proc. Natl. Acad. Sci. U.S.A. 104, 10240–10245.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Honey, C. J., and Sporns, O. (2008). Dynamical consequences of lesions in cortical networks. Hum. Brain Mapp. 29, 802–809.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Humphries, M., and Gurney, K. (2008). Network “small-world-ness”: a quantitative method for determining canonical network equivalence. PLoS ONE 3:e0002051. doi:10.1371/journal.pone.0002051

CrossRef Full Text

Humphries, M., Gurney, K., and Prescott, T. (2006). The brainstem reticular formation is a small-world, not scale-free, network. Proc. Biol. Sci. 273, 503–511.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Johnston, J. (1923). Further contributions to the study of the evolution of the forebrain. J. Comp. Neurol. 35, 337–482.

CrossRef Full Text

Jones, B. F., Groenewegen, H. J., and Witter, M. P. (2005). Intrinsic connections of the cingulate cortex in the rat suggest the existence of multiple functionally segregated networks. Neuroscience 133, 193–207.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Joyce, K. E., Laurienti, P. J., Burdette, J. H., and Hayasaka, S. (2010). A new measure of centrality for brain networks. PLoS ONE 5:e12200. doi:10.1371/journal.pone.0012200

CrossRef Full Text

Kammer, F., and Täubig, H. (2004). “Connectivity,” in Network Analysis. Lecture Notes in Computer Science, Vol. 3418, eds U. Brandes and T. Erlebach (Berlin: Springer), 143–177.

Kim, J., and Jung, M. (2006). Neural circuits and mechanisms involved in pavlovian fear conditioning: a critical review. Neurosci. Biobehav. Rev. 30, 188–202.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kim, H. J., and Kim, J. M. (2005). Cyclic topology in complex networks. Phys. Rev. E. Stat. Nonlin. Soft Matter Phys. 72, 036109.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kishi, T., Tsumori, T., Yokota, S., and Yasui, Y. (2006). Topographical projection from the hippocampal formation to the amygdala: a combined anterograde and retrograde tracing study in the rat. J. Comp. Neurol. 496, 349–368.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kleinberg, J. M. (1999). Authoritative sources in a hyperlinked environment. J. ACM 46, 604–632.

CrossRef Full Text

Kötter, R. (2002). Neuroscience Databases – A Practical Guide. Norwell, MA: Kluwer Academic Publishers.

Kötter, R., Reid, A. T., Krumnack, A., Wanke, E., and Sporns, O. (2007). Shapley ratings in brain networks. Front. Neuroinform. 1:2. doi:10.3389/neuro.11.002.2007

CrossRef Full Text

Kötter, R., and Stephan, K. E. (2003). Network participation indices: characterizing component roles for information processing in neural networks. Neural Netw. 16, 1261–1275.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

MacDonald, N. (1983). Trees and Networks in Biological Models. Chichester: John Wiley & Sons Ltd.

Meunier, M., Bachevalier, J., Murray, E. A., Málková, L., and Mishkin, M. (1999). Effects of aspiration versus neurotoxic lesions of the amygdala on emotional responses in monkeys. Eur. J. Neurosci. 11, 4403–4418.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Meynert, T. (1867). Der Bau der Gross-Hirnrinde und seine örtlichen Verschiedenheiten, nebst einem pathologisch-anatomischen Corollarium. Vierteljahrsschrift für Psychiatrie 2, 77–93, 125–170, 198–217.

Modha, D., and Singh, R. (2010). Network architecture of the long-distance pathways in the macaque brain. Proc. Natl. Acad. Sci. U.S.A. 107, 13485–13490.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Musen, M., Noy, N., O’Connor, M., Redmond, T., Rubin, D., Tu, S., et al. (2009). Protégé. Available at:

Nagyessy, L., Nepusz, T., Kocsis, L., and Bazso, F. (2006). Prediction of the main cortical areas and connections involved in the tactile function of the visual cortex by network analysis. Eur. J. Neurosci. 23, 1919–1930.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Newman, M., and Girvan, M. (2004). Finding and evaluating community structure in networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 69, 1–15.

CrossRef Full Text

Newman, M. E. J. (2004). Fast algorithm for detecting community structure in networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 69, 066133.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ozik, J., Hunt, B., and Ott, E. (2004). Growing networks with geographical attachment preference: Emergence of small worlds. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 69, 026108.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Palombi, O., Shin, J.-W., Watson, C., and Paxinos, G. (2006). Neuroanatomical affiliation visualization-interface system. Neuroinformatics 4, 299–317.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Paxinos, G., and Watson, C. (1986). The Rat Brain in Stereotaxic Coordinates, 2nd Edn. Amsterdam: Academic Press.

Paxinos, G., and Watson, C. (2007). The Rat Brain in Stereotaxic Coordinates, 6th Edn. Amsterdam: Elsevier Academic Press.

Pitkänen, A. (2000). “Connectivity of the rat amygdaloid complex,” in The Amygdala: A Functional Analysis, ed. J. P. Aggleton (New York: Oxford University Press), 31–115.

Poliakov, A., Hertzenberg, X., Moore, E., Corina, D., Ojemann, G., and Brinkley, J. (2007). Unobtrusive integration of data management with fMRI analysis. Neuroinformatics 5, 3–10.

Pubmed Abstract | Pubmed Full Text

Price, J. (1981). “Connectivity of the rat amygdaloid complex,” in The Amygdaloid Complex, INSERM Symposium No. 20 (Amsterdam: Elsevier/North-Holland Biomedical Press), 13–18.

Raghuraj, R. K., and Lakshminarayanan, S. (2006). “Alternate complexity measures and stability analysis of process and biological networks,” in Proceedings of the 11th APCChE Congres 2006 27th–30th August 2006, Kuala Lumpur.

Rosen, J. B., and Davis, M. (1990). Enhancement of electrically elicited startle by amygdaloid stimulation. Physiol. Behav. 48, 343–349.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52, 1059–1069.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Russchen, F. (1986). “Cortical and subcortical afferents of the amygdaloid complex,” in Excitatory Amino Acids and Epilepsy, eds R. Schwarcz and Y. Ben-Ari (New York: Plenum Press), 35–52.

Sah, P., Faber, E., Lopez de Armentia, M., and Power, J. (2003). The amygdaloid complex: Anatomy and physiology. Physiol. Rev. 83, 803–834.

Pubmed Abstract | Pubmed Full Text

Scannell, J. W., Burns, G. A. P. C., Hilgetag, C. C., O’Neil, M. A., and Young, M. P. (1999). The connectional organization of the corticothalamic system of the cat. Cereb. Cortex 9, 277–299.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schmitt, O., and Eipert, P. (2012). neuroVIISAS: approaching multiscale simulation of the rat connectome. Neuroinformatics 10, 243–267.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schmitt, O., Usunoff, K., Lazarov, N., Itzev, D., Eipert, P., Rolfs, A., et al. (2012). Orexinergic innervation of the extended amygdala and basal ganglia in the rat. Brain Struct. Funct. 217, 233–256.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schreiber, T. (2000). Measuring information transfer. Phys. Rev. Lett. 85, 461–464.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sporns, O. (2002). “Graph theory methods for the analysis of neural connectivity patterns,” in Neuroscience Databases, ed. R. Kötter (Boston: Kluwer Academic), 171–185.

Sporns, O. (2011). Networks of the Brain, Cambridge: The MIT Press.

Sporns, O., Honey, C. J., and Kötter, R. (2007). Identification and classification of hubs in brain networks. PLoS ONE 2:e1049. doi:10.1371/journal.pone.0001049

CrossRef Full Text

Sporns, O., and Kötter, R. (2004). Motifs in brain networks. PLoS Biol. 2:e369. doi:10.1371/journal.pbio.0020369

CrossRef Full Text

Sporns, O., Tononi, G., and Edelman, G. M. (2000). Connectivity and complexity: the relationship between neuroanatomy and brain dynamics. Neural Netw. 13, 909–922.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sporns, O., Tononi, G., and Edelman, G. M. (2002). Theoretical neuroanatomy and the connectivity of the cerebral cortex. Behav. Brain Res. 135, 69–74.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sporns, O., and Zwi, J. (2004). The small world of the cerebral cortex. Neuroinformatics 2, 145–162.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Stephan, K. E., Hilgetag, C. C., Burns, G. A. P. C., O’Neill, M. A., Young, M. P., and Kötter, R. (2000). Computational analysis of functional connectivity between areas of primate cerebral cortex. Philos. Trans. R. Soc. Lond. B Biol. Sci. 355, 111–126.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sugar, J., Witter, M., van Strien, N., and Cappaert, N. (2011). The retrosplenial cortex: intrinsic connectivity and connections with the (para)hippocampal region in the rat. An interactive connectome. Front. Neuroinform. 5:7. doi:10.3389/fninf.2011.00007

CrossRef Full Text

Swanson, L. (1992). Brain Maps: Structure of the Rat Brain, 1st Edn. Amsterdam: Elsevier Science Ltd.

Swanson, L., and Petrovich, G. (1998). What is the amygdala? Trends Neurosci. 21, 323–331.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tononi, G., and Sporns, O. (2003). Measuring information integration. BMC Neurosci. 4:31. doi:10.1186/1471-2202-4-31

CrossRef Full Text

Tononi, G., Sporns, O., and Edelman, G. (1994). A measure for brain complexity: relating functional segregation and integration in the nervous system. Proc. Natl. Acad. Sci. U.S.A. 91, 5033–5037.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tye, K. M., and Deisseroth, K. (2012). Optogenetic investigation of neural circuits underlying brain disease in animal models. Nat. Rev. Neurosci. 13, 251–266.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yamamoto, T., and Ueji, K. (2011). Brain mechanisms of flavor learning. Front. Syst. Neurosci. 5:76. doi:10.3389/fnsys.2011.00076

CrossRef Full Text

Young, A. W., Aggleton, J. P., Hellawell, D. J., Johnson, M., Broks, P., and Hanley, J. R. (1995). Face processing impairments after amygdalotomy. Brain 118, 15–24.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Young, M. (1992). Objective analysis of the topological organization of the primate cortical visual system. Nature 358, 152–155.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Young, M. (1993). The organization of neural systems in the primate cerebral cortex. Proc. Biol. Sci. 252, 13–18.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Young, M., Scannell, J., Burns, G., and Blakemore, C. (1994). Analysis of connectivity: neural systems in the cerebral cortex. Rev. Neurosci. 5, 227–250.

Pubmed Abstract | Pubmed Full Text

Zhang, S., Bodenreider, O., and Golbreich, C. (2006). Experience in reasoning with the foundational model of anatomy in owl dl. Pac. Symp. Biocomput. 2006, 200–211.

CrossRef Full Text

A. Appendix

The definitions of expressions, parameters, matrices, and simulation models (random graph models) used in the article are summarized in the following.

A.1. Basic Definitions

A.1.1. Node, vertex

The smallest subunit of a network. With regard to connectomes a node is a circumscribed or disjunctive region that contains neuron perikarya (sources of physiological action potential) and/or axonal terminals (targets of physiological action potentials).

A.1.2. Set of indexed nodes

The set of all indices of nodes is


A.1.3. Number of nodes

The number of nodes (regions, vertices) is


A.1.4. Edge

A directed edge (i, j) ∈ N × N is the line that connects vertices i and j with source i and target j. The set of directed edges E is


A.1.5. Edges

The number of edges (connections, links) ϵ is


A.1.6. Set of edges

The set of all not self-referencing edges is


A.1.7. Graph


A.1.8. Adjacency matrix

The adjacency matrix (connectivity matrix) A is


A.1.9. Weighted matrix

The weighted matrix W is


whereby wij is the weight of the edge (i, j) that connects i and j. 0 ≤ wij ≤ 1.

A.1.10. Path

A sequence of vertices (v1, …, vk) is a path from (v1 to vk) if ∀i ∈ {1, …, k − 1}: (v1, vi+1) ∈ E. The length of a path v1, …, vk is k − 1.

A.1.11. Distance matrix

The distance matrix D is



dij=d(i,j)=length of the shortestpath from i to j,if such a path exists,else(A11)

A.1.12. Degree all

Self-references of nodes are not considered for all three degree measures. deg(i) = degall(i)


A.1.13. Degree out


A.1.14. Degree in


A.1.15. Neighborhoods

Out-neighbors of i:


In-neighbors of i:


All neighbors of i:


A.2. Network Parameters

A.2.1. Leverage Lev(i)

The leverage centrality was introduced by Joyce et al. (2010).


A.2.2. Directed leverage Lev(i)


The weighted directed leverage Lev(i)w is analogue to Lev(i) with weighted degrees.

A.2.3. Communicability matrix G


The communicability is the sum of paths from i to j each weighted with 1k! for a path with length k.

A.2.4. Modularity measure

Let M = {M1, …, Mm} be a partition of N. Mi is a group, module, or cluster of vertices. With


the fraction of edges that fall within group MiN and


the fraction of ends of edges that are attached to vertices in group Mi, the modularity


whereby ai2 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


The method of Newman and Girvan (2004) was used.

A.2.5. Global efficiency (GE)


GE and GEw analog with d(i, j) and dw(i,j)

A.2.6. Directed global efficiency


A.2.7. Harmonic mean (HM)


The directed and weighted versions use the directed and weighted global efficiencies.

A.2.8. 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 local efficiency of the network is the average local efficiency of all nodes.

A.2.9. Directed local efficiency


A.2.10. Weighted directed local efficiency


Whereby ni = |Ni| and djk(Ni), respectively, djkwNi is the length of the shortest path between j and k that contains only neighbors of i.

A.2.11. Directed assortativity coefficient r


A.2.12. Directed and weighted assortativity coefficient rw


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.

A.2.13. Average path length = characteristic path length d¯

With P = {(i, j) ∈ N × N | d(i, j) < ∞}, the set of paths.


In the weighted case the distances d(i,j) are replaced by the weighted distances.

A.2.14. Average directed degree (deg)


A.2.15. Heterogeneity VC

Coefficient of variation (VC) of the degreeall parameter.


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.

A.2.16. Line density (Ld)

The line density is the same like the connectedness of Raghuraj and Lakshminarayanan (2006).


Without self-referencing edges.

A.2.17. Diameter (Diam)


A.2.18. Katz index

The Katz index (Katz status index, Katz centrality) is a measure for the direct and indirect input of a node.


The attenuation factor α 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 neuroVIISAS the Katz centrality is multiplied by the mean of the quotient degin(i)/CKatz(i) of all nodes with CKatz(i) > 0. Hence, the values lie in the same range as the indegrees.

A.2.19. Number of triangles


The maximum number of possible triangles that can be derived from a complete reciprocal triangle is 8.

A.2.20. Weighted number of triangles


Instead of the sum of triangles [t(i)] the sum of geometric means of edge weights of each triangle is calculated. The following example provides (wij · wjk · wik)1/3 as the summand:


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


A.2.22. Directed and weighted transitivity


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:


For the directed and weighted case:


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


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


Couti refers to the out-neighbors of i.


Cin(i) refers to the in-neighbors of i.


In the weighted case the aij are replaced by the wij.

A.2.25. Average cluster coefficient




A.2.26. Small-worldness S


d¯rand and Crand are the parameters C and d¯ of the Erdös Renyi randomizations.

A.2.27. Centrality


This centrality (degree centrality) is defined for an undirected network based on undirected degrees. A directed or weighted version is not available yet. For the calculation the directed network is transferred to an undirected one.

A.2.28. Circle length (LC)


A.2.29. Eccentricity out

Eccentricity out, the output eccentricity of the vertex i is the maximum distance from i to any vertex.


A.2.30. Eccentricity in

Eccentricity in, the input eccentricity of the vertex i is the maximum distance from any vertex to i.


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


the set of second neighbors of node i is:

C2(i)=1|N2(i)|(|N2(i)|-1)j,kN2(i)jkajk, if |N2(i)|>10, otherwise(A59)

In the weighted case the aij are replaced by wij.

A.2.32. Average neighbor degree

The non-weighted average neighbor degree NB(i) of node i is


A.2.33. Weighted average neighbor degree

The weighted average neighbor degree NB(i) of node i is


A.2.34. Variation coefficient of neighbor degree


The weighted case is analog.

A.2.35. 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+.


The weighted case is analog. 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.

A.2.36. Closeness centrality out CCout(i)

The closeness centrality out with indices of nodes from which node i can be reached RNOUT(i) = {jN{i}|d(i, j) < ∞}


A.2.37. Closeness centrality in CCin(i)

The closeness centrality in with indices of nodes which can be reached from node i RNIN(i) = {jN\{i}|d(j, i) < ∞}


A.2.38. Betweenness centrality (BC)


Where ρj,k is the number of shortest paths from j to k and ρj,k(i) is the number of shortest paths from j to k that pass through i. The directed and weighted definitions are the same.

A.2.39. Stress (S)


The directed and weighted definitions are the same.

A.2.40. Central point distance (CPD)


Where BCmax=iNmaxBCi is the maximum betweenness centrality. The directed and weighted versions use the directed and weighted betweenness centralities.

A.2.41. Participation coefficient

The partition M = {M1, …, Mm} is generated as described in the definition of modularity.


with x ∈ {in, out, all} and


(Number of edges from vertices of Mj to i).


(Number of edges from i to vertices of Mj).


(Number of edges between i and vertices of Mj ).


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.

A.2.42. Z-score/within module degree

Let Mi be the module containing node i. degx(i, Mi) x ∈ {in, out, all} is explained in the definition of the participation coefficient.


is the mean and


the standard deviation of the within module Mi degree distribution. Then the z-score is defined as


and analog


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.

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

A.2.44. Shapley rating ϕ

The Shapley rating is a measure that provides information about the loss of connectivity following the removal of a node.


Where SCCN^ is the set of strongly connected components of N^. The smaller the value is, the more important is the 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.

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

A.2.46. Input radiality (Radin)

The input radiality of a node Radin is


In the weighted case the weighted distances are used.

A.2.47. Output radiality (Radout)

The output radiality of a node Radout is


In the weighted case the weighted distances are used.

A.2.48. Centroid value (Cen)

With gout(i, j) = |{kN|d(i, k) < d(j, k) < ∞}| and gin(i, j) = |{kN|d(k, i) < d(k, j) < ∞}| which are the number of nodes closer to node i than to node j with regard to In- and Out-distance, the centroid value is defined in the following.

A.2.49. Output centroid value (Cenout)


A.2.50. Input centroid value (Cenin)


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.

A.2.51. Page rank centrality (PRC)

PRC(i) = ri where r is the solution of the linear system


with the damping factor α = 0.85, the identity matrix I, and the diagonal matrix B, whereby


In the weighted case the weight matrix W is used instead of A and the weighted version degoutwi of the outdegree.

A.2.52. 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 such paths.


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.


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


A.2.54. Subgraph centrality (SC)


The subgraph centrality of the network is the average subgraph centrality of its nodes.


A.2.55. Undirected cyclic coefficient (CyclC)

The undirected cyclic coefficient as published by Kim and Kim (2005).



disti(j,k)=length of the shortest path from j to k that does not contains i,if such a path exists,otherwise(A94)

A.2.56. Directed cyclic coefficient (CyclC)

A publication about the directed cyclic coefficient is unknown. The directed cyclic coefficient is implemented here as follows:


A.2.57. Directed weighted cyclic coefficient (CyclCw)

A publication about the directed weighted cyclic coefficient is unknown. The directed weighted cyclic coefficient is implemented here as follows:


with distiw(j,k) is the weighted version of disti(j, k) with the weighted path length.

A.2.58. Cyclic network coefficient (CyclC)

The cyclic coefficient of the network is the average cyclic coefficient of its nodes:


A.2.59. Independent cycle count (CyC)

The number of independent cycles is:


The relative number of independent cycles is:


A.2.60. Vulnerability (V)

The vulnerability V is the maximum relative decrease of the global efficiency removing a single node.


Where GE(i) is the global efficiency of the graph (N\{i},{(j, k) ∈ E | jik}) that originates by removal of node i and all edges adjacent to i. The weighted version is analog using the weighted global efficiencies.

A.2.61. Hubness

The hubness is the degree to which a node has output to authorities (Kleinberg, 1999).

A.2.62. Authoritativeness

The authoritativeness is the degree to which a node gets input from hubs (Kleinberg, 1999).

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

A.3.1. Erdös Renyi graph


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 Erdös Renyi random graph is binominal in terms of


A.3.2. 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 <k> closest neighbors using are wiring 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:


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


The degree distribution of a Barabasi-Albert network is scale-free following the power law distribution of the form:


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

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

A.3.6. Rewiring graph

The rewiring-models connects each target of an edge of a network to another target node.

A.3.7. Power law


Δ is the deviation (error) of an empirical distribution of degrees from the power law function. A small Δ value means that the empirical distribution is similar with the power law function.

Keywords: amygdala, connectome, tract-tracing, network analysis, stereotaxic atlas, visualization, rat brain, simulation

Citation: Schmitt O, Eipert P, Philipp K, Kettlitz R, Fuellen G and Wree A (2012) The intrinsic connectome of the rat amygdala. Front. Neural Circuits 6:81. doi: 10.3389/fncir.2012.00081

Received: 07 June 2012; Accepted: 24 October 2012;
Published online: 11 December 2012.

Edited by:

Luke R. Johnson, Uniformed Services University of the Health Sciences, USA

Reviewed by:

Astrid A. Prinz, Emory University, USA
Donald Rainnie, Emory University, USA

Copyright: © 2012 Schmitt, Eipert, Philipp, Kettlitz, Fuellen and Wree. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.

*Correspondence: Oliver Schmitt, Department of Anatomy, Gertrudenstr. 9, Institute for Biostatistics and Informatics in Medicine and Ageing Research, University of Rostock, Ernst-Heydemann-Street 8, 18057 Rostock, Germany. e-mail: