Research Article  Open Access
YongHo Choi, JeYeon Yun, BoHyun Kim, MinHo Lee, SaKwang Song, JongMin Lee, "GenderRelated and Hemispheric Effects in Cortical ThicknessBased Hemispheric Brain Morphological Network", BioMed Research International, vol. 2020, Article ID 3560259, 13 pages, 2020. https://doi.org/10.1155/2020/3560259
GenderRelated and Hemispheric Effects in Cortical ThicknessBased Hemispheric Brain Morphological Network
Abstract
Objective. The current study examined genderrelated differences in hemispheric asymmetries of graph metrics, calculated from a cortical thicknessbased brain structural covariance network named hemispheric morphological network. Methods. Using the T1weighted magnetic resonance imaging scans of 285 participants (150 females, 135 males) retrieved from the Human Connectome Project (HCP), hemispheric morphological networks were constructed per participant. In these hemispheric morphologic networks, the degree of similarity between two different brain regions in terms of the distributed patterns of cortical thickness values (the Jensen–Shannon divergence) was defined as weight of network edge that connects two different brain regions. After the calculation and summation of global and local graph metrics (across the network sparsity levels ), asymmetry indexes of these graph metrics were derived. Results. Hemispheric morphological networks satisfied smallworldness and global efficiency for the network sparsity ranges of . Betweengroup comparisons (female versus male) of asymmetry indexes revealed opposite directionality of asymmetries (leftward versus rightward) for global metrics of normalized clustering coefficient, normalized characteristic path length, and global efficiency (all ). For the local graph metrics, larger rightward asymmetries of cingulatesuperior parietal gyri for nodal efficiency in male compared to female, larger leftward asymmetry of temporal pole for degree centrality in female compared to male, and opposite directionality of interhemispheric asymmetry of rectal gyrus for degree centrality between female (rightward) and male (leftward) were shown (all ). Conclusion. Patterns of interhemispheric asymmetries for cingulate, superior parietal gyrus, temporal pole, and rectal gyrus are different between male and female for the similarities of the cortical thickness distribution with other brain regions. Accordingly, possible effect of genderbyhemispheric interaction has to be considered in future studies of brain morphology and brain structural covariance networks.
1. Introduction
Structural MRI (sMRI) has been used in attempts to construct group brain networks by detecting wholebrain morphological connectivity patterns based on the interregional morphological similarities or “brain structural covariance” across participants (interindividual) [1] or per individual (intraindividual) [2]. On the other hand, diffusion tensor imaging (DTI) and functional magnetic resonance imaging MRI (fMRI) have been commonly used for constructing the individual brain networks connected by way of the axonal pathways (in DTI) or coupled changes of brain functional activations (in fMRI), respectively [3, 4]. The brain structural covariance network partly reflects the patterns of brain white matterbased physical connectivity [5–7], coordinated oscillations in BOLD signal changes across the whole brain during the resting status (socalled resting state functional connectivity) [8, 9], and coordinated brain development and maturation [1, 10, 11]. For the construction of brain network, the sMRI could be preferred to fMRI or DTI because of its advantages in terms of easy access, high signaltonoise ratio, and relative insensitivity to artifacts.
The hemispheric asymmetry of the structural and functional networks has been studied because it could be an important aspect in understanding the organization of the human brain [12–14]. Interregional asymmetry studies were focused on identifying which regions and/or connections are stronger in one hemisphere than in the other. Further, representing the characteristics of the brain as a network can provide novel insights into the brain system [15]. For instance, a network of hemispheric brain asymmetry could assess the differences in the configuration (degrees of integration as well as segregation) of anatomical substrates for facilitating the information transfer between different brain regions [12]. A graphtheoretical analysis of brain networks has been applied successfully in investigating how the human brain is organized in healthy individuals and in patients with neuropsychological diseases such as Alzheimer’s disease (AD), schizophrenia [16–19], major depressive disorder [20], obsessivecompulsive disorder [2], and eating disorder [21], among others. Tian et al. investigated the difference in the topologies of the hemispheric functional networks in healthy righthanded adults. The researchers showed regions that exhibited hemisphericrelated differences in regions formerly observed to be functionally or structurally asymmetric [13]. Zhong et al. investigated how topological asymmetries evolve from adolescence to young adulthood and determined that rightward asymmetry in both global and local network efficiencies was consistently observed in adolescents and young adults and that the degree of asymmetry was significantly decreased in young adults [22].
Of note, gender is known as one of the key factors associated with not only the interregional asymmetry but also the hemispheric asymmetry of the brain networks. For instance, IturriaMedina et al. studied the difference between the topological organizations of the two hemispheric structural networks in healthy righthanded individuals and discovered that the left hemisphere presented more central and essential regions, whereas the right hemisphere was more efficient and interconnected [12]. Tian et al. showed that males tend to be more locally efficient in their right hemispheric networks and females in their left hemispheric networks, which suggested that the local efficiency of the hemispheric network could be associated with behavioral and cognitive differences between men and women [13]. Caeyenberghs and Leemans found that males have a greater global efficiency of the structural hemispheric networks than do females [23].
Compared to the interindividual brain structural covariance network, intraindividual brain covariance networks could be more convenient in finding the neural correlates of specific neuropsychiatric disorders (to be used for the machine learningbased disease classification) and of specific human behaviors (by way of the calculation of correlation coefficients). Accordingly, several methods have been developed to estimate the intraindividual version of brain structural covariance network [24–26]. Tijms et al. constructed an individual network using the intensity similarity between the patches on volume space, with different numbers of nodes for each individual. Thus, it was difficult to be used for group analysis because the sizes of the constructed individual networks were generally different [24]. Kong et al. proposed to construct a morphological brain network based on interregional morphological connectivity in terms of the Kullback–Leibler (KL) divergence among a regional morphological distribution and revealed longitudinal changes in morphological connectivity of the thalamus after longterm sleep deprivation [25, 27]. However, because the KL divergence is asymmetric, the average of a corresponding pair of values was used to calculate the morphological connectivity between two regions.
To the best of our knowledge, no study has reported hemisphererelated differences in the topological organization of brain structural (morphological) covariance networks. Furthermore, very little is known about the patterns of genderbyhemispheric interaction in terms of the global (i.e., balance of brain network integration) versus segregation and local (i.e., which brain regions well reflect the general trend of brain morphology across the brain) graph metrics. By applying the norm of Jensen–Shannon (JS) divergence for calculating the weights of network edges [28, 29], an intraindividual brain structural covariance network named “hemispheric morphological network”—that reflects the levels of similarities between different brain regions in terms of the cortical thickness distribution per participant—was constructed. We investigated (1) whether each hemispheric morphological network exhibits smallworld and high efficiency properties, (2) whether the hemispheric and gender effects exist in the morphological network at the global and regional scales, and (3) whether a gender difference in the hemispheric asymmetry exists in the morphological network at the global and regional scales.
2. Materials and Methods
The framework of this study design is shown in Figure 1. The details of the procedures are described below. All MRI data were analyzed using CIVET v2.1 pipeline (https://wiki.bic.mni.mcgill.ca/ServicesSoftware/CIVET), MATLAB (MATLAB R2016b, The MathWorks, Inc., Natick, Massachusetts, United States), and R [30].
2.1. Sample Characteristics and MRI Acquisition
The data used in this study were selected from the publicly available S900 Release of Human Connectome Project (HCP), WUMinn Consortium [31, 32]. The samples were selected from the HCP data according to these two criteria: (1) they had a handedness score greater than or equal to 50 and (2) they did not share family members between them. A total of 285 righthanded individuals between the ages of 22 and 36 years (150 females and 135 males) were selected. The experiments were performed in accordance with relevant guidelines and regulations, and the experimental protocol was approved by the Institutional Review Board (IRB) (IRB # 201204036; Title: Mapping the Human Connectome: Structure, Function, and Heritability). Written informed consent was obtained from all participants. The data analysis was performed in accordance with the ethical guidelines of the Hanyang University Ethics Committee. The structural T1weighted MRI was acquired on a 3T Siemens Skyra scanner using a threedimensional magnetization prepared rapid gradientecho (MPRAGE) sequence. The main MR parameters were the repetition time , echo time , inversion time , , field of view , and 0.7 mm isotropic voxels.
2.2. Cortical Surface Modeling and Extraction of Cortical Thickness Values
T1weighted MRIs were preprocessed by PreFreeSurfer using HCP Pipelines anatomical modules [33]. And preprocessed T1weighted MRIs were processed using CIVET v2.1 pipeline developed by the Montreal Neurological Institute (MNI). First, native MRI images were registered into a standardized stereotaxic space using an affine transformation [34]. The intensity nonuniformity resulting from the inhomogeneity in the magnetic field was corrected using N3 algorithms [35], and nonbrain tissue was excluded using the brain extraction tool (BET) algorithm [36]. Subsequently, gray matter (GM), white matter (WM), and cerebrospinal fluid were defined on the stereotaxic brain mask using an artificial neural net classifier [37]. The inner surfaces consisting of polygons of triangular components were extracted from the defined WM mask images and then expanded to the outer surfaces using the constrained Laplacianbased automated segmentation with proximities (CLASP) algorithm [38]. The cortical thickness was defined using the Euclidean distance between the linked vertices of the inner and outer surfaces, after applying an inverse transformation matrix to the cortical surfaces and reconstructing them in the native space [38, 39].
2.3. Construction of the Hemispheric Brain Morphological Network
The cortex was initially divided into 78 regions using the automated anatomical labeling (AAL) atlas [40] (see Table S1 in Supplementary Materials). The regions were then upsampled into 512 cortical regions with a similar number of vertices in each region [41, 42] because recent studies have highlighted the effect of region size on structural connectivity [43, 44]. We used the sphere model corresponding cerebral surface model to balance the region size in the AAL atlas. We applied means algorithm with different number of clusters to each region in order to make all the subregions of similar size. Network nodes were specified as 512 cortical regions, which contained an average of 149 vertices with a standard deviation of 22 vertices. Since we used the AAL atlas, the generated subregions maintain anatomical boundaries such as gyri and sulci. Regions within the upsampled AAL atlas are therefore approximately equal in size and maintaining anatomically constrained. Finally, left and righthemispheric morphological networks were constructed from each subject. Network edges between two regions were defined using the Jensen–Shannon divergence (JSD), as follows:
The parameters and are the distributions of the cortical thickness in each region, which were estimated using the kernel density estimation method [25, 45]. While the bandwidth is difficult to determine in kernel density estimation, it is associated with sample size [46, 47]. The mesh points that is sample size () was experimentally determined; therefore, the bandwidth was manually determined, because the brain regions were parcellated into a similar size. Because the estimated regional distributions provide a description of brain regions, similarity based on morphological distribution can present a reasonable way of quantifying relationship between two regions. The JSD was converted to a similarity measure (JSS), as follows:
Note that the JSS ranges from 0 to 1 for two probability distributions. Finally, each hemispheric network was binarized over sparsity ranges from 10% to 36% at 0.02 intervals. Sparsity is defined as the ratio of the total number of edges to the maximum possible number of edges in a network. The minimum sparsity calculated from the mean degree over all the nodes should be at least , or 10%, and the maximum sparsity was chosen to allow prominent smallworld properties in the brain network, at 36% [16, 48, 49]. Note that we did not covariate the constructed network matrix for any score of brain lateralization (i.e., the Edinburgh score).
2.4. Network Analysis: Global Graph Metrics
The global graph metrics were calculated for each hemispheric morphological network at a sparsity level of (with 0.02 interval) using the Brain Connectivity Toolbox (BCT) [50]. The global graph metrics included the (1) normalized characteristic path length, (2) normalized clustering coefficient, (3) smallworldness, (4) global efficiency, and (5) local efficiency.
First, the characteristic path length (a global graph metric reflecting the degree of network integration and information transfer) was derived as a harmonic mean of the shortest path lengths between two different nodes (; is the number of nodes comprising a given network (), , , ) comprising a network () [49]:
Second, the clustering coefficient (a graph metric of network segregation) was calculated by averaging the fraction of triangles around each node (or the fraction of node’s neighbors that were also connected each; is the number of nodes within a network () connected to the th node, and is the number of neighboring edges around the th node) across the whole network () [49]:
In addition, the normalized characteristic path length and normalized clustering coefficient values were derived by dividing the original or values using the averaged values of characteristic path length or clustering coefficient for a total of 100 random networks ( or , respectively; constructed by random shuffling of edges within the original network, the number of nodes and edges, and the degree (number of edges connected to each node) distribution preserved) [51].
Third, the smallworldness (a global graph metric for the balance of network integration and network segregation ) was defined as [49]
Fourth, the global efficiency (a global graph metric of network integration and information transfer) was defined as an average of the inverse shortest path length [52]:
Fifth, the local efficiency (a graph metric of local network segregation) was calculated by averaging the changed values of local clustering among the neighbors of th node ( is the number of nodes comprising a given network (), ) when the th node is deleted across the whole nodes within a network () [52]:
2.5. Network Analysis: Local Graph Metrics
The local graph metrics were calculated for each hemispheric morphological network at a sparsity level of (with 0.02 interval) using BCT. The (1) degree centrality and (2) nodal efficiency were calculated to explore the local network properties of the hemispheric morphological network.
First, the degree centrality of the th node was defined as the number of links ( is the presence (1) or absence (0) of edge connecting the th and th nodes within a network (); is the number of nodes comprising a given network (), , , ) connected to the th node of a network ():
Second, the nodal efficiency of the th node was defined as the global efficiency value computed on the neighborhood of the th node and is related to the clustering coefficient [53]:
2.6. Network Analysis: Integrated Graph Metrics
The integrated graph metrics were obtained because a specific sparsity level selection may not be sufficient to reveal the topological properties [13]. The integrated global graph metrics were defined as follows:
The parameter is the set of sparsity values with a 0.02 interval over a range of 0.1–0.36. The parameter is the number of elements in the set of sparsity values, and is the global network metric (, , , , and ) at sparsity . The integrated local graph metrics of node are defined similarly as follows:
The parameter is the set of sparsity values with an interval 0.02 over the network sparsity range of 0.1–0.36. The parameter is the number of elements in the set of sparsity values, and is the mean local network metric (, ) of a region at sparsity . The parameter is one of the 78 cortical regions within the AAL atlas.
2.7. Asymmetry Index of Global and Local Graph Metrics
For the global and local graph metrics, the asymmetry indexes (AIs) of the hemispheric morphological network were calculated as follows:
The parameters and are the integrated global and local graph metrics of the right and left hemispheric networks, respectively. The integrated global and local graph metrics of the hemispheric network were used as the summary of the graph metrics over the sparsity range. Note that positive means leftward asymmetry and negative rightward asymmetry.
2.8. Statistical Analyses
A twoway repeatedmeasures analysis of variance (ANOVA) was performed with hemisphere (left and right) as a repeatedmeasures factor, gender as a betweensubject factor, and age as covariate, to determine whether there were significant effects on any of the global graph metrics. A onesample test within each gender was performed to determine whether the global and local metrics in the hemispheric networks of each gender showed asymmetry. A twosample test between each gender was performed to identify whether the AIs of the global and local metrics were different. The significance threshold of the statistical analysis was set at for the global and local integrated metrics. The significance threshold was set at the false discovery rate (FDR) corrected value for the local graph metrics. All statistical analyses were performed in R.
3. Results
3.1. SmallWorldness and Global/Local Efficiencies of Hemispheric Brain Morphological Networks
For both male and female, the hemispheric brain morphological networks showed smallworldness since is larger than 1.5 over the entire range (see Figure 2) [49], in addition to the higher and lower than those of the matched random networks (see Figure 3) in the network sparsity range of . The results are consistent with those of previous hemispheric brain network studies [13, 14].
(a)
(b)
3.2. Gender and Hemispheric Effect on Global and Local Graph Metrics
The asymmetries of global graph metrics showed a significant hemisphere effect on (, ), a significant gender effect on (, ), and hemispheregender interaction effects on and (, ; , ) (see Table S2). The asymmetries of local graph metrics showed significant hemisphere effects (FDRcorrected ) on 38 nodes and significant hemispheregender interaction effects (FDRcorrected ) on 2 nodes (see Figure 4). For the degree, we observed a significant hemisphere effect (FDRcorrected ) on 33 nodes, a significant gender effect (FDRcorrected ) on 7 nodes, and a significant hemispheregender interaction effect (FDRcorrected ) on 2 nodes (see Figure 4). Significant hemispheric effects on two local graph metrics were observed on nearly the same nodes. Both global and local characteristics of morphological network showed not only hemispheric effect but also gender effect. Furthermore, there is interaction effect also in the global and local characteristics of morphological network.
3.3. The Asymmetry Index of the Global and Local Graph Metrics within Groups
Females showed significant rightward hemispheric asymmetries in (, ) and significant leftward hemispheric asymmetry in ( and ). Males showed significant rightward hemispheric asymmetry in (, ) (Table 1).

For asymmetries of nodal global efficiency, females showed significant asymmetries on 35 nodes (20 leftward asymmetries and 15 rightward asymmetries), and males showed significant asymmetries on 29 nodes (15 leftward asymmetries and 14 rightward asymmetries) at the FDRcorrected value < 0.05 (Figure 5); asymmetries of nodal global efficiencies are in the same direction on 28 nodes in both females and males. For asymmetries of degree centrality, females showed significant asymmetries on 28 nodes (12 leftward asymmetries and 16 rightward asymmetries), and males showed significant asymmetries on 31 nodes (16 leftward asymmetries and 15 rightward asymmetries) at the FDRcorrected value < 0.05 (see Figure 5). The asymmetries of the females and males are in the same direction in 26 nodes. While local graph metrics showed leftward and rightward asymmetries together in females and males, global graph metrics of the morphological network tended leftward in both females and males.
3.4. BetweenGroup Differences in Asymmetry Index of Global and Local Graph Metrics
For the asymmetries of global graph metrics, significant group differences were observed on , , and (, , ; , , ; , , and ) (see Table 1). The females showed significantly greater rightward asymmetry in and significantly greater leftward asymmetry in . Note that the degree of asymmetry in the females is greater. While the females showed rightward asymmetry, the males showed leftward asymmetry on .
For asymmetries of nodal global efficiency, the males showed greater rightward asymmetries on 2 nodes (median cingulate and paracingulate gyri, , FDRcorrected , mean of , mean of ; superior parietal gyrus, , FDRcorrected , mean of , mean of ) (see Figure 5). For asymmetries of degree centrality, the males showed leftward asymmetries, and the females showed rightward asymmetries on 1 node (the gyrus rectus, , FDRcorrected , mean of , mean of ) and showed greater leftward asymmetries on 1 node (temporal pole: middle temporal gyrus, , FDRcorrected , mean of , mean of ).
Opposite directionality of asymmetries (leftward versus rightward) was shown for global metrics of clustering coefficient, characteristic path length, and global efficiency. For the local metrics, larger rightward asymmetries of cingulatesuperior parietal gyri for nodal global efficiency in male compared to female, larger leftward asymmetry of temporal pole for degree centrality in female compared to male, and opposite directionality of interhemispheric asymmetry of rectal gyrus for degree centrality between female (rightward) and male (leftward) were shown.
4. Discussion
The current study examined genderrelated differences in hemispheric asymmetries of graph metrics, calculated from cortical thicknessbased brain structural covariance network named hemispheric morphological network. Hemispheric morphological networks satisfied smallworldness and global efficiency for the network sparsity ranges of . Betweengroup comparisons (female versus male) of asymmetry indexes revealed opposite directionality of asymmetries (leftward versus rightward) for global metrics of clustering coefficient, characteristic path length, and global efficiency. For the local metrics, larger rightward asymmetries of cingulatesuperior parietal gyri for nodal global efficiency in male compared to female, larger leftward asymmetry of temporal pole for degree centrality in female compared to male, and opposite directionality of interhemispheric asymmetry of rectal gyrus for degree centrality between female (rightward) and male (leftward) were shown. The overall results of this study indicate that brain network analysis using morphological features provides insights into the understanding of hemispheric asymmetry related to gender.
4.1. Economical SmallWorld Network Properties of the Hemispheric Brain Morphological Networks
The identification of economical smallworld properties of the human brain network could help in the understanding of the human brain [4, 15]. It has been shown that the whole brain network of healthy subjects is characterized by economical smallworld properties [3, 54–56]. In an individual morphological brain network, the presence of economical smallworld properties has been reported as convergent evidence [24, 25, 57]. Furthermore, the presence of economical smallworld network properties in hemispheric networks was reported using sMRI and fMRI [12, 13, 58]. Our results extend previous findings by showing that individual hemispheric morphological networks also demonstrate economical smallworld properties, indicating that information processing within each hemisphere could be of similar efficiency to that of the whole brain.
4.2. GenderRelated Effects in Asymmetry of Global Graph Metrics
The left hemisphere in females may be more efficient in the exchange of information in parallel. The in females shows leftward asymmetry. Although the left and right hemispheres of the human brain consistently communicate with each other, asymmetry of the human brain hemispheric network and its differences with respect to gender have been reported. For instance, IturriaMedina et al. showed that the right hemisphere is more efficient with rightward and than the left hemisphere in righthanded individuals [12]. Shang et al. showed that the functional left hemispheric network is more globally efficient than the functional right hemispheric network, and Sun et al. reported, using the structural hemispheric network, that the left hemispheric network is globally efficient [59, 60]. Gong et al. showed that the cortical anatomical network is more efficient globally and locally in females than in males [61]. Tian et al. reported that the two hemispheres and genders are not significantly different in transferring information between the brain regions; however, females and males are both more globally efficient in their right hemispheres [13]. These findings are different, depending on the modality and methodology used. Some aspects of asymmetry observed in this study are reflected observations based on the structural and functional networks. This study suggests that the left hemispheric morphological network is more globally efficient in females, in accordance with previous observations based on the structural and functional networks [59–61].
4.3. GenderRelated Effects in Asymmetry of Local Graph Metrics
It was observed that males and females have similar asymmetric patterns in nodal global efficiency and degree centrality. Significant leftward asymmetries were observed in the frontal region, precuneus, and temporal region, while significant rightward asymmetries were observed in the cingulate gyrus and occipital gyrus. These patterns have been observed to be structurally asymmetric in previous studies [62]. This asymmetry pattern along the frontooccipital axis is similar to that reported by Plessen et al. [63] and may be related to the Yakovlevian torque, where frontal/occipital bending in the human brain is present [64]. However, some regions related to behavior (properties) differences showed gender difference in asymmetry. A gender difference in asymmetries was observed in four brain regions (: the MCG, SPG; Degree: REC and TPOmid) in this study (see Figure 5). Most of these regions, such as the MCG, SPG, and TPOmid, have been reported as hubs in previous studies based on the structural and functional networks [65]. Given that the left hemisphere is dominant in language processing and the right hemisphere is dominant in spatial processing [66], the asymmetries of the local graph metrics may underlie advantageous verbal processing in females and advantageous spatial processing in males [67]. The nodal global efficiency in the MCG and SPG regions shows more rightward asymmetry in males than in females. In particular, the visuospatial processing regions, such as the right SPG, showed greater nodal global efficiency in males, which may directly contribute to the previous observation [68, 69]. These regions were also reported in the studies on the nodal global efficiency of the brain network [22, 23, 61]. For instance, Zhong et al. reported that the nodal global efficiency of the hemispheric structural network is more rightward asymmetric on the MCG and SPG. Caeyenberghs and Leemans observed rightward asymmetry of the nodal global efficiency in the MCG region using the structural network [23]. The rightward asymmetry of the nodal global efficiency in these regions is consistent with those found in previous studies using the structural network. The right hemispheric morphological network is more locally efficient in males than in females. The degree centrality in TPOmid region shows more leftward asymmetry in females than in males. The structural leftward asymmetry of the TPOmid region was studied [70, 71]. Price reported left hemispheric dominance for language in the middle temporal pole [72]. The significant leftward asymmetry of the degree centrality in the TPOmid regions is consistent with those found in previous studies. The degree centrality in the REC region shows rightward asymmetry in females and leftward asymmetry in males. The REC region is known as a nonfunctional gyrus [73]; hence, it is difficult to interpret gender difference from this region. For example, Belfi et al. reported that the REC region, a narrow strip of the cerebral cortex, has a larger volume in females; however, the gender difference could be in terms of psychological gender rather than biological gender. Therefore, higher femininity scores in gender assessments of males were associated with a larger volume of the REC region [74]. The opposing directions of asymmetry of the degree centrality in the REC regions between females and males might reflect the individual’s personality rather than the biological gender.
4.4. Limitations
There are still some limitations in this study. First, the selection of the brain atlas could affect the topological properties of the individual brain network [57]. In the future, it is important to validate this experiment with a different atlas. Second, we used cortical thickness as a morphological feature to construct a morphological brain network. Several recent studies have reported methodologies to construct networks using multiple features [75, 76]. In the future, it would be interesting to use multiple features to construct the brain network using our method. However, it is noticeable that an exact physiological explanation of this network is still difficult. Third, we observed significant hemisphere and gender effects in topological properties of the morphological brain network. Some aspects of these effects observed in this study are consistent with those based on structural and functional brain networks, while others are not [13, 61]. Because morphological similarity is related to the underlying axonal connectivities [77], a further study simultaneously evaluating the topologies of structural networks and of functional networks is expected [78–80]. Fourth, we observed significant asymmetries related to gender in topological properties of the morphological network. Brain asymmetry is closely related to lateralized behaviors. Thus, further analysis should be implemented for a better understanding of the basis of lateralized functions such as language and visuospatial processing. Finally, changes in brain asymmetries are closely related to the pathophysiology of various brain diseases such as schizophrenia and Alzheimer’s disease. Therefore, evaluating the topological organization of morphological brain networks within hemispheres is likely to improve our knowledge about the pathology of various brain diseases.
5. Conclusions
In this study, hemispheric morphological brain networks were constructed, and graphtheoretical approaches were used to examine the hemispheric asymmetry in females and males. First, the hemispheric morphological networks showed smallworld properties and a high efficiency, similar to those of structural and functional networks. Second, the hemispheric asymmetry within and between the gender groups was investigated. In the global scale, males tend to be locally efficient in the left hemispheric network, and females are more globally efficient in the left hemispheric network. In the local scale, the right asymmetry in the nodal global efficiency is greater around the cingulate gyrus and superior parietal gyrus in males than in females. The left asymmetry in degree centrality is greater around the temporal pole in females than in males. These findings may provide evidence for the topological difference in the hemispheric morphological network and for the behavioral differences related to gender. The overall results of this study indicate that brain network analysis using morphological features provides insights into the understanding of hemispheric asymmetry related to gender.
Data Availability
The brain MRI dataset analyzed during this study is available in the Human Connectome Project repository (http://www.humanconnectome.org/).
Conflicts of Interest
The authors have no conflicts of interest to declare.
Acknowledgments
This work was supported by Institute of Information & Communications Technology Planning & Evaluation (IITP) grant funded by the Korea government (MSIT) (No. 2020001373, Artificial Intelligence Graduate School Program (Hanyang University)). Data were provided by the Human Connectome Project, WUMinn Consortium (principal investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657), funded by the 16 NIH institutes and centers that support the NIH Blueprint for Neuroscience Research and by the McDonnell Center for Systems Neuroscience at Washington University.
Supplementary Materials
Table S1: abbreviations for the brain regions of interest. Table S2: hemisphere and gender effects on the global integrated measures of the hemispheric network revealed by twoway repeatedmeasures ANOVA. Twoway repeatedmeasures ANOVA was performed to investigate hemisphere and gender effects and interaction effects on the global integrated measures. The significance level for all analyses is . , , , , and denote the clustering coefficient, characteristic path length, smallworldness, global efficiency, and local efficiency, respectively. The significance () is indicated by the bold text and indicator (). (Supplementary materials)
References
 A. AlexanderBloch, J. N. Giedd, and E. Bullmore, “Imaging structural covariance between human brain regions,” vol. 14, no. 5, pp. 322–336, 2013. View at: Publisher Site  Google Scholar
 J. Y. Yun, P. S. W. Boedhoe, C. Vriend et al., “OUP accepted manuscript,” Brain, vol. 143, no. 2, pp. 684–700, 2020. View at: Publisher Site  Google Scholar
 P. Hagmann, M. Kurant, X. Gigandet et al., “Mapping human wholebrain structural networks with diffusion MRI,” PLoS One, vol. 2, no. 7, article e597, 2007. 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
 G. Ball, R. Beare, and M. L. Seal, “Charting shared developmental trajectories of cortical thickness and structural connectivity in childhood and adolescence,” Human Brain Mapping, vol. 40, no. 16, pp. 4630–4644, 2019. View at: Publisher Site  Google Scholar
 D. A. Hoagey, J. R. Rieck, K. M. Rodrigue, and K. M. Kennedy, “Joint contributions of cortical morphometry and white matter microstructure in healthy brain aging: a partial least squares correlation analysis,” Human Brain Mapping, vol. 40, no. 18, pp. 5315–5329, 2019. View at: Publisher Site  Google Scholar
 W. Sato and S. Uono, “The atypical social brain network in autism: advances in structural and functional MRI studies,” Current Opinion in Neurology, vol. 32, no. 4, pp. 617–621, 2019. View at: Publisher Site  Google Scholar
 L. G. Matthews, C. D. Smyser, S. Cherkerzian et al., “Maternal pomegranate juice intake and brain structure and function in infants with intrauterine growth restriction: a randomized controlled pilot study,” PLoS One, vol. 14, no. 8, article e0219596, 2019. View at: Publisher Site  Google Scholar
 F. Han, Y. Gu, G. L. Brown, X. Zhang, and X. Liu, “Neuroimaging contrast across the cortical hierarchy is the feature maximally linked to behavior and demographics,” Neuroimage, vol. 215, article 116853, 2020. View at: Publisher Site  Google Scholar
 R. J. Jirsaraie, A. N. Kaczkurkin, S. Rush et al., “Accelerated cortical thinning within structural brain networks is associated with irritability in youth,” Neuropsychopharmacology, vol. 44, no. 13, pp. 2254–2262, 2019. View at: Publisher Site  Google Scholar
 T. Qi, G. Schaadt, R. Cafiero, J. Brauer, M. A. Skeide, and A. D. Friederici, “The emergence of longrange language network structural covariance and language abilities,” Neuroimage, vol. 191, pp. 36–48, 2019. View at: Publisher Site  Google Scholar
 Y. IturriaMedina, A. P. Fernández, D. M. Morris et al., “Brain hemispheric structural efficiency and interconnectivity rightward asymmetry in human and nonhuman primates,” Cerebral Cortex, vol. 21, no. 1, pp. 56–67, 2011. View at: Publisher Site  Google Scholar
 L. Tian, J. Wang, C. Yan, and Y. He, “Hemisphere and genderrelated differences in smallworld brain networks: a restingstate functional MRI study,” NeuroImage, vol. 54, no. 1, pp. 191–202, 2011. View at: Publisher Site  Google Scholar
 M. Li, H. Chen, J. Wang et al., “Handedness and hemisphererelated differences in smallworld brain networks: a diffusion tensor imaging tractography study,” Brain Connectivity, vol. 4, no. 2, pp. 145–156, 2014. View at: Publisher Site  Google Scholar
 O. Sporns, “Structure and function of complex brain networks,” Dialogues in Clinical Neuroscience, vol. 15, no. 3, pp. 247–262, 2013. View at: Google Scholar
 D. S. Bassett, E. Bullmore, B. A. Verchinski, V. S. Mattay, D. R. Weinberger, and A. MeyerLindenberg, “Hierarchical organization of human cortical networks in health and schizophrenia,” The Journal of Neuroscience, vol. 28, no. 37, pp. 9239–9248, 2008. View at: Publisher Site  Google Scholar
 K. Supekar, V. Menon, D. Rubin, M. Musen, and M. D. Greicius, “Network analysis of intrinsic functional brain connectivity in Alzheimer’s disease,” PLoS Computational Biology, vol. 4, no. 6, article e1000100, 2008. View at: Publisher Site  Google Scholar
 C.Y. Lo, P.N. Wang, K.H. Chou, J. Wang, Y. He, and C.P. Lin, “Diffusion tensor tractography reveals abnormal topological organization in structural cortical metworks in Alzheimer’s disease,” The Journal of Neuroscience, vol. 30, no. 50, pp. 16876–16885, 2010. View at: Publisher Site  Google Scholar
 A. Griffa, P. S. Baumann, J. P. Thiran, and P. Hagmann, “Structural connectomics in brain diseases,” NeuroImage, vol. 80, pp. 515–526, 2013. View at: Publisher Site  Google Scholar
 T. Wang, K. Wang, H. Qu et al., “Disorganized cortical thickness covariance network in major depressive disorder implicated by aberrant hubs in largescale networks,” Scientific Reports, vol. 6, no. 1, 2016. View at: Publisher Site  Google Scholar
 E. Collantoni, P. Meneguzzo, E. Tenconi, R. Manara, and A. Favaro, “Smallworld properties of brain morphological characteristics in anorexia nervosa,” PLoS One, vol. 14, no. 5, p. e0216154, 2019. View at: Publisher Site  Google Scholar
 S. Zhong, Y. He, H. Shu, and G. Gong, “Developmental changes in topological asymmetry between hemispheric brain white matter networks from adolescence to young adulthood,” Cerebral Cortex, vol. 27, no. 4, pp. 2560–2570, 2017. View at: Publisher Site  Google Scholar
 K. Caeyenberghs and A. Leemans, “Hemispheric lateralization of topological organization in structural brain networks,” Human Brain Mapping, vol. 35, no. 9, pp. 4944–4957, 2014. View at: Publisher Site  Google Scholar
 B. M. Tijms, P. Series, D. J. Willshaw, and S. M. Lawrie, “Similaritybased extraction of individual networks from gray matter MRI scans,” Cerebral Cortex, vol. 22, no. 7, pp. 1530–1541, 2012. View at: Publisher Site  Google Scholar
 X. Z. Kong, Z. Liu, L. Huang et al., “Mapping individual brain networks using statistical similarity in regional morphology from MRI,” PLoS One, vol. 10, no. 11, article e0141840, 2015. View at: Publisher Site  Google Scholar
 H. J. Kim, J. H. Shin, C. E. Han et al., “Using individualized brain network for analyzing structural covariance of the cerebral cortex in Alzheimer’s patients,” Frontiers in Neuroscience, vol. 40, p. 394, 2016. View at: Publisher Site  Google Scholar
 X. Z. Kong, X. Wang, L. Huang et al., “Measuring individual morphological relationship of cortical regions,” Journal of Neuroscience Methods, vol. 237, pp. 103–107, 2014. View at: Publisher Site  Google Scholar
 D. M. Endres and J. E. Schindelin, “A new metric for probability distributions,” IEEE Transactions on Information Theory, vol. 49, no. 7, pp. 1858–1860, 2003. View at: Publisher Site  Google Scholar
 B. Fuglede and F. Topsoe, “JensenShannon divergence and Hilbert space embedding,” in International Symposium on Information Theory, 2004, ISIT 2004, Proceedings, p. 31, Chicago, IL, USA, 2004. View at: Publisher Site  Google Scholar
 R. Core Team, R: a language and environment for statistical computing, 2017.
 D. C. Van Essen, K. Ugurbil, E. Auerbach et al., “The Human Connectome Project: a data acquisition perspective,” NeuroImage, vol. 62, no. 4, pp. 2222–2231, 2012. View at: Publisher Site  Google Scholar
 D. C. Van Essen, S. M. Smith, D. M. Barch, T. E. J. Behrens, E. Yacoub, and K. Ugurbil, “The WUMinn Human Connectome Project: an overview,” NeuroImage, vol. 80, pp. 62–79, 2013. View at: Publisher Site  Google Scholar
 M. F. Glasser, S. N. Sotiropoulos, J. A. Wilson et al., “The minimal preprocessing pipelines for the Human Connectome Project,” Neuroimage, vol. 80, pp. 105–124, 2013. View at: Publisher Site  Google Scholar
 D. L. Collins, P. Neelin, T. M. Peters, and A. C. Evans, “Automatic 3D intersubject registration of MR volumetric data in standardized Talairach space,” Journal of Computer Assisted Tomography, vol. 18, no. 2, pp. 192–205, 1994. View at: Publisher Site  Google Scholar
 J. G. Sled, A. P. Zijdenbos, and A. C. Evans, “A nonparametric method for automatic correction of intensity nonuniformity in MRI data,” IEEE Transactions on Medical Imaging, vol. 17, no. 1, pp. 87–97, 1998. View at: Publisher Site  Google Scholar
 S. M. Smith, “Fast robust automated brain extraction,” Human Brain Mapping, vol. 17, no. 3, pp. 143–155, 2002. View at: Publisher Site  Google Scholar
 A. Zijdenbos, A. Evans, F. Riahi, J. Sled, J. Chui, and V. Kollokian, “Automatic quantification of multiple sclerosis lesion volume using stereotaxic space,” in Visualization in Biomedical Computing. VBC 1996, K. H. Höhne and R. Kikinis, Eds., vol. 1131 of Lecture Notes in Computer Science, pp. 439–448, Springer, Berlin, Heidelberg, 1996. View at: Publisher Site  Google Scholar
 J. S. Kim, V. Singh, J. K. Lee et al., “Automated 3D extraction and evaluation of the inner and outer cortical surfaces using a Laplacian map and partial volume effect classification,” NeuroImage, vol. 27, no. 1, pp. 210–221, 2005. View at: Publisher Site  Google Scholar
 K. Im, J.M. Lee, J. Lee et al., “Gender difference analysis of cortical thickness in healthy young adults with surfacebased methods,” NeuroImage, vol. 31, no. 1, pp. 31–38, 2006. View at: Publisher Site  Google Scholar
 N. TzourioMazoyer, B. Landeau, D. Papathanassiou et al., “Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI singlesubject brain,” NeuroImage, vol. 15, no. 1, pp. 273–289, 2002. View at: Publisher Site  Google Scholar
 A. M. Hermundstad, D. S. Bassett, K. S. Brown et al., “Structural foundations of restingstate and taskbased functional connectivity in the human brain,” Proceedings of the National Academy of Sciences, vol. 110, no. 15, pp. 6169–6174, 2013. View at: Publisher Site  Google Scholar
 Z. Zhang, E. Klassen, and A. Srivastava, “Robust comparison of kernel densities on spherical domains,” Sankhya A, vol. 81, pp. 144–171, 2018. View at: Publisher Site  Google Scholar
 D. S. Bassett, J. A. Brown, V. Deshpande, J. M. Carlson, and S. T. Grafton, “Conserved and variable architecture of human white matter connectivity,” Neuroimage, vol. 54, no. 2, pp. 1262–1279, 2011. View at: Publisher Site  Google Scholar
 P. Hagmann, L. Cammoun, X. Gigandet et al., “Mapping the structural core of human cerebral cortex,” PLoS Biology, vol. 6, no. 7, article e159, 2008. View at: Publisher Site  Google Scholar
 Z. I. Botev, J. F. Grotowski, and D. P. Kroese, “Kernel density estimation via diffusion,” The Annals of Statistics, vol. 38, no. 5, pp. 2916–2957, 2010. View at: Publisher Site  Google Scholar
 N. G. Ushakov and V. G. Ushakov, “On bandwidth selection in kernel density estimation,” Journal of Nonparametric Statistics, vol. 24, no. 2, pp. 419–428, 2012. View at: Publisher Site  Google Scholar
 S. Chen, “Optimal Bandwidth Selection for Kernel Density Functionals Estimation,” Journal of Probability and Statistics, vol. 2015, Article ID 242683, 21 pages, 2015. View at: Publisher Site  Google Scholar
 W. Liao, J. Ding, D. Marinazzo et al., “Smallworld directed networks in the human brain: multivariate Granger causality analysis of restingstate fMRI,” NeuroImage, vol. 54, no. 4, pp. 2683–2694, 2011. 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: Publisher Site  Google Scholar
 M. Rubinov and O. Sporns, “Complex network measures of brain connectivity: uses and interpretations,” NeuroImage, vol. 52, no. 3, pp. 1059–1069, 2010. View at: Publisher Site  Google Scholar
 S. Maslov and K. Sneppen, “Specificity and stability in topology of protein networks,” Science, vol. 296, no. 5569, pp. 910–913, 2002. View at: Publisher Site  Google Scholar
 V. Latora and M. Marchiori, “Efficient behavior of smallworld networks,” Physical Review Letters, vol. 87, no. 19, article 198701, 2001. 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, p. e17, 2007. View at: Publisher Site  Google Scholar
 D. S. Bassett and E. Bullmore, “Smallworld brain networks,” The Neuroscientist, vol. 12, no. 6, pp. 512–523, 2007. View at: Publisher Site  Google Scholar
 Y. He, Z. J. Chen, and A. C. Evans, “Smallworld anatomical networks in the human brain revealed by cortical thickness from MRI,” Cerebral Cortex, vol. 17, no. 10, pp. 2407–2419, 2007. View at: Publisher Site  Google Scholar
 J. Wang, “Graphbased network analysis of restingstate functional MRI,” Frontiers in Systems Neuroscience, vol. 4, p. 16, 2010. View at: Publisher Site  Google Scholar
 H. Wang, X. Jin, Y. Zhang, and J. Wang, “Singlesubject morphological brain networks: connectivity mapping, topological characterization and testretest reliability,” Brain and Behavior: A Cognitive Neuroscience Perspective, vol. 6, no. 4, pp. 1–21, 2016. View at: Google Scholar
 Y. Sun, Y. Chen, S. L. Collinson, A. Bezerianos, and K. Sim, “Reduced hemispheric asymmetry of brain anatomical networks is linked to schizophrenia: a connectome study,” Cerebral Cortex, vol. 27, pp. 602–615, 2015. View at: Publisher Site  Google Scholar
 Y. Sun, J. Li, J. Suckling, and L. Feng, “Asymmetry of hemispheric network topology reveals dissociable processes between functional and structural brain connectome in communityliving elders,” Frontiers in Aging Neuroscience, vol. 9, p. 361, 2017. View at: Publisher Site  Google Scholar
 W. Shang, T. Li, J. Xiang et al., “Hemispheric asymmetry of the functional brain connectome,” in Cognitive Systems and Signal Processing. ICCSIP 2016, F. Sun, H. Liu, and D. Hu, Eds., vol. 710 of Communications in Computer and Information Science, pp. 541–547, Springer, Singapore, 2017. View at: Publisher Site  Google Scholar
 G. Gong, P. RosaNeto, F. Carbonell, Z. J. Chen, Y. He, and A. C. Evans, “Age and genderrelated differences in the cortical anatomical network,” The Journal of Neuroscience, vol. 29, no. 50, pp. 15684–15693, 2009. View at: Publisher Site  Google Scholar
 T. Guadalupe, D. C. Glahn, P. M. Thompson et al., “Mapping cortical brain asymmetry in 17,141 healthy individuals worldwide via the ENIGMA consortium,” Proceedings of the National Academy of Sciences of the United States of America, vol. 115, no. 22, pp. E5154–E5163, 2018. View at: Google Scholar
 K. J. Plessen, K. Hugdahl, R. Bansal, X. Hao, and B. S. Peterson, “Sex, age, and cognitive correlates of asymmetries in thickness of the cortical mantle across the life span,” The Journal of Neuroscience, vol. 34, no. 18, pp. 6294–6302, 2014. View at: Publisher Site  Google Scholar
 J. J. Maller, R. Anderson, R. H. Thomson, J. V. Rosenfeld, Z. J. Daskalakis, and P. B. Fitzgerald, “Occipital bending (Yakovlevian torque) in bipolar depression,” Psychiatry Research: Neuroimaging, vol. 231, no. 1, pp. 8–14, 2015. View at: Publisher Site  Google Scholar
 M. P. van den Heuvel and O. Sporns, “Network hubs in the human brain,” Trends in Cognitive Sciences, vol. 17, no. 12, pp. 683–696, 2013. View at: Publisher Site  Google Scholar
 S. P. Springer and G. Deutsch, “Left brain, right brain: perspectives from cognitive neuroscience,” in A Series of Books in Psychology, W H Freeman/Times Books/Henry Holt & Co., 5th edition, 1998. View at: Google Scholar
 E. M. Weiss, G. Kemmler, E. A. Deisenhammer, W. W. Fleischhacker, and M. Delazer, “Sex differences in cognitive functions,” vol. 35, no. 4, pp. 863–875, 2003. View at: Publisher Site  Google Scholar
 E. Zarahn, G. Aguirre, and M. D’Esposito, “Replication and further studies of neural mechanisms of spatial mnemonic processing in humans,” vol. 9, no. 1, pp. 1–17, 2000. View at: Publisher Site  Google Scholar
 R. Everts, K. Lidzba, M. Wilke et al., “Strengthening of laterality of verbal and visuospatial functions during childhood and adolescence,” Human Brain Mapping, vol. 30, no. 2, pp. 473–483, 2009. View at: Publisher Site  Google Scholar
 K. E. Watkins, T. Paus, J. P. Lerch et al., “Structural asymmetries in the human brain: a voxelbased statistical analysis of 142 MRI scans,” Cerebral Cortex, vol. 11, no. 9, pp. 868–877, 2001. View at: Publisher Site  Google Scholar
 E. Luders, K. L. Narr, P. M. Thompson, D. E. Rex, L. Jancke, and A. W. Toga, “Hemispheric asymmetries in cortical thickness,” Cerebral Cortex, vol. 16, no. 8, pp. 1232–1238, 2006. View at: Publisher Site  Google Scholar
 C. J. Price, “The anatomy of language: contributions from functional neuroimaging,” Journal of Anatomy, vol. 197, no. 3, pp. 335–359, 2000. View at: Publisher Site  Google Scholar
 M. S. Joo, D. S. Park, C. T. Moon, Y. Il Chun, S. W. Song, and H. G. Roh, “Relationship between gyrus rectus resection and cognitive impairment after surgery for ruptured anterior communicating artery aneurysms,” Journal of Cerebrovascular and Endovascular Neurosurgery, vol. 18, no. 3, pp. 223–228, 2016. View at: Publisher Site  Google Scholar
 A. M. Belfi, A. L. Conrad, J. Dawson, and P. Nopoulos, “Masculinity/femininity predicts brain volumes in normal healthy children,” Developmental Neuropsychology, vol. 39, no. 1, pp. 25–36, 2009. View at: Publisher Site  Google Scholar
 J. Seidlitz, P. B. Jones, A. Raznahan et al., “Morphometric similarity networks detect microscale cortical organization and predict interindividual cognitive variation,” Neuron, vol. 97, no. 1, pp. 231–247.e7, 2017. View at: Publisher Site  Google Scholar
 X. Li, X. Wang, K. Yu, Q. Li, X. Zhang, and S. Li, “Individual morphological brain network construction based on multivariate Euclidean distances between brain regions,” Frontiers in Human Neuroscience, vol. 12, 2018. View at: Publisher Site  Google Scholar
 D. C. Van Essen, “A tensionbased theory of morphogenesis and compact wiring in the central nervous system,” vol. 385, no. 6614, pp. 313–318, 1997. View at: Publisher Site  Google Scholar
 R. Meuli, C. J. Honey, O. Sporns et al., “Predicting human restingstate functional connectivity from structural connectivity,” Proceedings of the National Academy of Sciences, vol. 106, no. 6, pp. 2035–2040, 2009. View at: Publisher Site  Google Scholar
 N. Langer, A. Pedroni, and L. Jäncke, “The problem of thresholding in smallworld network analysis,” PLoS One, vol. 8, no. 1, article e53199, 2013. View at: Publisher Site  Google Scholar
 M. P. van den Heuvel, S. C. de Lange, A. Zalesky, C. Seguin, B. T. T. Yeo, and R. Schmidt, “Proportional thresholding in restingstate fMRI functional connectivity networks and consequences for patientcontrol connectome studies: issues and recommendations,” NeuroImage, vol. 152, pp. 437–449, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 YongHo Choi 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.