Qu Liu1†
Xiao-fei Hu- 1Department of Neurology, Southwest Hospital, Third Military Medical University (Army Medical University), Chongqing, China
- 2Department of Radiology, Guizhou Provincial People’s Hospital, Guiyang, China
- 3Department of Centre for Disease Prevention and Control, Chengdu Military Region, Chengdu, China
- 4Department of Radiology, Southwest Hospital, Third Military Medical University (Army Medical University), Chongqing, China
This study investigated the degree of brain functional impairment in persistent somatoform pain disorder (PSPD) by examining changes in the patterns of brain functional hubs. Resting-state functional magnetic resonance imaging was performed in 21 PSPD patients with headache as the main symptom and 17 sex- and age-matched healthy controls. Degree centrality (DC) analysis as well as the connectivity among these hubs by functional connectivity (FC) analysis and Granger causality analysis (GCA) were performed to characterize abnormal brain networks in PSPD (Gaussian random field corrected: P < 0.001, Z > 3.09). The relationships between DC and connectivity and clinical parameters were also examined. DC values in the bilateral inferior occipital gyrus (IOG), bilateral calcarine fissure (CAL), and left paracentral lobule (PCL) and FC values of right IOG–left CAL, right IOG–right CAL, right IOG–left IOG, left CAL–right CAL, left CAL–left IOG, left CAL–left PCL, right CAL–left PCL, and left IOG–left PCL were lower in PSPD patients as compared to controls. A negative causal effect from the left CAL to the left paracentral lobule and a positive effect from the right CAL to the right IOG were observed in PSPD patients. Abnormal DC, FC, and signed-path coefficients in PSPD patients were negatively correlated with self-rating anxiety and depression scale scores. These results indicate that altered functional hubs and connectivity patterns in the somatosensory cortex may reflect emotional disturbance in PSPD patients.
Introduction
Persistent somatoform pain disorder (PSPD) is characterized by medically unexplained somatic symptoms (Smith et al., 2000). PSPD patients experience long-term, persistent, severe, and distressing pain (typically lasting for over 6 months) that cannot be attributed to a specific cause (World Health Organization [WHO], 1993; American Psychiatric Association [APA], 2013). Various studies have described emotional and cognitive impairment in PSPD patients (Luo et al., 2014, 2016). For instance, alexithymia and anxiety are very common symptoms in adolescents with PSPD (Burba et al., 2006). Others have reported that PSPD patients exhibited deficits in emotion perception but normal ability to recognize facial expressions (Schonenberg et al., 2014). Additionally, lower scores on the quality of life sub-scale of the 36-Item Short Form Survey were shown to be correlated with pain, depression, and anxiety scores in PSPD patients (Luo et al., 2014). However, none of these studies investigated the mechanism(s) underlying the impairment of brain function in PSPD patients.
Functional magnetic resonance imaging (fMRI) provides an effective means for studying PDSD-related cortical alterations. A task-based fMRI study showed that the prefrontal–temporal–limbic circuit was activated in PSPD patients in response to pinprick pain stimulation under negative emotional conditions, suggesting a potential role for emotion in pain processing and the pathophysiology of PSPD (Luo et al., 2016). Meanwhile, a resting-state (rs-)fMRI study reported changes in regional homogeneity (ReHo) in brain regions related to pain in PSPD patients during the resting state, including the default mode network, somatosensory and prefrontal cortices, and posterior cerebellum (Huang et al., 2016). A functional connectivity analysis showed altered co-activation within aberrant resting-state networks such as sensorimotor, default-mode, and salience networks, suggesting that PSPD patients experience large-scale reorganization (Zhao et al., 2017). However, these studies have focused on local spontaneous brain activity within selected brain regions based on a priori assumptions or by independent components analysis, and have not fully characterized the functional connectome of PSPD patient brains.
Degree centrality (DC) is a measure of network organization based on graph theory that has been applied to the identification of candidate functional hubs in several rsfMRI studies (Cole et al., 2010; Fransson et al., 2011; Di et al., 2013). This has allowed characterization of the functional relationships of a given node (voxel) within the full-brain connectivity matrix as opposed to relationships with specific nodes or networks (Tomasi and Volkow, 2011; Zuo et al., 2012). Thus, DC can reveal the complexity of the functional connectome of PSPD patients.
After representing the hubs of the brain functional connectome using DC, the subsequent analysis of the complexity and patterns the interactions among these functional hubs should been investigated with undirected and directed functional connectivity (FC). Directed FC can be estimated with traditional FC (Biswal et al., 1995), and Granger causality model could investigate positive and negative signed-path coefficients in an directed manner (Seth et al., 2015). Combining traditional FC and Granger causality analysis (GCA) get a more complete characterization of connectivity patterns among brain functional hubs in PSPD.
In this study, we tested the hypothesis that aberrations in brain functional hubs and their connectivity contribute to PSPD. Candidate functional hubs were identified based on DC, and FC analysis and GCA were used to investigate their interconnectivity. We also investigated the relationships among changes in brain function, clinical data, and neuropsychological performance. Our results provide insight into the neurological basis for the impairment of brain function in PSPD.
Materials and Methods
Participants
This study was approved by the Ethics Committee of Southwest Hospital (Chongqing, China), and all study subjects provided written, informed consent to participate in the study.
Consecutive PSPD patients were recruited from inpatients at the hospital and healthy controls (HCs) were recruited from the community between December 2016 and December 2017. A total of 21 patients were diagnosed with PSPD by an experienced specialist according to ICD-10 criteria; inclusion criteria for this group were as follows: (1) right-hand dominance (which was tested with the Edinburgh Handedness Inventory); (2) age between 18 and 65 years; and (3) clinical pain (headache as the main symptom) persistently for at least 6 months. All patients’ treatment and all the medication like antidepressants and pain-relieve drugs were taken after MRI scan. Additionally, 17 HCs matched in terms of age, sex, and handedness were recruited.
Exclusion criteria for all subjects were as follows: (1) presence of pain symptoms due to severe somatic disease; (2) pre-existing neurological or psychiatric disorder (including a history of seizures, global cognitive impairment, aphasia, neglect, substantial sensory disturbance, severe “major depressive disorders” in Statistical Manual for Mental Disorders (DSM)-V or claustrophobia); (3) existence of uncontrolled disease such as congestive heart failure, hypertension, cerebrovascular disease, and thyropathy; (4) metal clips in the brain; or (5) pneumonia at the time of enrollment.
Clinical Assessment
The Visual Analog Scale (VAS), Zung Self-Rating Anxiety Scale (SAS), and Zung Self-Rating Depression Scale (SDS) were used to assess pain characteristics, anxiety symptoms, and depression symptoms in each patient, respectively.
The VAS is a psychometric response scale that is the most widely used quantitative approach for evaluating pain (Reips and Funke, 2008). Clinical pain was scored from 0 (no pain) to 10 (extreme pain).
The SAS is a 20-item self-assessment devised to measure anxiety levels based on scores for cognitive, autonomic, motor, and central nervous system symptoms (Zung, 1971, 1974), with total raw scores ranging from 20 to 80. The raw score is then converted to an anxiety index (normal, 20–44; mild to moderate, 45–59; marked to severe, 60–74; and extreme, 75–80) for clinical interpretation of anxiety level.
The SDS is a self-administered survey for assessing the level of depression symptom (Zung et al., 1965), with scores ranging from 20 to 80. Depression scores were categorized as follows: normal range, 20–44; mildly depressed, 45–59; moderately depressed, 60–69; and severely depressed, ≥70.
MRI Data Acquisition
The MRI scan was performed on the same day that clinical data were obtained. Images were acquired using a 3.0 T Siemens Tim Trio whole-body MRI system (Siemens Medical Solutions, Erlangen, Germany). Subjects were instructed to stay awake and close their eyes, and to try not to think of anything. Foam padding and earplugs were used to reduce head motion and scanner noise. Imaging data were collected transversely using an echo-planar imaging sequence with the following settings: TR = 2000 ms, TE = 30 ms, flip angle = 90°, FOV = 192 mm × 192 mm, in-plane matrix = 64 × 64, thickness = 3 mm, voxel size = 3.0 mm × 3.0 mm × 3.0 mm. For each subject, a total of 240 volumes were acquired with a scan time of 480 s. T1-weighted structural images were collected by volumetric three-dimensional magnetization prepared with the following rapid-acquisition gradient-echo sequence (TR = 1900 ms, TE = 2.52 ms, flip angle = 9°, slice thickness = 1 mm, slices number = 176, FOV = 256 mm × 256 mm, matrix size = 256 × 256 and voxel size = 1 mm × 1 mm × 1 mm and sagittal scanning).
MRI Data Analysis
All preprocessing steps were performed using the Data Processing Assistant for Resting-State fMRI (DPARSF2.31), which is based on the Statistical Parametric Mapping (SPM8) program2 (Chao-Gan and Yu-Feng, 2010). Intracranial tissue segmentation was performed with Voxel Based Morphometry Toolbox 8 v.435 software. The main steps included individual 3D T1-weighted structural images co-registered to the mean of the realigned EPI images, and segmentation of intracranial tissue into gray matter, white matter, and cerebrospinal fluid (CSF), which automatically produced volume information for each brain tissue. Gray matter was smoothed with a full-width at half-maximum of 4 mm.
Prior to preprocessing of functional data, (1) the first 10 volumes were removed, (2) slice timing was performed, (3) re-alignment was performed to correct head motion. Based on the recorded motion correction estimates, the subjects with more than 2.0 mm in any direction (x, y, or z) or more than 2.0° at any angle were excluded from the study. The Friston 24-parameter model was used to regress out head motion effects (Yan et al., 2013). Other nuisance variables including white matter signal and CSF were regressed out. (4) Individual functional images were normalized in Montreal Neurological Institute (MNI) space for inter-subject comparison. (5) The resultant images were smoothed with a full-width at half-maximum of 4 mm. (6) De-trending was performed. (7) Data were bandpass filtered (0.01–0.08 Hz).
Based on the pre-processing, DC calculations were performed using DPARSF in a voxel-wise manner (Beucke et al., 2013; Zhou et al., 2016). First, the preprocessed functional data sets were subjected to voxel-based whole-brain correlation analysis. The time course of each voxel within the gray matter mask from each participant was correlated with the time course of every other voxel, which generated a correlation matrix. An undirected adjacency matrix was then obtained by thresholding each correlation at r > 0.25 (Beucke et al., 2013). Then, the DC was computed as the sum of the weights of the significant correlations for each voxel. Finally, by subtracting the mean DC across the entire brain and then dividing by the standard deviation of the whole-brain DC, these individual-level voxel-wise DC maps were standardized into a z-score. A smoothing kernel of 4 mm was applied. Peak MNI coordinates of candidate brain functional hubs identified by DC inter-group analyses and were used as seeds in subsequent analyses, and a sphere ROI was generated for each of the brain sites, the center of which corresponded to the peak voxel with a 6-mm radius. An independent samples t-test of DC maps between groups were performed by REST v.1.8 software3, with age and sex as covariates. Multiple comparisons correction was performed using a Gaussian random field (GRF) at P < 0.001, Z > 3.09 (Hayasaka et al., 2004).
Connectivity among the ROIs was analyzed using REST v.1.8 software, including FC and GCA analyses. For FC, correlation analysis of time course was performed between every two of the seed region for each subject. Fisher’s r-to-z transformation was applied to improve the normality of the FC maps. For GCA, signed-path coefficients between ROIs were computed in multivariate mode for subsequent parametric statistical analyses (Hamilton et al., 2011). The GCA was performed on the gray matter mask using REST tool-box (see text footnote 1). First, we used Granger approach that examined the time lagged effects between two nodes to infer the causal effects between regions (Chen et al., 2011). The signed-path coefficient generated using a time lag order of 1 TR (2 s) was used to estimate the probable excitatory or inhibitory effects of the directed physiological influence (Hamilton et al., 2011). The bivariate GCA accounted for the physiological probability of simultaneously bidirectional influences in the brain (Palaniyappan et al., 2013). In addition, the path coefficients were normally distributed and could be used in parametric statistical analysis for group level inference (Hamilton et al., 2011). Differences in FC and GCA between PSPD patients and HCs were computed with independent samples t-test using SPSS software (SPSS Inc., Chicago, IL, United States), with age and sex as covariates. Multiple comparisons correction was also performed using multiple comparisons by Bonferroni correction (P < 0.05/10 for FC and P < 0.05/20 for GCA).
Correlation Analysis
The relationships among DC, connectivity, clinical parameters, and neuropsychological test scores were examined by partial correlation analyses using SPSS software in all study participants while controlling for the effects of covariates used in the inter-group analyses of DC maps. P < 0.05 was considered statistically significant.
Results
Comparison of Demographic and Clinical Data
There were no significant differences between groups in terms of age, sex, and education. PSPD patients had higher SAS and SDS scores (all P < 0.001) (Table 1).
Table 1. Demographic characteristics and clinical assessment of participants.
DC Analysis
Compared to HCs, lower DC values were observed in the bilateral inferior occipital gyrus (IOG), bilateral calcarine fissure (CAL), and left paracentral lobule (PCL) of PSPD patients (GRF corrected: P < 0.001, Z > 3.09) (Table 2 and Figures 1-1,1-2).
Table 2. Brain regions with significant DC differences between the two groups.
Figure 1. DC value distribution of inter-group comparisons and correlations among the altered DC and PSPD-related parameters (SAS and SDS). (1) The significantly altered DC map in the PSPD group. (2) Comparison of DC value between the two groups (GRF corrected: P < 0.001, Z > 3.09). The color bar denotes the t-value. Error bars define the SEM. (A–J) Correlations among the altered DC and PSPD-related parameters (SAS and SDS). IOG, lateral occipital gyrus; CAL, calcarine fissure; PCL, paracentral gyrus; R, right; L, left; SAS, Self-Rating Anxiety Scale; SDS, Self-Rating Depression Scale.
Figure 2. FC pattern and GCA pattern of inter-group comparisons, and the altered FC and signed-path coefficients significantly correlated with SAS and SDS. (1) Comparison of FC z-scores between the two groups. *P < 0.05, #P < 0.05/10 (Bonferroni correction). (2) Comparison of signed-path coefficients between the two groups. *P < 0.05. Error bars define the SEM. (A–C,E–G) FC between right IOG and left CAL, right IOG and left IOG, left CAL and left IOG correlated with SAS and SDS. (D,H) Signed-path coefficients of the left CAL to the left PCL vs. SAS and SDS. IOG, lateral occipital gyrus; CAL, calcarine fissure; PCL, paracentral gyrus; R, right; L, left; SAS, Self-Rating Anxiety Scale; SDS, Self-Rating Depression Scale.
Connectivity Analysis
Results obtained with the independent t-test indicated that FC values of right IOG–left CAL (P = 0.033), right IOG–right CAL (P = 0.036), right IOG–left IOG (P < 0.001), left CAL–right CAL (P < 0.001), left CAL–left IOG (P = 0.044), left CAL–left PCL (P = 0.012), right CAL–left PCL (P = 0.003), and left IOG–left PCL (P = 0.006) were lower in PSPD patients than in HCs. After multiple comparisons correction, the FC of the right IOG–left IOG, left CAL–right CAL, and right CAL–left PCL was lower in PSPD patients than in HCs (P < 0.05/10) (Figure 2-1).
In the GCA analyses, a negative causal effect from the left CAL to the left paracentral lobule (PCL) and a positive effect from the right CAL to the right IOG were observed in PSPD patients (P < 0.05) (Figure 2-2); however, the independent t-test results were non-significant after multiple comparisons correction (P > 0.05/20) (Figure 2-2).
Correlation Analysis
A significant negative correlation was observed between DC values of the five above-mentioned brain regions (bilateral IOG, bilateral CAL, and left PCL) and SAS as well as SDS scores in all participants (P < 0.05) (Figures 1A–J).
FC values of right IOG–left CAL (r = -0.369, P = 0.025; r = −0.342, P = 0.038), right IOG–left IOG (r = −0.571, P < 0.001; r = −0.694, P < 0.001), and left CAL–left IOG (r = −0.343, P = 0.038; r = −0.345, P = 0.037) were negatively correlated with SAS and SDS scores (Figures 2A–C,E–G) in all participants. Signed-path coefficients from the left CAL to the left PCL were also negatively correlated with SAS and SDS scores (r = −0.343, P = 0.035; r = −0.415, P = 0.010; Figures 2D,H).
Discussion
In this study we used a novel approach combining DC, FC, and GCA to investigate changes in PSPD patient multiple brain regions and the whole-brain functional network connectome. We found that PSPD patients showed abnormal network DC in multiple brain regions as well as abnormal FC between these regions as expected. The GCA analysis further revealed disordered connectivity from the left CAL to the left PCL and from the right CAL to the right IOG. Moreover, significant negative correlations were found between all abnormal network DC values, most abnormal FC values, signed-path coefficients from the left CAL to the left PCL, and SAS as well as SDS scores. These findings provide novel insight into the large-scale functional reorganization that occurs in the PSPD patient brain.
In humans, the IOG (BA19) is a visual association area that functions in feature extraction, shape recognition, attention, and multimodal integration, while the CAL harbors the primary visual cortex (V1) (Thiebaut de Schotten et al., 2014). We found that PSPD patients showed abnormal network DC in the bilateral IOG and bilateral CAL and lower FC values of bilateral CAL–right IOG, bilateral IOG, bilateral CAL, and left CAL–left IOG. Furthermore, GCA showed a disruption of direct connectivity from the right CAL to the right IOG, which is where abnormal FC was detected. Thus, PSPD patients exhibit altered brain functional hubs and connectivity in the occipital lobe, which is presumed to be is related to visual information processing (Sveinbjornsdottir and Duncan, 1993). However, the occipital cortex is connected to structures associated with inhibition of the descending pain pathway in the rats animal experiments (Reis et al., 2010). Several studies have demonstrated decreased signals in the occipital lobes of patients with chronic pain disorder by positron emission tomography, electroencephalogram, and rs-fMRI (Cauda et al., 2009; Karibe et al., 2010; Klug et al., 2011). Decreased ReHo in the occipital lobe was found to be negatively correlated with migraine duration (Zhao et al., 2013) and was observed in PSPD patients (Huang et al., 2016). GCA in current study showed the disruption of direct connectivity from the right CAL to the right IOG also showed less neural impulses in the occipital lobe of PSPD patients than of HC. We speculate that altered functional hubs and connectivity in this brain region contributes to the pathological impairment of sensation and pain-related emotion in PSPD. However, the relationship between abnormal activation of the occipital lobe and pain-related brain dysfunction remains to be elucidated.
The PCL controls motor and sensory innervation and plays an important role in somatosensation (Jellinger, 2008). Many studies have reported that the pain pathways of the human brain are widely distributed and encompass the somatosensory cortex, which is a sensory/motivational association area and involved in affective/discriminative aspects of pain (Xie et al., 2009; Seifert and Maihofner, 2011; Quintero, 2013). Neuroscientific research has indicated that the effects of expectancy on the subjective pain experience are paralleled by changes in somatosensory areas which are afferent nociceptive brain areas, and that this is partly mediated by descending pain modulatory circuits (Peters, 2015). A meta-analysis of PET, fMRI, EEG, and MEG studies provides clarity regarding the regions including primary and secondary somatosensory, insular, anterior cingulate and thalamus found active during an acute pain experience, investigating that these areas are the fundamental core network of human nociceptive processing (Apkarian et al., 2005). PSPD is a mental disorder with pain that is unrelated to somatic injury but can cause significant somatosensory and emotional disturbance (Egloff et al., 2014). Decreased ReHo were found in the bilateral primary somatosensory cortex of PSPD patients (Huang et al., 2016) and altered co-activation in the primary somatosensory cortex (Zhao et al., 2017). We found that PSPD patients had abnormal network DC in the PCL as well as lower FC values of bilateral CAL–left PCL and left IOG–left PCL. Our results suggest that pain-related networks are in fact highly interactive and disrupt intra-network FC within sensory and cognitive systems.
However, GCA showed disordered connectivity from the left CAL to the left PCL, which is where abnormal FC was detected. GCA result revealed more neural impulses or output information from visual cortex to somatosensory cortex but less from somatosensory cortex to visual cortex of PSPD patients than of HC. The combined effect might cause the lower FC between the left CAL and the left PCL. But inconsistent between the undirected and directed connectivity may result from the distinction of mathematical theories (Seth et al., 2015). So this phenomenon requires a further investigation.
The significant negative correlations observed between all abnormal network DC values, most abnormal FC values, signed-path coefficients from the left CAL to the left PCL, and SAS as well as SDS scores. Our results suggest that emotion perception could be affected by sensory/discriminative pain processing, consistent with the notion that various states of patients could influence on pain perception via sensory systems, such as anxiety, sadness, and depression (Schonenberg et al., 2014). We speculate that these characterization on large-scale functional networks and correlation could explain pain-related emotional disturbance in PSPD. Thus, psychological interventions that enhance emotional awareness may be beneficial for PSPD patients.
This study had a few limitations. Firstly, since the population size was relatively small, inter-group comparisons in GCA were non-significant after stringent multiple corrections, which affected not only the statistical power but also the interpretation of the results, especially the sample size, age distribution, and gender bias, which may be the reason that the relationship between DC, FC and clinical assessment within PSPD group wasn’t found. So the result needed to be replicated in larger samples. What’s more, one of the main limitation is the relatively long TR = 2 s which is sustantilly longer than typical inter-neuron delays, although most functional MRI study use this TR. EEG could offer millisecond-range temporal resolution which is not possible with MRI (Britz et al., 2010). EEG is a widely utilized technique to measure oscillatory activity in different frequency bands originated by neuronal assemblies, and to link these oscillations with specific behavioral states (Britz et al., 2010; Xu et al., 2010, 2015). There is accumulating evidence that the spatial distribution of scale-free properties was compatible between EEG and fMRI (Lei et al., 2011; Xu et al., 2014). We are also considering the application of using simultaneous EEG-fMRI to get a more complete characterization of scale-free organization of brain in PSPD.
Conclusion
We took the novel approach of using a combined DC, FC, and GCA approach to investigate the changes in the PSPD patients brain functional hubs and their connectivity, and found altered large-scale functional networks involved occipital lobe and somatosensory cortices, which may mirror the underlying neural mechanism of the dysfunction of sensation and pain-related emotional disturbances in PSPD patients.
Author Contributions
QL conceived and executed the project and wrote the first draft of the manuscript. X-cZ performed the MRI data analysis. X-MJ contributed to the participant recruitment and clinical data acquisition. Z-hZ and X-fH contributed to the project conception, manuscript review, and critique.
Funding
This study was funded by grants from the National Natural Science Foundation of China (81701664) and the Technology Innovation Program in Southwest Hospital (SWH2016JCYB-30). The clinical novel technology program of military medicine and injury treatment of Southwest Hospital (SWH2016JSZD-02).
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 International Science Editing (http://www.internationalscienceediting.com) for editing this manuscript.
Footnotes
References
American Psychiatric Association [APA] (2013). Diagnostic and Statistical Manual of Mental Disorders (DSM-5). Washington, DC: American Psychiatric Pub.
Apkarian, A. V., Bushnell, M. C., Treede, R. D., and Zubieta, J. K. (2005). Human brain mechanisms of pain perception and regulation in health and disease. Eur. J. Pain 9, 463–484. doi: 10.1016/j.ejpain.2004.11.001
Beucke, J. C., Sepulcre, J., Talukdar, T., Linnman, C., Zschenderlein, K., Endrass, T., et al. (2013). Abnormally high degree connectivity of the orbitofrontal cortex in obsessive-compulsive disorder. JAMA Psychiatry 70, 619–629. doi: 10.1001/jamapsychiatry.2013.173
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
Britz, J., Van De Ville, D., and Michel, C. M. (2010). BOLD correlates of EEG topography reveal rapid resting-state network dynamics. NeuroImage 52, 1162–1170. doi: 10.1016/j.neuroimage.2010.02.052
Burba, B., Oswald, R., Grigaliunien, V., Neverauskiene, S., Jankuviene, O., and Chue, P. (2006). A controlled study of alexithymia in adolescent patients with persistent somatoform pain disorder. Can. J. Psychiatry 51, 468–471. doi: 10.1177/070674370605100709
Cauda, F., Sacco, K., Duca, S., Cocito, D., D’Agata, F., Geminiani, G. C., et al. (2009). Altered resting state in diabetic neuropathic pain. PLoS One 4:e4542. doi: 10.1371/journal.pone.0004542
Chao-Gan, Y., and Yu-Feng, Z. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4:13. doi: 10.3389/fnsys.2010.00013
Chen, G., Glen, D. R., Saad, Z. S., Hamilton, J. P., Thomason, M. E., Gotlib, I. H., et al. (2011). Vector autoregression, structural equation modeling, and their synthesis in neuroimaging data analysis. Comput. Biol. Med. 41, 1142–1155. doi: 10.1016/j.compbiomed.2011.09.004
Cole, M. W., Pathak, S., and Schneider, W. (2010). Identifying the brain’s most globally connected regions. NeuroImage 49, 3132–3148. doi: 10.1016/j.neuroimage.2009.11.001
Di, M. A., Zuo, X. N., Kelly, C., Grzadzinski, R., Mennes, M., Schvarcz, A., et al. (2013). Shared and distinct intrinsic functional network centrality in autism and attention-deficit/hyperactivity disorder. Biol. Psychiatry 74, 623–632. doi: 10.1016/j.biopsych.2013.02.011
Egloff, N., Camara, R. J., von Kanel, R., Klingler, N., Marti, E., and Ferrari, M. L. (2014). Hypersensitivity and hyperalgesia in somatoform pain disorders. Gen. Hosp. Psychiatry 36, 284–290. doi: 10.1016/j.genhosppsych.2014.01.011
Fransson, P., Åden, U., Blennow, M., and Lagercrantz, H. (2011). The functional architecture of the infant brain as revealed by resting-state FMRI. Cereb. Cortex 21, 145–154. doi: 10.1093/cercor/bhq071
Hamilton, J. P., Chen, G., Thomason, M. E., Schwartz, M. E., and Gotlib, I. H. (2011). investigating neural primacy in major depressive disorder: multivariate granger causality analysis of resting-state FMRI time-series data. Mol. Psychiatry 16, 763–772. doi: 10.1038/mp.2010.46
Hayasaka, S., Phan, K. L., Liberzon, I., Worsley, K. J., and Nichols, T. E. (2004). Nonstationary cluster-size inference with random field and permutation methods. NeuroImage 22, 676–687. doi: 10.1016/j.neuroimage.2004.01.041
Huang, T., Zhao, Z., Yan, C., Lu, J., Li, X., Tang, C., et al. (2016). Altered spontaneous activity in patients with persistent somatoform pain disorder revealed by regional homogeneity. PLoS One 11:e0151360. doi: 10.1371/journal.pone.0151360
Karibe, H., Arakawa, R., Tateno, A., Mizumura, S., Okada, T., Ishii, T., et al. (2010). Regional cerebral blood flow in patients with orally localized somatoform pain disorder: a single photon emission computed tomography study. Psychiatry Clin. Neurosci. 64, 476–482. doi: 10.1111/j.1440-1819.2010.02119.x
Klug, S., Stefanie, K., Anderer, P., Peter, A., Saletuzyhlarz, G., Gerda, S. Z., et al. (2011). Dysfunctional pain modulation in somatoform pain disorder patients. Eur. Arch. Psychiatry Clin. Neurosci. 261, 267–275. doi: 10.1007/s00406-010-0148-4
Lei, X., Xu, P., Luo, C., Zhao, J., Zhou, D., and Yao, D. (2011). fMRI functional networks for EEG source imaging. Hum. Brain Mapp. 32, 1141–1160. doi: 10.1002/hbm.21098
Luo, Y., Yan, C., Huang, T., Fan, M., Liu, L., Zhao, Z., et al. (2016). Altered neural correlates of emotion associated pain processing in persistent somatoform pain disorder: an fMRI study. Pain Pract. 16, 969–979. doi: 10.1111/papr.12358
Luo, Y. L., Heeramun-Aubeeluck, A., Huang, X., Ye, G., Wu, H., Sun, L., et al. (2014). Factors influencing quality of life in chinese patients with persistent somatoform pain disorder. Psychol. Health Med. 19, 744–752. doi: 10.1080/13548506.2013.878804
Palaniyappan, L., Simmonite, M., White, T. P., Liddle, E. B., and Liddle, P. F. (2013). Neural primacy of the salience processing system in schizophrenia. Neuron 79, 814–828. doi: 10.1016/j.neuron.2013.06.027
Peters, M. L. (2015). Emotional and cognitive influences on pain experience. Mod. Trends Pharmacopsychiatry 30, 138–152. doi: 10.1159/000435938
Reips, U. D., and Funke, F. (2008). Interval-level measurement with visual analogue scales in internet-based research: vas generator. Behav. Res. Methods 40, 699–704. doi: 10.3758/brm.40.3.699
Reis, G. M., Dias, Q. M., Silveira, J. W. S., Vecchio, F. D., Garcia-Cairasco, N., and Prado, W. A. (2010). Antinociceptive effect of stimulating the occipital or retrosplenial cortex in rats. J. Pain 11, 1015–1026. doi: 10.1016/j.jpain.2010.01.269
Schonenberg, M., Mares, L., Smolka, R., Jusyte, A., Zipfel, S., and Hautzinger, M. (2014). Facial affect perception and mentalizing abilities in female patients with persistent somatoform pain disorder. Eur. J. Pain 18, 949–956. doi: 10.1002/j.1532-2149.2013.00440.x
Seifert, F., and Maihofner, C. (2011). Functional and structural imaging of pain-induced neuroplasticity. Curr. Opin. Anaesthesiol. 24, 515–523. doi: 10.1097/ACO.0b013e32834a1079
Seth, A. K., Barrett, A. B., and Barnett, L. (2015). Granger causality analysis in neuroscience and neuroimaging. J. Neurosci. 35, 3293–3297. doi: 10.1523/JNEUROSCI.4399-14.2015
Smith, G. C., Clarke, D. M., Handrinos, D., Dunsis, A., and Mckenzie, D. P. (2000). Consultation-Liaison psychiatrists’ management of somatoform disorders. Psychosomatics 41, 481–489. doi: 10.1176/appi.psy.41.6.481
Sveinbjornsdottir, S., and Duncan, J. S. (1993). Parietal and occipital lobe epilepsy: a review. Epilepsia 34, 493–521. doi: 10.1111/j.1528-1157.1993.tb02590.x
Thiebaut de Schotten, M., Urbanski, M., Valabregue, R., Bayle, D. J., and Volle, E. (2014). Subdivision of the occipital lobes: an anatomical and functional MRI connectivity study. Cortex 56, 121–137. doi: 10.1016/j.cortex.2012.12.007
Tomasi, D., and Volkow, N. D. (2011). Functional connectivity hubs in the human brain. NeuroImage 57, 908–917. doi: 10.1016/j.neuroimage.2011.05.024
World Health Organization [WHO] (1993). The ICD-10 Classification of Mental and Behavioural Disorders: Diagnostic Criteria for Research. Geneva: World Health Organization.
Xie, Y. F., Huo, F. Q., and Tang, J. S. (2009). Cerebral cortex modulation of pain. Acta Pharmacol. Sin. 30, 31–41. doi: 10.1038/aps.2008.14
Xu, L., Chuan, Q., Peng, X., and Dezhong, Y. (2010). A parallel framework for simultaneous EEG/fMRI analysis: methodology and simulation. NeuroImage 52, 1123–1134. doi: 10.1016/j.neuroimage.2010.01.024
Xu, L., Wang, Y., Hong, Y., and Chen, A. (2015). Brain scale-free properties in awake rest and NREM Sleep: a simultaneous EEG/fMRI study. Brain Topogr. 28, 292–304. doi: 10.1007/s10548-014-0399-x
Xu, L., Yulin, W., Hong, Y., and Dante, M. (2014). Neuronal oscillations and functional interactions between resting state networks. Hum. Brain Mapp. 35, 3517–3528. doi: 10.1002/hbm.22418
Yan, C. G., Cheung, B., Kelly, C., Colcombe, S., Craddock, R. C., Di Martino, A., et al. (2013). A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics. NeuroImage 76, 183–201. doi: 10.1016/j.neuroimage.2013.03.004
Zhao, L., Liu, J., Dong, X., Peng, Y., Yuan, K., Wu, F., et al. (2013). Alterations in regional homogeneity assessed by fMRI in patients with migraine without aura stratified by disease duration. J. Headache Pain 14:85. doi: 10.1186/1129-2377-14-85
Zhao, Z., Huang, T., Tang, C., Ni, K., Pan, X., Yan, C., et al. (2017). Altered resting-state intra- and inter- network functional connectivity in patients with persistent somatoform pain disorder. PLoS One 12:e0176494. doi: 10.1371/journal.pone.0176494
Zhou, C., Hu, X., Hu, J., Liang, M., Yin, X., Lin, C., et al. (2016). Altered brain network in amyotrophic lateral sclerosis: a resting graph theory-based network study at voxel-wise level. Front. Neurosci. 10:e42816. doi: 10.3389/fnins.2016.00204
Zung, W. W. (1971). A rating instrument for anxiety disorders. Psychosomatics 12, 371–379. doi: 10.1016/s0033-3182(71)71479-0
Zung, W. W. (1974). The measurement of affects: depression and anxiety. Mod. Probl. Pharmacopsychiatry 7, 170–188. doi: 10.1159/000395075
Zung, W. W., Richards, C. B., and Short, M. J. (1965). Self-rating depression scale in an outpatient clinic. further validation of the sds. Arch. Gen. Psychiatry 13, 508–515.
Keywords: persistent somatoform pain disorder, emotional disturbance, resting-state functional MRI, degree centrality, Granger causality analysis
Citation: Liu Q, Zeng X-c, Jiang X-M, Zhou Z-h and Hu X-f (2019) Altered Brain Functional Hubs and Connectivity Underlie Persistent Somatoform Pain Disorder. Front. Neurosci. 13:415. doi: 10.3389/fnins.2019.00415
Received: 09 December 2018; Accepted: 11 April 2019;
Published: 30 April 2019.
Edited by:
Quanying Liu, California Institute of Technology, United StatesReviewed by:
Xiaoyun Liang, Australian Catholic University, AustraliaYali Jiang, Central South University, China
Copyright © 2019 Liu, Zeng, Jiang, Zhou and Hu. 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: Zhen-hua Zhou, ZXhwbG9pdGVyMDAxQDEyNi5jb20= Xiao-fei Hu, aGFycnl6b25ldG1tdUAxNjMuY29t
†These authors have contributed equally to this work