Jing Wei1
Shubao Wei1Rongxing Yang1Lu Yang1Qiong Yin1Huihui Li1Yuhong Qin2Yiwu Lei2Chao Qin1Jingqun Tang1
Shuguang Luo1*
Wenbin Guo3*- 1Department of Neurology, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 2Department of Radiology, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 3Department of Psychiatry, The Second Xiangya Hospital of Central South University, Changsha, China
Objective: Several networks in human brain are involved in the development of blepharospasm. However, the underlying mechanisms for this disease are poorly understood. A voxel-mirrored homotopic connectivity (VMHC) method was used to quantify the changes in functional connectivity between two hemispheres of the brain in patients with blepharospasm.
Methods: Twenty-four patients with blepharospasm and 24 healthy controls matched by age, sex, and education were recruited. The VMHC method was employed to analyze the fMRI data. The support vector machine (SVM) method was utilized to examine whether these abnormalities could be applied to distinguish the patients from the controls.
Results: Compared with healthy controls, patients with blepharospasm showed significantly high VMHC in the inferior temporal gyrus, interior frontal gyrus, posterior cingulate cortex, and postcentral gyrus. No significant correlation was found between abnormal VMHC values and clinical variables. SVM analysis showed a combination of increased VMHC values in two brain areas with high sensitivities and specificities (83.33 and 91.67% in the combined inferior frontal gyrus and posterior cingulate cortex; and 83.33 and 87.50% in the combined inferior temporal gyrus and postcentral gyrus).
Conclusion: Enhanced homotopic coordination in the brain regions associated with sensory integration networks and default-mode network may be underlying the pathophysiology of blepharospasm. This phenomenon may serve as potential image markers to distinguish patients with blepharospasm from healthy controls.
Introduction
Blepharospasm (BSP) is a common clinical type of focal dystonia characterized by involuntary blinking and eyelid spasms (Hallett et al., 2008). BSP is a chronic disease with an incidence rate of 4.2/10 million in general population (Steeves et al., 2012). This disease occurs mostly in adults and manifested with an increased frequency of blinking in the early stage and consistent closed-eyes or even functional blindness in the late stages. Furthermore, BSP will lead to poor quality of life, decline in employment, anxiety, and depression (Muller et al., 2002; Biuk et al., 2013; Valls-Sole and Defazio, 2016; Bedarf et al., 2017). However, the pathophysiology of BSP remains unclear.
Recent techniques, such as neuroimaging, facilitate the exploration of structural and functional abnormalities in patients with BSP. For example, Baker et al. reported regional abnormalities in the cortical and subcortical brain areas in BSP (Baker et al., 2003; Yang et al., 2013). A voxel-based morphometry study demonstrated increased gray matter in the putamen (Etgen et al., 2006) and primary sensorimotor cortex (Suzuki et al., 2011) in patients with BSP. Zhou et al. (Zhou et al., 2013) used the amplitude of a low-frequency fluctuation method and found abnormalities in the bilateral somatosensory regions in patients with BSP.
The above mentioned studies indicate that multiple abnormal nodes are involved in the pathophysiology of BSP. However, the reported findings are inconsistent in regard to specific networks. For example, a report with independent component analysis exhibited decreased connectivity within the sensorimotor network and the right fronto-parietal network but increased connectivity in the salience network in patients with BSP (Huang et al., 2017). By contrast, a functional connectivity study reported that patients with BSP showed altered functional connectivity at rest in widespread brain regions including basal ganglia, cerebellar, primary/secondary sensorimotor, and visual areas (Jochim et al., 2018). Therefore, BSP may be related to network abnormality, and dysfunction of these nodes is likely to be involved in the pathogenesis of BSP. However, whether or not these abnormalities result in simultaneous bilateral eyelid spasms in BSP remain ambiguous. Moreover, complex functional connectivities exist between two cerebral hemispheres, and their effects are complementary and coordinated. Both eyelids simultaneously experience spasms in BSP, and this event is related to abnormalities of functional interaction between two cerebral hemispheres. Thus, BSP might have also resulted from brain connection dysfunction. However, whether or not an abnormal functional interaction exists between two cerebral hemispheres in BSP should be confirmed.
Resting-state fMRI (rs-fMRI) has attracted increasing attention (Biswal et al., 1995), because it can detect patterns of coherent intrinsic activities of the brain and interactions between two hemispheres. Voxel-mirrored homotopic connectivity (VMHC) (Zuo et al., 2010) is a method of rs-fMRI used to analyze functional homotopy between two hemispheres. The synchrony in patterns of spontaneous activity between homotopic regions in each hemisphere is an important feature of the functional structure of the brain (Stark et al., 2008). The VMHC is designed to directly compare the interhemispheric resting-state functional connectivity. This process can also measure the correlations between blood oxygen level-dependent time series and reflect the communication pattern of information between two cerebral hemispheres. Thus, VMHC is extremely important for information integration of the brain. Recent studies used VMHC to show abnormal homotopic connection in patients with schizophrenia (Hoptman et al., 2012; Guo et al., 2014b, 2017a,b) and their unaffected siblings (Guo et al., 2014a), depression (Guo et al., 2013a,b), somatization disorder (Su et al., 2016), and Parkinson’s disease (Hu et al., 2015). Moreover, Anderson et al. reported that homotopic resting-state functional connectivity was disrupted in individuals with autism (Anderson et al., 2011), and this finding indicated that homotopic function is an important part of the brain function. In addition, excision of corpus callosum could damage the sensory, motor, and cognitive functions of human brain. This phenomenon implies that the coordination between two hemispheres plays an important role for human behavior (Luo et al., 2015). The abovementioned studies suggest that abnormal interhemispheric connectivity may be an important factor in the occurrence of BSP. Based on the onset pattern of simultaneous effect on the bilateral eyelids of patients with BSP and the presence of different networks and loop injuries, we hypothesized that hemisphere connectivity might play a key role in the pathogenesis of BSP. Using VMHC to compare the fMRI data of patients with BSP with those of healthy controls, we evaluated spatial heterogeneity of interhemispheric functional connectivity and clarified the pathogenesis of BSP. To date, findings in brain networks of patients with BSP are limited. To our knowledge, the current study is the first work to employ VMHC to investigate the resting-state functional connectivity of the two hemispheres in patients with BSP.
Patients with BSP usually demonstrate psychiatric symptoms, such as anxiety and depression (Muller et al., 2002; Yang et al., 2017). Self-rating Anxiety Scale (SAS) and Self-rating Depression Scale (SDS) were used to evaluate the symptoms of anxiety and depression, respectively, to control the confounding impact from these conditions. Prior to the rating, all patients were screened to determine whether they were suffering from anxiety or depression. Moreover, Jankovic Rating Scale (JRS-S) was used to assess the severity of eyelid spasm. Then, The VMHC values of different brain areas were correlated with the degree of JRS-S, the course of disease, and the SAS and SDS scores to examine the relationship between VMHC values and clinical symptoms. Finally, SVM was employed to examine whether abnormal VMHC values could discriminate patients from controls with high accuracy and specificity.
Materials and Methods
Subjects
We recruited 26 outpatients with BSP who visited the Department of Neurology in the First Affiliated Hospital of Guangxi Medical University between November 2012 and June 2014. The inclusion criteria for the BSP group were as follows: (1) met the criteria of BSP diagnosis according to the Clinical Guidelines of BSP (Defazio et al., 2013); (2) absence of structural changes with conventional MRI examination; (3) had not used botulinum toxin within 3 months prior to the study; (4) had not used medication for dystonia within 1 month prior to the study; and (5) right-handedness. Patients with the following conditions were excluded: (1) secondary BSP from other diseases, such as hepatolenticular degeneration and dry eye; and (2) history of neurological and psychiatric disorders.
The following information was collected from patients with BSP: sex, age, education level, course of disease, degree of illness severity, scores of SAS, and scores of SDS. JRS-S was used to assess the severity of BSP. SAS and SDS were used to examine the severity of anxiety and depression.
Twenty-four healthy volunteers from the community were recruited as age- sex-, and education level-matched controls. All individual controls were right-handed and had no history of severe neuropsychiatric diseases, medical illness, or family history of neurological or psychiatric disorders from the first-degree relatives.
All participants read and signed a consent form before the examination. The Ethics Committee of the First Affiliated Hospital, Guangxi Medical University approved this study.
Image Acquisition
Data were acquired using a German Siemens Trio Tim 3.0T scanner (Erlangen, Germany). Head fixers and earplugs were used for all subjects to reduce head movement and machine noise, respectively. The subjects were required to remain motionless, awake, and eye-closed during image acquisition. All patients underwent routine examination (T1W1 and T2W1) to exclude intracranial lesions. The following parameters were used for functional imaging: repetition time/echo time (TR/TE) = 2000/30 ms, 30 slices, 64 × 64 matrix, 90° flip angle, 24 cm FOV, 4 mm section thickness, 0.4 mm gap, and 250 volume (500 s).
Data Preprocessing
Images were preprocessed using the data assistant software (DPABI) (Yan et al., 2016). The first 10 time points were removed. All participants had less than 2 mm maximum displacement in the x, y, or z-axis and 2° of angular motion during data acquisition. The images were then normalized to the standard SPM81 echo planar imaging template with resampling to 3 × 3 × 3 mm3 voxels. The processed images were smoothed with an isotropic Gaussian kernel (full-width at half-maximum = 8 mm). Finally, the acquired data were subjected to temporal bandpass filtering (0.01–0.08 Hz) and linear detrending to reduce the effect of low-frequency drifts and physiological high-frequency noise. Spurious covariates and their temporal derivatives, including Friston-24 head motion parameters, white matter signals, and cerebrospinal fluid signals, were removed from the data using linear regression.
Interhemispheric Correlation and Statistical Analysis
The software REST2 was used to analyze VMHC, and the details of VMHC computation has been expounded in a previous study (Zuo et al., 2010). Individual VMHC maps were generated for each participant by computing Pearson correlation (Fisher z-transformed) between a given voxel and a corresponding voxel in the contralateral hemisphere. Correlation values were then Fisher z-transformed to improve the normality. The resultant values were applied for group comparisons and generation of the VMHC maps.
Individual-level VMHC maps were analyzed using group-level voxel wise t-test to determine regional group differences in VMHC. The significance level was set to p < 0.05 for multiple comparisons corrected by the Gaussian Random Field theory (voxel significance: p < 0.001, cluster significance: p < 0.05). Given that resting-state functional connectivity could be influenced by micro motions (Power et al., 2012), the frame wise displacement values were computed as a covariate for each subject in the group comparisons. Pearson correlation was used to evaluate the relationships between VMHC values with significant group differences and degree of JRS-S, SAS scores, SDS scores, and the course of disease after the normality of these variables being assessed.
The demographic and clinical information was compared with the two-sample t-test and χ2-test. The significance level was set to P < 0.05.
Classification Analysis Using SVM
SVM was conducted to evaluate the possibility of abnormal VMHC in these clusters to discriminate patients from healthy controls using the LIBSVM software package3 in MATLAB (Chang and Lin, 2011). The LIBSVM classifier was trained by providing examples of the form <x, c>, where x represents the VMHC values of these abnormal clusters, and c is the class label (c = +1 for patients with BSP and c = −1 for healthy controls). A sample set was divided into a training set and a test set for SVM to evaluate the classification performance on the unobserved data. The leave-one-pair-out method was applied for the LIBSVM classifier algorithm. We constructed a random SVM cluster based on the brain fMRI data of the subjects for classification and feature selection. To acquire the optimal sensitivity and specificity, default functional kernels of Gaussian radial basis and the grid search method were applied to optimize the parameters with the “leave-one-subject-out” method.
Results
Demographics and Clinical Characteristics of the Participants
A total of 24 patients were included in the further analysis (two patients were excluded due to excessive head movement). About 20.83 and 29.17% of patients with BSP had anxiety and depressive symptoms, respectively. Additionally, 3 (12.5%) patients with BSP had comorbidity of depression. Moreover, 19 patients with BSP (79.16%) had sensory tricks. These patients might temporarily improve eyelid spasms by wearing glasses, having material in the mouth, or touching the cheeks, forehead, or jaw. Half of the 24 patients (50.00%) experienced worsened eyelid spasms when talking. The patient and control groups did not differ significantly in sex, age, and education levels (Table 1).
TABLE 1. Clinical profile of the participants.
VMHC: Group Differences
Compared with the healthy control group, significantly high VMHC was found in the inferior temporal gyrus, interior frontal gyrus, posterior cingulate cortex and postcentral gyrus in the patient group (Figure 1 and Table 2).
FIGURE 1. Statistical maps showing VMHC differences between patients and controls. Red denotes increased VMHC, and color bars indicate T-values from t-tests between groups. VMHC, voxel-mirrored homotopic connectivity.
TABLE 2. Significant group differences in VMHC.
Correlation Analysis
No significant correlation was found between abnormal VMHC values and the severity of symptom, clinical course, and scores of SAS and SDS in the patients (P > 0.05).
SVM for Classification Analysis
SVM analysis was performed to determine whether or not abnormal VMHC values could satisfactorily discriminate patients with BSP from healthy controls. The results showed that the VMHC values in a single brain region could not discriminate patients with BSP from healthy controls with optimal sensitivity, specificity, and accuracy (Table 3), but the patients may be distinguished with high sensitivity, specificity, and accuracy (more than 80%) using a combination of two brain regions. Among the brain regions, the ability of discriminating patients with BSP from the healthy controls by the combination of the VMHC values in the inferior frontal gyrus and posterior cingulate cortex was optimal with an accuracy of 87.5% (42 of 48 in the 2 groups), a sensitivity of 83.33% (20 of 24 in the BSP group), and a specificity of 91.67% (22 of 24 in the control group). The combination of the VMHC values in the inferior temporal gyrus and postcentral gyrus showed an accuracy of 85.4% (41 of 48 in the 2 groups), a sensitivity of 83.33% (20 of 24 in the BSP group), and a specificity of 87.5% (21 of 24 in the control group; Figures 2, 3). The other combinations such as inferior frontal gyrus/inferior temporal gyrus or posterior cingulate cortex/postcentral gyrus had unsatisfactory sensitivity and specificity.
TABLE 3. Discriminate the patients from the controls by using the VMHC values of a single region with the SVM method.
FIGURE 2. Visualization of classification by support vector machine (SVM) using the combination of the VMHC values in the Inferior Frontal Gyrus and Posterior Cingulate Cortex. Left: SVM parameters selection result (3D visualization) [Grid Search Method]: Best c = 0.70711; Best g = 0.088388; Right: the visualization of classification with a combination of the VMHC values in the Inferior Frontal Gyrus and Posterior Cingulate Cortex.
FIGURE 3. Visualization of classification by support vector machine (SVM) using the combination of the VMHC values in the abnormal brain regions. Left: SVM parameters selection result (3D visualization) [Grid Search Method]: Best c = 0.70711; Best g = 11.3137; Right: the visualization of classification with a combination of the VMHC values in the Inferior Temporal Gyrus and Postcentral Gyrus.
Discussion
In this study, higher VMHC values were found in some brain regions, including the inferior temporal gyrus, interior frontal gyrus, posterior cingulate cortex, and postcentral gyrus in patients with BSP than those in healthy controls (P < 0.05). Moreover, SVM analysis revealed that the VMHC values in a single brain region could not discriminate patients with BSP from healthy controls, contrary to the combination of two brain regions with high sensitivity, specificity, and accuracy. Moreover, no significant correlation was found between abnormal VMHC values and the severity of symptoms, clinical course, and scores of SAS and SDS in the patient group (P > 0.05).
Previous studies indicated that BSP was related to abnormal function of the basal ganglia (Perlmutter et al., 1997; Federico et al., 1998; Esmaeli-Gutstein et al., 1999). Meanwhile, other studies (Martino et al., 2011; Zhou et al., 2013) have reported the abnormalities in the central posterior and posterior cingulate gyrus in patients with BSP. Our report is consistent with those from previous studies (Martino et al., 2011; Zhou et al., 2013). Patients with BSP demonstrated various patterns of sensory tricks, such as motor, imaginary, forcible, and reverse sensory tricks (Hallett, 2002; Ramos et al., 2014), and more likely to suffer from significant prepulse inhibition hand stimulation when presented sensory tricks (Gomez-Wong et al., 1998). One report indicated increased sensory afferents in the eyelid activity control in patients with BSP (Hallett, 2002). Another report implied that the changes in sensory input played an important role in sensory integration abnormalities which might result from increased sensory afferents (Patel et al., 2014). Thus, balancing the movement output and sensory afferents might be a valuable strategy to improve symptoms of eyelid spasm. Furthermore, this report also implied that abnormalities in the sensory center and sensory integration were involved in the pathogenesis of dystonia. Wong et al. reported that 8 of 17 patients with BSP (47.06%) showed sensory tricks (Gomez-Wong et al., 1998). In our patient group, 19 of the 24 BSP patients (79.16%) presented sensory tricks. The eyelid spasms could be alleviated when the patient wears glasses and touches his cheek, forehead, jaw, or other parts of his body. The VMHC of the postcentral gyrus was significantly increased. This phenomenon suggested that patients with BSP had to process more complex information and programs to control eyelid movement. Thus, the integration of sensory-motor information might play an important role in the pathogenesis of BSP. Furthermore, 12 of the 24 BSP patients (50.00%) demonstrated reverse sensory tricks. These patients experienced worse symptoms when they talked (Greene and Bressman, 1998). Interestingly, a stimulus or action can alleviate the symptoms in some patients but not effective in other patients with BSP or even worsening their symptoms (Wider et al., 2004; Martino et al., 2010; Ramos et al., 2014). Thus, a sensory trick has a very unique individual nature. However, the complex input mechanism of sensory integration is unclear, and more studies are needed to explore the related neurophysiological factors in patients with BSP.
The patients in our study presented increased connection from bilateral inferior frontal gyrus, where Brodmann areas 44 and 45 of the inferior frontal gyrus are the main parts of Broca area. Previous study suggested that the frontal lobe plays an integrated role in sensory information processing, and this concept supports the hypothesis that the Broca area is part of the sensory network (Okada et al., 2016). Recent studies showed that the human hyperdirect low-frequency interactions between the prefrontal cortex and subthalamic nucleus (STN) support the regulation of several related brain functions (Kelley et al., 2018), and the STN is the key point of the motor network (Alexander and Crutcher, 1990). In Parkinson’s disease, high frequencies of typical STN deep brain stimulation protocols are used to treat motor symptoms. Therefore, we can speculate that the STN can also become a therapeutic target in the future therapy of BSP.
In 2001, Raichle et al. (Raichle et al., 2001) first proposed default mode network (DMN). This network exhibits high levels of activity at rest. However, DMN becomes deactivated when specific goal-directed behavior is required. DMN consists of the posterior cingulate cortex, precuneus, medial prefrontal cortex, ventral anterior cingulate cortex, inferior parietal lobule, and several temporal lobes (Fransson, 2006). Zhou et al. found abnormalities in the medial prefrontal cortex and posterior cingulate cortex, and these abnormalities suggested significant differences in the DMN in patients with BSP (Zhou et al., 2013). The posterior cingulate cortex, which involves monitoring sensation, stereotactic positioning, and other functions, is the only region in the cingulate cortex that accepts thalamic occipital medial projection (Baleydier and Mauguiere, 1980). A bidirectional connection exists between the pulvinar thalamus and the secondary sensory cortex area (Brodmann areas 37, 39, and 40). The joint fibers connect to the ipsilateral occipital lobe, temporal lobe, frontal lobe, and contralateral brain region. An abnormal homozygous connection was found in the posterior cingulate cortex. This result suggested an abnormal loop between the posterior cingulate cortex and the thalamus and secondary sensory cortex. The posterior cingulate cortex is the intermediate connector between the two hemispheres. Moreover, the bilateral cerebral hemispheres are asymmetric. These findings suggest that the abovementioned changes may be compensatory reactions to the cortical center. Furthermore, the posterior cingulate cortex is an important node as a static state of DMN and might affect the functions and related connections of DMN. Damage in the bilateral temporal lobe of the monkey brain resulted in the degeneration of the fiber that connects the posterior cingulate cortex and the temporal lobe (Papez, 1995). In addition, the temporal gyrus is relevant to the information integration of vision. A study that experimented on color discrimination and visual contrast perception has demonstrated that the visual impairment of patients with BSP was not dependent on the severity of the disease or the course of the disease (Buttner et al., 1999). In the present study, increased VMHC values were found in the posterior cingulate cortex and inferior temporal gyrus, and these values suggested that patients with BSP might have abnormal visual-spatial sensory integration. However, whether or not the posterior cingulate cortex and inferior temporal gyrus are the exact key points of the visual spatial sensory integration remains unclear. The question needs to be confirmed by animal experiments. Our results provide evidence for the establishment of animal models in the future.
No correlation was found between SAS or SDS scores and increased VMHC values in the patients. Three of the 24 patients with BSP were diagnosed with depression in this study. This finding indicated the higher morbidity of depression in patients with BSP than in the normal population (The morbidity of depression in the normal population was approximately 2–6%) (Blazer et al., 1994; Dunlop et al., 2003; Ohayon and Schatzberg, 2003; Gu et al., 2013). Moreover, several scholars inferred that patients with focal dystonia easily suffer from anxiety and depressive symptoms (Voon et al., 2010). SAS and SDS are self-rating scales which might be affected by some confounders, such as educational levels, intelligence, illness duration, and social environment. No correlation was found between abnormal VMHC and anxiety and depression of BSP, which is consistent with a previous study (Stamelou et al., 2012). Furthermore, no correlation was found between abnormal VMHC and the course of disease and severity degrees of BSP symptoms. In previous studies, no difference was found in the age and severity of dystonia between patients with BSP with or without mood disorders (Fabbrini et al., 2010). This phenomenon suggested that mood disorders were not a direct response to focal dystonia (Lencer et al., 2009; Horovitz et al., 2012). Thus, patients with BSP may suffer from anxiety and depression, but these conditions do not interfere with the evolution of the disease. Previous studies reported that changes in the gray matter volume in patients with BSP were not related to illness duration (Martino et al., 2011; Horovitz et al., 2012) and severity degree (Martino et al., 2011) of BSP symptoms. However, another study has shown that abnormal gray matter density in patients with BSP was related to the course of disease (Suzuki et al., 2011). The conflicting results may be due to abnormal sensory motor plasticity (Abbruzzese and Berardelli, 2011) and neuronal remodeling uncertainty (Opavsky et al., 2006).
SVM analyses showed that the VMHC values in one single brain region could not discriminate patients with BSP from healthy controls with optimal sensitivity, specificity, and accuracy, contrary to a combination of two brain regions with high sensitivity, specificity, and accuracy. Meanwhile, the specificity was particularly remarkable, because every healthy control was correctly classified. Therefore, the combination of high VMHC values in the inferior frontal gyrus and posterior cingulate cortex, as well as in the inferior temporal gyrus and postcentral gyrus, may serve as a potential image marker to distinguish the patients with BSP from healthy controls. Our report is consistent with those of fMRI studies in patients with schizophrenia (Li et al., 2018; Wang et al., 2018). Therefore, early detection and improved accuracy of the diagnosis in patients with BSP may be achieved based on the fMRI imaging using the combination of inferior frontal gyrus and posterior cingulate cortex or the inferior temporal gyrus and postcentral gyrus.
This study has several limitations. Patients with BSP demonstrate a wide spectrum of symptoms with different severity degrees and illness durations. However, stratified analysis was not performed because of the small sample size. More patients should be recruited to confirm the present report. We will use stratified analysis to study the relationship between abnormal VMHC values and the symptoms or course of disease. Then, whether or not connection abnormality is a fundamental change or compensatory performance of BSP will be confirmed. The VMHC method reflects the synchrony in patterns of spontaneous activity between homotopic regions in each hemisphere and does not fully reflect the functional state of the whole brain. Therefore, we can simultaneously employ the regional homogeneity method to analyze local consistency changes in the whole brain regions.
Conclusion
Abnormalities of BSP were found in the inferior temporal gyrus, interior frontal gyrus, posterior cingulate cortex, and post-central gyrus. The results suggested that those related brain areas may be good candidate regions to explore the nature of BSP and highlight the significance of sensory integration and DMN in the pathophysiology of the disorder. Moreover, a combination of high VMHC values in two brain areas can serve as a potential image marker to distinguish patients with BSP from healthy controls.
Ethics Statement
This study was carried out in accordance with the recommendations of the First Affiliated Hospital of Guangxi Medical University. The protocol was approved by the the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University. All subjects gave written informed consent in accordance with the Declaration of Helsinki.
Author Contributions
SL and WG designed the study. JW, RY, SW, LY, HL, QY, YL, YQ, JT, and CQ collected the original imaging data. WG, JW, and SW managed and analyzed the imaging data. JW and SW wrote the first draft of the manuscript.
Funding
This study was supported by grants from the National Key R&D Program of China (2016YFC1307100), the National Natural Science Foundation of China (Grant Nos. 81571310 and 81771447), and the Guangxi Appropriate Technology for Medical and Health Research and Development Project (Grant No. S201415-05).
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.
Acknowledgments
We thank all individuals who served as the study participants.
Footnotes
- ^ http://www.fil.ion.ucl.ac.uk/spm
- ^ http://resting-fmri.sourceforge.net
- ^ http://www.csie.ntu.edu.tw/cjlin/libsvm
References
Abbruzzese, G., and Berardelli, A. (2011). Further progress in understanding the pathophysiology of primary dystonia. Mov. Disord. 26, 1185–1186. doi: 10.1002/mds.23707
Alexander, G. E., and Crutcher, M. D. (1990). Functional architecture of basal ganglia circuits: neural substrates of parallel processing. Trends Neurosci. 13, 266–271. doi: 10.1016/0166-2236(90)90107-L
Anderson, J. S., Druzgal, T. J., Froehlich, A., Dubray, M. B., Lange, N., Alexander, A. L., et al. (2011). Decreased interhemispheric functional connectivity in autism. Cereb. Cortex 21, 1134–1146. doi: 10.1093/cercor/bhq190
Baker, R. S., Andersen, A. H., Morecraft, R. J., and Smith, C. D. (2003). A functional magnetic resonance imaging study in patients with benign essential blepharospasm. J. Neuroophthalmol. 23, 11–15. doi: 10.1097/00041327-200303000-00003
Baleydier, C., and Mauguiere, F. (1980). The duality of the cingulate gyrus in monkey. Neuroanatomical study and functional hypothesis. Brain 103, 525–554. doi: 10.1093/brain/103.3.525
Bedarf, J. R., Kebir, S., Michelis, J. P., Wabbels, B., and Paus, S. (2017). Depression in blepharospasm: a question of facial feedback? Neuropsychiatr. Disease Treat. 13, 1861–1865. doi: 10.2147/NDT.S141066
Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409
Biuk, D., Karin, A. A., Matic, S., Barac, J., Benasic, T., and Stiglmayer, N. (2013). Quality of life in patients with blepharospasm. Coll. Antropol. 37, 29–33.
Blazer, D. G., Kessler, R. C., Mcgonagle, K. A., and Swartz, M. S. (1994). The prevalence and distribution of major depression in a national community sample: the National Comorbidity Survey. Am. J. Psychiatry 151, 979–986. doi: 10.1176/ajp.151.7.979
Buttner, T., Kuhn, W., Dietz, M., Muller, T., Postert, T., Przuntek, H., et al. (1999). Impaired visual function in focal idiopathic dystonia. Eur. Neurol. 41, 94–98. doi: 10.1159/000008010
Chang, C.-C., and Lin, C.-J. (2011). LIBSVM: a library for support vector machines. ACM Trans. Intell. Syst. Technol. 2, 1–27. doi: 10.1145/1961189.1961199
Defazio, G., Hallett, M., Jinnah, H. A., and Berardelli, A. (2013). Development and validation of a clinical guideline for diagnosing blepharospasm. Neurology 81, 236–240. doi: 10.1212/WNL.0b013e31829bfdf6
Dunlop, D. D., Song, J., Lyons, J. S., Manheim, L. M., and Chang, R. W. (2003). Racial/ethnic differences in rates of depression among preretirement adults. Am. J. Public Health 93, 1945–1952. doi: 10.2105/AJPH.93.11.1945
Esmaeli-Gutstein, B., Nahmias, C., Thompson, M., Kazdan, M., and Harvey, J. (1999). Positron emission tomography in patients with benign essential blepharospasm. Ophthalmic Plast. Reconstr. Surg. 15, 23–27. doi: 10.1097/00002341-199901000-00006
Etgen, T., Muhlau, M., Gaser, C., and Sander, D. (2006). Bilateral grey-matter increase in the putamen in primary blepharospasm. J. Neurol. Neurosurg. Psychiatry 77, 1017–1020. doi: 10.1136/jnnp.2005.087148
Fabbrini, G., Berardelli, I., Moretti, G., Pasquini, M., Bloise, M., Colosimo, C., et al. (2010). Psychiatric disorders in adult-onset focal dystonia: a case-control study. Mov. Disord. 25, 459–465. doi: 10.1002/mds.22983
Federico, F., Simone, I. L., Lucivero, V., Defazio, G., De Salvia, R., Mezzapesa, D. M., et al. (1998). Proton magnetic resonance spectroscopy in primary blepharospasm. Neurology 51, 892–895. doi: 10.1212/WNL.51.3.892
Fransson, P. (2006). How default is the default mode of brain function? Further evidence from intrinsic BOLD signal fluctuations. Neuropsychologia 44, 2836–2845. doi: 10.1016/j.neuropsychologia.2006.06.017
Gomez-Wong, E., Marti, M. J., Tolosa, E., and Valls-Sole, J. (1998). Sensory modulation of the blink reflex in patients with blepharospasm. Arch. Neurol. 55, 1233–1237. doi: 10.1001/archneur.55.9.1233
Greene, P. E., and Bressman, S. (1998). Exteroceptive and interoceptive stimuli in dystonia. Mov. Disord. 13, 549–551. doi: 10.1002/mds.870130329
Gu, L., Xie, J., Long, J., Chen, Q., Chen, Q., Pan, R., et al. (2013). Epidemiology of major depressive disorder in mainland china: a systematic review. PLoS One 8:e65356. doi: 10.1371/journal.pone.0065356
Guo, W., Liu, F., Chen, J., Wu, R., Li, L., Zhang, Z., et al. (2017a). Treatment effects of olanzapine on homotopic connectivity in drug-free schizophrenia at rest. World J. Biol. Psychiatry 1–9. doi: 10.1080/15622975.2017.1346280
Guo, W., Liu, F., Chen, J., Wu, R., Li, L., Zhang, Z., et al. (2017b). Family-based case-control study of homotopic connectivity in first-episode, drug-naive schizophrenia at rest. Sci Rep. 7:43312. doi: 10.1038/srep43312
Guo, W., Jiang, J., Xiao, C., Zhang, Z., Zhang, J., Yu, L., et al. (2014a). Decreased resting-state interhemispheric functional connectivity in unaffected siblings of schizophrenia patients. Schizophr. Res. 152, 170–175. doi: 10.1016/j.schres.2013.11.030
Guo, W., Xiao, C., Liu, G., Wooderson, S. C., Zhang, Z., Zhang, J., et al. (2014b). Decreased resting-state interhemispheric coordination in first-episode, drug-naive paranoid schizophrenia. Prog. Neuropsychopharmacol. Biol. Psychiatry 48, 14–19. doi: 10.1016/j.pnpbp.2013.09.012
Guo, W., Liu, F., Dai, Y., Jiang, M., Zhang, J., Yu, L., et al. (2013a). Decreased interhemispheric resting-state functional connectivity in first-episode, drug-naive major depressive disorder. Prog. Neuropsychopharmacol. Biol. Psychiatry 41, 24–29. doi: 10.1016/j.pnpbp.2012.11.003
Guo, W., Liu, F., Xue, Z., Gao, K., Liu, Z., Xiao, C., et al. (2013b). Decreased interhemispheric coordination in treatment-resistant depression: a resting-state fMRI study. PLoS One 8:e71368. doi: 10.1371/journal.pone.0071368
Hallett, M. (2002). Blepharospasm: recent advances. Neurology 59, 1306–1312. doi: 10.1212/01.WNL.0000027361.73814.0E
Hallett, M., Evinger, C., Jankovic, J., and Stacy, M. (2008). Update on blepharospasm: report from the BEBRF International Workshop. Neurology 71, 1275–1282. doi: 10.1212/01.wnl.0000327601.46315.85
Hoptman, M. J., Zuo, X. N., D’angelo, D., Mauro, C. J., Butler, P. D., Milham, M. P., et al. (2012). Decreased interhemispheric coordination in schizophrenia: a resting state fMRI study. Schizophr. Res. 141, 1–7. doi: 10.1016/j.schres.2012.07.027
Horovitz, S. G., Ford, A., Najee-Ullah, M. A., Ostuni, J. L., and Hallett, M. (2012). Anatomical correlates of blepharospasm. Transl. Neurodegener. 1:12. doi: 10.1186/2047-9158-1-12
Hu, X., Zhang, J., Jiang, X., Zhou, C., Wei, L., Yin, X., et al. (2015). Decreased interhemispheric functional connectivity in subtypes of Parkinson’s disease. J. Neurol. 262, 760–767. doi: 10.1007/s00415-014-7627-x
Huang, X.-F., Zhu, M.-R., Shan, P., Pei, C.-H., Liang, Z.-H., Zhou, H.-L., et al. (2017). Multiple neural networks malfunction in primary blepharospasm: an independent components analysis. Front. Hum. Neurosci. 11:235. doi: 10.3389/fnhum.2017.00235
Jochim, A., Li, Y., Gora-Stahlberg, G., Mantel, T., Berndt, M., Castrop, F., et al. (2018). Altered functional connectivity in blepharospasm/orofacial dystonia. Brain Behav. 8:e00894. doi: 10.1002/brb3.894
Kelley, R., Flouty, O., Emmons, E. B., Kim, Y., Kingyon, J., Wessel, J. R., et al. (2018). A human prefrontal-subthalamic circuit for cognitive control. Brain 141, 205–216. doi: 10.1093/brain/awx300
Lencer, R., Steinlechner, S., Stahlberg, J., Rehling, H., Orth, M., Baeumer, T., et al. (2009). Primary focal dystonia: evidence for distinct neuropsychiatric and personality profiles. J. Neurol. Neurosurg. Psychiatry 80, 1176–1179. doi: 10.1136/jnnp.2008.170191
Li, R. R., Lyu, H. L., Liu, F., Lian, N., Wu, R. R., Zhao, J. P., et al. (2018). Altered functional connectivity strength and its correlations with cognitive function in subjects with ultra-high risk for psychosis at rest. CNS Neurosci. Ther. doi: 10.1111/cns.12865 [Epub ahead of print].
Luo, C., Guo, X., Song, W., Zhao, B., Cao, B., Yang, J., et al. (2015). Decreased resting-state interhemispheric functional connectivity in parkinson’s disease. Biomed Res. Int. 2015:692684. doi: 10.1155/2015/692684
Martino, D., Di Giorgio, A., D’ambrosio, E., Popolizio, T., Macerollo, A., Livrea, P., et al. (2011). Cortical gray matter changes in primary blepharospasm: a voxel-based morphometry study. Mov. Disord. 26, 1907–1912. doi: 10.1002/mds.23724
Martino, D., Liuzzi, D., Macerollo, A., Aniello, M. S., Livrea, P., and Defazio, G. (2010). The phenomenology of the geste antagoniste in primary blepharospasm and cervical dystonia. Mov. Disord. 25, 407–412. doi: 10.1002/mds.23011
Muller, J., Kemmler, G., Wissel, J., Schneider, A., Voller, B., Grossmann, J., et al. (2002). The impact of blepharospasm and cervical dystonia on health-related quality of life and depression. J. Neurol. 249, 842–846. doi: 10.1007/s00415-002-0733-1
Ohayon, M. M., and Schatzberg, A. F. (2003). Using chronic pain to predict depressive morbidity in the general population. Arch. Gen. Psychiatry 60, 39–47. doi: 10.1001/archpsyc.60.1.39
Okada, K., Rogalsky, C., O’grady, L., Hanaumi, L., Bellugi, U., Corina, D., et al. (2016). An fMRI study of perception and action in deaf signers. Neuropsychologia 82, 179–188. doi: 10.1016/j.neuropsychologia.2016.01.015
Opavsky, R., Hlustik, P., and Kanovsky, P. (2006). Cortical plasticity and its implications for focal hand dystonia. Biomed. Pap. Med. Fac. Univ. Palacky Olomouc Czech. Repub. 150, 223–226. doi: 10.5507/bp.2006.031
Papez, J. W. (1995). A proposed mechanism of emotion. 1937. J. Neuropsychiatry Clin. Neurosci. 7, 103–112. doi: 10.1176/jnp.7.1.103
Patel, N., Jankovic, J., and Hallett, M. (2014). Sensory aspects of movement disorders. Lancet Neurol. 13, 100–112. doi: 10.1016/S1474-4422(13)70213-8
Perlmutter, J. S., Stambuk, M. K., Markham, J., Black, K. J., Mcgee-Minnich, L., Jankovic, J., et al. (1997). Decreased [18F]spiperone binding in putamen in idiopathic focal dystonia. J. Neurosci. 17, 843–850. doi: 10.1523/JNEUROSCI.17-02-00843.1997
Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018
Raichle, M. E., Macleod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682. doi: 10.1073/pnas.98.2.676
Ramos, V. F., Karp, B. I., and Hallett, M. (2014). Tricks in dystonia: ordering the complexity. J. Neurol. Neurosurg. Psychiatry 85, 987–993. doi: 10.1136/jnnp-2013-306971
Stamelou, M., Edwards, M. J., Hallett, M., and Bhatia, K. P. (2012). The non-motor syndrome of primary dystonia: clinical and pathophysiological implications. Brain 135, 1668–1681. doi: 10.1093/brain/awr224
Stark, D. E., Margulies, D. S., Shehzad, Z. E., Reiss, P., Kelly, A. M., Uddin, L. Q., et al. (2008). Regional variation in interhemispheric coordination of intrinsic hemodynamic fluctuations. J. Neurosci. 28, 13754–13764. doi: 10.1523/JNEUROSCI.4544-08.2008
Steeves, T. D., Day, L., Dykeman, J., Jette, N., and Pringsheim, T. (2012). The prevalence of primary dystonia: a systematic review and meta-analysis. Mov. Disord. 27, 1789–1796. doi: 10.1002/mds.25244
Su, Q., Yao, D., Jiang, M., Liu, F., Long, L., Dai, Y., et al. (2016). Decreased interhemispheric functional connectivity in insula and angular gyrus/supramarginal gyrus: significant findings in first-episode, drug-naive somatization disorder. Psychiatry Res. 248, 48–54. doi: 10.1016/j.pscychresns.2016.01.008
Suzuki, Y., Kiyosawa, M., Wakakura, M., Mochizuki, M., and Ishii, K. (2011). Gray matter density increase in the primary sensorimotor cortex in long-term essential blepharospasm. Neuroimage 56, 1–7. doi: 10.1016/j.neuroimage.2011.01.081
Valls-Sole, J., and Defazio, G. (2016). Blepharospasm: update on epidemiology, clinical aspects, and pathophysiology. Front. Neurol. 7:45. doi: 10.3389/fneur.2016.00045
Voon, V., Butler, T. R., Ekanayake, V., Gallea, C., Ameli, R., Murphy, D. L., et al. (2010). Psychiatric symptoms associated with focal hand dystonia. Mov. Disord. 25, 2249–2252. doi: 10.1002/mds.23250
Wang, S., Zhan, Y., Zhang, Y., Lyu, L., Lyu, H., Wang, G., et al. (2018). Abnormal long- and short-range functional connectivity in adolescent-onset schizophrenia patients: a resting-state fMRI study. Prog. Neuropsychopharmacol. Biol. Psychiatry 81, 445–451. doi: 10.1016/j.pnpbp.2017.08.012
Wider, C., Ghika, J., Bogousslavsky, J., and Vingerhoets, F. (2004). Segmental dystonia induced by wearing glasses with a ribbon: an unusual case of a reverse sensory geste. Mov. Disord. 19, 966–967.
Yan, C. G., Wang, X. D., Zuo, X. N., and Zang, Y. F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4
Yang, J., Luo, C., Song, W., Chen, Q., Chen, K., Chen, X., et al. (2013). Altered regional spontaneous neuronal activity in blepharospasm: a resting state fMRI study. J. Neurol. 260, 2754–2760. doi: 10.1007/s00415-013-7042-8
Yang, J., Shao, N., Song, W., Wei, Q., Ou, R., Wu, Y., et al. (2017). Nonmotor symptoms in primary adult-onset cervical dystonia and blepharospasm. Brain Behav. 7:e00592. doi: 10.1002/brb3.592
Zhou, B., Wang, J., Huang, Y., Yang, Y., Gong, Q., and Zhou, D. (2013). A resting state functional magnetic resonance imaging study of patients with benign essential blepharospasm. J. Neuroophthalmol. 33, 235–240. doi: 10.1097/WNO.0b013e31828f69e5
Keywords: blepharospasm, resting state, functional magnetic resonance imaging, voxel-mirrored homotopic connectivity, support vector machine
Citation: Wei J, Wei S, Yang R, Yang L, Yin Q, Li H, Qin Y, Lei Y, Qin C, Tang J, Luo S and Guo W (2018) Voxel-Mirrored Homotopic Connectivity of Resting-State Functional Magnetic Resonance Imaging in Blepharospasm. Front. Psychol. 9:1620. doi: 10.3389/fpsyg.2018.01620
Received: 02 June 2018; Accepted: 13 August 2018;
Published: 11 September 2018.
Edited by:
Jiongjiong Yang, Peking University, ChinaCopyright © 2018 Wei, Wei, Yang, Yang, Yin, Li, Qin, Lei, Qin, Tang, Luo and Guo. 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: Shuguang Luo, cm9iZXJ0NTgyNDNAc29odS5jb20= Wenbin Guo, Z3Vvd2VuYmluNzZAY3N1LmVkdS5jbg==