Construction of Individual Morphological Brain Networks with Multiple Morphometric Features

Introduction

As a sophisticated but well-organized system, the human brain is capable of processing multiple tasks with high efficiency (Sporns et al., 2004). The interconnected regions of the brain are apt to delineate the operation pattern of integration and segregation, which indicating the ability to rapidly organize specialized information from distributed brain regions functionally, while specialized processing to occur within densely interconnected groups of brain regions (Rubinov and Sporns, 2010; Sporns, 2013). Recently, the brain network has been applied to characterize the effects of age, gender (Damoiseaux et al., 2008; Gong et al., 2009b; Chen et al., 2011; Yan et al., 2011) and psychiatric disorders (Bassett et al., 2008; He et al., 2009; Van Den Heuvel et al., 2010; Alexander et al., 2012; Shi et al., 2013).

Brain networks are generally categorized into three modalities through different neuroimaging and electrophysiological techniques. Functional networks are formed using functional magnetic resonance imaging (fMRI) (Eguiluz et al., 2005; Bassett and Bullmore, 2006), electroencephalogram (EEG) (Micheloyannis et al., 2006; Stam et al., 2007), or magnetoencephalogram (MEG) (Stam, 2004). An anatomical network is built using diffusion tensor imaging (DTI) (Hagmann et al., 2007; Gong et al., 2009a), and a morphological network is constructed using structural MRI (He et al., 2007).

Despite the remarkable progress of studies on functional and anatomical connectomes, the morphological brain network issue was raised relatively late (Evans, 2013). He et al. (2007) used cortical thickness as the regional descriptor to realize the construction of a network by computing the correlations across regions across a large scale; their study is the first to document that the morphological brain network exhibits small-worldness and a scale-free degree distribution. This methodology has also been applied to different descriptors such as regional gray matter volume (Bassett et al., 2008) and surface area (Sanabria-Diaz et al., 2010). Moreover, the morphological brain network can be employed to explore the hierarchical and modular organizations of the cerebral cortex (Bassett et al., 2008; Chen et al., 2011).

Nevertheless, population-based studies on the morphological brain network result in a severe loss of information of inter-individual differences (Kanai and Rees, 2011). Fortunately, several novel approaches have been recently developed to extract structural information directly from T1-weighted MR images and obtain the interregional connectivities for a single subject (Raj et al., 2010; Zhou et al., 2011; Tijms et al., 2012; Kong et al., 2015). For instance, Raj et al. (2010) proposed an individual-based morphological network construction method by using Gibbs probability models. Similarly, in another study, the probability density function of local morphological features was used as the regional descriptor to build the network (Kong et al., 2015). In addition to these mathematical methods, Tijms et al. (2012) employed a cube-based approach. In this approach, after the intensity of voxels is concatenated into a feature vector for each cube, the inter-cube connectivities can be derived individually as the correlation coefficient. Notably, in this method, the size of the node in the brain network has shrunk from the atlas-based region to the voxel-based cube. Furthermore, Batalle et al. (2013) presented a method to normalize the atypical morphological networks, such as the cube partition (Tijms et al., 2012), to standardized brain networks with nodes identified using a parcellation scheme (i.e., AAL brain atlas).

Despite the breakthroughs in research on individual-based morphological brain networks, one particular aspect draws our attention. Previous studies have built their morphological networks by using a single cortical feature, such as cortical thickness or voxel intensity. However, there are multiple morphometric features can be extracted from the brain structural organization. Therefore, we extend previous studies by proposing a novel idea for realizing the formation of individual-based morphological brain network using a combination of multiple morphometric features. In the network construction, the interregional connections were computed as the correlation of feature vectors instead of one feature; each of these vectors comprised nine complementary morphometric features.

To test this idea, we applied it to a sample of 55 healthy participants (24 males) aged 20–29 years old. The network nodes were determined using a brain atlas with multiple anatomical regions of interest. The global and local topological properties of the network were computed in accordance with graph theory. The results of small-world configurations were compared with those of previous studies at a similar sparsity level. Furthermore, the nodal betweenness centrality and the spatial distribution of hubs were thoroughly investigated. Finally, the reproducibility of the method was evaluated.

Materials and Methods

Participants

Fifty-five right-handed healthy subjects (24 male ages ranging from 20 to 29 with mean = 22.92 and standard deviation = 2.89; 31 female ages ranging from 20 to 28 with mean = 21.71 and standard deviation = 2.19) were selected from the Open Access Series of Imaging Studies Database (http://www.oasis-brains.org/). All the subjects were questioned about their medical histories and use of psychoactive drugs before a trained physician began image acquisition. For details on the clinical and demographic information of the subjects, please refer to Marcus et al. (2007).

Image Acquisition

For each subject, three or four individual T1-weighted magnetization-prepared rapid gradient-echo images were acquired using a single Siemens 1.5 T Vision scanner with the following parameters: repetition time = 9.7 ms, echo time = 4.0 ms, inversion time = 20 ms, flip angle = 10°, sagittal orientation with 128 slices, and resolution = 1 × 1 × 1.25 mm. Multiple T1 images obtained for each subject were motion-corrected and then averaged to achieve an improved signal-to-noise ratio. For additional details on the post-processing information regarding the raw images, please refer to Marcus et al. (2007).

Measurements of Multiple Morphometric Features

All the images were preprocessed through FreeSurfer (version 5.3.0), which can be freely downloaded online (http://surfer.nmr.mgh.harvard.edu/). Studies have assessed FreeSurfer's performance with images acquired from various MRI scanners or sequences and received reliable results for cortical and subcortical measurements (Khan et al., 2008; Shi et al., 2013; Mulder et al., 2014). Details on FreeSurfer's processing pipeline can be found in Dale et al. (1999) and Fischl et al. (1999). In short, the raw images were resampled into 256 × 256 × 256 isotropic resolutions with a voxel size of 1 × 1 × 1 mm. The intensity bias was subsequently corrected in a volume-based stream, and the skull was stripped, followed by volumetric labeling and white matter segmentation (Figure 1A). In the surface-based stream, the surface of white matter was extracted as the inner surface, and then nudged to the gray–pial interface to generate the outer surface (Figure 1B). The FreeSurfer surface was created as a mesh with a non-uniform grid, also known as a vertex. Hence, the cortical morphometric features can be measured using every vertex between the inner and outer surfaces (Figure 1C). Finally, the built-in Desikan-Killiany cortical atlas (Desikan et al., 2006) was applied to obtain the regional measurements (Figure 1D). This atlas presents the parcellation of 34 regions for each hemisphere on the basis of the structural patterns of the gyrus and sulcus (Supplementary Table 1).

FIGURE 1

Figure 1. The general flowchart of the individual-based morphological brain network construction and analyses. The steps from (A–D) were accomplished using FreeSurfer. Intensity bias was corrected, and the skull was stripped (A) before the inner and outer surfaces of gray matter were extracted (B). All the involved morphometric features were first measured between inner and outer surfaces (C), and then the feature results were mapped to the Desikan–Killiany cortical atlas to obtain the regional feature measurements (D). For each region, nine morphometric features were concatenated as a feature vector to obtain the correlation coefficient between any two regions to frame the individual-based morphological connectivity (E). The connection matrix was then repeatedly binarized based on the sparsities (from 20 to 40% with a step size of 1%) to generate the network graphs (F), and then the network topological properties were analyzed in accordance with graph theory (G).

Nine morphometric features were initially employed in this study (Figure 1E): (1) the number of vertices in each region, (2) regional gray matter volume, (3) surface area, (4) mean cortical thickness, (5) standard deviation of cortical thickness, (6) mean curvature, (7) Gauss curvature, (8) curvature index, and (9) fold index. In particular, the first two features are tallied as the sum of the vertices and voxels in a region, respectively. Surface area was computed as the total of the areas of all vertices falling within a region (Panizzon et al., 2009). The cortical thickness is measured as the distance between the inner and outer cortical surfaces at each vertex (Fischl et al., 1999). The curvature is calculated as the reciprocal of the radius of the inscribed circle for each vertex, and the signs of the gyrus and sulcus are opposite (Li et al., 2014). The Gaussian curvature of a surface at a given point is the product of the principal curvatures, which are the eigenvalues of the shape operator at the point (Pienaar et al., 2008). They measure how the surface bends by different amounts in different directions at that point. The folding index is calculated as the number of cortices buried within the sulcal folds, as compared with the number of cortices on the visible outer cortex (Schaer et al., 2008). Notably, all the surface computations were performed in the native space, enabling the above-mentioned features to be measured without deformation.

Individual-Based Morphological Network Construction

To form the network, the regional descriptor was defined as a feature vector, which is the concatenation of the nine morphometric features in each region, as shown in Figure 1. Therefore, the interconnected matrix is generated by computing the Pearson correlation coefficient for each pair of feature vectors (Figure 1E). Significant differences exist between the orders of magnitude of each morphometric feature (10⁻² to 10³). Thus, the z-score values were computed for each feature as the standardized values before the correlation calculation. The z-score values were calculated from each morphometric feature by first subtracting the mean of all values from each individual value and then dividing each remainder by the standard deviation of all the values, as realized by the normalization function in the statistical analysis software SPSS v22.0 (SPSS Inc., Chicago IL, USA). The correlation matrix was then binarized to generate the unweighted and undirected networks (Figure 1F), because this matrix captures the underlying anatomical connection patterns, while providing a simple method for investigating the network (He et al., 2007). To maintain both positive and negative high correlations, the absolute values of matrixes were applied for binarization. Lastly, such network properties as small-world configurations and definition of hubs (Figure 1G) were analyzed in accordance with graph theory.

For binarization, sparsity was used to express the extent of thresholding, which is the ratio of existing connections to the total possible ones. However, the appropriate sparsity remains an open question (Zhu et al., 2012). Therefore, the network properties were evaluated as a function of sparsity ranging from 20 to 40% with a step size of 1%. The range was determined using the requirements stating that the minimum sparsity guarantees no isolated node in the network, and the maximum one ensures small-worldness. The other advantage of using a certain range of sparsity is that it allows us to investigate the network properties at the same level for all the subjects, as the same number of existing connections are found in each case. Additionally, the sparsity of 23% was highlighted to enable direct comparison with previous studies.

Analyses of Network Properties

The network properties were computed using the Brain Connectivity Toolbox (BCT) (Rubinov and Sporns, 2010) and Graph-theoretical Network Analysis (GRETNA) toolkit (Wang et al., 2015). SPSS was used for the single-feature normalization (see individual-based morphological network construction) and all the statistical analyses, while the spatial distribution of hubs was visualized using the BrainNet Viewer toolkit (Xia et al., 2013). The multiple comparison was corrected using the false discovery rate (FDR) proposed by Genovese et al. (2002).

Small-World Configurations

In 1998, Watts and Strogatz first proposed the notion of a small-world network, which exhibits a similar characteristic path length (L_p) but a higher clustering coefficient (C_p) than a random network (Watts and Strogatz, 1998). In particular, C_p represents the average of the clustering coefficients of all the nodes (C_i) in a network, where C_i expresses the likelihood that the neighborhoods of the nodes are connected. Specifically, C_i of a given node is computed as the proportion of connections among its neighbors which are actually realized compared with the number of all possible connections. Thus, C_p quantifies the extent of local cliquishness or the efficiency of a particular network's information transfer (Latora and Marchiori, 2001). In addition, L_p represents the characteristic path length of the graph, which is the average shortest path length among all pairs of nodes in the network. However, correct calculation of L_p based on the definition is unachievable because not every pair of nodes in a binary brain network is connected. To avoid this situation, L_p is measured as the “harmonic mean” distance between pairs of nodes proposed by Newman (2004). L_p quantifies the parallel information propagation ability, or the global efficiency of a network (Latora and Marchiori, 2001). Therefore, small-worldness can be demonstrated mathematically as

$\gamma =\frac{C_p^{observed}}{C_p^{random}}>1\text{ and}\lambda =\frac{L_p^{brain}}{L_p^{random}}\approx 1$

or merged into one formula

$\sigma =\frac{\gamma }{\lambda }>1,$

where the corresponding random networks consist of the same number of nodes and edges (Achard et al., 2006).

Spatial Distribution of Hubs

The betweenness centrality (BC) is defined as the number of shortest paths between any two nodes running through the given node, and measures the nodal ability of information flow throughout the network (Freeman, 1977). In particular, two considerations are adopted to determine the hubs. The first specifies the hubs for each subject; thus, the average of BC for the entire sparsity range (mspBC) is used to represent nodal BC for each subject. The next consideration involves realizing the identification of hubs at each sparsity; hence, the mean BC of all the subjects (msjBC) was employed to represent nodal BC at each sparsity. Overall, the hubs were selected with the same criterion, indicating that the nodes should achieve a higher BC (mspBC or msjBC) than the sum of the mean and standard deviation for the entire network, as demonstrated below.

$BC>mea{n}_{BC}+standarddeviatio{n}_{BC}$

The nodal degree was not used to determine the hubs because in contrast with BC, this degree only measures the connections linked with the node, rather than the shortest path.

Reproducibility

To evaluate the reproducibility of morphometric features and the present method, the intraclass correlation coefficient (ICC) was applied. ICC was first defined by Mcgraw and Wong (1996) as the fraction of the variance of chosen graphic property between subjects to the total variance, which is the summation variance of between and within subjects of that property.

$ICC=\frac{{\sigma }_{between}^2}{{\sigma }_{between}^2+{\sigma }_{within}^2}$

If the measurements of repeated scans are consistent for each subject, the ICC would be close to one. An ICC value above 0.75 is considered excellent, and one ranging from 0.6 to 0.75 is considered good (Bennett and Miller, 2010). ICC was computed using the “Reliability Analysis” function of SPSS.

Results

Before the formation of networks, nine of the morphometric features were firstly measured of its robustness by SPSS. The mean value of the entire brain was used to realize the ICC measurement. The number of vertices (ICC = 0.95 with p = 3.1 × 10⁻⁴), mean curvature (ICC = 0.80 with p = 2.4 × 10⁻²), surface area (ICC = 0.99 with p = 7.3 × 10⁻⁶), volume (ICC = 0.96 with p = 1.6 × 10⁻⁴) and standard deviation of thickness (ICC = 0.94 with p = 6.3 × 10⁻⁴) are significantly reliable. Robustness was also found in Gauss curvature (ICC = 0.68) and mean thickness (ICC = 0.73), but the results are less significant with p-value of 7.6 × 10⁻² and 5.3 × 10⁻², separately. However, the index of fold and curvature are failed to present acceptable reproducibility (ICC = 0.73 with p = 0.76 and ICC = 0.36 with p = 0.29, respectively). Therefore, only seven of morphometric features were employed in the network construction (Figure 1E, number of vertices, mean and Gauss curvature, surface area, volume, mean and standard deviation of thickness). The subsequent network analyses included the following: (1) assessment of small-world configurations (C_p, L_p, γ, λ, and σ) of the networks, (2) investigation of BC and spatial distribution of the networks' hubs, and (3) estimation of the method's reproducibility. For a rational demonstration, the sparsity applied ranged from 20 to 40% with a step size of 1% to extract all of the above-mentioned measurements for each individual-based brain network.

Small-World Configurations

Initially, apart from the subjects' respective brain networks, C_p and L_p were also derived from the corresponding random networks at each sparsity for each subject. Hence, the comparisons were conducted between the different subjects' brain networks and random networks (Figures 2A,B) or sparsities (Figures 2C–E). The comparison shows that C_p and L_p are profoundly higher than the random ones all over the sparsity range with the maximum p = 1.14 × 10⁻¹¹ (t-values ranged from 11.75 to 27.77 with FDR correction) and p = 6.15 × 10⁻³ (t-values ranged from 2.98 to 5.00 with FDR correction), as revealed using the independent two-sample t-test (Figures 2A,B, respectively). Likewise, γ was larger than one (max = 1.99, min = 1.37) throughout the whole sparsity range, while λ was close to one (max = 1.28, min = 1.02). Hence, as expected σ are found larger than one throughout the entire sparsity range (max = 1.56, min = 1.35) which confirms the existence of small-worldness. Overall, the individual networks exhibit significantly higher C_p than the random network, while remaining similar L_p. Thus, small-worldness is confirmed for each subject at every sparsity. Moreover, the increase and decrease of small-world configurations, as sparsity increased, is illustrated in Figure 2C. Such small-world configuration fluctuations are in accordance with the variation tendency documented in previous reports (He et al., 2007; Kong et al., 2015).

FIGURE 2

Figure 2. Small-world configurations of the individual-based morphological networks. (A,B) Represent the average C_p and L_p of sparsities for each subject. The error bar indicates the deviation caused by different sparsities. (C–F) Show the average C_p, L_p, γ, λ, and σ of the subjects for each sparsity (from 20 to 40% with a step size of 1%). The error bar indicates the deviation caused by different subjects.

Furthermore, the sparsity of 23% is highlighted for convenient comparison with previous studies, including those based on individual-based morphological brain networks (Tijms et al., 2012; Kong et al., 2015), those involving population-based morphological brain networks with various regional descriptors (He et al., 2007; Sanabria-Diaz et al., 2010; Zhu et al., 2012), and functional brain network studies (Zhang et al., 2011). As listed in Table 1, the population-based morphological networks and functional network exhibit smaller results than the current findings, whereas the individual-based morphological network studies present relatively similar results, particularly that of Kong et al. (2015). This finding suggests that the individual morphological brain networks may demonstrate a stronger integration and segregation because the inter-individual variability is highly reserved (Kanai and Rees, 2011).

TABLE 1

Table 1. Comparison of small world configurations between the present study and previous studies.

Spatial Distribution of Hubs

In addition to the assessment of small-world configurations, the spatial distribution of hubs was investigated for different subjects and sparsities. First, the effect of sparsities on hub identification was computed. The BC was averaged through all the subjects at each sparsity to investigate the extent of the sparsity range for each hub region (msjBC). A total of ten regions identified as hubs achieved a sparsity range of over 20% (Figure 3), including the bilateral entorhinal cortex (left 100.00%, right 85.71%), superior temporal gyrus (left 100.00%, right 100.00%), lateral occipital gyrus (left 95.24%, right 100.00%), frontal pole (left 71.43%, right 28.57%) and left caudal middle frontal gyrus (100.00%), left rostral anterior cingulate cortex (71.43%), left isthmus cingulate cortex (33.33%), left parahippocampal gyrus (33.33%). In addition, a total of 4 regions were regarded as hubs in a sparsity range of 4.76 to 14.29%, and the rest of the 52 regions were never observed presenting a hub at any sparsity. The specific extent of sparsity range for each node identified as a hub may be found in Supplementary Table 2.

FIGURE 3

Figure 3. Hubs with a sparsity range exceeding 20%, based on the average of subjects. L, left; R, right; A, anterior; P, posterior; CMF, caudal middle frontal gyrus; ENT, entorhinal cortex; FP, frontal pole; ISTC, isthmus cingulate; LOCC, lateral occipital gyrus; RAC, rostral anterior cingulate; PHG, parahippocampal; ST, superior temporal gyrus. The size of the node represents the proportion of sparsity; the largest node denotes 100% of the proportion. The different color of nodes denotes the different lobes: red-frontal lobe; green-temporal lobe; blue-occipital lobe; yellow-cingulate.

Thereafter, the influence of individual differences on hub determination was investigated. Likewise, the BC was averaged throughout the sparsity range for each subject to determine the number of subjects presenting each region as a hub (mspBC). A total of 14 nodes are hubs for over 20% of the subjects (Figure 4), including the superior temporal gyrus (left 54.55%, right 43.64%), caudal middle frontal (left 38.18%, right 27.27%), entorhinal cortex (left 38.18%, right 25.45%), lateral occipital gyrus (left 32.73%, right 38.18%), parahippocampal gyrus (left 23.64%, right 29.09%), and left isthmus cingulate cortex (27.27%), left rostral anterior cingulate cortex (23.64%), left frontal pole (23.64%), left inferior temporal gyrus (21.82%). A total of 53 regions are considered hubs in 1.82 to 18.18% of the subjects, and just one region comprised no subjects. The specific number of subjects presenting each node as a hub is found in Supplementary Table 3.

FIGURE 4

Figure 4. Hubs in over 20% of subjects based on the average of sparsities. L, left; R, right; A, anterior; P, posterior; CMF, caudal middle frontal gyrus; ENT, entorhinal cortex; FP, frontal pole; ISTC, isthmus cingulate; IT, inferior temporal gyrus; LOCC, lateral occipital gyrus; RAC, rostral anterior cingulate; PHG, parahippocampal; ST, superior temporal gyrus. The size of the node represents the proportion of subjects; the largest node denotes 54.55% of the proportion. The different color of nodes denotes the different lobes: red-frontal lobe; green-temporal lobe; blue-occipital lobe; yellow-cingulate.

The similarity of BC across subjects and sparsities was evaluated through the proximity matrix generated by SPSS (Figure 5). As illustrated in the figure (yellow: similar, blue: dissimilar), the inter-subject and inter-sparsity similarities of BC were observed, and the highly comparable BC was found in adjacent sparsities (close to diagonal) (Figure 5A). The individual differences in subjects significantly influence the BC value (Figure 5B). Clearly, the similarities were significantly high across sparsities. These findings are also in agreement with the above-mentioned results, which state that the hubs can be mostly retained for different sparsities, while hub identification varied profoundly for different subjects.

FIGURE 5

Figure 5. Similarity across sparsities and subjects. (A) Represents the similarity between different sparsities from 20 to 40% with the step size of 1%. (B) Exhibits the similarity between different subjects. The color bar shows the extent of similarity with blue as dissimilar and yellow as similar.

Furthermore, Figure 6 shows the visualization of hub distribution based on the average BC of all the subjects and the entire sparsity range. The results indicate that ten hubs were identified, comprising six heteromodal or unimodal associative regions and four paralimbic regions (Table 2). All these hub regions have been reported in at least one previous study.

FIGURE 6

Figure 6. Spatial distribution of hubs identified based on the average of sparsities and subjects. L, left; R, right; A, anterior; P, posterior; CMF, caudal middle frontal gyrus; ENT, entorhinal cortex; FP, frontal pole; ISTC, isthmus cingulate; LOCC, lateral occipital gyrus; RAC, rostral anterior cingulate; ST, superior temporal gyrus. The size of the node represents the value of BC; the largest node denotes 107.22 of BC. The different color of nodes denotes the different lobes: red-frontal lobe; green-temporal lobe; blue-occipital lobe; yellow-cingulate.

TABLE 2

Table 2. The hubs distribution based on the average BC of all the subjects and the entire sparsity range.

Reproducibility

The method's reproducibility was evaluated by measuring the ICC index of network properties for scans with acquisitions of two different time points in the same subjects. The ICC was investigated throughout the entire sparsity range. Only C_p, L_p, and BC were examined in this study because of the instability of random network generation processing.

The results indicated that C_p was highly reproducible (minimum ICC = 0.71, average ICC = 0.83), as plotted in Figure 7A. Moreover, the robustness of L_p (minimum ICC = 0.63, average ICC = 0.81) was fairly stable (Figure 7B), as the similarities to the BC (minimum ICC = 0.629, average ICC = 0.78) are shown in Figure 7C. Most of the results are significant, except for L_p and BC at sparsity of 40% (p = 5.6 × 10⁻² and p = 5.7 × 10⁻², separately). Hence, the reliability of the present method is well-documented in accordance with the standard proposed by Bennett and Miller (2010). The specific value of ICC for the three measurements across the entire sparsity range is shown in Supplementary Table 4.

FIGURE 7

Figure 7. Reproducibility of the method. (A–C) Represents the ICC of C_p, L_p, and BC as a function of sparsity, respectively. The hollow dot indicates the insignificant result.

Discussions

To our knowledge, this study is the first to collectively employ multiple morphometric features to evaluate the individual morphological brain network. First, small-worldness is confirmed at each sparsity (20%–40% with a step size of 1%) for all the subjects by investigating C_p, L_p, γ, λ, and σ (Figure 2). In addition, the small-world configurations are found to be compatible with previous reports at the sparsity of 23% (Table 1). Thereafter, subject diversity was assessed and was found to significantly influence hub identification, whereas the influence of variety between sparsities was shown to be minimal (Figures 3–5). In addition, the hubs determined by the average of the entire sparsity range and population were found to be highly consistent with previous findings (Figure 6 and Table 2). Finally, the method's reproducibility was ascertained by verifying the ICC of C_p, L_p, and BC throughout the sparsity range (Figure 7). Overall, the results suggest that the proposed method may offer a new way to build the morphological brain network individually with multiple morphometric features.

Interregional Morphological Similarity

The interregional morphological similarity has been observed and verified repeatedly in recent studies based on such morphometric features as cortical thickness and regional volume (Lerch et al., 2006; He et al., 2007; Bassett et al., 2008; Sanabria-Diaz et al., 2010). The interpretation of this observation implies that the brain regions exhibit covaring morphometric features that can form the connected structure, which may capture long-term neurobiological effects (Mechelli et al., 2005). However, the underlying mechanism of this covariance pattern remains elusive. Some conjectures have been debated in the scientific literature, including mutually trophic effects (Ferrer et al., 1995; Aid et al., 2007), environment-related plasticity (Maguire et al., 2000; Draganski et al., 2004; Mechelli et al., 2004), genetic influence (Schmitt et al., 2008), and normal development (Raz et al., 2005; Chen et al., 2011). The axonal tension theory has also been raised recently (Tijms et al., 2012; Kong et al., 2015), stating that the interconnected areas are becoming either thicker or thinner as a result of being pulled by a mechanical force (Van Essen, 1997).

Moreover, the relationships among morphological (M), functional (F), and anatomical (A) connectivities were investigated in previous studies. For instance, Reid et al. (2017) have demonstrated that the extent of accordance between F-M modalities varied remarkably across seed regions. Wang et al. (2015) have illustrated a tight coupling of F-M modalities in connectivity strength and network topologic organizations (i.e., modules, rich club, and motifs), further indicating that the neuropsychiatric disorders may considerably break this coupling. In addition, Gong et al. (2012) have found an A-M convergence in 35–40% of morphological connectivities, with the majority of the convergences observed in the positive morphological connectivities.

Further research in neurobiology would significantly assist the construction and exploration of morphological brain networks. In addition, the data on different modalities were typically used or analyzed separately because they were supposed to capture the distinct underlying phenomena. However, the correspondences in topological structures have been recently observed across modalities, suggesting that the combination of multimodal research will become increasingly important in neuroscience.

Spatial Distribution of Hubs

The spatial distribution of hubs (nodes with higher BC of average subjects and the entire sparsity range) in the present study was found to be strikingly similar to the results of functional studies. The consistent regions are the bilateral superior temporal gyrus, entorhinal cortex, and lateral occipital gyrus, and left isthmus cingulate, frontal pole (Table 2). Such accordance is quite thought-provoking, although the different criteria were employed to identify the hubs (BC, degree, and L_p). Hence, the findings may suggest that the individual-based morphological brain networks may exhibit an improved consistency with functional studies.

In addition, the hub spatial distribution of different subjects (averaged by sparsities) shows striking diversity in the present study (Supplementary Table 3). As illustrated in Figures 4, 5 (right), inter-individual variability is implied as a major influence on hub determination. Indeed, a number of studies have demonstrated unmistakable discrepancies in brain structures across individuals (Kanai and Rees, 2011). For example, genetic differences can lead to brain morphological changes (Thompson et al., 2001; Pol et al., 2006), as well as the variance in occupations (Maguire et al., 2000; Gaser and Schlaug, 2003) and skills (Mechelli et al., 2004). Moreover, cognition is found to be correlated with cerebral structures. For instance, individual differences in degrees of empathy were observed to be associated with gray matter density (Eres et al., 2015) and volume (Banissy et al., 2012). Psychiatric diseases such as epilepsy also involve unpredictable foci locations and numbers (Engel, 2006). Therefore, individual-based investigations are indispensable in clinical and scientific research.

Hence, investigation of the organization of morphological brain networks, along with functional and anatomical network organization on a case-by-case basis is both intriguing and indispensable. The individual-based accordance across modalities may also be an indicator in future cognitive and psychiatric research.

Methodology Issues and Future Research

Several methodological issues are found in the present study, which should be addressed and solved in future research.

First, the interregional connections were measured as the Pearson correlation coefficient in the present study, as in many earlier studies (He et al., 2007; Chen et al., 2011; Shi et al., 2013). However, this measurement shows that the observed connections between any two regions are actually the summation of the direct and indirect correlations. A widely used approach for eliminating the influence of other regions is to alternatively use partial correlation (Bassett et al., 2008). The number of variables should be less than a number of samples of each variable in partial correlation computation. However, in the present study, a total of 68 regions (i.e., variables) with nine morphometric features (i.e., samples) are used, which goes against the regular partial correlation. Consequently, in future research, the partial correlation computation must be adjusted to increase the accuracy of results addressing interregional connections.

Second, the networks were binarized with a sparsity range from 20 to 40% for topological analyses, as in previous studies (He et al., 2007; Kong et al., 2015). This binarization leads to the exclusion of 60–80% of connection information simply because such connections were deemed “not important” based on the correlation algorithms. The underlying mechanism of morphological connections remains unclear (see interregional morphological similarity). Therefore, whether or not the “not important” connections are really useless has not been confirmed (Barrat et al., 2004). As such, further exploration of brain networks as the full-connected weighted graph is essential to obtaining additional insights into the human brain.

Third, the nine morphometric features applied in the present study were far from painstakingly selected. However, the anatomical information including cortical thickness, gray matter volume, surface area, and curvature is varied (He et al., 2007; Sanabria-Diaz et al., 2010). In addition, redundant information may have existed in the three curvature-related measurements (mean curvature, Gauss curvature, and curvature index). Therefore, an optimal selection of morphometric features will be further explored.

Fourth, the present study only investigated the individual-based morphological connectivities, whereas the functional and anatomical connections were found to demonstrate a mutual relationship with the morphological connections in the present and previous studies (Sui et al., 2012; Mueller et al., 2013). Hence, in future research, the multimodal images should be employed in individual-based brain network analyses.

Conclusion

In this study, a new idea is proposed to construct an individual-based morphological brain network. To our knowledge, this study is the first to collectively use multiple morphometric features to form interregional connections. The results of network topological analyses have demonstrated this method's feasibility, and the verification of reproducibility has supported its excellent robustness. Our findings on hubs' spatial distribution have provided profound indications of individual differences that cannot and should not be overlooked. Moreover, the hubs obtained by averaging across subjects and sparsities have been shown to be largely compatible with an individual-based functional study, leading us to investigate whether individual-based multimodal brain networks share further similarities. The proposed method may offer a novel approach in investigating the cerebral organization individually. The interrelationship between modalities, combined with functional connectivities, is worthy of further exploration.

Ethics Statement

This study was carried out with the data of human subjects from the Open Access Series of Imaging Studies database, which is made available by the Washington University Alzheimer's Disease Research Center, Dr. Randy Buckner at the Howard Hughes Medical Institute (HHMI) at Harvard University, the Neuroinformatics Research Group (NRG) at Washington University School of Medicine, and the Biomedical Informatics Research Network (BIRN).

Author Contributions

The design of the work was proposed by WL and CY. The data processing was completed by WL, YN, and XZ. Drafting of the work was completed by WL and CY. Revision of the draft was conducted under the tutelage of FS. SW and QW provided the equipment and medical suggestions for the work. All the authors offered final approval for the version to be published. All the authors agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Funding

This research was partly supported by Beijing Nova Program (xx2016120), National Natural Science Foundation of China (81101107, 31640035 and 71661167001), Natural Science Foundation of Beijing Municipality (4162008), and Beijing Municipal Education Commission (PXM2017_014204_500012).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnana.2017.00034/full#supplementary-material

References

Achard, S., Salvador, R., Whitcher, B., Suckling, J., and Bullmore, E. (2006). A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J. Neurosci. 26, 63–72. doi: 10.1523/JNEUROSCI.3874-05.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Aid, T., Kazantseva, A., Piirsoo, M., Palm, K., and Timmusk, T. (2007). Mouse and rat BDNF gene structure and expression revisited. J. Neurosci. Res. 85, 525–535. doi: 10.1002/jnr.21139

PubMed Abstract | CrossRef Full Text | Google Scholar

Alexander, G. E., Bergfield, K. L., Chen, K., Reiman, E. M., Hanson, K. D., Lin, L., et al. (2012). Gray matter network associated with risk for Alzheimer's disease in young to middle-aged adults. Neurobiol. Aging 33, 2723–2732. doi: 10.1016/j.neurobiolaging.2012.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Banissy, M. J., Kanai, R., Walsh, V., and Rees, G. (2012). Inter-individual differences in empathy are reflected in human brain structure. Neuroimage 62, 2034–2039. doi: 10.1016/j.neuroimage.2012.05.081

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrat, A., Barthélemy, M., Pastor-Satorras, R., and Vespignani, A. (2004). The architecture of complex weighted networks. Proc. Natl. Acad. Sci. U.S.A. 101, 3747–3752. doi: 10.1073/pnas.0400087101

PubMed Abstract | CrossRef Full Text | Google Scholar

Bassett, D. S., and Bullmore, E. (2006). Small-world brain networks. Neuroscientist 12, 512–523. doi: 10.1177/1073858406293182

PubMed Abstract | CrossRef Full Text | Google Scholar

Bassett, D. S., Bullmore, E., Verchinski, B. A., Mattay, V. S., Weinberger, D. R., and Meyer-Lindenberg, A. (2008). Hierarchical organization of human cortical networks in health and schizophrenia. J. Neurosci. 28, 9239–9248. doi: 10.1523/JNEUROSCI.1929-08.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Batalle, D., Muñoz-Moreno, E., Figueras, F., Bargallo, N., Eixarch, E., and Gratacos, E. (2013). Normalization of similarity-based individual brain networks from gray matter MRI and its association with neurodevelopment in infants with intrauterine growth restriction. Neuroimage 83, 901–911. doi: 10.1016/j.neuroimage.2013.07.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Bennett, C. M., and Miller, M. B. (2010). How reliable are the results from functional magnetic resonance imaging? Ann. N.Y. Acad. Sci. 1191, 133–155. doi: 10.1111/j.1749-6632.2010.05446.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L., Sepulcre, J., Talukdar, T., Krienen, F. M., Liu, H., Hedden, T., et al. (2009). Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease. J. Neurosci. 29, 1860–1873. doi: 10.1523/JNEUROSCI.5062-08.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z. J., He, Y., Rosa-Neto, P., Gong, G., and Evans, A. C. (2011). Age-related alterations in the modular organization of structural cortical network by using cortical thickness from MRI. Neuroimage 56, 235–245. doi: 10.1016/j.neuroimage.2011.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Dale, A. M., Fischl, B., and Sereno, M. I. (1999). Cortical surface-based analysis: I. Segmentation and surface reconstruction. Neuroimage 9, 179–194. doi: 10.1006/nimg.1998.0395

PubMed Abstract | CrossRef Full Text | Google Scholar

Damoiseaux, J. S., Beckmann, C. F., Arigita, E. J., Barkhof, F., Scheltens, P., Stam, C. J., et al. (2008). Reduced resting-state brain activity in the “default network” in normal aging. Cereb. Cortex 18, 1856–1864. doi: 10.1093/cercor/bhm207

PubMed Abstract | CrossRef Full Text | Google Scholar

Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Draganski, B., Gaser, C., Busch, V., Schuierer, G., Bogdahn, U., and May, A. (2004). Neuroplasticity: changes in grey matter induced by training. Nature 427, 311–312. doi: 10.1038/427311a

PubMed Abstract | CrossRef Full Text | Google Scholar

Eguiluz, V. M., Chialvo, D. R., Cecchi, G. A., Baliki, M., and Apkarian, A. V. (2005). Scale-free brain functional networks. Phys. Rev. Lett. 94:018102. doi: 10.1103/PhysRevLett.94.018102

PubMed Abstract | CrossRef Full Text | Google Scholar

Engel, J. (2006). ILAE classification of epilepsy syndromes. Epilepsy Res. 70, 5–10. doi: 10.1016/j.eplepsyres.2005.11.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Eres, R., Decety, J., Louis, W. R., and Molenberghs, P. (2015). Individual differences in local gray matter density are associated with differences in affective and cognitive empathy. Neuroimage 117, 305–310. doi: 10.1016/j.neuroimage.2015.05.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, A. C. (2013). Networks of anatomical covariance. Neuroimage 80, 489–504. doi: 10.1016/j.neuroimage.2013.05.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrer, I., Blanco, R., Carulla, M., Condom, M., Alcántara, S., Olivé, M., et al. (1995). Transforming growth factor-α immunoreactivity in the developing and adult brain. Neuroscience 66, 189–199. doi: 10.1016/0306-4522(94)00584-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Fischl, B., Sereno, M. I., and Dale, A. M. (1999). Cortical surface-based analysis: II: inflation, flattening, and a surface-based coordinate system. Neuroimage 9, 195–207. doi: 10.1006/nimg.1998.0396

PubMed Abstract | CrossRef Full Text | Google Scholar

Freeman, L. C. (1977). A set of measures of centrality based on betweenness. Sociometry 40, 35–41. doi: 10.2307/3033543

CrossRef Full Text | Google Scholar

Gaser, C., and Schlaug, G. (2003). Brain structures differ between musicians and non-musicians. J. Neurosci. 23, 9240–9245.

PubMed Abstract | Google Scholar

Genovese, C. R., Lazar, N. A., and Nichols, T. (2002). Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage 15, 870–878. doi: 10.1006/nimg.2001.1037

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, G., He, Y., Chen, Z. J., and Evans, A. C. (2012). Convergence and divergence of thickness correlations with diffusion connections across the human cerebral cortex. Neuroimage 59, 1239–1248. doi: 10.1016/j.neuroimage.2011.08.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, G., He, Y., Concha, L., Lebel, C., Gross, D. W., Evans, A. C., et al. (2009a). Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography. Cereb. Cortex 19, 524–536. doi: 10.1093/cercor/bhn102

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, G., Rosa-Neto, P., Carbonell, F., Chen, Z. J., He, Y., and Evans, A. C. (2009b). Age-and gender-related differences in the cortical anatomical network. J. Neurosci. 29, 15684–15693. doi: 10.1523/JNEUROSCI.2308-09.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Hagmann, P., Cammoun, L., Gigandet, X., Meuli, R., Honey, C. J., Wedeen, V. J., et al. (2008). Mapping the structural core of human cerebral cortex. PLoS Biol. 6:e159. doi: 10.1371/journal.pbio.0060159

PubMed Abstract | CrossRef Full Text | Google Scholar

Hagmann, P., Kurant, M., Gigandet, X., Thiran, P., Wedeen, V. J., Meuli, R., et al. (2007). Mapping human whole-brain structural networks with diffusion MRI. PLoS ONE 2:e597. doi: 10.1371/journal.pone.0000597

PubMed Abstract | CrossRef Full Text | Google Scholar

Harriger, L., Van Den Heuvel, M. P., and Sporns, O. (2012). Rich club organization of macaque cerebral cortex and its role in network communication. PLoS ONE 7:e46497. doi: 10.1371/journal.pone.0046497

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Y., Chen, Z. J., and Evans, A. C. (2007). Small-world anatomical networks in the human brain revealed by cortical thickness from MRI. Cereb. Cortex 17, 2407–2419. doi: 10.1093/cercor/bhl149

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Y., Dagher, A., Chen, Z., Charil, A., Zijdenbos, A., Worsley, K., et al. (2009). Impaired small-world efficiency in structural cortical networks in multiple sclerosis associated with white matter lesion load. Brain 132, 3366–3379. doi: 10.1093/brain/awp089

PubMed Abstract | CrossRef Full Text | Google Scholar

Honey, C. J., Kötter, R., Breakspear, M., and Sporns, O. (2007). Network structure of cerebral cortex shapes functional connectivity on multiple time scales. Proc. Natl. Acad. Sci. U.S.A. 104, 10240–10245. doi: 10.1073/pnas.0701519104

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanai, R., and Rees, G. (2011). The structural basis of inter-individual differences in human behaviour and cognition. Nat. Rev. Neurosci. 12, 231–242. doi: 10.1038/nrn3000

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, A. R., Wang, L., and Beg, M. F. (2008). FreeSurfer-initiated fully-automated subcortical brain segmentation in MRI using large deformation diffeomorphic metric mapping. Neuroimage 41, 735–746. doi: 10.1016/j.neuroimage.2008.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, X. Z., Liu, Z., Huang, L., Wang, X., Yang, Z., Zhou, G., et al. (2015). Mapping individual brain networks using statistical similarity in regional morphology from MRI. PLoS ONE 10:e0141840. doi: 10.1371/journal.pone.0141840

PubMed Abstract | CrossRef Full Text | Google Scholar

Latora, V., and Marchiori, M. (2001). Efficient behavior of small-world networks. Phys. Rev. Lett. 87:198701. doi: 10.1103/PhysRevLett.87.198701

PubMed Abstract | CrossRef Full Text | Google Scholar

Lerch, J. P., Worsley, K., Shaw, W. P., Greenstein, D. K., Lenroot, R. K., Giedd, J., et al. (2006). Mapping anatomical correlations across cerebral cortex (MACACC) using cortical thickness from MRI. Neuroimage 31, 993–1003. doi: 10.1016/j.neuroimage.2006.01.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Yuan, X., Pu, F., Li, D., Fan, Y., Wu, L., et al. (2014). Abnormal changes of multidimensional surface features using multivariate pattern classification in amnestic mild cognitive impairment patients. J. Neurosci. 34, 10541–10553. doi: 10.1523/JNEUROSCI.4356-13.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Liang, M., Zhou, Y., He, Y., Hao, Y., Song, M., et al. (2008). Disrupted small-world networks in schizophrenia. Brain 131, 945–961. doi: 10.1093/brain/awn018

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo, C. Y., Wang, P. N., Chou, K. H., Wang, J., He, Y., and Lin, C. P. (2010). Diffusion tensor tractography reveals abnormal topological organization in structural cortical networks in Alzheimer's disease. J. Neurosci. 30, 16876–16885. doi: 10.1523/JNEUROSCI.4136-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Maguire, E. A., Gadian, D. G., Johnsrude, I. S., Good, C. D., Ashburner, J., Frackowiak, R. S., et al. (2000). Navigation-related structural change in the hippocampi of taxi drivers. Proc. Natl. Acad. Sci. U.S.A. 97, 4398–4403. doi: 10.1073/pnas.070039597

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcus, D. S., Wang, T. H., Parker, J., Csernansky, J. G., Morris, J. C., and Buckner, R. L. (2007). Open Access Series of Imaging Studies (OASIS): cross-sectional MRI data in young, middle aged, nondemented, and demented older adults. J. Cogn. Neurosci. 19, 1498–1507. doi: 10.1162/jocn.2007.19.9.1498

PubMed Abstract | CrossRef Full Text | Google Scholar

Mcgraw, K. O., and Wong, S. P. (1996). Forming inferences about some intraclass correlation coefficients. Psychol. Methods 1:30. doi: 10.1037/1082-989x.1.1.30

CrossRef Full Text | Google Scholar

Mechelli, A., Crinion, J. T., Noppeney, U., O'doherty, J., Ashburner, J., Frackowiak, R. S., et al. (2004). Neurolinguistics: structural plasticity in the bilingual brain. Nature 431, 757–757. doi: 10.1038/431757a

PubMed Abstract | CrossRef Full Text | Google Scholar

Mechelli, A., Friston, K. J., Frackowiak, R. S., and Price, C. J. (2005). Structural covariance in the human cortex. J. Neurosci. 25, 8303–8310. doi: 10.1523/JNEUROSCI.0357-05.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mesulam, M. M. (1998). From sensation to cognition. Brain 121, 1013–1052. doi: 10.1093/brain/121.6.1013

PubMed Abstract | CrossRef Full Text | Google Scholar

Micheloyannis, S., Pachou, E., Stam, C. J., Vourkas, M., Erimaki, S., and Tsirka, V. (2006). Using graph theoretical analysis of multi channel EEG to evaluate the neural efficiency hypothesis. Neurosci. Lett. 402, 273–277. doi: 10.1016/j.neulet.2006.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Mueller, S., Keeser, D., Samson, A. C., Kirsch, V., Blautzik, J., Grothe, M., et al. (2013). Convergent findings of altered functional and structural brain connectivity in individuals with high functioning autism: a multimodal MRI study. PLoS ONE 8:e67329. doi: 10.1371/journal.pone.0067329

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulder, E. R., de Jong, R. A., Knol, D. L., Van Schijndel, R. A., Cover, K. S., Visser, P. J., et al. (2014). Hippocampal volume change measurement: quantitative assessment of the reproducibility of expert manual outlining and the automated methods FreeSurfer and FIRST. Neuroimage 92, 169–181. doi: 10.1016/j.neuroimage.2014.01.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Newman, M. E. (2004). Fast algorithm for detecting community structure in networks. Phys. Rev. E 69:066133. doi: 10.1103/physreve.69.066133

PubMed Abstract | CrossRef Full Text | Google Scholar

Panizzon, M. S., Fennema-Notestine, C., Eyler, L. T., Jernigan, T. L., Prom-Wormley, E., Neale, M., et al. (2009). Distinct genetic influences on cortical surface area and cortical thickness. Cereb. Cortex 19, 2728–2735. doi: 10.1093/cercor/bhp026

PubMed Abstract | CrossRef Full Text | Google Scholar

Pienaar, R., Fischl, B., Caviness, V., Makris, N., and Grant, P. E. (2008). A methodology for analyzing curvature in the developing brain from preterm to adult. Int. J. Imaging Syst. Technol. 18, 42–68. doi: 10.1002/ima.20138

PubMed Abstract | CrossRef Full Text | Google Scholar

Pol, H. E. H., Schnack, H. G., Posthuma, D., Mandl, R. C., Baaré, W. F., Van Oel, C., et al. (2006). Genetic contributions to human brain morphology and intelligence. J. Neurosci. 26, 10235–10242. doi: 10.1523/JNEUROSCI.1312-06.2006

CrossRef Full Text | Google Scholar

Raj, A., Mueller, S. G., Young, K., Laxer, K. D., and Weiner, M. (2010). Network-level analysis of cortical thickness of the epileptic brain. Neuroimage 52, 1302–1313. doi: 10.1016/j.neuroimage.2010.05.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Raz, N., Lindenberger, U., Rodrigue, K. M., Kennedy, K. M., Head, D., Williamson, A., et al. (2005). Regional brain changes in aging healthy adults: general trends, individual differences and modifiers. Cereb. Cortex 15, 1676–1689. doi: 10.1093/cercor/bhi044

PubMed Abstract | CrossRef Full Text | Google Scholar

Reid, A. T., Hoffstaedter, F., Gong, G., Laird, A. R., Fox, P., Evans, A. C., et al. (2017). A seed-based cross-modal comparison of brain connectivity measures. Brain Struct. Funct. 222, 1131–1151. doi: 10.1007/s00429-016-1264-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanabria-Diaz, G., Melie-García, L., Iturria-Medina, Y., Alemán-Gómez, Y., Hernández-González, G., Valdés-Urrutia, L., et al. (2010). Surface area and cortical thickness descriptors reveal different attributes of the structural human brain networks. Neuroimage 50, 1497–1510. doi: 10.1016/j.neuroimage.2010.01.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaer, M., Cuadra, M. B., Tamarit, L., Lazeyras, F., Eliez, S., and Thiran, J. P. (2008). A surface-based approach to quantify local cortical gyrification. IEEE Trans. Med. Imaging 27, 161–170. doi: 10.1109/TMI.2007.903576

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitt, J., Lenroot, R., Wallace, G., Ordaz, S., Taylor, K., Kabani, N., et al. (2008). Identification of genetically mediated cortical networks: a multivariate study of pediatric twins and siblings. Cereb. Cortex 18, 1737–1747. doi: 10.1093/cercor/bhm211

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, F., Wang, L., Peng, Z., Wee, C. Y., and Shen, D. (2013). Altered modular organization of structural cortical networks in children with autism. PLoS ONE 8:e63131. doi: 10.1371/journal.pone.0063131

PubMed Abstract | CrossRef Full Text | Google Scholar

Sporns, O. (2013). Network attributes for segregation and integration in the human brain. Curr. Opin. Neurobiol. 23, 162–171. doi: 10.1016/j.conb.2012.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Sporns, O., Chialvo, D. R., Kaiser, M., and Hilgetag, C. C. (2004). Organization, development and function of complex brain networks. Trends Cogn. Sci. 8, 418–425. doi: 10.1016/j.tics.2004.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Stam, C. J. (2004). Functional connectivity patterns of human magnetoencephalographic recordings: a ‘small-world’ network? Neurosci. Lett. 355, 25–28. doi: 10.1016/j.neulet.2003.10.063

PubMed Abstract | CrossRef Full Text | Google Scholar

Stam, C. J., Jones, B. F., Nolte, G., Breakspear, M., and Scheltens, P. (2007). Small-world networks and functional connectivity in Alzheimer's disease. Cereb. Cortex 17, 92–99. doi: 10.1093/cercor/bhj127

PubMed Abstract | CrossRef Full Text | Google Scholar

Sui, J., Adali, T., Yu, Q., Chen, J., and Calhoun, V. D. (2012). A review of multivariate methods for multimodal fusion of brain imaging data. J. Neurosci. Methods 204, 68–81. doi: 10.1016/j.jneumeth.2011.10.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, P. M., Cannon, T. D., Narr, K. L., van Erp, T., Poutanen, V.-P., Huttunen, M., et al. (2001). Genetic influences on brain structure. Nat. Neurosci. 4, 1253–1258. doi: 10.1038/nn758

PubMed Abstract | CrossRef Full Text | Google Scholar

Tijms, B. M., Seriès, P., Willshaw, D. J., and Lawrie, S. M. (2012). Similarity-based extraction of individual networks from gray matter MRI scans. Cereb. Cortex 22, 1530–1541. doi: 10.1093/cercor/bhr221

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Den Heuvel, M. P., Mandl, R. C., Stam, C. J., Kahn, R. S., and Pol, H. E. H. (2010). Aberrant frontal and temporal complex network structure in schizophrenia: a graph theoretical analysis. J. Neurosci. 30, 15915–15926. doi: 10.1523/JNEUROSCI.2874-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Essen, D. C. (1997). A tension-based theory of morphogenesis and compact wiring in the central nervous system. Nature 385:313. doi: 10.1038/385313a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Wang, X., Xia, M., Liao, X., Evans, A., and He, Y. (2015). GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Front. Hum. Neurosci. 9:386. doi: 10.3389/fnhum.2015.00386

CrossRef Full Text | Google Scholar

Wang, L., Metzak, P. D., Honer, W. G., and Woodward, T. S. (2010). Impaired efficiency of functional networks underlying episodic memory-for-context in schizophrenia. J. Neurosci. 30, 13171–13179. doi: 10.1523/JNEUROSCI.3514-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Dai, Z., Gong, G., Zhou, C., and He, Y. (2015). Understanding structural-functional relationships in the human brain a large-scale network perspective. Neuroscientist 21, 290–305. doi: 10.1177/1073858414537560

PubMed Abstract | CrossRef Full Text | Google Scholar

Watts, D. J., and Strogatz, S. H. (1998). Collective dynamics of ‘small-world’ networks. Nature 393, 440–442. doi: 10.1038/30918

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, K., Taki, Y., Sato, K., Kinomura, S., Goto, R., Okada, K., et al. (2012). Age-related changes in topological organization of structural brain networks in healthy individuals. Hum. Brain Mapp. 33, 552–568. doi: 10.1002/hbm.21232

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, M., Wang, J., and He, Y. (2013). BrainNet Viewer: a network visualization tool for human brain connectomics. PLoS ONE 8:e68910. doi: 10.1371/journal.pone.0068910

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C., Gong, G., Wang, J., Wang, D., Liu, D., Zhu, C., et al. (2011). Sex-and brain size-related small-world structural cortical networks in young adults: a DTI tractography study. Cereb. Cortex 21, 449–458. doi: 10.1093/cercor/bhq111

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Wang, J., Wu, Q., Kuang, W., Huang, X., He, Y., et al. (2011). Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder. Biol. Psychiatry 70, 334–342. doi: 10.1016/j.biopsych.2011.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Wang, Y., Li, Y., Yap, P. T., Shen, D., and Initiative, A. S. D. N. (2011). Hierarchical anatomical brain networks for MCI prediction: revisiting volumetric measures. PLoS ONE 6:e21935. doi: 10.1371/journal.pone.0021935

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, W., Wen, W., He, Y., Xia, A., Anstey, K. J., and Sachdev, P. (2012). Changing topological patterns in normal aging using large-scale structural networks. Neurobiol. Aging 33, 899–913. doi: 10.1016/j.neurobiolaging.2010.06.022

PubMed Abstract | CrossRef Full Text | Google Scholar