Temporal and Intensity Coding of Pain in Human Cortex

Carlo A. Porro1, Valentina Cettolo1, Maria Pia Francescato1, and Patrizia Baraldi2

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
Abstract
Introduction
Methods
Results
Discussion
References

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.

    INTRODUCTION
Abstract
Introduction
Methods
Results
Discussion
References

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; 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).

However, several aspects need to be further clarified. To assess the physiological role of human cortical areas in pain perception, it is necessary to show that cortical neuronal assemblies display activity specifically related to the site, intensity, and time course of perceived pain. This is yet to be conclusively demonstrated. Preliminary PET and functional magnetic resonance imaging (fMRI) studies have shown hemodynamic changes related to pain intensity in different cortical areas (e.g., Davis 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).

To investigate the involvement of human cortical areas in coding temporal and intensive aspects of pain, we used a high-resolution fMRI technique (Frahm et al. 1993; 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

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). 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.

Twenty-nine subjects ("on-line" experiments: 13 ascorbic acid, 8 saline, and 8 controls) were asked to communicate, with both hands, the intensity of perceived pain (0-100) at 1-min intervals during the whole course of the experiments (both before and after stimulation), in response to a light touch of the unstimulated foot. Specifically, they were instructed to lift their hands and fingers twice, the first time to indicate the tens (1-10 lifted fingers = 10-100) and the second one to indicate the units. Closed hands meant 0. The whole motor sequence was usually completed in 3-4 s. Twenty-seven subjects were requested to draw the time profile of the perceived pain intensity immediately after the end of the experiment ("recall" experiments: 11 ascorbic acid, 8 saline and 8 controls). To this purpose, they were instructed to calculate time by paying attention to the rhythmical noises generated by the scanner (2 images were acquired in ~21 s, and a brief but distinguishable pause occurred between each 21-s epoch). In both the recall and on-line experiments the end of pain sensation, if any, was marked by a conventional sign. Apart from these brief motor activities, all subjects were instructed to maintain their head and body still and not to speak throughout the experiments. In subjects injected with ascorbic acid, no difference between the recall and on-line experiments was found in maximum pain intensity (51 ± 6 vs. 45 ± 3, mean ± SE: t = 0.80, P > 0.4), in the area under the pain intensity curve (930 ± 209 vs. 901 ± 108 arbitrary units: t = 0.12, P > 0.9), or in the spatial extent of activated clusters in any of the investigated regions [analysis of variance (ANOVA) group: F = 0.22, P > 0.6; group × region: F = 0.74, P > 0.6]. In both the recall and on-line groups, a linear relationship was found between the extent of neural clusters showing signal changes related to pain intensity and the individual peak pain intensity (r = 0.871 and r = 0.782 for positively correlated clusters; r = 0.836 and r = 0.759 for negatively correlated clusters in the recall and on-line experiments, respectively). Mean normalized fMRI signals, in images acquired at times corresponding to hand signaling in the on-line subjects, were also similar in the two groups, suggesting that the brief hand movements were not associated to marked changes of the pattern of activity within the investigated brain areas (a nonsignificant trend toward higher values was indeed found in the recall subjects: group: F = 3.89, P > 0.06; group × region: F = 0.91, P > 0.5). On these grounds, data from the two groups were pooled together for subsequent analyses.

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). 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.

Data analyses

Functional images were in-plane aligned with the AFNI software package (Cox 1996). 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).

For each subject and each plane, the first image of the functional data set was coregistered to the corresponding 2-D anatomic image. The boundary of cerebral cortex was outlined, excluding areas corresponding to visible artifacts or major blood vessels. Pixels displaying in the first functional image a signal intensity 2 SDs higher or lower than the mean of all pixels lying within the outlined cortical area, which were more likely to be associated with blood vessels or cortical sulci, were excluded from further analyses. Global drifts of signal intensity relative to baseline (mean percentage values at the end of the experimental period 0.57 ± 0.05%, unsigned) were linearly corrected before statistical tests.

The time course of the fMRI signal in each surviving pixel was compared with the individual (in the 40 subjects injected with ascorbic acid or saline) and/or the mean psychophysical pain intensity curve of the ascorbic acid group (all subjects), thus creating statistical maps based on Pearson's correlation coefficient r (possible range -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; lambda  = 0.754). The empiric alpha  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.

To perform regional analyses, boundaries of eight areas of interest (Fig. 1) were outlined onto the individual structural images according to anatomic landmarks, such as the pattern of cortical sulci and the position relative to the anterior and posterior commissures (Picard and Strick 1996; 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 this window]
[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:
[in this window] [in a new window]
 
TABLE 1. Mean centers of mass of clusters significantly correlated with the individual pain intensity curve

The normalized spatial extents of neural clusters were compared with ANOVA, with group as a between-subjects factor and region (the 8 identified areas) and sign (positively or negatively correlated clusters) as within-subjects factors, followed by specific contrasts if appropriate. Psychophysical and heart rate data and the intensity of fMRI signal changes were analyzed by ANOVA and by the Student's t-test. A value of P < 0.05 was assumed as significant.

    RESULTS
Abstract
Introduction
Methods
Results
Discussion
References

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.


View larger version (26K):
[in this window]
[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.

Mean heart rate value during the 4-min baseline period was 73 ± 2 beats/min. A significant increase to 82 ± 3 beats/min (t = 6.46, P < 0.0001) occurred during the minute immediately preceding the onset of stimulation. In almost all subjects, including the ascorbic acid group, this increase rapidly subsided, so that averaged heart rate in successive 5-min periods (0-5 min after stimulus: 73 ± 2; 6-10 min: 72 ± 2; 11-15 min: 71 ± 2) was not significantly different from baseline in any group (group × time, F = 0.50, P > 0.8).

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 -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).

Representative maps from two subjects are reported in Fig. 3. On average, 3.8 ± 0.7% (range 0.03-15.03%) of the total investigated area displayed significant changes; the spatial extent of the response was linearly related to the individual peak intensity of perceived pain (POS clusters: r = 0.85, P < 0.0001; NEG clusters: r = 0.74, P < 0.0001).


View larger version (90K):
[in this window]
[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).

On the basis of a priori hypotheses stemming from previous neurophysiological and functional imaging studies, regional analyses were performed on data from eight selected cortical areas (Fig. 1 and Table 1). Despite interindividual differences in their extent and spatial distribution, positively and negatively correlated clusters tended to concentrate in different regions. Clusters of positively correlated pixels were present in the majority of subjects in the posterior (SI) and/or anterior (MI) portions of the paracentral lobule, supplementary motor area (SMA), medial parietal, and anterior and posterior cingulate cortex. In all these regions, the spatial extent of the activated areas was linearly related to the individual maximal pain intensity. Clusters of negatively correlated pixels were found predominantly in the medial parietal, posterior and perigenual cingulate areas, and in medial prefrontal cortex (Table 1 and Fig. 4). No difference was detected between the spatial extent of the response in the contralateral hemisphere after injection into the right or left foot (side: F = 1.24, P > 0.25; side × region: F = 1.19, P > 0.30).


View larger version (39K):
[in this window]
[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.

The time profiles of fMRI signal changes of POS or NEG clusters in selected brain regions are shown in Fig. 5. In POS clusters, peak signal changes with respect to baseline (during the first 3 min after the s.c. ascorbic acid injection) tended to be higher in SI and in medial prefrontal cortex (on average, approximately +4.2%), but no significant difference was found in the investigated areas. In NEG clusters, peak signal changes were higher in medial prefrontal cortex (on average, -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 this window]
[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.

A small number of clusters (on average 0.5% of the investigated brain area, range 0.05-1.13%) whose signal time course was positively or negatively correlated (P < 0.0005) with the individual psychophysical pain intensity curve were identified in 13 of 16 subjects of the saline group. ANOVA showed that the spatial extent of POS and NEG clusters was significantly higher in the ascorbic acid than the saline group in any of the investigated regions (group: F = 12.03, P < 0.001; group × region: F = 1.44, P > 0.15).

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.

A further correlation analysis was done between the mean psychophysical curve of the ascorbic acid group (Fig. 2A) and the time course of signal intensity in the cortex of volunteers of all three groups (ascorbic acid, saline, and controls). Expectation of a potentially painful stimulus was common to all 56 subjects (see METHODS). However, the overall volume including clusters displaying significant (P < 0.0005) changes was larger in subjects injected with ascorbic acid than the other two groups (ANOVA + difference contrast, P < 0.0001). Regional analyses showed that positively correlated areas were significantly wider in the ascorbic acid than both other groups in SI, SMA, and anterior and posterior cingulate cortex, whereas volumes of negatively correlated clusters were larger in medial parietal, medial prefrontal, and perigenual cingulate cortex (Fig. 6).


View larger version (33K):
[in this window]
[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. bullet , significantly different from the control group, analysis of variance (ANOVA) + Duncan test, P < 0.05; triangle , significantly different from the saline group, ANOVA + Duncan test, P < 0.05.

    DISCUSSION
Abstract
Introduction
Methods
Results
Discussion
References

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; 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).

However, some limitations inherent in the adopted fMRI technique must also be underlined. The in-plane spatial resolution, which is in principle on the order of 2-3 mm, may be hampered by signals coming from large veins, draining from potentially wide brain areas (Lai et al. 1993; 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.

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).

On the other hand, it is worth noting that the identified neural populations displayed significant signal changes (although of lower intensity) even during anticipation of the impending stimulus, suggesting that their activity can be modulated by attentional and emotional factors.

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; 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).

Decreases in brain fMRI signals during acute pain were previously described (Apkarian 1995). 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).

Another interesting issue concerns the relationship between the spatial extent of neural changes and pain intensity. Different sources of evidence suggest that both the frequency of impulse discharges and the total number of central neurons activated are critical for encoding noxious events and discriminating their intensity (Price 1988). 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
Abstract
Introduction
Methods
Results
Discussion
References

0022-3077/98 $5.00 Copyright ©1998 The American Physiological Society