Skip to main content

Umbrella menu

  • SfN.org
  • eNeuro
  • The Journal of Neuroscience
  • Neuronline
  • BrainFacts.org

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
  • ALERTS
  • FOR AUTHORS
    • Preparing a Manuscript
    • Submission Guidelines
    • Fees
    • Journal Club
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
  • SfN.org
  • eNeuro
  • The Journal of Neuroscience
  • Neuronline
  • BrainFacts.org

User menu

  • Log in
  • Subscribe
  • My alerts

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • Subscribe
  • My alerts
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
  • ALERTS
  • FOR AUTHORS
    • Preparing a Manuscript
    • Submission Guidelines
    • Fees
    • Journal Club
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Behavioral/Systems/Cognitive

A Resilient, Low-Frequency, Small-World Human Brain Functional Network with Highly Connected Association Cortical Hubs

Sophie Achard, Raymond Salvador, Brandon Whitcher, John Suckling and Ed Bullmore
Journal of Neuroscience 4 January 2006, 26 (1) 63-72; DOI: https://doi.org/10.1523/JNEUROSCI.3874-05.2006
Sophie Achard
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Raymond Salvador
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brandon Whitcher
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John Suckling
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ed Bullmore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Article Figures & Data

Figures

  • Tables
  • Figure 1.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 1.

    Schematic of wavelet correlation analysis, thresholding, and functional network visualization. Top, fMRI time series recorded from each of 90 regions in each subject are decomposed using the MODWT, and the inter-regional correlation is estimated at each scale of the MODWT for each pair of regions in each subject; individual wavelet correlation matrices are then averaged over subjects at each scale to produce a set of six group mean wavelet correlation matrices. Middle, The wavelet correlation matrices are thresholded to generate binary matrices, each element of which is either black (if there is no significant connection between regions) or white (if there is). The stringency of the probabilistic thresholding operation is determined by the value of the correlation threshold R, as illustrated by applying three different thresholds (R = 0.3, 0.4, 0.5) to the scale 4 wavelet correlation matrix. Bottom, Thresholded matrices are visualized in anatomical space by locating each region according to its y and z centroid coordinates in Talairach space and drawing an edge between regions that are significantly connected.

  • Figure 2.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 2.

    Small-world properties of brain networks as a function of correlation threshold. a, As the correlation threshold R is increased, mean degree k monotonically decreases (the networks become more sparsely connected) at all scales of the wavelet transform: black lines, scale 1; red lines, scale 2; green lines, scale 3; dark blue lines, scale 4; light blue lines, scale 5; purple lines, scale 6. When the mean degree is less than the log of the number of regions [i.e., knet < log(n)], small-world properties are not estimable. b, The largest cluster size also tends to decrease as the correlation threshold is increased (i.e., the networks become progressively more fragmented as connections are eliminated at higher thresholds). c, The ratio γ = Cnet/Cran tends to increase as the correlation threshold is increased (i.e., compared with closely matched random networks, the brain networks demonstrate progressively greater clustering at higher thresholds). d, The ratio λ = Lnet/Lran is only modestly increased as a function of correlation threshold (i.e., compared with closely matched random networks, the brain networks demonstrate approximately equivalent path lengths over all thresholds). e, The ratio σ = γ/λ, a scalar summary of small-worldness, therefore tends to increase as a function of increasing the correlation threshold. Note that the evidence for small-world properties is clearest at high thresholds in scales 4 and 5 (collectively corresponding to the frequency interval 0.01–0.06Hz), but there is some evidence for small-world properties (σ > 1) over a range of thresholds in all scales.

  • Figure 3.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 3.

    Anatomical map of a small-world human brain functional network created by thresholding the scale 4 wavelet correlation matrix representing functional connectivity in the frequency interval 0.03–0.06 Hz. a, Four hundred five undirected edges, ∼10% of the 4005 possible inter-regional connections, are shown in a sagittal view of the right side of the brain. Nodes are located according to the y and z coordinates of the regional centroids in Talairach space. Edges representing connections between nodes separated by a Euclidean distance <7.5 cm are red; edges representing connections between nodes separated by Euclidean distance >7.5 cm are blue. b, Short-distance connections, predominantly in the posterior cortex, are shown separately in red. c, Long-distance connections (e.g., between the frontal cortex and regions of the parietal and temporal association cortex) are shown separately in blue.

  • Figure 4.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 4.

    Topological map of a small-world human brain functional network created by thresholding the scale 4 wavelet correlation matrix representing functional connectivity in the frequency interval 0.03–0.06 Hz. The regions have been located by multidimensional scaling of the (binary) thresholded matrix so that the distance between them in the space of this plot approximates the path length between them; network hubs are clustered centrally, and less well connected regions are located peripherally. Long-distance connections subtending a Euclidean distance between regional centroids >7.5 cm are drawn as black lines, and short-distance connections are drawn as gray lines; regional labels and color codes are as in Figure 6. See Table 2 for the list of abbreviations.

  • Figure 5.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 5.

    Degree distribution of a small-world brain functional network. a, Histogram of regional degree ki distribution. b, Plot of the log of the cumulative probability of degree, log(P(ki)), versus log of degree, log(ki). The plus sign indicates observed data, the solid line is the best-fitting exponentially truncated power law, the dotted line is an exponential, and the dashed line is a power law.

  • Figure 6.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 6.

    Relationship between local clustering and mean physical distance of connections to brain regions. The scatter plot of Euclidean distance Di (y-axis) versus clustering coefficient Ci (x-axis) is shown. The regions are labeled with the abbreviations in Table 2 and color-coded as follows: black, association cortex; red, paralimbic/limbic cortex; green, primary sensory or motor cortex. The fitted lines are shown for regression of distance on clustering for neocortical (black) and limbic/paralimbic (red) regions; the dotted lines indicate the network mean values of distance and clustering coefficient. Distance and clustering were negatively correlated over neocortical (but not limbic) regions; the unimodal association cortex tended to have high clustering and short mean connection distances, whereas the heteromodal association cortex tended to have the opposite pattern of low clustering and long mean connection distances.

  • Figure 7.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 7.

    Resilience of the human brain functional network (rightcolumn) compared with random (left column) and scale-free (middle column) networks. Top row, Size of the largest connected cluster in the network (scaled to maximum; y-axis) versus the proportion of total nodes eliminated (x-axis) by random error (dashed line) or targeted attack (solid line). Bottom row, Global mean path length (Lnet; y-axis) versus the proportion of total nodes eliminated (x-axis) by random error (dashed line) or targeted attack (solid line). The size of the largest connected cluster in the brain functional network is more resilient to targeted attack and about equally resilient to random error compared with the scale-free network.

Tables

  • Figures
    • View popup
    Table 1.

    Wavelet scale dependency of functional connectivity and small-world parameters for an entire human brain network

    Scale Hz rRLnetCnet λ γ σ
    1 0.23-0.45 0.12 0.13 2.9 0.534 1.28 1.81 1.42
    2 0.11-0.23 0.21 0.2 2.6 0.566 1.12 2.14 1.92
    3 0.06-0.11 0.39 0.39 2.69 0.555 1.16 2.22 1.91
    4 0.03-0.06 0.45 0.44 2.49 0.525 1.08 2.38 2.19
    5 0.01-0.03 0.44 0.35 2.4 0.554 1.04 2.39 2.30
    6 0.007-0.01 0.41 0.17 2.65 0.515 1.15 2.15 1.88
    • Scales 1-6 of the MODWT denote progressively lower-frequency intervals (Hz). r is the mean inter-regional correlation, and R is the correlation threshold. Lnet and Cnet are the mean path length and clustering coefficient, respectively, of the thresholded network. The λ and γ are ratios of brain network path length and clustering coefficient, respectively, to comparable random network metrics. The equation σ = γ/λ is a scalar measure of “small-worldness.”

    • View popup
    Table 2.

    Regions of an entire human brain functional network ranked in order of increasing path length

    Abbreviation Region Class LikiCiDi ΔLnet
    PCUN Precuneus Association 1.95 22.5 0.31 57.6 2.59
    MTG Middle temporal gyrus Association 2.03 16.5 0.46 76.6 2.44
    MFG Middle frontal gyrus Association 2.10 18.0 0.35 72.6 3.55
    LING Lingual gyrus Association 2.11 21 0.54 55.8 2.07
    MOG Middle occipital gyrus Association 2.14 20.5 0.55 57.0 2.33
    SPG Superior parietal gyrus Association 2.14 15.5 0.43 61.7 2.23
    PoCG Postcentral gyrus Primary 2.15 20 0.47 62.0 1.34
    ITG Inferior temporal gyrus Association 2.18 14 0.52 95.9 1.58
    PreCG Precentral gyrus Primary 2.18 20 0.44 62.4 2.49
    DCG Dorsal cingulate gyrus Paralimbic 2.21 9.5 0.46 42.2 2.08
    SFGdor Superior frontal gyrus (dorsal) Association 2.23 13 0.48 72.2 1.32
    CAL Calcarine cortex Primary 2.23 14.5 0.73 37 0.93
    STG Superior temporal gyrus Association 2.24 17 0.43 81.9 2.57
    FFG Fusiform gyrus Association 2.30 16.5 0.63 63.7 3.77
    SMA Supplementary motor area Association 2.39 11 0.62 49.4 1.13
    CUN Cuneus Association 2.41 11 0.86 34.3 0.49
    IPL Inferior parietal lobule Association 2.42 9.5 0.56 72.4 0.54
    ORBsup Orbitofrontal cortex (superior) Paralimbic 2.42 9.5 0.37 67.2 3.82
    IFGtriang Inferior frontal gyrus (triangular) Association 2.44 9.5 0.49 71.6 0.61
    SOG Superior occipital gyrus Association 2.49 14 0.74 45.0 0.28
    THA Thalamus Subcortical 2.51 8 0.51 42.8 0.83
    ORBinf Orbitofrontal cortex (inferior) Paralimbic 2.61 6 0.51 57.1 1.14
    ORBmid Orbitofrontal cortex (middle) Paralimbic 2.66 6 0.57 59.5 -0.15
    SMG Supramarginal gyrus Association 2.70 5.5 0.70 77 -0.27
    PCG Posterior cingulate gyrus Paralimbic 2.73 4.5 0.68 30.4 -0.35
    IOG Inferior occipital gyrus Association 2.75 6.5 0.92 58.7 -0.40
    PCL Paracentral lobule Association 2.77 6 0.93 26.6 -0.43
    ROL Rolandic operculum Association 2.77 7.5 0.74 62.1 -0.45
    IFGoperc Inferior frontal gyrus (opercular) Association 2.78 5.5 0.62 65.0 -0.46
    SFGmed Superior frontal gyrus (medial) Association 2.8 7 0.50 49.6 0.49
    INS Insula Paralimbic 2.94 5.5 0.63 68.6 -0.64
    ANG Angular gyrus Association 3.03 5 0.0.20 79.1 -0.89
    HES Heschl's gyrus Primary 3.07 3.5 0.92 60.4 -1.21
    TPOsup Temporal pole (superior) Paralimbic 3.14 5 0.42 70.7 -0.38
    TPOmid Temporal pole (middle) Paralimbic 3.2 3 0.08 80.3 -0.46
    REC Rectus gyrus Paralimbic 3.2 3 0.5 36 -1.53
    HIP Hippocampus Limbic 3.3 3.5 0.15 34.5 -0.25
    CAU Caudate Subcortical 3.35 3 1 33.6 -1.94
    PHG Parahippocampal gyrus Paralimbic 3.38 2.5 0 28 -1.08
    ORBsupmed Orbitofrontal cortex (superior-medial) Paralimbic 3.46 2.5 0 26.2 -2.09
    • Major “hubs” of the network are listed first. Li, Ci, and ki are the characteristic path length, clustering coefficient, and degree, respectively, of the ith region. Di is the mean Euclidean distance of connections to that region. ΔLnet is the percentage of change in global mean path length when the network is attacked by eliminating the ith region and its connections. The regions are classified as the primary, association, paralimbic, or limbic cortex as described by Mesulam (2000). Network hubs, defined as those regions with characteristic path length less than global network mean path length (2.49), are listed in bold; clustering coefficients are italicized for those hubs in which they are less than the global mean clustering coefficient (0.525).

Back to top

In this issue

The Journal of Neuroscience: 26 (1)
Journal of Neuroscience
Vol. 26, Issue 1
4 Jan 2006
  • Table of Contents
  • About the Cover
  • Index by author
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
A Resilient, Low-Frequency, Small-World Human Brain Functional Network with Highly Connected Association Cortical Hubs
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Citation Tools
A Resilient, Low-Frequency, Small-World Human Brain Functional Network with Highly Connected Association Cortical Hubs
Sophie Achard, Raymond Salvador, Brandon Whitcher, John Suckling, Ed Bullmore
Journal of Neuroscience 4 January 2006, 26 (1) 63-72; DOI: 10.1523/JNEUROSCI.3874-05.2006

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
A Resilient, Low-Frequency, Small-World Human Brain Functional Network with Highly Connected Association Cortical Hubs
Sophie Achard, Raymond Salvador, Brandon Whitcher, John Suckling, Ed Bullmore
Journal of Neuroscience 4 January 2006, 26 (1) 63-72; DOI: 10.1523/JNEUROSCI.3874-05.2006
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Appendix
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

  • Influence of Reward on Corticospinal Excitability during Movement Preparation
  • Identification and Characterization of a Sleep-Active Cell Group in the Rostral Medullary Brainstem
  • Gravin Orchestrates Protein Kinase A and β2-Adrenergic Receptor Signaling Critical for Synaptic Plasticity and Memory
Show more Behavioral/Systems/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(JNeurosci logo)
(SfN logo)

Copyright © 2021 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.