1 Dipartimento Scienze e Tecnologie Biomediche, Università di Udine, I-33100 Udine; and 2 Dipartimento Scienze Biomediche, Università di Modena, I-41100 Modena, Italy
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Porro, Carlo A., Valentina Cettolo, Maria Pia Francescato, and Patrizia Baraldi. Temporal and intensity coding of pain in human cortex. J. Neurophysiol. 80:3312-3320, 1998. We used a high-resolution functional magnetic resonance imaging (fMRI) technique in healthy right-handed volunteers to demonstrate cortical areas displaying changes of activity significantly related to the time profile of the perceived intensity of experimental somatic pain over the course of several minutes. Twenty-four subjects (ascorbic acid group) received a subcutaneous injection of a dilute ascorbic acid solution into the dorsum of one foot, inducing prolonged burning pain (peak pain intensity on a 0-100 scale: 48 ± 3, mean ± SE; duration: 11.9 ± 0.8 min). fMRI data sets were continuously acquired for ~20 min, beginning 5 min before and lasting 15 min after the onset of stimulation, from two sagittal planes on the medial hemispheric wall contralateral to the stimulated site, including the cingulate cortex and the putative foot representation area of the primary somatosensory cortex (SI). Neural clusters whose fMRI signal time courses were positively or negatively correlated (P < 0.0005) with the individual pain intensity curve were identified by cross-correlation statistics in all 24 volunteers. The spatial extent of the identified clusters was linearly related (P < 0.0001) to peak pain intensity. Regional analyses showed that positively correlated clusters were present in the majority of subjects in SI, cingulate, motor, and premotor cortex. Negative correlations were found predominantly in medial parietal, perigenual cingulate, and medial prefrontal regions. To test whether these neural changes were due to aspecific arousal or emotional reactions, related either to anticipation or presence of pain, fMRI experiments were performed with the same protocol in two additional groups of volunteers, subjected either to subcutaneous saline injection (saline: n = 16), inducing mild short-lasting pain (peak pain intensity 23 ± 4; duration 2.8 ± 0.6 min) or to nonnoxious mechanical stimulation of the skin (controls: n = 16) at the same body site. Subjects did not know in advance which stimulus would occur. The spatial extent of neural clusters whose signal time courses were positively or negatively correlated with the mean pain intensity curve of subjects injected with ascorbic acid was significantly larger (P < 0.001) in the ascorbic acid group than both saline and controls, suggesting that the observed responses were specifically related to pain intensity and duration. These findings reveal distributed cortical systems, including parietal areas as well as cingulate and frontal regions, involved in dynamic encoding of pain intensity over time, a process of great biological and clinical relevance.
The involvement of several cortical regions in pain mechanisms is supported by behavioral and neurophysiological observations in experimental animals and by clinical and functional imaging studies in humans (Backonja 1996 Subjects
Fifty-six right-handed volunteers (24 males and 32 females) aged 19-39 yr (mean, 24.3) were included in the study after informed consent. All of them were students or trainees at the School of Medicine, University of Udine, and were free from systemic diseases and painful conditions. To exclude intracranial pathologies and to reduce the potential stress of a novel and restricted environment, each subject underwent a brief training session in the MR system, during which a preliminary set of anatomic images was acquired.
Experimental protocol
The experimental protocol was approved by the Ethics Committee of the University of Udine. Each participant was studied only once and was subjected to only one of three kinds of stimulation on the dorsum of one foot: subcutaneous injection of a commercial buffered ascorbic acid solution (0.5 ml, 20%, pH 6.87; ascorbic acid, n = 24), subcutaneous injection of 0.5 ml of saline (saline, n = 16), or innocuous touching the skin with a needle for 20 s (controls, n = 16). Side of stimulation was balanced across individuals. Subjects were aware of which foot would be stimulated and that the stimulation might be painful. However, they did not know in advance which stimulus would actually occur. Heart rate was continuously monitored in 40 subjects. All volunteers were instructed to rate pain intensity with a numerical 0-100 scale, where 0 denoted "no pain" and 100 denoted "the maximum imaginable intensity of pain." It was explained that the estimate concerned the sensory intensity of pain and not unpleasantness (Price 1988 Data acquisition
The study was performed with a 1.5-T Siemens MR scanner with quadrature head coils. Head motion was minimized by an adjustable padded head holder and orthopedic collars. Field homogeneity was adjusted by a global shimming procedure for each subject, line widths being ~20 Hz. In each volunteer, fMRI data sets were continuously acquired by a FLASH gradient-echo sequence (TR/TE = 63/40 ms, FA = 40°; FOV = 230-245 mm; 128 × 128 acquisition matrix, interpolated to 256 × 256; scan time, 21 s for 2 images) over ~20 min, beginning 5 min before and lasting 15 min after the onset of stimulation, from two sagittal planes (5 mm thick, interplane gap 1 mm) on the medial hemispheric wall contralateral to the stimulated foot. The imaged planes included the cingulate cortex and the putative foot representation area of SI (Fox et al. 1987 Data analyses
Functional images were in-plane aligned with the AFNI software package (Cox 1996
Psychophysical and physiological data
All twenty-four subjects subcutaneously injected with a dilute ascorbic acid solution into the dorsum of one foot reported the occurrence of localized pain with stinging and burning quality, occasionally radiating at the distal anterior portion of the leg during peaks of intensity. Pain lasted 11.9 ± 0.8 min (range 4-15 min), maximum intensity being 48 ± 3 (range 20-85; Fig. 2A). Volunteers receiving subcutaneous saline at the same body site experienced pain lasting 2.8 ± 0.6 min, maximal intensity being 23 ± 4; both values were significantly lower (P < 0.001) than the ascorbic acid group. No subject in the control group reported pain.
Cortical fMRI signal changes related to the individual pain intensity curve
In all 24 subjects of the ascorbic acid group, clusters of pixels whose signal time course was positively or negatively correlated (P < 0.0005) with the individual psychophysical pain intensity curve were identified by cross-correlation analyses. Mean percentage fMRI signal differences from baseline during the first 3 min after the onset of stimulation, corresponding in most subjects to peak pain intensities, were approximately +3.9% for positively correlated (POS) and
Control analyses
Additional analyses were performed to test whether the observed fMRI signal changes were related to factors other than pain intensity and duration. First, a correlation analysis was done between the signal time course of POS and NEG pixels in the ascorbic acid group and the individual heart rate profile during the experimental period. The results showed that, on average, only 10 ± 4% of the POS pixels and 10 ± 4% of the NEG pixels were identified with both the pain intensity and the heart rate curves.
Time is a fundamental dimension of conscious experience, and this study yields the first direct demonstration of a link between the time course of perceived pain intensity and the activity of discrete cortical populations during noxious somatic stimulation lasting several minutes. Neural clusters encoding pain intensity over time were found not only in the somatotopically appropriate portion of the primary somatosensory cortex but also in premotor, cingulate, and medial prefrontal regions. Interestingly, the predominant sign of the response (positive or negative correlations between the fMRI signals and pain intensity) differed in various cortical regions, and the spatial extent of the identified clusters was in several areas linearly related to peak pain intensity across subjects. We hypothesize that these cortical neural populations may contribute to the sensory-discriminative aspects of pain perception and to the graded involvement of attentional, limbic, and motor networks during prolonged pain.
Advantages and limitations of the adopted fMRI technique
fMRI allows individual-based analyses and is potentially characterized by a better spatial and temporal resolution than conventional PET techniques (Cohen and Bookheimer 1994 Specificity of the observed fMRI changes
Several lines of evidence suggest that the changes of activity of the identified cortical clusters after noxious input were specifically related to time-varying changes in pain intensity and not to aspecific arousal or emotional reactions caused by anticipation or presence of pain, or to the activation of nonnociceptive afferents. First, in subjects experiencing pain lasting several minutes discrete, and not widespread, activated clusters were found. Second, the spatial extent of cortical areas showing significant changes was related to the individual maximal intensity of perceived pain. Third, the mean time course of the response was different from that of heart rate, a vegetative index of arousal. Indeed, only a minor fraction of the neural clusters encoding pain intensity could also be identified by cross-correlation with the heart rate profile. Fourth, cortical activity changes of similar entity and time courses were not found after nonnoxious or mild, short-lasting noxious stimulation in subjects who were unaware of the kind of impending stimulus. Were the observed fMRI signal changes, in subjects receiving ascorbic acid, related to stimulation of nonnociceptive slowly adapting receptors or to pain experience per se, irrespective of its magnitude and duration, the mean extent of the cortical response would have been similar in the saline and ascorbic acid groups. It is also unlikely that factors such as memory load or silent counting contributed to the observed results because they had conceivably a different time course than pain intensity (namely, either remained unchanged or increased over time after stimulation).
Cortical circuits showing changes of activity time-related to the perceived pain intensity
Different cortical areas, which are likely to subserve different functions, showed fMRI signal changes significantly related to pain intensity. The activation of the primary somatosensory cortex is consistent with a role for this region in sensory-discriminative aspects of pain (Backonja 1996
INTRODUCTION
Abstract
Introduction
Methods
Results
Discussion
References
; Jones and Derbyshire 1996
; Kenshalo and Willis 1991
; Talbot et al. 1991
). Experimental findings and clinical data suggest that parietal areas, including the primary somatosensory cortex (SI), are mainly involved in the sensory-discriminative aspects (Melzack and Casey 1968
), whereas frontolimbic networks subserve the affective dimension of pain experience (see Albe-Fessard et al. 1985
; Backonja 1996
; Price 1988
). This view is supported by the functional characteristics of neurons projecting to, or situated in SI, many of which show small receptive fields and code the intensity of noxious stimuli, whereas frontal and cingulate neurons usually respond to stimulation of wide body regions (Kenshalo et al. 1988
; Sikes and Vogt 1992
; Snow et al. 1992
). A recent positron emission tomography (PET) study showed indeed that selective manipulation of pain unpleasantness by hypnotic suggestions changes the regional cerebral blood flow (rCBF) response in anterior cingulate but not in SI (Rainville et al. 1997
).
; Svensson et al. 1996), but the relationship between the temporal aspects of perceived pain and cortical activity is unknown. So far, most functional imaging studies in humans provided integrated rCBF measurements during brief periods of noxious stimulation, inducing relatively constant levels of pain intensity, whereas under natural circumstances pain often varies over time. Time-related changes of the activity of CNS cell assemblies during prolonged pain were described in experimental animals (Porro and Cavazzuti 1993
, 1996
).
; Porro et al. 1996
) mainly sensitive to local changes in blood oxygenation accompanying neuronal activation (the so-called BOLD effect) (Cohen and Bookheimer 1994
; Kwong et al. 1992
; Ogawa et al. 1990
). Subcutaneous injection of ascorbic acid was chosen as an ethically acceptable noxious stimulus, inducing a prolonged pain sensation with minimal tissue damage.1 It was recently demonstrated that fMRI may reveal sustained signal changes during prolonged sensory stimulation (Bandettini et al. 1997
), thus allowing the monitoring the activity of cortical populations over time.
METHODS
Abstract
Introduction
Methods
Results
Discussion
References
). No hint was given on the expected intensity and duration of pain. Experiments lasted ~20 min; a 4-min baseline period was followed by a "warning" signal (cleaning the skin with an antiseptic solution) and, 1 min later, by delivery of the stimulus.
). Fifty-seven images per plane were analyzed; the first image of each series was discarded. Anatomical images from the same planes were acquired with a two-dimensional (2-D) T1-weighted sequence (TR/TE = 600/15 ms, 256 × 256). To detect the course of major blood vessels in the two planes of interest and to allow transformation in stereotactic coordinates, respectively, a 2-D FLASH sequence (TR/TE = 32/10 ms, FA = 50°, 256 × 256) and a three-dimensional Turbo-FLASH sequence (TR/TE = 10/4 ms, TI = 300 ms, FA = 10°, 192 × 256) were also run in each subject.
). Twelve volunteers showing clear motion artifacts were not included in the study. In the 56 investigated subjects, mean estimated motion parameters for the x and y directions (expressed in pixel fractions, unsigned) were 0.175 ± 0.017 and 0.384 ± 0.029, respectively; mean rotational parameter was 0.334 ± 0.037°. They were not significantly different in the three groups (ascorbic acid, saline, and controls).
1 to 1) (Bandettini et al. 1993
). Given the relatively long acquisition time for functional images and the sampling period of psychophysical data, no shift was introduced to account for the time lag of the hemodynamic response (Friston et al. 1994
; Kwong et al. 1992
). A value of |r | = 0.45 (corresponding to z = 3.56) and a cluster size = 4 pixels were assumed as significance thresholds, yielding a nominal P < 0.0005 (2-tailed), corrected for the number of comparisons, according to Xiong et al. (1995)
(maximum number of investigated pixels = 18,903;
= 0.754). The empiric
level, estimated with the mean psychophysical curve of the ascorbic acid group in 57 images acquired with the same sequence in phantoms (n = 20), was also <0.05%. To account for variations in brain size and to allow interindividual comparisons, the spatial extents of identified areas were normalized by calculating differences between the individual brains and the standard Talairach brain and correcting values accordingly.
; Talairach and Tournoux 1988
; Vogt et al. 1995
) and then projected onto the aligned functional images. Pixels corresponding to cortical sulci or visible blood vessels were not included in the outlined regions. The foot and leg representation areas of the primary somatosensory cortex (SI) and of the primary motor cortex (MI) were defined as the portions of the paracentral lobule located posterior or anterior, respectively, to the medial root of the central sulcus, which was identifiable in all subjects. The cingulate cortex was divided into three different subareas, namely perigenual (anteroventral), anterior, and posterior cingulate (Devinsky et al. 1995
; Vogt et al. 1995
); the boundary between the anterior and posterior portions was arbitrarily set at the level of the posterior commissure. It is underlined that, given the interindividual variability in the structural and functional parcellation of cortical areas relative to sulcal anatomy (e.g., Rademacher et al. 1993
), the identification of regional boundaries and the assignment of cytoarchitectural areas to the eight regions of interest (Table 1) is by necessity approximate. Functional maps were transformed in the Talairach space along with anatomic images and coordinates of the centers of mass of the identified clusters were calculated for each region of interest.
View larger version (40K):
[in a new window]
FIG. 1.
Schematic drawing of the medial hemispheric wall, modified from the Talairach and Tournoux stereotaxic atlas, showing the location and extent of the 8 cortical areas of interest. AC, anterior cingulate; MI, primary motor cortex; MP, medial parietal cortex (precuneus); MPF, medial prefrontal cortex; PC, posterior cingulate; PGC, perigenual cingulate; SI, primary somatosensory cortex; SMA, supplementary motor area. The horizontal and vertical lines correspond to the bicommissural line and to the vertical planes passing through the anterior and posterior commissure, respectively.
View this table:
TABLE 1.
Mean centers of mass of clusters significantly correlated with the individual pain intensity curve
RESULTS
Abstract
Introduction
Methods
Results
Discussion
References
View larger version (26K):
[in a new window]
FIG. 2.
Time profiles of the mean (±SE, n = 24) intensity of perceived pain (A) and of changes of functional magnetic resonance imaging (fMRI) signal intensity in the cortical populations displaying increases (B) or decreases (C) significantly related to the individual psychophysical pain intensity curve, in subjects injected with subcutaneous ascorbic acid. Arrows: "warning" signal (cleaning the skin with an antiseptic solution). Vertical lines indicate time of injection. fMRI signals were normalized in every subject and for each pixel by dividing the actual value of signal intensity in each image by the mean intensity value of the same pixel in the 1st 10 images acquired during the baseline period.
3.5% for negatively correlated (NEG) clusters; these signal changes gradually declined toward baseline along with decreases in the intensity of perceived pain (Fig. 2, B and C). In the same neural populations, significant signal changes (approximately +1.4% for POS and
1% for NEG clusters: t = 3.51 and t = 4.78, P < 0.01 and P < 0.001, respectively) occurred even during the 40-s period after the "warning" signal, before stimulation onset (Fig. 2, B and C). However, mean fMRI signals during peak pain intensity were significantly different both from baseline and from the postwarning period (ANOVA + difference contrast: t = 7.12, P < 0.0001 and t = 11.89, P < 0.0001 for POS and NEG clusters, respectively).
View larger version (90K):
[in a new window]
FIG. 3.
Statistical maps from 2 representative subjects, showing the spatial distribution of neural clusters displaying positive (red-yellow) or negative (green-blue) correlations with the individual pain intensity curve. The blue arrow on top points to a cluster in the postcentral gyrus (SI).
View larger version (39K):
[in a new window]
FIG. 4.
Maximal pain intensity plotted against the normalized spatial extent (see METHODS) of neural clusters showing positive or negative correlations with the individual curve of perceived pain intensity, in subjects injected with ascorbic acid. Lines show linear best fits.
4.6%) than all other regions, with the exception of SI and MI (P < 0.05, corrected for the number of comparisons).
View larger version (25K):
[in a new window]
FIG. 5.
Time profiles of the mean changes of fMRI signal intensity in neural clusters displaying increases (SI, anterior cingulate) or decreases (perigenual cingulate, medial prefrontal cortex) significantly related to the individual psychophysical pain intensity curve, in subjects subcutaneously injected with ascorbic acid. Arrows: "warning" signal (cleaning the skin with an antiseptic solution). Vertical lines indicate time of injection.
View larger version (33K):
[in a new window]
FIG. 6.
Comparison of the normalized spatial extent of regions displaying fMRI signal changes significantly related to the mean pain intensity profile of the ascorbic acid group, in volunteers subjected to different kinds of stimuli. Pos and Neg: clusters of pixels displaying positive or negative correlations, respectively. , significantly different from the control group, analysis of variance (ANOVA) + Duncan test, P < 0.05;
, significantly different from the saline group, ANOVA + Duncan test, P < 0.05.
DISCUSSION
Abstract
Introduction
Methods
Results
Discussion
References
; Jones and Derbyshire 1994
). Thus, fMRI maps may reflect the dynamic behavior of neural populations over time and the different pattern of response of adjacent cortical networks. Individual-based analyses, which do not require extensive spatial smoothing, also allow to quantify the extent of activated areas and relate it to psychophysical data. The interindividual variability of the response may in part be responsible for conflicting results obtained in PET studies when group analyses were performed (Jones and Derbyshire 1994
; Porro and Cavazzuti 1996
).
; Segebarth et al. 1994
). In this study, care was taken to exclude from the analyses pixels corresponding to visible blood vessels or cortical sulci, where large veins are often encountered (e.g., Yousry et al. 1996
). Although intraparenchimal small veins and venules obviously contribute to the observed signal changes, the activation of adjacent cortical areas such as MI and SI may be clearly differentiated with similar fMRI techniques (Porro et al. 1996
). The temporal resolution of fMRI is otherwise limited by the time constant of the hemodynamic response (Friston et al. 1994
; Kwong et al. 1992
) and by the hardware and software configuration of the scanner. The main disadvantage of the adopted fMRI technique was that signals from only two planes could be imaged at a reasonable temporal resolution. We focused on the hemisphere contralateral to the injected site, where most rCBF changes are known to occur during brief noxious input (Casey et al. 1994
, 1996
; Coghill et al. 1994
; Hsieh et al. 1995a
; Jones et al. 1991
; Talbot et al. 1991
). However, additional fMRI studies characterized by higher acquisition speed (Kwong 1995
) will allow to include areas not investigated in the present work, such as the anterior insula and secondary somatosensory cortex (Coghill et al. 1994
; Dong et al. 1994
; Jones et al. 1991
; Talbot et al. 1991
) and the ipsilateral hemisphere, yielding a more comprehensive view of brain regions involved in temporal and intensity coding of pain.
; Kenshalo and Willis 1991
), at least at suprathreshold intensity levels. However, pain intensity coding appears to occur in parallel cortical channels receiving inputs either from the "lateral" ascending system (parietal cortex) or the "medial" system (e.g., anterior cingulate) (Albe-Fessard et al. 1985
; Craig 1994
). Information regarding pain intensity may be conveyed to anterior cingulate by medial thalamic neurons (Bushnell and Duncan 1989
), and, indeed, a linear relationship between peak hemodynamic changes and pain intensity during brief transcutaneous electrical stimulation of the median nerve was reported in a recent fMRI work in the posterior part of the anterior cingulate cortex (Davis et al. 1997
). The anterior cingulate includes different subregions and may be associated with several functions, including affective aspects of pain (Devinsky et al. 1995
; Picard and Strick 1996
; Rainville et al. 1997
). Interestingly, the mean coordinates of sites showing significant changes in this study closely correspond to foci of activation found in PET studies during acute noxious stimulation (Casey et al. 1994
; Jones et al. 1991
; Talbot et al. 1991
), a perceived illusion of pain (Craig et al. 1996
) or chronic neuropathic pain (Hsieh et al. 1995b
). The presence of activated clusters in primary motor and medial premotor (SMA) areas confirms and extends the results of previous PET studies (e.g., Casey et al. 1996
; Coghill et al. 1994
; Hsieh et al. 1995a
; Peyron et al. 1998
) and it may be related to the organization (or inhibition) of motor reactions to noxious input (see also Hsieh et al. 1994
).
). Although their exact physiological significance is yet to be clarified, it is worth mentioning that several areas showing mainly negatively correlated clusters, such as the ventromedial prefrontal, medial parietal, and posterior cingulate cortex, also displayed reduced rCBF during acute (Coghill et al. 1994
; Vogt et al. 1996
) and chronic (Hsieh et al. 1995b
) pain as well as during active visual tasks relative to passive viewing of the same stimuli (Shulman et al. 1997
). These hemodynamic changes may reflect shifts of activity in structures involved in monitoring the environment or the self (Damasio 1994
). Neural changes in the medial prefrontal-perigenual cingulate cortex are of particular interest, given the extensive connections of these regions with components of the limbic system (Devinsky et al. 1995
; Fuster 1989
) and their likely involvement in mood regulation (Drevets et al. 1997
; George et al. 1995
; LeDoux 1996
; Mayberg 1997
). Interestingly, pain intensity and unpleasantness may closely covary during different kinds of pain (Price 1988
; Wade et al. 1996
). Our results support the hypothesis that encoding of pain intensity may be a major factor in the regulation of emotional and autonomic responses (Duncan et al. 1992
).
). At the spinal cord level, it has been shown by 2-deoxyglucose autoradiography in the rat that nonnoxious stimulation is associated with more delimited increases in glucose metabolic rates than noxious input; moreover, higher intensities of noxious stimulation induce wider metabolic changes (Coghill et al. 1993
). Our demonstration of a relationship between the spatial extents of the identified clusters in specific cortical regions and peak intensities of perceived pain, and of a greater extent of the response in the group experiencing higher pain levels, is in line with the hypothesis that neural recruitment mechanisms are important for pain intensity coding.
![]() |
ACKNOWLEDGEMENTS |
---|
The authors thank P. E. di Prampero and M. Bazzocchi for valuable support to the project and L. Bon, M. E. Diamond, D. Mostofsky, P. Nichelli, and G. Pagnoni for helpful discussions.
Supported by grants from Consiglio Nazionale delle Ricerche and Ministero dell'Università e della Ricera Scientifica e Technologia, Italy, to C. A. Porro.
![]() |
FOOTNOTES |
---|
1 A preliminary report on the results from 18 subjects was presented at the 3rd International Meeting on Human Brain Mapping (C. A. Porro et al., Neuroimage 5: S224 1997).
Address for reprint requests: C. A. Porro, Dip. Scienze e Tecnologie Biomediche, Università di Udine, P.le Kolbe 5, I-33100 Udine, Italy.
Received 3 June 1998; accepted in final form 5 August 1998.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|