Common brain activations for painful and non-painful aversive stimuli

Aversion: painful and non-painful stimuli

Identification of potentially harmful stimuli is necessary for the well-being and self-preservation of all organisms. Organisms with relatively simple nervous systems (e.g. worms, fruit flies) display motivated approach and avoidance behaviours to rewarding and aversive stimuli, respectively, implying the existence of some evolutionarily conserved mechanisms [1, 2]. Aversive stimuli are those which an organism will generally expend energy to minimize or avoid [3]; in this context, aversion is operationally opposite to reward [4]. However, the strength of aversive stimuli and the context in which they occur can produce a variety of psychophysical (e.g. negative emotion, pain) and behavioural (e.g. reduced behaviour following punishment, avoidance) responses. While recent work has suggested the existence of a common aversion-related network of brain regions responsible for the basic processing of aversive stimuli [5], this work focused only on studies involving non-painful stimuli and studies including painful stimuli were not considered; as such, it is unclear if those results extend to pain-related processing.

Pain-associated brain activity

Pain, which typically results from activating the nociceptive system (but can also involve non-nociceptive mechanisms, such as in neuropathic pain), is experienced across mammals and is critical for survival [6]. Studies in humans and non-human animals have generally supported the notion that pain is processed differentially in the brain according to affective (e.g. amygdala, anterior insula, hippocampus) and sensory (e.g. somatosensory cortices, posterior insula) dimensions (e.g. [7–9]; though see also [10] for a review on the influential 3-dimension theory of pain). Nonetheless, the assumption that this network (sometimes referred to as the ‘Pain Matrix’) is specifically activated by painful stimuli has been questioned [11–13]. Using fMRI in humans, Mouraux et al. (2011) uncovered strong support for the notion that the typical regions of the Pain Matrix are largely involved in salience processing [12]. They showed that multimodal non-painful aversive stimuli and painful stimuli activate similar regions in the MCC, insula, thalamus, and sensory cortex and that the BOLD signals in these regions correlated largely with the perceived saliency of the stimulus (regardless of modality or stimulus type).

Studies suggest shared regions for pain and non-pain aversion

While much work has identified pain as a uniquely important experience (e.g. [10, 14]), numerous studies using non-painful aversive stimuli (e.g. unpleasant sounds, sights, etc.) have implicated many of the same cortical and subcortical regions [5], suggesting that the processing of various painful and non-painful aversive stimuli require many of the same neurobiological substrates. In this regard, human studies have been key to understanding the role of cortical regions (e.g. prefrontal and insular cortices; e.g. [15, 16]). Alternately, studies in animals have highlighted the importance of subcortical areas such as the periaqueductal grey, hypothalamus, bed nucleus of the stria terminalis, nucleus accumbens/ventral striatum and ventral tegmental area [17–20]. While prior work has identified a network of regions involved in non-painful aversion-related processing [5], it nonetheless remains unclear which, if any, of those identified areas are also involved in processing painful stimuli.

Systematic translational analysis of aversion-related circuitry

The present hypothesis is that there exists a core aversion-related circuit involved in processing aversive stimuli regardless of whether they are painful or non-painful. In an analogical sense, this network would be similar to the basic underlying (e.g. mesocorticolimbic) circuitry identified in the field of reward [21–23]. Prior meta-analyses in humans have outlined core regions associated with pain processing [6, 24], and some animal work has even suggested the existence of an overlapping pain and non-pain-related aversion network [25]. Nonetheless, no investigations have used both human and animal data to directly explore the possibility of a shared network for pain- and non-pain-related processing.

To this end, a translational cross-species approach was used to identify the core components of the potential aversion-related network. More specifically, our first aim was to compare brain activations to the passive reception of painful aversive stimuli in healthy adults (using a meta-analysis of human imaging data; i.e. functional magnetic resonance imaging, fMRI, or positron emission tomography, PET) to those in rodents (using a systematic review of studies including markers of cellular activation and available imaging studies). Secondly, we aimed to compare the results on pain-related processing to those gathered previously on the processing of passive non-painful aversive stimuli in humans (meta-analysis) and animals (systematic review) [5].

Our main hypothesis is that aversive stimuli, regardless of origin (e.g. sensory modality) or perception (e.g. painful or non-painful), are processed largely by a common network of brain regions. However, some areas may be more (or uniquely) involved in different aspects of pain- and non-pain-related aversive processing. The use of a meta-analytical approach allows for the clear distinction of areas which have been identified reliably across numerous studies – in comparison to individual studies which may have low power and a higher probability of reporting false positive activations [26]. The incorporation of animal studies allows for a cross-species comparison and ensures that especially subcortical areas, which may be important for aversion-related processing, are identified. Studying these areas in humans has proven difficult given limitations in optimal imaging resolution and the correct interpretation of subcortical activations (or the lack thereof) [27]. Importantly, this translational approach allowed for the direct comparison of the overlap between areas identified in pain and non-pain aversion studies.

Pain-related activation in humans (meta-analysis) and rodents (systematic review)

Results of the meta-analysis revealed a general pain-related brain circuitry involving the bilateral insula, mid cingulate cortex (MCC), postcentral gyrus (primary and secondary somatosensory cortices), precentral gyrus (motor cortex), secondary/supplementary motor area (SMA), and thalamus (Thal). Additional extent-based clusters, extending from regions with peak activations, were also noted in the anterior cingulate cortex (ACC), posterior cingulate cortex (PCC), dorsomedial prefrontal cortex (DMPFC), bilateral operculum, bilateral supramarginal gyri, right ventrolateral orbitofrontal cortex (VLOFC), right rostral temporal gyrus (RTG), right hippocampal/parahippocampal area (Hipp/Parahipp), inferior frontal gyrus, dorsal striatum (DS), cerebellar cortex, and midbrain and rostral pons (Figure 1 and Additional file 1: Table S1A).

Non-pain-related aversive brain activation in humans and other animals

As published previously by Hayes & Northoff (2011), the meta-analysis results of human neuroimaging studies using passive non-painful aversive stimuli implicated brain circuitry involving the amygdala (Amyg), anterior insula (AI), ventrolateral orbitofrontal cortex (VLOFC), hippocampus (Hipp), and parahippocampal gyrus (Parahipp), dorsal striatum (DS), rostral temporal gyri (RTG), and thalamus (Thal). Extent-based clusters were also noted in the anterior and middle cingulate cortex (ACC and MCC), dorsomedial prefrontal cortex (DMPFC), secondary motor area (SMA), and midbrain (Figure 1 and Additional file 1: Table S1B) [5].

Comparison between pain- and non-pain-related activations in humans and animals

Conjunction and contrast analyses in humans

A conjunction analysis of the human meta-analysis results for pain- and non-pain-related aversive stimuli revealed a common network of brain areas including: MCC, posterior cingulate cortex (PCC), AI/claustrum, right VLOFC, DMPFC, right RTG, SMA, thalamus, right Hipp/Parahipp, dorsal striatum, and midbrain (see Figure 2, Table 3). We used the conservative minimum statistic resulting in the intersection of 5900 voxels at the cluster-based FWE corrected level [28]. While all clusters are displayed in Table 3, only those greater than 10 voxels were considered in the discussion. The percentage of overlapping voxels in the present study accounted for 24% of pain-related activations and 35% of aversion-related activations.

Cluster	x	y	z	voxels	Volume (mm3)
1	16	-3	-1	4019	32152
2	1	9	39	1418	11344
3	-33	17	2	400	3200
4	-2	-23	36	20	160
5	7	-38	-21	20	160
6	-21	12	-3	5	40
7	-24	-6	-9	3	24
8	-16	15	-4	2	16
9	4	-25	38	2	16
10	9	-30	-20	2	16
11	-12	18	8	1	8
12	-18	-28	0	1	8
13	8	0	-6	1	8
14	-32	-2	-10	1	8
15	-34	4	-12	1	8
16	-12	12	8	1	8
17	-20	6	4	1	8
18	18	16	0	1	8
19	4	-28	-20	1	8

Results of MKDA analysis: Overlap for painful and non-painful stimuli highlighted in the bottom row of Figure 2. All results are family-wise error rate whole-brain corrected at p < 0.05. Only clusters >10 voxels were considered.

In addition to the activation of this common network, the presentation of painful compared to non-painful stimuli (painful > non-painful) resulted in additional unique activations in the primary and secondary sensorimotor cortices, posterior cingulate cortex (PCC), the bilateral operculum (including bilateral SMG), cerebellum, and the rostral midbrain (in a region consistent with the VTA; see [29]) and rostral pons (see Figure 2).

Alternately, non-painful > painful aversive stimuli resulted in unique activations in the amygdalae and left hemisphere VLOFC, RTG, and Hipp/Parahipp. In addition, extended activation clusters unique to non-painful stimulus presentation were noted for the more anterior portions of the ACC (pregenual and subgenual ACC), the DMPFC, a more anterior portion of the PCC/posterior MCC, an area encompassing the hypothalamus and dorsal (DS) and ventral striatal (VS) areas, left SMA, and bilateral VLOFC (Figure 2).

Comparison to non-human animals

While a conjunction analysis was not performed across non-human animal studies (given that precise coordinates related to neural/cellular activation are less often given), the pattern of activation seen in humans was reflected in these studies. The general pattern of brain activation between humans and animals appears to be maintained, as is generally reflected by the relative percentage of studies reporting such activations (see Tables 3 and 3). For instance, a large percentage of animal studies have reported cingulate activations in pain-related studies (Table 1) which is mirrored by the large activation in the pain (but not the non-pain) meta-analysis; in comparison, a majority of animal studies have reported amygdala activations in response to non-painful aversive stimuli (Table 2), which is reflected by the large activation in the non-pain (but not the pain) meta-analysis results.

Additional subcortical regions of importance (which may be too difficult to detect or differentiate precisely using human brain imaging techniques) have also been noted in animals (see Tables 1, 2). Most seem to be involved in general aversion-related processing (i.e. BNST, Hab, Hyp, LC, NAc, NTS, PAG, PBN, septal area, anterior raphe nuclei) while this investigation suggests that some may be involved more in pain- (i.e. IC, SC, pretectal area, red nucleus, retrosplenial area) or non-pain-related (i.e. basolateral and central nuclei of the amygdalae, VTA) processing. Results for overlapping and non-overlapping activations across species are illustrated in Figure 3.

Areas activated by painful and non-painful aversive stimuli are similar across species

Results from the pain-related meta-analysis (Figure 1, Additional file 1: Table S1A) and animal literature review (Table 1, Additional file 1: Table S2A) are consistent with prior human [30, 31] and animal studies [32, 33]. Similar regions were noted here, including: bilateral insula, MCC, primary and secondary motor and somatosensory cortices, VLOFC, Hipp/Parahipp, DS, Thal and midbrain. It is also worth noting that our results are nearly identical to another recent human meta-analysis on pain [24]. Additional, mainly subcortical, areas identified across animal studies included the Amyg, BNST, Hab, Hyp, NAc, NTS, PBN, PAG, SC, IC, and septal nuclei. Brain imaging techniques may be limited by low subcortical resolution (see [34] for a brief discussion of fMRI capabilities), however some fMRI studies focusing on this level (e.g. PAG, PBN, VTA) have corroborated their involvement in human pain processing (e.g. [35]).

As reported previously [5], non-painful aversive stimuli result in consistent activations in humans (Figure 1, Additional file 1: Table S1B) and other animals (Table 2, Additional file 1: Table S2B). The main activations to non-painful aversive stimuli across species include: Amyg, ACC, VLOFC, DMPFC, secondary motor cortex, Hipp/Parahipp, DS, RTG, Thal, and midbrain. Additional subcortical activations, identified most consistently in animal studies (although see human imaging studies for instance by [36] and [37]), included the BNST, Hab, Hyp, NTS, NAc, PAG, PBN and septal nuclei. These results are also consistent with prior studies across species looking at aversion-related concepts such as fear [38, 39], threat [40, 41] and social punishment [42, 43]).

Together, these findings suggested the involvement of similar regions in pain- and non-pain-related aversive processing (Figure 1). However, the specific substrates involved remain unclear as does whether the apparent spatial overlap involved in both processing types reflects an actual overlap in function. Iannetti and Mouraux (2010) have argued that the so-called Pain Matrix (reflecting a pain-specific network of brain regions; see [44, 45] for early uses of the term) may be misleading. While the Pain Matrix is activated reliably across studies (as underscored here and in other meta-analsyses, e.g. [24]), and the magnitude of its response is highly correlated with the intensity of pain perception [46, 47], there is, nonetheless, presently no evidence that any single brain region or clearly defined network is devoted to the processing of pain.

Iannetti and Mouraux (2010) (and other authors, e.g. [13]) subsequently suggested that, instead of pain-specific processing, this system may be better regarded as a multimodal network related to the detection of saliency [11]. The notion that this network is largely multimodal is supported by prior work demonstrating activations in these regions regardless of the sensory modality in which the aversive stimuli were presented [5, 12]. As well, this notion is further supported by the present findings that painful stimuli (which are all tactile in nature) show overlapping activations with mainly non-tactile, non-painful, aversive stimuli. In particular, in a series of fMRI experiments, Mouraux et al. (2011) showed that most of the activations related to painful stimuli appear to also be involved in non-painful processes [12].

In another recent example, fMRI work by Moulton et al., (2011) showed similar overlapping between activations related to non-painful (i.e. unpleasant pictures) and painful (i.e. noxious heat) aversive stimuli in the cerebellum; these signals were correlated with many of the same regions noted in the present study (e.g. Hyp, ACC, Hipp/Parahipp) [48]. Interestingly, it is difficult to assess the cerebellum’s potential involvement in general aversion-related processing given that it is often overlooked in imaging studies and most animal studies use painful shock-tone pairing to investigate its role. As such, while its role in general conditioning mechanisms is well established (e.g. [49]), fewer studies have looked at non-pain-related aversive stimuli (though note research showing its involvement in conditioned taste aversion; e.g. [50]).

As the frequent occurrence of similar brain activity appeared to be further supported by the initial cross-species results above, the second aim of this study was to specifically compare the potential overlap between pain- and non-pain-related activations in humans and other animals.

Aversion-related stimuli activate a common core network

The conjunction analysis of pain- and non-pain-related aversion studies in humans (Figure 2, Table 3) and the systematic review in animals (Tables 1 & 2; Additional file 1: Tables S2A & S2B) revealed a common core aversion-related network consisting of MCC, PCC, AI/claustrum, right VLOFC, DMPFC, SMA, right Hipp/Parahipp, Thal, DS, and midbrain (including an area encompassing the PAG). Additional, mainly subcortical, areas were noted in animals including the Amyg, BNST, Hab, Hyp, NAc, septal nuclei, NTS, and PBN. This core, cross-species, aversion-related network (illustrated in Figure 3) suggests that passively received aversive stimuli are, at least in part, processed in a similar way using similar neurobiological substrates.

Identified consistently here, both the cingulate (especially the MCC) and bilateral AI have been implicated in many functions though there is still much debate on whether a basic role for these regions exists. For instance, research in humans and other animals suggests they are involved in detecting and processing errors [51, 52], in reward-related processing [53, 54], and in adaptive decision making [55, 56]. In fact, given the roles that these, and other, regions play across various networks (e.g. salience, interoceptive, resting state, and valuative), it is possible that they are active participants in many underlying processes. For instance, various lines of evidence in humans and primates have implicated the anterior portion of the MCC as being especially key for the integration of negative affect, pain, and cognitive control (as reviewed in [57]). Moreover, despite our attempt to reduce the impact of cognitive processes by limiting our studies to those involving passive stimuli, it is important to note that implicit cognitive processing (e.g. brain activity related to preparatory escape, emotional regulation) may contribute to the activations noted across species.

Differential activations related to painful and non-painful aversive stimuli

Although a core network of brain regions for the processing of painful and non-painful aversive stimuli was identified, the areas of activation noted in the human meta-analysis do not overlap completely. In fact, there are a number of regions which appear to be unique to painful or non-painful stimuli. Specifically, painful stimuli resulted in unique activations in the primary and secondary sensory cortices, the PCC, bilateral operculum/mid and posterior insula, and cerebellum. Non-painful aversive stimuli resulted in unique activations in the AMYG, L VLOFC, L RTG, and L Hipp/Parahipp as well as additional extended activations in the pACC, DMPFC, HYP, posterior MCC, DS/VS, and bilateral VLOFC. Similar subcortical regions were identified in animal studies of both pain- and non-pain-related processing (e.g. NTS, PBN, BNST). Nonetheless, from the animal studies included here, a few areas are noted in only pain- (i.e. red nucleus, pretectal area, IC/SC) or non-pain- (i.e. VTA) related studies.

The differences may in fact be unique, thus providing unique neural signatures for painful and non-painful aversive stimuli. However, there are alternative explanations for these apparently unique activations which seem more tenable. Firstly, the differential weighting of activations (i.e. that some regions are more involved) may help explain the absence of activations in human studies on pain (e.g. amygdala) compared to non-pain aversion. Secondly, as the analysis of imaging studies on pain and non-pain aversion were limited to the passive reception period, this is equivalent to taking a temporal snapshot of brain activity. While intended to be similar across pain and non-pain studies, this snapshot may in fact reflect unique brain processing for two types of stimuli differentiated by their timing properties (i.e. painful stimuli result in fast aversive responding; non-painful stimuli are slower and more variable).

Strengths and limitations

The greatest strength of the current work is its translational nature. The inclusion of animal studies has two main advantages. Firstly, they help to support the findings from human imaging and add insight regarding subregional differences and underlying mechanisms. Secondly, they underscore the involvement of many subcortical regions which are generally underreported in imaging studies (as noted previously). We believe using this approach somewhat offsets the potential selection biases which may be found across human imaging studies (e.g. lowering statistical thresholds for a priori regions) and animal studies (e.g. looking only for activity in a priori regions; also see below).

The strict and narrow criteria used in both human and animal studies allowed for a clearer interpretation of results (e.g. the use of passive and acute aversive stimuli only; the exclusion of studies/subjects using explicit cognitive tasks; see Methods section for all criteria). These criteria were used to isolate, as clearly as the present methods allow, the period of brain activation during which acute aversive stimuli are present (e.g. most neuroimaging studies look at periods around 5–10 seconds; the animal studies included here extract the brains as soon as possible following stimulus presentation) – thus, attempting to separate this period from others (e.g. anticipation, termination). It is in this sense that we have attempted to identify a network associated with aversion-related processing (see also [26] for further discussion on using meta-analyses to identify functionally related brain regions). Nonetheless, it is worth pointing out again (as discussed briefly above in the Differential weighting and Differential temporal dynamics sections) that the inference of a temporal relationship between regional activations relies heavily on the inclusion, and exclusion, of appropriate studies. Ultimately, the identification of such networks through meta-analyses and systematic review should be used as the basis for testing future hypotheses regarding co-activation.

This approach also shed some light on one inherent and important limitation of meta-analyses – particularly those using human imaging studies. While some neuroimaging studies do report subcortical activations in aversion-related processing (e.g. [37, 97, 98]), their relative scarcity means some subcortical regions may not be noted in the final meta-analysis results. This absence of activation likely also extends to highly variable cortical regions. The corollary is that meta-analysis results underscore the most consistent nodes of activation across studies (with the coordinates being more informative than the size or shape of the clusters per se; see also Differential weighting discussion above), while regions not identified may still be active (and even essential) components – findings that are made clearer through animal studies looking directly at brain tissue.

Although the results of the animal studies outlined in Tables 1 and 2 (i.e. listing the percentage of reported brain activations) should be considered illustrative due to reporting and researcher interest biases and the lack of a whole-brain approach (as insisted upon for the imaging data), most studies investigated at least 5 brain regions. In fact, only 9 of the 42 non-pain-related studies [89, 90, 99–106] and 13 of the 32 pain-related studies [87, 107–118] focused on 4 or less regions. In addition, none of the pain-related studies focused solely on the cingulate (a key node identified in both human and animal data), and only 3 studies from the non-pain-related aversion studies [100, 105, 106] focused solely on the amygdala (perhaps the single best described aversion-related region). Although the issue of selection bias and the reporting of positive data (the so-called file-drawer problem) cannot be fully accounted for, taken together, the animal and human data allow for a more confident interpretation regarding the inclusion of brain areas involved in aversion-related processing.

Painful and non-painful aversion-related brain activation in humans

Literature search: We identified all imaging studies – positron emission tomography (PET) and functional magnetic resonance imaging (fMRI) – published from 2000 to August 2011 with PubMed (http://www.pubmed.gov) and Web of Science (http://apps.webofknowledge.com; though no additional studies were found here) searches. Keywords included "aversion", "aversive", “avoidance”, “punishment”, “reinforcement” (to capture some studies focusing on reward, but also using an independently analysed aversive control condition), “fear”, “anger”, “disgust”, “sadness”, “negative emotion”, “pain”, “nociceptive”, “unpleasant”, "positron emission tomography", and "functional magnetic resonance imaging" and others. Furthermore, we searched the reference list of articles and reviews, including meta-analyses (for example see [6, 119, 120]). The data regarding brain activations associated with non-painful aversive stimuli were reported previously in [5].

Inclusion and exclusion criteria: Our main goal was to compare the basic brain activity of painful- and non-painful aversive stimuli. Therefore, we included only those studies which used the passive presentation of acute aversive stimuli (e.g. the viewing of unpleasant pictures; exposure to painful stimuli) without active responses. Without behavioural measures, aversive experiences were determined subjectively and often supported through physiological measures such as electrodermal activity. Designs whose contrasts did not include specific comparisons relevant to the current analysis (i.e. involving the passive reception of painful or non-painful aversive stimuli, independent of explicit cognitive processes, memory, or attention) were excluded. In addition, only studies that reported coordinates from whole-brain analysis were included (although some studies discussed region-of-interest data, those coordinates were not included here). These criteria lead to the inclusion of studies involving both exteroceptive (e.g. pictures, shock) and interoceptive (e.g. rectal distension) aversive stimuli.

Although the related search terms were included for completeness, studies reporting the responses to specific negative emotions (e.g. sadness, anger) were excluded given their social nature as were other stimuli which may involve ambiguous interpretations such as those involving empathy, disgust, or physical contamination. For related examples involving such analyses showing results consistent with the present study, see [121, 122]. Other studies looking explicitly at social aspects of aversion, such as social exclusion, were also excluded given that they typically require complex behavioural responses and involve other potentially confounding processes such as empathy and theory of mind. Moreover, it is important to note that most studies involving conditioned fearful stimuli in both humans and animals were excluded to avoid issues related to recent learning effects and complex conditioning designs, although it should be noted that the aversive stimuli used are virtually identical to those included here and the results from such studies are largely similar to those reported here [123, 124]. Finally, studies involving manipulations in homeostatic states (e.g. hunger, thirst) were also not included given their non-acute nature (e.g. usually involving forced deprivation over time). Although we aimed to cast as large an initial net as possible for both the human and other animal studies, we must concede that, given the conceptual complexity of aversion-related terms, we may have overlooked some studies which do not employ standard terminologies.

We screened all the articles for Talairach or Montreal Neurological Institute (MNI) coordinates and tabulated the reported regional foci. We included the data of healthy subjects only. Studies including individuals with psychiatric illnesses, or a history thereof, those with volumetric abnormalities or brain injuries, those taking any medications or illicit drugs, and those belonging to a group that may result in a sample bias (e.g. war veterans) were excluded. In addition, a large number of studies were excluded due to the absence of coordinates, identification of coordinate systems, and/or incomplete statistical information. These criteria resulted in the selection of 34 non-pain-related studies which included 44 contrasts, and 27 pain-related aversion studies which included 32 contrasts (see Additional file 1 for study references included in the meta-analysis). Regions were labelled macroanatomically by the probabilistic Harvard-Oxford atlas. The nomenclature for the cingulate by Vogt (2005) was used here [125].

Multilevel kernel density analysis (MKDA) meta-analytic technique: The MKDA meta-analytic approach has been covered in depth elsewhere [26, 126]. Briefly, MKDA is a coordinate-based meta-analytic method which determines the activation probability of each voxel and contiguous voxel clusters (to create voxel-based and cluster-based study comparison maps, respectively) across the brain. Compared with other meta-analysis methods, MKDA prevents any single study reporting a large number of activations from biasing the results (i.e. the study is the unit of analysis) and weights contrasts based on the quality of the study (e.g. random vs. fixed effects) and the sample size. It should be noted that study-specific statistical thresholding is not taken into account, however, because weighting activation peaks by their respective Z-scores introduces the presupposition that the scores across studies are comparable, which is untrue (e.g. scores from studies with small sample sizes would be inflated relative to the larger population; [26] for further discussion and examples). All results are reported in MNI space; co-ordinates from studies using Talairach space were converted to MNI space in the MKDA software using the Brett transform (http://imaging.mrc-cbu.cam.ac.uk/downloads/MNI2tal/tal2mni.m). Peaks from each study were convolved with a spherical kernel of 10 mm radius; doing so ensures that multiple nearby peaks are not counted as multiple activations (additionally helping to reduce the impact of studies employing lower statistical thresholds and which typically report a greater number of activated voxels). This approach also ensures that single studies do not drive the results of the meta-analysis. The threshold for significance was determined using a Monte Carlo simulation with 3000 iterations (5000 iterations did not alter the results) and a null hypothesis which assumes that the activated regions within each map are not spatially consistent (i.e. that the cluster centers are randomized throughout the grey matter). The voxel size was 2 x 2 x 2 mm (i.e. 1 voxel = 8 mm³) and cluster sizes were all greater than 10 voxels (> 80 mm³). We have reported peak voxel-wise activations as well as peak cluster-wise activations (which include contiguous voxels significant at p < 0.001; whole-brain FWE corrected). All results are family-wise error rate whole-brain corrected at p < 0.05. Analyses were performed in Matlab 2009a (Mathworks, Naticks, MA) using MKDA software created by Tor Wager (http://www.columbia.edu/cu/psychology/tor/).

Painful and non-painful aversion-related brain activation in animals

Literature search: PubMed (http://www.pubmed.gov) and Web of Science (http://apps.webofknowledge.com; again, no additional studies were found here) searches identified rodent studies related to pain-related (~700 studies) and non-pain-related (~350 studies) aversion published in English from 2000 to August 2011. Keywords were similar to those above but also included specific terms related to animal studies such as: “immediate early genes”, “IEG”, “c-fos”, “rat”, “mice”, “monkey”, “mammal”, “avoidance”, “fear”, “threat”, “electrophysiology”. Of the total studies identified, only those clearly showing altered brain metabolism (e.g. increased/decreased c-Fos or blood oxygenated level dependent activity, or BOLD) were included in the systematic review (i.e. 34 pain and 42 non-pain aversion studies). Due to lack of methodological instruments, absence of precise standardized coordinate systems, and the wide range of experimental procedures, we did not conduct the same rigorous meta-analysis in animals as in humans. As with the human data, the non-human animal data regarding the presentation of passive non-painful aversive stimuli were reported previously in [5].

Inclusion and exclusion criteria: We looked at the following metabolic indexes of non-human animal brain activity: immediate early gene activation (e.g. c-Fos or Fos-like expression), BOLD activity in fMRI, [14 C]-2-deoxyglucose, and [14 C]-iodoantipyrine. Each of these indexes has previously been related to neural activity and/or metabolism. Considering the broad spectrum of animal models of aversion-related behaviour (e.g. formalin-induced nociception, foot shock, conditioned taste aversion etc.), we looked at all those data that report clear effects in brain activity between control animals and those exposed to painful and/or non-painful aversive stimuli. It is important to note that some methods must necessarily use between-subjects measures (e.g. counting IEG-stained cells in brain slices) whereas others use mainly within-subjects measures (e.g. neuroimaging).

Studies involving adolescent animals, chronic exposure to aversive stimuli, and exposure to drugs of abuse and those having a direct effect on aversion- or reward-related brain circuitry, were excluded (although non-drug-exposed controls were included where appropriate). This was done in order to avoid confounding issues related to neurodevelopment and drug interactions and/or drug-induced changes in brain structure or function unrelated to the acute aversive treatment. Studies using electrical/chemical lesions or other irreversible alterations (e.g. the use of knock-out or transgenic rodents, animals bred for psychiatric disorder-related phenotypes; although one study using rats bred for high or low anxiety levels, though not pre-exposed to anxiolytic stimuli, was included [127] as the results following an aversive probe largely converged regardless of anxiety levels) were also excluded for clarity.

Any comparison between human and animal data raises the question of homology of brain regions. Since they show analogous anatomy, comparisons of subcortical regions are not often an issue (see also [128]). In contrast, the issue of homology becomes more problematic in the case of cortical regions that show both anatomical and terminological differences between humans and animals. Nonetheless, even areas which may be considered largely ‘higher-order’ or evolutionarily more recent, such as the prefrontal cortex, may show strong structural and functional homologies between primates and other mammals, such as rodents [129, 130]. Concerning cortical regions, we relied on criteria of homology as established by various authors [131–133]. In addition, it should be noted that due to the strict inclusion/exclusion criteria nearly all of the animal studies involve rodents with the exception of one study using monkeys [102]. As such, the generalizability to other mammals should be considered cautiously. Nonetheless, beyond the studies included here, studies involving primates were also considered in the discussion to add support to the translatability of the results (e.g.) [53, 66].