fMRI Measures of Evoked Acute Pain - Model-based and Model-free Analysis
Model-based
We reported a temporal and spatial segregation of the brain response to acute painful thermal stimulation [12] (Figure 1). We interpreted the spatial maps associated with an early and late response to noxious stimuli as pathways involving emotional (early) and discriminative (late) processing. A closer look at such results in light of new understanding of brain networks seem to indicate that the main structures of the early response (anterior cingulate/anterior insula) are associated with the salience network while the late response seems to elicit activity in sensory/discriminative structures.
A follow up study by Chen and colleagues [11] of evoked (innocuous) brush and (noxious) thermal stimulation demonstrated clear fMRI signals for noxious heat that do not represent the standard hemodynamic response and that also show a complex response. Innocuous brush stimulation, however, provides a hemodynamic response close to the canonical one. They, however, did not model the signals with more than one hemodynamic response to delineate other brain structures potentially associated with the different signals. Nevertheless, application of sensory stimuli, be it noxious or innocuous, clearly raised the attention of the subjects and seems to be demonstrated in Chen’s paper by a response that peaks about 6 seconds after the initiation of the stimulus. Only for noxious heat, it is possible to observe a second, delayed response appearing about 15 seconds post initiation of the stimulus.
The multiphasic response appears to be independent of the type of noxious stimuli, Ibinson et al. [13] demonstrated a response to noxious electrical stimulation that was multiphasic and similar to the observed one with heat pain.
The data suggest that stimuli that are noxious further induce an additional response that it is delayed with respect to a first one that seems to be independent of the nature of the stimulus. In a recent study, [14] we aimed at assessing test-retest activation to painful stimuli, we utilized a two-phase model to characterize brain activity. Our results indicated that for graded pain stimuli, there’s a differential response in the anterior insula (Figure 2) according to the intensity of the pain stimulus; the early phase, which seems to capture activity in the salience network, displays a graded activity according to the intensity of the nociceptive stimulus suggesting that the salient network encodes pain intensity. Interestingly, the late phase seems to be present in the anterior insula only when the nociceptive stimulus is clearly painful (7/10).
These results were obtained from analyzing fMRI experiments with a traditional general linear model (GLM) approach. Such approach tends to emphasize brain areas with temporal dynamics that better match the model of choice. Given the complexity of the observed response to noxious stimuli in the brain, a model-free approach was utilized to provide an unbiased picture of brain dynamics not obtainable with traditional fMRI analyses.
Model-free Analysis
In order to further explore brain responses to noxious stimuli without imposing a specific response, we performed an independent component analysis, typically used in resting state analysis, of the evoked pain in 10 healthy subjects that were scanned while receiving 4 noxious stimuli of 46°C for 25 seconds with resting periods of 30 seconds at 35°C [15]. fMRI scans were pre-processed using standard approaches and a concatenated ICA analysis (fMRI Melodic Tool) was done with a fixed number of components (70). Following the ICA decomposition of the response of the subjects, brain networks were identified based on published results that characterize different networks and their corresponding connection [16]. Figure 3 displays identified brain networks out of the 70 components obtained in our ICA analysis. The figure further displays a matrix of correlations between the different networks; from the diagonal up absolute correlations are color coded, from the diagonal down partial correlation between networks are displayed (partial correlation takes into account potential correlations of specific networks with the rest of networks) (the figure was generated with FSL’s FSLNets tool). Some immediate, known dynamics are visible in the matrix; for example, the salience network (SN) and the posterior default mode network (pDMN) indicate a strong anticorrelation (r=-0.7) [5,6] and at the same time, the SN and the sensorimotor area representing the hand (rSMN) are strongly correlated. The dynamics indicated here between the DMN, SN and CEN have been extensively studied in resting state experiments [5,6] and are further emphasized here as a result of the noxious stimuli driving the signals across different networks.
We further inspected the temporal responses of specific networks based on the results of Figure 3. In Figure 4, the temporal response in 3 networks of interest; DMN, SN, and SMN are presented as well as a depiction of the explanatory variable of the noxious stimulus used in standard GLM analysis. DMN provides a signal that indicates de-activation when the stimulus was on as observed in other pain and cognitive experiments. The SMN signal seems to better match the GLM model.
The SN presents a striking signal in which it is quickly activated at the beginning of the stimulus returning to baseline and re-activated when the stimulus ends. Furthermore, the initial response is attenuated over time as the subjects become familiar with a painful but not threatening stimulus, suggesting adaptation to the stimuli over time. In a previous, similar study [17] we recorded subjects’ ratings of the 4 noxious stimuli (Figure 5) and, although we did not detect differences between stimuli, there was a trend similar to what we observe in Figure 4 with the second salience peak. Interestingly, the stimulus “offset” response seems to remain constant across all stimuli. We suggest that offset of the noxious stimuli results in pain relief and that, although a subject could have become familiar with the onset of pain, the relief when the stimulus finishes is constant.
fMRI of Evoked Pain in Chronic Pain Patients - Model-free Analysis
We have previously reported a study in pediatric patients suffering from complex regional pain syndrome (CRPS) [18]. In the study, patients were recruited after being diagnosed and scanned before and after a non-pharmacological treatment. The fMRI study was carried out performing evoked stimuli of the affected area producing allodynia (brush stimulation) or hyperalgesia (cold stimulation). Two fMRI scans for each type of stimulation were acquired; during each fMRI scan, the stimulus was delivered for 15 seconds twice with 30 seconds resting time. Pre-processing of fMRI data was performed as described in the original report [18]. Here, we analyzed the brush stimulation using an ICA approach as described above. Spatial maps that identified the SN and DMN were selected and time courses extracted based on a group-concatenated ICA analysis using melodic. Figure 6 depicts the average temporal signal for patients pre- and post-treatment. For the SN, there is a clear initial signal that coincided with the onset and offset of the stimulus before treatment. For post-treatment, the initial SN response was similar to the pre-treatment signal but was mostly absent once the stimulus ended. Similarly, the DMN presented a strong signal when subjects were stimulated pre-treatment and was attenuated after treatment. VAS values pre- and post-treatment were 4.8±0.4 and 1.0±0.3, respectively. We interpret these signals as an indication in the pre-treatment phase that the initial brushing elicited a painful event that was salient upon its termination (pain relief). However, after recovering, although brushing provoked an initial salient response, it was not painful, and its termination did not represent a salient event. Furthermore, the DMN was not suppressed by the salient network as it was determined that this was not a painful event.
The significant attenuation of DMN’s deactivation post-treatment, compared to its response before treatment, may support the notion that this network’s activity is modulated by the saliency of the stimulus. The experimental temporal resolution of the fMRI datasets prevents us from elucidating the dynamical processing that takes place, and we can only speculate that the self-monitoring network (DMN) perceives an event (noxious stimulation) and the salience network assesses and determines is saliency inducing a proportionate disengagement of the DMN while prompting the central executive network to plan for action.
Overall Evaluation of Modeling Evoked fMRI Responses
Converging lines of evidence suggest that the response to evoked painful stimuli does not simply follow the standard GLM approach in fMRI of modeling the hemodynamic response. A complex pattern is observed that requires careful modeling to properly capture the signals and their associated brain structures.
Dissecting the Evoked Response to Acute Nociception in Chronic Pain Patients
Based on the results of Figure 4, we generated 2 explanatory variables (EV) to reproduce the signal observed in the salience network (Figure 4 blue line) for the onset (Figure 7, blue) and offset phases (Figure 7, black) of noxious evoked stimuli. We then constructed EVs for the CRPS patients’ brush response utilizing the onset and offset salience EVs generated form healthy subjects’ data (ref [12]). A standard GLM-based analysis of the brush response for pre- and post-treatment brush scans using FSL feat tools was performed. Differences in activation for the 3 EVs (onset, offset, and standard EV based on stimulus presentation) comparing pre- and post-treatment scans were evaluated and results are presented in Figure 8: Panel A. BRUSH corresponds to the standard “block” EV. For the EV representing the initial salience response (“Salience On” in Figure 8, Panel B), there were no significant differences in activation in pre- vs. post-treatment as evidenced from the group maps and from the comparison map. The results for the EV representing the offset of the stimulus (“Salience Off” in Figure 8, Panel C) displayed significant differences. Comparing pre- vs. post-treatment; the standard BRUSH results seem to indicate that cingulate and frontal areas are more active pre-treatment than post-treatment. Inspecting the Salience-On comparison between pre- and post-treatment indicates non-significant differences in activation, suggesting that although subjects have recovered from chronic pain, their brains remain alert to stimuli of the one affected area. On the other hand, comparison of activation of the Salience-Off response between pre- and post- demonstrates a significant difference due to a smaller response in the post-treatment scan. This result suggests that non-noxious stimuli do not elicit a salience response once the stimulus is finished. Although patients recovered from chronic pain, stimulation of the affected area elicited a salient response, alerting that an area previously sensitive to pain, has been touched, but upon realization that pain is no longer elicited, once the stimulus is stopped, the SN does not find such event as salient.
Summary
We present evidence that brain response to phasic noxious stimuli is a complex inter-dynamic process involving several networks resulting in a fMRI response that is not well characterized utilizing standard models. Using a model-free approach the temporal behavior of the different network provides a picture in which the sequential involvement of these can be delineated and characterized. Accordingly, our results seem to indicate that current approaches do not characterize the brain response to pain in an accurate way resulting in a distorted perception that the brain response mainly represents activation of the salience network. Further studies will provide a better understanding of the brain dynamics of pain processing; specifically related to properly identifying and explaining the reasons for the duration of the onset/offset of the salience response.