Gender-Related Differences in the Dysfunctional Resting Networks of Migraine Suffers

Abstract

Background

Migraine shows gender-specific incidence and has a higher prevalence in females. However, little is known about gender-related differences in dysfunctional brain organization, which may account for gender-specific vulnerability and characteristics of migraine. In this study, we considered gender-related differences in the topological property of resting functional networks.

Methodology/Principal Findings

Data was obtained from 38 migraine patients (18 males and 20 females) and 38 healthy subjects (18 males and 20 females). We used the graph theory analysis, which becomes a powerful tool in investigating complex brain networks on a whole brain scale and could describe functional interactions between brain regions. Using this approach, we compared the brain functional networks between these two groups, and several network properties were investigated, such as small-worldness, network resilience, nodal centrality, and interregional connections. In our findings, these network characters were all disrupted in patients suffering from chronic migraine. More importantly, these functional damages in the migraine-affected brain had a skewed balance between males and females. In female patients, brain functional networks showed worse resilience, more regions exhibited decreased nodal centrality, and more functional connections revealed abnormalities than in male patients.

Conclusions

These results indicated that migraine may have an additional influence on females and lead to more dysfunctional organization in their resting functional networks.

About the Authors:

Jixin Liu

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Wei Qin

* E-mail: [email protected] (JT); [email protected] (FL); [email protected] (WQ)

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Jiaofen Nan

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Jing Li

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Kai Yuan

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Ling Zhao

Affiliation: The 3rd Teaching Hospital, Chengdu University of Traditional Chinese Medicine, Chengdu, China

Fang Zeng

Affiliation: The 3rd Teaching Hospital, Chengdu University of Traditional Chinese Medicine, Chengdu, China

Jinbo Sun

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Dahua Yu

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Minghao Dong

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Peng Liu

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Karen M. von Deneen

Affiliation: Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China

Qiyong Gong

Affiliation: Huaxi MR Research Center (HMRRC), Department of Radiology, The Center for Medical Imaging, West China Hospital of Sichuan University, Chengdu, Sichuan, People's Republic of China

Fanrong Liang

* E-mail: [email protected] (JT); [email protected] (FL); [email protected] (WQ)

Affiliation: The 3rd Teaching Hospital, Chengdu University of Traditional Chinese Medicine, Chengdu, China

Jie Tian

* E-mail: [email protected] (JT); [email protected] (FL); [email protected] (WQ)

Affiliations Life Sciences Research Center, School of Life Sciences and Technology, Xidian University, Xi'an, Shaanxi, China, Institute of Automation, Chinese Academy of Sciences, Beijing, China

Introduction

Migraine is a debilitating condition, characterized by severe headaches and nausea, which leads to disability, lost productivity, and a decreased overall quality of life [1]–[3]. Recently, researchers have shown an increased interest in the damage and dysregulation of the central nervous system (CNS) occurring in patients with headaches and migraine. In particular, Filippi et al. (2008), May (2009) and Chiapparini et al. (2010) provided summaries of the main investigations in the modern imaging field [4]–[6]. Based on neuroimaging findings, several brain regions in patients with migraine (PM) have been reported to show structural and functional alteration, including the insula (INS), anterior cingulate cortex (ACC), cerebellum, thalamus, prefrontal cortex (PFC), and temporal cortex [7]–[11]. These brain abnormalities were considered as the secondary effects of migraine, and may be the results of dysfunction from infrequent to frequent cortical overstimulation associated with headaches [4], [6]. Neuroimaging findings provide useful information on brain activity outside of headache attacks, and are just beginning to learn about CNS damage and cortical reorganization in migraine.

Recently, differences in the incidence of migraine between men and women have been observed [12], [13]. Stewart et al. (2008) pointed out that migraine had a cumulative lifetime incidence of 43% in women and 18% in men [12], which was about three times more common in women than in men [13]. As we know, gender is an important influential element in brain neural systems [14]. Anatomical and functional network organization in the human brain is associated with gender [15]–[17]. Yan et al. (2010) investigated gender differences in the structural brain networks and found that females had greater local efficiencies than males did [16]. Tian et al. (2011) demonstrated that there existed gender-related differences in the organization of brain functional networks [17]. Additionally, the gender difference is also an influential factor for behavior and cognitive performance [18], such as for the mental rotation test [19], emotion recognition test [20], and social sensitivity test [21]. However, to date, there has been little discussion on the cortical network organization differences between men and women with migraine. Due to the different incidence of this disease, there is the high possibility that migraine may have a different influence on the organization of cortical networks between males and females.

Previous studies focused on the local variation of regional cerebral blood flow (positron emission tomography, PET) and structural damage (voxel-based morphometry, VBM) [7], [11], [22], [23] in migraine. However, physiological effects of brain injury were best evaluated throughout all of the cortical networks, rather than assessed locally [24]. In this study, we used resting-state functional MRI to consider resting networks which can be conceived as neurocognitive entities that incorporate both local and global processes [25], [26]. Furthermore, differences in resting networks between patients and control subjects could imply differences in an underlying disease-specific brain abnormality [27]. The recent application of graph theory analysis (GTA) has become a powerful tool to investigate complex brain networks on a whole brain scale by defining a graph as a set of nodes (brain regions) and edges (functional connections) [28]–[30]. A graph represents the cortical network, where it describes the basis of cognitive processing for distributed functional interactions between brain regions [31]. Cortical network analysis may accelerate our understating of a dysregulated connectome in the PM.

Here, we investigated the topological properties of brain networks by using GTA and hypothesized that that there was a different connection model which arose from gender differences in PM. To test our aim, we constructed functional brain networks to characterize the interregional relationships between brain regions in PM and healthy controls (HC) respectively. We shed light on an integrative systems' perspective and detected the patients' abnormal network characteristics and their associations with gender.

Materials and Methods

All research procedures were approved by the West China Hospital Subcommittee on Human Studies and were conducted in accordance with the Declaration of Helsinki.

2.1 Participants

Thirty-eight right-handed migraine patients (20 females, 32.3±8.7 years (mean age ± SD); 18 males, 33±7.9 years (mean age ± SD)) were recruited for the study (Table 1 and Table 2). All patients fulfilled the ICHD-II criteria for migraine without aura. According to Detsky et al. (2006), five questions were asked: (1) “Is it a pulsating headache?” (2) “Does it last between 4 and 72 hours without medication?” (3) “Is it unilateral?” (4) “Is there nausea?” and (5) “Is the headache disabling?” [32]. All patients answered “yes” to 4 or more of the 5 questions. In all patients, migraine attack frequency was 6.7±2.1 days/month, migraine attack duration was 14.6±7.1 hours, the migraine duration was 11.3±6.7 years, and the average pain intensity was 5.3±1.6 for the past four weeks. No patients had macroscopic brain T2-visible lesions on MRI scans. No patients had any history of drug abuse. All patients had been free from a typical migraine attack for at least 1 week prior to MRI examination. After scanning, all subjects reported that they did not experience any headaches or migraines during the experiment.

[Figure omitted. See PDF.]

Table 1. Demographic characteristics of subjects.

https://doi.org/10.1371/journal.pone.0027049.t001

[Figure omitted. See PDF.]

Table 2. Demographic characteristics of subjects.

https://doi.org/10.1371/journal.pone.0027049.t002

Thirty-eight age-, education- and gender-matched, healthy, right-handed controls (age 32.6±6.9 years) were recruited from the local community. All HC had no history of major medical illnesses, head trauma, and neuropsychiatric disorders, and had not used prescription medications within the last month. To exclude abnormal neurological findings, all HC were screened by a neurologist specialized in headaches. All subjects gave written, informed consent after the experimental procedures had been fully explained.

2.2 fMRI Data Acquisition

The fMRI experiments were carried out in a 3T GE scanner. We collected resting state scans of 205 continuous EPI functional volumes (TE 30 ms, TR 2 s, matrix 64×64, FOV 240 mm, flip angle 90°) for HC and patient groups. Subjects were instructed to keep their eyes closed, to not think about anything, and to stay awake during scanning. According to the subjects' reports, they all stayed awake during all scans. For spatial normalization and localization, a set of T1-weighted high-resolution anatomical images were acquired (TE 3.39 ms, TR 2.7 s, matrix 256×256, FOV 256 mm, flip angle 7°, in plane resolution 1×1 mm, slice thickness 1 mm).

2.3 Data preprocessing

Image preprocessing was carried out using SPM5 (http://www.fil.ion.ucl.ac.uk/spm). All datasets were initially corrected for temporal offsets using sinc interpolation and head movement-related effects using a six-parameter spatial transformation [33]. To minimize movement artifacts, individuals with an estimated maximum displacement in any direction larger than 1.5 mm or head rotation larger than 1.5° were discarded from the study. No data were excluded under this criterion. The mean EPI image for each subject was then co-registered to a corresponding T1-weighted high-resolution image volume and subsequently spatially normalized to the Montreal Neurological Institute (MNI) echoplanar imaging (EPI) template image and resampled to 2-mm isotropic voxels. In order to avoid local correlations, the spatially normalized data were not spatially smoothed in this study.

2.4 Construction of the cortical brain network

A key issue in characterizing the brain topological network is the construction of the functional connection matrix. In this study, automated anatomically labeled (AAL) template images were used [34], which divided the whole brain into 90 regions of interest (ROIs)(Table 3). These 90 brain regions were considered as a set of nodes in our network analysis. Then, the fMRI time series of all regions was filtered using a bandpass filter (0.01–0.08 Hz) to remove the effects of magnetic field drift and high-frequency noise. After that, the mean time courses from deep white matter, ventricles and the 6 rigid-body motion parameters were regressed out from the filtered time series. Finally, we computed the mean time series of each seed region and obtained a 90*90 matrix of the Pearson correlation coefficients between all possible connections of node pairs [35]–[37].

[Figure omitted. See PDF.]

Table 3. Cortical and subcortical regions defined by the AAL template image in standard stereotaxic space.

https://doi.org/10.1371/journal.pone.0027049.t003

As mentioned in previous studies [36]–[38], we invested the brain's topological properties by way of binarized graphs in which each correlation matrix was thresholded and converted to the adjacency matrix. The adjacency matrix is a means of representing which nodes of the graph are adjacent to other vertices with 1 dedicating the existing edge and 0 indicating the absence of an edge between the two nodes. A brain functional connection could be represented as an undirected edge if the correlation coefficient between the two nodes is larger than a correlation threshold corresponding to a sparsity value which was defined as the total number of edges in a network divided by the maximum possible number of edges [39]. In this way, the resulting graphs would be comprised of the same number of connections, which make cortical networks in the HC and patient groups to have the same wiring cost [39], [40]. We thresholded each correlation matrix repeatedly over a wide range of thresholds (15% to 25%) and then estimated the network properties at each threshold value. To consider the abnormal connectivity and cortical hubs between the HC and patient groups, the network pattern at threshold sparsity (S)  = 15% was shown as being typical in the following analysis. This value was the lowest threshold that could avoid generating isolate node in the networks (the cortical networks are fully connected), and it could also minimize the number of spurious edges in the networks [41], [42]. This threshold was employed in several previous GTA studies [15]–[17], [40], [43], [44].

2.5 Small-world properties analysis

The discovery of the small-world network, in which most nodes are not neighbors of one another but can be attained via a small number of steps, has offered new insights into various aspects of the network properties [45], [46]. Small-world structure, which has a highly efficient neuronal architecture, was frequently found in the brain network organization [38], [41]. Several network metrics were calculated to assess small-world properties. The key parameters of a small-world network were the clustering coefficient of a network C and the mean minimum path length L. The clustering coefficient of a node 0<Ci<1 is a ratio that defines the proportion of possible connections that actually exist between the nearest neighbors of a node [36], [37], [45]. The clustering coefficient of a network C is the average over each node's clustering coefficient:

where N is the total number of nodes in the network, is the number of edges connecting the neighbors of node j, and is the number of neighbors of node j. The minimum path length is the average of the shortest path lengths over each possible pair of vertices [36], [37], [45]:

where is the shortest path length between the node and the node, and the path length is defined as the number of edges included in the path. Corresponding parameters for a random graph of C and L with the same number of nodes were also calculated, as denoted by and . The original degree distribution is preserved when random networks are created with a random reconnection of the edges in the network of origin at least ten times to destroy the neighboring structures [47]. Random networks have a small average shortest path length, but with limited local interconnections resulting in a small and . A graph is considered small-world if its average clustering coefficient C is significantly higher than a random graph constructed on the same number of nodes, and if the graph has a small average shortest path length [36]. We examined the ratio and the ratio in our resting networks. The ratio could be summarized for small-world networks as typically being >1 [36].

Betweenness centrality

Based on anatomical evidence, Hagmann et al. (2008) pointed out that the brain region within the structural core was composed of connector hubs that linked all of the major structural modules [48]. These cortical hubs play a critical role in integrating diverse informational sources and balancing the opposing pressure to evolve segregated, specialized pathways [49]. To calculate the relative importance of a vertex in a graph, we applied the betweenness centrality measurement in our analysis. Nodes that occur on many shortest paths between other nodes have higher betweenness than those that do not [50]:

where is the betweenness centrality of a node i, is the shortest path number from node j to node k, and is the shortest path number from node j to node k that passes through node i. estimates the significance of a brain region over information flow between other brain regions in the cortical networks.

2.6 Network robustness analysis

To assess the resilience of the brain networks to targeted attack in the pathological conditions of a migraine, we used the approach described by Albert and Barabasi (2002) [51]. Targeted attack of the network was stimulated by selecting the node with the largest betweenness centrality. We then removed it and assessed the size of the remaining largest connected cluster. We repeated this process and incrementally eliminated additional nodes in decreasing order of their network betweenness centrality across all subjects [36], [42], [51]. Finally, we obtained network robustness measures in the HC and PM group respectively and estimated their differences.

2.7 Interregional correlations analysis

While graph visualization provided strong hints of gender-related dysregulated brain regions in PM, we further evaluated the abnormal interregional correlations of the resting networks in the male and female PM. The correlation coefficient between the two regions was preserved as the weights of the edges. Then, Fisher's z transformation was performed to convert the correlation coefficients to a z value [40].

2.8 Statistical analysis

To investigate whether there were significant differences in the topological properties of the resting networks between HC and PM, a two-way ANOVA was performed on each network metric (small-worldness, betweenness centrality, network robustness and interregional correlations) to model gender (male vs. female) and disease state (patients vs. controls) effects simultaneously. The test evaluated the hypothesis that gender, disease state and interaction effects were all the same against the alternative that they were not all the same. We paid more attention to the significant interactions between gender and disease state, and performed a test to determine which pairs of effects were significantly different. In addition, we performed partial correlation analyses between the changes of network organization and migraine duration across all subjects while controlling for age.

Results

3.1 Small-world property

In the current study, we constructed resting functional networks by using GTA between gender (male vs. female) and disease state (patients vs. controls). Small-world network properties were obtained at different network sparsity values (Fig. 1 A and B) from 15% to 25% in 0.01 increments. For the significant interactions, we used a multiple comparison test to determine which pairs of effects were significantly different. For males, small-world network organization was significant altered at 20%, 23% and 24% between PM and HC (Fig. 1 A, p<0.05, FDR corrected). For females, there were significant small-worldness differences to the migraine compared to HC at 15%, 16% and 17% (Fig. 1 B, p<0.05, FDR corrected). The horizontal star line in Fig.1 indicates the significant difference in the network parameters between the HC and PM groups among all of the thresholds.

[Figure omitted. See PDF.]

Figure 1. The mean and standard deviation of small-worldness from the resting networks at different sparsity values.

The colored lines show small-worldness differences among (A) male PM, (B) female PM and matched controls. The horizontal stars indicate the significant difference of the small-worldness between HC and PM (p<0.05, FDR corrected).

https://doi.org/10.1371/journal.pone.0027049.g001

The sparsity S = 15% was selected as being typical in the following network analysis in all subjects (see Materials and Methods). The results exhibited that the small-worldness in female PM was positively correlated with the duration of migraine by applying a linear partial correlation model controlling for the patients' age (Fig. 2, r = 0.53, p = 0.02). We did not find any significant correlation in male PM (p>0.05).

[Figure omitted. See PDF.]

Figure 2. Correlation between the small-worldness of female PM resting networks and the migraine duration while controlling for patients' age at threshold S = 15%.

https://doi.org/10.1371/journal.pone.0027049.g002

3.2 Network robustness

The network robustness evaluation could provide us with a precise quantitative measure of global impairment for PM. Our results showed that the resting networks in HC were resilient to targeted attack than in PM (Fig.3). Specifically, the mean size of the largest connected component in the resting networks of male HC was decreased to 60% by targeted attack when the top 33% nodes with high betweenness values were removed, but was 30% of the nodes in male PM (Fig. 3 A). For females, the brain networks in HC were disintegrated to the same extent when the top 34% nodes were deleted, but was 27% of the nodes in female PM (Fig. 3 B). Comparing topological properties of the male PM resting network to that of the female PM, gender-related differences were found in our results. As can be seen from Figure 3, the network robustness revealed different degrees of loss in male and female PM. Significant network breakdowns (p<0.05) were found in male and female PM due to targeted node attacks (horizontal dotted line in Fig.3).

[Figure omitted. See PDF.]

Figure 3. Network robustness analysis for targeted attacks.

Graph presents the mean size of the largest connected component versus the proportion of total nodes eliminated by targeted attacks among (A) male PM, (B) female PM and matched controls. The horizontal dotted line indicates the significant difference between HC and PM (p<0.05).

https://doi.org/10.1371/journal.pone.0027049.g003

3.3 Nodal centrality

Based on each nodal betweenness centrality, our aim was to assess the prominent nodal centrality differences between the HC and PM resting networks, and investigate the gender-related differences in local impairments between male and female PM networks. Five ROIs showed a significant interaction between gender (male vs. female) and disease state (patients vs. controls) from the whole brain ANOVA analysis (Fig. 4, p<0.05, FDR corrected). As shown in Figure 4, these brain regions included the precentral gyrus (PreCG), dorsolateral superior frontal gyrus (SFGdor), orbital inferior frontal gyrus (ORBinf), anterior cingulate gyrus (ACG), and parahippocampal gyrus (PHG). Furthermore, there was a significant decrease in PM females compared to HC in all of these brain regions (Table 4). In male groups, only the ORBinf exhibited a significant decrease in nodal centrality in PM as compared to HC. In addition, the nodal centrality of the PreCG and ACG presented a significant correlation with the duration of migraine (Table 4) in female PM.

[Figure omitted. See PDF.]

Figure 4. Significant interactions (nodal centrality differences) between gender (male vs. female) and disease state (patients vs. controls) effects.

A two-way ANOVA was performed (p<0.05, corrected). This figure was visualized with the BrainNet viewer.

https://doi.org/10.1371/journal.pone.0027049.g004

[Figure omitted. See PDF.]

Table 4. The nodal centrality (betweeness) across gender and disease state.

https://doi.org/10.1371/journal.pone.0027049.t004

3.4 Interregional correlations

Several pairs of connections were significantly altered in PM (p<0.05, FDR corrected) (Fig. 5 and Fig. 6). In the male PM group, significantly increased connections were found in the PFC, SMG, amygdala (AMYG), HIP, inferior parietal lobule and temporal lobes (Fig. 5). In the female PM group, the network of patients was more disorganized compared with matched controls, and more regions showed significant increases in interregional correlations (Fig. 6), including the orbital part of the PFC, posterior cingulate gyrus (PCG), PHG, cuneus (CUN), putamen (PUT), caudate (CAU), parietal lobule, temporal lobes, and occipital cortex (Fig. 6). No significantly decreased connections were found in PM.

[Figure omitted. See PDF.]

Figure 5. Significant differences in the intensity of the functional connection in male PM resting networks.

Ten connections (red lines) showed increased intensity in the male PM resting networks (migraine > HC, p<0.05, FDR corrected). This figure was visualized with the BrainNet viewer.

https://doi.org/10.1371/journal.pone.0027049.g005

[Figure omitted. See PDF.]

Figure 6. Significant differences in the intensity of the functional connections in the female PM resting networks.

Ten connections (red lines) showed increased intensity in the female PM resting networks (migraine > HC, p<0.05, FDR corrected). This figure was visualized with the BrainNet viewer.

https://doi.org/10.1371/journal.pone.0027049.g006

Discussion

In this study, we viewed resting functional networks of migraine from an integration perspective. Compared with HC, we found abnormal topological organization in PM brain networks, including the small-world property, resilience, nodal centrality, and interregional connections. Furthermore, gender-related differences of these network parameters in migraine patients were also observed.

Small-world models are useful for connectivity studies of nervous systems because they have high clustering and a short path length which confers the capability for both specialized or modular processing in local neighborhoods, as well as including distributed or integrated processing over the entire network [36]. In the present study, salient small-worldness was found in the resting networks of HC, consistent with several previous cortical network studies [52]. However, both male and female PM resting networks showed different small-world scalars over a wide range of sparsity values as compared with HC, demonstrating that more variable topological properties existed in these two groups (Fig. 1). Previous studies demonstrated that a migraine was associated with a significant gray matter reduction in several brain areas [23], [53]–[56]. While subtle brain lesions led to abnormal brain function [57]–[59], long-term headaches could alter functional brain organization of migraine patients [60]. As mentioned above, our results indicated that the small-world property was disrupted in PM, which may be the consequence of frequent migraine-related nociceptive input in patients' daily lives.

The small-worldness measurement of the resting network for HC, male and female PM could provide valuable insights into the characteristics of abnormal organization in a migraine patient's brain. It still needs further analysis to conclude the global damage for the resting network in migraine and its gender-related differences. Achard et al. (2006) suggested that healthy brain organization had advantages in terms of robustness to a selective attack on cortical hubs, and it could confer resilience against the loss of network functionality in the face of a pathological attack [36]. As compared with HC, male and female PM resting networks tended to fragment rapidly when responding to targeted attacks (Fig. 3). These findings might be interpreted as the results of long-term and high-frequency headache attacks, which give rise to lesions and to migraine-related loss of functional integrity, thus leading to unfavorable influences on the functional connections within the resting functional networks. Interestingly, we noted differences between male and female PM in their network robustness. Resting networks in female PM exhibited the worst tolerance against targeted node removal among the three groups. This indicated that the functional networks in females may be more vulnerable to disruption of the developmental aberration resulting from migraine.

Many studies have investigated the central pain network involved in nociceptive processing in the human brain, and several cortical and subcortical brain regions exhibited altered neural activity in response to pain in different experimental conditions [23], [61], [62]. In our study, there was a significant nodal centrality decrease in the PreCG, SFGdor, ORBinf, ACG, and PHG (Fig. 4 and Table 4). These brain regions were largely considered to be involved in pain processing [63]. Our results were consistent with many previous studies of migraine on brain structure and function [22], [23], [54], [63], [64]. Furthermore, we found differences between male and female PM in the distribution of abnormal nodal centrality as compared to matched controls (Table 4). While male PM nodal centrality only decreased in the orbital part of the PFC, more regions showed less nodal centrality in female PM. From our observations, it could be inferred that migraine may affect males and females differently in the organizational patterns of brain functional networks in PM. Moreover, nodal centrality is a useful measure for the relative importance of a brain region to integrate diverse global information into the resting networks [16], [17], [40]. More regions showed significantly decreased centrality in females indicating more brain damage done in female PM, which may have a strong possibility of resulting in abnormal information integration during the experience and the anticipation of pain [6].

Figure 5 and Figure 6 present abnormal interregional correlations in male and female PM. Compared with matched controls, there were many functional connections showing a significant increase in the intensity of the interregional correlations, mainly in the SI, SII, SMG, PFC, striatum, HIP, PHG, AMYG, occipital cortices, and temporal cortices. During noxious stimuli, most of these regions were significantly activated to mediate the unpleasant-affective dimension of pain and the motivation to escape from noxious events in central pain processing [4], [65]–[67]. According to previous studies [4]–[6], these dysfunctional interregional correlations observed between brain regions within the pain processing networks could be understood as a brain injury and may be the result of the secondary effect of having migraine. It is important to note that abnormal functional connectivities in the chronically migraine-inflicted brain have a skewed balance between males and females. More dysfunctional connections were found in the female PM resting networks (Fig. 6). This finding was consistent with our results above that migraine had a distinction between male and female patients, where the resting networks of female PM had more functional abnormalities than in males. Recently, Kruit et al. (2004) pointed out that women, but not men, with migraine were associated with increased risk of white matter abnormalities [68]. There exist gender differences in the underlying brain network organization [15]–[17], [69]. Studies showed that white matter in women was more efficiently utilized [15], [70]. Gong et al. (2009) utilized diffusion tensor imaging tractography to construct cortical structure networks, and that the cortical networks of women had both higher overall global and local efficiency [15]. Gur et al. (1999) found a stronger association between white matter volume and cognitive performance in women [70]. According to our results, this may have a potential correlation between migraine and female brain characteristics. It is possible that the distinct features of female cortical networks may likely suffer more brain damage in females with migraine. We inferred that women's cortical networks may have a higher risk for migraine, which may explain the different prevalence of migraine between males and females.

There is a methodological limitation to the present study. While a graph represents functional connection relationships between brain regions, GTA characterizes the functional brain network by constructing a voxel-based network with each voxel as a network node [47], [49], or by constructing a region-based network with nodes corresponding to anatomically-defined ROIs [71]. Several studies also pointed out that different parcellation strategies of graph analytical techniques may result in a considerable variation in the exact values defining key parameters of network organization [72], [73]. In our study, functional brain networks were constructed at a regional level by dividing the brain into 90 regions. We do not know if the dysfunctional network organization of PM we found will change significantly with different parcellation strategies. To test the reproducibility of our results, future studies needs to consider the effect that a specific parcellation approach has on graph analytical findings.

In summary, by exploring the topological properties of the resting brain, we investigated the brain organization of migraine that included a gender component. We observed that functional damages showed a skewed balance between males and females in their network robustness, nodal centrality, and functional connections. The resting networks of females were more vulnerable to migraine. Our findings may help us better understand gender differences in PM during the resting state. Gender-related differences should be considered when designing experiments or interpreting results in future migraine studies.

Acknowledgments

We would like to thank Qizhu Wu and Junran Zhang for valuable technical assistance in conducting this research.

Author Contributions

Conceived and designed the experiments: J. Liu WQ QG FL JT. Performed the experiments: J. Liu J. Li JN KY JS FZ LZ. Analyzed the data: J. Liu. Contributed reagents/materials/analysis tools: FZ LZ DY MD PL. Wrote the paper: J. Liu KMD.

Citation: Liu J, Qin W, Nan J, Li J, Yuan K, Zhao L, et al. (2011) Gender-Related Differences in the Dysfunctional Resting Networks of Migraine Suffers. PLoS ONE6(11): e27049. https://doi.org/10.1371/journal.pone.0027049

References

1. Lipton RB, Hamelsky SW, Kolodner KB, Steiner TJ, Stewart WF (2000) Migraine, quality of life, and depression: a population-based case-control study. Neurology 55: 629–635.RB LiptonSW HamelskyKB KolodnerTJ SteinerWF Stewart2000Migraine, quality of life, and depression: a population-based case-control study.Neurology55629635

2. Schwedt TJ, Dodick DW (2009) Advanced neuroimaging of migraine. Lancet Neurol 8: 560–568.TJ SchwedtDW Dodick2009Advanced neuroimaging of migraine.Lancet Neurol8560568

3. Terwindt GM, Ferrari MD, Tijhuis M, Groenen SMA, Picavet HSJ, et al. (2000) The impact of migraine on quality of life in the general population. Neurology 55: 624–629.GM TerwindtMD FerrariM. TijhuisSMA GroenenHSJ Picavet2000The impact of migraine on quality of life in the general population.Neurology55624629

4. Chiapparini L, Ferraro S, Grazzi L, Bussone G (2010) Neuroimaging in chronic migraine. Neurol Sci 31: 19–22.L. ChiappariniS. FerraroL. GrazziG. Bussone2010Neuroimaging in chronic migraine.Neurol Sci311922

5. Filippi M, Rocca MA (2008) Headache and migraine. Neurol Sci 29: 336–338.M. FilippiMA Rocca2008Headache and migraine.Neurol Sci29336338

6. May A (2009) New insights into headache: an update on functional and structural imaging findings. Nat Rev Neurol 5: 199–209.A. May2009New insights into headache: an update on functional and structural imaging findings.Nat Rev Neurol5199209

7. Matharu MS, Bartsch T, Ward N, Frackowiak RSJ, Weiner R, et al. (2004) Central neuromodulation in chronic migraine patients with suboccipital stimulators: a PET study. Brain 127: 220–230.MS MatharuT. BartschN. WardRSJ FrackowiakR. Weiner2004Central neuromodulation in chronic migraine patients with suboccipital stimulators: a PET study.Brain127220230

8. Cao Y, Aurora SK, Nagesh V, Patel SC, Welch KMA (2002) Functional MRI-BOLD of brainstem structures during visually triggered migraine. Neurology 59: 72–78.Y. CaoSK AuroraV. NageshSC PatelKMA Welch2002Functional MRI-BOLD of brainstem structures during visually triggered migraine.Neurology597278

9. Welch KMA, Cao Y, Aurora S, Wiggins G, Vikingstad EM (1998) MRI of the occipital cortex, red nucleus, and substantia nigra during visual aura of migraine. Neurology 51: 1465–1469.KMA WelchY. CaoS. AuroraG. WigginsEM Vikingstad1998MRI of the occipital cortex, red nucleus, and substantia nigra during visual aura of migraine.Neurology5114651469

10. Afridi SK, Giffin NJ, Kaube H, Friston KJ, Ward NS, et al. (2005) A positron emission tomographic study in spontaneous migraine. Arch Neurol 62: 1270–1275.SK AfridiNJ GiffinH. KaubeKJ FristonNS Ward2005A positron emission tomographic study in spontaneous migraine.Arch Neurol6212701275

11. Afridi SK, Matharu MS, Lee L, Kaube H, Friston KJ, et al. (2005) A PET study exploring the laterality of brainstem activation in migraine using glyceryl trinitrate. Brain 128: 932–939.SK AfridiMS MatharuL. LeeH. KaubeKJ Friston2005A PET study exploring the laterality of brainstem activation in migraine using glyceryl trinitrate.Brain128932939

12. Stewart WF, Wood C, Reed ML, Roy J, Lipton RB (2008) Cumulative lifetime migraine incidence in women and men. Cephalalgia 28: 1170–1178.WF StewartC. WoodML ReedJ. RoyRB Lipton2008Cumulative lifetime migraine incidence in women and men.Cephalalgia2811701178

13. Stovner LJ, Zwart JA, Hagen K, Terwindt GM, Pascual J (2006) Epidemiology of headache in Europe. Eur J Neurol 13: 333–345.LJ StovnerJA ZwartK. HagenGM TerwindtJ. Pascual2006Epidemiology of headache in Europe.Eur J Neurol13333345

14. Sowell ER, Peterson BS, Kan E, Woods RP, Yoshii J, et al. (2007) Sex differences in cortical thickness mapped in 176 healthy individuals between 7 and 87 years of age. Cereb Cortex 17: 1550–1560.ER SowellBS PetersonE. KanRP WoodsJ. Yoshii2007Sex differences in cortical thickness mapped in 176 healthy individuals between 7 and 87 years of age.Cereb Cortex1715501560

15. Gong G, Rosa-Neto P, Carbonell F, Chen ZJ, He Y, et al. (2009) Age-and gender-related differences in the cortical anatomical network. J Neurosci 29: 15684–15693.G. GongP. Rosa-NetoF. CarbonellZJ ChenY. He2009Age-and gender-related differences in the cortical anatomical network.J Neurosci291568415693

16. Yan C, Gong G, Wang J, Wang D, Liu D, 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.C. YanG. GongJ. WangD. WangD. Liu2011Sex-and brain size–related small-world structural cortical networks in young adults: a DTI tractography study.Cereb Cortex21449458

17. Tian L, Wang J, Yan C, He Y (2011) Hemisphere-and gender-related differences in small-world brain networks: A resting-state functional MRI study. NeuroImage 54: 191–202.L. TianJ. WangC. YanY. He2011Hemisphere-and gender-related differences in small-world brain networks: A resting-state functional MRI study.NeuroImage54191202

18. Hamilton C (2009) Cognition and sex differences. 51: Basingstoke,UK: Palgrave Macmillan. 965–965.C. Hamilton2009Cognition and sex differences.51Basingstoke,UKPalgrave Macmillan965965

19. Shepard RN, Metzler J (1971) Mental rotation of three-dimensional objects. Science 171: 701–703.RN ShepardJ. Metzler1971Mental rotation of three-dimensional objects.Science171701703

20. McClure EB (2000) A meta-analytic review of sex differences in facial expression processing and their development in infants, children, and adolescents. Psychol Bull 126: 424–453.EB McClure2000A meta-analytic review of sex differences in facial expression processing and their development in infants, children, and adolescents.Psychol Bull126424453

21. Baron-Cohen S, O'Riordan M, Stone V, Jones R, Plaisted K (1999) Recognition of faux pas by normally developing children and children with Asperger syndrome or high-functioning autism. J Autism Dev Disord 29: 407–418.S. Baron-CohenM. O'RiordanV. StoneR. JonesK. Plaisted1999Recognition of faux pas by normally developing children and children with Asperger syndrome or high-functioning autism.J Autism Dev Disord29407418

22. Schmitz N, Admiraal-Behloul F, Arkink EB, Kruit MC, Schoonman GG, et al. (2008) Attack frequency and disease duration as indicators for brain damage in migraine. Headache 48: 1044–1055.N. SchmitzF. Admiraal-BehloulEB ArkinkMC KruitGG Schoonman2008Attack frequency and disease duration as indicators for brain damage in migraine.Headache4810441055

23. Valfrè W, Rainero I, Bergui M, Pinessi L (2008) Voxel-based morphometry reveals gray matter abnormalities in migraine. Headache 48: 109–117.W. ValfrèI. RaineroM. BerguiL. Pinessi2008Voxel-based morphometry reveals gray matter abnormalities in migraine.Headache48109117

24. Carter AR, Astafiev SV, Lang CE, Connor LT, Rengachary J, et al. (2009) Resting inter-hemispheric fMRI connectivity predicts performance after stroke. Ann Neurol 67: 365–375.AR CarterSV AstafievCE LangLT ConnorJ. Rengachary2009Resting inter-hemispheric fMRI connectivity predicts performance after stroke.Ann Neurol67365375

25. Fox MD, Raichle ME (2007) Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci 8: 700–711.MD FoxME Raichle2007Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging.Nat Rev Neurosci8700711

26. Sporns O, Tononi G (2002) Classes of network connectivity and dynamics. Complexity 7: 28–38.O. SpornsG. Tononi2002Classes of network connectivity and dynamics.Complexity72838

27. Fox MD, Greicius M (2010) Clinical applications of resting state functional connectivity. Fron Syst Neurosci 4: 19.MD FoxM. Greicius2010Clinical applications of resting state functional connectivity.Fron Syst Neurosci419

28. He Y, Evans A (2010) Graph theoretical modeling of brain connectivity. Curr Opin Neurol 23: 341–350.Y. HeA. Evans2010Graph theoretical modeling of brain connectivity.Curr Opin Neurol23341350

29. Wang J, Zuo X, He Y (2010) Graph-based network analysis of resting-state functional MRI. Frontiers Syst Neurosci 4: 16.J. WangX. ZuoY. He2010Graph-based network analysis of resting-state functional MRI.Frontiers Syst Neurosci416

30. Bullmore E, Bassett DS (2011) Brain graphs: graphical models of the human brain connectome. Annu Rev Clin Psychol 7: 113–140.E. BullmoreDS Bassett2011Brain graphs: graphical models of the human brain connectome.Annu Rev Clin Psychol7113140

31. Bullmore E, Sporns O (2009) Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci 10: 186–198.E. BullmoreO. Sporns2009Complex brain networks: graph theoretical analysis of structural and functional systems.Nat Rev Neurosci10186198

32. Detsky ME, McDonald DR, Baerlocher MO, Tomlinson GA, McCrory DC, et al. (2006) Does this patient with headache have a migraine or need neuroimaging? JAMA 296: 1274–1283.ME DetskyDR McDonaldMO BaerlocherGA TomlinsonDC McCrory2006Does this patient with headache have a migraine or need neuroimaging?JAMA29612741283

33. Friston KJ, Ashburner J, Frith CD, Poline JB, Heather JD, et al. (1995) Spatial registration and normalization of images. Hum Brain Mapp 3: 165–189.KJ FristonJ. AshburnerCD FrithJB PolineJD Heather1995Spatial registration and normalization of images.Hum Brain Mapp3165189

34. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, et al. (2002) Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15: 273–289.N. Tzourio-MazoyerB. LandeauD. PapathanassiouF. CrivelloO. Etard2002Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain.Neuroimage15273289

35. Liu Y, Liang M, Zhou Y, He Y, Hao Y, et al. (2008) Disrupted small-world networks in schizophrenia. Brain 131: 945–961.Y. LiuM. LiangY. ZhouY. HeY. Hao2008Disrupted small-world networks in schizophrenia.Brain131945961

36. Achard S, Salvador R, Whitcher B, Suckling J, Bullmore E (2006) A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J Neurosci 26: 63–72.S. AchardR. SalvadorB. WhitcherJ. SucklingE. Bullmore2006A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs.J Neurosci266372

37. Salvador R, Suckling J, Coleman MR, Pickard JD, Menon D, et al. (2005) Neurophysiological architecture of functional magnetic resonance images of human brain. Cereb Cortex 15: 1332–1342.R. SalvadorJ. SucklingMR ColemanJD PickardD. Menon2005Neurophysiological architecture of functional magnetic resonance images of human brain.Cereb Cortex1513321342

38. He Y, Chen ZJ, Evans AC (2007) Small-world anatomical networks in the human brain revealed by cortical thickness from MRI. Cereb Cortex 17: 2407–2419.Y. HeZJ ChenAC Evans2007Small-world anatomical networks in the human brain revealed by cortical thickness from MRI.Cereb Cortex1724072419

39. Achard S, Bullmore E (2007) Efficiency and cost of economical brain functional networks. PLoS Comput Biol 3: e17.S. AchardE. Bullmore2007Efficiency and cost of economical brain functional networks.PLoS Comput Biol3e17

40. Yao Z, Zhang Y, Lin L, Zhou Y, Xu C, et al. (2010) Abnormal Cortical Networks in Mild Cognitive Impairment and Alzheimer's Disease. PLoS Comput Biol 6: e1001006.Z. YaoY. ZhangL. LinY. ZhouC. Xu2010Abnormal Cortical Networks in Mild Cognitive Impairment and Alzheimer's Disease.PLoS Comput Biol6e1001006

41. He Y, Chen Z, Evans A (2008) Structural insights into aberrant topological patterns of large-scale cortical networks in Alzheimer's disease. J Neurosci 28: 4756–4766.Y. HeZ. ChenA. Evans2008Structural insights into aberrant topological patterns of large-scale cortical networks in Alzheimer's disease.J Neurosci2847564766

42. Bernhardt BC, Chen Z, He Y, Evans AC, Bernasconi N (2011) Graph-Theoretical Analysis Reveals Disrupted Small-World Organization of Cortical Thickness Correlation Networks in Temporal Lobe Epilepsy. Cereb Cortex: In Press. BC BernhardtZ. ChenY. HeAC EvansN. Bernasconi2011Graph-Theoretical Analysis Reveals Disrupted Small-World Organization of Cortical Thickness Correlation Networks in Temporal Lobe Epilepsy.Cereb Cortex: In Press

43. Gong G, He Y, Concha L, Lebel C, Gross DW, et al. (2009) Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography. Cereb Cortex 19: 524–536.G. GongY. HeL. ConchaC. LebelDW Gross2009Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography.Cereb Cortex19524536

44. Lo CY, Wang PN, Chou KH, Wang J, He Y, et al. (2010) Diffusion tensor tractography reveals abnormal topological organization in structural cortical networks in Alzheimer's disease. J Neurosci 30: 16876–16885.CY LoPN WangKH ChouJ. WangY. He2010Diffusion tensor tractography reveals abnormal topological organization in structural cortical networks in Alzheimer's disease.J Neurosci301687616885

45. Watts DJ, Strogatz SH (1998) Collective dynamics of ‘small-world’networks. Nature 393: 440–442.DJ WattsSH Strogatz1998Collective dynamics of ‘small-world’networks.Nature393440442

46. Reijneveld JC, Ponten SC, Berendse HW, Stam CJ (2007) The application of graph theoretical analysis to complex networks in the brain. Clin Neurophysiol 118: 2317–2331.JC ReijneveldSC PontenHW BerendseCJ Stam2007The application of graph theoretical analysis to complex networks in the brain.Clin Neurophysiol11823172331

47. Hayasaka S, Laurienti PJ (2010) Comparison of characteristics between region-and voxel-based network analyses in resting-state fMRI data. Neuroimage 50: 499–508.S. HayasakaPJ Laurienti2010Comparison of characteristics between region-and voxel-based network analyses in resting-state fMRI data.Neuroimage50499508

48. Hagmann P, Cammoun L, Gigandet X, Meuli R, Honey CJ, et al. (2008) Mapping the structural core of human cerebral cortex. PLoS Biol 6: 1479–1493.P. HagmannL. CammounX. GigandetR. MeuliCJ Honey2008Mapping the structural core of human cerebral cortex.PLoS Biol614791493

49. Buckner RL, Sepulcre J, Talukdar T, Krienen FM, Liu H, 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.RL BucknerJ. SepulcreT. TalukdarFM KrienenH. Liu2009Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease.J Neurosci2918601873

50. Freeman LC (1977) A set of measures of centrality based on betweenness. Sociometry 40: 35–41.LC Freeman1977A set of measures of centrality based on betweenness.Sociometry403541

51. Albert R, Barabási AL (2002) Statistical mechanics of complex networks. Rev mod Phys 74: 47–97.R. AlbertAL Barabási2002Statistical mechanics of complex networks.Rev mod Phys744797

52. Bassett DS, Bullmore E (2006) Small-world brain networks. Neuroscientist 12: 512–523.DS BassettE. Bullmore2006Small-world brain networks.Neuroscientist12512523

53. May A, Gaser C (2006) Magnetic resonance-based morphometry: a window into structural plasticity of the brain. Curr Opin Neurol 19: 407–411.A. MayC. Gaser2006Magnetic resonance-based morphometry: a window into structural plasticity of the brain.Curr Opin Neurol19407411

54. Rocca MA, Ceccarelli A, Falini A, Colombo B, Tortorella P, et al. (2006) Brain gray matter changes in migraine patients with T2-visible lesions: a 3-T MRI study. Stroke 37: 1765–1770.MA RoccaA. CeccarelliA. FaliniB. ColomboP. Tortorella2006Brain gray matter changes in migraine patients with T2-visible lesions: a 3-T MRI study.Stroke3717651770

55. Schmidt-Wilcke T, Leinisch E, Straube A, K mpfe N, Draganski B, et al. (2005) Gray matter decrease in patients with chronic tension type headache. Neurology 65: 1483–1486.T. Schmidt-WilckeE. LeinischA. StraubeN. K mpfeB. Draganski2005Gray matter decrease in patients with chronic tension type headache.Neurology6514831486

56. May A (2008) Chronic pain may change the structure of the brain. Pain 137: 7–15.A. May2008Chronic pain may change the structure of the brain.Pain137715

57. Salat DH, Kaye JA, Janowsky JS (2002) Greater orbital prefrontal volume selectively predicts worse working memory performance in older adults. Cereb Cortex 12: 494–505.DH SalatJA KayeJS Janowsky2002Greater orbital prefrontal volume selectively predicts worse working memory performance in older adults.Cereb Cortex12494505

58. Driscoll I, Hamilton DA, Petropoulos H, Yeo RA, Brooks WM, et al. (2003) The aging hippocampus: cognitive, biochemical and structural findings. Cereb Cortex 13: 1344–1351.I. DriscollDA HamiltonH. PetropoulosRA YeoWM Brooks2003The aging hippocampus: cognitive, biochemical and structural findings.Cereb Cortex1313441351

59. Wen W, Zhu W, He Y, Kochan NA, Reppermund S, et al. (2011) Discrete neuroanatomical networks are associated with specific cognitive abilities in old age. J Neurosci 31: 1204–1212.W. WenW. ZhuY. HeNA KochanS. Reppermund2011Discrete neuroanatomical networks are associated with specific cognitive abilities in old age.J Neurosci3112041212

60. Baliki MN, Geha PY, Apkarian AV, Chialvo DR (2008) Beyond feeling: chronic pain hurts the brain, disrupting the default-mode network dynamics. J Neurosci 28: 13981403.MN BalikiPY GehaAV ApkarianDR Chialvo2008Beyond feeling: chronic pain hurts the brain, disrupting the default-mode network dynamics.J Neurosci2813981403

61. Apkarian AV, Bushnell MC, Treede RD, Zubieta JK (2005) Human brain mechanisms of pain perception and regulation in health and disease. Eur J Pain 9: 463–484.AV ApkarianMC BushnellRD TreedeJK Zubieta2005Human brain mechanisms of pain perception and regulation in health and disease.Eur J Pain9463484

62. Tracey I (2005) Nociceptive processing in the human brain. Curr Opin Neurobiol 15: 478–487.I. Tracey2005Nociceptive processing in the human brain.Curr Opin Neurobiol15478487

63. Kim JH, Suh SI, Seol HY, Oh K, Seo WK, et al. (2008) Regional grey matter changes in patients with migraine: a voxel-based morphometry study. Cephalalgia 28: 598–604.JH KimSI SuhHY SeolK. OhWK Seo2008Regional grey matter changes in patients with migraine: a voxel-based morphometry study.Cephalalgia28598604

64. Schmidt-Wilcke T, G n bauer S, Neuner T, Bogdahn U, May A (2008) Subtle grey matter changes between migraine patients and healthy controls. Cephalalgia 28: 1–4.T. Schmidt-WilckeS. G n bauerT. NeunerU. BogdahnA. May2008Subtle grey matter changes between migraine patients and healthy controls.Cephalalgia2814

65. Do Hyung Kang JHS, Kim YC (2010) Neuroimaging Studies of Chronic Pain. Kor J Pain 23: 159–165.JHS Do Hyung KangYC Kim2010Neuroimaging Studies of Chronic Pain.Kor J Pain23159165

66. Price DD (2000) Psychological and neural mechanisms of the affective dimension of pain. Science 288: 17691772.DD Price2000Psychological and neural mechanisms of the affective dimension of pain.Science28817691772

67. Treede RD, Kenshalo DR, Gracely RH, Jones AKP (1999) The cortical representation of pain. Pain 79: 105–111.RD TreedeDR KenshaloRH GracelyAKP Jones1999The cortical representation of pain.Pain79105111

68. Kruit MC, van Buchem MA, Hofman PAM, Bakkers JTN, Terwindt GM, et al. (2004) Migraine as a risk factor for subclinical brain lesions. JAMA 291: 427–434.MC KruitMA van BuchemPAM HofmanJTN BakkersGM Terwindt2004Migraine as a risk factor for subclinical brain lesions.JAMA291427434

69. Gong G, He Y, Evans AC (2011) Brain Connectivity: Gender Makes a Difference. Neuroscientist: In Press. G. GongY. HeAC Evans2011Brain Connectivity: Gender Makes a Difference.Neuroscientist: In Press

70. Gur RC, Turetsky BI, Matsui M, Yan M, Bilker W, et al. (1999) Sex differences in brain gray and white matter in healthy young adults: correlations with cognitive performance. J Neurosci 19: 4065–4072.RC GurBI TuretskyM. MatsuiM. YanW. Bilker1999Sex differences in brain gray and white matter in healthy young adults: correlations with cognitive performance.J Neurosci1940654072

71. Liu J, Liang J, Qin W, Tian J, Yuan K, et al. (2009) Dysfunctional connectivity patterns in chronic heroin users: An fMRI study. Neurosci Lett 460: 72–77.J. LiuJ. LiangW. QinJ. TianK. Yuan2009Dysfunctional connectivity patterns in chronic heroin users: An fMRI study.Neurosci Lett4607277

72. Zalesky A, Fornito A, Harding IH, Cocchi L, Yücel M, et al. (2010) Whole-brain anatomical networks: Does the choice of nodes matter? Neuroimage 50: 970–983.A. ZaleskyA. FornitoIH HardingL. CocchiM. Yücel2010Whole-brain anatomical networks: Does the choice of nodes matter?Neuroimage50970983

73. Fornito A, Zalesky A, Bullmore ET (2010) Network scaling effects in graph analytic studies of human resting-state FMRI data. Fron Syst Neurosci 4: 22.A. FornitoA. ZaleskyET Bullmore2010Network scaling effects in graph analytic studies of human resting-state FMRI data.Fron Syst Neurosci422