- Open Access
S1 is associated with chronic low back pain: a functional and structural MRI study
Molecular Painvolume 9, Article number: 43 (2013)
A fundamental characteristic of neural circuits is the capacity for plasticity in response to experience. Neural plasticity is associated with the development of chronic pain disorders. In this study, we investigated 1) brain resting state functional connectivity (FC) differences between patients with chronic low back pain (cLBP) and matched healthy controls (HC); 2) FC differences within the cLBP patients as they experienced different levels of endogenous low back pain evoked by exercise maneuvers, and 3) morphometric differences between cLBP patients and matched HC. We found the dynamic character of FC in the primary somatosensory cortex (S1) in cLBP patients, i.e., S1 FC decreased when the patients experienced low intensity LBP as compared with matched healthy controls, and FC at S1 increased when cLBP patients experienced high intensity LBP as compared with the low intensity condition. In addition, we also found increased cortical thickness in the bilateral S1 somatotopically associated with the lower back in cLBP patients as compared to healthy controls. Our results provide evidence of structural plasticity co-localized with areas exhibiting FC changes in S1 in cLBP patients.
Chronic low back pain (cLBP) is one of the most common reasons for all physician visits in the USA and is a leading cause of disability [1, 2]. The etiology of cLBP is heterogeneous ; non-specific cLBP, which is characterized by a lack of recognizable pathology, represents the majority of cLBP patients [3–5]. The financial cost associated with LBP care is staggering; however, treatment for cLBP has achieved limited success . To develop more effective treatments, it is crucial to understand the underlying neurobiology of cLBP.
Recent neuroimaging studies [7–11] have found significant differences in brain functions when comparing patients with cLBP and matched healthy controls (HC). Compared to HC, cLBP patients showed augmented activation in pain related brain regions [9, 12] during administration of comparable experimental pain, differences in activation during emotional decision-making tasks , reduced deactivation in several key default mode network (DMN) brain regions during a simple visual attention task , and increased high-frequency BOLD oscillations circumscribed mainly to the mPFC and parts of the DMN . More recently, using arterial spin labeling, our group has found  that provoked increases in endogenous LBP ratings were positively associated with statistically significant increases in regional cerebral blood flow in a widespread network of cortical areas, including the bilateral medial and dorsolateral prefrontal cortices (mPFC and DLPFC), superior parietal lobules, S1 and S2, and unilaterally in the right insula in cLBP patients. In addition, we  also found that compared with healthy controls, patients demonstrated stronger default mode network connectivity to the pregenual anterior cingulate cortex, left inferior parietal lobule, and right insula. Patients' baseline chronic low back pain intensity was significantly correlated positively with connectivity strength between the DMN and right insula. Taken together, results from previous studies indicate that chronic low back pain is associated with both sensory and overall changes in brain function.
In addition to changes in brain function, multiple morphometric imaging studies report structural changes in the brains of patients with cLBP [16–21]. For instance, using voxel based morphometry (VBM) methods, Apkarian and colleagues  found that cLBP patients showed 5-11% less neocortical gray matter volume than control subjects and that the decreased volume was related to pain duration. These studies indicate that cLBP is also associated with structural pathological changes in the brain.
Few studies have been performed to simultaneously investigate the functional and structural changes in a single study. In one of the few such studies, Seminowicz and colleagues  found that cLBP is associated with decreased cortical thickness and abnormal activity during attention-demanding tasks in the DLPFC. Most intriguingly, after successful treatment, both cortical thickness and functional connectivity (FC) of the DLPFC were normalized. In another study, Baliki and colleagues  found that when low back pain persisted for one year, the gray matter density decreased significantly in brain regions including the striatum and insula, and left sensorimotor cortex. Additionally, they found that increased FC of the nucleus accumbens with the medial prefrontal cortex during the subacute stage of back pain predicted pain persistence. In one study from our group on fibromyalgia (FM) , we found that FM patients displayed a distinct overlap between decreased cortical thickness, brain volumes and measures of functional regional coherence in the rostral anterior cingulate cortex. The morphometric changes were more pronounced with longer exposure to FM pain.
In this study, we investigated 1) brain resting state FC differences between cLBP patients and matched HC; 2) FC differences within the cLBP patients as they experienced different levels of endogenous low back pain (before and after exercise maneuvers), and 3) morphometric differences between cLBP patients and matched HC. There are multiple methods to perform FC; here, we applied a whole-brain voxel-by-voxel hub FC method described in previous studies [25, 26]. The basic analytic strategy was to compute an estimate of the FC of each voxel within the brain by performing Pearson correlation between each voxel and all other voxels of the whole brain . One advantage of this method, we believe, is that rather than focusing on the FC between different brain regions/networks, it focuses on the synchrony of each brain region, which could facilitate the identification of the key regions disrupted in cLBP.
We believe that a combined investigation of both structure and FC changes between cLBP patients and HCs, as well as the dynamic change of FC when cLBP patients were experiencing different endogenous low back pain intensities, will shed new light on our understanding of pathophysiological mechanism, adaptation and reorganization of brain structure and function in chronic pain conditions.
Materials and methods
18 cLBP patients (age = 36.1 ± 9.9, 6 males) and 18 healthy controls, matched for age and gender, completed the study. The Institutional Review Board at Massachusetts General Hospital approved the study and all subjects gave written informed consent.
All patients were clinically diagnosed with nonspecific cLBP with a duration of at least six months by a clinical evaluation, including the use of X-ray/MRI reports, when available. Only those patients meeting Quebec Low Back Pain Task Force classification criteria for Classes I or II were enrolled . For instance, patients reporting radicular pain and/or numbness below the knee were excluded. Subjects were also excluded if they reported major systemic diseases, or history of head injury or coma. cLBP patients were asked to rate their pain using a visual analog scale (0 no pain, 10 maximum imaginable pain). Depressive symptoms were assessed using the Beck Depression Inventory (BDI-II) for all patients who participated in the study . All questionnaires were administered immediately prior to brain scanning. Demographics, clinical assessments and characteristics for cLBP patients and HCs are presented in Table 1. Healthy controls, matched on gender, age and race, were recruited through the flyers/emails in the community. All subjects passed the health subjects screening procedures, indicating that they did not have low back pain.
Medication use per self-report was limited to non-steroidal anti-inflammatory drugs (NSAIDs, e.g., ibuprofen, Motrin, Advil, and Naproxen) and acetaminophen (e.g., Tylenol). Additional non-pharmacological methods of self-reported pain management included chiropractic massages, physical therapy, exercises, and acupuncture.
MRI data acquisition and pain exacerbation procedures
All MRI data were acquired with a 3T Siemens whole-body scanner with echo-planar imaging capability using a 32-channel radio-frequency head coil at the Martinos Center for Biomedical Imaging. During the resting state fMRI scan, subjects were asked to keep their eyes open and look at a darkened screen for 6 minutes. The BOLD fMRI scan acquisition included 47 slices with slice thickness of 3 mm, TR 3000 ms, TE 30 ms, and a 3×3 mm in-plane spatial resolution. T1-weighted MPRAGE type structural images were acquired using the following parameters: voxel size 1.2×1.2×1.2 mm, TR, 2.2 s, TE, 1.54 ms, flip angle: 7 degrees, slices: 144; field of view: 230. Before and immediately after each 6 minute scan, subjects were asked to rate the intensity of their LBP using the 0–10 pain scale.
After the first resting state scan, cLBP patients were taken out of the scanner to perform exercises for a period lasting up to 10 minutes to exacerbate their endogenous lower back pain. These exercises were tailored to each patient based on their report of which movements exacerbated their pain. They included lumbar flexion and extension exercises, sit-ups, or lumbar rotation exercises, where the subject rotated his or her body from side to side at a self-selected speed. The patients were required to perform their exercise(s) slowly, and at maximum flexion/extension, such that they might feel the exercise trigger low back pain. During the screening, all subjects were asked to confirm that they could perform these exercises. If at the end of the first resting state scan the patient’s cLBP pain rating was too strong (≥ 7 in 0–10 scale) and the patient was reluctant to perform exercises to enhance their pain experience, they were asked to wait for 10 minutes in a comfortable position that they chose before starting the second half of the scanning session.
During the exercises, cLBP patients were asked intermittently to report their level of pain using the 0–10 pain scale; the exercises were repeated until subjects reported an increase in pain of approximately 3 points on the pain scale. Once this level of pain was achieved, subjects were placed back in the scanner to repeat the same MRI scans that were acquired before the exercise maneuvers. For healthy controls, structural and resting state scans were only collected once.
Whole brain voxel-wise functional connectivity analysis
To summarize the procedures, data pre-processing included deleting the first four volumes to remove non steady-state images, slice-timing correction, motion correction, and co-registration to the standard brain atlas (MNI152). In addition, the images were corrected for linear trends over each run, low-pass filtered at 0.08 Hz and spatially smoothed with a 4 mm FWHM Gaussian kernel. Nine nuisance parameters (six motion parameters, whole brain, white matter and CSF) were also extracted. Finally, the data were down-sampled into 4 mm isotropic resolution to reduce the computational burden.
Pearson correlation coefficients were calculated between time series data of each voxel and every other voxel of the whole brain. Cross-correlation measures were weighted by the number of strongly correlated links (i.e., correlation above a threshold of r > 0.25) across the entire brain to every given voxel to determine the degree of connectivity, and a z-transformation was applied to prepare for group comparison [25, 26]. The z-transformation did not affect the topography of the connectivity map, but rather normalized each individual’s connectivity map. Since global signal removal might induce artificial negative correlations [29, 30], further analyses were restricted to positive correlations.
Paired and unpaired two-sample t-tests were then used for group analysis to compare FC changes within patients with different endogenous LBP intensities, and between HCs and cLBP patients respectively using FSL. The results were cluster corrected for multiple comparisons and a threshold was set at Z > 2.3 and p < 0.05 with 10 contiguous voxels and spatially smoothed with 1mm spherical kernel for display. Mean Z values over the S1 region (Z > 2.3) were also extracted from individual subjects and plotted for different groups.
To identify networks related to the change in endogenous pain intensity, the difference in subjective pain ratings before and after the maneuver (high pain – low pain) was used as a covariate of interest for regression analysis (Z > 2.3 and p < 0.05 with 10 contiguous voxels) across the whole brain.
Brain structure analysis using FreeSurfer
Automated cortical surface reconstruction was performed on the T1 MPRAGE scans in Free Surfer version 5.1 (http://surfer.nmr.mgh.harvard.edu) using previously described methods [31–33]. Any inaccuracies in the reconstruction of white and pial surfaces of individual subjects were manually corrected before calculating cortical thickness following instructions provided by the software developers. The cortical thickness measure was computed as the distance between the pial and white matter surfaces at each point across the cortical mantle.
Group analyses were performed by resampling each subject’s data to the FreeSurfer average atlas, distributed as a part of FreeSurfer. Cortical thickness maps were smoothed using a Gaussian kernel with a FWHM of 10 mm. Vertex wise analyses of cortical thickness were performed using FreeSurfer. For both cLBP patients and HCs, a linear model of the cortical thickness was calculated at each vertex on the surface. Two-sample t-tests were used to compare the cortical thickness between cLBP patients and HCs. Based on previous structural studies [16, 18, 22, 23], we identified three bilateral a priori regions of interest (ROIs), including the lateral prefrontal cortex (LPFC), primary somatosensory cortex (S1), and insula, and restricted our analysis to a mask of these ROIs. The mask was created by merging the FreeSurfer parcellation labels for the rostral and caudal middle frontal areas, insula and the top portion of the postcentral gyrus. Based on the somatosensory homunculus as defined by Penfield  we included only the portion of the somatosensory cortex that include the cortical representation of the lower back (the top third of S1) in our mask using a method provided by FreeSurfer. For a priori ROIs, a vertex wise threshold was set at p < 0.05 corrected for multiple comparisons in FreeSurfer using Monte Carlo permutations with 5000 iterations within the mask.
In addition, we also conducted a whole brain cortical thickness comparison between cLBP patients and HCs using a vertex wise threshold (p < 0.05) corrected for multiple comparisons in FreeSurfer using Monte Carlo permutations with 5000 iterations. To further explore the relationship between brain cortical thickness changes and duration of cLBP, we also applied a regression analysis to investigate the relationship between cortical thickness and the duration of LBP within cLBP patients.
The automated segmentation procedure for labeling different brain structures and extracting their volumetric measures is described in detail by Fischl et al. . This procedure assigns a neuroanatomical label to each voxel in an MRI volume based on probabilistic information automatically estimated from a manually labeled training set, including both gray and white matter. This technique has been shown to be comparable in accuracy to manual labeling [35, 36]. The automatic segmentations were also visually inspected for accuracy. In the present study, we focus on the volumetric difference in the top one third of primary somatosensory cortex for the following reasons:1) previous studies have suggested the involvement of S1 in both the experimental presentation of noxious stimuli and in pathological pain states such as cLBP [11, 37–39]; in particular, studies have suggested the possibility of a reorganization of S1 in cLBP patients  in the top portion of the postcentral gyrus, and 2) FC results from the present study showed dynamic FC changes at this region (see Results section for more details).
Using the automatic labeling of brain structures described previously, we obtained volumetric measures of the entire S1 cortical region as well as specific portions of S1, divided perpendicular to its long axis. For the present analysis, we divided S1 into 3 sections as determined by the length of the long axis. A two sample t-test was then applied to compare the volume differences between the patients and controls in the bilateral top third of S1. As an exploratory control, we also applied the same analysis on the total S1 volume, as well as bilateral middle and bottom thirds of S1 separately.
A total of 18 cLBP patients, and age- and gender-matched HCs, completed the study.
Please see Table 1 for more details on patient/control characteristics. One patient had strong chronic pain at baseline and thus did not perform any exercises. After lying down for 10 minutes, the patient felt a reduction in low back pain. The patient received the exact same set of scan procedures before and after the 10 minute rest period, comparable to the healthy control condition. This data was included in the data analysis (see Methods for details).
In this study, the LBP intensity was recorded before and after each resting state fMRI scan. We used the average of pre- and post-scan ratings to represent the LBP intensity level for the particular resting scan. All patients, except one who did not perform the exercises, reported an increase in low back pain intensity from pre- to post-maneuver. In our data analysis, we define the scanning period during which patients’ low back pain ratings are lower as the low pain (LP) condition, and the other as the high pain (HP) condition. The average LBP intensity was 3.8 ± 2.5 for LP condition, and 6.7 ± 2.0 for HP condition. No back pain was reported for healthy controls.
Functional connectivity analysis results
Whole brain voxel-by-voxel hub FC in HCs and cLBP patients during the LP condition was compared using a two-sample t-test. Results showed that HCs had significantly greater FC in the right primary somatosensory and motor areas (S1 and M1) located in the upper portion of the pre- and postcentral gyrus compared with patients with cLBP LP condition (Figure 1). The opposite contrast showed that cLBP patients had greater FC in the left fusiformgyrus, occipital gyrus, right posterior cingulate cortex, and inferior parietal gyrus (Table 2). When we compared HCs with cLBP HP condition, HCs showed significant greater FC in the left superior frontal gyrus. No regions showed significant differences in the opposite comparison (Table 2).
A paired t-test comparing cLBP LP and cLBP HP conditions in patients showed that patients in the HP condition had greater FC in bilateral S1 and M1, and left superior frontal cortex, where as patients in the cLBP LP condition showed greater FC in the right inferior parietal lobule, cuneus, and middle occipital gyrus (Table 2 and Figure 1).
To test the association between the FC and subjective cLBP rating change, we applied a regression analysis between the FC difference of cLBP HP versus cLBP LP (HP-LP) and the corresponding low back pain rating differences. The results showed significant positive correlations between FC and LBP rating changes at the left insula, precuneus, amygdala and fusiform (Table 3 and Figure 2).
Volumetric analysis showed that there were no significant differences between cLBP patients and controls with regard to total gray volume (p = .78), subcortical gray volume (p = 0.52) and intracranial volume (p =0.87).
Direct comparison between cLBP patients and the control group showed that the cortical thickness measure of the bilateral postcentral gyri was significantly greater in cLBP patients. No other brain regions showed significant differences between the two groups (Table 4 and Figure 3).
To further explore the relationship between S1 morphometric changes and cLBP, we compared the volume of the top third of the postcentral gyrus bilaterally. The results showed a significant difference in volume between the cLBP patients (mean ± SD, 5416 ± 932 mm3) and matched HCs (4781 ± 783 mm3) (LBP patients > HC) bilaterally in the top third of the postcentral gyrus (p = 0.033). When age and gender were included in the model as covariates, the p value remained significant (p = 0.17).
In a further test of volumetric differences, we compared the whole volume, middle and bottom third of the postcentral gyrus bilaterally. The results showed there were no significant differences between the cLBP patients and controls in total volume (16977 ± 2868 mm3vs 16198 ± 2739 mm3, p = .37), middle third (5696 ± 946 mm3 vs 5554 ± 1080 mm3, p = .411) and bottom third of S1 (5863 ± 1029 mm3 vs 5863 ± 1198 mm3, p = 1.0), which highlights the specificity of the changes in the top third of S1 to cLBP patients.
Voxelwise whole brain regression analysis between the cortical thickness and the duration of LBP within cLBP patients no brain regions survived the threshold we set.
In this study, we investigated the FC and structural difference between cLBP patients and matched HCs. The results of this study identify not only that FC differs between cLBP patients and healthy controls at S1 (the area corresponding to the lower back), but also that FC varies with endogenous low back pain intensity (high intensity pain condition showed stronger FC). In addition, we also found differences in brain structure (cortical thickness and volume) between the cLBP patients and HCs at S1. The combined neuroanatomical and functional changes of S1 imply that this brain region in particular may play an important role in the pathophysiology of cLBP.
Recently, the spontaneous brain activity fluctuations observed in resting state fMRI, have drawn the attention of brain imaging investigators [41–43]. It is believed that low-frequency components of the spontaneous functional MR imaging signal can provide information about the intrinsic functional and anatomical organization of the brain . In this study, we applied a data driven method [25, 26] that calculates the connectivity of each voxel with other voxels of the whole brain. The advantage of this method is that it can elucidate key regions (rather than networks) that show FC changes among different conditions. In this study, we found that the HCs demonstrated significantly greater connectivity bilaterally at S1 than cLBP patients during the LP condition (less severe endogenous low back pain). More interestingly, we found that this decreased FC at S1 during the LP condition increased as endogenous low back pain intensity increased, suggesting the dynamic character of FC at S1. The negative correlation between FC and LBP rating changes at S1 indicates that the rate of increase in FC is more dramatic when pain intensity starts to increase.
The involvement of S1 in pain processing has been known for a long time. Brain imaging studies [11, 37–39, 45–47] have found that S1 is activated during the presentation of noxious stimuli as well as in association with pathological pain states such as chronic low back pain. Although arguable, some investigators believe that S1 plays a prominent and highly modulated role in the sensory aspects of pain, including the localization and discrimination of pain intensity . Thus, the dynamic FC changes observed at S1, in association with the varying endogenous pain intensity levels observed in our study, may represent the important role that S1 plays in monitoring different pain intensities.
In addition to FC changes, we also found increased cortical thickness in the superior section of S1 in the cLBP patients. A further analysis of brain volume in the top one third of the postcentral gyrus (S1), but not the whole S1 volume, also demonstrated that S1 was significantly larger in cLBP patients compared to HCs. This result is consistent with a recent study  in which the authors found increased brain gray matter(GM) density in chronic low back pain patients using a relative large cohort of patients (47 cLBP patients). It is noted that the location of the significant FC changes and cortical thickness differences in S1 is also similar to the recent study .
These results are also consistent with previous animal studies that suggest S1 is highly plastic during both development and in adulthood . This plasticity/reorganization can occur in various situations including peripheral lesions and passive sensory experiences, including chronic pain [21, 40, 49, 50]. In addition, our results are also in line with previous morphometric studies [51–53] that found exercise, learning and stimulation can change brain structures that are associated with the task. For instance, Teutsch and colleagues  found that after 8 consecutive days of identical 20-minute pain sequences, repeated painful stimulation resulted in stimulation dependent substantial increases of gray matter in S1 in healthy subjects (average age 26).
In a recent study of chronic pain (migraine) patients , investigators found that patients had greater cortical thickness at S1 compared to the control group. In another study , investigators found a significant increase in the activation response in areas of the ventral postcentral gyrus (POG) in stroke patients relative to controls. The same ventral POG areas showed a significant increase in cortical thickness in the patients.
Gray matter differences could result from changes in cell volume, synaptic densities, blood flow, or interstitial fluid [45, 56]. In a recent study on healthy subjects (average age 25) , Erpelding and colleagues found strong correlation between greater thermal and pain sensitivity and cortical thickness of the primary somatosensory cortex, indicating that individuals who are highly sensitive to pain have thicker cortical gray matter. Studies [57–60] also suggest that cLBP patients exhibit hypersensitivity in the central nervous system as indicated by lower perception thresholds, lower pain thresholds, lower pain tolerance values, and reduced habituation compared to healthy controls. Thus, increased cortical thickness and increased volume (top one third) in the primary somatosensory cortex in cLBP patients is in line with the central sensitization of chronic pain patients and may represent some level of compensation for the constant experience of pain. One question that remains to be answered is whether the increased changes observed in the primary somatosensory cortex are the consequence of cLBP or a cause in the development of cLBP. Theoretically, individuals with thicker/larger somatosensory cortices may be more sensitive to pain, and thus be more vulnerable to developing cLBP . A study including pre-pain brain structure measurements would be able to answer this question.
Although consistent with a recent study with large cohort of patients , our finding of increased cortical thickness at S1 is different from previous studies [18, 19, 22, 23], in which the authors found a significant decrease in gray matter in the somatosensory cortex (S1). We believe this difference may arise from differences between subgroups of the cLBP population, an observation which has been documented in other chronic pain populations that show differential changes in brain structure (e.g., chronic head pain) . In the previous studies [18, 22, 23] that found somatosensory gray matter decreases, all patients had spinal disk pathology changes; in contrast, subjects in our study were characterized as nonspecific cLBP patients, representing the majority of cLBP patients [4, 5, 61].
In other studies, a discogenic back pain subgroup may have injured the nerve surrounding the affected region and prevented the function/movement of the related regions, which caused the gray matter decrease in the corresponding somatosensory region . In nonspecific cLBP, however, this may not be the case [61, 63], as pain is not attributable to a recognizable specific pathology. For these patients, central vulnerable characteristics that have been associated with high sensitivity to pain , such as increased cortical thickness, may be unique to this subgroup of cLBP patients.
Another potential reason that we see conflicting observations (cortical thickening versus thinning) between studies is the age of study population. In a previous study , Schweinhardt and colleagues found vestibulodynia (PVD) patients had significantly higher gray matter densities in pain modulatory and stress-related areas, including the parahippocampal gyrus/hippocampus and basal ganglia when compared with healthy controls. In a subsequent study from the same group , the authors found that whereas older fibromyalgia patients (average age 56) had the commonly observed gray matter decreases, younger patients (average age 43) showed exclusively gray matter increases relative to age-matched controls in areas including the insula, basal ganglia, and ventrolateral prefrontal cortex, which provides further support of this hypothesis in a different chronic pain population.
In our study, the average age of cLBP patients was 36; the average ages in the previous studies were 50 in the Schmidt-Wilcke et al. study , 46 in the Seminowicz et al. study ; 48  (n.b., the average age of healthy controls in this study is 39, which is lower than the patient population) and 43 in the Baliki et al. studies (the patients in this study have a total duration of only one year of chronic pain, which is much shorter than studies from other groups). Thus, we speculate that the cortical thickness increase observed in S1 in our study may be due to the relatively young age of the patients.
In this study, we also found significant positive correlations between FC changes and LBP rating changes at the left insula and amygdala. The insula is one of the key regions in pain processing [11, 39, 66, 67]. The amygdala is a key region involved in emotion processes such as stress and anxiety [68, 69]. Increased low back pain intensity may be associated with more anxiety and stress. Our result may indicate that the increased LBP intensity can be reflected by brain FC increases in brain regions associated with both pain and emotional processes.
There are some potential limitations in this study. We did not include medication in the model, and although it is unlikely for medication to affect our results, we cannot completely rule out the effect of medication. It is important to note, however, that we excluded all patients using opioids from this study, as a previous study  found that administering oral morphine daily for 1 month can cause anatomical changes in the brain.
Another potential limitation is the order effects between the high endogenous LBP and low endogenous LBP conditions. One challenge of cLBP studies is that once LBP is provoked, it is hard to control without any pharmacological intervention. In this study, we used exercise to provoke the patients’ LBP; thus, the high pain condition tended to follow the low pain condition. However, previous studies have suggested high test-retest reliability across different sessions [44, 71], indicating that the significant findings observed in this study are unlikely due to order effects. Future studies are needed to address the potential confounding order effects in this study.
Finally, we did not collect second scan data for healthy controls. However, previous studies have suggested that resting state FC measurements have moderate to high reliability [72, 73]. In a previous study , we investigated the FC change before and after leg lift maneuver, we found that in cLBP patients, the correlation between these connectivity patterns is significantly associated with the presence of ongoing clinical pain, but not in healthy controls who performed the same maneuver, indicating the specificity of FC change and the clinical pain intensity change. Thus, we believe the difference observed in cLBP patients before and after maneuver tend to be related to the low back pain changes rather than the random changes between the two scan sessions.
In summary, our study showed that compared with healthy controls, cLBP patients have increased cortical thickness and increased volume in the top one third of S1, and decreased FC during low intensity LBP conditions at S1. Interestingly, we also found that the FC at S1 increases when patients experience more endogenous back pain. We believe our results provide the evidence of structural plasticity co-localized with areas exhibiting dynamic FC characteristics in S1 in cLBP patients.
Hart LG, Deyo RA, Cherkin DC: Physician office visits for low back pain. Frequency, clinical evaluation, and treatment patterns from a U.S. National survey. Spine 1995, 20: 11–19. 10.1097/00007632-199501000-00003
Chou R, Shekelle P: Will this patient develop persistent disabling low back pain? Jama 2010, 303: 1295–1302. 10.1001/jama.2010.344
Ehrlich GE: Low back pain. Bull World Health Organ 2003, 81: 671–676.
Savigny P, Watson P, Underwood M, et al.: Low back pain: early management of persistent non-specific low back pain. London: National Collaborating Centre for Primary Care and Royal College of General Practitioners (Clinical guideline 88); 2009.
Chou R, Qaseem A, Snow V, Casey D, Cross JT Jr, Shekelle P, Owens DK: Diagnosis and treatment of low back pain: a joint clinical practice guideline from the American college of physicians and the American pain society. Ann Intern Med 2007, 147: 478–491. 10.7326/0003-4819-147-7-200710020-00006
Bogduk N: Pharmacological alternatives for the alleviation of back pain. Expert Opin Pharmacother 2004, 5: 2091–2098. 10.1517/14656522.214.171.1241
Tagliazucchi E, Balenzuela P, Fraiman D, Chialvo DR: Brain resting state is disrupted in chronic back pain patients. Neurosci Lett 2010, 485: 26–31. 10.1016/j.neulet.2010.08.053
Baliki MN, Baria AT, Apkarian AV: The cortical rhythms of chronic back pain. J Neurosci 2011, 31: 13981–13990. 10.1523/JNEUROSCI.1984-11.2011
Kobayashi Y, Kurata J, Sekiguchi M, Kokubun M, Akaishizawa T, Chiba Y, Konno S-i, Kikuchi S-i: Augmented cerebral activation by lumbar mechanical stimulus in chronic Low back pain patients: an fMRI study. Spine 2009, 34: 2431–2436. 2410.1097/BRS.2430b2013e3181b2431fb2476 10.1097/BRS.0b013e3181b1fb76
Apkarian AV, Baliki MN, Geha PY: Towards a theory of chronic pain. Prog Neurobiol 2009, 87: 81–97. 10.1016/j.pneurobio.2008.09.018
Wasan AD, Loggia ML, Chen LQ, Napadow V, Kong J, Gollub RL: Neural correlates of chronic Low back pain measured by arterial spin labeling. Anesthesiology 2011, 115: 364–374. 10.1097/ALN.0b013e318220e880
Giesecke T, Gracely RH, Grant MA, Nachemson A, Petzke F, Williams DA, Clauw DJ: Evidence of augmented central pain processing in idiopathic chronic low back pain. Arthritis Rheum 2004, 50: 613–623. 10.1002/art.20063
Apkarian AV, Sosa Y, Krauss BR, Thomas PS, Fredrickson BE, Levy RE, Harden RN, Chialvo DR: Chronic pain patients are impaired on an emotional decision-making task. Pain 2004, 108: 129–136. 10.1016/j.pain.2003.12.015
Baliki MN, Chialvo DR, Geha PY, Levy RM, Harden RN, Parrish TB, Apkarian AV: Chronic pain and the emotional brain: specific brain activity associated with spontaneous fluctuations of intensity of chronic back pain. J Neurosci 2006, 26: 12165–12173. 10.1523/JNEUROSCI.3576-06.2006
Loggia ML, Kim J, Gollub RL, Vangel MG, Kirsch I, Kong J, Wasan AD, Napadow V: Default mode network connectivity encodes clinical pain: an arterial spin labeling study. Pain 2013, 154: 24–33. 10.1016/j.pain.2012.07.029
Apkarian AV, Sosa Y, Sonty S, Levy RM, Harden RN, Parrish TB, Gitelman DR: Chronic back pain is associated with decreased prefrontal and thalamic gray matter density. J Neurosci 2004, 24: 10410–10415. 10.1523/JNEUROSCI.2541-04.2004
Flor H: Cortical reorganisation and chronic pain: implications for rehabilitation. J Rehabil Med 2003, 41(Suppl):66–72.
Schmidt-Wilcke T, Leinisch E, Gänssbauer S, Draganski B, Bogdahn U, Altmeppen J, May A: Affective components and intensity of pain correlate with structural differences in gray matter in chronic back pain patients. Pain 2006, 125: 89–97. 10.1016/j.pain.2006.05.004
Baliki MN, Schnitzer TJ, Bauer WR, Apkarian AV: Brain morphological signatures for chronic pain. PLoS One 2011, 6: e26010. 10.1371/journal.pone.0026010
May A: Chronic pain may change the structure of the brain. Pain 2008, 137: 7–15. 10.1016/j.pain.2008.02.034
Henry DE, Chiodo AE, Yang W: Central nervous system reorganization in a variety of chronic pain states: a review. Pm R 2011, 3: 1116–1125. 10.1016/j.pmrj.2011.05.018
Seminowicz DA, Wideman TH, Naso L, Hatami-Khoroushahi Z, Fallatah S, Ware MA, Jarzem P, Bushnell MC, Shir Y, Ouellet JA, Stone LS: Effective treatment of chronic low back pain in humans reverses abnormal brain anatomy and function. J Neurosci 2011, 31: 7540–7550. 10.1523/JNEUROSCI.5280-10.2011
Baliki MN, Petre B, Torbey S, Herrmann KM, Huang L, Schnitzer TJ, Fields HL, Apkarian AV: Corticostriatal functional connectivity predicts transition to chronic back pain. Nat Neurosci 2012. 10.1038/nn.3153
Jensen K, Srinivasan P, Spaeth R, Tan Y, Kosek E, Petzke F, Carville S, Fransson P, Marcus H, Williams SC, et al.: Overlapping structural and functional brain changes in patients with long-term exposure to fibromyalgia. Arthritis & Rheumatism 2013. In press
Buckner RL, Sepulcre J, Talukdar T, Krienen FM, Liu H, Hedden T, Andrews-Hanna JR, Sperling RA, Johnson KA: Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease. J Neurosci 2009, 29: 1860–1873. 10.1523/JNEUROSCI.5062-08.2009
Sepulcre J, Liu H, Talukdar T, Martincorena I, Yeo BT, Buckner RL: The organization of local and distant functional connectivity in the human brain. PLoS Comput Biol 2010, 6: e1000808. 10.1371/journal.pcbi.1000808
Werneke MW, Hart DL: Categorizing patients with occupational low back pain by use of the Quebec task force classification system versus pain pattern classification procedures: discriminant and predictive validity. Phys Ther 2004, 84: 243–254.
Beck AT, Ward CH, Mendelson M, Mock J, Erbaugh J: An inventory for measuring depression. Arch Gen Psychiatry 1961, 4: 561–571. 10.1001/archpsyc.1961.01710120031004
Fox MD, Snyder AZ, Vincent JL, Corbetta M, Van Essen DC, Raichle ME: The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci USA 2005, 102: 9673–9678. 10.1073/pnas.0504136102
Saad ZS, Gotts SJ, Murphy K, Chen G, Jo HJ, Martin A, Cox RW: Trouble at rest: how correlation patterns and group differences become distorted after global signal regression. Brain Connect 2012, 2: 25–32. 10.1089/brain.2012.0080
Dale AM, Fischl B, Sereno MI: Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 1999, 9: 179–194. 10.1006/nimg.1998.0395
Fischl B, Salat DH, Busa E, Albert M, Dieterich M, Haselgrove C, van der Kouwe A, Killiany R, Kennedy D, Klaveness S, et al.: Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron 2002, 33: 341–355. 10.1016/S0896-6273(02)00569-X
Fischl B, Sereno MI, Dale AM: Cortical surface-based analysis. II: inflation, flattening, and a surface-based coordinate system. Neuroimage 1999, 9: 195–207. 10.1006/nimg.1998.0396
Penfield W, Rasmussen T: The cerebral cortex of Man. New York: Macmillan; 1955.
Salat DH, Buckner RL, Snyder AZ, Greve DN, Desikan RS, Busa E, Morris JC, Dale AM, Fischl B: Thinning of the cerebral cortex in aging. Cereb Cortex 2004, 14: 721–730. 10.1093/cercor/bhh032
Kuperberg GR, Broome MR, McGuire PK, David AS, Eddy M, Ozawa F, Goff D, West WC, Williams SC, van der Kouwe AJ, et al.: Regionally localized thinning of the cerebral cortex in schizophrenia. Arch Gen Psychiatry 2003, 60: 878–888. 10.1001/archpsyc.60.9.878
Bushnell MC, Duncan GH, Hofbauer RK, Ha B, Chen JI, Carrier B: Pain perception: is there a role for primary somatosensory cortex? Proc Natl Acad Sci USA 1999, 96: 7705–7709. 10.1073/pnas.96.14.7705
Kong J, White NS, Kwong KK, Vangel MG, Rosman IS, Gracely RH, Gollub RL: Using fMRI to dissociate sensory encoding from cognitive evaluation of heat pain intensity. Hum Brain Mapp 2006, 27: 715–721. 10.1002/hbm.20213
Kong J, Loggia ML, Zyloney C, Tu P, Laviolette P, Gollub RL: Exploring the brain in pain: activations, deactivations and their relation. Pain 2010, 148: 257–267. PMID: 20005043 10.1016/j.pain.2009.11.008
Flor H, Braun C, Elbert T, Birbaumer N: Extensive reorganization of primary somatosensory cortex in chronic back pain patients. Neurosci Lett 1997, 224: 5–8. 10.1016/S0304-3940(97)13441-3
Biswal B, Yetkin FZ, Haughton VM, Hyde JS: Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med 1995, 34: 537–541. 10.1002/mrm.1910340409
Raichle ME, Mintun MA: Brain work and brain imaging. Annu Rev Neurosci 2006, 29: 449–476. 10.1146/annurev.neuro.29.051605.112819
Fox MD, Raichle ME: Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci 2007, 8: 700–711. 10.1038/nrn2201
Van Dijk KR, Hedden T, Venkataraman A, Evans KC, Lazar SW, Buckner RL: Intrinsic functional connectivity as a tool for human connectomics: theory, properties, and optimization. J Neurophysiol 2010, 103: 297–321. 10.1152/jn.00783.2009
Erpelding N, Moayedi M, Davis KD: Cortical thickness correlates of pain and temperature sensitivity. Pain 2012, 153: 1602–1609. 10.1016/j.pain.2012.03.012
Teutsch S, Herken W, Bingel U, Schoell E, May A: Changes in brain gray matter due to repetitive painful stimulation. Neuroimage 2008, 42: 845–849. 10.1016/j.neuroimage.2008.05.044
Davis KD, Moayedi M: Central Mechanisms of Pain Revealed Through Functional and Structural MRI. J Neuroimmune Pharmacol 2013, 8(3):518–534. 10.1007/s11481-012-9386-8
Ung H, Brown JE, Johnson KA, Younger J, Hush J, Mackey S: Multivariate classification of structural MRI data detects chronic Low back pain. Cereb Cortex 2012. 10.1093/cercor/bhs1378
Feldman DE, Brecht M: Map plasticity in somatosensory cortex. Science 2005, 310: 810–815. 10.1126/science.1115807
Eto K, Wake H, Watanabe M, Ishibashi H, Noda M, Yanagawa Y, Nabekura J: Inter-regional contribution of enhanced activity of the primary somatosensory cortex to the anterior cingulate cortex accelerates chronic pain behavior. J Neurosci 2011, 31: 7631–7636. 10.1523/JNEUROSCI.0946-11.2011
Gaser C, Schlaug G: Brain structures differ between musicians and non-musicians. J Neurosci 2003, 23: 9240–9245.
Draganski B, Gaser C, Busch V, Schuierer G, Bogdahn U, May A: Neuroplasticity: changes in grey matter induced by training. Nature 2004, 427: 311–312. 10.1038/427311a
May A, Hajak G, Ganssbauer S, Steffens T, Langguth B, Kleinjung T, Eichhammer P: Structural brain alterations following 5 days of intervention: dynamic aspects of neuroplasticity. Cereb Cortex 2007, 17: 205–210.
DaSilva AF, Granziera C, Snyder J, Hadjikhani N: Thickening in the somatosensory cortex of patients with migraine. Neurology 2007, 69: 1990–1995. 10.1212/01.wnl.0000291618.32247.2d
Schaechter JD, Moore CI, Connell BD, Rosen BR, Dijkhuizen RM: Structural and functional plasticity in the somatosensory cortex of chronic stroke patients. Brain 2006, 129: 2722–2733. 10.1093/brain/awl214
Gage FH: Neurogenesis in the adult brain. J Neurosci 2002, 22: 612–613.
Puta C, Schulz B, Schoeler S, Magerl W, Gabriel B, Gabriel HH, Miltner WH, Weiss T: Enhanced sensitivity to punctate painful stimuli in female patients with chronic low back pain. BMC Neurol 2012, 12: 98. 10.1186/1471-2377-12-98
Neziri AY, Curatolo M, Limacher A, Nuesch E, Radanov B, Andersen OK, Arendt-Nielsen L, Juni P: Ranking of parameters of pain hypersensitivity according to their discriminative ability in chronic low back pain. Pain 2012, 153: 2083–2091. 10.1016/j.pain.2012.06.025
Kleinbohl D, Holzl R, Moltner A, Rommel C, Weber C, Osswald PM: Psychophysical measures of sensitization to tonic heat discriminate chronic pain patients. Pain 1999, 81: 35–43. 10.1016/S0304-3959(98)00266-8
Flor H, Diers M, Birbaumer N: Peripheral and electrocortical responses to painful and non-painful stimulation in chronic pain patients, tension headache patients and healthy controls. Neurosci Lett 2004, 361: 147–150. 10.1016/j.neulet.2003.12.064
Balague F, Mannion AF, Pellise F, Cedraschi C: Non-specific low back pain. Lancet 2012, 379: 482–491. 10.1016/S0140-6736(11)60610-7
Taylor KS, Anastakis DJ, Davis KD: Cutting your nerve changes your brain. Brain 2009, 132: 3122–3133. 10.1093/brain/awp231
Staud R: Evidence for shared pain mechanisms in osteoarthritis, low back pain, and fibromyalgia. Curr Rheumatol Rep 2011, 13: 513–520. 10.1007/s11926-011-0206-6
Schweinhardt P, Kalk N, Wartolowska K, Chessell I, Wordsworth P, Tracey I: Investigation into the neural correlates of emotional augmentation of clinical pain. Neuroimage 2008, 40: 759–766. 10.1016/j.neuroimage.2007.12.016
Ceko M, Fitzcharles MA, Bushnell MC, Schweinhardt P: Cerebral grey and white matter changes in fibromyalgia depend on patients’ age. In 14th World Congress on Pain. Milan, Italy; 2012.
Craig AD: Pain mechanisms: labeled lines versus convergence in central processing. Annu Rev Neurosci 2003, 26: 1–30. 10.1146/annurev.neuro.26.041002.131022
Tracey I, Mantyh PW: The cerebral signature for pain perception and its modulation. Neuron 2007, 55: 377–391. 10.1016/j.neuron.2007.07.012
Davidson RJ, McEwen BS: Social influences on neuroplasticity: stress and interventions to promote well-being. Nat Neurosci 2012, 15: 689–695. 10.1038/nn.3093
Kim MJ, Loucks RA, Palmer AL, Brown AC, Solomon KM, Marchante AN, Whalen PJ: The structural and functional connectivity of the amygdala: from normal emotion to pathological anxiety. Behav Brain Res 2011, 223: 403–410. 10.1016/j.bbr.2011.04.025
Younger JW, Chu LF, D'Arcy NT, Trott KE, Jastrzab LE, Mackey SC: Prescription opioid analgesics rapidly change the human brain. Pain 2011, 152: 1803–1810. 10.1016/j.pain.2011.03.028
Faria AV, Joel SE, Zhang Y, Oishi K, van Zjil PC, Miller MI, Pekar JJ, Mori S: Atlas-based analysis of resting-state functional connectivity: evaluation for reproducibility and multi-modal anatomy-function correlation studies. Neuroimage 2012, 61: 613–621. 10.1016/j.neuroimage.2012.03.078
Zuo XN, Kelly C, Adelstein JS, Klein DF, Castellanos FX, Milham MP: Reliable intrinsic connectivity networks: test-retest evaluation using ICA and dual regression approach. Neuroimage 2010, 49: 2163–2177. 10.1016/j.neuroimage.2009.10.080
Shehzad Z, Kelly AM, Reiss PT, Gee DG, Gotimer K, Uddin LQ, Lee SH, Margulies DS, Roy AK, Biswal BB, et al.: The resting brain: unconstrained yet reliable. Cereb Cortex 2009, 19: 2209–2229. 10.1093/cercor/bhn256
This work was supported by KO1AT003883 (NCCAM), R21AT004497 (NCCAM), R03AT218317 (NIDA), R01AT006364 (NCCAM) to Jian Kong, R01AT005280 (NCCAM) to Randy Gollub, K23DA020681 (NIDA) and Arthritis Investigator Award to Ajay Wasan, PO1-AT002048 (NCCAM) to Bruce Rosen, M01-RR-01066 and UL1 RR025758-01 for Clinical Research Center Biomedical Imaging Core from National Center for Research Resources (NCRR), and P41RR14075 for Center for Functional Neuroimaging Technologies from NCRR.
There is no conflict of interest to claim for all authors.
JK: experiment design, data collection, data analysis, manuscript preparation. RS: data collection, data analysis, manuscript preparation. HYW: data analysis, manuscript preparation. AC: data collection, data analysis, manuscript preparation. AC: data analysis, manuscript preparation. KJ: manuscript preparation. YT: data analysis. HSL: data analysis. DHW: data analysis. ML: manuscript preparation. VN: experiment design. JS: experiment design, manuscript preparation. AW: experiment design, manuscript preparation. RG: experiment design, manuscript preparation. All authors read and approved the final manuscript.