Front. Neurosci.Frontiers in NeuroscienceFront. Neurosci.1662-453XFrontiers Media S.A.10.3389/fnins.2018.00900NeuroscienceOriginal ResearchRegional Homogeneity and Multivariate Pattern Analysis of Cervical Spondylosis Neck Pain and the Modulation Effect of TreatmentChenJun1*WangZengjian2TuYiheng3LiuXian1JorgensonKristen3YeGuoxi1LinChenlin1LiuJianhua4*ParkJoel3LangCourtney3LiuBo1*KongJian31Department of Radiology, Guangdong Provincial Hospital of Chinese Medicine, Guangzhou, China2Department of Maternal and Child Health, School of Public Health, Sun Yat-sen University, Guangzhou, China3Department of Psychiatry, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States4Key Laboratory for Studying Regularities and Mechanism of Acu-moxibustion, Department of Acu-moxibustion, Guangdong Provincial Hospital of Chinese Medicine, Guangzhou, China
Edited by: Xi-Nian Zuo, Institute of Psychology (CAS), China
Reviewed by: Lijun Bai, Xi’an Jiaotong University, China; Jiliang Fang, Guang’anmen Hospital, China Academy of Chinese Medical Sciences, China
*Correspondence: Jun Chen, junesums@qq.com Jianhua Liu, jyhf08@sina.com Bo Liu, lbgdhtcm@163.com
This article was submitted to Brain Imaging Methods, a section of the journal Frontiers in Neuroscience
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.
Objects: We investigated brain functional alteration in patients with chronic cervical spondylosis neck pain (CSNP) compared to healthy controls (HCs) and the effect of intervention.
Methods: 104 CSNP patients and 96 matched HCs were recruited. Patients received 4 weeks of treatment. Resting-state fMRI and Northwick Park Neck Pain Questionnaire (NPQ) were collected before and after treatment. Resting state regional homogeneity (rs-ReHo) and multivariate pattern analysis (MVPA) were applied to (1) investigate rs-ReHo differences between CSNP patients and controls and the effect of longitudinal treatment and (2) classify CSNP patients from HCs and predict clinical outcomes before treatment using MVPA.
Results: We found that (1) CSNP patients showed decreased rs-ReHo in the left sensorimotor cortex and right temporo-parietal junction (rTPJ), and rs-ReHo at the rTPJ significantly increased after treatment; (2) rs-ReHo at rTPJ was associated with NPQ at baseline, and pre- and post-treatment rs-ReHo changes at rTPJ were associated with NPQ changes in CSNP patients; and (3) MVPA could discriminate CSNP patients from HCs with 72% accuracy and predict clinical outcomes with a mean absolute error of 19.6%.
Conclusion: CSNP patients are associated with dysfunction of the rTPJ and sensorimotor area.
Significance: rTPJ plays on important role in the pathophysiology and development of CSNP.
CSNP patients showed decreased ReHo in the left S1/M1 and right TPJ.
ReHo at rTPJ increased after treatment in CSNP patients.
ReHo at rTPJ and S1/M1 can classify CSNP patients and HCs.
ReHo at baseline can predict NPQ changes after treatment.
Introduction
Cervical spondylosis (CS) is a common age-related chronic disease characterized by stiffness and neck and upper back pain. Chronic CS neck pain (CSNP) has become an important public health and social issue because of its high prevalence, unsatisfactory treatment options, large medical burden, and reduction in quality of life (Binder, 2007; Hoy et al., 2010; Baime, 2016). Recent studies have demonstrated that patients with cervical spondylotic myelopathy (CSM), the later stage of CS, have cerebral cortex reorganization as a result of spinal cord compression (Holly et al., 2007; Dong et al., 2008; Holly, 2009; Duggal et al., 2010). However, the characteristics of cortical activity in CSNP patients without spinal cord compression have not been reported.
Previous studies showed that patients with spinal cord compression in CSM are associated with increased activation in the primary motor and premotor cortices and activation loss in the sensory cortex (Kowalczyk et al., 2012; Zhou et al., 2014; Zhou F. et al., 2015; Zhou F. Q. et al., 2015). Single-photon emission computed tomography (SPECT) and positron emission tomography (PET) studies found significant abnormal hypoperfusion and hypometabolism in brain regions such as the bilateral parieto-occipital cortex in chronic neck pain patients with whiplash syndrome (Otte et al., 1995, 1996, 1997a,b; Bicik et al., 1998; Lorberboym et al., 2002; Guez, 2006; Linnman et al., 2009). Neuroimaging studies also found changes in brain structure and function in chronic back pain (CBP) patients (Wand et al., 2011; Wasan et al., 2011; Kong et al., 2013; Loggia et al., 2013; Mao et al., 2014; Yu et al., 2014a; Pijnenburg et al., 2015; Hotz-Boendermaker et al., 2016).
Recent studies also found improvements in brain functional activity and connectivity in patients with spinal cord decompression following intervention (Holly et al., 2007; Dong et al., 2008; Holly, 2009; Duggal et al., 2010; Tam et al., 2010; Hrabalek et al., 2014; Green et al., 2015; Hrabalek et al., 2015; Tan et al., 2015; Bhagavatula et al., 2016). Dong and colleagues found that abnormalities in the cortical sensorimotor recruitment pattern in CSM patients gradually disappeared and were accompanied by behavioral improvements after surgery (Dong et al., 2008). Effective longitudinal treatment can also modulate the aberrant functional neural activity and brain structure of chronic low back pain patients (Lorimer Moseley, 2005; Baliki et al., 2008; Seminowicz et al., 2011; Younger et al., 2011; Li et al., 2014; Ceko et al., 2015; Louw et al., 2015; Shi et al., 2015; Vrana et al., 2015; Braden et al., 2016). Nonetheless, there is a lack of research evaluating the modulation effect of effective treatment on the brain activity of CSNP patients.
In recent years, machine learning has gained popularity in the translational neuroimaging community for it allows researchers to identify brain signatures of patients with psychiatric and neurologic diseases and make accurate classifications (Arbabshirani et al., 2017). A widely applied machine learning approach, multivariate pattern analysis (MVPA) uses pattern classifiers to capture the information encoded by spatially correlated voxels in fMRI and has proven to be more sensitive to the functional organization of the cortex (e.g., general linear model; Haxby, 2012). This technique has been used to explore the neurophysiology of chronic pain disorders, including migraine (Chong et al., 2017), chronic low back pain (Ung et al., 2014) and fibromyalgia (Lopez-Sola et al., 2017).
In this study, we explored (1) the differences in resting state brain activity between CSNP patients and matched healthy controls (HCs), (2) how longitudinal treatment can modulate brain abnormalities in CSNP patients using regional homogeneity (ReHo), and (3) if MVPA approaches can classify CSNP patients and HCs and predict clinical outcomes before treatment. We hypothesized that (1) CSNP patients would be associated with altered activity in specific brain regions (e.g., sensorimotor cortex or other regions) compared with HCs, (2) symptom reduction after an effective treatment could modulate the altered brain activities of CSNP patients, and (3) CSNP patients and HCs could be discriminated using MVPA.
Materials and MethodsParticipants
This study was approved by the Ethics Committee at the Second Affiliated Hospital of Guangzhou University of Chinese Medicine. The experiment was performed in accordance with approved guidelines. Participants were recruited from the outpatient clinics of the Second Affiliated Hospital, Guangzhou University of Traditional Chinese Medicine. This study was registered on http://www.chictr.org.cn/index.aspx (ChiCTR-IPR-14005282). Written informed consent was obtained from each participant. Subjects were recruited from July 2010 to December 2012.
The inclusion criteria were as follows: (1) a confirmed diagnosis of CSNP in accordance with the diagnostic criteria published by the Chinese Association of Rehabilitation (2010) and in accordance with the International Classification of Diseases, 10th edition (ICD-10) codes1; (2) male or female 18–35 years of age; (3) no signs of musculoskeletal system abnormalities upon physical examination; (4) not receiving any treatment within the last 7 days; and (5) informed consent obtained. Patients were excluded if they met any of the following criteria: (1) history of neck trauma; (2) vertebral body or spinal canal cancer, tuberculosis, or severe osteoporosis; (3) history of neck surgery or presence of congenital malformation of the cervical vertebrae; (4) pregnant or breastfeeding (females); and (5) presence of a severe systemic disease such as tumors, diabetes mellitus, kidney disease, or digestive system disease.
Additionally, pain-free age and gender-matched individuals were recruited for this study as HCs. Each subject underwent a medical history evaluation and physical examination.
Interventions
All CSNP patients received acupuncture treatments for one month at the Second Affiliated Hospital of Guangzhou University of Chinese Medicine. The CSNP patients were assigned to one of five acupuncture treatment groups based on the acupuncture points applied. The acupoints selected were bilateral Bailao (EX-HN15) in acupuncture group 1 (AG1), bilateral EX-HN15 and Hegu (LI4) in AG2, bilateral EX-HN15 and Zhongzhu (TE3) in AG3, bilateral LI4 in AG4, and bilateral Zusanli (ST36) in AG5.
Each patient received about 30 min of acupuncture treatment approximately once every three days for a total of 10 sessions over a 4-week period. The needles were inserted to a depth of 20–30 mm, perpendicular to the surface of the skin. Electroacupuncture was performed to stimulate the relevant acupoints for 30 min using a continuous-wave, 1-Hz frequency, and a comfortable strength until deqi sensation (soreness, numbness, distention, and heaviness) was obtained.
The Northwick Park Neck Pain Questionnaire (NPQ) was administered to CSNP patients to assess neck pain and disability (Leak et al., 1994). The NPQ scores were assessed before the first treatment as baseline data and again after the final treatment by the patients themselves following thorough instruction from the researchers.
Image Acquisition
All participants underwent a structural and functional MRI scan using a 3.0T MR system (Magnetom Verio, Siemens, Germany) at the Department of Radiology at the Second Affiliated Hospital of Guangzhou University of Chinese Medicine. Cushions were used to support participants’ heads within the coil to prevent excessive head motion. Participants were asked to stay still with their eyes closed. All images were acquired parallel to the anterior-commissure-posterior-commissure line with an auto-align technique. The functional data were collected using an echo-planar imaging sequence: 31 axial slices; repetition time (TR) = 2000 ms; echo time (TE) = 30 ms; slice thickness = 3.5 mm; gap = 0.35 mm; flip angle (FA) = 90°; matrix = 64 × 64; field of view (FOV) = 224 mm × 224 mm, 240 time points. 3D structural images were acquired using a T1-weighted MP-RAGE sequence: 176 sagittal slices; TR = 1900 ms; TE = 2.27 ms; inversion time = 900 ms; slice thickness = 1.0 mm; no gap; FA = 9°; matrix = 256 × 256; FOV = 256 mm × 256 mm.
Data Preprocessing
Data preprocessing was carried out using Data Processing Assistant for Resting-State fMRI (DPARSF2), which is based on Statistical Parametric Mapping (SPM3) and the toolbox for Data Processing and Analysis of Brain Imaging (DPABI4). The preprocessing procedures included: time alignment across slices, motion correction, within-subject registration between T1 and EPI images, T1 segmentation, and normalization of the BOLD fMRI datasets to register them to Montreal Neurologic Institute (MNI) space with voxels resampled at 3 mm × 3 mm × 3 mm. Participants were excluded if head motion exceeded 1.5 mm and 1.5°. Finally, the waveform of each voxel was passed through a band-pass filter (0.01–0.08 Hz) to reduce the effects of low-frequency drift and high-frequency physiological noise. Finally, several nuisance signals were regressed out from each voxel’s time course, including 24-parameter head-motion profiles, mean white matter (WM) and cerebrospinal fluid (CSF) time series.
Regional Homogeneity (ReHo)
ReHo was calculated by the Kendall’s coefficient of concordance (KCC) (Kendall, 1990) using the REST toolkit5. Individual resting state ReHo (rs-ReHo) maps were generated by assigning each voxel a value corresponding to the KCC of its time series with its 26 nearest neighboring voxels (Zang et al., 2004). The individual rs-ReHo maps were standardized by their own mean KCC. Then, a Gaussian kernel with a full-width at half-maximum of 4 mm was used to smooth the images in order to reduce noise and residual differences.
Group analysis was calculated with a random effects model using the Data Processing and Analysis of Brain Imaging toolbox (DPABI6). We first compared the rs-ReHo differences between CSNP patients and HCs using two-sample t-tests. Then, a partial correlation analysis was used to investigate the relationship between the rs-ReHo of these aberrant brain regions and NPQ scores in the CSNP patients, controlling for age and gender as covariates using SPSS.
In order to explore whether clinical improvement after an effective treatment could affect the brain’s intrinsic neural activity in those aforementioned aberrant brain regions, we defined all brain regions that showed significant rs-ReHo changes between CSNP patients and HCs as regions of interest (ROI) and compared the pre- and post-treatment rs-ReHo changes in CSNP patients using a paired t-test within the ROIs as well as the whole brain. Then, partial correlation analysis was used to investigate the relationship between pre- and post-treatment rs-ReHo changes within the ROIs and corresponding NPQ score changes in the CSNP patients, controlling for age and gender as covariates using SPSS. A threshold of voxel-wise p < 0.005 and p < 0.05 family-wise error corrected at cluster level was applied for whole brain analyses, and a small volume correction method (voxel-wise p < 0.005 and p < 0.05 FWE corrected) was used for all ROIs (derived from a comparison between the CSNP patients and healthy controls before and after acupuncture treatments).
Multivariate Pattern Analysis (MVPA)
The ROIs defined in the ReHo analyses were used for the MVPA analyses with two objectives: (1) classify CSNP patients and HCs and (2) predict treatment responses using baseline (pre-treatment) ReHo values. To avoid the risk of overfitting, all analyses were based on fivefold cross-validation (CV) (Pereira et al., 2009).
In the first step, machine learning models were trained to classify CSNP patients and HCs using ReHo values within the ROIs. A support vector machine (SVM) classifier was used and the implementation of SVM was based on LIBSVM (Chang and Lin, 2011). To quantify the performance of the SVM classifier, classification accuracy, sensitivity and specificity were calculated. Here, sensitivity and specificity represent the proportion of patients correctly classified and the proportion of HCs correctly classified, respectively. To further assess the performance of the classifier and evaluate the significance of classification accuracy, we ran permutation testing. In each test, we randomly permuted the class labels of the data prior to training. Fivefold CV was then performed on the permuted dataset and the procedure was repeated 10,000 times. If the classifier trained on real class labels had an accuracy exceeding the 95% confidence interval generated from the accuracies of the classifiers trained on randomly relabeled class labels, this classifier was considered to be well-performing.
In the second step, we aimed to predict treatment responses (percentage change of NPQ scores) using baseline ReHo values. Here, we used a novel fMRI prediction approach, namely sliced inverse regression (SIR), to extract the most predictive information from ROIs and subsequently predict responses. SIR is powerful because it can work efficiently regardless of the linearity of the relationship between fMRI data and their labels. The details of this implementation of SIR-based fMRI prediction can be found in Tu et al. (2017). Briefly, in the first step, we performed SIR on the features Φ ∈ Rn×p and corresponding behavioral data ψ ∈ Rn (i.e., percentage changes of NPQ), resulting in the effective dimension reduction (e.d.r.) direction β ∈ Rp×K and SIR-derived feature set S = Φβ ∈ Rn×K. In the second step, we used a support vector regression (SVR) ψ = C(Stα;t = 1, … ,n) with α ∈ RK to estimate the predicted value and the contribution of each variable in the SIR reduced space to predict ψ. In the third step, we transformed α back to the original space by multiplying β, resulting in the weights for all features 𝜃 = βα. To quantify the performance of prediction, we used the prediction-outcome correlation which is defined as the correlation between the actual and predicted differences of NPQ scores, as well as mean absolute error (MAE) (Tu et al., 2016), which is defined as
MAE=1N∑n=1N|ψn−ψ^n|.
where ψn and ψ^n denote, respectively, the actual and predicted differences of NPQ scores (1-Post NPQ/Pre NPQ) for patient n, and N is the number of CSNP patients. The number of slices for SIR was set to 5 and the number of directions was set to 4. Performance measures were assessed using permutation testing (as mentioned in the classification session).
Results
A total of 104 right-handed CSNP patients (45 males and 59 females, age 24.90 ± 1.98 years, neck pain duration 36.49 ± 24.93 months) and 96 HCs (50 males and 46 females, age 24.80 ± 1.52 years) were recruited. There were no significant differences in age and gender between CSNP patients and HCs (P > 0.05).
Compared with baseline scores, all treatment groups (AG1, AG2, AG3, AG4, and AG5) showed improvement in NPQ scores (P < 0.01). We found no significant differences among AG1, AG2, AG3, AG4, and AG5 in NPQ score improvement (P > 0.05). The average NPQ score across all patients was 24.1 ± 9.2 before treatment and 12.7 ± 9.0 at the end of treatment. There was a significant decrease in NPQ scores after 1 month of treatment (P < 0.01).
Based on these clinical findings, we merged all acupuncture treatment groups (AAG) in the rs-ReHo analyses to investigate how an effective treatment can modulate rs-ReHo.
ReHo Results
To explore the neural pathophysiology of neck pain, we first compared all CSNP patients to HCs using a two-sample t-test. Compared to HCs, CSNP patients showed decreased rs-ReHo in the right temporo-parietal junction (rTPJ) and left sensorimotor cortex (left post-central gyrus and precentral gyrus) (Figure 1 and Table 1). Regression analysis indicated that the rs-ReHo in the rTPJ was positively associated with NPQ scores at baseline in CSNP patients (Figure 1).
Altered resting state ReHo in CSNP patients and the modulation effect of clinical treament. The left sensorimotor cortex and rTPJ showed reduced rs-ReHo in CSNP patients compared to healthy controls. The rs-ReHo value in rTPJ was associated with neck pain intensity (NPQ) at baseline in CSNP patients.
Brain regions with significant ReHo difference between CSNP patients and health controls. rTPJ, right temporo-parietal junction.
Brain region
Brodmann area
Voxels
Coordinates (MNI space)
Peak value
rTPJ/temporal cortex
22/13/21
117
60 -48 6
-3.950
left sensorimotor cortex
3/4/2
84
-51 -27 42
-4.378
We then compared post- and pre-treatment rs-ReHo differences using paired t-tests. Results revealed that after longitudinal clinical treatment, CSNP patients showed increased rs-ReHo in the rTPJ (peak 60, -48, 6, 13 voxels). There were no significant differences between the pre- and post-treatment rs-ReHo in other ROIs and non-ROIs at the threshold we set. Regression analysis indicated that the pre- and post-treatment rs-ReHo changes at the rTPJ were positively associated with corresponding NPQ changes in CSNP patients.
MVPA Results
Classification accuracy for separating CSNP patients and HCs is summarized in Figure 2, upper panel. Using the ReHo values within the rTPJ, left postcentral gyrus, and precentral gyrus, the classifier achieved an accuracy of 72.0% (76.0% for CSNP and 67.7% for HCs; p < 0.0001, permutation testing).
The MVPA results. Upper panel: The classification accuracy of discriminating CSNP patients and HCs. Lower panel: The prediction-outcome correlation of predicting treatment responses using baseline fMRI data.
Similar to a previous study (Doehrmann et al., 2013), two subjects were removed as outliers from the prediction analysis since their NPQ changes were more than 3 SDs higher than the average of our samples. The correlation between actual and predicted treatment response was 0.30 (p = 0.002; Figure 2, lower panel). We obtained an MAE of 19.6% for predicting treatment responses (p = 0.004, permutation testing).
Discussion
In this study, we investigated rs-ReHo differences between CSNP patients and HCs as well as how symptom relief after treatment can modulate rs-ReHo. We found that CSNP patients were associated with reduced rs-ReHo in brain regions including the rTPJ and left sensorimotor cortex as compared with HCs. In addition, we found that effective longitudinal treatment (acupuncture treatment) can normalize (increase) decreased rs-ReHo in CSNP patients. Moreover, rs-ReHo in the rTPJ was associated with neck pain intensity, as measured by NPQ scores at baseline. Following treatment, rs-ReHo changes at the rTPJ were associated with alleviation of neck pain intensity (i.e., decrease in NPQ scores) in CSNP patients. We applied MVPA for discriminating CSNP patients from HCs with an accuracy of 72% as well as predicting treatment responses using baseline data with an error of 19.6%.
ReHo has been successfully used to detect local abnormalities in subjects with different disorders such as chronic pain (Wang et al., 2014, 2015; Ke et al., 2015; Huang et al., 2016; Zhang et al., 2016) and psychiatric disorders (Liang et al., 2013; Gao et al., 2014; Wang et al., 2015; Yang et al., 2015, 2016). Compared to other resting state functional connectivity (rs-FC) methods such as seed-based rs-FC analysis or independent component analysis, ReHo can target key regions involved in the pathophysiology of neck pain rather than elucidating the connectivity between and among different regions (Yu et al., 2014b).
Neck pain accompanied by neck stiffness and discomfort is the earliest symptom of CS. Patients who do not receive appropriate treatment in the initial stage of CS later develop a complex series of clinical symptoms including spinal cord compression (i.e., CSM). Recently, studies have investigated brain functional changes in CSM (Kowalczyk et al., 2012; Zhou et al., 2014; Zhou F. et al., 2015; Zhou F. Q. et al., 2015) and found patients with CSM are associated with cortical reorganization in the primary motor cortex (Holly et al., 2007; Duggal et al., 2010; Tam et al., 2010; Kowalczyk et al., 2012; Bhagavatula et al., 2016), primary somatosensory cortex (Duggal et al., 2010; Bhagavatula et al., 2016), and sensorimotor cortex (Dong et al., 2008; Zhou et al., 2014; Tan et al., 2015; Zhou F. Q. et al., 2015).
Consistent with previous findings, we found CSNP patients were associated with reduced rs-ReHo in the sensorimotor cortex compared to HCs. This result is also consistent with previous studies that report alterations in sensorimotor cortices in conjunction with clinical deterioration in CSM patients (Dong et al., 2008; Tan et al., 2015). Some studies have found that sensorimotor cortices normalized with the recovery of function following surgical spinal cord decompression in CSM patients (Holly et al., 2007; Dong et al., 2008; Duggal et al., 2010; Tam et al., 2010; Green et al., 2015; Tan et al., 2015; Bhagavatula et al., 2016). In our study, we did not find any significant rs-ReHo changes in the sensorimotor cortex after treatment. We speculate these differences might be due to the different stages of CS investigated in the two studies. Patients with CSM often show motor and sensory deficits, while CSNP patients often show discomfort of the neck including chronic neck pain and stiffness.
In the present study, decreased ReHo in the sensorimotor cortex was mainly localized to S1. Previous studies have indicated that S1 plays a crucial role in nociceptive processing pathways in the human brain (Bushnell et al., 1999; Ploner et al., 1999). In our previous study we also found that chronic low back pain patients have increased cortical thickness and increased volume in S1 and decreased voxel-by-voxel functional connectivity during low pain conditions at S1 compared to HCs (Kong et al., 2013). Our results provide further evidence of rs-ReHo changes in the sensorimotor cortex in CSNP patients, and we speculate this may reflect aberrant central processing of neck pain afferentiation.
We also found reduced rs-ReHo in the rTPJ, including the posterior end of the right superior temporal gyrus and right inferior parietal lobule, in CSNP patients compared to HCs. Our findings were in line with previous studies that reported abnormal structure and function in the TPJ for various chronic pain conditions. It is worth noting that the literature sometimes uses the superior temporal gyrus and inferior parietal lobule instead of the TPJ. For instance, Friebel et al. (2011) used coordinate-based meta-analysis to find that patients with chronic persistent neuropathic pain were associated with activation of the right inferior parietal lobe, a component of the TPJ.
The TPJ is an approximate term that generally refers to the regions including the posterior end of the superior temporal gyrus, the inferior parietal lobule, and the lateral occipital gyrus in the human brain (Mars et al., 2012; Krall et al., 2015). Studies have suggested that the TPJ is an area in the mirror neuron system, which maintains a coherent representation of our own, as well as others’, body (Iacoboni and Dapretto, 2006; Tsakiris et al., 2008; Rizzolatti and Sinigaglia, 2010) and is also involved in imitative motor learning and control (Buccino and Riggio, 2006; Fogassi, 2011).
The mirror neuron system also participates in pain perception (Hojat and Cohen, 2012) and includes the brain regions associated with pain experience, such as the cingulate and insular cortices (Ramachandran and Oberman, 2006). Studies have suggested that the TPJ receives information from the thalamus and limbic system and projects to many other brain areas associated with pain perception and modulation, such as the lateral anterior prefrontal cortex, ventral prefrontal cortex, anterior insula, posterior cingulate, and anterior medial prefrontal cortex (Mars et al., 2012). In our study, we found that effective longitudinal acupuncture treatments can increase reduced rs-ReHo in the rTPJ, and this normalization is associated with a reduction in neck pain intensity in CSNP patients. More studies are needed to explore the role of the TPJ in the development of CSNP.
Machine learning techniques have been widely applied in translational neuroimaging studies to provide a basis for identifying neuropathological features of different diseases and show potential clinical utility beyond current clinical diagnostic categories (Pereira et al., 2009; Arribas et al., 2010; Orru et al., 2012; Zeng et al., 2012; Abi-Dargham and Horga, 2016; Woo et al., 2017). The earliest work identified brain signatures for neurological diseases, particularly Alzheimer’s Disease and related dementia (Kippenhan et al., 1992). The studies have been extended to other neurological and psychiatric conditions, such as Parkinson’s disease (Tang et al., 2010), major depression (Zeng et al., 2012), schizophrenia (Shen et al., 2010), and ADHD (Hart et al., 2014). Most studies have focused on identifying brain signatures for discriminating patients from HCs and consequently established a meaningful neurophysiological basis for the disorder of interest. Another machine learning model uses a patient’s brain measurements to predict treatment outcomes. Most focus on depression and anxiety disorders (Doehrmann et al., 2013; Hahn et al., 2015; Whitfield-Gabrieli et al., 2016) as well as schizophrenia (Sarpal et al., 2016) and Parkinson’s disease (Ye et al., 2016).
To the best of our knowledge, the present study is the first attempt to classify CSNP patients from healthy controls using resting-state fMRI and a state-of-the-art machine learning approach. Therefore, it provides a possible method of using objective measures to diagnose CSNP patients. More importantly, we used baseline resting-state fMRI to predict clinical symptom changes for a longitudinal treatment with an error of 19.6% and prediction-outcome correlation of 0.30. The present approach offers a useful tool for clinicians when determining the effectiveness of acupuncture treatment before executing treatment procedures.
Limitations
It is worth noting that the aim of this study was to investigate the brain’s neural plasticity in CSNP patients, how an effective treatment can modulate altered neural plasticity, and the association between ReHo changes and clinical symptom changes. This paradigm has been used in previous studies from our group (Li et al., 2016) and from other groups (Zhou and Jin, 2008; Rodriguez-Raecke et al., 2009; Sato et al., 2009; Fang et al., 2012). Although our results may provide some novel insights into the mechanisms underlying effective treatment, the lack of sham treatment as a control has significantly limited our ability to infer an underlying longitudinal acupuncture treatment mechanism. Additional controls should be applied in future studies to elucidate the specific effect of the corresponding treatment.
Conclusion
We found that CSNP patients are associated with increased regional coherences at the right temporo-parietal junction and left sensorimotor cortex. Additionally, the alleviation of symptoms after treatment can normalize the reduced rs-ReHo in CSNP patients, and this normalized rs-ReHo value is also associated with a reduction in neck pain intensity. We were also able to use MVPA to discriminate CSNP patients from HCs with an accuracy of 72% as well as predict treatment responses using baseline data with an error of 19.6%. Our results demonstrate the crucial role of the rTPJ in the pathophysiology and development of non-specific neck pain.
Author Contributions
BL and JL proposed the study and were the guarantors. JC, ZW, YT, KJ, JP, CL, and JK analyzed the data and/or wrote the original and revised drafts of the manuscript. All authors contributed to the design and/or interpretation of the study.
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.
Funding. This work was supported by grants from the National Natural Science Foundation of China under Grant Nos. 81673888 and 81603681, the National Key Basic Research and Development Program (973) under Grant No. 2010CB530503, the Science and Technology Project of Guangdong Province under Grant Nos. 2010B080701025, 2014A020212595, and 2016A020215201, the Youth Elite Project of Guangzhou University of Chinese Medicine under Grant No. 20120105, and the Special Fund for the Research of Traditional Chinese Medicine Science and Technology of Guangdong Provincial Hospital of Chinese Medicine under Grant No. NoYN2015QN18.
ReferencesAbi-DarghamA.HorgaG. (2016). The search for imaging biomarkers in psychiatric disorders.221248–1255. 10.1038/nm.419027783066ArbabshiraniM. R.PlisS.SuiJ.CalhounV. D. (2017). Single subject prediction of brain disorders in neuroimaging: promises and pitfalls.145137–165. 10.1016/j.neuroimage.2016.02.07927012503ArribasJ. I.CalhounV. D.AdaliT. (2010). Automatic Bayesian classification of healthy controls, bipolar disorder, and schizophrenia using intrinsic connectivity maps from FMRI data.572850–2860. 10.1109/TBME.2010.208067920876002BaimeM. J. (2016). In chronic nonspecific neck pain, adding Alexander Technique lessons or acupuncture to usual care improved pain.164:JC29. 10.7326/ACPJC-2016-164-6-02926974729BalikiM. N.GehaP. Y.JabakhanjiR.HardenN.SchnitzerT. J.ApkarianA. V. (2008). A preliminary fMRI study of analgesic treatment in chronic back pain and knee osteoarthritis.4:47. 10.1186/1744-8069-4-4718950528BhagavatulaI. D.ShuklaD.SadashivaN.SaligoudarP.PrasadC.BhatD. I. (2016). Functional cortical reorganization in cases of cervical spondylotic myelopathy and changes associated with surgery.40:E2. 10.3171/2016.3.FOCUS163527246485BicikI.RadanovB. P.SchaferN.DvorakJ.BlumB.WeberB. (1998). PET with 18fluorodeoxyglucose and hexamethylpropylene amine oxime SPECT in late whiplash syndrome.51345–350. 10.1212/WNL.51.2.345BinderA. (2007). The diagnosis and treatment of nonspecific neck pain and whiplash.4379–89.BradenB. B.PipeT. B.SmithR.GlaspyT. K.DeatherageB. R.BaxterL. C. (2016). Brain and behavior changes associated with an abbreviated 4-week mindfulness-based stress reduction course in back pain patients.6:e00443. 10.1002/brb3.44326925304BuccinoG.RiggioL. (2006). The role of the mirror neuron system in motor learning.385–15.BushnellM. C.DuncanG. H.HofbauerR. K.HaB.ChenJ. I.CarrierB. (1999). Pain perception: is there a role for primary somatosensory cortex?967705–7709. 10.1073/pnas.96.14.770510393884CekoM.ShirY.OuelletJ. A.WareM. A.StoneL. S.SeminowiczD. A. (2015). Partial recovery of abnormal insula and dorsolateral prefrontal connectivity to cognitive networks in chronic low back pain after treatment.362075–2092. 10.1002/hbm.2275725648842ChangC.-C.LinC.-J. (2011). LIBSVmml: a library for support vector machines.21–27. 10.1145/1961189.1961199Chinese Association of Rehabilitation (2010). Beijing: Chinese Association of Rehabilitation.ChongC. D.GawN.FuY.LiJ.WuT.SchwedtT. J. (2017). Migraine classification using magnetic resonance imaging resting-state functional connectivity data.37828–844. 10.1177/033310241665209127306407DoehrmannO.GhoshS. S.PolliF. E.ReynoldsG. O.HornF.KeshavanA. (2013). Predicting treatment response in social anxiety disorder from functional magnetic resonance imaging.7087–97. 10.1001/2013.jamapsychiatry.522945462DongY.HollyL. T.Albistegui-DuboisR.YanX.MarehbianJ.NewtonJ. M. (2008). Compensatory cerebral adaptations before and evolving changes after surgical decompression in cervical spondylotic myelopathy.9538–551. 10.3171/SPI.2008.10.083119035745DuggalN.RabinD.BarthaR.BarryR. L.GatiJ. S.KowalczykI. (2010). Brain reorganization in patients with spinal cord compression evaluated using fMRI.741048–1054. 10.1212/WNL.0b013e3181d6b0ea20200344FangZ.NingJ.XiongC.ShulinY. (2012). Effects of electroacupuncture at head points on the function of cerebral motor areas in stroke patients: a PET study.2012:902413. 10.1155/2012/90241322956979FogassiL. (2011). The mirror neuron systemml: how cognitive functions emerge from motor organization.7766–75. 10.1016/j.brainresbull.2008.01.01618394524FriebelU.EickhoffS. B.LotzeM. (2011). Coordinate-based meta-analysis of experimentally induced and chronic persistent neuropathic pain.581070–1080. 10.1016/j.neuroimage.2011.07.02221798355GaoW.JiaoQ.LuS.ZhongY.QiR.LuD. (2014). Alterations of regional homogeneity in pediatric bipolar depression: a resting-state fMRI study.14:222. 10.1186/s12888-014-0222-y25095790GreenA.CheongP. W.Fook-ChongS.TiruchelvarayanR.GuoC. M.YueW. M. (2015). Cortical reorganization is associated with surgical decompression of cervical spondylotic myelopathy.2015:389531. 10.1155/2015/38953126609437GuezM. (2006). Chronic neck pain. An epidemiological, psychological and SPECT study with emphasis on whiplash-associated disorders.77 preceding 13–33. 16544560HahnT.KircherT.StraubeB.WittchenH. U.KonradC.StrohleA. (2015). Predicting treatment response to cognitive behavioral therapy in panic disorder with agoraphobia by integrating local neural information.7268–74. 10.1001/jamapsychiatry.2014.174125409415HartH.ChantilukeK.CubilloA. I.SmithA. B.SimmonsA.BrammerM. J. (2014). Pattern classification of response inhibition in ADHD: toward the development of neurobiological markers for ADHD.353083–3094. 10.1002/hbm.2238624123508HaxbyJ. V. (2012). Multivariate pattern analysis of fMRI: the early beginnings.62852–855. 10.1016/j.neuroimage.2012.03.01622425670HojatM.CohenM. J. (2012). “Physicians’ perception of pain as related to empathy, sympathy and the mirror-neuron system,” in , 1st Edn, edsIncayawarM.ToddK. H. (Oxford: Oxford University Press), 63–74.HollyL. T. (2009). Management of cervical spondylotic myelopathy with insights from metabolic imaging of the spinal cord and brain.22575–581. 10.1097/WCO.0b013e3283325ea719741530HollyL. T.DongY.Albistegui-DuBoisR.MarehbianJ.DobkinB. (2007). Cortical reorganization in patients with cervical spondylotic myelopathy.6544–551. 10.3171/spi.2007.6.6.517561743Hotz-BoendermakerS.MarcarV. L.MeierM. L.BoendermakerB.HumphreysB. K. (2016). Reorganization in secondary somatosensory cortex in chronic low back pain patients.41E667–E673. 10.1097/BRS.000000000000134827244113HoyD. G.ProtaniM.DeR.BuchbinderR. (2010). The epidemiology of neck pain.24783–792. 10.1016/j.berh.2011.01.01921665126HrabalekL.HlustikP.HokP.CechakovaE.WanekT.OtrubaP. (2015). [Influence of cervical spondylotic spinal cord compression on cerebral cortical adaptation. radiological study].82404–411. 26787180HrabalekL.HlustikP.HokP.WanekT.OtrubaP.CechakovaE. (2014). [Effects of spinal cord decompression in patients with cervical spondylotic myelopathy oncortical brain activations].93530–535. 25418940HuangT.ZhaoZ.YanC.LuJ.LiX.TangC. (2016). Altered spontaneous activity in patients with persistent somatoform pain disorder revealed by regional homogeneity.11:e0151360. 10.1371/journal.pone.015136026977802IacoboniM.DaprettoM. (2006). The mirror neuron system and the consequences of its dysfunction.7942–951. 10.1038/nrn202417115076KeJ.QiR.LiuC.XuQ.WangF.ZhangL. (2015). Abnormal regional homogeneity in patients with irritable bowel syndrome: a resting-state functional MRI study.271796–1803. 10.1111/nmo.1269226403620KendallM. G. (1990). Oxford: Oxford University Press.KippenhanJ. S.BarkerW. W.PascalS.NagelJ.DuaraR. (1992). Evaluation of a neural-network classifier for PET scans of normal and Alzheimer’s disease subjects.331459–1467. 1634935KongJ.SpaethR. B.WeyH. Y.CheethamA.CookA. H.JensenK. (2013). S1 is associated with chronic low back pain: a functional and structural MRI study.9:43. 10.1186/1744-8069-9-4323965184KowalczykI.DuggalN.BarthaR. (2012). Proton magnetic resonance spectroscopy of the motor cortex in cervical myelopathy.135461–468. 10.1093/brain/awr32822180462KrallS. C.RottschyC.OberwellandE.BzdokD.FoxP. T.EickhoffS. B. (2015). The role of the right temporoparietal junction in attention and social interaction as revealed by ALE meta-analysis.220587–604. 10.1007/s00429-014-0803-z24915964LeakA. M.CooperJ.DyerS.WilliamsK. A.Turner-StokesL.FrankA. O. (1994). The northwick park neck pain questionnaire, devised to measure neck pain and disability.33469–474. 10.1093/rheumatology/33.5.469LiJ.ZhangJ. H.YiT.TangW. J.WangS. W.DongJ. C. (2014). Acupuncture treatment of chronic low back pain reverses an abnormal brain default mode network in correlation with clinical pain relief.32102–108. 10.1136/acupmed-2013-01042324280949LiZ.LiuM.LanL.ZengF.MakrisN.LiangY. (2016). Altered periaqueductal gray resting state functional connectivity in migraine and the modulation effect of treatment.6:20298. 10.1038/srep2029826839078LiangM. J.ZhouQ.YangK. R.YangX. L.FangJ.ChenW. L. (2013). Identify changes of brain regional homogeneity in bipolar disorder and unipolar depression using resting-state FMRI.8:e79999. 10.1371/journal.pone.007999924324588LinnmanC.AppelL.SoderlundA.FransO.EnglerH.FurmarkT. (2009). Chronic whiplash symptoms are related to altered regional cerebral blood flow in the resting state.1365–70. 10.1016/j.ejpain.2008.03.00118486506LoggiaM. L.KimJ.GollubR. L.VangelM. G.KirschI.KongJ. (2013). Default mode network connectivity encodes clinical pain: an arterial spin labeling study.15424–33. 10.1016/j.pain.2012.07.02923111164Lopez-SolaM.WooC. W.PujolJ.DeusJ.HarrisonB. J.MonfortJ. (2017). Towards a neurophysiological signature for fibromyalgia.15834–47. 10.1097/j.pain.000000000000070727583567LorberboymM.GiladR.GorinV.SadehM.LamplY. (2002). Late whiplash syndrome: correlation of brain SPECT with neuropsychological tests and P300 event-related potential.52521–526. 10.1097/00005373-200203000-0001711901329Lorimer MoseleyG. (2005). Widespread brain activity during an abdominal task markedly reduced after pain physiology education: fMRI evaluation of a single patient with chronic low back pain.5149–52. 10.1016/S0004-9514(05)70053-215748125LouwA.PuenteduraE. J.DienerI.PeoplesR. R. (2015). Preoperative therapeutic neuroscience education for lumbar radiculopathy: a single-case fMRI report.31496–508. 10.3109/09593985.2015.103837426395827MaoC. P.ZhangQ. L.BaoF. X.LiaoX.YangX. L.ZhangM. (2014). Decreased activation of cingulo-frontal-parietal cognitive/attention network during an attention-demanding task in patients with chronic low back pain.56903–912. 10.1007/s00234-014-1391-624986218MarsR. B.SalletJ.SchuffelgenU.JbabdiS.ToniI.RushworthM. F. (2012). Connectivity-based subdivisions of the human right “temporoparietal junction area”: evidence for different areas participating in different cortical networks.221894–1903. 10.1093/cercor/bhr26821955921OrruG.Pettersson-YeoW.MarquandA. F.SartoriG.MechelliA. (2012). Using support vector machine to identify imaging biomarkers of neurological and psychiatric disease: a critical review.361140–1152. 10.1016/j.neubiorev.2012.01.00422305994OtteA.EttlinT.FierzL.Mueller-BrandJ. (1996). Parieto-occipital hypoperfusion in late whiplash syndrome: first quantitative SPET study using technetium-99m bicisate (ECD).2372–74. 10.1007/BF017369938586106OtteA.EttlinT. M.NitzscheE. U.WachterK.HoegerleS.SimonG. H. (1997a). PET and SPECT in whiplash syndrome: a new approach to a forgotten brain?63368–372. 9328255OtteA.Mueller-BrandJ.FierzL. (1995). Brain SPECT findings in late whiplash syndrome.345:1513. 10.1016/S0140-6736(95)91075-1OtteA.Mueller-BrandJ.NitzscheE. U.WachterK.EttlinT. M. (1997b). Functional brain imaging in 200 patients after whiplash injury.38:1002. 9189162PereiraF.MitchellT.BotvinickM. (2009). Machine learning classifiers and fMRI: a tutorial overview.45S199–S209. 10.1016/j.neuroimage.2008.11.00719070668PijnenburgM.BrumagneS.CaeyenberghsK.JanssensL.GoossensN.MarinazzoD. (2015). Resting-state functional connectivity of the sensorimotor network in individuals with nonspecific low back pain and the association with the sit-to-stand-to-sit task.5303–311. 10.1089/brain.2014.030925557846PlonerM.SchmitzF.FreundH. J.SchnitzlerA. (1999). Parallel activation of primary and secondary somatosensory cortices in human pain processing.813100–3104. 10.1152/jn.1999.81.6.310010368426RamachandranV. S.ObermanL. M. (2006). Broken mirrors: a theory of autism.29562–69. 10.1038/scientificamerican1106-62RizzolattiG.SinigagliaC. (2010). The functional role of the parieto-frontal mirror circuit: interpretations and misinterpretations.11264–274. 10.1038/nrn280520216547Rodriguez-RaeckeR.NiemeierA.IhleK.RuetherW.MayA. (2009). Brain gray matter decrease in chronic pain is the consequence and not the cause of pain.2913746–13750. 10.1523/JNEUROSCI.3687-09.200919889986SarpalD. K.ArgyelanM.RobinsonD. G.SzeszkoP. R.KarlsgodtK. H.JohnM. (2016). Baseline striatal functional connectivity as a predictor of response to antipsychotic drug treatment.17369–77. 10.1176/appi.ajp.2015.1412157126315980SatoM.InubushiM.ShigaT.HirataK.OkamotoS.KamibayashiT. (2009). Therapeutic effects of acupuncture in patients with rheumatoid arthritis: a prospective study using (18)F-FDG-PET.23311–316. 10.1007/s12149-009-0238-419337783SeminowiczD. A.WidemanT. H.NasoL.Hatami-KhoroushahiZ.FallatahS.WareM. A. (2011). Effective treatment of chronic low back pain in humans reverses abnormal brain anatomy and function.317540–7550. 10.1523/JNEUROSCI.5280-10.201121593339ShenH.WangL.LiuY.HuD. (2010). Discriminative analysis of resting-state functional connectivity patterns of schizophrenia using low dimensional embedding of fMRI.493110–3121. 10.1016/j.neuroimage.2009.11.01119931396ShiY.LiuZ.ZhangS.LiQ.GuoS.YangJ. (2015). Brain network response to acupuncture stimuli in experimental acute low back pain: an fMRI study.2015:210120. 10.1155/2015/21012026161117TamS.BarryR. L.BarthaR.DuggalN. (2010). Changes in functional magnetic resonance imaging cortical activation after decompression of cervical spondylosis: case report.67E863–E864. 10.1227/01.NEU.0000374848.86299.1720657323TanY.ZhouF.WuL.LiuZ.ZengX.GongH. (2015). Alteration of regional homogeneity within the sensorimotor network after spinal cord decompression in cervical spondylotic myelopathy: a resting-state fMRI study.2015:647958. 10.1155/2015/64795826605335TangC. C.PostonK. L.EckertT.FeiginA.FruchtS.GudesblattM. (2010). Differential diagnosis of parkinsonismml: a metabolic imaging study using pattern analysis.9149–158. 10.1016/S1474-4422(10)70002-8TsakirisM.CostantiniM.HaggardP. (2008). The role of the right temporo-parietal junction in maintaining a coherent sense of one’s body.463014–3018. 10.1016/j.neuropsychologia.2008.06.00418601939TuY.TanA.BaiY.HungY. S.ZhangZ. (2016). Decoding subjective intensity of nociceptive pain from pre-stimulus and post-stimulus brain activities.10:32. 10.3389/fncom.2016.0003227148029TuY. H.FuZ. N.TanA.HuangG.HuL.HungY. S. (2017). A novel and effective fMRI decoding approach based on sliced inverse regression and its application to pain prediction.273373–384. 10.1016/j.neucom.2017.07.045UngH.BrownJ. E.JohnsonK. A.YoungerJ.HushJ.MackeyS. (2014). Multivariate classification of structural MRI data detects chronic low back pain.241037–1044. 10.1093/cercor/bhs37823246778VranaA.Hotz-BoendermakerS.StampfliP.HanggiJ.SeifritzE.HumphreysB. K. (2015). Differential neural processing during motor imagery of daily activities in chronic low back pain patients.10:e0142391. 10.1371/journal.pone.014239126569602WandB. M.ParkitnyL.O’ConnellN. E.LuomajokiH.McAuleyJ. H.ThackerM. (2011). Cortical changes in chronic low back pain: current state of the art and implications for clinical practice.1615–20. 10.1016/j.math.2010.06.00820655796WangP.DuH.ChenN.GuoJ.GongQ.ZhangJ. (2014). Regional homogeneity abnormalities in patients with tension-type headache: a resting-state fMRI study.30949–955. 10.1007/s12264-013-1468-625098351WangY.ZhangX.GuanQ.WanL.YiY.LiuC. F. (2015). Altered regional homogeneity of spontaneous brain activity in idiopathic trigeminal neuralgia.112659–2666. 10.2147/NDT.S9487726508861WasanA. D.LoggiaM. L.ChenL. Q.NapadowV.KongJ.GollubR. L. (2011). Neural correlates of chronic low back pain measured by arterial spin labeling.115364–374. 10.1097/ALN.0b013e318220e88021720241Whitfield-GabrieliS.GhoshS. S.Nieto-CastanonA.SayginZ.DoehrmannO.ChaiX. J. (2016). Brain connectomics predict response to treatment in social anxiety disorder.21680–685. 10.1038/mp.2015.10926260493WooC. W.ChangL. J.LindquistM. A.WagerT. D. (2017). Building better biomarkers: brain models in translational neuroimaging.20365–377. 10.1038/nn.447828230847YangH.LiL.PengH.LiuT.YoungA. H.AngstJ. (2016). Alterations in regional homogeneity of resting-state brain activity in patients with major depressive disorder screening positive on the 32-item hypomania checklist (HCL-32).20369–76. 10.1016/j.jad.2016.05.00427280965YangX. Y.SunJ.LuoJ.ZhongZ. X.LiP.YaoS. M. (2015). Regional homogeneity of spontaneous brain activity in adult patients with obsessive-compulsive disorder before and after cognitive behavioural therapy.188243–251. 10.1016/j.jad.2015.07.04826378734YeZ.RaeC. L.NombelaC.HamT.RittmanT.JonesP. S. (2016). Predicting beneficial effects of atomoxetine and citalopram on response inhibition in Parkinson’s disease with clinical and neuroimaging measures.371026–1037. 10.1002/hbm.2308726757216YoungerJ. W.ChuL. F.D’ArcyN. T.TrottK. E.JastrzabL. E.MackeyS. C. (2011). Prescription opioid analgesics rapidly change the human brain.1521803–1810. 10.1016/j.pain.2011.03.02821531077YuR.GollubR. L.SpaethR.NapadowV.WasanA.KongJ. (2014a). Disrupted functional connectivity of the periaqueductal gray in chronic low back pain.6100–108. 10.1016/j.nicl.2014.08.01925379421YuR.GollubR. L.VangelM.KaptchukT.SmollerJ. W.KongJ. (2014b). Placebo analgesia and reward processing: integrating genetics, personality, and intrinsic brain activity.354583–4593. 10.1002/hbm.2249624578196ZangY.JiangT.LuY.HeY.TianL. (2004). Regional homogeneity approach to fMRI data analysis.22394–400. 10.1016/j.neuroimage.2003.12.03015110032ZengL. L.ShenH.LiuL.WangL.LiB.FangP. (2012). Identifying major depression using whole-brain functional connectivity: a multivariate pattern analysis.1351498–1507. 10.1093/brain/aws05922418737ZhangJ.SuJ.WangM.ZhaoY.YaoQ.ZhangQ. (2016). Increased default mode network connectivity and increased regional homogeneity in migraineurs without aura.17:98. 27771875ZhouF.GongH.LiuX.WuL.LukK. D.HuY. (2014). Increased low-frequency oscillation amplitude of sensorimotor cortex associated with the severity of structural impairment in cervical myelopathy.9:e104442. 10.1371/journal.pone.010444225111566ZhouF.WuL.LiuX.GongH.LukK. D.HuY. (2015). Characterizing thalamocortical disturbances in cervical spondylotic myelopathy: revealed by functional connectivity under two slow frequency bands.10:e0125913. 10.1371/journal.pone.012591326053316ZhouF. Q.TanY. M.WuL.ZhuangY.HeL. C.GongH. H. (2015). Intrinsic functional plasticity of the sensory-motor network in patients with cervical spondylotic myelopathy.5:9975. 10.1038/srep0997525897648ZhouY.JinJ. (2008). Effect of acupuncture given at the HT 7, ST 36, ST 40 and KI 3 acupoints on various parts of the brains of Alzheimer’ s disease patients.339–17. 10.3727/03601290880386118618672741