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

**|**HHS Author Manuscripts**|**PMC3854964

Formats

Article sections

Authors

Related links

Med Image Comput Comput Assist Interv. Author manuscript; available in PMC 2013 December 6.

Published in final edited form as:

Med Image Comput Comput Assist Interv. 2012; 15(0 2): 196–203.

PMCID: PMC3854964

NIHMSID: NIHMS522835

Johnson J. GadElkarim,^{1,}^{2} Dan Schonfeld,^{1} Olusola Ajilore,^{2} Liang Zhan,^{3} Aifeng Zhang,^{2} Jamie D. Feusner,^{4} Paul M. Thompson,^{3} Tony J. Simon,^{5} Anand Kumar,^{2} and Alex D. Leow^{2,}^{6}

jgadel2/at/uic.edu, alexfeuillet/at/gmail.com

See other articles in PMC that cite the published article.

We propose a framework for quantifying node-level community structures between groups using anatomical brain networks derived from DTI-tractography. To construct communities, we computed hierarchical binary trees by maximizing two metrics: the well-known modularity metric (Q), and a novel metric that measures the difference between inter-community and intra-community path lengths. Changes in community structures on the nodal level were assessed between generated trees and a statistical framework was developed to detect local differences between two groups of community structures. We applied this framework to a sample of 42 subjects with major depression and 47 healthy controls. Results showed that several nodes (including the bilateral precuneus, which have been linked to self-awareness) within the default mode network exhibited significant differences between groups. These findings are consistent with those reported in previous literature, suggesting a higher degree of ruminative self-reflections in depression.

Modern imaging techniques have allowed us to construct anatomical (structural) and functional networks of the human brain and a new term - the human “connectome” [1] has been recently coined to describe their mathematical properties. A brain connectivity network consists of nodes which represent gray matter regions and edges connecting nodes which represent white matter fibers. Many tools and metrics have been adapted from the field of graph-theory to help characterize and analyse the human connectome [1-2]. It has been shown that brain networks are sparse (networks having the number of nodes in the order of the number of edges) [3], thus allowing the investigation of their community structure properties in which nodes may be grouped into different communities [4, 5] (in our context, “*community structure*” is a property for sparse networks, which allows us to extract the topological organization of a network and to partition it into a set of non-overlapping *communities*, sometimes also referred to as *modules* or *clusters*).

Studying the structural organization of the brain communities has lately shed light on alterations related to aging [6, 7] as well as neuropsychiatric disorders such as schizophrenia [8-12]. In this paper, we proposed a novel metric, with which we extracted community structures in the form of top-down binary trees (also known as dendrograms). This is followed by the computation of a node “consistency metric”, designed to quantify differences between trees at each node. To our best knowledge, we are the first group to propose a recipe for quantifying node-level community differences between two groups of brain networks. Lastly, statistical analyses were conducted to detect significant community alterations on the nodal level for group studies. We tested the proposed method and compared the binary trees generated from the structural brain networks of depressed versus normal control subjects. Significant changes were found on the nodal level in the community structure in depression relative to controls.

We scanned 47 healthy subjects (20 male/27 female; age: 59.74±14.8)) and 42 subjects with major depression as defined by DSM-IV (18 male/24 female; age: 57.64±13.4). A Philips 3.0 T Achieva scanner supplied with 8-channel SENSE head-coil located at the University of Illinois Medical Center Advanced Imaging Center was used to acquire the brain MRI data. High resolution T1-weighted images were acquired with MPRAGE sequence (FOV=240 mm; TR/TE=8.4/3.9 ms; flip angle 8^{°}; voxel size 1.1×1.1×1.1mm). Diffusion weighted (DW) images were acquired using SS-SE-EPI sequences (FOV = 240 mm; resolution 0.83×0.83 mm; TR/TE=6994/71 ms; flip angle = 90^{°}, 32 gradient directions, b= 700 s/mm^{2} and one minimally DW scan: b0 image). Parallel imaging was also used with a SENSE factor of 2.5 to reduce scan time to ~ 4 min.

Structural brain networks were generated using a pipeline which integrates multiple image analysis techniques. First, DW images were eddy current corrected using the automatic image registration (AIR) tool embedded in DtiStudio software (http://www.mristudio.org) by registering all DW images to their corresponding b0 images with 12-parameter a ne transformations. This was followed by the computation of diffusion tensors and deterministic tractography using the DtiStudio program. T1-weighted images were used to generate label maps using the Freesurfer software (http://surfer.nmr.mgh.harvard.edu). Brain networks formed by the 68 cortical regions were generated using an in-house program in Matlab by counting the number of fibers connecting each pair of nodes.

Measuring community structure has been a scientifically important task in network science. To this end, we note that nodes belonging to the same community should exhibit stronger connections than those in different communities. A well-known metric for such task is the Q modularity proposed in [13] and is mathematically defined as:

(1)

where Q is a function of a graph G, m is the total number of edges, A* _{ij}* = 1 if an edge links nodes i and j and 0 otherwise, δ(i, j) = 1 if nodes i and j are in the same community and 0 otherwise, and k

(2)

where N* _{i}* denotes the number of nodes in a community C

(3)

Here, notice that can be considered the weighted characteristic path length (a measure of global network integration) defined on the sub-network formed by C* _{i}*. Thus, maximizing

There exist two main types of hierarchical clustering techniques: the agglomerative and the divisive. The agglomerative method starts bottom-up. At the first level, every node is a community by itself. In each subsequent step (level), two communities that are considered closest (with respect to a similarity measure) are merged into one, until all nodes are in one community. The divisive method follows the opposite or a top-down approach. It starts with all nodes belonging to one single community, and at each step each community is further split into several sub-communities (unlike the agglomerative method, the divisive clustering usually uses a dis-similarity measure). It should be noted that one may choose to stop the clustering process at a pre-specified level (stopping criteria can be a given number of communities or the optimization of a quality function [14]). Here, to maximize Q or *Ψ ^{PL}*, we used a top-down hierarchical clustering approach, as it has been suggested that agglomerative methods may fail to find correct communities in networks where the community structures are known a priori [13]. At the first level, brain regions were randomly assigned to one of two communities, and optimal assignment was determined by maximizing Q or

In order to assess local differences between communities in different groups, we first construct the mean binary tree for the control group (by maximizing Q or *Ψ ^{PL}*, with respect to the mean connectivity matrix for this group). All individual subjects’ trees will then be compared to the mean control tree. To this end, a local metric is needed to quantify how two trees differ. One candidate metric is via the consistency metric assessment method introduced in [16], which yields a consistency vector

(4)

where N* _{pq}* denotes the number of common nodes between the two communities C

(5)

where N* _{c}* denotes the number of common nodes between the two communities C

To statistically test group differences in community structures at the nodal level, all individual subjects’ trees (both depressed and control) are compared to the mean control reference tree, thus yielding 47 *V* or *V ^{m}* vectors in the healthy group, and 42 in the depressed group. Group differences on local community structures can now be assessed using node-wise 2-sample Mann-Whitney-Wilcoxon (MWW) tests for

We first generated and compared the mean community structures of the normal controls in this study, using both Q and *Ψ ^{PL}* (figure 1). Note that in the first level of the binary tree (partitioning the brain into two modules), the constructed mean tree using the Q function exhibits an exact left/right separation while the tree generated using

To detect community structure abnormalities in depression, we constructed the Hotelling T-squared statistics as described in the method section for all 16 communities in the reference tree. Results showed a significantly lower consistency (p = 0.0079) in depression versus control for the community containing the right precuneus, superior parietal gyrus, inferior parietal gyrus, inferior temporal gyrus, and isthmus cingulate (figures 2 and and3).3). On a node level, conducting 2-sample T-tests or the MWW-test on *V ^{m}* or

Hotelling’s T-Squared statistics are overlaid on the 16 communities of the reference tree, showing regional group differences in the community structure between groups. The most significant differences are in the right precuenus and superior parietal **...**

4-level binary trees showing mean community structures via the proposed metric *Ψ*^{PL} for the healthy controls (left) and the depressed subjects (right). Notice that the right precuneus and right superior parietal gyrus belong to the same community **...**

In this paper we presented a new approach to construct hierarchical community structures from weighted structural brain networks. Our proposed metric took into account the path lengths within and between communities, and we used simulated annealing to construct top-down four-level binary trees. Moreover, we are the first to develop a statistical framework that allows for the quantification of node-level community structural differences between two groups of brain networks, via the consistency vector metric. This framework was validated on a sample of depressed patients. Results revealed local community differences in parts of the default mode network (especially the right precuneus), which has been previously linked to self-awareness and shown to exhibit abnormalities in depression [17, 18].

1. Sporns O, Tononi, Ktter R. The human connectome: A structural description of the human brain. PLoS Computational Biology. 2005;1:245–251. [PMC free article] [PubMed]

2. Rubinov M, Sporns O. Complex network measures of brain connectivity: Uses and interpretations. J of NeuroImage. 2010;52:1059–1069. [PubMed]

3. Bassett BS. Small World Brain Networks. The Neuroscientist. 2006;12(6):512–523. [PubMed]

4. Meunier D, Lambiotte R, Fornito A, Ersche KD, Bullmore E. Hierarchical modularity in human brain functional networks. Front. Neuroinformatics. 2009b;3:37. [PMC free article] [PubMed]

5. Telesford Q, Simpson S, Burdette J, Hayasaka S, Laurienti P. The Brain as a Complex System: Using Network Science as a Tool for Understanding the Brain. Brain Connectivity. 2011;1(4):295–308. [PMC free article] [PubMed]

6. Meunier D, Achard S, Morcom A, Bullmore E. Age-related changes in modular organization of human brain functional networks. J of Neuroimage. 2009a;44:715723. [PubMed]

7. Fair DA, Cohen AL, Power JD, Dosenbach NU, Church JA, Miezin FM, Schlaggar BL, Petersen SE. Functional brain networks develop from a “local to distributed” organization. PLoS Comput Biol. 2009;5(5):e1000381. [PMC free article] [PubMed]

8. Liu Y, Liang M, Zhou Y, He Y, Hao Y, Song M, Yu C, Liu H, Liu Z, Jiang T. Disrupted small-world networks in schizophrenia. Brain. 2008;131:945–961. [PubMed]

9. Bassett DS, Bullmore E, Verchinski BA, Mattay VS, Weinberger DR, Meyer-Lindenberg A. Hierarchical organization of human cortical networks in health and Schiz-ophrenia. J of Neuroscience. 2008;28(37):9239–9248. [PMC free article] [PubMed]

10. Lynall ME, Bassett DS, Kerwin R, McKenna PJ, Kitzbichler M, Muller U, Bullmore E. Functional connectivity and brain networks in schizophrenia. J of Neuroscience. 2010;30(28):9477–9487. [PMC free article] [PubMed]

11. Van Den Heuvel MP, Mandl RCW, Stam CJ, Kahn RS, Hulsho Pol HE. Aber-rant frontal and temporal complex network structure in schizophrenia: A graph theoretical analysis. J of Neuroscience. 2010;30(47):15915–15926. [PubMed]

12. Alexander-Bloch A, Lambiotte R, Roberts B, Giedd J, Gogtay N, Bullmore E. The discovery of population differences in network community structure: New methods and applications to brain functional networks in schizophrenia. J of NeuroImage. 2012;59(4):3889–3900. [PMC free article] [PubMed]

13. Newman MEJ, Girvan M. Finding and evaluating community structure in networks. Phys. Rev. E. 2004;69:026113. [PubMed]

14. Fortunato S. Community detection in graphs. Physical Reports. 2010;486(3-5):75–175.

15. Kirkpatrick S, Gelatt CD, Vecchi MP. Optimization by Simulated Annealing. Science. 1983;220(4598):671680. [PubMed]

16. Steen M, Hayasaka S, Joyce K, Laurienti P. Assessing the consistency of community structure in complex networks. Physical Review E. 2011;84:016111. [PMC free article] [PubMed]

17. Zhu X, Wang X, Xiao J, Liao J, Zhong M, Wang W, Yao S. Evidence of a Dissociation Pattern in Resting-State Default Mode Network Connectivity in First-Episode, Treatment-Naive Major Depression Patients. Biol. Psychiatry. 2011 [Epub ahead of print] [PubMed]

18. Lou HC, Luber B, Crupain M, Keenan JP, Nowak M, Kjaer TW, Sackeim HA, Lisanby SH. Parietal cortex and representation of the mental Self. Proc. of the National Academy of Sciences of the United States of America. 2004;101 [PubMed]

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