Complex Brain Networks Graph Theoretical Analysis of Structural and Functional Systems

REVIEWS
Complex brain networks: graph

theoretical analysis of structural and
functional systems
Ed Bullmore* and Olaf Sporns
Abstract | Recent developments in the quantitative analysis of complex networks, based

largely on graph theory, have been rapidly translated to studies of brain network organization.
The brains structural and functional systems have features of complex networks such as
small-world topology, highly connected hubs and modularity both at the whole-brain
scale of human neuroimaging and at a cellular scale in non-human animals. In this article, we
review studies investigating complex brain networks in diverse experimental modalities
(including structural and functional MRI, diffusion tensor imaging, magnetoencephalography and electroencephalography in humans) and provide an accessible introduction to the
basic principles of graph theory. We also highlight some of the technical challenges and key
questions to be addressed by future developments in this rapidly moving field.
Graph theory
A branch of mathematics that
deals with the formal
description and analysis of
graphs. A graph is defined
simply as a set of nodes
(vertices) linked by connections
(edges), and may be directed
or undirected. When describing
a real-world system, a graph
provides an abstract
representation of the systems
elements and their interactions.
*University of Cambridge,
Behavioural & Clinical
Neurosciences Institute,
Department of Psychiatry,
Addenbrookes Hospital,
Cambridge, CB2 2QQ, UK.
GlaxoSmithKline,
Clinical Unit Cambridge,
Addenbrookes Hospital,
Cambridge, CB2 2QQ, UK.
Department of Psychological
and Brain Sciences, Indiana
University, Bloomington,
Indiana 47405, USA.
Correspondence to E.B.
email: [email protected]
doi:10.1038/nrn2575
Published online
4 February 2009
We have known since the nineteenth century that the

neuronal elements of the brain constitute a formidably
complicated structural network1,2. Since the twentieth
century it has also been widely appreciated that this
anatomical substrate supports the dynamic emergence
of coherent physiological activity, such as phase-locked
high-frequency electromagnetic oscillations, that can
span the multiple spatially distinct brain regions that
make up a functional network3,4. Such networks are
thought to provide the physiological basis for information processing and mental representations59. In this
article, we focus on graph theoretical approaches to the
analysis of complex networks that could provide a powerful new way of quantifying the brains structural and
functional systems (BOX 1).
Since the mid 1990s, developments in our understanding of the physics of complex systems1012 have led
to the rise of network science13 as a transdisciplinary
effort to characterize network structure and function.
In this body of literature, complexity arises in the macroscopic behaviour of a system of interacting elements that
combines statistical randomness with regularity 14. The
ascendancy of network science has been driven by the
growing realization that the behaviour of complex systems be they societies, cells or brains is shaped
by interactions among their constituent elements. The
increasing availability and tractability of large, highquality data sets on a wide range of complex systems1517
has led to a fundamental insight: substantively different

complex systems often share certain key organizational
principles, and these can be quantitatively characterized
by the same parameters (BOX 2). In other words, many
complex systems show remarkably similar macroscopic
behaviour despite profound differences in the microscopic details of the elements of each system or their
mechanisms of interaction.
One example of an apparently ubiquitous macroscopic behaviour in complex systems is the small-world
phenomenon18 (BOX 3). Recently, small-world architectures have been found in several empirical studies of
structural and functional brain networks1922 in humans
and other animals, and over a wide range of scales in
space and time; small-worldness is naturally therefore
a key topic for this Review. However, discovering that
brain networks are small-world networks is only the
first step towards a comprehensive understanding of
how these networks are structurally organized and how
they generate complex dynamics. In network science,
methodological advances allow us to quantify other
topological properties of complex systems such as
modularity 23, hierarchy 24, centrality 25 and the distribution of network hubs26,27 many of which have already
been measured in brain networks. There have also been
significant efforts to model the development or evolution of complex networks28, to link network topology to
network dynamics, and to explore network robustness
186 | maRcH 2009 | VOlume 10
www.nature.com/reviews/neuro
2009 Macmillan Publishers Limited. All rights reserved
REVIEWS
Box 1 | Structural and functional brain networks
Anatomical parcellation
Histological or
imaging data
Recording sites
Time series data
Structural brain network
Functional brain network

Sensorimotor
Premotor
Parietal
Prefrontal
Occipital
Inferior temporal
Orbitofrontal
Temporal pole
Graph theoretical analysis
Structural and functional brain networks can be explored using graph theory through the following four steps (see the figure):
Define the network nodes. These could be defined as electroencephalography or multielectrode-array electrodes, or as
Nature Reviews | Neuroscience
anatomically defined regions of histological, MRI or diffusion tensor imaging data.
Complex network
An informal description of a
network with certain
topological features, such as
high clustering,
small-worldness, the presence
of high-degree nodes or hubs,
assortativity, modularity or
hierarchy, that are not typical
of random graphs or regular
lattices. Most real-life networks
are complex by this definition,
and analysis of complex
networks therefore forms an
important methodological tool
for systems biology.
Adjacency matrix
An adjacency matrix indicates
the number of edges between
each pair of nodes in a graph.
For most brain networks, the
adjacency matrix is specified
as binary that is, each
element is either 1 (if there is
an edge between nodes) or 0
(if there is no edge). For
undirected graphs the
adjacency matrix is
symmetrical.
Estimate a continuous measure of association between nodes. This could be the spectral coherence or Granger causality
measures between two magnetoencephalography sensors, or the connection probability between two regions of an
individual diffusion tensor imaging data set, or the inter-regional correlations in cortical thickness or volume MRI
measurements estimated in groups of subjects.
Generate an association matrix by compiling all pairwise associations between nodes and (usually) apply a threshold to
each element of this matrix to produce a binary adjacency matrix or undirected graph.
Calculate the network parameters of interest in this graphical model of a brain network and compare them to the
equivalent parameters of a population of random networks.
Each step entails choices that can influence the final results and must be carefully informed by the experimental question.
At step 1, parcellation schemes can use prior anatomical criteria or be informed by the functional connectivity profiles of
different regions. Several such parcellation schemes may be available and can affect network measures147. In most magnetoencephalography and electroencephalography studies, network nodes are equivalent to individual electrodes or sensors,
but networks could also be based on reconstructed anatomical sources. However, some reconstruction algorithms will
determine the brain location of each source by minimizing the covariance between sensors, which has major effects on the
configuration of functional networks. At step 2, a range of different coupling metrics can be estimated, including measures
of both functional and effective connectivity. A crucial issue at step 3 is the choice of threshold used to generate an
adjacency matrix from the association matrix: different thresholds will generate graphs of different sparsity or connection
density, and so network properties are often explored over a range of plausible thresholds. Finally, at step 4 a large number of
network parameters can be quantified (BOX 2). These must be compared with the (null) distribution of equivalent parameters
estimated in random networks containing the same number of nodes and connections. Statistical testing of network
parameters may best be conducted by permutation- or resampling-based methods of non-parametric inference given the
lack of statistical theory concerning the distribution of most network metrics.
Most graph theoretical network studies to date have used symmetrical measures of statistical association or functional
connectivity such as correlations, coherence and mutual information to construct undirected graphs. This approach
could be generalized to consider asymmetrical measures of causal association or effective connectivity such as Granger
causal148,149 or dynamic causal66 model coefficients to construct directed graphs. It is also possible to avoid the
thresholding step (step 3) by analysing weighted graphs that contain more information than the simpler unweighted and
undirected graphs that have been the focus of attention to date. Structural brain network image is reproduced from
ReF. 59. Functional brain network image is reproduced, with permission, from ReF. 70 (2006) Society for Neuroscience.
naTuRe ReVIeWS | NeuroscieNce
VOlume 10 | maRcH 2009 | 187

REVIEWS
and vulnerability topics that are likely to become
increasingly relevant in relation to neuroscience.
In this article, we describe and discuss the expanding
interface between systems neuroscience and the physics
of complex networks. We review the existing empirical
data on topological and dynamical properties of structural and functional brain networks, and ask what this
literature tells us about how structural networks shape
functional brain dynamics. Space limitations prevent us
from providing coverage of all animal models, in vitro
Box 2 | Network measures

A network is defined in graph theory as a set of nodes or vertices and the edges or lines between them. Graph topology
can be quantitatively described by a wide variety of measures, some of which are discussed here. It is not yet established
which measures are most appropriate for the analysis of brain networks. The figure shows a schematic diagram of a brain
network drawn as a directed (left) and an undirected (right) graph; both structural and functional networks can be either
directed or undirected (BOX 1).
Node degree, degree distribution and assortativity

The degree of a node is the number of connections that link it to the rest of the network this is the most fundamental
network measure and most other measures are ultimately linked to node degree. The degrees of all the networks nodes
form a degree distribution15. In random networks all connections are equally probable, resulting in a Gaussian and
symmetrically centred degree distribution. Complex networks generally have non-Gaussian degree distributions, often
with a long tail towards high degrees. The degree distributions of scale-free networks follow a power law90. Assortativity is
the correlation between the degrees of connected nodes. Positive assortativity indicates that high-degree nodes tend to
connect to each other.
clustering coefficient and motifs
If the nearest neighbours of a node are also directly connected to each other they form a cluster. The clustering coefficient
quantifies the number of connections that exist between the nearest neighbours of a node as a proportion of the
maximum number of possible connections18. Random networks have low average clustering whereas complex networks
have high clustering (associated with high local efficiency of information transfer and robustness). Interactions between
neighbouring nodes can also be quantified by counting the occurrence of small motifs of interconnected nodes150. The
distribution of different motif classes in a network provides information about the types of local interactions that the
network can support48.
Path length and efficiency
Path length is the minimum number of edges that must be traversed to go from one node to another. Random and
complex networks have short mean path lengths (high global efficiency of parallel information transfer) whereas regular
lattices have long mean path lengths. Efficiency is inversely related to path length but is numerically easier to use to
estimate topological distances between elements of disconnected graphs.
connection density or cost
Connection density is the actual number of edges in the graph as a proportion of the total number of possible edges and is
the simplest estimator of the physical cost for example, the energy or other resource requirements of a network.
Hubs, centrality and robustness
Hubs are nodes with high degree, or high centrality. The centrality of a node measures how many of the shortest paths
between all other node pairs in the network pass through it. A node with high centrality is thus crucial to efficient
communication151. The importance of an
individual node to network efficiency can be
assessed by deleting it and estimating the
efficiency of the lesioned network.
Robustness refers either to the structural
integrity of the network following deletion of
nodes or edges or to the effects of
Threshold
perturbations on local or global network
states.
Modularity
Many complex networks consist of a number
of modules. There are various algorithms that
estimate the modularity of a network, many of
them based on hierarchical clustering23. Each
module contains several densely
interconnected nodes, and there are relatively
few connections between nodes in different
modules. Hubs can therefore be described in
terms of their roles in this community
structure27. Provincial hubs are connected
mainly to nodes in their own modules,
whereas connector hubs are connected to
nodes in other modules.
Motif
Path
Clustering
Connector hub
Community 2
Community 1
Provincial hub

188 | maRcH 2009 | VOlume 10
REVIEWS
Box 3 | Random, scale-free and small-world networks
In random graphs each pair of nodes has an equal probability, p, of being connected152.
Large random graphs have Gaussian degree distributions (BOX 2). It is now known that
most graphs describing real-world networks significantly deviate from the simple
random-graph model.
Some networks (including the Internet and the World Wide Web) have degree
distributions in the form of a power law: that is, the probability that a node has degree k
is given as Prob(k) ~ k. In biological systems, the degree exponent often ranges
between 2 and 3, and the very gradual (heavy-tail) power law decay of the degree
distribution implies that the network lacks a characteristic scale hence such
networks are called scale-free networks. Barabsi and Albert90 demonstrated that
scale-free networks can originate from a process by which each node that is added to
the network as it grows connects preferentially to other nodes that already have high
degree. Scale-free networks are unlikely if the attachment of connections is subject to
physical constraints or associated with a cost15. Therefore, physically embedded
networks, in which nodes have limited capacity for making connections, often do not
have pure power law degree distributions but may instead demonstrate exponentially
truncated power law degree distributions, which are associated with a lower
probability of very high degree nodes.
Originally described in social networks153, the small-world property combines high
levels of local clustering among nodes of a network (to form families or cliques) and
short paths that globally link all nodes of the network. This means that all nodes of a
large system are linked through relatively few intermediate steps, despite the fact that
most nodes maintain only a few direct connections mostly within a clique of
neighbours. Small-world organization is intermediate between that of random
networks, the short overall path length of which is associated with a low level of local
clustering, and that of regular networks or lattices, the high-level of clustering of which
is accompanied by a long path length18. A convenient single-number summary of
small-worldness is thus the ratio of the clustering coefficient to the path length after
both metrics have been standardized by comparing their values to those in equivalent
random networks154. Evidence for small-world attributes has been reported in a wide
range of studies of genetic, signalling, communications, computational and neural
networks. These studies indicate that virtually all networks found in natural and
technological systems have non-random/non-regular or small-world architectures
and that the ways in which these networks deviate from randomness reflect their
specific functionality.
Connectivity
In the brain, connectivity can
be described as structural,
functional or effective.
Structural connectivity denotes
a network of anatomical links,
functional connectivity denotes
the symmetrical statistical
association or dependency
between elements of the
system, and effective
connectivity denotes directed
or causal relationships
between elements.
Microcircuit
A neuronal network composed
of specific cell types and
synaptic connections, often
arranged in a modular
architecture and capable of
generating functional outputs.
Connectome
The complete description of
the structural connections
between elements of a nervous
system.
preparations and human studies that have contributed to

this endeavour: we necessarily focus on what we consider
to be some representative examples of graph theoretical
research in brain networks, with an emphasis on studies
of the human brain. We thus consider the implications of
complex brain networks for an understanding of neuropsychiatric disorders and conclude with some general
remarks about the evolution of scale-invariant topology
in brain networks and some key future challenges for this
emerging field.
Structural brain networks

Topological and physical properties of structural networks. The anatomical configuration of brain networks,
ranging from inter-neuronal connectivity to inter-regional
connectivity, has long been a focus of empirical neuroscience. network analysis, and in particular graph theory
(BOX 4), offers new ways to quantitatively characterize
anatomical patterns. according to graph theory, structural brain networks can be described as graphs that are
composed of nodes (vertices) denoting neural elements
(neurons or brain regions) that are linked by edges representing physical connections (synapses or axonal
projections). although graph theory emphasizes topological connectivity patterns, the topological and physical
distances between elements in brain networks are often

intricately related29. neurons and brain regions that are
spatially close have a relatively high probability of being
connected, whereas connections between spatially remote
neurons or brain regions are less likely 3032. longer axonal
projections are more expensive in terms of their material
and energy costs33. It has been suggested that the spatial
layout of neurons or brain regions34,35 is economically
arranged to minimize axonal volume. Thus, conservation of wiring costs is likely to be an important selection
pressure on the evolution of brain networks.
Mapping structural networks in animal models. Tracing
individual neuronal processes has been a long-standing
technological challenge, and few structural brain networks have been mapped at cellular resolution. Serial
reconstruction of electron microscopy sections allowed
the complete connection matrix of the nematode
Caenorhabditis elegans to be visualized36. currently,
this is the only nervous system to have been comprehensively mapped at a cellular level, and it was the first
to be described as a small-world network18. In larger
brains, combinations of physiological and anatomical
techniques have allowed patterns in neuronal connections to be identified, leading to the identification of
neuronal microcircuits37 and the formulation of probabilistic connection rules30. Small-world properties have
been demonstrated in biologically accurate models of
cellular networks in the reticular formation of the vertebrate brainstem38. Reconstruction of cellular networks in
the mammalian neocortex from multielectrode activity
recordings has revealed several highly nonrandom features of connectivity 39, including a tendency for synaptic
connections to be reciprocal and clustered. a promising
approach for mapping connectivity involves the stochastic expression of several fluorescent proteins40, and this
may ultimately deliver a complete map of the cellular
interconnections of an entire brain41.
Histological dissection and staining, degeneration
methods and axonal tracing have been used to map cerebral white matter connections. The pathways identified by these methodologies have formed the basis for
the systematic collation of species-specific anatomical
connection matrices, including those for the macaque
visual cortex 42 and the cat thalamocortical system43.
network analyses of such data sets demonstrated high
clustering of functionally related areas with short average
path lengths4446, hallmarks of a small-world architecture.
clusters identified by network analysis map on to known
functional subdivisions of the cortex 43. long-distance
cortical projections facilitate short-path communication despite increasing axonal volume47. The cortical
networks of several mammalian species also consistently
demonstrated an overabundance of motif classes associated with network modularity and functionally diverse
circuitry 48,49.
Mapping structural networks in the human brain.
Several attempts have been made to map the structural networks of the human brain, also known as the
human connectome50, at the scale of brain regions. One
VOlume 10 | maRcH 2009 | 189

REVIEWS
Box 4 | Origins of graph theory
Diffusion tensor imaging
(DTI). An MRI technique that
takes advantage of the
restricted diffusion of water
through myelinated nerve
fibres in the brain to map the
anatomical connectivity
between brain areas.
Diffusion spectrum imaging

An MRI technique that is
similar to DTI, but with the
added capability of resolving
multiple directions of diffusion
in each voxel of white matter.
This allows multiple groups of
fibres at each location,
including intersecting fibre
pathways, to be mapped.
Cortical parcellation
A division of the continuous
cortical sheet into discrete
areas or regions; Brodmanns
division of the cortex into areas
defined by their
cytoarchitectonic criteria is the
most famous but not the only
parcellation scheme.
Neuronographic
measurements
Recordings of epileptiform
electrical activity at specific
sites in the cortex following
topical application of a
pro-convulsive drug to a distant
cortical site; rapid propagation
of electrical activity from
stimulation to recording sites
implies that the sites are
anatomically connected.
Functional MRI
(fMRI). The detection of
changes in regional brain
activity through their effects on
blood flow and blood
oxygenation (which, in turn,
affect magnetic susceptibility
and tissue contrast in magnetic
resonance images).
Electroencephalography
(eeG). A technique used to
measure neural activity by
monitoring electrical signals
from the brain, usually through
scalp electrodes. eeG has good
temporal resolution but
relatively poor spatial
resolution.
Magnetoencephalography
(MeG). A method of measuring
brain activity by detecting
minute perturbations in the
extracranial magnetic field that
are generated by the electrical
activity of neuronal
populations.
Graph theory is rooted in the physical world155. In 1736, Euler showed that it was impossible to traverse the city of
Knigsbergs seven bridges across the river Pregel exactly once and return to the starting point. To prove this conjecture,
Euler represented the problem as a graph, and his original publication156 is generally taken to be the origin of a new
branch of mathematics called graph theory. In the middle of the nineteenth century, the analysis of electrical circuits and
the exploration of chemical isomers led to the discovery of additional graph theoretic concepts149. Today, graph theory
pervades many areas of science.
Significant progress in graph theory has come from the study of social networks157. One prominent experiment153
tracked paths of acquaintanceship across a large social network and found that even very large networks could be
traversed, on average, in a small number of steps. Although this small-world phenomenon quickly captured the public
imagination, its origins remained obscure until its association with specific types of connectivity was demonstrated
(BOX 3). The dual discoveries of small-world18 and scale-free90 networks launched the modern era of graph theory, which
now extends into biology and neuroscience. Neural systems have long been described as sets of discrete elements linked
by connections. Nonetheless, graph theory has essentially only been applied to neuroscience in the past 10 years.
study derived structural connection patterns from crosscorrelations in cortical thickness or volume across
individual brains, which might indirectly indicate the
presence of corticocortical pathways51,52. Graph analysis
revealed small-world attributes and the existence of local
communities of brain regions. a more detailed analysis53
of the modularity or community structure of connection
data sets derived from cortical thickness correlations
revealed significant overlap between anatomical network
modules and functional systems in the cortex.
Human brain structural networks have also been
mapped using diffusion imaging and tractography. a map
of 7090 cortical and basal brain grey matter areas was
constructed using diffusion tensor imaging (DTI) and analysed using methods from graph theory 54,55. The network
exhibited high clustering and short path length, and contained motif classes similar to those identified from tracttracing data48. Several areas, including the precuneus, the
insula, the superior parietal cortex and the superior frontal
cortex, were found to have high betweenness centrality
and thus to constitute putative hubs. another study that
mapped connections between 78 cortical regions using
DTI also identified several hub regions, including the
precuneus and the superior frontal gyrus56.
Due to limitations in the model that is used to infer
fibre bundle orientation, DTI has difficulty detecting
crossing fibre bundles. Diffusion spectrum imaging can overcome this limitation by reconstructing multiple diffusion
directions in each voxel57 and was used to build cortical
connection matrices between 5004,000 homogeneously
distributed regions of interest 58. again, network analyses
revealed the small-world architecture of the cortical network. a more extensive analysis of 998 region-of-interest
networks obtained from 5 participants 59 identified
structural modules interconnected by highly central hub
regions. When considering multiple network measures
(including node degree, connection strength and centrality), a particular set of brain regions located predominantly in the posterior medial cortex, including portions
of the posterior cingulate and the precuneus, was highly
and densely interconnected, forming a structural core59.
although they differ in terms of their experimental methodology and cortical parcellation, most of these
studies reveal highly clustered large-scale cortical
networks, with most pathways existing between areas
that are spatially close and functionally related. These

clusters or modules are interlinked by specialized hub
regions, ensuring that overall path lengths across the
network are short. most studies identified hubs among
parietal and prefrontal regions, providing a potential
explanation for their well-documented activation
by many cognitive functions. Particularly notable is
the prominent structural role of the precuneus 55,56,59,
a region that is homologous to the highly connected
posteromedial cortex in the macaque60. The precuneus
is involved in self-referential processing, imagery and
memory 61, and its deactivation is associated with anaesthetic-induced loss of consciousness62. an intriguing
hypothesis suggests that these functional aspects can
be explained on the basis of its high centrality in the
cortical network.
Functional brain networks

although analysing structural networks helps us to
understand the fundamental architecture of interregional connections, we must also consider functional
networks directly to elucidate how this architecture supports neurophysiological dynamics. Despite considerable
heterogeneity in the methodological approaches, there is
an encouraging degree of convergence between studies
of functional brain networks. The first such study used a
set of neuronographic measurements of the propagation of
epileptiform activity following localized applications
of strychnine to the macaque cortex 63. This demonstrated a pattern of functional connections between
cortical areas that was consistent with a small-world
network. as we discuss below, these findings have been
extended by studies based on functional MRI (fmRI),
electroencephalography (eeG), magnetoencephalography
(meG) or multielectrode array (mea) data.
although such studies based on graph theory are
the focus of this Review, we note that other methods to
investigate brain functional systems have recently been
developed, including mathematical models of effective connectivity between regions. effective connectivity models, such as structural equation modelling 64,65,
dynamic causal modelling 66 or Granger causality 67,
involve estimating the causal influence that each element of a system exerts on the behaviour of other
elements. Thus, measures of effective connectivity
190 | maRcH 2009 | VOlume 10
REVIEWS
between multiple regions can be used to generate a
directed graph, which can then be topologically described
using graph theory. However, the functional network
studies reviewed below have all been based on undirected graphs, derived from simpler measures of functional connectivity or symmetrical statistical association
between brain regions. The key point is that, in principle,
graph theory could be applied to an association matrix of
either functional or effective connectivity measures, to
generate either undirected or directed graphs, respectively, although all neuroscientific studies to date have in
fact been based on measures of functional connectivity.
Multielectrode array
(MeA). A technique for
simultaneously measuring the
electrical activity of local
neuronal populations or single
neurons, usually in tissue slices
or cell cultures in vitro.
Association matrix
A matrix that represents the
strength of the association
between each pair of nodes in
a graph. Association between
nodes can be quantified by
many continuously variable
metrics, such as correlation or
mutual information. either
functional or effective
connectivity measures can be
used to construct an
association matrix.
Blood oxygen
level-dependent (BOLD)
signals
Changes in magnetic
susceptibility and MRI tissue
contrast that are indirectly
indicative of underlying
changes in spontaneous or
experimentally controlled brain
activation.
Default-mode network
A set of brain regions, including
medial frontal and posterior
cingulate areas of the cortex,
that are consistently
deactivated during the
performance of diverse
cognitive tasks.
Mapping functional networks using fMRI. The first

graph theoretical study of fmRI data measured the
partial correlations of resting-state blood oxygen leveldependent (BOLD) signals between 90 cortical and subcortical regions and reported small-world properties of
the resulting whole-brain networks68. almost simultaneously, another study reported small-world properties
of functional networks derived from a set of activated
voxels in fmRI data; this voxel-level network was also
reported to have a scale-free degree distribution69.
Subsequently, small-world properties, with parameter
values similar to those previously reported in topological studies of cat and macaque anatomical connectivity
matrices50, were confirmed in a low-frequency (0.03
0.06 Hz) whole-brain network derived from wavelet
correlations between regional mean time series70. The
high-degree nodes or hubs of this network were mostly
regions of multimodal association cortex, and the degree
distribution was an exponentially truncated power law 70.
Other studies have explored the community structure of
fmRI networks using a hierarchical cluster analysis68,71,72
and shown that functionally and/or anatomically related
brain regions are more densely interconnected, with relatively few connections between functional clusters, again
echoing prior work on anatomical connectivity matrices.
The high density of connections between functionally
related regions increases the clustering coefficient of
the graph, whereas the long-range connections between
different modules or clusters, even though they are relatively few in number, keep the path length low. Thus, the
small-world architecture of a brain functional network
is closely related to its modularity 72.
There are several other metrics for quantifying smallworld architecture in brain functional networks. Studies
in statistical physics73,74 have shown that path length is
inversely related to the global efficiency of a network for
the transfer of information between nodes by multiple
parallel paths, and that global efficiency is easier to estimate than path length when studying sparse networks.
Furthermore, the clustering coefficient can be regarded as
a measure of the local efficiency of information transfer,
or of the robustness of the network to deletion of individual nodes. The structural network of the macaque brain
was found to have high global and local efficiency and to
be sparsely connected73,74. Thus, the macaque cortex has
economical small-world properties: it has high global
efficiency of parallel information transfer and high local
fault tolerance for relatively low connection density.
These concepts were translated to the analysis of

resting-state fmRI data acquired from young and elderly adults75, using the wavelet correlation (a measure
of the association between time series in a specific frequency band) to estimate the functional connectivity
between regional BOlD time series endogenously oscillating in the frequency interval 0.060.1 Hz (for a more
detailed review of the rationale for wavelet analysis in
fmRI, see ReF. 76). In younger adults, functional brain
networks demonstrated small-world properties over a
broad range of connection densities or costs. Relatively
sparse networks were associated with maximum cost
efficiency. The older age group also showed evidence of
small-world properties, but had significantly reduced
cost efficiency: they had to be relatively over-connected
to provide the efficiency of parallel information transfer seen in a younger brain network. The suggestion
that aging is associated with changes in the economical small-world properties of brain functional networks
converges with studies that show differences in attentional and default-mode networks between children and
young adults77. normal processes of brain maturation
and senescence might thus be reflected in quantifiable
changes in functional network topology.
Mapping functional networks using electrophysiological techniques. When comparing the results of fmRI
studies to results obtained using electrophysiological
techniques (eeG, meG or mea), many aspects of the
data clearly differ. fmRI has good spatial resolution (on
the order of millimetres) but poor temporal resolution
(on the order of seconds), restricting the measurable
bandwidth to approximately 0.0010.5 Hz, and fmRI
measures activation-related haemodynamics rather than
neuronal activity per se. all electrophysiological methods measure neuronal activity more directly and have
better temporal resolution, with bandwidths typically of
1100 Hz, but they often have worse spatial resolution
(on the order of millimetres or centimetres for meG and
eeG) or less complete anatomical coverage (in the case
of mea) than fmRI. another point of difference is that
the nodes of a network derived from fmRI data will be
anatomically localized regions or voxels of the image,
whereas the nodes of a network derived from meG or
eeG data could be the surface sensors or recording electrodes. However, we can compare the topologies of networks derived from these different data sets by invoking
a general operating principle of complex network analysis: microscopically distinct systems can be informatively
compared in terms of their macroscopic organization
using graph theory.
Functional connectivity between pairs of electrodes
has been estimated using a measure of generalized synchronization78, and then thresholded to generate functional networks, in several studies of eeG or meG data
sets. This has shown that small-world topology is represented at many frequency intervals79 and can be related
to cognitive performance and normal aging 80. an alternative approach used the wavelet correlation to estimate
frequency-dependent functional connectivity between
meG sensors, again revealing that many topological and
VOlume 10 | maRcH 2009 | 191

REVIEWS
dynamic properties of brain functional networks were
conserved across frequencies81. This timescale invariance of wavelet correlations and the brain functional
network parameters derived from them is a theoretically
predictable corollary of the long-range autocorrelations
of neurophysiological time series8284.
Metastable dynamics
Transitions between marginally
stable network states; these
transitions can occur
spontaneously or as a result of
weak external perturbations.
Resting state
A cognitive state in which a
subject is quietly awake and
alert but does not engage in or
attend to a specific cognitive
or behavioural task.
Conservation of functional network properties. There

have been fewer studies of functional networks in nonhuman species. In anaesthetized rats, fmRI has been
used to demonstrate small-world and modular properties of whole-brain functional networks85. a recent
study considered functional coupling between pairs of
cortical neurons, measured by multiunit electrodes, in
the visual cortex of the anaesthetized cat 86. most of the
properties of these neuronal-interaction networks could
be accounted for by the pairwise correlations between
electrodes87,88, and the networks had small-world properties with some highly connected (hub) neurons. The
small-world organization of functional networks at a cellular level has also been described on the basis of mea
recordings from in vitro cultures of cortical networks89.
These results indicate that small-worldness might be a
conserved property of brain functional networks over
different species and spatial scales.
Preliminary evidence suggests that other topological properties of brain functional networks, such as the
degree distribution, might also be conserved over spatial
and temporal scales. many large networks exhibit scalefree power law degree distributions90 indicative of the
existence of highly connected nodes (BOX 3). Some studies of functional brain networks carried out at high spatial
resolution (single voxels in fmRI) have provided evidence
of scale-free organization69,91. However, in other studies of
functional networks derived from fmRI and meG data
sets over a wide range of frequency intervals70,81, the
empirical degree distribution conforms to an exponentially truncated power law, implying that the probability of a highly connected node or hub is greater than
in an equivalent random network but less than would
be expected in a scale-free network. Truncated power
law degree distributions have also been reported from
analysis of anatomical networks derived from structural
mRI data in humans55, and from analysis of functional
networks derived from mea data in cats86. These distributions are typical of physically embedded networks,
such as the global air-transportation network, in which
the maximum degree is limited by physical considerations such as the finite capacity of any node to receive
connections27. The reasons for the differences between
these findings and reports of scale-free properties are
currently unknown; however, it is notable that pure
power law scaling of the degree distributions of human
brain functional networks has only been reported by
voxel-level analysis, whereas exponentially truncated
power laws have been reported by region-level analysis.
The form of the degree distribution could be affected
by spatial interpolation, and other pre-processing steps
applied before the construction of functional networks
and the standardization of such methodological procedures will be important in elucidating the impact of
anatomical resolution on the extent to which power

law scaling of the degree distribution is truncated at
high degree.
Structurefunction relations in brain networks

How do functional brain networks emerge from structural brain connectivity? Structural maps indicate that
each neural node maintains a specific pattern of structural
connections with other nodes. Different nodes often have
different functionalities, such as specific response preferences to sensory stimuli. From a network perspective, the
functionality of an individual neural node is partly determined by the pattern of its interconnections with other
nodes in the network92: nodes with similar connection
patterns tend to exhibit similar functionality 36. although
functional properties are expressed locally, they are the
result of the action of the entire network as an integrated
system. Structural connectivity places constraints on
which functional interactions occur in the network.
Structural and functional connectivity in cellular networks undergo dynamic changes. The degree to which
synaptic connectivity is modified in the adult brain is
highly debated. although some evidence suggests that
mammalian cellular networks are continually remodelled93, other evidence indicates that most synaptic
spines are stable94. changes in neuronal connectivity
necessitate homeostatic mechanisms to ensure functional stability 95. multielectrode recording data suggest
that cellular functional networks exhibit transient synchronization96 and metastable dynamics97. These changes
occur within seconds, and it seems possible that even
more rapid transitions and network reconfigurations
may take place33. These observations of relatively slow
structural modifications accompanied by faster changes
in functional linkages pose major unresolved questions
regarding functional stability in neural circuits.
It is currently unknown whether large-scale cortical
networks in the adult brain undergo structural modifications on fast timescales. most of the changes that
have been observed were associated with aging, disease
progression or experience-dependent plasticity. By contrast, patterns of functional connectivity between brain
regions undergo spontaneous fluctuations and are
highly responsive to perturbations, such as those that
are induced by sensory input or cognitive tasks, on a
timescale of hundreds of milliseconds. These rapid reconfigurations do not affect the stability of global topological characteristics81,98. On longer timescales of seconds
to minutes, correlations between spontaneous fluctuations in brain activity 99101 form functional networks that
are particularly robust. For example, a set of posterior
medial, anterior medial and lateral parietal brain regions
comprises the default mode network102,103.
The persistence of functional networks associated
with the brains resting state provides an opportunity
to investigate how much of the pattern of functional
connections is determined by underlying structural
networks. Observations from a single cortical slice104,
structural imaging of fibre bundles linking components
of the default network105, and direct comparisons of
structural and functional connectivity in the same cohort
192 | maRcH 2009 | VOlume 10
REVIEWS
Structural brain network
Anatomical
tract tracing
Coupling matrix and neural dynamics

Simulated
electrophysiological
time series
Simulated BOLD
time series
Functional brain network
Mutual information
Cross-correlation
Figure 1 | computational modelling of structural and functional brain networks.

Computational models have been used to demonstrate how dynamic patterns arise as a
result of interactions between anatomically connected neural units. Shown is how such
a model is generated and used. A structural brain network derived from anatomical data
serves as a matrix of coupling coefficients that link neuronal nodes, the activities of which
unfold through time. This time evolution is governed by physiologically motivated dynamic
equations. In the example shown, the surface of the macaque cortex was subdivided into
47 areas (nodes) and a structural brain network linking these nodes was compiled from
anatomical tract-tracing data. The dynamic equations were derived from a model of large
neuronal populations, the parameters of which were set to physiological values109. Data
from computer simulations then yield functional brain networks. Such networks are
derived from measures of association between the simulated time series for example, an
information theoretic measure such as the mutual information (computed on
voltagetime data) or cross-correlations in neural activity that are computed from
simulated blood oxygen level-dependent (BOLD) data. These matrices can then be
thresholded to yield binary networks from which network measures can be derived. The
fact that both structural and functional networks are completely specified in the model
facilitates their comparative analysis. The structural brain network panel is reproduced,
with permission, from ReF. 109 (2007) National Academy of Sciences. The rest of the
figure is modified, with permission, from ReF. 158 (2009) Academic Press.
of participants63,106,107 suggest that structural connections

are highly predictive of functional connections. Indirect
interactions can account for additional functional linkages. Such indirect connections can lead to discrepancies between structural and functional connectivity;
however, current evidence suggests that topological
parameters are generally conserved between structural
and functional networks.
computational models offer a complementary
method to investigate structurefunction relations in
brain networks. Several studies have used empirically
derived structural brain networks to specify the wiring
between neural units, which then results in spatially

patterened dynamic interactions (FIG. 1). The covariance structure of the endogenous activity generated by
these neural units yields a functional network that can
be generated by computer simulations or in some cases
analytically 108. This strategy has been used to determine
the effects of structural topology on functional networks
and dynamics. Simulation studies of large-scale cortical networks demonstrated the emergence of complex
spatiotemporal structure in neural correlations at multiple timescales109, as well as realistic patterns of modelled BOlD resting-state correlations that depend on
the topology 109 and time delays110 in the structural coupling matrix. The topology of structural and functional
networks was identical when functional connectivity
was estimated from long time samples, but functional
networks estimated on shorter time samples or at higher
frequencies were less strongly constrained by the structural wiring diagram109. likewise it has been shown
that the modularity of structural networks can determine the hierarchical organization of functional networks111,112 and may be important for generating diverse
and persistent dynamic patterns113. In a computational
model of phase synchronization between coupled neurons, patterns of local versus global synchronization
were found to depend strongly on the balance between
high clustering and short path length in the connectivity of the network114. computational studies have also
begun to determine the effects of functional activity on
structural topology. It has been shown that an initially
random wiring diagram can evolve through synaptic
plasticity to a functional state characterized by a smallworld topology of the most strongly connected nodes
and by self-organized critical dynamics115.
These findings indicate that the brains structural
and functional networks are intimately related and
share common topological features, such as modules
and hubs (FIG. 2). although most studies provide support for the idea that structural networks determine
some aspects of functional networks, especially at low
frequencies or over long time periods, it is less clear how
the structural topology both supports the emergence of
fast and flexibly reconfigured functional networks and
is itself remodelled by function-related plasticity on a
slower timescale.
Clinical and translational aspects

Since the work of pioneers such as Wernicke, meynert
and Dejerine, it has been appreciated that many neurological and psychiatric disorders can be described as dysconnectivity syndromes116. The emergence of symptoms
or functional impairment in these disorders can be theoretically related to the disruption or abnormal integration
of spatially distributed brain regions that would normally
constitute a large-scale network subserving function.
Using network properties as diagnostic markers. One
application of complex-network theory in this context is
to provide new measures to quantify differences between
patient groups and appropriate comparison groups.
Several studies have reported that the parameters of
VOlume 10 | maRcH 2009 | 193

REVIEWS
Cellular functional network
Module 1
Hub node
Module 2
PIP
46
VP
DP
CITd
VOT
VIP
LIP
MT
V2
PO
V4
V3
PITv
TF
FEF
AITv
V3A
V1
TH
FST
V4t
Whole-brain structural network
PITd
CITv
Figure 2 | cellular and whole-brain networks demonstrate consistent topological

Reviews
| Neuroscience
features. The top panel shows a cellular functional networkNature
constructed
from
multielectrode-array recordings made in the anaesthetized cat; each node (represented by a
circle) corresponds approximately to one neuron and the connections represent high
functional connectivity between neurons86. The different coloured nodes constitute
separate clusters or modules. The plots in each circle illustrate cellular responses to
stimuli of different orientations, and the circle size corresponds to the degree (number of
functional connections) of each node. The bottom panel shows a whole-brain structural
network constructed from histological data on the macaque cortex; each node
corresponds to a brain area and the connections represent axonal projections between
areas49. The network has two main modules, shown here with yellow and grey circles
corresponding to mostly dorsal and ventral visual regions, respectively. Both networks
exhibit the small-world attributes of high clustering and short path length (see BOX 3);
both have an exponentially truncated power law degree distribution (see BOX 2),
associated with the existence of high-degree hubs (V4 in the anatomical network); and
both have a community structure characterized by sparse connectivity between modules
(each module is enclosed by stippled lines) and linked by hubs (nodes circled in red).
AITv, anterior inferotemporal ventral area; CITd, central inferotemporal dorsal area;
CITv, central inferotemporal ventral area; DP, dorsal preluneate area; FEF, frontal eye
field; FST, floor of superior temporal area; LIP, lateral intraparietal area; MT, middle
temporal area; PIP, posterior intraparietal area; PITd, posterior inferotemporal dorsal
area; PITv, posterior inferotemporal ventral area; PO, parieto-occipital area; TF, area TF;
TH, area TH ;V14, visual cortical areas 14; VIP, ventral intraparietal area; VOT, ventral
occipitotemporal area; VP, ventral posterior area. The top panel is reproduced, with
permission, from ReF. 86 (2008) Oxford University Press. The bottom panel is
reproduced from ReF. 49.
brain networks derived from fmRI, eeG or structural

mRI data are altered in patients with schizophrenia or
alzheimers disease (aD).
In an fmRI study of aD, clustering was significantly
reduced at a global level (in whole-brain networks
operating at frequencies below 0.1 Hz) and at a local
level (in both hippocampi), and global clustering was
able to discriminate aD patients from age-matched
comparison subjects with high specificity and sensitivity, implying that the loss of small-world network
properties might provide a clinically useful diagnostic
marker 117. In a comparable eeG study, path length in
beta band (1535 Hz) functional networks was significantly increased in patients with aD. Importantly, the
variability of cognitive function across both control and
aD groups was negatively correlated with path length,
providing direct evidence that functional-network
topology can be related to variation in cognitive performance118. a third meG study of resting-state functional networks confirmed a degradation of small-world
attributes in patients with aD and suggested that this
effect is due to disease-related changes at highly connected network hubs119. However, in a fourth study that
used between-subject covariation in regional measures
of cortical thickness to infer anatomical networks from
a large structural mRI data set, global clustering was
significantly increased in patients with aD and there
were abnormalities in the topological configuration
of crucial, high-centrality nodes in regions of the
multimodal association cortex 120.
Schizophrenia has been investigated by comparable
approaches. The economical small-world properties of
low-frequency functional networks derived from fmRI
data were shown to be impaired in patients with schizophrenia121. Studies that used eeG to measure synchronization likelihood or nonlinear interactions between
cortical nodes found that the clustering and path-length
parameters of functional networks were closer to their
values in random graphs for patients with schizophrenia
than for comparison subjects122,123. anatomical networks
derived from inter-regional covariation of grey matter density in structural mRI data showed differences
in the hierarchy and assortativity of multimodal and
transmodal cortical subnetworks in healthy volunteers,
suggesting that these major divisions of the cortex may
have developed according to different growth rules or
evolved to meet different selection criteria. Topological
abnormalities in people with schizophrenia included
a reduced hierarchy of multimodal cortex (FIG. 3). The
schizophrenic groups networks also exhibited relatively
long physical distances between connected regions,
compatible with inefficient axonal wiring 124.
Thus, there is convergent evidence from methodologically disparate studies that both aD and schizophrenia are associated with abnormal topological
organization of structural and functional brain networks. However, there are also inconsistencies between
existing studies for example, clustering is reportedly
increased in the structural networks120 but decreased in
the functional networks of patients with aD117 which
might be attributable to the clinical heterogeneity of the
194 | maRcH 2009 | VOlume 10
REVIEWS
a Healthy volunteers
b People with schizophrenia

37
Frontal
Temporal
Parietal
Occipital
High clustering
Low clustering
44'
45'
46
12
37
20
36'
Degree
8
22
34'
44
21'
19
36
46'
47'
44
22
21
10
43'
46
39
43
39'
22
42'
47'
40
47
20
8'
9'
40
44'
10'
40'
20'
43
45'
39'
45
10'
20
37'
39
40'
19'
37'
19'
47
10
8'
45
36
43'
19
9'
21'
22'
Figure 3 | Disease-related disorganization of brain anatomical networks derived

from structural Mri data. In both parts, the nodes (circles) represent cortical regions
Nature
Reviews
| Neuroscience
and the connections represent high correlation in grey matter
density
between
nodes.
The nodes are arranged vertically by degree and are separated horizontally for clarity of
representation. The numbers indicate approximate Brodmann area, and the prime
symbols () denote left-sided regions. The clustering coefficient of each node, a measure
of its local connectivity, is indicated by its size: nodes with high clustering are larger.
a | The brain anatomical network of the healthy volunteers has a hierarchical
organization characterized by low clustering of high-degree nodes24. b | The equivalent
network constructed from MRI data on people with schizophrenia shows loss of this
hierarchical organization high-degree nodes are more often highly clustered. Figure is
reproduced, with permission, from ReF. 124 (2008) Society for Neuroscience.
Assortativity
A measure of the tendency for
nodes to be connected to
other nodes of the same or
similar degree.
Endophenotype
A quantifiable biological
marker of the genetic risk for a
neuropsychiatric disorder.
patient groups as well as to the differences in imaging and analytic methods. Some of these differences
may perhaps be resolved by studies combining network
measurements on structural and functional neuroimaging data acquired on the same patients. It also seems
likely that evidence for network abnormalities in other
neuropsychiatric disorders and conditions (such as epilepsy 125127, attention-deficit hyperactivity disorder 128 or
spinal cord injury 129) will accumulate as the disorders
are increasingly investigated from this perspective.
Understanding the pathogenesis and treatment of brain
disorders from a network perspective. many psychiatric
disorders are highly heritable and are likely to represent the clinical outcome of aberrations in the formation of large-scale networks in utero or during early
postnatal life. measures of network topology may be
worth investigating as intermediate phenotypes, or

endophenotypes , that indicate the genetic risk for a
neuropsychiatric disorder; however, network metrics
have not yet been adopted for this purpose. a study of
healthy twin pairs has shown that classical small-world
metrics on brain functional networks derived from eeG
data have high heritability 130, a necessary prerequisite for
their candidacy as disease endophenotypes. another
study of graph theoretical measures of anatomical networks derived from inter-regional correlations in cortical-thickness mRI measurements on a sample of normal
twins, singletons and singleton siblings of twins showed
that genetically determined frontoparietal networks had
small-world properties131. network metrics are arguably
more attractive as intermediate phenotypes than local
measures of brain (dis)organization, because computational models of network development are often available to test mechanistic hypotheses for how an observed
profile of anatomical or functional dysconnectivity in a
mature network might have been generated by earlier
developmental abnormalities24,132.
another example of how empirical and computational
approaches can be usefully combined is provided by studies that have lesioned anatomical or functional network
models for example, by deleting nodes or connections
to explore how acute and focal damage could affect the
overall performance of brain networks70,133,134. networks
can be lesioned by random deletion of nodes or edges,
or by targeted attack on the highest-degree nodes in the
network. The vulnerability of the network to damage is
assessed by comparing its topological or dynamical behaviour after the lesioning to its intact behaviour. Different
network topologies confer different vulnerabilities to the
effects of random or targeted attack. For example, scalefree networks are robust to random error but highly
vulnerable to deletion of the network hubs. Brain functional networks with an exponentially truncated power
law degree distribution were found to be less vulnerable
to attack than scale-free networks70. In an anatomically
informed computational model, deletion of hub nodes
produced widespread disruptions of functional connectivity 53,134 that were consistent with effects reported in
focal human brain lesions135,136. computational lesioning
of network models was also used to explore the functional consequences of a gradual and precisely specified
disease process: the elimination of long-range projections
and the sprouting of short-range connections in a model
of epileptogenesis in the rat dentate gyrus137. The topology of the normal or non-epileptic dentate gyrus became
relatively over-connected and dynamically hyperexcitable
as a result of cellular changes previously described in relation to temporal lobe epilepsy. Other studies of models
of temporal lobe epilepsy have shown loss of small-world
topology in cellular networks during hypersynchronized
bursting 138 and have shown that variation of small-world
topological and synaptic properties of a computational
model can cause transitions between normal, bursting
and seizing behaviours139.
It is also conceivable that network analysis can be
used to further our understanding of the therapeutic
effects of pharmacological or psychological therapies.
VOlume 10 | maRcH 2009 | 195

REVIEWS
Dopaminergic drugs can modulate measures of functional
connectivity in animal and human fmRI140,141 and meG
recordings142, suggesting that drug effects might be quantifiable in terms of altered functional network topology.
This has been confirmed directly in a study which demonstrated that a dopamine D2 receptor antagonist impaired
the economical small-world properties of human brain
fmRI networks75. Future work might include efforts
using graph theoretical measures to quantify how therapeutically effective treatments remediate topologically
sub-optimal network configurations in patients.
Conclusions and prospects

It is clear that certain aspects of the organization of complex brain networks are highly conserved over different
scales and types of measurement, across different species
and for functional and anatomical networks. The archetypal brain network has a short path length (associated
with high global efficiency of information transfer), high
clustering (associated with robustness to random error),
a degree distribution compatible with the existence of
hubs, and a modular community structure. Furthermore,
anatomical networks are sparsely connected, especially
between nodes in different modules, and the wiring
length (the physical distance that connections span) is
close to minimal. This profile of topological and geometric properties is typical not just of brain networks but
also of many other complex networks, including transport systems and intracellular signalling pathways17,73.
Why might this be so?
a parsimonious hypothesis is that many spatially
embedded complex networks have evolved to optimize
the same set of competitive selection criteria high efficiency of information transfer between nodes at low connection cost or to achieve an optimal balance between
functional segregation and integration that yields high
complexity dynamics14. If wiring cost was exclusively
prioritized the network would be close to a regular lattice, whereas if efficiency was the only selection criterion
the network would be random. The existence of a few
long-range anatomical connections can deliver benefits
Cajal, S. R. Histology of the Nervous System of Man

and Vertebrates (Oxford Univ. Press, New York, 1995).
Swanson, L. W. Brain Architecture (Oxford Univ. Press,
Oxford, 2003).
3.
Singer, W. Neuronal synchrony: a versatile code for the
definition of relations? Neuron 24, 4965 (1999).
4.
Fries, P. A mechanism for cognitive dynamics: neuronal
communication through neuronal coherence. Trends
Cogn. Sci. 9, 474480 (2005).
5.
Bressler, S. L. Large-scale cortical networks and
cognition. Brain Res. Brain Res. Rev. 20, 288304
(1995).
6.
Mesulam, M. M. From sensation to cognition. Brain
121, 10131052 (1998).
7.
McIntosh, A. R. Towards a network theory of cognition.
Neural Netw. 13, 861870 (2000).
8.
Friston, K. Beyond phrenology: what can
neuroimaging tell us about distributed circuitry?
Annu. Rev. Neurosci. 25, 221250 (2002).
9.
Buzski, G. Rhythms of the Brain (Oxford Univ. Press,
New York, 2006).
10. Strogatz, S. H. Exploring complex networks. Nature
410, 268277 (2001).
11. Albert, R. & Barabsi, A. L. Statistical mechanics of
complex networks. Rev. Mod. Phys. 74, 4797
(2002).
1.
2.
12.
13.
14.
15.
16.
17.
18.
in terms of efficiency and could arguably account for

the evolution of economical small-world properties in
brain networks at all scales20,51. This hypothesis needs to
be more directly explored and tested, perhaps using evolutionary algorithms in computational models of brain
network selection48.
a key issue for the future will be to consolidate our
understanding of how functional networks interact with
their structural substrates. at low frequencies, or over long
time periods, there are reasons to expect that functional
networks should be highly isomorphic with underlying
structural networks84,109. But clearly function can be adaptive over much shorter timescales than structure. We need
to understand more about the non-stationarity or metastability 143 of brain functional networks. How does functional network topology change over time? Do functional
networks exist in a dynamically critical state at some or all
frequency intervals144146? What constraints on the itinerancy of network dynamics are imposed anatomically and
how does the long-term history of functional activity in a
network feed back on the development and remodelling
of the anatomical connections between nodes?
a related question concerns how the parameters of
complex brain networks relate to cognitive and behavioural functions. One can make an intuitively reasonable claim that high clustering favours locally specialized
processing whereas short path length favours globally
distributed processing; but the empirical evidence is currently almost non-existent. This will probably be a key
focus of future work that might be combined with further studies of clinical disorders or cohorts at different
stages of normal development.
The emerging field of complex brain networks raises
a number of interesting questions and provides some of
the first quantitative insights into general topological
principles of brain network organization. The fundamental growth in the statistical mechanics of complex
networks, and the power and elegance of graph theoretical analysis, suggests that this approach will play an
increasingly important part in our efforts to comprehend
the physics of the brains connectome.
A scholarly review of the early literature on the

physics of complex networks, with an emphasis on
various types of scale-free and small-world
connectivity.
Boccaletti, S., Latora, V., Moreno, Y., Chavez, M. &
Hwang, D.-U. Complex networks: structure and
dynamics. Phys. Rep. 424, 175308 (2006).
Brner, K., Sanyal, S. & Vespignani, A. Network
science. Annu. Rev. Inform. Sci. Technol. 41, 537607
(2007).
Tononi, G., Sporns, O. & Edelman, G. M. A measure
for brain complexity: relating functional segregation
and integration in the nervous system. Proc. Natl
Acad. Sci. USA 91, 50335037 (1994).
Amaral, L. A. N., Scala, A., Barthelemy, M. & Stanley,
H. E. Classes of small-world networks. Proc. Natl Acad.
Sci. USA 97, 1114911152 (2000).
Amaral, L. A. N. & Ottino, J. M. Complex networks.
Augmenting the framework for the study of complex
systems. Eur. Phys. J. B 38, 147162 (2004).
Barabsi, A. L. & Oltvai, Z. N. Network biology:
understanding the cells functional organization.
Nature Rev. Genet. 5, 101113 (2004).
Watts, D. J. & Strogatz, S. H. Collective dynamics of
small-world networks. Nature 393, 440442
(1998).
196 | maRcH 2009 | VOlume 10
19.
20.
21.
22.
23.
24.
25.
26.
This seminal paper on small-world networks

demonstrated their ubiquitous occurrence in
natural, social and technological systems.
Sporns, O., Chialvo, D., Kaiser, M. & Hilgetag, C. C.
Organization, development and function of complex
brain networks. Trends Cogn. Sci. 8, 418425 (2004).
Bassett, D. S. & Bullmore, E. T. Small world brain
networks. Neuroscientist 12, 512523 (2006).
Reijneveld, J. C., Ponten, S. C., Berendse, H. W. &
Stam, C. J. The application of graph theoretical
analysis to complex networks in the brain. Clin.
Neurophysiol. 118, 23172331 (2007).
Stam, C. J. & Reijneveld, J. C. Graph theoretical
analysis of complex networks in the brain. Nonlin.
Biomed. Phys. 1, 3 (2007).
Girvan, M. & Newman, M. E. J. Community structure
in social and biological networks. Proc. Natl Acad. Sci.
USA 99, 78217826 (2002).
Ravasz, E. & Barabsi, A. L. Hierarchical organization
in complex networks. Phys. Rev. E Stat. Nonlin. Soft
Matter Phys. 67, 026112 (2003).
Barthlemy, M. Betweenness centrality in large complex
networks. Eur. Phys. J. B 38, 163168 (2004).
Guimer, R. & Amaral, L. A. N. Functional cartography
of complex metabolic networks. Nature 433, 895900
(2005).
REVIEWS
27. Guimer, R., Mossa, S., Turtschi, A. & Amaral, L. A.
The worldwide air transportation network: anomalous
centrality, community structure and cities global roles.
Proc. Natl Acad. Sci. USA 102, 77947799 (2005).
28. Kashtan, N. & Alon, U. Spontaneous evolution of
modularity and network motifs. Proc. Natl Acad. Sci.
USA 102, 1377313778 (2005).
29. Laughlin, S. B. & Sejnowski, T. J. Communication in
neuronal networks. Science 301, 18701874 (2003).
30. Braitenberg, V. & Schz, A. Statistics and Geometry of
Neuronal Connectivity (Springer, Berlin, 1998).
31. Hellwig, B. A quantitative analysis of the local
connectivity between pyramidal neurons in layers 2/3
of the rat visual cortex. Biol. Cybern. 82, 111121
(2000).
32. Averbeck, B. B. & Seo, M. The statistical
neuroanatomy of frontal networks in the macaque.
PLoS Comput. Biol. 4, e1000050 (2008).
33. Cherniak, C. Component placement optimization in
the brain. J. Neurosci. 14, 24182427 (1994).
34. Chklovskii, D. B., Schikorski, T. & Stevens, C. F. Wiring
optimization in cortical circuits. Neuron 34, 341347
(2002).
35. Klyachko, V. A. & Stevens, C. F. Connectivity
optimization and the positioning of cortical areas.
36. White, J. G., Southgate, E., Thomson, J. N. & Brenner, S.
The structure of the nervous system of the nematode
Caenorhabditis elegans. Philos. Trans. R. Soc. Lond. B
Biol. Sci. 314, 1340 (1986).
37. Silberberg, G., Grillner, S., LeBeau, F. E. N., Maex, R.
& Markram, H. Synaptic pathways in neural
microcircuits. Trends Neurosci. 28, 541551 (2005).
38. Humphries, M. D., Gurney, K. & Prescott, T. J. The
brainstem reticular formation is a small-world, not scalefree, network. Proc. Biol. Sci. 273, 503511 (2006).
39. Song, S., Sjstrm, P. J., Reigl, M., Nelson, S. &
Chklovskii, D. B. Highly nonrandom features of synaptic
connectivity in local cortical circuits. PLoS Biol. 3, e68
(2005).
Presented recordings from multiple cortical
neurons that revealed the small-world topology of
cellular functional networks.
40. Livet, J. et al. Transgenic strategies for combinatorial
expression of fluorescent proteins in the nervous
system. Nature 450, 5662 (2007).
41. Lichtman, J. W., Livet, J. & Sanes, J. R. A technicolour
approach to the connectome. Nature Rev. Neurosci. 9,
417422 (2008).
42. Felleman, D. J. & van Essen, D. C. Distributed
hierarchical processing in the primate cerebral cortex.
Cereb. Cortex 1, 147 (1991).
43. Scannell, J. W., Burns, G. A. P. C., Hilgetag, C. C.,
ONeil, M. A. & Young, M. P. The connectional
organization of the cortico-thalamic system of the cat.
Cereb. Cortex 9, 277299 (1999).
44. Sporns, O., Tononi, G. & Edelman, G. M. Theoretical
neuroanatomy and the connectivity of the cerebral
cortex. Cereb. Cortex 10, 127141 (2000).
One of the first papers to describe small-world
topological properties, and to investigate the
relationship between topology and complex
dynamics, in brain networks.
45. Hilgetag, C. C., Burns, G. A., ONeill, M. A., Scannell,
J. W. & Young, M. P. 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, 91110 (2000).
46. Sporns, O. & Zwi, J. The small world of the cerebral
cortex. Neuroinformatics 2, 145162 (2004).
47. Kaiser, M. & Hilgetag, C. C. Nonoptimal component
placement, but short processing paths, due to longdistance projections in neural systems. PLoS Comput.
Biol. 2, e95 (2006).
A comprehensive analysis of the relationship
between economical wiring and small-world
topology of brain networks, and its evolutionary
significance.
48. Sporns, O. & Ktter, R. Motifs in brain networks. PLoS
Biol. 2, 19101918 (2004).
49. Sporns, O., Honey, C. J. & Ktter, R. Identification and
classification of hubs in brain networks. PLoS ONE 2,
e1049 (2007).
50. Sporns, O., Tononi, G. & Ktter, R. The human
connectome: a structural description of the human
brain. PLoS Comp. Biol. 1, e42 (2005).
This review article argued for the fundamental
importance of structural connectivity in cognitive
neuroscience and proposed an effort to
systematically collect data on structural
connections in the human brain.
51. He, Y., Chen, Z. J. & Evans, A. C. Small-world

anatomical networks in the human brain revealed by
cortical thickness from MRI. Cereb. Cortex 17,
24072419 (2007).
This study was the first to derive a structural
network of the human brain on the basis of
correlations in cortical grey matter thickness
measured using MRI.
52. Wright, I. C. et al. Supra-regional brain systems and
the neuropathology of schizophrenia. Cereb. Cortex 9,
366378 (1999).
53. Chen, Z. J., He, Y., Rosa-Neto, P., Germann, J. &
Evans, A. C. Revealing modular architecture of human
brain structural networks by using cortical thickness
from MRI. Cereb. Cortex 18, 23742381 (2008).
54. Iturria-Medina, Y. et al. Characterizing brain anatomical
connections using diffusion weighted MRI and graph
theory. Neuroimage 36, 645660 (2007).
55. Iturria-Medina, Y., Sotero, R. C., Canales-Rodriguez,
E. J., Aleman-Gomez, Y. & Melie-Garcia, L. Studying
the human brain anatomical network via diffusionweighted MRI and graph theory. Neuroimage 40,
10641076 (2008).
56. Gong, G. et al. Mapping anatomical connectivity
patterns of human cerebral cortex using in vivo
diffusion tensor imaging tractography. Cereb. Cortex
20 Jun 2008 (doi:10.1093/cercor/bhn102).
57. Wedeen, V. J., Hagmann, P., Tseng, W. Y., Reese, T. G. &
Weisskoff, R. M. Mapping complex tissue architecture
with diffusion spectrum magnetic resonance imaging.
Magn. Reson. Med. 54, 13771386 (2005).
58. Hagmann, P. et al. Mapping human whole-brain
structural networks with diffusion MRI. PLoS ONE 2,
e597 (2007).
59. Hagmann, P. et al. Mapping the structural core of
human cerebral cortex. PLoS Biol. 6, e159 (2008).
This paper demonstrated the existence of modules,
hubs and a structural core in the human anatomical
network derived from DTI.
60. Parvizi, J., Van Hoesen, G. W., Buckwalter, J. &
Damasio, A. Neural connections of the posteromedial
cortex in the macaque. Proc. Natl Acad. Sci. USA 103,
15631568 (2006).
61. Cavanna, A. E. & Trimble, M. R. The precuneus: a
review of its functional anatomy and behavioural
correlates. Brain 129, 564583 (2006).
62. Alkire, M. T., Hudetz, A. G. & Tononi, G. Consciousness
and anesthesia. Science 322, 876880 (2008).
63. Stephan, K. E. et al. Computational analysis of
functional connectivity between areas of primate
cerebral cortex. Philos. Trans. R. Soc. Lond. B Biol. Sci.
355, 111126 (2000).
64. McIntosh, A. R. et al. Network analysis of cortical
visual pathways mapped with PET. J. Neurosci. 14,
655666 (1994).
65. Bullmore, E. T. et al. How good is good enough in path
analysis of fMRI data? Neuroimage 17, 573582
(2002).
66. Friston, K. J., Harrison, L. & Penny, W. Dynamic causal
modelling. Neuroimage 19, 12731302 (2003).
67. Brovelli, A. et al. Beta oscillations in a large-scale
sensorimotor cortical network: directional influences
revealed by Granger causality. Proc. Natl Acad. Sci.
USA 101, 98499854 (2004).
68. Salvador, R. et al. Neurophysiological architecture of
functional magnetic resonance images of human brain.
Cereb. Cortex 15, 13321342 (2005).
69. Eguluz, V. M., Chialvo, D. R., Cecchi, G. A., Baliki, M.
& Apkarian, A. V. Scale-free brain functional networks.
Phys. Rev. Lett. 94, 018102 (2005).
70. Achard, S., Salvador, R., Whitcher, B., Suckling, J. &
Bullmore, E. T. A resilient, low-frequency, small-world
human brain functional network with highly connected
association cortical hubs. J. Neurosci. 26, 6372
(2006).
This paper presented one of the first detailed
analyses of small-world brain functional networks
derived from human fMRI data.
71. Ferrarini, L. et al. Hierarchical functional modularity in
the resting-state human brain. Hum. Brain Mapp.
1 Oct 2008 (doi:10.1002/hbm.20663).
72. Meunier, D., Achard, S., Morcom, A. & Bullmore, E.
Age-related changes in modular organization of
human brain functional networks. Neuroimage 44,
715723 (2008).
73. Latora, V. & Marchiori, M. Efficient behaviour of
small-world networks. Phys. Rev. Lett. 87, 198701
(2001).
The first formulation of the economical small-world
concept and its key parameters: topological
efficiency and connection cost.
74. Latora, V. & Marchiori, M. Economic small-world

behavior in weighted networks. Eur. Phys. J. B 32,
249263 (2003).
75. Achard, S. & Bullmore, E. T. Efficiency and cost of
economical brain functional networks. PLoS Comput.
Biol. 3, e17 (2007).
76. Bullmore, E. T. et al. Wavelets and functional magnetic
resonance imaging of the human brain. Neuroimage
23, S234S249 (2004).
77. Fair, D. A. et al. Development of distinct cortical
networks through segregation and integration. Proc.
Natl Acad. Sci. USA 104, 1350713512 (2007).
78. Stam, C. J. & van Dijk, B. W. Synchronization
likelihood: an unbiased measure of generalized
synchronization in multivariate data sets. Physica D
163, 236251 (2002).
79. Stam, C. J. Functional connectivity patterns of human
magnetoencephalographic recordings: a small-world
network? Neurosci. Lett. 355, 2528 (2004).
80. Micheloyannis, S. et al. The influence of ageing on
complex brain networks: a graph theoretical analysis.
Hum. Brain Mapp. 30, 200208 (2009).
81. Bassett, D. S., Meyer-Lindenberg, A., Achard, S.,
Duke, T. & Bullmore, E. T. Adaptive reconfiguration of
fractal small-world human brain functional networks.
This study provides evidence for fractal or
scale-invariant small-world networks across
multiple frequency ranges and for their
reconfiguration during cognitive tasks.
82. Linkenkaer-Hansen, K., Nikouline, V. V., Palva, J. M. &
Ilmoniemi, R. J. Long-range temporal correlations and
scaling behavior in human brain oscillations.
J. Neurosci. 21, 13701377 (2001).
83. Maxim, V. et al. Fractional Gaussian noise, functional
MRI and Alzheimers disease. Neuroimage 25,
141158 (2005).
84. Achard, S., Bassett, D. S., Meyer-Lindenberg, A. &
Bullmore, E. T. Fractal connectivity of long memory
networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys.
77, 036104 (2008).
85. Schwarz, A., Gozzi, A. & Bifone, A. Community
structure and modularity in networks of correlated
brain activity. Magn. Reson. Imaging 26, 914920
(2008).
86. Yu, S., Huang, D., Singer, W. & Nikolic, D. A small world
of neuronal synchrony. Cereb. Cortex 18, 28912901
(2008).
This paper was one of the first to apply graph
theoretical techniques to map the topology of
functionally characterized cortical neuronal
circuits.
87. Schneidman, E., Still, S., Berry, M. J. & Bialek, W.
Network information and connected correlations.
Phys. Rev. Lett. 91, 238701 (2003).
88. Schneidman, E., Berry, M. J., Segev, R. & Bialek, W.
Weak pairwise correlations imply strongly correlated
network states in a neural population. Nature 440,
10071012 (2006).
89. Bettencourt, L. M., Stephens, G. J., Ham, M. I. &
Gross, G. W. Functional structure of cortical neuronal
networks grown in vitro. Phys. Rev. E Stat. Nonlin. Soft
Matter Phys. 75, 021915 (2007).
90. Barabsi, A. L. & Albert, R. Emergence of scaling in
random networks. Science 286, 509512 (1999).
This landmark paper was the first to describe the
scale-free organization of many complex networks
and proposed a simple growth rule for their
formation.
91. Van den Heuvel, M. P., Stam, C. J., Boersma, M. &
Hulshoff Pol, H. E. Small-world and scale-free
organization of voxel-based resting-state functional
connectivity in the human brain. Neuroimage 43,
528539 (2008).
92. Passingham, R. E., Stephan, K. E. & Ktter, R. The
anatomical basis of functional localization in the
cortex. Nature Rev. Neurosci. 3, 606616 (2002).
93. Alvarez, V. A. & Sabatini, B. L. Anatomical and
physiological plasticity of dendritic spines. Annu. Rev.
Neurosci. 30, 7997 (2007).
94. Grutzendler, J., Kasthuri, N. & Gan, W.-B. Long-term
dendritic spine stability in the adult cortex. Nature
420, 812816 (2002).
95. Marder, E. & Goaillard, J.-M. Variability, compensation
and homeostasis in neuron and network function.
Nature Rev. Neurosci. 7, 563574 (2006).
96. Harris, K. D., Csicsvari, J., Hirase, H., Dragoi, G. &
Buzski, G. Organization of cell assemblies in the
hippocampus. Nature 424, 552556 (2003).
97. Sasaki, T., Matsuki, N. & Ikegaya, Y. Metastability of
active CA3 networks. J. Neurosci. 27, 517528 (2007).
VOlume 10 | maRcH 2009 | 197

REVIEWS
98. Valencia, M., Martinerie, J., Dupont, S. & Chavez, M.
Dynamic small-world behaviour in functional brain
networks unveiled by an event-related networks
approach. Phys. Rev. E Stat. Nonlin. Soft Matter Phys.
77, 050905 (2008).
99. Raichle, M. E. et al. A default mode of brain function.
Using quantitative metabolic and haemodynamic
measures, this paper first proposed the existence
of an organized pattern of resting or default-mode
brain activity.
100. Gusnard, D. A. & Raichle, M. E. Searching for a
baseline: functional imaging and the resting human
brain. Nature Rev. Neurosci. 2, 685694 (2001).
101. Fox, M. D. & Raichle, M. E. Spontaneous fluctuations
in brain activity observed with functional magnetic
resonance imaging. Nature Rev. Neurosci. 8, 700711
(2007).
102. Greicius, M. D., Krasnow, B., Reiss, A. L. & Menon, V.
Functional connectivity in the resting brain: a network
analysis of the default mode hypothesis. Proc. Natl
Acad. Sci. USA 100, 253258 (2003).
103. Fox, M. D. et al. The human brain is intrinsically
organized into dynamic, anticorrelated functional
networks. Proc. Natl Acad. Sci. USA 102, 96739678
(2005).
104. Koch, M. A., Norris, D. G. & Hund-Georgiadis, M. An
investigation of functional and anatomical connectivity
using magnetic resonance imaging. Neuroimage 16,
241250 (2002).
105. Greicius, M., Supekar, K., Menon, V. & Dougherty,
R. F. Resting-state functional connectivity reflects
structural connectivity in the default mode network.
Cereb. Cortex 19, 7278 (2008).
106. Honey, C. J. et al. Predicting human resting-state
functional connectivity from structural connectivity.
Proc. Natl Acad. Sci. USA (in the press).
107. Park, C. H., Kim, S. Y., Kim, Y.-H. & Kim, K. Comparison
of the small-world topology between anatomical and
functional connectivity in the human brain. Physica A
387, 59585962 (2008).
108. Galn, R. F. On how network architecture determines
the dominant patterns of spontaneous neural activity.
PLoS ONE 3, e2148 (2008).
109. Honey, C. J., Ktter, R., Breakspear, M. & Sporns, O.
Network structure of cerebral cortex shapes functional
connectivity on multiple time scales. Proc. Natl Acad.
Sci. USA 104, 1024010245 (2007).
This paper used a large-scale computational model
to relate topological features of structural and
functional brain networks at multiple timescales.
110. Ghosh, A., Rho, Y., McIntosh, A. R., Ktter, R. & Jirsa,
V. K. Cortical network dynamics with time delays
reveals functional connectivity in the resting brain.
Cogn. Neurodyn. 2, 115120 (2008).
111. Zhou, C., Zemanova, L., Zamora, G., Hilgetag, C. C. &
Kurths, J. Hierarchical organization unveiled by
functional connectivity in complex brain networks.
Phys. Rev. Lett. 97, 238103 (2006).
112. Mller-Linow, M., Hilgetag, C. C. & Htt, M.-T.
Organization of excitable dynamics in hierarchical
biological networks. PLoS Comput. Biol. 4, e1000190
(2008).
113. Kaiser, M., Grner, M. & Hilgetag, C. C. Criticality of
spreading dynamics in hierarchical cluster networks
without inhibition. New J. Phys. 9, 110 (2007).
114. Percha, B., Dzakpasu, R., Zochowski, M. & Parent, J.
Transition from local to global phase synchrony in
small world neural network and its possible
implications for epilepsy. Phys. Rev. E Stat. Nonlin.
Soft Matter Phys. 72, 031909 (2005).
115. Siri, B., Quoy, M., Delord, B., Cessac, B. & Berry, H.
Effects of Hebbian learning on the dynamics and
structure of random networks with inhibitory
and excitatory neurons. J. Physiol. (Paris) 101,
136148 (2007).
116. Catani, M. & fftyche, D. H. The rises and falls of
disconnection syndromes. Brain 128, 22242239
(2005).
117. Supekar, K., Menon, V., Rubin, D., Musen, M. &
Greicius, M. D. Network analysis of intrinsic functional
brain connectivity in Alzheimers disease. PLoS
Comput. Biol. 4, e1000100 (2008).
118. Stam, C. J., Jones, B. E., Nolte, G., Breakspear, M. &
Scheltens, P. Small-world networks and functional
connectivity in Alzheimers disease. Cereb. Cortex 17,

9299 (2007).
This paper was one of the first to use graph theory
to demonstrate disease-related differences in brain
functional network topology.
119. He, Y., Chen, Z. & Evans, A. C. Structural insights into
aberrant topological patterns of large-scale cortical
networks in Alzheimers disease. J. Neurosci. 28,
81488159 (2008).
120. Stam, C. J. et al. Graph theoretical analysis of
magnetoencephalographic functional connectivity in
Alzheimers disease. Brain 24 Oct 2008 (doi:10.1093/
brain/awn262).
121. Liu, Y. et al. Disrupted small-world networks in
schizophrenia. Brain 131, 945961 (2008).
122. Micheloyannis, S. et al. Small-world networks and
disturbed functional connectivity in schizophrenia.
Schizophr. Res. 87, 6066 (2006).
123. Rubinov, M. et al. Small-world properties of nonlinear
brain activity in schizophrenia. Hum. Brain Mapp. 10
Dec 2007 (doi:10.1002/hbm.20517).
124. Bassett, D. S. et al. Hierarchical organization of
human cortical networks in health and schizophrenia.
J. Neurosci. 28, 92399248 (2008).
125. Ponten, S. C., Bartolomei, F. & Stam, C. J. Small-world
networks and epilepsy: graph theoretical analysis of
intracerebrally recorded mesial temporal lobe
seizures. Clin. Neurophysiol. 118, 918927 (2007).
126. Kramer, M. A., Kolaczyk, E. D. & Kirsch, H. E.
Emergent network topology at seizure onset in
humans. Epilepsy Res. 79, 173186 (2008).
127. Schindler, K. A., Bialonski, S., Horstmann, M. T.,
Elger, C. E. & Lehnertz, K. Evolving functional network
properties and synchronizability during human
epileptic seizures. Chaos 18, 033119 (2008).
128. Wang, L. et al. Altered small-world brain functional
networks in children with attention-deficit/
hyperactivity disorder. Hum. Brain Mapp. 24 Jan
2008 (doi:10.1002/hbm.20530).
129. De Vico Fallani, F. et al. Cortical functional connectivity
networks in normal and spinal cord injured patients:
evaluation by graph analysis. Hum. Brain Mapp. 28,
13341346 (2007).
130. Smit, D. J., Stam, C. J., Posthuma, D., Boomsma, D. I.
& de Geus, E. J. Heritability of small-world networks
in the brain: a graph theoretical analysis of restingstate EEG functional connectivity. Hum. Brain Mapp.
29, 13681378 (2008).
131. Schmitt, J. E. et al. Identification of genetically
mediated cortical networks: a multivariate study of
pediatric twins and siblings. Cereb. Cortex 18,
17371747 (2008).
132. Sporns, O. Small-world connectivity, motif composition
and complexity of fractal neuronal connections.
Biosystems 85, 5564 (2006).
133. Kaiser, M., Robert, M., Andras, P. & Young, M. P.
Simulation of robustness against lesions of cortical
networks. Eur. J. Neurosci. 25, 31853192 (2007).
134. Honey, C. J. & Sporns, O. Dynamical consequences of
lesions in cortical networks. Hum. Brain Mapp. 29,
802809 (2008).
135. He, B. J., Shulman, G. L., Snyder, A. Z. & Corbetta, M.
The role of impaired neuronal communication in
neurological disorders. Curr. Opin. Neurol. 20,
655660 (2007).
136. He, B. J. et al. Breakdown of functional connectivity in
frontoparietal networks underlies behavioral deficits
in spatial neglect. Neuron 53, 905918 (2007).
137. Dyhrfjeld-Johnsen, J. et al. Topological determinants
of epileptogenesis in large-scale structural and
functional models of the dentate gyrus derived from
experimental data. J. Neurophysiol. 97, 15661587
(2007).
This paper used biologically realistic computational
modelling to study the effects of epileptogenic
cellular changes on the topology and dynamics of
functional networks in the rat hippocampus.
138. Srinivas, K. V., Jain, R., Saurav, S. & Sikdar, S. K.
Small-world network topology of hippocampal
neuronal network is lost, in an in vivo glutamate injury
model of epilepsy. Eur. J. Neurosci. 25, 32763286
(2007).
139. Netoff, T. I., Clewley, R., Arno, S., Keck, T. & White,
J. A. Epilepsy in small-world networks. J. Neurosci.
24, 80758083 (2004).
198 | maRcH 2009 | VOlume 10
140. Honey, G. D. et al. Dopaminergic drug effects on

physiological connectivity in a human cortico-striatothalamic system. Brain 126, 17671781 (2003).
141. Schwarz, A. J., Gozzi, A., Reese, T., Heidbreder, C. A. &
Bifone, A. Pharmacological modulation of functional
connectivity: the correlation structure underlying the
phMRI response to d-amphetamine modified by
selective dopamine D3 receptor antagonist SB277011A.
Magn. Reson. Imaging 25, 277811277820 (2007).
142. Stoffers, D., Bosboom, J. L., Wolters, E. Ch., Stam, C. J.
& Berendse, H. W. Dopaminergc modulation of corticocortical functional connectivity in Parkinsons disease:
an MEG study. Exp. Neurol. 213, 191195 (2008).
143. Bressler, S. & Kelso, J. A. S. Cortical coordination
dynamics and cognition. Trends Cogn. Sci. 5, 2636
(2001).
144. Barahona, M. & Pecora, L. M. Synchronization in smallworld systems. Phys. Rev. Lett. 89, 054101 (2002).
145. Shin, C. W. & Kim, S. Self-organized criticality and
scale-free properties in emergent functional neural
networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys.
74, 045101 (2006).
146. Kitzbichler, M., Smith, M., Sorensen, C. & Bullmore, E.
Broadband criticality of human brain network
synchronization. PLoS Comput. Biol. (in the press).
147. Wang, J. et al. Parcellation-dependent small-world brain
functional networks: a resting-state fMRI study. Hum.
Brain Mapp. 22 Jul 2008 (doi:10.1002/hbm.20623).
148. Roebroeck, A., Formisano, E. & Goebel, R. Mapping
directed influence over the brain using Granger
causality and fMRI. Neuroimage 25, 230242 (2005).
149. Bressler, S. L., Tang, W., Sylvester, C., Shulman, G. &
Corbetta, M. Top-down control of human visual cortex
by frontal and parietal cortex in anticipatory visual
spatial attention. J. Neurosci. 28, 1005610061
(2008).
150. Milo, R. et al. Network motifs: simple building blocks
of complex networks. Science 298, 824827 (2002).
151. Freeman, L. C. A set of measures of centrality based
on betweenness. Sociometry 40, 3541 (1977).
152. Erds, P. & Rnyi, A. On the evolution of random
graphs. Publ. Math. Inst. Hung. Acad. Sci. 5, 1761
(1960).
153. Milgram, S. The small world problem. Psychol. Today
1, 6167 (1967).
154. Humphries, M. D. & Gurney, K. Network small-worldness: a quantitative method for determining canonical
network equivalence. PLoS ONE 3, e0002051 (2008).
155. Harary, F. Graph Theory (Perseus, Reading,
Massachusetts, 1969).
156. Euler, L. Solutio problematis ad geometriam situs
pertinentis. Commentarii Academiae Scientiarum
Imperialis Petropolitanae 8, 128140 (1736).
157. Wasserman, S. & Faust, K. Social Network Analysis:
Methods and Applications (Cambridge Univ. Press,
1994).
158. Sporns, O. in Diffusion MRI: from Quantitative
Measurement to InVivo Neuroanatomy (eds JohansenBerg, H. & Behrens, T.) 309332 (Academic, London,
2009).
Acknowledgements
E.B. was supported by a Human Brain Project grant from the

National Institute of Mental Health and the National Institute
of Biomedical Imaging & Bioengineering. The Behavioural &
Clinical Neurosciences Institute in the University of Cambridge
is supported by the Wellcome Trust and the Medical Research
Council (UK). O.S. was supported by the JS McDonnell
Foundation.
Competing interests statement
The authors declare competing financial interests: see web

version for details.
FURTHER INFORMATION
Ed Bullmores homepage: http://www.neuroscience.cam.
ac.uk/directory/profile.php?edbullmore
Olaf Sporns homepage: http://www.indiana.edu/~cortex
Brain Connectivity Toolbox: http://www.brainconnectivity-toolbox.net/
BrainwaveR Toolbox: http://www.nitrc.org/projects/
brainwaver/
All liNks Are Active iN tHe oNliNe PDf

Complex Brain Networks Graph Theoretical Analysis of Structural and Functional Systems

Uploaded by

Copyright:

Available Formats

Complex Brain Networks Graph Theoretical Analysis of Structural and Functional Systems

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Complex Brain Networks Graph Theoretical Analysis of Structural and Functional Systems

Uploaded by

Copyright:

Available Formats

REVIEWS

Complex brain networks: graph

Abstract | Recent developments in the quantitative analysis of complex networks, based

We have known since the nineteenth century that the

has led to a fundamental insight: substantively different

186 | maRcH 2009 | VOlume 10

Time series data

Structural brain network

Functional brain network

Graph theoretical analysis

naTuRe ReVIeWS | NeuroscieNce

VOlume 10 | maRcH 2009 | 187

Box 2 | Network measures

Node degree, degree distribution and assortativity

Nature Reviews | Neuroscience

preparations and human studies that have contributed to

Structural brain networks

distances between elements in brain networks are often

naTuRe ReVIeWS | NeuroscieNce

VOlume 10 | maRcH 2009 | 189

Diffusion spectrum imaging

that are spatially close and functionally related. These

Functional brain networks

190 | maRcH 2009 | VOlume 10

Mapping functional networks using fMRI. The first

These concepts were translated to the analysis of

naTuRe ReVIeWS | NeuroscieNce

VOlume 10 | maRcH 2009 | 191

Conservation of functional network properties. There

anatomical resolution on the extent to which power

Structurefunction relations in brain networks

192 | maRcH 2009 | VOlume 10

Coupling matrix and neural dynamics

Functional brain network

Figure 1 | computational modelling of structural and functional brain networks.

of participants63,106,107 suggest that structural connections

between neural units, which then results in spatially

Clinical and translational aspects

naTuRe ReVIeWS | NeuroscieNce

VOlume 10 | maRcH 2009 | 193

Figure 2 | cellular and whole-brain networks demonstrate consistent topological

brain networks derived from fmRI, eeG or structural

194 | maRcH 2009 | VOlume 10

b People with schizophrenia

Figure 3 | Disease-related disorganization of brain anatomical networks derived

worth investigating as intermediate phenotypes, or

naTuRe ReVIeWS | NeuroscieNce

VOlume 10 | maRcH 2009 | 195

Conclusions and prospects

Cajal, S. R. Histology of the Nervous System of Man

in terms of efficiency and could arguably account for

A scholarly review of the early literature on the

196 | maRcH 2009 | VOlume 10

This seminal paper on small-world networks

51. He, Y., Chen, Z. J. & Evans, A. C. Small-world

naTuRe ReVIeWS | NeuroscieNce

74. Latora, V. & Marchiori, M. Economic small-world

VOlume 10 | maRcH 2009 | 197

connectivity in Alzheimers disease. Cereb. Cortex 17,

198 | maRcH 2009 | VOlume 10