Original Research ARTICLE
The Abnormality of Topological Asymmetry in Hemispheric Brain Anatomical Networks in Bipolar Disorder
- 1College of Information and Computer, Taiyuan University of Technology, Taiyuan, China
- 2Faculty of Human Health Science, Graduate School of Medicine, Kyoto University, Kyoto, Japan
- 3Translational Medicine Research Center, Shanxi Medical University, Taiyuan, China
Convergent evidences have demonstrated a variety of regional abnormalities of asymmetry in bipolar disorder (BD). However, little is known about the alterations in hemispheric topological asymmetries. In this study, we used diffusion tensor imaging to construct the hemispheric brain anatomical network of 49 patients with BD and 61 matched normal controls. Graph theory was then applied to quantify topological properties of the hemispheric networks. Although small-world properties were preserved in the hemispheric networks of BD, the degrees of the asymmetry in global efficiency, characteristic path length, and small-world property were significantly decreased. More changes in topological properties of the right hemisphere than those of left hemisphere were found in patients compared with normal controls. Consistent with such changes, the nodal efficiency in patients with BD also showed less rightward asymmetry mainly in the frontal, occipital, parietal, and temporal lobes. In contrast to leftward asymmetry, significant rightward asymmetry was found in supplementary motor area of BD, and attributed to more deficits in nodal efficiency of the left hemisphere. Finally, these asymmetry score of nodal efficiency in the inferior parietal lobule and rolandic operculum were significantly associated with symptom severity of BD. Our results suggested that abnormal hemispheric asymmetries in brain anatomical networks were associated with aberrant neurodevelopment, and providing insights into the potential neural biomarkers of BD by measuring the topological asymmetry in hemispheric brain anatomical networks.
The human brain is structurally and functionally asymmetrical or lateralized (Watkins et al., 2001; Toga and Thompson, 2003). Even subtle perturbations to anatomical asymmetries between two hemispheres, such as gray matter volume (Sigalovsky et al., 2006), cortical thickness (Rimol et al., 2010), or white matter (WM) integrity (Park et al., 2004), can cause disturbances in cognitive and emotion processes. Studies have shown that aberrant brain region asymmetries are highly correlated with disturbed functions such as executive function (Yin et al., 2013), emotion (Schulte et al., 2012), and language (O’Donoghue et al., 2017). Hence, abnormal anatomical asymmetries have been observed in a variety of neurological and psychiatric disorders, including schizophrenia (Highley et al., 1999), depression (Kwon et al., 1996), attention-deficit hyperactivity disorder (Miller et al., 2006), and Bipolar disorder (BD) (Beyer et al., 2005; Bruno et al., 2008). As one of the most distinct syndromes in psychiatry, BD is mainly characterized as episodic elevations in emotion and disturbances in cognition (Belmaker, 2004; Saunders and Goodwin, 2010). Convergent studies on BD have showed abnormal asymmetries in WM (Bruno et al., 2008; Kafantaris et al., 2009; Wessa et al., 2009; Mahon et al., 2013), for example, less WM volume within the left frontal lobes, the rightward WM in orbital frontal. These results indicated that the alteration in WM asymmetries have been proposed as a key factor in the manifestation of BD symptoms.
As an imaging method, diffusion tensor imaging (DTI) can reconstruct the major WM tracts faithfully (Dae-Jin et al., 2013) and has been proved to be a promising tool for assessing WM abnormalities. Recently, using the DTI tractography and graph theory, connectome studies demonstrated abnormal topological properties in BD (Leow et al., 2013; Forde et al., 2015; Collin et al., 2016, 2017; Spielberg et al., 2016; O’Donoghue et al., 2017). The patients with BD showed decrease in global integration (longer characteristic path length, smaller global efficiency), increase in functional segregation (larger clustering coefficient and local efficiency), and loss of small-world property (the balance between local integration and functional segregation) (Leow et al., 2013; Forde et al., 2015; Collin et al., 2016; Spielberg et al., 2016; O’Donoghue et al., 2017). These studies have mainly focused on the WM topological properties in whole network rather than hemispheric network. Analyzing the hemispheric anatomical networks and further determination of the status of the anatomical network asymmetries might benefit the understanding of the underlying nature of alteration in the brain of BD, and potentially help to elucidate the etiology of the disorder. However, the hemispheric asymmetries of anatomical network in patients with BD were remained unclear.
In present study, we adopted the DTI deterministic tractography method and graph theory to investigate the abnormality of hemispheric asymmetries in brain anatomical networks in BD. In particular, we focused on global graph measures, including small-world property, global and local efficiency, and regional parameters to evaluate (1) the abnormal hemispheric asymmetries in brain anatomical networks in patients with BD and (2) whether the abnormal hemispheric asymmetries in network organization were related to clinical features of BD.
Materials and Methods
Data were selected from the UCLA Consortium for Neuropsychiatric Phenomics LA5c Study, and the study was approved by the UCLA Institutional Review Board. The data were obtained via a public database, openfMRI (Poldrack and Gorgolewski, 2015). About 49 BD patients and 61 age- and gender-matched normal subjects were selected for further analyzing. All subjects were right-handed. There are more details available in openfMRI1 (ds000030). The detailed demographics and clinical features of the patients with BD and normal controls are described in Table 1. Patient symptoms were evaluated using the 17-item Hamilton Depression Rating Scale (HAMD) (Hamilton, 1960) and the Young Mania Rating Scale (YMRS) (Young et al., 1978).
Data Acquisition and Preprocessing
Structural MRI data were acquired on 3T Siemens Trio scanners located at the Ahmanson-Lovelace Brain Mapping Center (Siemens version syngo MRB15) and the Staglin Center for Cognitive Neuroscience (Siemens version syngo MRB17) at UCLA. A high-resolution 3D echoplanar imaging was acquired with the following parameters: TR = 1.9 s, TE = 2.26 ms, FOV = 250 mm, matrix = 256 × 256, sagittal plane, slice thickness = 1 mm, 176 slices. Diffusion weighted imaging (DWI) data were collected using an echo-planar sequence with parameters: 64 directions, 2 mm slices, TR/TE = 9000/93 ms, 1 average, 96 × 96 matrix, 90° flip angle, axial slices, b = 1000 s/mm2.
This study used the MATLAB toolbox named PANDA to perform data preprocessing and the construction of the brain network (Cui et al., 2013). Specifically, data preprocessing approaches included correction for simple head movements and eddy current distortions using affine transformation to the b0 image (Jenkinson et al., 2002). After data preprocessing, the seven independent components of the diffusion tensor were estimated and from which fractional anisotropy (FA, a DTI measurement) was calculated. Subsequently, the deterministic fiber tracking algorithms were used to reconstruct fiber paths (Mori et al., 1999). The fiber tracking procedure started from the deep WM regions and terminated if two consecutive moving directions had a crossing angle above 35° or the FA was out of the threshold range (0.1∼1).
Network Construction and Analysis
In this study, the method of constructing the WM network was described in Gong et al. (Gong et al., 2009). Based on the automated anatomical labeling (AAL) parcellation scheme (Tzouriomazoyer et al., 2002), the brain was divided into 90 regions (45 in each hemisphere). Each region was defined as one node in the anatomical network. A linear transformation was applied locally within each subject’s DTI image correlated with the T1-weighted image to coregister them to the b0 image with DTI space followed by applying a nonlinear transformation to map to the ICBM152 T1 template [Montreal Neurological Institute (MNI)]. Then, the subject-specific AAL mask was weaved from the MNI space to the DTI native space with the corresponding inverse transformation, such that separate labeling values were maintained via nearest-neighbor interpolation (Gong et al., 2009; Cui et al., 2013). The FA between two regions was defined as the network edge (Shu et al., 2011; Bai et al., 2012). Prior to constructing the network, the connection between these two regions was adopted if the fiber number (FN) between the two regions was larger than 3 (Shu et al., 2011). It was helpful to reduce the influence of pseudo-connections owing to possible noise effects on whole-brain tractography. For each subject, we eliminated the inter-hemispheric connections and then obtained two weighted 45 × 45 hemispheric brain networks, one for the left hemisphere and the other for the right hemisphere.
The network architecture was then investigated at both global and regional levels for the constructed WM networks. The small-world property suggests the architecture of networks with higher local clustering and equivalent characteristic path length compared with the random network (Watts and Strogatz, 1998). In this work, eight network properties were used to analyze the topological organization of the WM networks. The clustering coefficient (Cp) of a network is thus the average of clustering coefficients across nodes and is a measure of functional segregation. The characteristic path length (Lp) of a network is the average shortest path length between all pairs of nodes in the network and is the most commonly used measure of global integration. The normalized clustering coefficients (γ), γ = Cp/Crand and the normalized characteristic path length (λ), λ = Lp/Lrand, Crand and Lrand represent indices derived from matched random networks (100 matched random networks were selected). The small-world property of a network can be characterized by both γ and λ, indicating a balance between integration and segregation. In a small-world network, the Cp is significantly higher than that of random networks (γ greater than 1), while the Lp is comparable to random networks (λ close to 1). The global efficiency (Eg), reflecting the efficiency of whole network information transmission (integration), is the inverse relation of Lp. The local efficiency (Eloc), reflected the efficiency of the network segregation. Regional properties were described in terms of nodal efficiency, Enodal(i) (Achard and Bullmore, 2007). It measures the information transmission ability of node i in the network. A node with high Enodal(i) indicates great interconnectivity with other regions in the network. The detailed descriptions of these graph theory parameters can be found in a previous work (Bullmore and Sporns, 2009; Rubinov and Sporns, 2010). The graph theory analysis was performed with GRETNA2 software.
The degree of hemispheric network asymmetry was estimated by the asymmetry score (Iturriamedina et al., 2011; Ratnarajah et al., 2013) using the following formulation: AS(X) = 100 × [X(R) − X(L)]/[X(R) + X(L)] where X(R) and X(L) represent the network properties of the right and left hemisphere, respectively. The asymmetry score AS(X) helps us look at the differences between the right and left hemisphere. Notably, a positive value of AS(X) represents a rightward asymmetry, while a negative value of AS(X) indicates a leftward asymmetry.
All statistical analyses were performed using SPSS 19 software (SPSS, Inc., Chicago, IL, United States). To determine whether there was any significant group difference in age and education, this work performed the separate two-tailed t-tests to analyze data. We used a χ2 test to analyze the gender data. To assess the group differences in hemispheric network properties, we used a general linear model (GLM) that was performed with hemisphere as a within-subject factor, group as a between-subject factor, and a group-by-hemisphere interaction. Studies have shown that hemispheric asymmetries are related to age (Cabeza, 2002; Dolcos et al., 2002) and gender (McGlone, 1980; Good et al., 2001). Hence, we set the age, gender, and age-by-gender interaction as covariates in the GLM. Further post hoc tests, including two-sample t-test for group differences and paired t-test for hemisphere differences, were used if any difference survived a threshold of P < 0.05. To determine whether the AS of the network properties within each group was significantly different from zero, one-sample t-test was performed on the asymmetry score. Meanwhile, a univariate analysis with covariance of age, gender, and age-by-gender interaction was performed on the AS to assess group differences. The threshold of P < 0.05 was considered to be significant for global properties. In particular, for the regional properties, the Bonferroni-correction was performed on that threshold (P < 0.05).
This work also studied the relationships between the hemispheric asymmetry scores and symptom severity of patients with BD. Considering age, gender, and age-by-gender as the covariates, we used both Pearson and Spearman correlation to analyze the relationship between network properties and BD symptom severity. Pearson is used to measure the linear relationship between two consecutive variables. Spearman does not require the distribution of primitive variables and belongs to a non-parametric statistical method. It was unnecessary to correct multiple comparisons because the aforementioned analyses were exploratory in nature. Hence, a significant relationship was considered at an uncorrected P-value of 0.05.
Global Properties of Hemispheric Networks
Group and Hemispheric Differences
As shown in Figure 1, all four hemispheric brain anatomical networks (2 hemispheres × 2 groups) exhibited prominent features of small-world property, as expressed by a larger γ (γ > 1) and a smaller λ (λ ≈ 1). Statistical analysis results showed significantly topological changes in the global properties in both groups and hemispheres. Except for local efficiency Eloc (F1,106 = 2.372, P = 0.127), the rest of the six network properties exhibited significant group differences between normal controls and patients with BD. There was a reduced global integration (decreased global efficiency and increased characteristic path length) and increased small-world property in patients with BD. Furthermore, we observed a significant group-by-hemisphere interaction on global efficiency Eg (F1,106 = 9.311, P = 0.003), the characteristic path length Lp (F1,106 = 7.323, P = 0.008), the normalized clustering coefficients γ (F1,106 = 30.107, P < 0.001) and the small-world property σ (F1,106 = 32.163, P < 0.001). Post hoc analysis indicated that this interaction resulted from a significant rightward trend in global integration (P = 0.001) and a significant leftward trend in the characteristic path length (P = 0.003), the normalized clustering coefficients (P < 0.001), and small-world architecture (P < 0.001) in normal controls and a symmetrical trend in patients with BD.
FIGURE 1. Further statistical analysis of network properties. Bars represent mean ± SD. Each horizontal line and associated number represent the P-value of a t-test (paired t-test for hemispheric difference at each group and two-sample t-test for group difference). Significant differences are marked with asterisk (∗∗∗P < 0.001). BD, bipolar disorder; NC, normal control.
The asymmetry score was helpful for us to evaluate the differences between the right and left hemisphere for the network properties. Additionally, the group differences in asymmetry scores would directly reflect the abnormal of hemispheric lateralization of topological organization in BD, and supplemented the result of the group-by-hemisphere interaction. Table 2 summarizes the statistical analysis results of the asymmetry scores of the global network properties for the two groups. Significant differences in hemispheric asymmetry were only observed in normal controls and disappeared in patients with BD. Normal controls showed more globally efficient in the right hemisphere than the left hemisphere [AS(Eg), t52 = 3.852, P < 0.001]. Additionally, the characteristic path length [AS(Lp), t52 = -3.852, P < 0.001], the normalized clustering coefficients [AS(γ), t52 = -6.447, P < 0.001], and the small-world property [AS(σ), t52 = -6.578, P < 0.001] showed leftward hemispheric asymmetries in normal controls. When comparing the asymmetry scores between two groups, we observed significant group differences in asymmetry scores of Eg (F1,106 = 8.268, P = 0.005), Lp (F1,106 = 8.268, P = 0.005), γ (F1,106 = 33.684, P < 0.001), and σ (F1,106 = 35.657, P < 0.001), which agreed with the significant group-by-hemisphere interaction on Eg, Lp, γ, and σ revealed by the GLM model. This result indicated that the rightward asymmetries of global integration and leftward asymmetries of small-world feature were observed only in normal controls but disappeared in patients with BD.
Regional Properties of the Hemispheric Networks
Hemispheric and Group Differences
The statistical results of nodal efficiency differences are summarized in Figure 2. Using Bonferroni-correction, we observed five regions exhibited significant hemispheric differences (P < 0.05) (Figure 2A). Among these five brain regions, regions with significant leftward advantage in nodal efficiency mainly included the anterior cingulate gyrus (ACG) and the inferior parietal lobule (IPL), whereas regions with significant rightward nodal efficiency were predominantly located at the supramarginal gyrus (SMG), the angular gyrus (ANG), and the rolandic operculum (ROL). Moreover, regions with significant group differences (BD < CN, P < 0.05, Bonferroni-corrected) in nodal efficiency included the amygdala (AMYG), IPL, putamen (PUT), and temporal pole (middle) (TPOmid) (Figure 2B). The PUT region could not be marked in Figure 2B. Furthermore, we found significant group-by-hemisphere interactions in the ROL (F1,106 = 12.930, P < 0.001), the lingual gyrus (LING) (F1,106 = 19.011, P < 0.001), the superior occipital gyrus (SOG) (F1,106 = 21.221, P < 0.001), the IPL (F1,106 = 36.225, P < 0.001), the SMG (F1,106 = 21.617, P < 0.001), the inferior temporal gyrus (ITG) (F1,106 = 15.900, P < 0.001), and the ACG (F1,106 = 10.547, P = 0.002, Figure 2C). Post hoc analysis indicated that these interaction effects resulted from significantly reduced rightward hemispheric asymmetry in nodal efficiency in patients with BD (Figure 3). Specially, we revealed significant group-by-hemisphere interaction in the supplementary motor area (SMA) (F1,106 = 34.907, P < 0.001). Post hoc analysis showed that this interaction in SMA was attributed to significantly rightward hemispheric asymmetry in nodal efficiency in patients with BD (Figure 3).
FIGURE 2. Cortical regions with significant differences (A) hemispheric difference, (B) group difference, and (C) group-by-hemispheric interaction on the nodal efficiency. The color bar denotes F-values. Significant difference was identified with a threshold value of P < 0.05 (Bonferroni-corrected). Significantly different regions were overlaid on surface maps provided by BrainNet Viewer software (Xia et al., 2013).
FIGURE 3. The nodal efficiency of normal controls and bipolar disorder in several regions including PUT, ROL, and SFGdor. Error bar indicates standard error. ∗∗∗P < 0.001, ∗∗P < 0.01.
The asymmetry score of nodal efficiency indicated that the region was leftward or rightward in each group subject. Consistent with the prior GLM hemispheric results, both groups exhibited significant hemispheric asymmetry (P < 0.05, Bonferroni-corrected) in nodal efficiency. For the normal controls shown in Figure 4A, the nodal efficiency with rightward asymmetry covered the inferior occipital gyrus (IOG), the orbitofrontal gyrus, the middle part (ORBmid), Cuneus (CUN), and ANG, SMG, SOG, and ROL. In contrast, the leftward asymmetric nodes were mainly located at the superior frontal gyrus, medial part (SFGmed), SMA, and Heschl (HES). For patients with BD shown in Figure 4B, the regions with significant leftward asymmetries in nodal efficiency involved the ITG, LING, SFGmed, ACG, middle occipital gyrus (MOG), HES, and IPL regions. Regions with significant rightward asymmetries in nodal efficiency were located at the ANG, Pallidium (PAL), CAL, and SMA. Additionally, the group differences in asymmetry score can find the reasons resulted from group-by-hemisphere interaction. When comparing the group differences in the asymmetry score, significant differences (P < 0.05, Bonferroni-corrected) were observed in the SMG, SOG, ROL, ITG, LING, IPL, ACG, and SMA. Specifically, patient with BD showed significant less rightward asymmetry in the ROL, SOG, and SMG, and more leftward asymmetry in the IPL, ACG, ITG, and LING, which were attributed to significantly more reductions nodal efficiency in right hemispheric. Only nodal efficiency of SMA in left hemisphere significantly reduced in patients with BD when compared with normal controls, and resulting in reversed hemispheric asymmetry.
FIGURE 4. The asymmetry scores of nodal efficiency for the hemispheric network with 45 nodes. Each bar indicates the mean asymmetry scores (AS). Error bars represent standard error. For both groups of (A) normal controls and (B) patients with BD, regions with significant AS are marked with asterisk (∗∗∗P < 0.001, Bonferroni-corrected). All the regions in each group were ascending sorted according to the statistical t-values. Similarly, the color bar indicated the statistical t-values. The spatial distributions of cortical regions with significant AS in both groups were also presented with BrainNet Viewer software (Xia et al., 2013). Cortical regions with significant group difference were presented with yellow background.
Relationship Between Hemispheric Asymmetry and BD Symptom Severity
Interestingly, we found the asymmetry scores of regional properties were correlated with BD symptom severity no matter using the Pearson or Spearman method. As shown in Figure 5, the nodal efficiency in ROL showed a significant and positive correlation with YMRS [AS(ROL), Pearson: r = 0.24, P = 0.05; Spearman: r = 0.29, P = 0.02], and the nodal efficiency in IPL showed a significant and negative correlation with YMRS [AS(IPL), Pearson: r = -0.31, P = 0.02; Spearman: r = -0.25, P = 0.04].
FIGURE 5. Partial correlation coefficients between the nodal asymmetry scores and clinical features of patients with BD.
This current study employed DTI tractography and graph theory to examine the hemispheric asymmetries in brain WM networks in patients with BD. The hemispheric networks of patients were observed to have abnormal small-world property and reduced in global integration. Significant group-by-hemisphere interaction was revealed in the global efficiency, characteristic path length, and small-world property, which was attributed to significantly reduced global integration and increased small-word characteristic of right hemisphere in patients with BD compared with those in normal controls. Specifically, we found atypical asymmetric nodal efficiency in several regions mostly located at the parietal, temporal, and occipital areas. Furthermore, we revealed that the network properties were significantly correlated with symptom severity in BD. Our findings might provide a potential neural biomarkers of that the altered hemispheric asymmetries in brain anatomical networks for clinical presentation of BD.
Hemisphere-Related Differences in Small-World Properties
In present work, hemispheric anatomical networks were used to evaluate the abnormal topological properties in patients with BD. Consistent with previous studies on whole-brain anatomical networks (Puetz et al., 2016; O’Donoghue et al., 2017), the hemispheric networks of both normal controls and BD patients preserved significant small-world property, that was significantly more clustered than random networks and had approximately the same characteristic path length as random networks (Watts and Strogatz, 1998). Moreover, we observed significant group differences in topological properties between normal controls and patients with BD. The smaller global efficiency, longer path length, and larger small-world property were shown in patients with BD compared with normal controls, suggesting a deficit in global integration (Munarini, 2013). Consistently, previous studies on whole-brain networks (Collin et al., 2016; Spielberg et al., 2016; Roberts et al., 2018) also showed smaller global efficiency and longer characteristic path length in patients with BD. Overall, these results suggested that the deficits in global integration were common in brain networks in patients with BD. Importantly, these abnormal global integrations were hemisphere-independent, which might be mainly due to the deficits of intra-hemispheric connections. Contemporary theories suggest that the complex clinical presentation of BD can be described as an aberration in the efficiency of information exchange between separate neural networks in the human brain (Vargas et al., 2013). Along this notion, our findings suggested the hemisphere-independent anatomical network with significantly smaller global efficiency and longer characteristic path length provides evidence for the hemispheric anatomical networks in patients with BD as a disconnection syndrome (Dae-Jin et al., 2013; Collin et al., 2016), especially for the right hemisphere.
The currently observed rightward asymmetry in global efficiency for normal controls suggests that the right hemisphere is intra-connected in a better integrated way, allowing for more efficient communication at the hemispheric level. Consistently, rightward asymmetry in network efficiency has been reported in healthy adults(Iturria-Medina et al., 2008; Sun et al., 2015). Interestingly, such rightward advantage in global integration in normal controls was absent in patients with BD. Instead, a roughly symmetrical pattern of global integration at the global level was found. The absence of rightward asymmetry in global integration was mainly due to the decrease of global integration in right hemisphere. According to lateralization of brain cognation theories, such absence of rightward asymmetry in global integration might underlie the BD dysfunctions in attention, visuospatial abilities (Cullen et al., 2016), and emotion regulation (Schwartz et al., 1975), which are considered to be dominantly processed in the right hemisphere. The decreased global integration in the right hemisphere was supported by the WM destruction in BD, specifically in the cingulum, corpus callosum, and superior longitudinal fasciculus (Ho et al., 2017). Our findings extend earlier work and provide network evidence that patients with BD have abnormal asymmetries in hemispheric networks and deficits in global integration in the functional networks of the right hemisphere.
Hemisphere-Related Differences in Nodal Properties
Compared with normal controls, patients with BD exhibited reduced nodal efficiency in several regions (including the AMYG, the IPL, the PUT, and the TPOmid). Previous DTI studies (Rosso et al., 2007; Passarotti et al., 2012) also reported that patients with BD exhibited degenerated WM connectivity in these regions. For example, the disturbed anatomical connectivity (decreased FA) in the AMYG (Mcintosh et al., 2008) and TPOmid (Barnea-Goraly et al., 2009) may represent the reduction of nodal efficiency observed in this study. Haller et al. (2011) found that the gray matter of the PUT in patients with BD was significantly smaller than that in normal controls. The aberrations in nodal efficiency provided regional evidence to underline the neurobiological basis of BD.
In addition to the group differences in nodal efficiency, significant group-by-hemisphere interactions were found in ROL, SMG, ITG, SOG, and LING, attributing to significantly reduced hemispheric asymmetry in patients with BD compared with normal controls. The nodal efficiency of patients with BD showed more decreased in right hemisphere than left hemisphere relative to normal controls. Consistently, previous studies on BD revealed disrupted connectivity in the right hemisphere including ROL (Gernsbacher and Kaschak, 2003) and SMG (Wang et al., 2016), ITG, SOG, and LING (Green, 2006; Bearden et al., 2015). The right ROL and SMG regions were demonstrated involving emotional regulation (Silani et al., 2013). Considering one of significant symptoms of patients with BD was emotional regulation (Pavuluri et al., 2006, 2007), the reduction of nodal efficiency in these regions were proposed to associate with the impaired in emotional regulation in BD (Pavuluri et al., 2007). Moreover, it has been demonstrated that the right SOG, ITG, and LING regions were as associated with the visuospatial processing (Green, 2006; Bearden et al., 2015). The reduction of nodal efficiency in these occipital and temporal regions might be associated with the deficits in visuospatial functions (Green, 2006; Bearden et al., 2015).
Moreover, we also found that attention cognitive function regions including the ACG and IPL exhibited significant group-by-hemisphere interaction and group differences in asymmetry score. Because the nodal efficiency of right hemisphere in BD patients showed more reduction than those in left hemisphere, the ACG and IPL exhibited significant leftward hemispheric asymmetry. Consistently, previous studies have reported decreased connectivity in these two regions in patients with BD. For example, significantly decreased connectivity was found in the right ACG (Anand et al., 2009; Wang et al., 2009). Barnea-Goraly et al. (2009) reported that adolescents with BD had lower FA values than normal controls in the IPL region. Moreover, studies have showed that alterations in two regions are closely to the brain function in patients with BD. Patients with BD was demonstrated with decrease activation in the right ACG relative to normal controls during response inhibition task (Passarotti et al., 2010). The gray matter volume in the right parietal cortical regions correlated positively with the better inhibitory control in BD patients (Haldane et al., 2008). We further proposed that the abnormal asymmetries of nodal efficiency in patients with BD were associated with abnormal attention or inhibition function.
Interestingly, the patients with BD exhibited the significant rightward hemispheric asymmetry only in the SMA region, whereas a reversed asymmetry (leftward) in normal controls. The SMA region was revealed significant difference between the BD group and normal controls (Khadka et al., 2013). In addition, convergent studies (Caligiuri et al., 2003, 2004) have reported that patients with BD exhibited greater activities in right SMA during motor performance than those in left SMA. These results suggested that abnormalities in asymmetries of nodal efficiency in SMA influenced processing of motor function in BD (Puri et al., 2010; Rashid et al., 2014). These findings might provide potential neural biomarkers of the altered asymmetries in nodal efficiency for clinical presentation of BD.
Relationship Between Regional Asymmetries and BD Symptom Severity
An interesting finding of the current study was that within patients with BD, the topological properties in hemispheric networks were related to BD symptom severity. We found that nodal efficiency was associated with the YMRS, which reflected the severity of the current manic episode; a higher YMRS resulted in a more severe episode. A positive relationship between the asymmetry score of nodal efficiency in the ROL and the YMRS was revealed, indicating that the more connectivity in the right ROL, the more severe symptoms of BD. Recently, Gao et al. (2017) found that connectivity in the right ROL was positively associated with the classification of BD. Moreover, a negative relationship between the asymmetry score of nodal efficiency in the IPL and YMRS was observed. We suggested that the negative correlation might result from significantly decreased connectivity in the right IPL. That is, the larger YMRS score the less connected communication in the right IPL region. Similarly, Barnea-Goraly et al. (2009) reported that adolescents with BD had lower FA values than normal subjects in inferior parietal region. Our observation suggested that abnormal properties of hemispheric asymmetries may underlie the dysfunctions existed in patients with BD.
Using the DTI deterministic tractography method and graph theory, the current work evaluated the hemispheric effects on brain anatomical networks in patients with BD. The results revealed abnormalities in hemispheric asymmetries in patients with BD compared with those in normal controls. For the global network properties, the hemispheric asymmetry in global efficiency was significantly decreased and small-world property was significantly increased. Compared with normal controls, the nodal efficiency of patients also showed decreased rightward asymmetry mainly in the frontal, occipital, parietal, and temporal lobes. Exceptionally, the SMA region in patients with BD showed increased rightward asymmetries, attributing to a significantly reduction of the efficiency in left SMA. The asymmetry score of nodal efficiencies in the IPL and ROL exhibited correlations with clinical features of BD. These observations highlight that the altered hemispheric asymmetries in brain anatomical networks and the potential of brain hemispheric network measures as neural biomarkers for clinical presentation of BD. Our findings suggested that abnormal asymmetries in properties of hemispheric networks may underlie the dysfunctions in emotion and attention in patients with BD.
The research of this paper still has some limitations. In this study, the subjects was chosen from one site with small number of BD patient, limiting the statistical results of hemispheric effects and may bring about the type I error on brain network and region alterations. In order to evaluate accurately alteration in hemispheric asymmetry in patients with BD, the subjects with large number or from multi-site will be considered in the future work. Additionally, studies have reported that both gender and age are potential factors linked to brain asymmetry. Hence, we will further examine gender and age effects on hemispheric asymmetries in brain anatomical networks in BD.
After receiving a thorough explanation, all participants gave written informed consent according to the procedures approved by the Taiyuan University of Technology.
BW completed the entire study of the experiment and writing. TL, MZ, SZ, YN, and XW revised the manuscript. TY, RC, and JX provided advice and guidance. DL provided the research ideas.
This work was supported by the National Natural Science Foundation of China (Grant Nos. 61503272, 61503273, 61873178, and 61741212), the project funded by the China Postdoctoral Science Foundation (Grant No. 2016 M601287), the Shanxi Provincial Foundation for Returned Scholars, China (Grant No. 2016–037), the Scientific Research Foundation for Returned Overseas Chinese Scholars (Grant No. 20151098-2), the Key Research and Development Project in Shanxi (Grant No. 201603D111014).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer YC and handling Editor declared their shared affiliation, at the time of the review.
Anand, A., Li, Y., Wang, Y., Lowe, M. J., and Dzemidzic, M. (2009). Resting state corticolimbic connectivity abnormalities in unmedicated bipolar disorder and unipolar depression. Psychiatry Res. 171, 189–198. doi: 10.1016/j.pscychresns.2008.03.012
Bai, F., Shu, N., Yuan, Y., Shi, Y., Yu, H., Wu, D., et al. (2012). Topologically convergent and divergent structural connectivity patterns between patients with remitted geriatric depression and amnestic mild cognitive impairment. J. Neurosci. 32, 4307–4318. doi: 10.1523/JNEUROSCI.5061-11.2012
Barnea-Goraly, N., Chang, K. D., Karchemskiy, A., Howe, M. E., and Reiss, A. L. (2009). Limbic and corpus callosum aberrations in adolescents with bipolar disorder: a tract-based spatial statistics analysis. Biol. Psychiatry 66, 238–244. doi: 10.1016/j.biopsych.2009.02.025
Bearden, C. E., Hoffman, K. M., and Cannon, T. D. (2015). The neuropsychology and neuroanatomy of bipolar affective disorder: a critical review. Bipolar Disord. 3, 106–150. doi: 10.1034/j.1399-5618.2001.030302.x
Beyer, J. L., Taylor, W. D., Macfall, J. R., Kuchibhatla, M., Payne, M. E., Provenzale, J. M., et al. (2005). Cortical white matter microstructural abnormalities in bipolar disorder. Neuropsychopharmacology 30, 2225–2229. doi: 10.1038/sj.npp.1300802
Bruno, S., Cercignani, M., and Ron, M. A. (2008). White matter abnormalities in bipolar disorder: a voxel-based diffusion tensor imaging study. Bipolar Disord. 10, 460–468. doi: 10.1111/j.1399-5618.2007.00552.x
Caligiuri, M. P., Brown, G. G., Meloy, M., Eyler, L. T., Kindermann, S. S., Eberson, S., et al. (2004). A functional magnetic resonance imaging study of cortical asymmetry in bipolar disorder. Bipolar Disord. 6, 183–196. doi: 10.1111/j.1399-5618.2004.00116.x
Caligiuri, M. P., Brown, G. G., Meloy, M. J., Eberson, S. C., Kindermann, S. S., Frank, L. R., et al. (2003). An fMRI study of affective state and medication on cortical and subcortical brain regions during motor performance in bipolar disorder. Psychiatry Res. 123, 171–182. doi: 10.1016/S0925-4927(03)00075-1
Collin, G., Scholtens, L. H., Kahn, R. S., Mhj, H., and van den Heuvel, M. P. (2017). Affected anatomical rich club and structural-functional coupling in young offspring of schizophrenia and bipolar disorder patients. Biol. Psychiatry 82, 746–755. doi: 10.1016/j.biopsych.2017.06.013
Collin, G., van den Heuvel, M. P., Abramovic, L., Vreeker, A., de Reus, M. A., van Haren, N. E., et al. (2016). Brain network analysis reveals affected connectome structure in bipolar I disorder. Hum. Brain Mapp. 37, 122–134. doi: 10.1002/hbm.23017
Cullen, B., Ward, J., Graham, N. A., Deary, I. J., Pell, J. P., Smith, D. J., et al. (2016). Prevalence and correlates of cognitive impairment in euthymic adults with bipolar disorder: a systematic review. J. Affect. Disord. 205, 165–181. doi: 10.1016/j.jad.2016.06.063
Dae-Jin, K., Bolbecker, A. R., Josselyn, H., Olga, R., Olaf, S., Hetrick, W. P., et al. (2013). Disturbed resting state EEG synchronization in bipolar disorder: a graph-theoretic analysis. Neuroimage Clin. 2, 414–423. doi: 10.1016/j.nicl.2013.03.007
Dolcos, F., Rice, H. J., and Cabeza, R. (2002). Hemispheric asymmetry and aging: right hemisphere decline or asymmetry reduction. Neurosci. Biobehav. Rev. 26, 819–825. doi: 10.1016/S0149-7634(02)00068-4
Forde, N. J., O’Donoghue, S., Scanlon, C., Emsell, L., Chaddock, C., Leemans, A., et al. (2015). Structural brain network analysis in families multiply affected with bipolar I disorder. Psychiatry Res. 234, 44–51. doi: 10.1016/j.pscychresns.2015.08.004
Gao, S., Osuch, E. A., Wammes, M., Theberge, J., Jiang, T. Z., Calhoun, V. D., et al. (2017). Discriminating bipolar disorder from major depression based on kernel SVM using functional independent components. Paper Presented at the IEEE International Workshop on Machine Learning for Signal Processing, Aalborg. doi: 10.1109/MLSP.2017.8168110
Gong, G., He, Y., Concha, L., Lebel, C., Gross, D. W., Evans, A. C., et al. (2009). Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography. Cereb. Cortex 19, 524–536. doi: 10.1093/cercor/bhn102
Good, C. D., Johnsrude, I., Ashburner, J., Henson, R. N., Friston, K. J., and Frackowiak, R. S. (2001). Cerebral asymmetry and the effects of sex and handedness on brain structure: a voxel-based morphometric analysis of 465 normal adult human brains. Neuroimage 14, 685–700. doi: 10.1006/nimg.2001.0857
Haldane, M., Cunningham, G., Androutsos, C., and Frangou, S. (2008). Structural brain correlates of response inhibition in bipolar disorder I. J. Psychopharmacol. 22, 138–143. doi: 10.1177/0269881107082955
Haller, S., Xekardaki, A., Delaloye, C., Canuto, A., Lövblad, K. O., Gold, G., et al. (2011). Combined analysis of grey matter voxel-based morphometry and white matter tract-based spatial statistics in late-life bipolar disorder. J. Psychiatry Neurosci. 36, 391–401. doi: 10.1503/jpn.100140
Highley, J. R., Mcdonald, B., Walker, M. A., Esiri, M. M., and Crow, T. J. (1999). Schizophrenia and temporal lobe asymmetry. A post-mortem stereological study of tissue volume. Br. J. Psychiatry 175, 127–134. doi: 10.1192/bjp.175.2.127
Ho, N. F., Li, Z., Ji, F., Wang, M., Kuswanto, C. N., Sum, M. Y., et al. (2017). Hemispheric lateralization abnormalities of the white matter microstructure in patients with schizophrenia and bipolar disorder. J. Psychiatry Neurosci. 42, 242–251. doi: 10.1503/jpn.160090
Iturriamedina, Y., Pérez, F. A., Morris, D. M., Canalesrodríguez, E. J., Haroon, H. A., García, P. L., et al. (2011). Brain hemispheric structural efficiency and interconnectivity rightward asymmetry in human and nonhuman primates. Cereb. Cortex 21, 56–67. doi: 10.1093/cercor/bhq058
Iturria-Medina, Y., Sotero, R. C., Canales-Rodríguez, E. J., Alemán-Gómez, Y., and Melie-García, L. (2008). Studying the human brain anatomical network via diffusion-weighted MRI and graph theory. Neuroimage 40, 1064–1076. doi: 10.1016/j.neuroimage.2007.10.060
Jenkinson, M., Bannister, P., Brady, M., and Smith, S. (2002). Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825–841. doi: 10.1006/nimg.2002.1132
Kafantaris, V., Kingsley, P., Ardekani, B., Saito, E., Lencz, T., Lim, K., et al. (2009). Lower orbital frontal white matter integrity in adolescents with bipolar I disorder. J. Am. Acad. Child Adolesc. Psychiatry 48, 79–86. doi: 10.1097/CHI.0b013e3181900421
Khadka, S., Meda, S. A., Stevens, M. C., Glahn, D. C., Calhoun, V. D., Sweeney, J. A., et al. (2013). Is aberrant functional connectivity a psychosis endophenotype? A resting state functional magnetic resonance imaging study. Biol. Psychiatry 74, 458–466. doi: 10.1016/j.biopsych.2013.04.024
Kwon, J. S., Youn, T., and Jung, H. Y. (1996). Right hemisphere abnormalities in major depression: quantitative electroencephalographic findings before and after treatment. J. Affect. Disord. 40, 169–173. doi: 10.1016/0165-0327(96)00057-2
Leow, A., Ajilore, O., Zhan, L., Arienzo, D., Gadelkarim, J., Zhang, A., et al. (2013). Impaired inter-hemispheric integration in bipolar disorder revealed using brain network analyses. Biol. Psychiatry 73, 183–193. doi: 10.1016/j.biopsych.2012.09.014
Mahon, K., Burdick, K. E., Ikuta, T., Braga, R. J., Gruner, P., Malhotra, A. K., et al. (2013). Abnormal temporal lobe white matter as a biomarker for genetic risk of bipolar disorder. Biol. Psychiatry 73, 177–182. doi: 10.1016/j.biopsych.2012.07.033
Mcintosh, A. M., Muñoz, M. S., Lymer, G. K., Mckirdy, J., Hall, J., Sussmann, J. E., et al. (2008). White matter tractography in bipolar disorder and schizophrenia. Biol. Psychiatry 64, 1088–1092. doi: 10.1016/j.biopsych.2008.07.026
Miller, S. R., Miller, C. J., Bloom, J. S., Hynd, G. W., and Craggs, J. G. (2006). Right hemisphere brain morphology, attention-deficit hyperactivity disorder (ADHD) subtype, and social comprehension. J. Child Neurol. 21, 139–144. doi: 10.1177/08830738060210021901
Mori, S., Crain, B. J., Chacko, V. P., and van Zijl, P. C. (1999). Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging. Ann. Neurol. 45, 265–269. doi: 10.1002/1531-8249(199902)45:2<265::AID-ANA21>3.0.CO;2-3
O’Donoghue, S., Kilmartin, L., O’Hora, D., Emsell, L., Langan, C., Mcinerney, S., et al. (2017). Anatomical integration and rich-club connectivity in euthymic bipolar disorder. Psychol. Med. 47, 1609–1623. doi: 10.1017/S0033291717000058
Park, H., Westin, C., Kubicki, M., Maier, S. E., Niznikiewicz, M., Baer, A., et al. (2004). White matter hemisphere asymmetries in healthy subjects and in schizophrenia: a diffusion tensor MRI study. Neuroimage 23, 213–223. doi: 10.1016/j.neuroimage.2004.04.036
Passarotti, A. M., Ellis, J., Wegbreit, E., Stevens, M. C., and Pavuluri, M. N. (2012). Reduced functional connectivity of prefrontal regions and amygdala within affect and working memory networks in pediatric bipolar disorder. Brain Connect. 2, 320–334. doi: 10.1089/brain.2012.0089
Passarotti, A. M., Sweeney, J. A., and Pavuluri, M. N. (2010). Neural correlates of response inhibition in pediatric bipolar disorder and attention deficit hyperactivity disorder. Psychiatry Res. 181, 36–43. doi: 10.1016/j.pscychresns.2009.07.002
Pavuluri, M. N., O’Connor, M. M., Harral, E., and Sweeney, J. A. (2007). Affective neural circuitry during facial emotion processing in pediatric bipolar disorder. Biol. Psychiatry 62, 158–167. doi: 10.1016/j.biopsych.2006.07.011
Pavuluri, M. N., Schenkel, L. S., Aryal, S., Harral, E. M., Hill, S. K., Herbener, E. S., et al. (2006). Neurocognitive function in unmedicated manic and medicated euthymic pediatric bipolar patients. Am. J. Psychiatry 163, 286–293. doi: 10.1176/appi.ajp.163.2.286
Puetz, V. B., Parker, D., Kohn, N., Dahmen, B., Verma, R., and Konrad, K. (2016). Altered brain network integrity after childhood maltreatment: a structural connectomic DTI-study. Hum. Brain Mapp. 38, 855–868. doi: 10.1002/hbm.23423
Puri, B. K., Agour, M., Gunatilake, K. D., Fernando, K. A., Gurusinghe, A. I., and Treasaden, I. H. (2010). Reduction in left supplementary motor area grey matter in adult female fibromyalgia sufferers with marked fatigue and without affective disorder: a pilot controlled 3-T magnetic resonance imaging voxel-based morphometry study. J. Int. Med. Res. 38, 1468–1472. doi: 10.1177/147323001003800429
Rashid, B., Damaraju, E., Pearlson, G. D., and Calhoun, V. D. (2014). Dynamic connectivity states estimated from resting fMRI Identify differences among Schizophrenia, bipolar disorder, and healthy control subjects. Front. Hum. Neurosci. 8:897. doi: 10.3389/fnhum.2014.00897
Ratnarajah, N., Rifkin-Graboi, A., Fortier, M. V., Chong, Y. S., Kwek, K., Saw, S. M., et al. (2013). Structural connectivity asymmetry in the neonatal brain. Neuroimage 75, 187–194. doi: 10.1016/j.neuroimage.2013.02.052
Rimol, L. M., Hartberg, C. B., Nesvåg, R., Fennemanotestine, C., Hagler, D. J., Pung, C. J., et al. (2010). Cortical thickness and subcortical volumes in schizophrenia and bipolar disorder. Biol. Psychiatry 117, 41–50. doi: 10.1016/j.biopsych.2010.03.036
Roberts, G., Perry, A., Lord, A., Frankland, A., Leung, V., Holmespreston, E., et al. (2018). Structural dysconnectivity of key cognitive and emotional hubs in young people at high genetic risk for bipolar disorder. Mol. Psychiatry 23, 413–421. doi: 10.1038/mp.2016.216
Rosso, I. M., Killgore, W. D., Cintron, C. M., Gruber, S. A., Tohen, M., and Yurgelun-Todd, D. A. (2007). Reduced amygdala volumes in first-episode bipolar disorder and correlation with cerebral white matter. Biol. Psychiatry 61, 743–749. doi: 10.1016/j.biopsych.2006.07.035
Schulte, T., Mülleroehring, E. M., Sullivan, E. V., and Pfefferbaum, A. (2012). White matter fiber compromise contributes differentially to attention and emotion processing impairment in alcoholism, HIV-infection, and their comorbidity. Neuropsychologia 50, 2812–2822. doi: 10.1016/j.neuropsychologia.2012.07.042
Shu, N., Liu, Y., Li, K., Duan, Y., Wang, J., Yu, C., et al. (2011). Diffusion tensor tractography reveals disrupted topological efficiency in white matter structural networks in multiple sclerosis. Cereb. Cortex 21, 2565–2577. doi: 10.1093/cercor/bhr039
Sigalovsky, I. S., Fischl, B., and Melcher, J. R. (2006). Mapping an intrinsic MR property of gray matter in auditory cortex of living humans: a possible marker for primary cortex and hemispheric differences. Neuroimage 32, 1524–1537. doi: 10.1016/j.neuroimage.2006.05.023
Silani, G., Lamm, C., Ruff, C. C., and Singer, T. (2013). Right supramarginal gyrus is crucial to overcome emotional egocentricity bias in social judgments. J. Neurosci. 33, 15466–15476. doi: 10.1523/JNEUROSCI.1488-13.2013
Spielberg, J. M., Beall, E. B., Hulvershorn, L. A., Altinay, M., Karne, H., and Anand, A. (2016). Resting state brain network disturbances related to hypomania & depression in medication-free bipolar disorder. Neuropsychopharmacology 41, 3016–3024. doi: 10.1038/npp.2016.112
Sun, Y., Chen, Y., Collinson, S. L., Bezerianos, A., and Sim, K. (2015). Reduced hemispheric asymmetry of brain anatomical networks is linked to schizophrenia: a connectome study. Cereb. Cortex 27, 602–615. doi: 10.1093/cercor/bhv255
Tzouriomazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., 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. doi: 10.1006/nimg.2001.0978
Vargas, C., López-Jaramillo, C., and Vieta, E. (2013). A systematic literature review of resting state network—functional MRI in bipolar disorder. J. Affect. Disord. 150, 727–735. doi: 10.1016/j.jad.2013.05.083
Wang, F., Kalmar, J. H., He, Y., Jackowski, M., Chepenik, L. G., Edmiston, E. E., et al. (2009). Functional and structural connectivity between the perigenual anterior cingulate and amygdala in bipolar disorder. Biol. Psychiatry 66, 516–521. doi: 10.1016/j.biopsych.2009.03.023
Wang, Y., Zhong, S., Jia, Y., Sun, Y., Wang, B., Liu, T., et al. (2016). Disrupted resting-state functional connectivity in nonmedicated bipolar disorder. Radiology 280, 529–536. doi: 10.1148/radiol.2016151641
Watkins, K. E., Paus, T., Lerch, J. P., Zijdenbos, A., Collins, D. L., Neelin, P., et al. (2001). Structural asymmetries in the human brain: a voxel-based statistical analysis of 142 MRI scans. Cereb. Cortex 11, 868–877. doi: 10.1093/cercor/11.9.868
Wessa, M., Houenou, J., Leboyer, M., Chanraud, S., Poupon, C., Martinot, J. L., et al. (2009). Microstructural white matter changes in euthymic bipolar patients: a whole-brain diffusion tensor imaging study. Bipolar Disord. 11, 504–514. doi: 10.1111/j.1399-5618.2009.00718.x
Yin, X., Han, Y., Ge, H., Xu, W., Huang, R., Zhang, D., et al. (2013). Inferior frontal white matter asymmetry correlates with executive control of attention. Hum. Brain Mapp. 34, 796–813. doi: 10.1002/hbm.21477
Keywords: bipolar disorder, diffusion tensor imaging, graph theory, hemispheric asymmetry, structural connectivity
Citation: Wang B, Li T, Zhou M, Zhao S, Niu Y, Wang X, Yan T, Cao R, Xiang J and Li D (2018) The Abnormality of Topological Asymmetry in Hemispheric Brain Anatomical Networks in Bipolar Disorder. Front. Neurosci. 12:618. doi: 10.3389/fnins.2018.00618
Received: 03 July 2018; Accepted: 14 August 2018;
Published: 03 September 2018.
Edited by:Jing Sui, Institute of Automation (CAS), China
Reviewed by:Yue Cui, Chinese Academy of Sciences, China
Xiang-Zhen Kong, Max-Planck-Institut für Psycholinguistik, Netherlands
Copyright © 2018 Wang, Li, Zhou, Zhao, Niu, Wang, Yan, Cao, Xiang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dandan Li, email@example.com