Research Article  Open Access
A Comparative Study of Theoretical Graph Models for Characterizing Structural Networks of Human Brain
Abstract
Previous studies have investigated both structural and functional brain networks via graphtheoretical methods. However, there is an important issue that has not been adequately discussed before: what is the optimal theoretical graph model for describing the structural networks of human brain? In this paper, we perform a comparative study to address this problem. Firstly, largescale cortical regions of interest (ROIs) are localized by recently developed and validated brain reference system named Dense Individualized Common Connectivitybased Cortical Landmarks (DICCCOL) to address the limitations in the identification of the brain network ROIs in previous studies. Then, we construct structural brain networks based on diffusion tensor imaging (DTI) data. Afterwards, the global and local graph properties of the constructed structural brain networks are measured using the stateoftheart graph analysis algorithms and tools and are further compared with seven popular theoretical graph models. In addition, we compare the topological properties between two graph models, namely, stickinessindexbased model (STICKY) and scalefree gene duplication model (SFGD), that have higher similarity with the real structural brain networks in terms of global and local graph properties. Our experimental results suggest that among the seven theoretical graph models compared in this study, STICKY and SFGD models have better performances in characterizing the structural human brain network.
1. Introduction
The human brain is intrinsically organized into distinct largescale functional networks, and cognitive functions arise from the dynamic interactions of distributed brain areas operating in these networks [1]. New advances in neuroimaging techniques have shown the possibility of systematic exploring the human brain formal complex network perspective. Graph theory provides a theoretical framework in which the topological properties of the brain networks can be examined such as centrality, clustering, efficiency, hierarchy, modularity, robustness, smallworldness, and synchronizability [2], and it can reveal important information about both the global and local organizations of the human brain networks. The improved characterization of brain networks achieved via graphtheoretical methods provides not only parsimonious accounts of normal cognitive processes [3], but also novel insights into psychiatric and neurological disorders such as Alzheimer’s disease [4, 5], multiple sclerosis [6], and attentiondeficit disorder [7].
Many complex systems show remarkably similar macroscopic behaviors despite profound differences in the microscopic details of the elements of each system or their mechanisms of interaction [2]. In this paper, we focus on an important issue that has not been adequately addressed before: what graph models can best possibly describe the structures of brain networks. In this way, the model provides a possible explanation for a key challenge for systems neuroscience: how to understand the complex network organization of the brain on the basis of neuroimaging data [8]. For instance, previous studies have demonstrated the smallworldness [9] and scale free properties of human brain networks [10].
In brain network studies, the methods for identifying regions of interest (ROIs), that is, the network nodes, can be generally classified into four categories. The first group is manual labeling by experts based on experience and domain knowledge. This method is widely used; however, it may not be reproducible due to both intersubject and intrasubject variations [11]. The second group is datadriven methods, which clusters ROIs from the brain image itself [12, 13]. However, it might be sensitive to the clustering parameters used in many datadriven approaches. The third group of methods identifies activated brain regions as ROIs by taskbased functional magnetic resonance imaging (taskbased fMRI), and it is regarded as the benchmark approach for ROI identification. However, taskbased fMRI itself has limitations such as being timeconsuming and expensive. Additionally, different patterns may be shown in groupbased activation maps from an individual’s activation map [14]. In short, it remains quite challenging to accurately localize ROIs for each individual by using standard analysis of taskbased fMRI data [15]. The last group methods are cortical parcellation based on image/surface registration, whose limitations have been comprehensively discussed in [16, 17].
In this paper, we apply our recently developed brain reference system, named Dense Individualized and Common Connectivitybased Cortical Landmarks (DICCCOL) [18] which discovers 358 consistent and corresponding ROIs across subjects based on diffusion tensor imaging (DTI) data, to localize the ROIs for each participated subject. DICCCOL possess intrinsically established structural and functional correspondences (universal), while their locations and sizes are determined in each individual’s space (individualized). With the identified brain ROIs, the structural brain network is constructed for each subject based on the corresponding DTI data. Then, through the large network analysis tool GraphCrunch2 [19], we evaluate the fitness of seven popular theoretical graph models for describing the real brain network by measuring both the global and local graph properties of the constructed brain networks and compare them with those graph models. The graph models are as follows: (1) ErdősRényi random graph (ER) [20]; (2) ErdősRényi random graph with the same degree distribution as the input data (ERDD) [19, 21]; (3) geometric random graph (GEO) [22, 23]; (4) geometric gene duplication model (GEOGD) [24]; (5) scalefree BarabásiAlbert preferential attachment model (SF) [25]; (6) scalefree gene duplication model (SFGD) [26]; (7) stickinessindexbased model (STICKY) [27]. Those graph models will be explained in detail in Section 2.
Our experimental results suggest that the SFGD model fits the real brain network the best in terms of global and local graph properties. We also demonstrate that the real brain network also has the STICKY property. In summary, the SFGD model, combined with its STICKY property, can best describe graph properties of the real brain network and also can indirectly describe the mechanism of structural network from biological properties. Importantly, the results are consistent across populations.
2. Materials and Methods
2.1. Data Acquisition and Preprocessing
DTI datasets for 104 healthy subjects including three age groups of adolescents (28), young adults (53), and elderly normal brains (23) were acquired on a 3T GE Signal magnetic resonance imaging (MRI) scanner. Acquisition parameters for the scans were as follows: 256 × 256 matrix, 3 mm slice thickness, 240 mm^{2} field of view (FOV), 50 slices, 15 diffusion weighted imaging (DWI) volumes, and value = 1000. The preprocessing of DTI data included brain skull removal, motion correction, and eddy current correction [28]. After preprocessing, fiber tracking was reconstructed via MEDINRIA [29]. Then, The grey matter (GM)/white matter (WM) cortical surface was reconstructed according to the brain tissue segmentation map based on the DTI data [30].
2.2. Structural Brain Network Construction
The structural brain network of each subject is represented as an unweighted undirected graph in which is the set of nodes and is the set of edges, where is the number of nodes in the network. The nodes in are identified via DICCCOL system [18], and the elements in are measured for each ROI pair independently. Specifically, DICCCOL employs a novel datadriven strategy to discover, from DTI datasets, dense and common cortical landmarks [18]. The basic idea is that we optimize the localizations of each DICCCOL landmark in individual subjects by maximizing the groupwise consistency of their white matter fiber connectivity patterns. We obtain 358 DICCCOL ROIs for each subject and regard them as network nodes. White matter fibers obtained via the deterministic DTI tractography are projected onto the reconstructed cortical surface and the number of fibers connecting the ROI pair is used to measure the structural connectivity between the ROI pairs. Thus, a structural connectivity matrix is obtained for each subject. A predefined threshold is applied on the connectivity matrix to remove the noise and errors in fiber tracking. Note that is the same for all the subjects. After thresholding the structural connectivity matrix, we obtained the adjacency matrix for each subject. In the adjacency matrix, two nodes and are connected if ; otherwise, . The selfloops are currently ignored in the constructed structural brain networks.
Following the procedure described above, we constructed structural brain networks for 28 adolescents, 53 adults and 23 elders. Figure 1 shows exemplar brain networks in different groups.
(a) Adolescent
(b) Adult
(c) Elder
2.3. Graph Models
In this study, we compared the topological properties of the constructed structural brain networks with those of 7 typical graph models via GraphCrunch2 [19], a tool for complex network analysis. GraphCrunch2 implements the following network models: (1) ErdősRényi random graph model (ER): ER graph is generated by using the Library of Efficient Data Types and Algorithms (LEDA) based random graph generator [19]. ER graph can provide a rigorous definition of what it means for a property to hold for almost all graphs, or be used in the probabilistic method to prove the existence of graphs satisfying various properties [20]. (2) ErdősRényi random graph with the same degree distribution as the input data (ERDD): ERDD graph is generated by using the “stubs method” [19]. In brief, according to the degree distribution of the realworld network being modeled, the number of “stubs” (to be filled by edges) is assigned to each node in the model network. After that, edges are created between randomly picked pairs of nodes. At last, the number of “stubs” left available at the corresponding “endnodes” of the edge is decreased by one [19, 21]. (3) Geometric random graph model (GEO): in GEO graph, nodes correspond to uniformly randomly distributed points in a metric space and if the corresponding points are close enough in the metric space, according to some distance norm, the edges are created between pairs of nodes [22, 23]. (4) Geometric gene duplication model (GEOGD): GEOGD graph is the extension of geometric random graph, in which the principles of gene duplications and mutations are incorporated [19]. Every model determines the principle by which the network is grown from a small seed network and adds new nodes intended to model gene duplications and mutations [24]. (5) Scalefree BarabásiAlbert preferential attachment model (SF) [25]: the most important characteristic in a scalefree network is the relative commonness of vertices with a degree that greatly exceeds the average. The highestdegree nodes are named “hubs.” The scalefree property strongly correlates with the robustness of network. The hierarchy allows for a faulttolerant behavior [26]. (6) Scalefree gene duplication model (SFGD): SFGD model is an evolution of SF model. The tolerance of SFGD to damage is determined by the scalefree nature of its multifractal distribution. However, they present novel properties: multifractal features are inherited in a model of growing networks [26]. (7) Stickinessindexbased model (STICKY): STICKY graph is based on stickiness indices that summarize node connectivities and the complexities of normalized degree of nodes in networks. The stickiness framework produces a convenient, parameterfree random network [27].
According to the input structural brain network, GraphCrunch2 repeatedly generates a number of instances of the defined graph models and measures the global and local graph properties that infer the similarity/dissimilarity between the input graph and each of the generated model graph instances. In our experiments, the number of repeatedly generated instances is 10. The global and local graph properties implemented in GraphCrunch2 will be introduced in the next section.
2.4. Topological Properties of Graphs
We calculated both global and local topological graph properties for the constructed structural brain networks and the 7 graph models. Statistical analysis was conducted to explore which graph model can best describe the constructed brain networks. In this paper, the following global properties are examined: (1) Pearson correlation coefficients between the degree distributions [31]; (2) average shortest path length difference ratio; (3) average clustering coefficient difference ratio. The local graph properties include the relative graphlet frequency (RGF) distance [23] and the graphlet degree distribution (GDD) agreement [32]. We briefly provide the definitions of these measurements below.
2.4.1. Global Graph Properties
Define the degree of a node as the number of neighbors it has in the network. The topological property of a graph can be obtained in terms of the degree distribution [2]. The degree distribution is one of indirect measures that reflect the network robustness to insult [8]. We compute Pearson correlation coefficients between the degree distributions of the graph model and the real brain network as the similarities between them in the aspect of degree distributions. Shortest path length plays an important role in characterizing the internal structure of a graph. Shortest path length of a node pair is defined as the path between two nodes in a graph that the number of its constituent edges is minimized. The average of the shortest path length, also known as characteristic path length, is the mean of the shortest paths over all node pairs [2] and is the most widely used measure to functional integration [8]. The clustering coefficient of node is defined as , where is the number of edges between neighbors of . The average clustering coefficient is the mean of over all nodes [2]. In this paper, we use average shortest path length difference ratio (PathDiff) and average clustering coefficient difference ratio (ClustDiff) to measure the difference between graph models and the real brain networks: where and Path_{real} are the average of shortest path of models and real brain networks, respectively and and are the average of clustering coefficient of models and real brain networks, respectively. Smaller difference indicates higher similarity between two networks.
2.4.2. Local Graph Properties
In addition to above global properties, we also measured the local graph properties to evaluate the fitness of the graph models to the real brain networks. In this paper, the local graph properties include RGF distance and GDD agreement.
Graphlet degree distribution agreement (GDD agreement) [32] is a similarity measure between topologies of two networks based on graphlet degree vector distributions. GraphCrunch2 computes 2–5 nodes graphlets [19]. By calculating the fitness of each of the 73 GDDs of the networks under comparison, GDD agreement contains 73 similarity constraints. High GDD agreement between two networks indicates that they are similar [19]. In general, GDD agreement is a local heuristic metric for measuring network structure. It imposes 73 highly structured constraints, and thus it increases the chances that two networks are truly similar if they are similar with respect to this measure [32].
Relative graphlet frequency distance (RGF distance) [23] is a measure that compares the frequencies of appearance of all 2 to 5 node graphlets in two networks [19]. Since there are 30 possible graphlets on up to 5 nodes, RGF distance includes 30 similarity constraints by examining the fit of 30 graphlets frequencies between two networks. The similarity between two graphs only depends on the differences between relative frequencies of graphlets. Smaller RGF distance indicates higher similarity between networks [19].
3. Results and Discussion
3.1. Topological Graph Properties
3.1.1. Global Graph Properties
Figure 2(a) shows the Pearson correlation coefficients between the degree distributions of the constructed brain networks and the 7 graph models for 28 subjects in the adolescent group. Note that the ERDD model needs to be excluded because it has the same degree distribution as the input graph. It is seen that the Pearson correlation of degree distribution is approaching to 1. Figure 2(a) indicates that the STICKY model and SFGD model have a high correlation with the real structural brain networks. The difference ratio of average shortest path length and average clustering coefficient are shown in Figures 2(b) and 2(c), respectively. It is seen that the SFGD and STICKY models have lower average shortest path length difference ratio and average clustering coefficient difference ratio.
(a)
(b)
(c)
In general, the comparison of global topological properties between real structural brain networks and the 7 graph models demonstrates the superiority of SFGD and STICKY graph models in characterizing the structural brain networks. SFGD and STICKY graph models share common features, namely, “inheritance” and “variance”: when the new nodes are created, the son nodes inherit the connectivities of their parent nodes (“duplication” and “stickiness”). Meanwhile, the son nodes have new connectivities with other nodes which are not linked with their parent nodes “divergence.” These characteristics could indirectly describe the biologically important properties of the mechanism of structural brain network. In the process of brain development, the new neurons are created and contain both functional and structural features. In this way, the specific regions of the brain which have specific function are formed, and the functions of these regions are improved and enhanced, such as the visual association area, motor speech area, and olfactory area. Meanwhile, the new neurons brought additional connectivities with other neurons in different regions; therefore, they enhanced the cooperation and coordination ability of different regions in human brain. The “stickiness” indicates that the important neurons which have large degree in structural brain network will maintain the original functional and structural features during brain development. On the other hand, the SFGD network shares common features with other scalefree networks, and its tolerance to damage is determined by the scalefree nature of its multifractal distribution. Therefore, it also means that, during development, the brain not only maintains the functional integrity, but also is with increased stability. The SFGD model has the lowest average shortest path length difference ratio indicates that “divergence” feature enhanced the interactions among different regions during brain development.
3.1.2. Local Graph Properties
The RGF distance and GDD agreement of the 7 graph models, in comparison with the real structural brain networks for the adolescent group with 28 subjects, are shown in Figures 3(a) and 3(b), respectively. It is seen that the GEOGD model has the lowest GDD agreements when compared with the real brain networks, and there is no significant difference among the rest of the graph models. It is also shown that STICKY, SF, and SFGD graph models have lower RGF distance, indicating their superior performance in characterizing the real structural brain networks. It is notable that GraphCrunch2 only calculates 2–5 nodes graphlets of the networks, as mentioned previously. In real brain networks, the connectivity patterns are complex and there exist graphlets with size over 5 nodes. Nevertheless, the local graph properties of RGF distance and GDD agreement demonstrate that STICKY, SF, and SFGD graph models have higher similarity to the real structural brain networks, especially for the STICKY model. This result reflects that the “stickiness” feature plays an important role in brain development, and it might be the main reason that brain maintains regional and group integrity.
(a)
(b)
In summary, the results of global and local topological properties indicate that the SFGD and STICKY graph model fit the brain network data better; that is, the SFGD and STICKY models have better performance in the description of structural brain networks. Since SFGD and STICKY are not in conflict [27], we can also regard the stickiness as a property of the real brain network.
3.2. SFGD versus Sticky
The global and local topological properties presented in the previous section indicate that SFGD and STICKY graph models have comparable performance in describing the real structural brain networks. In this section, we present a further comparison study between SFGD and STICKY by evaluating two topological properties of complex graph which were not implemented in GraphCrunch2, namely, (1) smallworldness [9]; (2) global efficiency [2]. Smallworldness has been reported as one of the most important properties of both the functional and structural human brain networks in previous studies. Smallworldness networks are defined as networks that are significantly more clustered than random networks, yet have approximately the same characteristic path length as random networks. More generally, smallworld networks should be simultaneously highly integrated and segregated [33]. This property is often analyzed by considering the fraction of nodes in the network that have a particular number of connections going into them. The average inverse shortest path length is a related measure known as the global efficiency. Some authors have claimed that the global efficiency may be a superior measure of integration [34]. Unlike the characteristic path length, the global efficiency can be more meaningfully computed on disconnected networks. Zero efficiency corresponds to infinite length of paths between disconnected nodes [8]. Figure 4 shows the mean of smallworldness and global efficiency in the adolescent group. It indicates that the SFGD model has lower discrepancy compared with the real brain networks in both smallworldness and global efficiency. This result to some extent demonstrates that SFGD graph model is relatively more suitable in characterizing the human brain networks than STICKY graph model.
(a)
(b)
3.3. Topological Graph Properties across Groups
As a reproducibility study, we compared the topological properties of the brain networks across three different groups (adolescents, adults, and elders). The results are shown in Tables 1, 2, 3, 4, and 5. The trend is almost the same and the difference across the three groups is relatively small. It is also seen that the SFGD model has higher correlation coefficients of degree distributions, lower average shortest path length difference ratio, and average clustering coefficient difference ratio in all the three groups. It indicates that SFGD model is the most similar network to the real structural brain networks. The STICKY model also has higher similarity to real brain networks in terms of GDD agreement and RGF distance.





4. Conclusion
In this paper, we compared the global and local graph properties, as well as the topological properties of structural human brain networks to 7 representative graph models. The objective is to explore which graph model can best describe structural human brain networks.
Our experimental results demonstrated that SFGD graph model in general has the best performance in characterizing the structural networks of the human brain, followed by STICKY graph model. SFGD graph model is based on the hypothesis of evolution by duplications and divergence of the genes which produce proteins [26]. It reproduces the topological properties of the proteinprotein interaction networks (PIN) with noticeable accuracy. Our experimental results also showed that the STICKY model has high GDD agreement, high Pearson correlation coefficients of the degree distributions, and low RGF distance and average shortest path length difference ratio when compared with SFGD graph model. The STICKY model was also widely used in the investigation of proteinprotein interaction networks, and the stickiness index can essentially capture the abundance and popularity of binding domains on a protein.
Since SFGD and STICKY are not in conflict [27], we can also regard the stickiness as a property of the real brain network. In SFGD model, networks revolute by duplication of nodes, and as a node is duplicated, it inherits most of the neighbors (interactions) of its parent nodes but gains some new neighbors as well [26]. Combined with stickiness index, a duplicated node would inherit its parents’ stickiness index along with many of the parents’ neighbors, and it would gain new neighbors in proportion to its inherited stickiness index and stickiness indices of the nodes already in the network [27]. The characteristics of “duplication,” “divergence,” and “stickiness” may explain the mechanism of structural brain network from biologically important properties shaped by natural selection. The “duplication” and “stickiness” make sure that the new neurons are created and contain both the functional and structural features. On the other hand, the “divergence” ensures that the new neurons have new connectivities with other neurons in different regions in order to enhance the cooperation and coordination ability of different regions in the human brain. This provides explanation why both STICKY and SFGD models are suitable for describing the human brain networks; that is, the human brain networks are similar to SFGD model and have STICKY property at the same time.
References
 M. D. Fox, A. Z. Snyder, J. L. Vincent, M. Corbetta, D. C. van Essen, and M. E. Raichle, “The human brain is intrinsically organized into dynamic, anticorrelated functional networks,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 27, pp. 9673–9678, 2005. View at: Publisher Site  Google Scholar
 E. Bullmore and O. Sporns, “Complex brain networks: graph theoretical analysis of structural and functional systems,” Nature Reviews Neuroscience, vol. 10, no. 3, pp. 186–198, 2009. View at: Publisher Site  Google Scholar
 S. L. Bressler and V. Menon, “Largescale brain networks in cognition: emerging methods and principles,” Trends in Cognitive Sciences, vol. 14, no. 6, pp. 277–290, 2010. View at: Publisher Site  Google Scholar
 Y. He, Z. Chen, and A. Evans, “Structural insights into aberrant topological patterns of largescale cortical networks in Alzheimer's disease,” Journal of Neuroscience, vol. 28, no. 18, pp. 4756–4766, 2008. View at: Publisher Site  Google Scholar
 C. J. Stam, B. F. Jones, G. Nolte, M. Breakspear, and P. Scheltens, “Smallworld networks and functional connectivity in Alzheimer's disease,” Cerebral Cortex, vol. 17, no. 1, pp. 92–99, 2007. View at: Publisher Site  Google Scholar
 Y. He, A. Dagher, Z. Chen et al., “Impaired smallworld efficiency in structural cortical networks in multiple sclerosis associated with white matter lesion load,” Brain, vol. 132, no. 12, pp. 3366–3379, 2009. View at: Publisher Site  Google Scholar
 L. Wang, C. Zhu, Y. He et al., “Altered smallworld brain functional networks in children with attentiondeficit/hyperactivity disorder,” Human Brain Mapping, vol. 30, no. 2, pp. 638–649, 2009. View at: Publisher Site  Google Scholar
 E. Bullmore, A. Barnes, D. S. Bassett et al., “Generic aspects of complexity in brain imaging data and other biological systems,” NeuroImage, vol. 47, no. 3, pp. 1125–1134, 2009. View at: Publisher Site  Google Scholar
 D. J. Watts and S. H. Strogatz, “Collective dynamics of “smallworld” networks,” Nature, vol. 393, no. 6684, pp. 440–442, 1998. View at: Google Scholar
 R. Cohen and S. Havlin, “Scalefree networks are ultrasmall,” Physical Review Letters, vol. 90, no. 5, Article ID 058701, 2003. View at: Google Scholar
 K. Amunts, A. Malikovic, H. Mohlberg, T. Schormann, and K. Zilles, “Brodmann's areas 17 and 18 brought into stereotaxic space: where and how variable?” NeuroImage, vol. 11, no. 1, pp. 66–84, 2000. View at: Publisher Site  Google Scholar
 C. F. Beckmann, M. DeLuca, J. T. Devlin, and S. M. Smith, “Investigations into restingstate connectivity using independent component analysis,” Philosophical Transactions of the Royal Society B, vol. 360, no. 1457, pp. 1001–1013, 2005. View at: Publisher Site  Google Scholar
 Y. Zang, T. Jiang, Y. Lu, Y. He, and L. Tian, “Regional homogeneity approach to fMRI data analysis,” NeuroImage, vol. 22, no. 1, pp. 394–400, 2004. View at: Publisher Site  Google Scholar
 K. Li, L. Guo, D. Zhu, X. Hu, J. Han, and T. Liu, “Individual functional ROI optimization via maximization of groupwise consistency of structural and functional profiles,” Neuroinformatics, vol. 10, no. 3, pp. 225–242, 2012. View at: Publisher Site  Google Scholar
 K. Li, L. Guo, C. Faraco et al., “Visual analytics of brain networks,” NeuroImage, vol. 61, no. 1, pp. 82–97, 2012. View at: Publisher Site  Google Scholar
 X. Hu, D. Zhu, P. Lv, K. Li et al., “Finegranularity functional interaction signatures for characterization of brain conditions,” Neuroinformatics, vol. 11, no. 3, pp. 301–317, 2013. View at: Publisher Site  Google Scholar
 T. Liu, “A few thoughts on brain ROIs,” Brain Imaging and Behavior, vol. 5, no. 3, pp. 189–202, 2011. View at: Publisher Site  Google Scholar
 D. Zhu, K. Li, L. Guo et al., “DICCCOL: dense individualized and common connectivitybased cortical landmarks,” Cerebral Cortex, vol. 23, no. 4, pp. 786–800, 2013. View at: Publisher Site  Google Scholar
 O. Kuchaiev, A. Stevanović, W. Hayes, and N. Pržulj, “GraphCrunch 2: software tool for network modeling, alignment and clustering,” BMC Bioinformatics, vol. 12, article 24, 2011. View at: Publisher Site  Google Scholar
 P. Erdos and A. Renyi, “On the evolution of random graphs,” Publicationes Mathematicae, vol. 6, pp. 290–297, 1959. View at: Google Scholar
 M. Molloy and B. Reed, “A critical point of random graphs with a given degree sequence,” Random Structures and Algorithms, vol. 6, no. 23, pp. 161–180, 1995. View at: Google Scholar
 M. Penrose, Random Geometric Graphs, Oxford University Press, New York, NY, USA, 2003.
 N. Przulj, D. G. Corneil, and I. Jurisica, “Modeling interactome: scalefree or geometric?” Bioinformatics, vol. 20, no. 18, pp. 3508–3515, 2004. View at: Publisher Site  Google Scholar
 N. Przulj, O. Kuchaiev, A. Stevanovic, and W. Hayes, “Geometric evolutionary dynamics of protein interaction networks,” in Proceedings of the Pacific Symposium on Biocomputing, pp. 178–189, Stanford, Calif, USA, 2010. View at: Google Scholar
 A.L. Barabasi and R. Albert, “Emergence of scaling in random networks,” Science, vol. 286, no. 5439, pp. 509–512, 1999. View at: Publisher Site  Google Scholar
 A. Vazqueza, A. Flamminia, A. Maritana, and A. Vespignani, “Modeling of protein interaction networks,” Complexus, vol. 1, no. 1, pp. 38–44, 2003. View at: Google Scholar
 N. Przulj and D. J. Higham, “Modelling proteinprotein interaction networks via a stickiness index,” Journal of the Royal Society Interface, vol. 3, no. 10, pp. 711–716, 2006. View at: Publisher Site  Google Scholar
 D. Zhang, L. Guo, G. Li et al., “Automatic cortical surface parcellation based on fiber density information,” in Proceedings of the 7th IEEE International Symposium on Biomedical Imaging (ISBI '10), pp. 1133–1136, Xi'an, China, April 2010. View at: Publisher Site  Google Scholar
 http://wwwsop.inria.fr/asclepios/software/MedINRIA/.
 T. Liu, H. Li, K. Wong, A. Tarokh, L. Guo, and S. T. C. Wong, “Brain tissue segmentation based on DTI data,” NeuroImage, vol. 38, no. 1, pp. 114–123, 2007. View at: Publisher Site  Google Scholar
 C. Bishop, Pattern Recognition and Machine Learning, Springer, Cambridge, Mass, USA, 2006.
 N. Przulj, “Biological network comparison using graphlet degree distribution,” Bioinformatics, vol. 23, no. 2, pp. e177–e183, 2007. View at: Publisher Site  Google Scholar
 M. D. Humphries and K. Gurney, “Network “smallworldness”: a quantitative method for determining canonical network equivalence,” Plos ONE, vol. 3, no. 4, Article ID e2051, 2008. View at: Publisher Site  Google Scholar
 S. Achard and E. Bullmore, “Efficiency and cost of economical brain functional networks,” Plos Computational Biology, vol. 3, no. 2, article e17, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Xiaojin Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.