Abstract

Evidence from both human and animal studies has demonstrated a key role for brainstem centers in the control of ascending nociceptive input. Nuclei such as the rostral ventromedial medulla and periaqueductal gray (PAG) are able to both inhibit and facilitate the nociceptive response. It has been proposed that altered descending modulation may underlie many of the chronic pain syndromes (both somatic and visceral). We used functional magnetic resonance imaging to image the neural correlates of visceral and somatic pain within the brainstem. Ten healthy subjects were scanned twice at 3 tesla, during which they received matched, moderately painful, electrical stimuli to either the midline lower abdomen or rectum. Significant activation was observed in regions consistent with the PAG, nucleus cuneiformis (NCF), ventral tegmental area/substantia nigra, parabrachial nuclei/nucleus ceruleus, and red nucleus bilaterally to both stimuli. Marked spatial similarities in activation were observed for visceral and somatic pain, although significantly greater activation of the NCF (left NCF, p = 0.02; right NCF, p = 0.01; Student's paired t test, two-tailed) was observed in the visceral pain group compared with the somatic group. Right PAG activity correlated with anxiety during visceral stimulation (r = 0.74; p < 0.05, Pearson's r, two-tailed) but not somatic stimulation. We propose that the differences in NCF and right PAG activation observed may represent a greater nocifensive response and greater emotive salience of visceral over somatic pain.

Introduction

Since 1906, we have known that spinal cord excitability is directly influenced by descending, tonically active signals that are mostly inhibitory in function. Work by Sherrington (1906) showed that nociceptive reflexes were enhanced after the spinal cord was transected. This supraspinal descending modulation was demonstrated graphically by World War I soldiers with severe injuries who reported disproportionately low pain sensation (Beecher, 1946), presumably because of mechanisms of distraction or stress-induced analgesia. Distraction from pain decreases both the intensity and affective components of the pain experience (Miron et al., 1989; Bantick et al., 2002). Therefore, cortical and subcortical circuitry exists to allow pain modulation in a top-down manner. Over the last 35 years, an accumulation of evidence has demonstrated key roles for brainstem centers in this control of nociception. In a pivotal study, Reynolds (1969) highlighted the role of the periaqueductal gray (PAG) in the nocifensive response: electrical stimulation of the PAG in rats allowed abdominal surgery without the use of general anesthesia. This inhibition of the normal pain response through stimulation of discrete brainstem centers such as the PAG [stimulus-produced analgesia (SPA)] was subsequently described in humans (Boivie and Meyerson, 1982; Baskin et al., 1986).

In animals, SPA has also been demonstrated in response to stimulation of the rostral ventromedial medulla (RVM) (Basbaum et al., 1976; Fields et al., 1977), which is situated centrally around the pontomedullary junction. Exposure of this region to both higher intensities of electrical stimulation and higher concentration of stimulatory neurotransmitters can induce analgesia (Zhuo and Gebhart, 1990; Gebhart, 2004). The RVM has a dual role in pain control: it is able to both inhibit and facilitate the nociceptive input. The latter is achieved through low-intensity electrical stimulation or low concentrations of stimulatory neurotransmitters. Facilitation of pain via the RVM has been implicated in the development of central sensitization and secondary hyperalgesia (Urban and Gebhart, 1997, 1999). It has been proposed therefore that disordered descending influences from the brainstem on nociceptive afferent information may underlie abnormal pain perception in both functional pain disorders, such as irritable bowel syndrome and fibromyalgia, and neuropathic pain (Porreca et al., 2002).

To establish whether aberrant brainstem pain processing underlies human chronic somatic or visceral pain syndromes, noninvasive methods of investigating brainstem function need to be developed. Technical difficulties in functionally imaging the brainstem, such as poor spatial resolution, local field inhomogeneity-induced signal losses, and image distortion have precluded any pain study in humans. For a review of these issues and of which studies to date have identified brainstem activation using neuroimaging methods, see Tracy and Iannetti (2005). Only a few studies to date have successfully identified regional brainstem activity in individual subjects during pain, fine touch, swallow, motor activities and taste (Komisaruk et al., 2002; Tracey et al., 2002; Topolovec et al., 2004). The aims of this study were to compare the regional brainstem activation between visceral and somatic pain. We therefore performed high-resolution functional magnetic resonance imaging (fMRI) in healthy subjects during both rectal and lower abdominal electrical stimulation matched for pain intensity and unpleasantness.

Materials and Methods

Subjects. Ten healthy, right-handed subjects (five female) participated in the study. The median age was 30 years of age (range, 21-33 years). None reported any abnormal gastrointestinal or somatic symptoms on a personal history and bowel symptom questionnaire (Talley et al., 1989). Clinical depression was excluded with the use of the Beck's Depression Inventory (Beck et al., 1961). None of the subjects took any medications likely to interfere with the blood oxygen level-dependent response. All subjects gave informed consent, and the Oxfordshire Clinical Research Ethics Committee approved the study.

Stimuli. Electrical stimuli were used to induce both somatic and visceral pain. Purpose-built electrical stimulating catheters were designed for safe use within the high field strength magnet (Standard Instruments, Karlsruhe, Germany). These contained no ferromagnetic material and had two bipolar ring electrodes 2 cm apart at the distal tip. The catheters were connected via a magnet-compatible extension cable to a Digitimer (Welwyn Garden City, UK) DS7A constant current stimulator. The DS7A automatically varies the voltage depending on the impedance in the circuit to provide a constant, predefined current. A Digitimer DG2 trigger generator was used to provide 6 s trains of electrical pulses delivered at a frequency of 5 Hz.

To induce visceral pain, the catheters were inserted into the rectum, 10 cm from the anal canal. For somatic pain, the catheters were attached to the lower abdomen, in the midline at the dermatome level of T12.

Study design. Before the functional scan, increasing 6 s stimuli were delivered to the subject in a stepwise manner to define the current required to induce a reliable pain intensity rating of 5 out of 10. The subjects were instructed that a rating of 0 meant no pain and a rating of 10 represented excruciating pain. This current was used for the remainder of the experiment. During the functional scan, 20 6-s stimuli were delivered to the subject with an average 60 s interstimulus interval. The subjects rated each stimulus via a numerical rating scale projected onto a screen that was visualized via prism glasses. An arrow on the scale could be moved via a button box. Unpleasantness was also rated (0, not unpleasant; 10, unbearably unpleasant) along with urge sensation (0, no urge; 10, intense urge). All scales were presented to the subjects during the interstimulus interval. The order of the visceral and somatic scans was randomized. Subjects were then asked to retrospectively rate their average anxiety levels during the stimuli (0, no anxiety; 10, highly anxious).

Data acquisition. For the functional data, subjects were scanned in a 3 tesla human MRI system (Oxford Magnet Technology, Oxford, UK) with the use of a Magnex SGRAD MK III head coil (Magnex Scientific, Oxford, UK). To minimize the image distortion induced through B0 magnetic field inhomogeneities in the vicinity of the brainstem, coronaloblique slices were acquired in the plane of the anterior wall of the fourth ventricle, graphically shown in Figure 1. Twenty-four 2 mm slices with an in-plane resolution of 3 × 3 mm were obtained with the following parameters: repetition time (TR), 3 s; echo time (TE), 30 ms; flip angle, 90°; field of view, 192 × 192 mm; matrix, 64 × 64. These parameters were determined after a series of pilot functional scans. They resulted in the most reliable brainstem activation across subjects at a resolution that was high enough to be able to identify individual brainstem nuclei. The coverage was limited to the brainstem and a segment of cortex extending from the superior colliculus posteriorly to the hypothalamus anteriorly but incorporated the full width of the brain. To minimize B0 field inhomogeneity-induced effects, particular attention was given to optimal shimming before the commencement of the functional imaging run. A whole-brain structural proton-density turbo spin echo (PDTSE) scan was acquired in an equivalent coronal-oblique plane on a Siemens (Munich, Germany) Magnetom 1.5 T scanner. The parameters used were as follows: turbo factor, 9; effective TE, 12 ms; TR, 6 s; slice thickness, 2 mm; in-plane resolution, 0.75 × 0.75 mm; three averages. To quantitatively evaluate the intensity of activation within individual brainstem nuclei, region of interest (ROI) masks were produced. To facilitate mask production, an additional structural, axially acquired PDTSE brainstem scan was obtained with maximal brainstem contrast (Fig. 2). Parameters for this scan were as follows: turbo factor, 3; effective TE, 12 ms; TR, 5 s; slice thickness, 3 mm; in-plane resolution, 0.8 × 0.8 mm; six averages.

Graphical representation of the slice orientation during the functional scan. For clarity, the underlying image is a T1-weighted structural scan, and the number of slices was reduced from 24 to 8. A coronal-oblique orientation was used with slices aligned with the ventral wall of the fourth ventricle.

Image analysis. Image analysis, aimed at delineating significant brain activation via changes in the blood oxygen level-dependent signal (Ogawa et al., 1992), was performed on each subject's functional data set using the Expert Analysis Tool (www.fmrib.ox.ac.uk/fsl/) of the Centre for Functional Magnetic Resonance Imaging of the Brain (FMRIB) (Smith et al., 2001). Before the statistical analysis, brain segmentation was performed with the use of FMRIB brain-extraction tool (BET) (Smith, 2002). All scans were scrutinized after BET, procedure and additional manual segmentation was performed if required, thus optimizing the accuracy of any registrations. The data were motion corrected with the use of MCFLIRT (Jenkinson et al., 2002), spatial smoothing was performed with a Gaussian kernel of 4 mm full-width at half-maximum, intensity normalization was performed with a single scaling factor, and high-pass temporal filtering was performed with a Gaussian-weighted least-squares straight-line fit and a high-pass cutoff filter of 60 s. The statistical analysis was performed with the improved linear model of FMRIB (Woolrich et al., 2001). A model of the relevant applied stimuli was thus designed and convolved with the hemodynamic response function. This convolved model was then fitted to the four-dimensional data set to demonstrate areas of brain activation.

Axial slices through a PDTSE sequence structural scan of the brainstem. These were designed to provide maximal resolution and contrast within the brainstem to aid accurate region of interest mask formation. The substantia nigra (SN) and PAG are seen in lighter contrast. The red nucleus (RN), dentate nucleus (DN), decussation of the corticospinal tracts (DCST), and inferior olivary nucleus (ION) are labeled.

Each voxel was analyzed against the convolved model with a resultant parameter estimate (PE) image. The parameter estimate is proportional to the signal change arising from the stimulus. Cluster thresholding with significance estimation defined by Gaussian Random Field Theory was used to identify clusters of activated voxels (z-score > 2.3; p < 0.01). The functional data set was coregistered onto the subject's high-resolution whole-brain PDTSE scan, which was then registered onto a standard brain [Montreal Neurological Institute (MNI) 152 brain] to facilitate group analysis.

A mixed-effects group analysis was performed for each of the two groups in standard brain space with Local Analysis of Mixed Effects (of the FMRIB). This incorporates variance within session and across time (fixed effects) and cross session variances (random effects). Cluster thresholding was performed with a z-threshold of 2.3 and corrected p value of < 0.01 (Worsley et al., 1992; Friston et al., 1994).

Left and right ROI masks were then defined for the PAG, RVM, nucleus cuneiformis (NCF), dorsolateral pons [the parabrachial nucleus (PBN) and nucleus ceruleus (NC)], and the ventral tegmental area (VTA). Because no stereotaxic atlas exists for the brainstem, the ROIs were defined by careful comparison with a detailed brainstem atlas (Duvernoy, 1995) using both T1-weighted structural scans and the axially acquired PDSE structural scan with optimal contrast within the brainstem (Fig. 2). The ROIs were then transformed into standard MNI space, and this transformed mask was applied to the individual's activation, also within standard MNI space. This method allows us therefore to take into account any subtle variations in brainstem anatomy between subjects. The mean PE was then calculated for each ROI for each subject's data set. Significant differences between the somatic and visceral mean group PE for each ROI were tested with a two-tailed paired Student's t test (not corrected for multiple comparisons). Because the stimuli used were midline, we would expect any significant effects in the ROI analysis to be bilateral. This being the case, we contend that this reduces the chances of the bilateral effects we found being falsely reported as positive. Also these brainstem ROIs are part of a well recognized ascending and descending pain network and are not independent and often coactivate. Therefore, if the typical Bonferroni's correction were applied, it would be excessively stringent given its requirement for needing no a priori hypotheses and total independence of behavior between regions of interest. An independent analysis of the likelihood of finding bilateral significant differences of right and left NCF parameter estimates between the two stimulation paradigms by chance supports our reasoning and validity of results (data not shown). Correlation analysis of the mean parameter estimate between these regions was performed using Pearson's r correlation analysis.

Imaging

A template, stereotaxic brain, such as the atlas by Talairach and Tournoux (1988), does not exist for the brainstem. Identification of the specific regions activated therefore depended on careful identification of the areas of the brainstem with the use of a detailed accurate atlas (Duvernoy, 1995). Of the 10 subjects studied, activation within the brainstem (mesencephalon, pons, and/or medulla) occurred in seven subjects during both visceral and somatic pain (Table 1). In the remaining three subjects, brainstem activation was observed in just one of the two scans (two visceral and one somatic). The group activation maps for all 10 of the subjects are shown in Figure 4. A similar spatial pattern of activation can be seen for both sensory modalities. Bilaterally, the posterior insula cortex and thalamus (mediodorsal and ventrolateral thalamic nuclei) were activated in both groups, regions frequently activated in pain imaging studies.

The group psychophysical data. The stimuli were well matched, with no significant differences in either intensity or unpleasantness rating between visceral or somatic stimulation (p > 0.05; Student's t test, paired, two-tailed). Electrical stimulation of the rectum induced a mild urge sensation. Error bars represent SE.

In addition, there was concordance between both groups in the bilateral activation of the red nucleus, NCF, and PAG. An additional area of activation occurred bilaterally in the dorsolateral pons, anterolateral to the fourth ventricle. This region is consistent with locus ceruleus medially and parabrachial nuclei laterally and was common to both visceral and somatic pain. Both sensory modalities activated a region on the right side of the ventral mid pons. This area includes but extends beyond (mainly anteriorly) adrenergic center A5 in the visceral group but is situated more anteriorly in the somatic group. This peak of activation overlies the right pontine nuclei (PN). Left-sided PN activity was also seen in the visceral group. The PN act as relay points for the corticoponto-cerebellar pathway. These areas of activation are consistent with the projection from the primary and secondary sensory cortices (Brodal, 1968; Leergaard et al., 2004). A region in the inferior, anterior mesencephalon was activated in both sensory modalities bilaterally. This region encompasses the VTA and extends laterally into the substantia nigra. A region in the left dorsolateral region of the superior medulla was activated solely during visceral pain. This region contains the nucleus tractus solitarius NTS), gracile nucleus (GN), and dorsoreticular nucleus, key regions in visceral and nociceptive processing. In several regions, areas of activation overlap. Despite this, separate peaks of activation are present in each area (Fig. 4, red nuclei and nucleus cuneiformis). In these cases, the voxel with the peak z-score has been used to identify the location of the individual nuclei. In addition, some of the nuclei discussed lie adjacent to each other (such as the PAG and NCF). Thus, even with higher resolution, clear demarcations between all nuclei would not be possible.

Group activation maps for somatic (first column) and visceral (second column) pain. The sagittal slices show columns of activation within the brainstem in each group. Axial slices (a-e) correspond to the plane indicated in the sagittal slice. A reduced field of view was used during functional scanning; thus, the activation map is limited in its anteroposterior plane. The two groups have a similar spatial pattern of activation. Regions commonly activated in whole-brain pain-imaging studies are significantly activated: thalamus (Thal) and posterior insula (Ins) bilaterally. Activation was also seen in the globus pallidus (GP) bilaterally in both groups. Significantly, activated brainstem nuclei included the red nucleus (RN), NCF, PAG, VTA [which extended laterally into the substantia nigra (SN)], and the dorsolateral pons (DLPons) bilaterally. Bilateral activation of the PN occurred in the visceral group but was limited to the right side in the somatic group. Activation in the region of the RVM was also seen in both groups. In the visceral group alone, a region of activation occurred in the left dorsolateral medulla (DLM). A small area of motion artifact (MA) occurred around the central canal in the somatic group. The axial slices in the visceral group have been magnified and juxtaposed to drawings at the same anatomical location [modified from Duvernoy (1995) with permission].

Because the stimuli used were midline, we would expect any significant effects in the ROI analysis to be bilateral. This is the case for NCF parameter estimate comparison between visceral and somatic pain and correlations of dorsolateral pons and RVM activity during visceral pain. Likewise, the correlation of PAG activity with anxiety during visceral pain shows a trend toward a bilaterally significant result. We contend that this dramatically reduces the chances that these effects are being falsely reported as positive. These values have therefore not been corrected for multiple comparisons because use of Bonferroni's correction would be excessively stringent.

Discussion

The results from this study demonstrate that it is feasible to perform an fMRI of the brainstem (from medulla to mesencephalon) in individual subjects, thus supporting previous studies (Komisaruk et al., 2002; Tracey et al., 2002; Topolovec et al., 2004). Detailed attention to accurate brain segmentation and registration has allowed us to further these observations by performing group analyses. These have shown bilateral regional brainstem activation in areas consistent with the PAG, RVM, NCF, and the dorsolateral pons. Furthermore, the pattern of spatial activation is similar for the two sensory modalities, visceral and somatic. These results are consistent with animal stimulation and tracing studies (Traub et al., 1996; Rodella et al., 1998; Millan, 2002; Monnikes et al., 2003) and the marked viscerosomatic convergence that occurs at the spinal cord level (Ruch, 1946). The main spatial difference between the two sensory modalities is activation of the dorsolateral medulla (including the gracile nucleus) in the visceral group only. We also identified differences in behavior of individual regions: greater quantitative activation bilaterally in the NCF and a correlation of anxiety with right PAG activity during visceral pain alone.

Nucleus cuneiformis

Generally, we do not feel peristalsis (passage of stool or gas), yet in rodents visceral afferents are excited by non-noxious, physiological stimuli (Andrew and Blackshaw, 2001), and subliminal visceral stimulation in humans results in cortical activation, as seen with fMRI (Kern and Shaker, 2002). Therefore, information must be reaching the dorsal horn or brainstem (via the vagus nerve). Why then do we not, in general, perceive the normal functioning of our gastrointestinal tracts? One possible explanation is that visceral input could be under a persistent inhibitory control. Undoubtedly, constant perception of gut activity would act as a considerable distraction to daily activity. Our results have shown that the NCF, which is known to have nocifensive functions, is activated to a significantly greater degree during visceral pain than somatic pain. Therefore, visceral pain may result in a heightened, descending inhibitory control from brainstem centers such as the NCF, thus preventing sufficient input to the cortex to allow conscious perception of afferent information from the bowel. The NCF, a nucleus of the brainstem reticular formation contains functionally distinct classes of neurons that are able to enhance nociception (“ON cells”) and inhibit nociception (“OFF cells”) (Haws et al., 1989), characteristics it shares with the PAG and RVM (Fields et al., 1983; Heinricher et al., 1987). It has strong neuronal connections with the PAG, RVM, PFC, hypothalamus, and amygdala (Edwards and de Olmos, 1976; Zemlan and Behbehani, 1988; Bernard et al., 1989; Sesack et al., 1989), and c-Fos immunoreactivity in this region has been induced by colorectal distension in rats (Monnikes et al., 2003). Therefore, it is tactically positioned to receive emotive, autonomic, and executive inputs, in concert with the PAG and RVM, and can therefore influence nociceptive processing in a facilitatory or inhibitory manner. Indeed, a recent human fMRI study of mechanical stimulation to experimentally induced secondary hyperalgesia demonstrated a significantly greater left NCF activation compared with mechanical stimulation in the control stimulation condition, thus illustrating its potential facilitatory role in the development of central sensitization (Zambreanu et al., 2005). The absence of a correlation with anxiety ratings suggests that the increased NCF activity seen in the visceral group is not related to the heightened anxiety observed.

Periaqueductal gray

Subjects experienced higher anxiety during the visceral stimulation compared with somatic stimulation. Furthermore, anxiety ratings correlated with right PAG activation during visceral but not somatic pain (Fig. 6). Anxiety and the stress response are coordinated through the PAG together with other regions such as the hypothalamus and amygdala. William James in 1884 proposed that emotions are dictated by and reliant on bodily sensations (James, 1884). Modern proponents of these theories comment that the physical sensation of anxiety, or any emotional response, is dependent on sensations from the viscera (Damasio, 1996). This has presumably led to the aphorisms “gut feelings,” “butterflies in the stomach,” and “a visceral response.” These emotively induced visceral sensations may be attributable to changes in autonomic outflow and resultant increased gut tone (Almy et al., 1949; Whorwell et al., 1992). However, the consequential visceral sensation induced by these changes in gut tone will transmit emotive meaning. Therefore, the correlation of visceral, but not somatic, pain-induced PAG activity with anxiety ratings could represent the neural correlates of visceral-based perception of anxiety [i.e. peripheral visceral signals preferentially (over somatic inputs) connect with the PAG to contribute to the sensation of anxiety perceived by the subject].

These theories are speculative, alternative explanations are possible. For example, pain perception can be increased through high anxiety (Ploghaus et al., 2001) or decreased through stress-induced analgesia (Fields, 2000). Thus, the correlation of right PAG activity with visceral pain may reflect an emotively/autonomically induced descending modulation of the noiceptive input. It is impossible to ascertain whether this modulation is facilitatory or inhibitory in our results as the pain intensity of the stimuli was locked.

Dorsolateral medulla

An additional difference in the observed activations was a region of activation in the left dorsolateral medulla in the visceral group only (Fig. 4). This area includes important, nociceptive nuclei: the NTS and the GN. The NTS receives a predominantly visceral input (nociceptive and non-nociceptive) from the vagal nerve and spinal afferents. Its most caudal aspect, extending the length of the medulla, receives input from the visceral organs and vasculature, with the superior pole receiving gustatory neurons. Electrical stimulation of the nucleus can induce both inhibition (Aicher and Randich, 1990) and facilitation (via the vagal nerve) (Randich and Gebhart, 1992) of nociception in rodents. It therefore would be expected to contribute to the measured activation in response to nociceptive visceral stimulation, as found in our study. The GN receives a mainly visceral nociceptive input from the dorsal column pathway (Willis et al., 1999). It is thought to represent an ascending part of an amplification loop that, when activated, heightens visceral nociceptive perception possibly, although facilitatory influences of the RVM (Palecek, 2004). Another nucleus found in this region of activation is the dorsal reticular nucleus of the medulla. This receives both visceral and somatic input and is both pronociceptive and nocifensive (Lima and Almeida, 2002). Studies suggest that this area is crucial for the development of diffuse noxious inhibitory controls (Villanueva and Le Bars, 1995). Future high spatial-resolution studies will aim to resolve these structures functionally to determine their separate contributions to visceral nociceptive processing.

Ventral tegmental area

These results also suggest a role for the VTA in both visceral and somatic pain processing. This region has been studied extensively in studies of addiction and reward in human and animal studies and represents part of the mesolimbic reward circuitry. It has close connections with the limbic system, especially the cingulate cortex (Oades and Halliday, 1987), and rodent studies have shown that connectivity exists with the PAG and RVM (Kirouac et al., 2004). Electrical stimulation induces an antinociceptive effect in rats (Sotres-Bayon et al., 2001). Furthermore, nondopaminergic neurons exist within the VTA that respond to aversive, not rewarding stimuli (Ungless et al., 2004), suggesting a role of the VTA in the affective/aversive dimension of pain processing.

RVM: dorsolateral pons correlation of activation

A significant correlation between left and right dorsolateral pons (incorporating the PBN and NC) and RVM activity occurred during visceral pain. The spinoparabrachial pathway provides a major tract for ascending nociceptive information (Hunt and Mantyh, 2001). More than 80% of lamina 1 dorsal horn projection neurons to the PBN express the neurokinin-1 (NK-1) receptor (Todd et al., 2000). Selective ablation of lamina 1 neurons with NK-1 receptors results in a marked attenuation of capsaicin-induced thermal and mechanical primary and secondary hyperalgesia in rodents (Mantyh et al., 1997; Suzuki et al., 2002). Inactivation of the RVM can abolish the development of secondary hyperalgesia in rodents (Urban and Gebhart, 1999). A spino-bulbo-spinal loop to and from lamina 1 to the PBN and back to the RVM to either facilitate or inhibit nociceptive input at the level of the dorsal horn has been postulated (Suzuki et al., 2004). This loop could explain the correlation observed in our study. Alternatively, this observed correlation may be attributable to similar activity of the NC and RVM. The NC also has descending inhibitory functions, electrical stimulation induces antinociception in the rat (Jones and Gebhart, 1986). Therefore, descending nociceptive inhibition may be occurring in parallel with the RVM.

At both an individual (Table 1) and group (Fig. 4) level, we observed areas of activation within the brainstem, the locality of which correlate with regions known to be involved in pain processing in animals (Traub et al., 1996; Rodella et al., 1998; Millan, 2002; Monnikes et al., 2003) and human studies (Boivie and Meyerson, 1982; Baskin et al., 1986; Tracey et al., 2002). The spatial pattern of activation was very similar for visceral and somatic pain; however, greater quantitative bilateral NCF activity during visceral pain compared with somatic was found. Furthermore, correlation of right PAG activity with anxiety was observed only during visceral pain. We propose that these findings may underlie some of the perceptual and behavioral differences observed between the two sensory modalities. Noninvasive tools are essential if we are to gain additional understanding of the complex interactions of the human brainstem nuclei (Iannetti and Tracey, 2005), not only in pain, but for all sensory, autonomic, and motor functions in humans in health and disease.

Footnotes

P.D. and L.C. are funded through Grant AR41622 from the National Institute of Arthritis and Musculoskeletal and Skin Diseases-National Institutes of Health (Bethesda, MD). R. G. W. is funded through Advanced Training Fellowship 067037 from The Wellcome Trust. Q.A. is funded by the Medical Research Council and Cancer Research Campaign (London, UK). I.T. is funded by the Higher Education Funding Council for England and Medical Research Council of Great Britain and Northern Ireland (FMRIB Centre). We acknowledge the Royal Society for an equipment grant (RSRG 24647). We thank Andreas Schoenfeld of Standard Instruments (Karlsruhe, Germany) for technical support. We also thank Prof. Derek Jewell for his guidance and support.