Changes of Functional Brain Networks in Major Depressive Disorder: A Graph Theoretical Analysis of Resting-State fMRI

Abstract

Recent developments in graph theory have heightened the need for investigating the disruptions in the topological structure of functional brain network in major depressive disorder (MDD). In this study, we employed resting-state functional magnetic resonance imaging (fMRI) and graph theory to examine the whole-brain functional networks among 42 MDD patients and 42 healthy controls. Our results showed that compared with healthy controls, MDD patients showed higher local efficiency and modularity. Furthermore, MDD patients showed altered nodal centralities of many brain regions, including hippocampus, temporal cortex, anterior cingulate gyrus and dorsolateral prefrontal gyrus, mainly located in default mode network and cognitive control network. Together, our results suggested that MDD was associated with disruptions in the topological structure of functional brain networks, and provided new insights concerning the pathophysiological mechanisms of MDD.

About the Authors:

Ming Ye

Contributed equally to this work with: Ming Ye, Tianliang Yang

Affiliations Faculty of Computer and Information Science, Southwest University, Chongqing, China, College of Electronic and Information Engineering, Southwest University, Chongqing, China

Tianliang Yang

Contributed equally to this work with: Ming Ye, Tianliang Yang

Affiliations Faculty of psychology, Southwest University, Chongqing, China, Key Laboratory of Cognition and Personality (Ministry of Education), Chongqing, China

Peng Qing

Affiliation: Faculty of Computer and Information Science, Southwest University, Chongqing, China

Xu Lei

* E-mail: [email protected] (LGY); [email protected] (LX)

Affiliations Faculty of psychology, Southwest University, Chongqing, China, Key Laboratory of Cognition and Personality (Ministry of Education), Chongqing, China

Jiang Qiu

Affiliations Faculty of psychology, Southwest University, Chongqing, China, Key Laboratory of Cognition and Personality (Ministry of Education), Chongqing, China

Guangyuan Liu

* E-mail: [email protected] (LGY); [email protected] (LX)

Affiliation: College of Electronic and Information Engineering, Southwest University, Chongqing, China

Introduction

Major depressive disorder (MDD) is one of the most prevalent mental disorders, and it is characterized by persistent feeling of sadness, low self-esteem, sleep disturbances, and withdrawal from pleasurable activities [1]. Although much efforts have been made in the treatments of MDD, up to 80% of patients still suffer from a relapse [2].

Previous studies suggest that MDD not only has relationships with the regional deficits, but also with the abnormal functional integration of distributed brain regions [3,4,5,6,7]. Using typical regional measures including regional homogeneity and the (fractional) amplitude of characteristic low-frequency fluctuations, a number of brain regions with abnormal activities in the resting-state have been identified associated with MDD, such as parahippocampal gyrus [4,5], prefrontal cortex [8,9,10], cingulate gyrus [11,12], fusiform gyrus [5,11], and thalamus [13]. Moreover, disruptions in functional connectivity have been observed between specific region pairs in MDD through functional connectivity analyses (including seed-based connectivity, independent component analyses and network homogeneity), such as increased medial prefrontal cortex and anterior cingulate cortex connectivity [6], increased subgenual cingulate-thalamic connectivity [3,14], and reduced bilateral dorsal lateral prefrontal cortex and right superior parietal lobule connectivity [15]. In addition, several neural networks have been observed to mediate depressive disorders based on functional connectivity analyses, such as the default mode network (DMN) and cognitive control network (CCN) [7,8,15,16].

Despite the above-mentioned advances, recent studies have investigated the MDD patients in their whole-brain networks using graph theory. Graph theory is a mathematical technique which can model the brain as a complex network represented graphically by sets of nodes and edges [17,18]. Based on graph theory, numerous studies of healthy subjects have concluded that human brain networks have many special organizational principles, including small-worldness (an optimal brain network organization characterized by high efficiency of information transfer at a low cost), and modularity (an optimal partition of a brain network into smaller functional communities of modules) [19,20]. Furthermore, disruptions in network organization have been related to many neuropsychiatric disorders. For patients with MDD, several studies have reported topological changes in human brain connectome, including a loss of the small-world network [19] and a significant reorganization of the community structure [14,21,22]. However, these studies were limited in examining the brain network organization of MDD patients at meso-level, which could be described by the modularity of a network [14,19,21,22]. In addition, the application of graph theory on MDD has provided conflicting results. Two functional neuroimaging studies found that MDD resulted in decreased path length, but no change in clustering coefficient [14,21]. However, another study showed a prominent changes of the community structure, but no significant differences in both path length and clustering coefficient [22]. These inconsistent findings may result from differences in methodology, study population (Chinese vs German cohort), and variability in the clinical definition of MDD.

In the present study, we employed resting-state fMRI and graph theoretical analyses to further explore the MDD-related disruptions in functional brain networks. Given previous evidence of abnormal functional connectivity in MDD, combined with findings of disrupted network structure, we predicted that MDD disrupted the topological organization of human brain connectome, such as significant anomalies in small-worldness and modularity. In addition, considering the importance of DMN and CCN in mediating MDD, we also hypothesized that the DMN and CCN module would show significant disruptions in MDD patients. Finally, the relationships between inter-group differences in topology properties and individual clinical variables (Hamilton Depression Scale, HAMD) were also investigated.

Materials and Methods

Subjects

Eighty-four subjects were recruited, including 42 MDD patients and 42 sex-, age-, and education-matched healthy controls (Table 1). The age of MDD patients and healthy controls both ranged from 18 to 60 years. The dataset reported here was randomly selected from our ongoing project, which examined the occurrence and development of depression. The diagnosis of MDD was made according to the Structured Clinical Interview of the DSM-IV by experienced psychiatrists from the First Affiliated Hospital of Chongqing Medical University [23]. All healthy controls were carefully screened for a current or lifetime diagnosis of any Axis I or II disorder using the Structured Clinical Interview of the DSM-IV Non-Patient Edition and Structural Clinical Interview for DSM-IV Axis II Personality Disorder. Organic or neurologic disorders were examined based on personal histories and complete physical examinations. The severity of depression was measured through the 17-item Hamilton Rating Scale for Depression (HAMD). To be suitable for this study, all patients were re-examined by a psychiatry expert after an initial outpatient assessment. Inclusion criteria for all patients were (1) patients met the judgment of MDD defined by DSM-IV; (2) 17-item HAMD scores were larger than 24 [24]; (3) patients were medication- naïve; and (4) diagnosed as depressive disorder, not bipolar disorder. Exclusion criteria included the presence of (1) other Axis I psychiatric disorders and symptoms; (2) a history of organic brain disorder, neurological disorders or cardiovascular diseases; (3) pregnancy or any physical illness as assessed by personal history and laboratory analysis; and (4) the inability to undergo an fMRI scan. Healthy controls met the exclusion criteria above and had no psychiatric illness history or any family history of major psychiatric disease in their first-degree relatives. This study was approved by the Ethics Committee of Southwest University, and the written informed consent was obtained for each subject.

[Figure omitted. See PDF.]

Table 1. Demographic and clinical characteristics of the study samples.

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

Image Acquisition and Preprocessing

A total of 242 volumes of resting-state functional images were obtained for each subject using an echo planar imaging (EPI) sequence through a 3T Siemens Trio scanner (TR/TE = 2000/30ms, flip angle = 90°, acquisition matrix = 64×64, field of view = 220×220mm², axial slices = 32, and thickness/gap = 3/1mm).

Functional data preprocessing was carried out using SPM8 (http://www.fil.ion.ucl.ac.uk/spm). The entire process included removal of the first 10 volumes, slice timing correction, realignment to the first volume for head-motion correction, filtering (0.01~0.08 Hz), normalization to the EPI template with a resampling voxel size of 3×3×3 mm³, smoothing with a 6mm full-width at half-maximum Gaussian kernel. No subjects were excluded because all the head motions were <2 mm or 2°.

Construction of Functional Brain Network

To construct a functional brain network, we firstly employed the automated anatomical labeling (AAL) template [25] to parcellate the brain into 90 regions of interest (ROIs). Secondly, the time series was acquired on each ROI by averaging the signals of all voxels within that area and then linearly regressing out the influences of head motion and global signal. Thirdly, by calculating the Pearson correlation coefficients in the residual time courses between each pair of ROIs, a correlation matrix was obtained for each subject. To improve the normality, the correlation map was Fisher transformed (r-z) [26]. Finally, the absolute z values were converted into a binary connection matrix to make a graphic model of a brain network. That is, if the absolute z_ij (Fisher r-to-z of the Pearson correlation coefficient) of a pair of brain regions, i and j, exceeded a given threshold T (Fisher r-to-z), an edge was said to exist; otherwise it did not exist.

The degree of each node, D_nod, is the number of connections that link it to the rest of the network. The total number of edges in a network, divided by the maximum possible number of edges (N²-N)/2, is called the connection density or the cost of the network [27]. Given that there was no accurate way to choose a threshold in studies of brain networks [28], so the functional brain networks were constructed over the whole value of costs (0.03~0.50) at the interval of 0.01. Because a similar trend for between-group differences was observed over the range of 0.03~0.50, and the biggest difference between MDD patients and healthy controls was found when cost is 0.21 for global measures, only results using a cost of 0.21 were reported for regional nodal analyses [29,30].

Functional Brain Network Analysis

For functional brain networks, both the global network parameters and the regional nodal parameters were calculated to characterize the global topological organization of functional network and regional properties of each node [18,31,32,33]. The global network parameters include global efficiency (E_global), local efficiency (E_local) and modularity (Q). The regional nodal parameters include nodal degree (D_nod), nodal efficiency (E_nod), and nodal betweenness (N_b).

Global Network Parameters

For a given graph G with N nodes, the global efficiency (E_global) is defined by the inverse of the harmonic mean of the minimum path length between each pair of nodes:where L_i,j is the shortest path length from node i to node j.

The local efficiency (E_local) of a graph G is defined as:where E_global(i) is the global efficiency of G_i. In particular, E_global measures the overall communication efficiency of the network, and E_local measures the local cliquishness of the network.

A brain network can be considered as a small-world network, as it meets the following criteria [19,34]: E_global(G_regular) < E_global(G_real) < E_global(G_random) and E_local(G_random) < E_local(G_real) < E_local(G_regular). E_global(G_regular), E_global(G_real) and E_global(G_random) represent global efficiencies of regular networks, real networks and random networks, respectively. And E_local(G_regular), E_local(G_real) and E_local(G_random) represent the local efficiencies of regular networks, real networks and random networks, respectively.

Modularity refers to the formation of the local modules that nodes in the same module are closely connected to each other while nodes in different modules are sparsely connected. For the study of the distribution of network modules by group, mean group functional matrices were calculated by averaging the N×N (N = 90 in the present study) absolute connection matrix of all the subjects within the group (MDD group and healthy group) [29]. There are several algorithms which are used in quantifying the partition in terms of module separation. For example, how well a partition differentiates subsets of nodes densely connected [32,35]. The modularity, Q, of a graph G can be quantified as the proportion of G’s edges that fall within modules, subtracted by the proportion that would be expected due to random chance alone, for a given partition of nodes into module [31,32]. Q can be defined as:where m is the total number of edges; A_ij = 1 if an edge links i and j, and 0 otherwise; δ(M_i, M_j) is 1 if i and j are in the same module and 0 otherwise, and ensures that only intra-modular edges are added to the sum; P_ij is the probability that there would be an edge between i and j, given a random graph comparable to G. The value of P_ij could be estimated by:where K_i is the total number of edges connecting node i. We include this information in the null model because it affects the expected proportion of intra-modular edges.

Nodal Centrality Parameters

The nodal degree of a node i is defined as:where e_ij is the (i,j)th element in the formerly obtained binary, undirected network.

The nodal efficiency of a node i is defined as:

The nodal betweenness of a node i is defined as:where δ_jk is the number of the shortest path length between node j and node k, and δ_jk(i) is the number of the shortest path length between node j and node k, which pass through node i. N_b measures the influence of node i on the flow of information between other nodes in the network.

Statistical Analysis

Statistical comparisons of topological properties (both the global network parameters and the regional nodal parameters) between MDD patients and healthy controls were performed by using the two-sample two-tailed t-tests. Moreover, the relationships between topological measurements and HAMD scores in MDD patients were assessed using Pearson Correlation. In addition, the false discovery rate (FDR) was used for multiple comparison corrections [36].

Results

Altered Small-world properties in MDD Patients

In this study, we constructed functional brain networks on the global scale for MDD patients and healthy controls. Graph theory analyses revealed that the global efficiencies of functional brain networks in MDD patients and healthy controls were greater than regular networks but less than random networks, and the local efficiencies were greater than random networks but less than regular networks over the whole range of 0.03~0.50 (Fig 1). These results were typical features of small-worldness and were compatible with previous studies of small-world brain networks [14,21].

[Figure omitted. See PDF.]

Fig 1. Small-world properties of functional brain networks.

(A) global efficiency and (B) local efficiency over the whole range of 0.03~0.50 for random (green), regular (blue), and real networks (MDD patients: red; Healthy controls: black). Error bars corresponded to standard error of the mean. Purple stars indicated where the difference between MDD patients and healthy controls was significant (p<0.05). On average, both healthy controls and MDD patients showed small-world properties. Additionally, MDD patients showed higher local efficiencies at the range of 0.15~0.30 (p<0.05).

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

Despite the common small-world properties, higher local efficiencies were found in MDD patients in comparison with healthy controls at the range of 0.15~0.30, whereas there were no significant differences in global efficiencies (Fig 1).

Altered Modularity in MDD Patients

Functional brain networks of MDD patients and healthy controls exhibited typical features of modular structure. Specifically, functional brain networks for both groups were significantly more modular than random networks with the same degree distribution over the range of 0.03~0.37 (Fig 2). In addition, a two-sample two-tailed t-test revealed that MDD patients showed increased modularity compared with healthy controls over the range of 0.14~0.22 (Fig 2). The brain networks were decomposed into 5 basic modules in healthy controls, but 6 basic modules in MDD patients (Fig 3).

[Figure omitted. See PDF.]

Fig 2. Modularity of functional brain networks.

Functional brain networks of MDD patients (red) and healthy controls (black) showed larger modularity than random networks (green) at the whole range of 0.03~0.37. Additionally, increased modularity in MDD patients was observed over the range of 0.14~0.22 (p<0.05). Error bars corresponded to standard error of the mean. Purple stars indicated where the difference between MDD patients and healthy controls was significant (p<0.05).

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

[Figure omitted. See PDF.]

Fig 3. The modular organization of group averaged functional brain network at cost = 0.21.

(A) Healthy controls had five modules. (B) MDD patients had six modules. Module 1~Module 5 were the same network for healthy controls and MDD patients, but Module 6 was a new one for MDD patients.

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

Altered Regional Nodal Properties in MDD Patients

Significant differences were found on nearly the same nodes (Fig 4, Table 2), except for right gyrus rectus and right thalamus whose nodal betweenness showed no significant differences between MDD patients and healthy controls (p<0.05, FDR corrected). Specifically, MDD patients revealed increased nodal centralities in many brain areas, including right gyrus rectus, right hippocampus, bilateral amygdala, right fusiform gyrus, bilateral middle temporal gyrus, and bilateral thalamus compared with healthy controls. Additionally, decreased nodal centralities in MDD patients were observed in bilateral dorsolateral prefrontal gyrus and bilateral anterior cingulate gyrus.

[Figure omitted. See PDF.]

Fig 4. Regions exhibited significant between-group differences in regional nodal parameters.

The blue color represented the higher values of regional nodal centralities in healthy controls, and the red color represented the higher values of regional nodal centralities in MDD patients (p<0.05, FDR corrected).

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

[Figure omitted. See PDF.]

Table 2. Regions showing disrupted nodal centralities in MDD patients as compared with healthy controls (cost = 0.21).

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

Relationships between Network Parameters and HAMD

We found no significant correlations between global network metrics (E_global and E_local) and HAMD scores in MDD patients. However, significant positive correlation between modularity (Q) and HAMD scores was observed (r = 0.36, p = 0.018, cost = 0.21; Fig 5). Additionally, there were significant correlations between regional nodal parameters and HAMD scores among many brain regions (Table 3). Specifically, nodal centralities were significantly and negatively correlated with HAMD scores in bilateral dorsolateral frontal gyrus and bilateral anterior cingulate gyrus. Furthermore, significant positive correlations were observed between HAMD scores and nodal centralities in right hippocampus, bilateral amygdala, right fusiform gyrus, bilateral thalamus, and bilateral middle temporal gyrus.

[Figure omitted. See PDF.]

Fig 5. Scatter plot with trend line showed the relationship between modularity and HAMD scores for cost = 0.21 in MDD patients.

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

[Figure omitted. See PDF.]

Table 3. Pearson correlation coefficients between regional nodal properties and HAMD scores of MDD patients (cost = 0.21).

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

Discussion

In this paper, we investigated the topological architecture of functional brain networks in MDD patients. The main findings were as follows: 1) functional brain networks of MDD patients showed increased local efficiency and modularity; 2) many local brain regions, mainly located in default mode network (DMN) and cognitive control network (CCN), were significantly affected by MDD: increased nodal centralities were observed in right gyrus rectus, right hippocampus, bilateral amygdala, right fusiform gyrus, bilateral middle temporal gyrus and bilateral thalamus, whereas decreased nodal centralities were found in bilateral dorsolateral prefrontal gyrus and bilateral anterior cingulate gyrus. In sum, these significant findings could expand our understanding of neurophysiologic mechanisms related to MDD from a network perspective.

Disrupted Small-world Properties in MDD Patients

Our results demonstrated that functional brain networks of both MDD patients and healthy controls showed prominent small-world property, which was in line with previous studies on MDD [14,21]. Networks with small-world properties ensure a higher information-processing efficiency for both locally specialized and globally integrated processing [17]. Despite the common small-world properties, the local efficiencies were significantly higher in MDD patients than healthy controls. Local efficiency is the measure of local network connectivity, so the increase of local efficiency in MDD may represent disrupted information processing among distant brain areas [37,38,39,40]. Contrasted with our results, one fMRI study showed a significantly increased global efficiency, but unchanged local efficiency[14]. The discrepancies may be attributed to the differences in the MDD populations. Zhang et al. (2011) recruited MDD patients with mean HAMD scores≧18, whereas our study recruited MDD patients with HAMD≧24. Patients at different stages may manifest different symptoms with distinct neuronal correlates [41].

MDD-related Changes of Nodal Characteristics

MDD-related increases in nodal centralities were mainly located in right hippocampus, bilateral amygdala, right fusiform gyrus, bilateral middle temporal gyrus and bilateral thalamus, most of which were components of DMN [15,42]. Increasing evidence indicated that abnormal activities in hippocampus, amygdala, fusiform gyrus, middle temporal gyrus, and thalamus were associated with disruptions in emotional functions [14,43,44]. Specifically, thalamus were thought to be involved in emotional perception, and amygdala, fusiform gyrus and middle temporal gyrus played key roles in the neural responses to negative stimuli [45,46,47,48,49]. Disruptions in these regions may result in negative bias in interpersonal feedback and somatic complaints [7,50]. In addition, hippocampus disruptions were also closely associated with negative bias in MDD patients [38,45]. Disruptions in hippocampus would increase memory sensitivity to negative stimuli in MDD patients [38,45]. Furthermore, we also found that MDD disrupted the functioning of the gyrus rectus, a critical region involved in the regulation of mood and cognition [47]. MDD patients showed abnormal brain activities in the gyrus rectus during specific tasks or in a resting state [48,49]. Particularly, Fitzgerald et al. (2008) showed that MDD patients showed increased activations in right gyrus rectus when exposed to positive stimuli. And more notably, we also found that there were significant positive correlations between the nodal centralities and HAMD scores in right hippocampus, bilateral amygdala, right fusiform gyrus, bilateral middle temporal gyrus and bilateral thalamus. This indicated that if the severity of depression state increased across patients, the nodal centralities of these regions would become higher. Moreover, because all of these regions were key regions implicated in the pathophysiology of MDD, the properties of these regions could predict the depressive state to some extent.

MDD-related decreases in nodal centralities were mainly found in bilateral anterior cingulate gyrus and bilateral dorsolateral prefrontal cortex, which were closely related to cognitive control network [15,16]. Abnormalities in anterior cingulate gyrus and dorsolateral prefrontal cortex were not surprising since they were both closely associated with emotion regulation [48,51]. The dorsolateral prefrontal gyrus has been suggested a key role in the cognitive control functioning, including allocation of attention [52,53]. Disruptions in this region may lead to more focused attention on negative aspects of one’s self in MDD patients. Furthermore, the anterior cingulate gyrus played an important role in the generation of negative mood states [43,54]. Anomalies of this region were associated with high emotional involvement [54,55]. In addition, significant negative correlations were observed between the nodal centralities and HAMD scores in bilateral anterior cingulate gyrus and bilateral dorsolateral prefrontal cortex. This implied that as the disease progresses, the emotion regulation in MDD patients would become worse. Moreover, all of these results further strengthened our understanding of the pathophysiological mechanism of MDD.

Disruptions in Modularity

Resting-state functional brain network has been shown as a modular organization [56]. The highly modularized structures in MDD patients and healthy controls of the current study provided further evidence. Key circuits associated with main brain functions were consistently observed in previous studies, such as cognitive control network (CCN), auditory system, visual system, default mode network (DMN) and subcortical system [32,57,58]. The current results were compatible with these findings: Module 1 was mainly related to cognitive control; Module 2 was primarily involved in subcortical regions; Module 3 and Module 5 were mainly associated with auditory and visual functions. Although Module 4 primarily involved in DMN was also identified, there were different brain regions in healthy controls and MDD patients. Additionally, a new module (Module 6) associated with cognitive control function was observed in MDD patients. From above analyses, we found that the CCN and DMN showed disruptions in MDD patients. The DMN has been shown to be associated with self-referential processing, so disruptions in DMN may be associated with the increases of recall and rehearsal of negative life events among MDD patients [42,59]. And the CCN has been postulated to play an important role in patients’ cognitive control of emotion regulation deficit [43,60]. Therefore, the DMN and CCN may represent important neural substrates of MDD. In addition, increased modularity in MDD patients indicated that there were relatively less inter-modular edges and more intra-modular edges, which may also be associated with the disruptions in emotion regulation by decreasing communications between the DMN and CCN [16]. Different from our results, although Lord et al. (2012) observed a significant reorganization of the community structure in MDD, they did not find the differences in the modularity. The inconsistencies may be due to more severely depressed subjects in our study (their mean HAMD scores were 15.8 compared with 24.5 in our study) [22]. In addition, culture differences (Chinese vs German populations) may also affect the classification of depression.

Limitations and Further Considerations

There are several issues needed to be addressed. Firstly, our results of increased local efficiencies in MDD patients were different from one fMRI study which showed comparable local efficiencies between MDD patients and healthy controls [22]. Variations in clinical characteristics of MDD may be one important factor to account for these discrepancies. Different clinical features and depression severity may result in different brain activations patterns and further cause different brain network topological properties [41,61]. Future studies, which use MDD patients with different clinical characteristics, may give us a more complete understanding of brain abnormalities. Secondly, the clinical diagnosis of MDD has some limitations and is easily affected by subjective knowledge and experience. Hence, the objective diagnosis of MDD has high clinical value in preventing severe disease. In addition, the present study suggested that the measurement of topological properties was a preferential candidate for diagnosing MDD. However, different topological properties were found between functional and structural brain networks [14,40]. Thus, investigating the topological properties of human connectome in combination with functional and structural neuroimaging is also important. Moreover, the ability to diagnose MDD could be further enhanced using multi-modality MRI. Thirdly, large studies have indicated that the cerebellum was closely associated with higher-order functions, including emotion regulation and cognitive processing, and have suggested that the cerebellum should be included in the pathophysiological models of MDD [62,63]. However, because the AAL template we chose did not contain the cerebellum, the current study analyzed the functional brain network without cerebellum. In future, the cerebellum should be involved in the analyses by selecting the more comprehensive template or making a special analysis on the cerebellum. Finally, the application of graph theory to investigate network disruptions in MDD patients is still in its preliminary stage, and there are still many issues to be deepening. Until now, MDD patients have been investigated widely using graph theory. Next, we should form a hypothesis-generating framework to understand the relationships between abnormal topological properties and MDD.

Conclusion

Using graph theory, we found that the MDD patients showed increased local efficiency and modularity. Additionally, brain regions mainly implicated in emotional and cognitive function were significantly affected by MDD. These results provide new and important insights into the neural correlates of MDD. Further work could be conducted to examine how the topological structure of functional brain networks is altered with different clinical depressive symptoms and severity.

Acknowledgments

The work was conducted in the Chongqing Medical University. We thank Dr. Peng Xie for support in data acquisition.

Author Contributions

Conceived and designed the experiments: MY PQ. Performed the experiments: MY PQ JQ GYL. Analyzed the data: MY PQ TLY. Contributed reagents/materials/analysis tools: MY PQ. Wrote the paper: MY TLY.

Citation: Ye M, Yang T, Qing P, Lei X, Qiu J, Liu G (2015) Changes of Functional Brain Networks in Major Depressive Disorder: A Graph Theoretical Analysis of Resting-State fMRI. PLoS ONE 10(9): e0133775. https://doi.org/10.1371/journal.pone.0133775

References

1. Belmaker R, Agam G. Major depressive disorder. New England Journal of Medicine. 2008;358:55–68. pmid:18172175

2. Gotlib IH, Hamilton JP. Neuroimaging and Depression Current Status and Unresolved Issues. Current Directions in Psychological Science. 2008;17:159–163.

3. Greicius MD, Flores BH, Menon V, Glover GH, Solvason HB, Kenna H, et al. Resting-state functional connectivity in major depression: abnormally increased contributions from subgenual cingulate cortex and thalamus. Biological psychiatry. 2007;62:429–437. pmid:17210143

4. Liu F, Guo W, Liu L, Long Z, Ma C, Xue Z, et al. Abnormal amplitude low-frequency oscillations in medication-naive, first-episode patients with major depressive disorder: a resting-state fMRI study. Journal of affective disorders. 2013;146:401–406. pmid:23116810

5. Wu QZ, Li DM, Kuang WH, Zhang TJ, Lui S, Huang XQ, et al. Abnormal regional spontaneous neural activity in treatment‐refractory depression revealed by resting‐state fMRI. Human brain mapping. 2011;32:1290–1299. pmid:20665717

6. Zhu X, Wang X, Xiao J, Liao J, Zhong M, Wang W, et al. Evidence of a dissociation pattern in resting-state default mode network connectivity in first-episode, treatment-naive major depression patients. Biological psychiatry. 2012;71:611–617. pmid:22177602

7. Wang L, Hermens D, Hickie I, Lagopoulos J. A systematic review of resting-state functional-MRI studies in major depression. Journal of affective disorders. 2012;142:6–12. pmid:22858266

8. Guo W, Liu F, Zhang J, Zhang Z, Yu L, Liu J, et al. Abnormal default-mode network homogeneity in first-episode, drug-naive major depressive disorder. PloS one. 2014;9:e91102. pmid:24609111

9. Liu C-H, Ma X, Li F, Wang Y-J, Tie C-L, Li S-F, et al. Regional homogeneity within the default mode network in bipolar depression: a resting-state functional magnetic resonance imaging study. PloS one. 2012;7:e48181. pmid:23133615

10. Wang L, Dai W, Su Y, Wang G, Tan Y, Jin Z, et al. Amplitude of low-frequency oscillations in first-episode, treatment-naive patients with major depressive disorder: a resting-state functional MRI study. PloS one. 2012;7:e48658. pmid:23119084

11. Guo W-b, Liu F, Chen J-d, Gao K, Xue Z-m, Xu X-j, et al. Abnormal neural activity of brain regions in treatment-resistant and treatment-sensitive major depressive disorder: a resting-state fMRI study. Journal of psychiatric research. 2012;46:1366–1373. pmid:22835912

12. Yao Z, Wang L, Lu Q, Liu H, Teng G. Regional homogeneity in depression and its relationship with separate depressive symptom clusters: a resting-state fMRI study. Journal of affective disorders. 2009;115:430–438. pmid:19007997

13. Wang L, Li K, Zhang Q, Zeng Y, Dai W, Su Y, et al. Short-term effects of escitalopram on regional brain function in first-episode drug-naive patients with major depressive disorder assessed by resting-state functional magnetic resonance imaging. Psychological medicine. 2014;44:1417–1426. pmid:23942213

14. Zhang J, Wang J, Wu Q, Kuang W, Huang X, He Y, et al. Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder. Biological psychiatry. 2011;70:334–342. pmid:21791259

15. Wei X, Shen H, Ren J, Li X, Xu X, Yang R, et al. Altered Resting-State Connectivity in College Students with Nonclinical Depressive Symptoms. PloS one. 2014;9:e114603. pmid:25502215

16. Sheline YI, Price JL, Yan Z, Mintun MA. Resting-state functional MRI in depression unmasks increased connectivity between networks via the dorsal nexus. Proceedings of the National Academy of Sciences. 2010;107:11020–11025.

17. Bullmore E, Sporns O. Complex brain networks: graph theoretical analysis of structural and functional systems. Nature Reviews Neuroscience. 2009;10:186–198. pmid:19190637

18. Wang J, Zuo X, He Y. Graph-based network analysis of resting-state functional MRI. Frontiers in systems neuroscience. 2010;4.

19. Achard S, Bullmore E. Efficiency and cost of economical brain functional networks. PLoS computational biology. 2007;3:e17. pmid:17274684

20. He Y, Evans A. Graph theoretical modeling of brain connectivity. Current opinion in neurology. 2010;23:341–350. pmid:20581686

21. Leistedt SJ, Coumans N, Dumont M, Lanquart JP, Stam CJ, Linkowski P. Altered sleep brain functional connectivity in acutely depressed patients. Human brain mapping. 2009;30:2207–2219. pmid:18937282

22. Lord A, Horn D, Breakspear M, Walter M. Changes in community structure of resting state functional connectivity in unipolar depression. PloS one. 2012;7:e41282. pmid:22916105

23. First MB, Spitzer RL, Gibbon M, Williams JB. User's guide for the Structured clinical interview for DSM-IV axis I disorders SCID-I: clinician version: American Psychiatric Pub; 1997.

24. Williams JB. A structured interview guide for the Hamilton Depression Rating Scale. Archives of general psychiatry. 1988;45:742–747. pmid:3395203

25. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage. 2002;15:273–289. pmid:11771995

26. Liu Y, Liang M, Zhou Y, He Y, Hao Y, Song M, et al. Disrupted small-world networks in schizophrenia. Brain. 2008;131:945–961. pmid:18299296

27. Latora V, Marchiori M. Economic small-world behavior in weighted networks. The European Physical Journal B-Condensed Matter and Complex Systems. 2003;32:249–263.

28. Bassett DS, Bullmore E. Small-world brain networks. The neuroscientist. 2006;12:512–523. pmid:17079517

29. Baggio HC, Sala‐Llonch R, Segura B, Marti MJ, Valldeoriola F, Compta Y, et al. Functional brain networks and cognitive deficits in Parkinson's disease. Human brain mapping. 2014;35:4620–4634. pmid:24639411

30. Bullmore ET, Bassett DS. Brain graphs: graphical models of the human brain connectome. Annual review of clinical psychology. 2011;7:113–140. pmid:21128784

31. Alexander-Bloch AF, Gogtay N, Meunier D, Birn R, Clasen L, Lalonde F, et al. Disrupted modularity and local connectivity of brain functional networks in childhood-onset schizophrenia. Frontiers in systems neuroscience. 2010;4:147. pmid:21031030

32. Sun Y, Yin Q, Fang R, Yan X, Wang Y, Bezerianos A, et al. Disrupted Functional Brain Connectivity and Its Association to Structural Connectivity in Amnestic Mild Cognitive Impairment and Alzheimer’s Disease. PloS one. 2014;9:e96505. pmid:24806295

33. Rubinov M, Sporns O. Complex network measures of brain connectivity: uses and interpretations. Neuroimage. 2010;52:1059–1069. pmid:19819337

34. Latora V, Marchiori M. Efficient behavior of small-world networks. Physical review letters. 2001;87:198701. pmid:11690461

35. Tagliazucchi E, Von Wegner F, Morzelewski A, Brodbeck V, Borisov S, Jahnke K, et al. Large-scale brain functional modularity is reflected in slow electroencephalographic rhythms across the human non-rapid eye movement sleep cycle. Neuroimage. 2013;70:327–339. pmid:23313420

36. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological). 1995;57:289–300.

37. Hamilton JP, Furman DJ, Gotlib IH. The neural foundations of major depression: Classical approaches and new frontiers. Neurobiology of depression. 2011:57–73.

38. Hamilton JP, Siemer M, Gotlib IH. Amygdala volume in major depressive disorder: a meta-analysis of magnetic resonance imaging studies. Molecular psychiatry. 2008;13:993–1000. pmid:18504424

39. Sacher J, Neumann J, Fünfstück T, Soliman A, Villringer A, Schroeter ML. Mapping the depressed brain: a meta-analysis of structural and functional alterations in major depressive disorder. Journal of affective disorders. 2012;140:142–148. pmid:21890211

40. Singh MK, Kesler SR, Hadi Hosseini S, Kelley RG, Amatya D, Hamilton JP, et al. Anomalous gray matter structural networks in major depressive disorder. Biological psychiatry. 2013;74:777–785. pmid:23601854

41. Sahay A, Hen R. Adult hippocampal neurogenesis in depression. Nature neuroscience. 2007;10:1110–1115. pmid:17726477

42. Sheline YI, Barch DM, Price JL, Rundle MM, Vaishnavi SN, Snyder AZ, et al. The default mode network and self-referential processes in depression. Proceedings of the National Academy of Sciences. 2009;106:1942–1947.

43. Bora E, Fornito A, Pantelis C, Yücel M. Gray matter abnormalities in major depressive disorder: a meta-analysis of voxel based morphometry studies. Journal of affective disorders. 2012;138:9–18. pmid:21511342

44. Hamilton JP, Etkin A, Furman DJ, Lemus MG, Johnson RF, Gotlib IH. Functional neuroimaging of major depressive disorder: a meta-analysis and new integration of baseline activation and neural response data. American Journal of Psychiatry. 2012;169:693–703. pmid:22535198

45. Disner SG, Beevers CG, Haigh EA, Beck AT. Neural mechanisms of the cognitive model of depression. Nature Reviews Neuroscience. 2011;12:467–477. pmid:21731066

46. Hamilton JP, Gotlib IH. Neural substrates of increased memory sensitivity for negative stimuli in major depression. Biological psychiatry. 2008;63:1155–1162. pmid:18281017

47. Ballmaier M, Toga AW, Blanton RE, Sowell ER, Lavretsky H, Peterson J, et al. Anterior cingulate, gyrus rectus, and orbitofrontal abnormalities in elderly depressed patients: an MRI-based parcellation of the prefrontal cortex. American Journal of Psychiatry. 2004;161:99–108. pmid:14702257

48. Fitzgerald PB, Laird AR, Maller J, Daskalakis ZJ. A meta‐analytic study of changes in brain activation in depression. Human brain mapping. 2008;29:683–695. pmid:17598168

49. Lisiecka D, Meisenzahl E, Scheuerecker J, Schoepf V, Whitty P, Chaney A, et al. Neural correlates of treatment outcome in major depression. The International Journal of Neuropsychopharmacology. 2011;14:521–534. pmid:21205435

50. Cullen KR, Gee DG, Klimes-Dougan B, Gabbay V, Hulvershorn L, Mueller BA, et al. A preliminary study of functional connectivity in comorbid adolescent depression. Neuroscience letters. 2009;460:227–231. pmid:19446602

51. Drevets WC, Price JL, Bardgett ME, Reich T, Todd RD, Raichle ME. Glucose metabolism in the amygdala in depression: relationship to diagnostic subtype and plasma cortisol levels. Pharmacology Biochemistry and Behavior. 2002;71:431–447.

52. Jin C, Gao C, Chen C, Ma S, Netra R, Wang Y, et al. A preliminary study of the dysregulation of the resting networks in first-episode medication-naive adolescent depression. Neuroscience letters. 2011;503:105–109. pmid:21871534

53. Koechlin E, Summerfield C. An information theoretical approach to prefrontal executive function. Trends in cognitive sciences. 2007;11:229–235. pmid:17475536

54. Mayberg HS, Liotti M, Brannan SK, McGinnis S, Mahurin RK, Jerabek PA, et al. Reciprocal limbic-cortical function and negative mood: converging PET findings in depression and normal sadness. 2014;156:675–682.

55. Ebert D, Ebmeier KP. The role of the cingulate gyrus in depression: from functional anatomy to neurochemistry. Biological psychiatry. 1996;39:1044–1050. pmid:8780840

56. Meunier D, Lambiotte R, Fornito A, Ersche KD, Bullmore ET. Hierarchical modularity in human brain functional networks. Frontiers in neuroinformatics. 2009;3:37. pmid:19949480

57. He Y, Wang J, Wang L, Chen ZJ, Yan C, Yang H, et al. Uncovering intrinsic modular organization of spontaneous brain activity in humans. PloS one. 2009;4:e5226. pmid:19381298

58. Meunier D, Achard S, Morcom A, Bullmore E. Age-related changes in modular organization of human brain functional networks. NeuroImage. 2009;44:715–723. pmid:19027073

59. Hamilton JP, Furman DJ, Chang C, Thomason ME, Dennis E, Gotlib IH. Default-mode and task-positive network activity in major depressive disorder: implications for adaptive and maladaptive rumination. Biological psychiatry. 2011;70:327–33. pmid:21459364

60. Drevets WC. Neuroimaging and neuropathological studies of depression: implications for the cognitive-emotional features of mood disorders. Current opinion in neurobiology. 2001;11:240–249. pmid:11301246

61. Duman RS, Nakagawa S, Malberg J. Regulation of adult neurogenesis by antidepressant treatment. Neuropsychopharmacology. 2001;25:836–844. pmid:11750177

62. Guo W, Liu F, Liu J, Yu M, Zhang Z, Liu G, et al. Increased cerebellar-default-mode-network connectivity in drug-naive major depressive disorder at rest. Medicine. 2015;94:e560. pmid:25738471

63. Liu L, Zeng L-L, Li Y, Ma Q, Li B, Shen H, et al. Altered cerebellar functional connectivity with intrinsic connectivity networks in adults with major depressive disorder. PloS one. 2012;7:e39516. pmid:22724025