Home | About | Journals | Submit | Contact Us | Français |

**|**Front Neuroinformatics**|**v.4; 2010**|**PMC2859882

Formats

Article sections

- Abstract
- Introduction
- Materials and Methods
- Results
- Summary and Discussion
- Conflict of Interest Statement
- References

Authors

Related links

Front Neuroinformatics. 2010; 4: 1.

PMCID: PMC2859882

Edited by: Claus C. Hilgetag, Jacobs University Bremen, Germany

Reviewed by: Steven Bressler, Florida Atlantic University, USA; David Meunier, University of Cambridge, UK

*Correspondence: Gorka Zamora-López, Interdisciplinary Center for Dynamics of Complex Systems, University of Potsdam, Komplex II – Golm (Haus 28) , Karl-Liebknecht-Str. 24, D-14476 Potsdam, Germany. e-mail: ed.madstop-inu.dlnga@akrog

Received 2009 April 6; Accepted 2010 February 2.

Copyright © 2010 Zamora-López, Zhou and Kurths.

This is an open-access article subject to an exclusive license agreement between the authors and the Frontiers Research Foundation, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are credited

This article has been cited by other articles in PMC.

Sensory stimuli entering the nervous system follow particular paths of processing, typically separated (segregated) from the paths of other modal information. However, sensory perception, awareness and cognition emerge from the combination of information (integration). The corticocortical networks of cats and macaque monkeys display three prominent characteristics: (i) modular organisation (facilitating the segregation), (ii) abundant alternative processing paths and (iii) the presence of highly connected hubs. Here, we study in detail the organisation and potential function of the cortical hubs by graph analysis and information theoretical methods. We find that the cortical hubs form a spatially delocalised, but topologically central module with the capacity to integrate multisensory information in a collaborative manner. With this, we resolve the underlying anatomical substrate that supports the simultaneous capacity of the cortex to segregate and to integrate multisensory information.

The mammalian nervous system is responsible for collecting and processing of information, and for providing adaptive responses which permit the organism to survive in a permanently changing environment. Sensory neurones encode environmental information into electrical signals which propagate in a “bottom-up” manner through different processing stages (Kandel et al., 2000; Bear et al., 2006). Each level provides responses of increasing complexity and at different time scales, e.g. reflex arcs, emotional responses and more elaborate cognitive responses. Information of the same modality (e.g. visual, auditory, somatosensory, etc.) traverses the body together, typically separated from the processing paths of other modalities. This permits that particular regions of the cortex specialise in detecting different features of the sensory stimuli, e.g. orientation, velocity and colour of the visual input; or frequency and pitch of the auditory stimuli. However, in order to generate a coherent perception of the reality, the brain needs to combine (integrate) this multisensory information at some place (Robertson, 2003) and during some time (Fahle, 1993; Singer and Gray, 1995; Engel and Singer, 2001). For that, the paths of information need to converge.

It has been argued that the functional capacity of the NS to balance between segregation (specialisation) and integration might be facilitated by its structural organisation (Sporns and Tononi, 2001). Analysis of the connectivity between regions of the cerebral cortex in macaque monkeys and cats has revealed the following characteristics: (i) clustered organisation of the cortical areas (Scannell and Young, 1993; Scannell et al., 1995; Hilgetag et al., 2000; Hilgetag and Kaiser, 2004) (see Figure Figure2),2), (ii) a large density of connections, and (iii) a broad degree distribution containing highly connected areas which are referred as *hubs* (Zemanová et al., 2006; Sporns et al., 2007; Zhou et al., 2007; Hagmann et al., 2008). While the organisation of cortical areas into clusters permits the segregated processing of information of different modality, the large number of connections involves that sensory information is highly accessible to all cortical areas, regardless of its modal origin. A detailed analysis of the corticocortical communication substrate has revealed the central role of the cortical hubs, by facilitating the communication between the different sensory modalities (Zamora-López et al., 2009).

Whether the cortical hubs act as passive transmitters of information, or they perform a more active function is a relevant question that we try to answer in the present paper. We start by summarising principles of complex network analysis and information theory in Section “Materials and Methods”. Section “Topological Capacity of Integration” contains a thorough application of graph theoretical measures which reveal that the cortical hubs form an additional module, expressed as a higher hierarchical level. In Section “Functional Capacity of Integration”, we challenge the intuitively assigned integrative properties of this central module by means of dynamical and information theoretical measures. On the one hand, we find that only simultaneous lesion of particular hubs leads to a dynamical segregation of the sensory modules (visual, auditory, somatosensory-motor and frontolimbic). On the other hand, the same hubs form a dynamical cluster after simultaneous excitation of primary sensory areas, a clear sign of their integrative capacities.

We first introduce basic concepts of graph theory. A network is an abstract manner to represent different aspects of a real system, providing it with a form (topology) which can be mathematically tractable. A network (*N*, *L*), composed of *N* nodes interconnected by *L* links, is described by an adjacency matrix **A** with entries *A _{ij}* =1 when there is a link pointing from node

The distance *d _{ij}* between two nodes

$${C}_{B}(i)={\displaystyle \sum _{s\ne i,t\ne i}^{N}\frac{{\text{\sigma}}_{st}(i)}{{\text{\sigma}}_{st}}}$$

(1)

where σ* _{st}*(

The topological similarity of two nodes can be characterised as the number of common neighbours they share. In the extreme case, two nodes are topologically identical if both have the same set of connections. The neighbourhood of node *i* is defined as the set of nodes it connects with, Г(*i*) ={*j* :* A _{ij}* =1}. In graphs without multiple links the size of the neighbourhood |Г(

$$MI(i,j)=\frac{|\text{\Gamma}(i)\cap \text{\Gamma}(j)|}{|\text{\Gamma}(i)\cup \text{\Gamma}(j)|}=\frac{{\displaystyle {\sum}_{n,m=1}^{N}{A}_{in}{A}_{jm}}}{k(i)+k(j)-{\displaystyle {\sum}_{n,m=1}^{N}{A}_{in}{A}_{jm}}}$$

(2)

Now, *MI*(*i*,*j*)=1 only if *i* and *j* are connected exactly to the same nodes, Г(*i*)= Г(*j*), and *MI*(*i*,*j*)=0 if they have no common neighbours.

Graph theoretical measures help understand the topological organisation of networks. Equally relevant is to uncover the features which are characteristic to the underlying system and the fundamental properties of its development. In this sense, the question is not whether a graph measure takes a specific numerical value, but whether this value distinguishes the empirical network _{emp} from others of similar characteristics. For that, the formulation of appropriate null-models is required. A typical such null case is to generate surrogate networks with the same size *N*, number of links *L* and degree distribution *p*(*k*) as in _{emp}. The *link switching method* (Katz and Powell, 1957; Holland and Leinhardt, 1977; Rao et al., 1996; Kannan et al., 1999; Roberts, 2000) consists of the following iterative process: starting from _{emp}, at each iteration two links are chosen at random (*i*_{1}→ *j*_{1}) and (*i*_{2}→ *j*_{2}). The links are rewired as (*i*_{1}→ *j*_{2}) and (*i*_{2}→ *j*_{1}) provided that the new links do not already exist and do not introduce self-loops, i.e. *i*→ *i*. Repeating the process sufficient times the resulting surrogate network conserves the initial degree distribution but any higher order structure is destroyed.

The classical textbook illustration of the cerebral cortex as a surface (grey matter) which can be subdivided into functional or cytoarchitectonic regions is only a limited picture. Additionally, long-range fibres link the cortical areas via the white matter forming a complex network which is neither regular nor completely random. This intricated structure enhances the richness and complexity of information processing capabilities of the cerebral cortex. In this paper we focus on the analysis of the cortical connectivity of the cat because it is, up to date, the most complete and reliable dataset of this kind.

The dataset of the corticocortical connections within the cortex of cats was created after an extensive collation of literature reporting anatomical tract-tracing experiments (Scannell and Young, 1993; Scannell et al., 1995). It consists of a parcellation into 53 cortical areas and 826 fibres of axons between them as summarised in Figure Figure2.2. The connections are weighted according to the axonal density of the projections. After application of data mining methods (Scannell and Young, 1993; Hilgetag and Kaiser, 2004), the network was found to be organised into four distinguishable clusters which closely follow functional subdivisions: visual (V), auditory (A), somatosensory-motor (SM) and frontolimbic (FL).

In order to perform significance tests of the graph measures, an ensemble of 1000 surrogate networks has been created following the *link switching method* (see Section “Graph Analysis”). All the resulting networks have the same size *N*=53, the same number of links *L*=826 and the same degree distribution as the corticocortical network of the cat. To assure that any further internal structure is destroyed, each surrogate network is the product of 10× *L* iterations. In the following, this set will be referred as the rewired ensemble {_{1n}}, and the original corticocortical network of the cat as _{cat}. The ensemble average of graph measures applied on the surrogate set {_{1n}} will be considered as the expected values.

Information theory has been very successful to describe transmission of information, encoding and channel capacity. At the root of this success lies the original idea of Shannon to apply concepts of statistical physics to represent the nature of communication. Consider a system *A* with *M* possible states. That is, a measurement made on *A* yields the values *a*_{1},*a*_{2},…,*a _{M}*, with a probability

$$H(A)=-{\displaystyle \sum _{i=1}^{M}p\left({a}_{i}\right)\text{log\hspace{0.05em}}p\left({a}_{i}\right).}$$

(3)

The entropy can be interpreted as the *amount of surprise one should feel upon reading the result of a measurement* (Faser and Swinney, 1986). It vanishes when the system has only one accessible state because the value *a* is always obtained, i.e. there is no surprise. *H*(*A*) is maximum when all the states are equally likely, i.e. there are no preferred states.

The *statistical dependence* between two systems *x*_{1} and *x*_{2} is quantified by their mutual information:

$$MI\left({x}_{1},{x}_{2}\right)=H\left({x}_{1}\right)+H\left({x}_{2}\right)-H\left({x}_{1},{x}_{2}\right).$$

(4)

By definition, the joint entropy is *H*(*x*_{1},*x*_{2})≤ *H*(*x*_{1})+ *H*(*x*_{2}). The equality is only fulfilled if *x*_{1} and *x*_{2} are statistically *independent*, hence *MI*(*x*_{1},*x*_{2})=0, and otherwise *MI*(*x*_{1},*x*_{2})>0.

In a series of papers Tononi and Sporns proposed a particular measure of integration (Tononi and Sporns, 1994; Tononi et al., 1996, 1998). Given a system *X* composed of *N* subsystems *x*_{1}, integration is defined as:

$$I(X)={\displaystyle \sum _{i=1}^{N}}H\left({x}_{i}\right)-H(X)$$

(5)

where *H*(*x*_{i}) is the entropy of one subsystem and *H*(*X*)= *H*(*x*_{1},*x*_{2},…,*x _{n}*) is the joint entropy of the system considered as a whole.

The *steady-state* of a linear system whose *N* subsystems **x** =(*x*_{1},*x*_{2},…,*x _{n}*) are driven by a Gaussian noise ξ =(ξ

$$\text{x}=g{\text{A}}^{\text{t}}\text{x}+\text{\xi}.$$

(6)

In practical terms the variable *x _{i}* might be interpreted as the

$$I(X)=\frac{1}{2}\text{log}\left[\frac{{\displaystyle {\prod}_{i=1}^{N}{\nu}_{i}}}{\text{|}COV(X)\text{|}}\right].$$

(7)

This expression shows that *I*(*X*) of the linear system is properly normalised and is independent of system size *N*. The covariance matrix can be analytically computed by solving the system such that $\text{x}={\scriptscriptstyle \frac{1}{1-g{A}^{t}}}\cdot \text{\xi}=Q\cdot \text{\xi}$, and averaging over the states produced by successive values of ξ one finds: *COV*(*X*)= x · x* ^{t}*=(

To compare *I*(*X*) of different systems, the matrix *A ^{t}* needs to be adequately normalised because application of the same coupling strength

Finally, a proper coupling strength *g* needs to be chosen. For that, we have estimated the covariance matrices of the cat cortical network under different coupling strengths (Figure (Figure4).4). They are similar to the correlation patterns arising from more complex models (Zemanová et al., 2006; Honey et al., 2007; Zhou et al., 2006, 2007). This similarity indicates the validity of the simple linear System (6) for the exploratory purposes here intended. All networks considered in Section “Functional Capacity of Integration” are normalised by their first pole and a coupling strength of *g*=0.5 is applied. Unless otherwise stated, the noise level is set to ξ* _{i}*=1.0.

The entropy of a subset of systems *S* *X* can be obtained by first computing *COV*(*X*) as indicated above, and then extracting the covariance submatrix *COV*(*S*) out of *COV*(*X*) by considering only the elements *x _{i}*

$$I(S)={\displaystyle \sum _{{x}_{j}\in S}^{}H({x}_{j})-H(S)=\frac{1}{2}\mathrm{log}\left[\frac{{\prod}_{j=1}^{{N}_{s}}{\nu}_{j}}{\left|COV\left(S\right)\right|}\right]}.$$

(8)

In order to characterise the connectional organisation of the nervous system and to understand its functional implications, the complex network approach has been applied in the recent years, particularly at the level of the cerebral cortex. This analysis has revealed several organisation properties, e.g. the clustering of cortical areas according to their sensory modality (visual, auditory, somatosensory-motor and frontolimbic). Recently, it has been reported that communication paths between cortical areas in different sensory modules are not random, but mediated by the hubs of the network (Zamora-López et al., 2009). In this section we present a more detailed graph analysis aiming to characterise the potential function of the cortical hubs.

The betweenness centrality *C _{B}*(ν) quantifies the relevance of a node

The centrality of a node usually correlates with its degree, hence, it is trivial to find out that precisely the hubs have larger centrality. Drawing any further conclusion requires performing a proper significance test. For comparison, the average *C _{B}*(ν) of the nodes in all the 1000 rewired networks of the surrogate ensemble {

This significance test permits us to uncover the most likely candidates to be a hub of the network, not only in terms of their number of links, but considering their contribution for the corticocortical communications. The hubs found here are potential candidates to perform high level integration because they have access to the information of different modalities. However, with the current results we can only affirm with certainty that the hubs are useful for the *transmission* of information from one modality to another. Concluding whether they perform any further function or not, it requires a more careful analysis.

A relevant question is now whether the cortical hubs are functionally independent of each other, i.e. each hub has a specialised function, or they perform some collaborative function. A graph measure to characterise the relation between the hubs of a network is the *rich-club phenomenon*. The k-density ϕ(*k*), is defined as the internal density of links between the nodes with degree larger than *k*′:

$$\text{\varphi}({k}^{\prime})=\frac{{L}_{{k}^{\prime}}}{{N}_{{k}^{\prime}}\left({N}_{{k}^{\prime}}-1\right)},$$

(9)

where ${N}_{{k}^{\prime}}$ is the number of nodes with *k*(*v*)≥ *k*^{′} and ${L}_{{k}^{\prime}}$ is the number of links between them. Notice that ϕ(*k*) is an increasing function of *k*. As $\text{\varphi}(0)={\scriptscriptstyle \frac{L}{N(N-1)}}$ is the density of the network, after the nodes with low degrees are removed the remaining reduced network contains more links per node. Thus, a plain measure of ϕ(*k*) is not very informative because hubs have a higher intrinsic chance of being connected to each other. Again, a conclusive interpretation requires the comparison to random networks with the same degree distribution. The question is then whether ϕ(*k*) of the real network grows faster or slower with *k* than the expected k-density ϕ_{1n}(*k*) out of the surrogate networks {_{1n}}. If ϕ(*k*) grows faster than ϕ_{1n}(*k*), it means that the hubs are more connected than expected and form a dense module (a *rich-club*). On the contrary, if ϕ(*k*) grows slower than ϕ_{1n}(*k*), the hubs are more independent of each other than expected.

In Figure Figure6A6A the k-density ϕ_{cat}(*k*) of _{cat} is presented together with the ensemble average ϕ_{1n}(*k*). For low degrees, ϕ_{cat}(*k*) follows very close the expectation, but for degrees *k*(*v*)>15, ϕ_{cat}(*k*) starts to grow faster showing that the hubs of the network form a *rich-club*. The largest difference occurs for *k*=23, comprising of 11 cortical hubs from all the four sensory systems (Figure (Figure6B).6B). Compared to the internal density of the four modules of the network, we find that the hubs form an even denser module (Table (Table11).

A central assumption in systems neuroscience is that the function of brain regions are specified by their afferents and efferents (Passingham et al., 2002). Under this assumption, it is to be expected that cortical areas of similar function, i.e. specialised in the processing of same modal information, should display a similar pattern of projections. In the case of the cortical hubs, it has been shown in the previous section that they form a tightly connected module. Whether this module could be regarded as a functional module, at least from a topological point of view, is the goal of the following analysis. The matching index *MI*(*v*,*v*′) is a graph measure to estimate the topological similarity of two nodes, by counting the number of common neighbours of *v* and *v*′ (Section “Graph Analysis”). In order to compare the values obtained for different pairs, the measure is normalised such that *MI*(*v*,*v*′)=1 only if all the neighbours of node *v* are also all the neighbours of *v*′. See the example in Figure Figure11.

We have computed the matching index for all pairs of cortical areas and the result is shown in matrix form, Figure Figure7A.7A. Visual inspection reveals the modular organisation of the network. This is reflected by the fact that *MI*(*v*,*v*′) is typically larger if both *v* and *v*′ belong to the same anatomical module, than if they belong to different modules. To highlight this difference, in Figure Figure7B7B the distribution of the matching values is shown: when the areas belong to the same module (internal matching), or to different modules (external matching). The external matching has a broad skewed distribution but peaking near *MI*=0.15. The internal matching displays a more constrained distribution with maximum at approximately *MI*=0.55. In Table Table11 the average matching of the network is compared to the average internal matching for each of the anatomical modules V, A, SM and FL. The internal averages are always larger than the global average despite the broad deviations, confirms the expected functional cohesiveness of the modules; not only in terms of their internal density of connections, but also in terms of their common connectivity.

As pointed out, the distribution of external matching is skewed and contains some larger values up to *MI*≈ 0.6. We find that most of these larger values are contributed precisely by the links between the cortical hubs which lie in different modules. We have recomputed the distribution of external matching, but ignoring the matching between the cortical hubs (solid line in Figure Figure7C).7C). The distribution decays now faster than in Figure Figure7B.7B. Finally, the distribution of the internal matching for the 11 hubs forming the rich-club is displayed (dotted line of Figure Figure7C).7C). It appears clearly separated from that of the distribution of external matching and peaking near *MI*=0.55. Its average is 0.52±0.10, comparable to, or larger than, the internal matching of the anatomical modules, Table Table1.1. These observations support the idea that the cortical hubs form a functional module on their own, as the anatomical modules do.

The two structural properties of the cortical hubs here presented, (*i*) hubs are densely connected with each other and (*ii*) they are functionally interrelated in terms of their inputs and outputs, extend the current understanding of cortical networks by uncovering that the multisensory hubs form yet another module which lies at a higher level in the hierarchical organisation. In the complex networks literature one finds two types of hierarchical topologies. The model after Arenas et al. (2006) considers hierarchies as the agglomeration of modules, say, small modules join to form larger modules, Figure Figure8A.8A. Another type of hierarchy after Ravasz and Barabási (2003) can be regarded as a tree-like fractal structure which produces modular networks with scale-free degree distribution. At each level, there is a central community connecting to all the modules at the same level, and to all modules in the hierarchies below. Such centralised patterns are repeated through different scales, Figure Figure88B.

The organisation that we uncover here is none of these two, but it might be regarded as a combination of them. Notice that in the model by Arenas et al. (2006), the small communities are randomly linked to each other such that their union forms a larger community. In the present case, the inter-community links are not random, but centralised. Therefore, the highest hierarchical level is formed by a *partial overlap* of the underlying modules. See Figure Figure8C8C for a schematic representation.

The functional implications of the topological findings described in this section, necessarily arise from intuitive interpretation of the intrinsic relationship between structure and function in neural systems. To provide a more solid ground to these intuitive interpretations, in the following section we challenge them by means of dynamical and information theoretical measures. We focus in a very simple dynamical model which has the benefit of being analytically solvable, although its validity for our purposes is confirmed by comparison to the dynamical output of more complex models, see Section “Information Theory and Integration”.

The structural organisation described in the previous section supports the idea that the cortical hubs might be responsible for combining the multisensory information hence facilitate the emergence of a global (integrated) perception. In this sense, we aim for a definition of integration which characterises *the capacity of one or more nodes to receive information of different character and combine it to produce new useful information*. Certainly, this definition involves crucial theoretical problems, e.g. what the *character* of information is, or what are the rules under which information is *combined*. Nevertheless, within a networked system, the nodes with a capacity to integrate information should obey certain measurable conditions. We propose the following:

1) Accessibility to information: A node can perform an integrative function only if it has general access to the information contained within the system.

2) Sharing of information: Two or more nodes can perform integrative function in a collaborative manner only if they are sufficiently connected with each other.

3) Segregation after selective damage: If a node has an integrative function, its removal should lead to a decrease of the integrative capabilities of the whole system.

From the structural point of view, the hubs listed in Figure Figure6B6B obey these three conditions. They are the most central areas and they are densely connected to each other. Besides, robustness studies (Kaiser et al., 2007) have shown that intentional lesion of the highly connected cortical areas largely affect the communication within the network. In the following, we introduce a framework to characterise the integrative function of the hubs by means of dynamical systems and information theory. Additionally, we perform a probabilistic analysis of the composition of the dynamical core, rather than a deterministic one. The reason is that even if the corticocortical networks of the cat is the most complete and reliable dataset of its kind up to date, it is not free of experimental errors. For example, some of the real connections might still be absent in the data. We aim to discriminate those hubs which, grouped together, possess a larger potential to integrate multisensory information from those groups which might have lesser capacities. For that, we arbitrarily choose all the areas with output degree *k*_{o}(ν)≥ 20 as potential members of the integrator module giving rise to a set of *N _{S}*=19 areas:

*S*_{hubs}= {20a, 7, AES, EPp, 6l, 6m, 5Am, 5Al, 5Bm, 5Bl, SSSAi, SSAo, PFCL, Ia, Ig, CGa, CGp, 35, 36}.

The statistical analysis consists in measuring the integrative capacities of all the 524,097 combinations of sizes *N _{S}*=1 to

Consider the linear System (6) with **Â** being the transposed and normalised adjacency matrix of the cat _{eat}. All areas are driven by a small Gaussian noise level ξ* _{i}*=1.0 and coupled by

The results depicted in Figure Figure9A9A show that *I ^{e}*(

Within a networked system the removal of critical nodes should lead to a decrease in its integrative capacities. In the following, we study the impact of targeted lesions of the corticocortical network of the cat, _{cat}. For all the possible subsets *S* composed of hubs in *S*_{hubs}, we perform a lesion to the network by simultaneously removing the nodes *x _{i}*

Recall that integration *I*(*X*) as defined in Eq. 7 is an extension of the mutual information for more than two systems. It represents the limit case in which the statistical dependence among all the elements *x _{i}* in the system

$${I}_{\text{\mathcal{P}}}\left(X\right)={\displaystyle \sum _{j=1}^{n}H\left({S}_{j}\right)}-H\left(X\right).$$

(10)

Note that when *n*= *N*, then *I*_{}(*X*)= *I*(*X*).

Considering the partition _{4}= {V, A, SM, FL} and the corticocortical network of the cat, then ${I}_{{\mathcal{P}}_{4}}\left({\mathcal{G}}_{cat}\right)=0.292$. The modular integration of each lesioned network _{S} is computed for the partition _{4}. Notice that (*a*) the nodes are also removed from the partition and (*b*) every _{S} is adequately normalised by its largest eigenvalue such that the measured observables are comparable across realisations (see Section “Information Theory and Integration”). The results in Figure Figure9C9C permit us again to discriminate between subsets of hubs whose simultaneous removal lead to a large segregation of the network, while removal of other subsets has barely no effect. For example, among all the possible lesions of size *N _{S}*=10, some trigger a large segregation of the modules ${I}_{{\mathcal{P}}_{4}}\left({\mathcal{G}}_{S}\right)\sim 0.05$ while other lesions do even increase their dynamical dependence: ${I}_{{\mathcal{P}}_{4}}\left({\mathcal{G}}_{S}\right)\sim 0.35>{I}_{{\mathcal{P}}_{4}}\left({\mathcal{G}}_{\text{cat}}\right)$.

Selecting only those subsets whose lesion leads to a larger segregation of the modules, i.e. ${I}_{{\mathcal{P}}_{4}}\left({\mathcal{G}}_{S}\right)$ lies among 10% of the minimal modular integration for each size *N _{S}* (red dots in Figure Figure9A),9A), a co-participation matrix

In summary, both the multiple lesion and the multisensory excitation analysis performed in this section lead to the identification of the same cortical hubs as responsible for the integration of multisensory information in the corticocortical network of the cat. Moreover, this set largely coincides with the top hierarchical level found by the graph analysis in Section “Topological Capacity of Integration”, corroborating the integrative function assigned to the hubs by intuitive interpretation of their topological characteristics.

In this paper we have analysed the modular and hierarchical organisation of the corticocortical network of the cat and its relationship to the intrinsic necessities of the brain to simultaneously segregate and integrate multisensory information. From the topological point of view, we have extended the current understanding of cortical organisation with the finding that the cortical hubs form a central module on top of the cortical hierarchy; which is expressed as the partial overlap of the four anatomical modules (visual, auditory, somatosensory-motor and frontolimbic). By means of dynamical and information theoretical measures, we have corroborated its capacity to integrate multisensory information, i.e. after simultaneous excitation of visual, auditory and somatosensory primary areas, a particular set of hubs becomes statistically dependent forming a dynamical cluster. Additionally, the simultaneous lesion of these hubs leads to a largest decrease in the integrative capacities of the network. Both structural and functional results indicate that visual areas 7 and AES, auditory area EPp and frontolimbic areas Ia, Ig, CGp, 35 and 36 are the most likely candidates to form the top hierarchical module. The participation of somatosensory-motor areas is less clear, although area 6m is the strongest candidate of them. Visual area 20a and somatosensory-motor areas 5Al and 5Bl are also potential candidates.

The modular and hierarchical organisation here detected agrees with the behaviour observed in dynamical simulations of cortical networks. The resting state dynamics are typically governed by the formation of dynamical clusters which closely relate to the anatomical modules, but the influence of the hierarchical organisation is also expressed. In Zemanová et al. (2006) and Zhou et al. (2006, 2007) it was shown that the correlation between the dynamical clusters is mediated by the cortical hubs. In Honey et al. (2007) the centrality of the hubs was found to oscillate in time. Simulation of excitable dynamics on hierarchical networks (Müller-Linow et al., 2008) has shown that the dynamical behaviour of the cortical network of the cat may be dominated either by the modular structure or by the hubs, depending on the time scales.

The separation of modal information paths is a relevant characteristic of organisation in the nervous system that permits simultaneous (parallel) processing of sensory input and detection of its features. Cortical regions containing neurones specialised in similar function, e.g. in processing information of the same sensory modality, lie geographically close to each other (Figure (Figure10A).10A). However, a coherent perception and the emergence of mental states such as awareness and consciousness require that information is integrated at different levels: the binding of sensory features into entities, the combination of entities with memories (personal experiences) into events, etc. While experimental techniques have led to a deep understanding about the basis of sensory perception, the nature of integration and the localisation of brain regions involved in it, is still under the subject of debate. As stated by Fuster (2003), simple extrapolation of the principles of sensory organisation do not lead to the identification of the substrate for cognition.

Several models have proposed that high-level functions are represented by distributed, interactive and overlapping networks of neurones, which transcend any of the traditional subdivisions of the cortex by structural (cytoarchitecture) or functional criteria (Damasio, 1989; Fuster, 2003, 2006; Tononi, 2004). During the recent years increasing experimental evidence has confirmed this hypothesis and the networked perspective has gained the favour against the assumption of a single brain region fully responsible for integration (Stam and Reijneveld, 2007; Bullmore and Sporns, 2009; Knight, 2009). The anatomical networked connectivity may serve as the basis in which localised and distributed functional networks rapidly emerge and dissolve governed by coordination dynamics according to the sensory stimulation and the ongoing activity (Bressler and Kelso, 2001).

As a further evidence, our results resolve the anatomical organisation substrate that supports the capacity of the cerebral cortex to simultaneously segregate and integrate information. In the light of this organisation, it could be envisioned that multisensory integration emerges from the collaborative function of the cortical hubs. While *early* sensory cortical regions perform specialised processing of the sensory input, the hubs of the network may work together to combine the multisensory information. A relevant organisation difference is that the cortical hubs form a module which is densely connected by axonal paths through the white matter, but is *geographically delocalised* (Figure (Figure1010B).

The current paper focuses in the corticocortical connectivity of cats because it is, up to date, the most complete and reliable dataset of its kind. Hence, it is the most suitable for a detailed and statistically consistent analysis. The main limitation is that it comprises of interconnection between cortical areas in only one cerebral hemisphere. Because of the known inter-hemisfere differences in many mammals, particularly in humans, it will be very valuable in the future to acquire the connectivity within and between both hemispheres in animal and human models. Based on current literature in which the cortical networks of the macaque and cat models display similar features, we expect that the general organisation principles here exposed to be valid in a wide range of mammals.

An interesting challenge is now to explain the emergence of this modular and hierarchical organisation in terms of evolution and development, in particular how the delocalised cluster of hubs could have evolved if, apparently, areas of similar function tend to be grouped close to each other. Very likely, the balancing between short wiring requirements (leading to minimisation of energy costs) and short processing paths allowing for robustness and fast responses (Kaiser and Hilgetag, 2006) plays a major role. It would also be of relevance to find out whether similar hierarchical patterns are repeated across smaller scales within the cortex, i.e. the interconnections between cortical columns and micro-columns. This would imply an underlying fractal-like complex architecture which can emerge from simple rules of assembly during development.

Finally, we should remind that current non-invasive techniques such as EEG, MEG and fMRI reveal only the presence of brain activity. They permit to identify which brain regions are *associated* with certain experimental condition. However, at the current stage it is very difficult, if not impossible, to understand *what* is exactly an activated region doing. Is it filtering a signal? Is it integrating information? Is an activation detected only because that particular region contains memories which are being retrieved and passed to other regions for processing? In our opinion, it would be highly interesting to further develop concepts of information theory as the modular and local capacity of integration here presented which applied to the time series of regional activity might help understand the particular function of individual brain regions within a given experimental task.

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.

We thank Lucia Zemanová, Claus-C. Hilgetag and Werner Sommer for valuable discussions. Gorka Zamora-López and Jürgen Kurths are supported by the Deutsche Forschungsgemeinschaft, research group FOR 868 (contract No. KU 837/23-1) and by the BioSim network of excellence (contract No. LSHB-CT-2004-005137 and No. 65533).

- Anthonisse J. M. (1971). The rush in a directed graph. Technical Report BN9/71. Amsterdam, Stichting Mathematisch Centrum
- Arenas A., Díaz-Guilera A., Pérez-Vicente C. (2006). Synchronization reveals topological scales in complex networks. Phys. Rev. Lett. 96, 114102.10.1103/PhysRevLett.96.114102 [PubMed] [Cross Ref]
- Bear M. F., Connors B. W., Paradiso M. A. (2006). Neuroscience: Exploring the Brain. Baltimore, MD, Lippincott Williams and Wilkins
- Bressler S. L., Kelso J. A. S. (2001). Cortical coordination dynamics and cognition. Trends Cogn. Sci. 5, 26–3610.1016/S1364-6613(00)01564-3 [PubMed] [Cross Ref]
- Bullmore E., Sporns O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 10, 1 [PubMed]
- Cover T. M., Thomas J. A. (1991). Elements of Information Theory. New York, Wiley; 10.1002/0471200611 [Cross Ref]
- Damasio A. R. (1989). Time-locked multiregional retroactivation: a systems-level proposal for the neuronal substrates of recall and recognition. Cognition 33, 25–6210.1016/0010-0277(89)90005-X [PubMed] [Cross Ref]
- Engel A. K., Singer W. (2001). Temporal binding and neural correlates of sensory awareness. Trends Cogn. Sci. 5, 16–2510.1016/S1364-6613(00)01568-0 [PubMed] [Cross Ref]
- Fahle M. (1993). Figure-ground discrimination from temporal information. Proc. R. Soc. Lond., B, Biol. Sci. 254, 199–20310.1098/rspb.1993.0146 [PubMed] [Cross Ref]
- Faser A. M., Swinney H. L. (1986). Independent coordinates for strange attractors from mutual information. Phys. Rev. A 33, 2318–2321 [PubMed]
- Freeman L. C. (1977). A set of measures of centrality based upon betweenness. Sociometry 40, 35–4110.2307/3033543 [Cross Ref]
- Fuster J. M. (2003). Cortex and Mind: unifying cognition. New York, Oxford University Press
- Fuster J. M. (2006). The cognit: a network model of cortical representation. Int. J. Psychophysiol. 60, 125–13210.1016/j.ijpsycho.2005.12.015 [PubMed] [Cross Ref]
- Graben P., Zhou C. S., Thiel M., Kurths J., editors. (eds). (2007). Lectures in Supercomputational Neuroscience. Berlin, Springer-Verlag
- Hagmann P., Cammoun L., Gigandet X., Meuli R., Honey C. J., Wedeen V. J., Sporns O. (2008). Mapping the structural core of human cerebral cortex. PLoS Biol. 6, e159.10.1371/journal.pbio.0060159 [PMC free article] [PubMed] [Cross Ref]
- Hilgetag C. C., Burns G. A., O'neill M. A., Scannell J. W., Young M. P. (2000). Anatomical connectivity defines the organization of clusters of cortical areas in the macaque monkey and the cat. Philos. Trans. R. Soc. Lond., B, Biol. Sci. 355, 91–11010.1098/rstb.2000.0551 [PMC free article] [PubMed] [Cross Ref]
- Hilgetag C. C., Kaiser M. (2004). Clustered organisation of cortical connectivity. Neuroinformatics 2, 353–36010.1385/NI:2:3:353 [PubMed] [Cross Ref]
- Holland P. W., Leinhardt S. (1977). The statistical analysis of local structure in social networks. In Sociological Methodology, Heise D. R., editor. , ed. (San Francisco, Jossey-Bass; ), pp. 1–45
- Honey C. J., Kötter R., Breakspear M., Sporns O. (2007). Network structure of cerebral cortex shapes functional connectivity on multiple time scales. Proc. Natl. Acad. Sci. U.S.A. 104, 10240–1024510.1073/pnas.0701519104 [PubMed] [Cross Ref]
- Kaiser M., Hilgetag C. C. (2006). Nonoptimal component placement, but short processing paths, due to long-distance projections in neural systems. PLoS Comput. Biol. 2, e95.10.1371/journal.pcbi.0020095 [PubMed] [Cross Ref]
- Kaiser M., Martin R., Andras P., Young M. P. (2007). Simulation of robustness against lesions of cortical networks. Eur. J. Neurosci. 25, 3185–319210.1111/j.1460-9568.2007.05574.x [PubMed] [Cross Ref]
- Kandel E. R., Schwartz J. H., Jessell T. M. (2000). Principles of Neural Science. New York, McGraw-Hill
- Kannan R., Tetali P., Vempala S. (1999). Simple Markov-chain algorithms for generating bipartite graphs and tournaments. Random Struct. Algorithms 14, 293–30810.1002/(SICI)1098-2418(199907)14:4<293::AID-RSA1>3.0.CO;2-G [Cross Ref]
- Katz L., Powell J. H. (1957). Probability distributions of random variables associated with a structure of the sample space of sociometric investigations. Ann. Math. Stat. 28, 442–44810.1214/aoms/1177706972 [Cross Ref]
- Knight R. T. (2009). Neural networks debunk phrenology. Science 316, 1578.10.1126/science.1144677 [PubMed] [Cross Ref]
- Kötter R., Sommer F. T. (2000). Global relationship between anatomical connectivity and activity propagation in the cerebral cortex. Philos. Trans. R. Soc. Lond., B, Biol. Sci. 355, 127–13410.1098/rstb.2000.0553 [PMC free article] [PubMed] [Cross Ref]
- Müller-Linow M., Hilgetag C. C., Hütt M. T. (2008). Organization of excitable dynamics in hierarchical biological networks. PLoS Comput. Biol. 4, e1000190.10.1371/journal.pcbi.1000190 [PMC free article] [PubMed] [Cross Ref]
- Newman M. E. J. (2003). The structure and function of complex networks. SIAM Rev. 45, 167–25610.1137/S003614450342480 [Cross Ref]
- Papoulis A. (1991). Probability, Random Variables and Stochastic Processes. New York, McGraw-Hill
- Passingham R. E., Stephan K. E., Kötter R. (2002). The anatomical basis of functional localization in the cortex. Nat. Rev. 3, 606. [PubMed]
- Rao A. R., Jana R., Bandyopadhyay S. (1996). A Markov chain Monte Carlo method for generating random (0, 1)-matrices with given marginals. Sankhya A 58, 225–242
- Ravasz E., Barabási A.-L. (2003). Hierarchical organization in complex networks. Phys. Rev. E 67, 026112.10.1103/PhysRevE.67.026112 [PubMed] [Cross Ref]
- Roberts J. M. (2000). Simple methods for simulating sociomatrices with given marginal totals. Soc. Netw. 22, 273–28310.1016/S0378-8733(00)00026-5 [Cross Ref]
- Robertson L. C. (2003). Binding, spatial attention and perceptual awareness. Nat. Rev. 4, 93–10210.1038/nrn1030 [PMC free article] [PubMed] [Cross Ref]
- Scannell J. W., Blakemore C., Young M. P. (1995). Analysis of connectivity in the cat cerebral cortex. J. Neurosci. 15, 1463–1483 [PubMed]
- Scannell J. W., Young M. P. (1993). The connectional organization of neural systems in the cat cerebral cortex. Curr. Biol. 3, 191–20010.1016/0960-9822(93)90331-H [PubMed] [Cross Ref]
- Shannon C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423
- Singer W., Gray C. M. (1995). Visual feature integration and the temporal correlation hypothesis. Annu. Rev. Neurosci. 18, 555–58610.1146/annurev.ne.18.030195.003011 [PubMed] [Cross Ref]
- Sporns O., Honey C. J., Kötter R. (2007). Identification and classification of hubs in brain networks. PLoS One 10, e1049.10.1371/journal.pone.0001049 [PMC free article] [PubMed] [Cross Ref]
- Sporns O., Tononi G. M. (2001). Classes of network connectivity and dynamics. Complexity 7, 28–3810.1002/cplx.10015 [Cross Ref]
- Stam C. J., Reijneveld J. C. (2007). Graph theoretical analysis of complex networks in the brain. Nonlinear Biomed. Phys. 1, 3.10.1186/1753-4631-1-3 [PMC free article] [PubMed] [Cross Ref]
- Tononi G. (2004). An information integration theory of consciousness. BMC Neurosci. 5, 42.10.1186/1471-2202-5-42 [PMC free article] [PubMed] [Cross Ref]
- Tononi G., McIntosh A. R., Russell D. P., Edelman G. M. (1998). Functional clustering: identifying strongly interactive brain regions in neuroimaging data. Neuroimage 7, 133.10.1006/nimg.1997.0313 [PubMed] [Cross Ref]
- Tononi G., Sporns O. (1994). A measure for brain complexity: relating functional segregation and integration in the nervous system. Proc. Natl. Acad. Sci. U.S.A. 91, 5033.10.1073/pnas.91.11.5033 [PubMed] [Cross Ref]
- Tononi G., Sporns O., Edelman G. M. (1996). A complexity measure for selective matching of signals by the brain. Proc. Natl. Acad. Sci. U.S.A. 93, 3422.10.1073/pnas.93.8.3422 [PubMed] [Cross Ref]
- Young M. P., Hilgetag C.-C., Scannell J. W. (2000). On imputing function to structure from the behavioural effects of brain lesions. Philos. Trans. R. Soc. Lond., B, Biol. Sci. 355, 147–16110.1098/rstb.2000.0555 [PMC free article] [PubMed] [Cross Ref]
- Zamora-López G., Zhou C. S., Kurths J. (2009). Graph analysis of cortical networks reveals complex anatomical communication substrate. Chaos 19, 015117.10.1063/1.3089559 [PubMed] [Cross Ref]
- Zemanová L., Zhou C. S., Kurths J. (2006). Structural and functional clusters of complex brain networks. Phy. D 224, 202–21210.1016/j.physd.2006.09.008 [Cross Ref]
- Zhou C. S., Zemanová L., Zamora-López G., Hilgetag C. C., Kurths J. (2006). Hierarchical organization unveiled by functional connectivity in complex brain networks. Phys. Rev. Lett. 97, 238103.10.1103/PhysRevLett.97.238103 [PubMed] [Cross Ref]
- Zhou C. S., Zemanová L., Zamora-López G., Hilgetag C. C., Kurths J. (2007). Structure-function relationship in complex brain networks expressed by hierarchical synchronization. New J. Phys. 9, 178.10.1088/1367-2630/9/6/178 [Cross Ref]

Articles from Frontiers in Neuroinformatics are provided here courtesy of **Frontiers Media SA**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |