Abstract
Background
Chronic pain conditions are highly prevalent in patients with mild traumatic brain injury. Supraspinal diffuse axonal injury is known to dissociate brain functional connectivity in these patients. The effect of this dissociated state on supraspinal pain network is largely unknown. A functional magnetic resonance imaging study was conducted to compare the supraspinal pain network in patients with mild traumatic brain injury to the gender and age-matched healthy controls with the hypothesis that the functional connectivities of the medial prefrontal cortices, a supraspinal pain modulatory region to other pain-related sensory discriminatory and affective regions in the mild traumatic brain injury subjects are significantly reduced in comparison to healthy controls.
Results
The mild traumatic brain injury group (N = 15) demonstrated significantly (P < 0.01, cluster threshold > 150 voxels) less activities in the thalamus, pons, anterior cingulate cortex, insula, dorsolateral prefrontal cortex, and medial prefrontal cortices than the healthy control group (N = 15). Granger Causality Analyses (GCA) indicated while the left medial prefrontal cortices of the healthy control group cast a noticeable degree of outward (to affect) causality inference to multiple pain processing related regions, this outward inference pattern was not observed in the mild traumatic brain injury group. On the other hand, only patients’ bilateral anterior cingulate cortex received multiple inward (to be affected) causality inferences from regions including the primary and secondary somatosensory cortices and the inferior parietal lobe. Resting state functional connectivity analyses indicated that the medial prefrontal cortices of the mild traumatic brain injury group demonstrated a significantly (P < 0.01, F = 3.6, cluster size > 150 voxels) higher degree of functional connectivity to the inferior parietal lobe, premotor and secondary somatosensory cortex than the controls. Conversely, the anterior cingulate cortex of the healthy group demonstrated significantly (P < 0.01, F = 3.84, cluster size > 150 voxels) less degree of functional connectivities to the inferior parietal lobe and secondary somatosensory cortex than their mild traumatic brain injury counterparts.
Conclusions
In short, the current study demonstrates that patients with mild traumatic brain injury and headaches appear to have an altered state of supraspinal modulatory and affective functions related to pain perception.
Keywords
Introduction
Chronic pain conditions such as persistent headache are highly prevalent in patients with mild traumatic brain injury (MTBI). This patient population was known to have a state of microscopic diffuse central axonal injury affecting supraspinal functional connectivities.1,2 To what extent this impaired functional connectivity may affect supraspinal pain processing both at evoked and resting states is largely unknown. As in other neurological diseases, understanding the underlying neurological functional changes can facilitate treatment development. Based on previous studies, the supraspinal pain processing network is known to involve: (1) the thalamus (TH) and pons, which relate sensory afferent signals to other supraspinal regions; (2) the sensory discriminatory regions including the primary and secondary somatosensory cortices (SSC1 and SSC2), and the inferior parietal lobe (IPL); (3) the affective regions such as the anterior cingulate cortex (ACC) and insula (IN); and(4) the modulatory regions involving the dorsolateral prefrontal cortex (DLPFC), and various regions of the prefrontal cortices (PFCs).3,4 The IN is also implicated in assessing the magnitude of pain.3,5,6 Furthermore, the IPL is also known to be involved in spatial discriminatory functions of pain perception.7–9
Chronic pain conditions can occur as a result of maladaptation in the supraspinal pain processing and functional connectivity.3,10 In the case of MTBI, a recent study with cranial pressure pain threshold assessments suggested that the occurrence of chronic pain in MTBI could be attributed to an elevated supraspinal affective pain state and/or a lack of supraspinal modulatory functions in this patient population. 11 However, further confirmatory studies are required to support this assertion and provide guidance for treatment development.
Here, the authors hypothesize that patients with MTBI suffer from a state of altered supraspinal modulatory and affective response to pain. To assess this hypothesis, a study with functional magnetic resonance imaging (fMRI) was conducted to compare the supraspinal resting state functional connectivity and response to evoked heat pain (HP) in patients with MTBI-related headache with gender and age-matched healthy controls.
Methods
With institutional human subject committee approval, subjects (all Veterans) who attended the traumatic brain injury (TBI) clinic were consented, screened, and enrolled based on the following inclusion criteria: male or female age between 18 and 60; history of MTBI and established diagnosis of posttraumatic headache based on the ICHD-212,13 diagnostic criteria including:
Headache, no typical characteristics known, fulfilling criteria C and D Head trauma that includes the following:
either no loss of consciousness or loss of consciousness of <30 min duration Glasgow Coma Scale (GCS) ≥13 symptoms and/or signs diagnostic of concussion as discussed in the below diagnostic criteria of MTBI Headache occurs within seven days after head trauma Headache persists for >three months after head trauma
Cohort of patients with mild traumatic brain injury related headache (MTBI-HA).
HA: headache; M-VAS: Mechanical Visual Analog Scale Score; M: male; F: female; B: blast related; NB: non-blast related; FALL: falling accident; BHT: blunt head trauma; MCA: motorcycle accident; IED: improvised explosive device.
Daily headache medication.
As needed headache medication.
Prescanning assessments
Prior to each fMRI session, each subject was asked to rate the average intensity of their persistent daily headaches on a M-VAS. 14 An elastic band consisting of 13 increments was used to consistently mark the location of the thermal threshold measurement and HP stimulation. The band was extended from the medial malleolus to the medial tibial plateau and the testing/stimulation location was marked at the medial aspect of the left calf between the sixth and seventh marking of an elastic band. A Thermal Sensory Analyzer (Medoc Advanced Medical Systems, MN) was used to assess the non-noxious and noxious thermal thresholds including cool and warm, cold, and heat pain thresholds via a thermode with a contact area of 30 × 30 mm. This method of peripheral sensory testing has been well established in literature and commonly applied to pain-related studies.18–22 After the initial thermal sensory threshold assessments, a 6-sec HP stimulus at a subject specific threshold was applied at the marked location. The subjects were then asked to rate the intensities of HP on a M-VAS. During the HP-fMRI session, the HP threshold temperature with the intensity score ≥30 was delivered as the HP stimulus to the subject. This method of subject specific HP delivery minimizes any unnecessary distraction to the subjects during the scanning and thus optimizing the correlation of supraspinal responses with the stimulation paradigms.16,17,23
Neuroimaging data acquisition
Patients were asked not to take any headache-related medications on the day of scanning. Head movement during scanning were minimized by instructing the subjects to hold the head still during the scanning, applying padding between the subjects’ head and the head coil, and having subjects wear a cervical collar to minimize both lateral and axial head movements.23,24 A 5-min resting fMRI data were collected via a 1.5T GE scanner with T2*-weighted EPI-sequence (TE = 30 ms, TR = 2.0 s, α = 90°, TH = 4 mm, 32 slices, FOV = 220 × 220 mm2, MA = 64 × 64). Evoked HP fMRI was conducted with the same pulse sequence and intermittent 6 sec of subject threshold specific HP stimuli delivered via a fMRI-compatible Peltier probe at the medial aspect of the left mid-calf with various intervals (20–40 sec) of a baseline temperature at 32℃. This established stimulation paradigm was repeated 20 times to complete the session. 16 Anatomical scans were obtained with rapid gradient-echo (MP RAGE) samplings (176 slices T1, 256 × 256 and 1 cm slice thickness).
Data analysis
A simple paired t-test analysis was conducted for the thermal sensory threshold data, and within- and between-group random effect analyses, and GCA were conducted for the fMRI data in the Brain Voyager platform using protocols described by Goebel et al. 25 Using Self-organizing group level Independent Component Analysis (SogICA),26,27 between-group (MTBI > healthy) pain-related clusters were adopted as seeded regions for subsequent ANCOVA resting state functional connectivity analyses.
fMRI data processing and analysis
Preprocessing of functional data
First, raw functional data in the Dicom format was converted into Brain Voyager’s internal “FMR” data format. Standard preprocessing steps including head motion correction, slice scan time correction, drift removal, and spatial smoothing with Gaussian filter (FWHM = 5 mm) were performed for all functional data sets.
Preprocessing of anatomical data
The anatomical data of each subject was also converted into Brain Voyager’s “VMR” data format. Intensity inhomogeneities correction was performed, and the data were then transformed into AC-PC Talairach standard space after being resampled to 1-mm resolution. To avoid quality loss due to successive data sampling, the three spatial transformations were combined and applied backward in one step. In addition, to form a single 4 × 4 transformation matrix, the two affine transformations, iso-voxel scaling, and AC–PC transformation were concatenated. Followed by application of the inverted spatial transformation matrix, the coordinates of each voxel in the target (Talairach) space was affined back to the original space so that data points could be sampled with sinc interpolation in the original three dimensional (3D) space.
Brain segmentation
The brain was further segmented from surrounding head tissue using an automatic “brain peeling” tool as described by Goebel et al. 25 so that a 3-D image of the brain could be visualized.
Cortex segmentation
For cortical functional data analysis, additional brain segmentation was carried out for the gray/white matter boundary using automatic segmentation routines. 28 Applying an analysis of intensity histograms, the white/gray matter border was segmented via a region-growing method. Additional smoothing steps were subsequently applied using a “bridge removal” algorithm in which each segmented hemisphere was smoothed to topologically correct mesh representations. 28 Repeated small morphing steps were also conducted until the central sulcus were visible and areas of activities including those in the sulcus could be visualized. The inflated cortical meshes were used to sample the functional data so that subsequent cortex-based data analysis in a mesh time course format could be carried out at each vertex (node) for each run of each subject.
Normalization of functional data
To transform the functional data into Talairach space, the functional time series data on each subject were first coregistered with the corresponding 3D anatomical data set. The coregistered data were then transformed into the Talairach space resulting in normalized 4D volume time course data ready for both within- and between-group analyses.
Within- and between-group general linear model analysis
First, a protocol file representing the onset and duration of the events for the stimulation conditions was derived. To account for the hemodynamic delay and dispersion, each of the predictors was derived by convoluting an appropriate box-car waveform with a double-gamma hemodynamic response function 29 to extract brain regions with both positively (activation) and negatively (deactivation) correlated blood oxygen level dependent responses. In addition, to address the issues of multiple comparisons, a False Discovery Rate correction was applied. For both within- and between-group random effect analyses, any supraspinal regions with significant (P < 0.01 and cluster threshold > 150 voxels) activation and deactivation were recorded. These cluster size and statistical significance cutoffs were used in previous studies using similar HP stimulation paradigms.16,17,23 The spatial coordinates of significant clusters (VOI) were first complied in a text file allowing anatomical naming via Talairach Client (http://www.talairach.org/client.html). The naming of the regions was further confirmed with BV Tutor probabilistic anatomical map. For 3D cortical labeling, Patches of Interest (POI) were first selected using BVQX “POI Analysis Tool.” The POI’s were then converted to VOI’s to extrapolate the details (including spatial coordinates) of the defined region and classify them using the above-mentioned VOI naming methodology.
Granger Causality Analysis
GCA was conducted to explore the causal interaction (inference) of seeded region(s) of interest (SROI) extracted from between-group random effect analysis according to an established analytical procedure.16,17 A SROI was used to estimate effective connectivity among clusters in the group with the Brain Voyager’s GCA plug-in. The result of the analysis was noted as either positive values signifying significant outward (to affect) inference from the SROI to the targeted regions or negative values representing inward (to be affected) inference to the SROI. 30 In addition, clusters information including coordinates, sizes, and Brodmann areas were converted by the Talairach Client into a text format after verifying the data with Brain Tutor.31,32 The resulting text was imported to a spreadsheet, and the network of inference was mapped onto a spatial representation of the brain network involved in acute thermal pain processing within each group.
All applied fMRI data processing steps, analytical approaches, and selected cluster threshold were similar to previous published HP-related studies.16,17,23
Resting state functional connectivity
Using Self-organizing group level Independent Component Analysis (SogICA),26,27 significant between-group (MTBI > healthy) pain-related clusters in the medial prefrontal cortices (MPFCs) and ACC were adopted as seeded regions for subsequent ANCOVA resting state functional connectivity analyses.
Results
Patient and healthy subjects demographic data
A total of 15 patients with MTBI and 15 healthy controls were enrolled in the study. The patient cohort consisted of 11 male subjects with an average age (years old ± SD) of 34.6 ± 8.7 and 4 female subjects with an average age of 31.2 ± 4.0. The average ages of the male and female healthy subjects were 33.6 ± 9.2 and 29.5 ± 5.5, respectively. Although three patients (Subject #s 1, 4, and 5 in Table 1) reported an additional head trauma event prior to the headache causing MTBI incidence, the remainder recalled only a single MTBI event leading to the symptom of persistent headache after the injury. In the healthy control group, none of the subjects reported any significant head trauma meeting the diagnostic criteria of either MTBI or posttraumatic chronic headache as stipulated in the study enrollment criteria. In the patient group, the average duration of headache was around 70 months, and the average daily headache intensity on a 0–100 Mechanical Analogue Scale (M-VAS ± SD) was 57.7 ± 16.4 (Table 1). Nine MTBI patients reported blast-related injury due to improvised explosive devices and the remaining six subjects sustained non-blast-related head trauma. The average headache intensities (M-VAS ± SD) for blast-related and non-blast-related injury were 62.2 ± 15.8 and 50.8 ± 16.2, respectively. No statistical significant difference (P = 0.20) in the headache intensity was found between the two injury mechanisms in a two-tail t-test analysis. All anatomical brain MRI showed no gross pathology.
Thermal sensory thresholds analysis
The cool, warm, cold, and heat pain thresholds (℃ ± SD) for the healthy subjects were 24.3 ± 3.0, 39.2 ± 3.9, 3.0 ± 4.5, and 47.6 ± 2.2, respectively, and for the patients were 23.2 ± 2.7, 37.8 ± 2.7, 1.5 ± 1.9, and 46.1 ± 3.0, respectively. No significant difference in thermal sensory thresholds was found between the two groups.
Within-group fMRI data analysis
MTBI within-group analysis.
ACC: Anterior Cingulate Cortex.
Healthy within-group analysis.
SSC: Secondary Somatosensory Cortex.
Between-group fMRI data analysis
The between-group random effect analysis (MTBI group > Healthy group) indicated the MTBI group demonstrated significantly (P < 0.01, cluster threshold > 150 voxels) less activity in the thalamus, pons, insula, DLPFC, MPFCs, and ACC (Figure 1(a)) than the healthy subject group when both groups were exposed to a similar intensity and duration of HP stimuli (see Table 4 and Figure 1(b)).
(a) Between-group brain activity differences in response to heat pain stimulation at the seeded regions (medial prefrontal cortex and anterior cingulate cortex) marked with crosshair showing significantly (P < 0.01, cluster threshold > 150 voxels) less activities (blue color in the reference Z-scale) in the mild traumatic brain injury patients in comparison to their healthy counterparts. A: Anterior: P: Posterior; R: Right; L: Left. (b) An overall cortical projection of between-group (mild traumatic brain injury patients minus healthy controls) differences in response to heat pain stimulation with brain areas showing significantly (P < 0.01, cluster threshold > 150 voxels) less activities (blue color in the Z-scale reference). PFCs: Medial Prefrontal Cortices; IPL: inferior Parietal Lobe; ACC: Anterior Cingulate Cortex; TH: Thalamus; IN: Insula. Between-group analysis (MTBI > healthy). ACC(B): bilateral anterior cingulate cortices; SSC2: secondary somatosensory cortex; DLPFC: dorsolateral prefrontal cortex; PFCs: prefrontal cortices; BA: brodmann area; n/a: non-applicable.
Granger Causality Analysis

(a and b). Granger Causality Analysis of left Medial Prefrontal Cortices (PFCs) and right Anterior Cingulate Cortex (ACC) as seeded region of interest (SROI) shown in (a) and (b), respectively. Healthy controls demonstrated significant outward (to affect) causality inference from the PFCs and no inward (to be affected) causality inference to ACC, whereas, patients with MTBI demonstrated significant inward inference to ACC and minimal outward inference from the left PFCs. The direction (>>>>) of the arrows indicates the causality relationship between the SROI and other brain regions. Green arrow indicates outward inference originating from SROI at either left (LT) or right (RT) hemisphere to the left (LT) hemisphere, or inference from the left (LT) hemisphere to the SROI at either left (LT) or right (RT) hemisphere. Orange arrow indicates outward inference originating from SROI at either left (LT) or right (RT) hemisphere to the right (LT) hemisphere, or inference from the right (RT) hemisphere to the SROI at either left (LT) or right (RT) hemisphere. SSC1: primary somatosensory cortex; SSC2: secondary somatosensory cortex; IPL: inferior parietal lobe; TH: thalamus; IN: insula; PM: premotor cortex; M: motor cortex; DLPFC: dorsolateral prefrontal cortex; AMG: amygdala; red and blue boxes indicate regions of activation and deactivation, respectively, in the within-group analysis.
Resting state functional connectivity
With the left MPFCs as the seeded region, the MTBI group demonstrated significantly (P < 0.01, F = 3.6, cluster size > 150 voxels) higher functional connectivities to the IPL, premotor, and secondary somatosensory cortex (SSC2). On the other hand, the ACC of the healthy group demonstrated significantly (P < 0.01, F = 3.84, cluster size > 150 voxels) less functional connectivities to the IPL and SSC2 than their MTBI counterparts at the resting state (see Figures 3 and 4).
Resting state functional connectivity difference with the left medial prefrontal cortex (seeded region) of the Healthy Controls (N = 15) demonstrating more significant (P < 0.01) functional connectivity to the left secondary somatosensory cortex (SSC2) and right inferior parietal lobe (IPL) than patients with mild traumatic brain injury (N = 15). Resting state functional connectivity difference with the right anterior cingulate cortex (seeded region) of the healthy controls (N = 15) demonstrating less significant (P < 0.01) functional connectivities to the right secondary somatosensory cortex (SSC2) and inferior parietal lobe (IPL) than patients with mild traumatic brain injury (N = 15).

Discussion
Several observations from the current study are worthy of discussion. First, in the within-group analysis, the MTBI patients demonstrated supraspinal activation mainly in regions associated with pain affective response (ACC), but no corresponding response from the modulatory regions when they were exposed to a short period of HP stimuli. On the other hand, the healthy subjects demonstrated robust responses from the supraspinal modulatory regions including the motor and PFCs when they were subjected to the same intensity and duration of noxious stimuli. Further within-group functional causal connectivity assessments of the modulatory response in the two groups indicated the MTBI subjects’ MPFCs cast a much lower number of outward causality inferences to other pain-related supraspinal regions than the healthy subjects. On the other hand, the ACC of the MTBI subjects received a much higher degree of inward causality inference from the sensory and discriminatory pain processing regions in comparison to the healthy subjects. This contrast in functional causality connectivities between the MTBI and the healthy subjects suggests that the supraspinal modulatory functions of MTBI subjects on the affective response to pain are diminished in comparison to the healthy controls. This observation is in line with previous studies in human experimental pain and patients with chronic pain states. In both cases, pain can occur when there is a supraspinal imbalance/dissociation between modulatory and affective responses to pain.16,33,34 Second, in the between-group comparison, the MTBI patients were found to be less reactive in the supraspinal pain modulatory response involving the MPFCs in comparison to the age and gender matched healthy controls. These findings highly correlate with the results of a previous human study demonstrating deficits in prefrontal cortical activities required for cognitive task 35 in the MTBI patients, and other animal studies demonstrating functional deficits in PFCs after TBI resulting in pain regulatory, mood and cognitive dysfunctions.36,37,38 Thus, these current findings further support an emerging assertion that chronic headache is a form of central pain conditions due to damage in brain functions involved in supraspinal pain processing. 39 In addition, increasing evidence suggests that either peripherally or centrally originated pain can become centralized through maladaptive responses within the supraspinal pain network, causing pain perception and mood alterations. 40 Therefore, chronic pain conditions such as intractable persistent headache should be considered as a brain disease. Given that the sample size in the current study is similar to previous published functional imaging studies involving HP, the authors believe these statistically significant findings are highly relevant to the supraspinal impairments caused by the traumatic injury incidences. 3
Second, no significant differences were found between the two groups in their sensory threshold assessments, suggesting that the ascending pain pathway in the patient group was not significantly affected as most of the axonal injury occurs in the supraspinal region. With diffuse axonal injury in the major cortical white matter tracts, MTBI patients are known to have functional deficit in fine motor skill, attention, mood, and memory.41–46 Correlating with these functional deficits is elevated motor evoked potential thresholds found in this patient population. 47 These abnormal electrophysiological findings in the MTBI patients suggest a deficiency in the cortical excitability and conductivity in brain areas associated with pain modulation/adaptation. This recently acquired understanding in MTBI symptom-related neuronal morphological and functional changes provides a feasible direction for treatment development. Already, recent studies have demonstrated that non-invasively stimulating supraspinal regions known to have pain modulatory functions may minimize the headache symptoms in patients with MTBI and other associated post-concussive symptoms.48,49 Thus, applicable non-invasive neuromodulatory or functionally restoring treatment options such as transcranial magnetic stimulation should be considered in this patient population.50,51
While both PTSD and depression are common comorbid conditions in MTBI, the impact of these conditions on pain perception has not been thoroughly defined. A previous study demonstrated that MTBI patients with major depression showed significantly more activities in bilateral amygdalae, thalamus, and prefrontal cortical activities during face matching of fearful stimuli than MTBI patients without depression. Interestingly, while some studies suggested MTBI subjects with major depression demonstrated similar supraspinal functions as those with depression independent of MTBI 52 or PTSD without MTBI, 53 other studies have demonstrated a positive correlation among white matter deficits, default mode network connectivity deficiency with the occurrence of PTSD, and depressive symptoms in the MTBI patients.54,55 To minimize these confounding issues, the current study excluded patients with significant depression and PTSD symptoms meeting the clinical diagnostic criteria. Thus, to fully understand the impact of these comorbid conditions in pain, future studies correlating the severity of these comorbid symptoms to the severities of supraspinal pain processing alteration should be conducted.
The current study consists of several limitations worthy of further discussion. First, the heterogeneity in the mechanisms of injury may minimize the generalization of the observed results to any particular mechanisms of injury. Although the mechanisms of injury may differ, recent studies suggest both blast and non-blast MTBI patients suffer from similar neuropsychological impairments, 56 and there are no significant differences between the two main injury mechanisms in various neuropsychological impairments.57,58 Other studies also suggest TBI itself, independent of injury mechanisms and combat exposure intensity, is a primary driver of adverse outcomes 59 and symptoms such as PTSD and depression are correlated with MTBI independent of the mechanisms of injury. 60 On the other hand, the severity of headache appears to correlate with the severity of depression and PTSD.61–65 Thus, future studies correlating the severity of headache or pain to supraspinal prefrontal cortical activities and psychological symptoms such as depression and PTSD can further enhance our current understanding in their relationship. Second, the chronicity of the headache problem in the studied patient population was quite extensive with the average duration of headache at 70 months and a duration range of 36 to 108 months. This rather wide range of headache durations and an overall high chronicity of the condition precluded the generalization of the result to any acute/sub-acute injury settings. Thus, aside from using the seven-day onset International Classification of Headache Disorder (ICHD-2) diagnostic criteria, future studies may consider assessing the effect of specific increments of chronicity on supraspinal pain processing and modulation. Furthermore, memory dysfunction is a common occurrence in MTBI patient population. Thus, although the investigators conducted due diligence in verifying the patients’ reported date of injury and the onset of headache by carefully reviewing their medical records, minor discrepancies could not be completely excluded. In addition, while the current study demonstrates patients with MTBI appear to suffer from a state of supraspinal pain modulatory functional deficiency, it does not definitively establish the causal effect of the altered state in the development of their headache. Nevertheless, as a first step in characterizing the supraspinal functional deficit associated with MTBI, the current study results do demonstrate with MTBI, supraspinal pain processing and modulation are being affected. Further studies comparing military personnel with non-traumatic headache and/or trauma without headaches are required to further verify the results of these initial mechanistic findings.
Conclusions
In short, patients with MTBI appear to suffer from an altered state of supraspinal modulatory and affective functions related to pain perception. Further assessments are required to correlate the specificity of the current findings to the development of chronic persistent headaches.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors thank the funding support of VA Rehabilitation and Research Development Award.
